A method is proposed which can be used to obtain the fibesr content distribution of compression molded long fiber-reinforced thermoplastic sheet for nonisothermal state. The fiber is modelled to be a sphere. Once the one-dimensional unsteady state heat conduction equation in solved, the mean temperature in defined across the thickness direction. The viscosity of matrix is determined with the mean temperature. Using the obtained viscosity, two-dimensional sheet0like part compression molding is simulated with the finite element method. Comparison with experiments shows that the method accurately predicts the distribution.