3-DIMENSIONAL FINITE-ELEMENT ANALYSIS OF FREE-EDGE STRESSES AND DELAMINATION OF COMPOSITE LAMINATES

Authors
Citation
Hty. Yang et Cc. He, 3-DIMENSIONAL FINITE-ELEMENT ANALYSIS OF FREE-EDGE STRESSES AND DELAMINATION OF COMPOSITE LAMINATES, Journal of composite materials, 28(15), 1994, pp. 1394-1412
Citations number
22
Categorie Soggetti
Materials Sciences, Composites
ISSN journal
00219983
Volume
28
Issue
15
Year of publication
1994
Pages
1394 - 1412
Database
ISI
SICI code
0021-9983(1994)28:15<1394:3FAOFS>2.0.ZU;2-K
Abstract
A 24-node, 64-degree-of-freedom, solid hexahedronal finite element is used to analyze the interlaminar free edge stresses and delamination o f composite laminated plates. The difficulty due to the requirement of large storage space which is frequently encountered in the modeling o f contoured delamination problems is resolved by using the preconditio ned conjugate gradient method. One of the advantages of this method is that there is no need to store the fill-in elements in the stiffness matrix decomposition process, which leads to a significant reduction i n Storage space required, thereby making the numerical analysis of the complicated contoured delamination problem feasible as well as effect ive. To evaluate the accurateness of the present formulation, numerica l algorithm and computer program, free edge stress analyses are perfor med for a [0/90], cross-ply and a [+/-45], angle-ply laminated plate. Both plates are subjected to an axial tensile strain loading. The pres ent three-dimensional results are found to be in good agreement with t hose obtained using an alternative quasi-three-dimensional methods. To further demonstrate the applicability of the present development, the strain energy release rate of a [+/-25/90], graphite-epoxy rectangula r laminated plate subjected to tensile strain is studied. Two initial delaminations are assumed to be located in the midspan of the two oppo site free edges and also at the middle of the thickness. The delaminat ion contour is assumed to be a straight line segment joined by two cir cular arcs. The strain energy release rates along the delamination fro nt are calculated using the modified crack closure method. As a step b eyond the solution using quasi-3D assumption, the present results show the non-uniform distribution of mode I strain energy release rate G(I ) along the boundary of the delamination front which fluctuates to a m aximum value when approaching both ends at the free edge of the lamina te.