Abstract: A three-dimensional non-linear progressive damage model was developed to predict the response of composite laminates with a single bolt . The analysis took into account the contact at the bolt-hole interface by Augmented Lagrangian Contact Algorithm. The objectives of the study were to investigate the onset and development of thefibre damage , mat rix damage , fibre-mat rix shear , and their relationship. From the load-displacement curves drawnfrom the calculation , the initial bearing damage strength was predicted. The effect s of stacking sequences and geomet rical parameters on the bearing st rength were determined experimentally. Good agreement between experimentalresult s and numerical predictions is observed.