Numerical micromechanical investigations of the mechanical behavior and damage evolution of glass fiber reinforced composites are presented. A program code for the automatic generation of 3D micromechanical unit cell models of composites with damageable elements is developed, and used in the numerical experiments. The effect of the statistical variability of fiber strengths, viscosity of the polymer matrix as well as the interaction between the damage processes in matrix, fibers and interface are investigated numerically. It is demonstrated that fibers with constant strength ensure higher strength of a composite at the pre-critical load, while the fibers with randomly distributed strengths lead to the higher strength of the composite at post-critical loads. In the case of randomly distributed fiber strengths, the damage growth in fibers seems to be almost independent from the crack length in matrix, while the influence of matrix cracks on the beginning of fiber cracking is clearly seen for the case of the constant fiber strength. Competition between the matrix cracking and interface debonding was observed in the simulations: in the areas with intensive interface cracking, both fiber fracture and the matrix cracking are delayed. Reversely, in the area, where a long matrix crack is formed, the fiber cracking does not lead to the interface damage. (C) 2009 Elsevier Ltd. All rights reserved.