Abstract: Based on the 2D model , the user subroutine USFLD has been programmed for ABAQUS to simulatethe damage evaluation in the composite laminate. The numerical result s show that the present method can predict thebear st rength and demonst rate the failure mode of the bolt-loaded laminate very well. Af ter analyzing the result s , asfor the cross-ply laminate , the shear non-linearity effect is not evident , both the non-linear analysis and the linear analysis can be applied to predict the st rength and the failure mode of the cross-ply laminate , but as for the angle-plylaminate , the difference between the above two methods is very large , and the non-linear analysis is more fit forthis laminate. The clearance between the bolt and the hole affect s the ultimate st rength of the laminate to a great extent . The bear st rength decreases with the clearance increasing. At the same time , the clearance can affect the dist ribution and degree of damage in the bolt-loaded laminate.