We describe the calculation of thermal expansion coefficients of $alpha$-quartz, muscovite and albite using a MD simulation method. The selection of interatomic potentials is important for the MD calculation, and we used the 2-body interatomic potential function. The coefficients are calculated using a differential operation of the temperature dependence of the lattice constant obtained from the NPT-ensemble molecular dynamics simulation. Reasonable agreement is found between the analytical results and measured data.