Due to the large scatters in Sheet Molding Compound (SMC) mechanical properties, the uncertainty in the prediction of mechanical responses of SMC composite structures is significantly high as compared to that for continuous fiber composite structures. Another challenge is the size effect resulting from the scatter in mechanical strengths, i.e., the measured material flexural strength is much higher than its tensile strength. As a result, classical approaches that rely on the mean tensile strength can substantially underpredict the flexural responses in Finite Element Analysis (FEA). This work investigates the use of a material model with bimodal Weibull distributions for its strength properties in Probabilistic FE (PFE) simulations. A Random Search-based optimization procedure was developed to obtain model parameters which are difficult to determine. The material model was verified by PFE simulations of standard experiments, and then examined in simulations of open-hole tension and disk bending tests.