Abstract: Based on continuous damage theory, the crack-tip stress singularity of composite fiber-reinforced resin laminate with a Mode I crack was analyzed by finite element method and a nonlinear and multi-scalar damage model. Considering the geometrical non-linearity, a full coupling method iterating the damage and displacement simultaneously is applied to computation. The numerical results indicate that FEM is efficient to analyze stress singularity, and the stress singularity always exists in the hardening phase of material and decreases with the increment of damage in crack tip. The relationships between singularity and material properties are presented.