A flow model of an anisotropic viscous flow is developed for the analysis of compression molding of sheet molding compounds(SMC) and penalty finite element formulation is presented to accurately predict flow front progressions. The SMC, which is substantially reinforced with fibers that are randomly oriented in planes, is assumed anisotropic Newtonian fluid which has different viscosities in extension and shear. Numerical examples are solved and the results are compared with those of experiments for the verification of the proposed flow model. With proposed model, initial charge shape is properly calculated so that no knit line develops for final L-shaped charge. Also the effects of various process conditions such as charge size, height, mold temperature, viscosity, and friction are studied.