This paper deals with the computational modeling of delamination and the pr
ediction of delamination growth in laminated composites. In the analysis of
post-buckled delaminations, an important parameter is the distribution of
the local strain-energy release rate along the delamination front. A study
using virtual crack closure technique is made for three-dimensional finite-
element models of circular delaminations embedded in woven and non-woven co
mposite laminates. The delamination is embedded at different depths along t
he thickness direction of the laminates. The issue of symmetry boundary con
ditions is discussed. It is found that fibre orientation of the plies in th
e delaminated part play an important role in the distribution of the local
strain-energy release rate. This implies that the popular use of quarter mo
dels in order to save computational effort is unjustified and will lead to
erroneous results. Comparison is made with experimental results and growth
of the delamination front with fatigue cycling is predicted. A methodology
for the prediction of delamination areas and directions using evolution cri
teria derived from test coupon data is also described. It is found that evo
lution criteria based on components of the strain-energy release rate predi
ct the rate of delamination growth much better than evolution criteria base
d on the total strain energy release rate. (C) 2001 Elsevier Science Ltd. A
ll rights reserved.