Research Papers

Predictions of Delamination Growth for Quasi-Static Loading of Composite Laminates

[+] Author and Article Information
Jiawen Xie

Department of Aerospace Engineering,
University of Michigan,
Ann Arbor, MI 48109-2140
e-mail: jwxie@umich.edu

Anthony M. Waas

Felix Pawlowski Collegiate Professor, Emeritus
University of Michigan,
Ann Arbor, MI 48109-2140
Boeing Egtvedt Endowed Chair & Chair
William E. Boeing Department of Aeronautics
and Astronautics,
University of Washington,
Seattle, WA 98195-2400
e-mail: awaas@aa.washington.edu

1Corresponding author.

Contributed by the Applied Mechanics Division of ASME for publication in the JOURNAL OF APPLIED MECHANICS. Manuscript received March 13, 2015; final manuscript received May 14, 2015; published online June 9, 2015. Assoc. Editor: Daining Fang.

J. Appl. Mech 82(8), 081004 (Aug 01, 2015) (12 pages) Paper No: JAM-15-1134; doi: 10.1115/1.4030684 History: Received March 13, 2015; Revised May 14, 2015; Online June 09, 2015

This paper presents an exact, two-dimensional (2D) quasi-static elastic analysis of predelaminated composite panels subject to arbitrary transverse pressure loads. The piecewise linear spring model and the shear bridging model are, respectively, used to simulate the normal contact and shear frictional behavior between the interfaces of the existing delamination. This general contact model can be further reduced to the “friction-free model” and the “constrained model” by assigning extreme values to the spring stiffnesses. The analysis yields a closed-form solution for the 2D displacement and stress fields. To predict the delamination propagation, different propagation criteria as suggested in the literature are used. Calculated load–displacement responses and delamination threshold loads are in good agreement with existing experimental data. The results are further compared against simple fracture models and a model that uses a modified classical laminated plate theory for the predelaminated composite with fracture criteria, showing an overestimation of the delamination threshold loads. The 2D elasticity theory that is formulated can be used with confidence to study other multilayered structures with multiple delaminations and subject to arbitrary loading profiles.

Copyright © 2015 by ASME
Your Session has timed out. Please sign back in to continue.


Abrate, S., 1998, Impact on Composite Structures, Cambridge University Press, New York.
Davies, G. A. O., and Robinson, P., 1992, “Predicting Failure by Debonding/Delamination,” AGARD 74th Structures and Materials Meeting, Debonding/Delamination of Composites, Patras, Greece, May 24–29, Paper No. 5.
Schoeppner, G. A., and Abrate, S., 2000, “Delamination Threshold Loads for Low Velocity Impact on Composite Laminates,” Composites, Part A, 31(9), pp. 903–915. [CrossRef]
Olsson, R., Donadon, M. V., and Falzon, B. G., 2006, “Delamination Threshold Load for Dynamic Impact on Plates,” Int. J. Solids Struct., 43(10), pp. 3124–3141. [CrossRef]
Kruntcheva, M. R., Collop, A. C., and Thom, N. H., 2005, “Effect of Bond Condition on Flexible Pavement Performance,” J. Transp. Eng., 131(11), pp. 880–888. [CrossRef]
Hu, X., and Walubita, L. F., 2010, “Effects of Layer Interfacial Bonding Conditions on the Mechanistic Responses in Asphalt Pavements,” J. Transp. Eng., 137(1), pp. 28–36. [CrossRef]
Abdul-Baqi, A., and Van der Giessen, E., 2001, “Indentation-Induced Interface Delamination of a Strong Film on a Ductile Substrate,” Thin Solid Films, 381(1), pp. 143–154. [CrossRef]
Sjoblom, P. O., Hartness, J. T., and Cordell, T. M., 1988, “On Low-Velocity Impact Testing of Composite Materials,” J. Compos. Mater., 22(1), pp. 30–52. [CrossRef]
Razi, H., and Kobayashi, A. S., 1993, “Delamination in Cross-Ply Laminated Composite Subjected to Low-Velocity Impact,” AIAA J., 31(8), pp. 1498–1502. [CrossRef]
Zhang, X., Bianchi, F., and Liu, H., 2012, “Predicting Low-Velocity Impact Damage in Composites by a Quasi-Static Load Model With Cohesive Interface Elements,” Aeronaut. J., 116(1186), pp. 1367–1381. http://aerosociety.com/News/Publications/Aero-Journal/Online/939/Predicting-lowvelocity-impact-damage-in-composites-by-a-quasistatic-load-model-with-cohesive-interface-elements
Pagano, N. J., 1969, “Exact Solution for Composite Laminates in Cylindrical Bending,” J. Compos. Mater., 3(3), pp. 398–411. [CrossRef]
Pagano, N. J., 1970, “Exact Solutions for Rectangular Bidirectional Composites and Sandwich Plates,” J. Compos. Mater., 4(1), pp. 20–34. [CrossRef]
Xie, J., and Waas, A. M., 2013, “2D Elastodynamic Solution for the Impact Response of Laminated Composites,” ASME J. Appl. Mech., 81(4), p. 041015. [CrossRef]
Qiao, P., and Chen, F., 2012, “On the Improved Dynamic Analysis of Delaminated Beam,” J. Sound Vib., 331(5), pp. 1143–1163. [CrossRef]
Brewer, J. C., and Lagace, P. A., 1988, “Quadratic Stress Criterion for Initiation of Delamination,” J. Compos. Mater., 22(12), pp. 1141–1155. [CrossRef]
Luo, H., and Hanagud, S., 2000, “Dynamics of Delaminated Beams,” Int. J. Solids Struct., 37(10), pp. 1501–1519. [CrossRef]
Rybicki, E. F., and Kanninen, M. F., 1977, “A Finite Element Calculation of Stress Intensity Factors by a Modified Crack Closure Integral,” Eng. Fract. Mech., 9(4), pp. 931–938. [CrossRef]
Davies, G. A. O., and Zhang, X., 1995, “Impact Damage Prediction in Carbon Composite Structures,” Int. J. Impact Eng., 16(1), pp. 149–170. [CrossRef]
Davies, G. A. O., Robinson, P., Robson, J., and Eady, D., 1997, “Shear Driven Delamination Propagation in Two Dimensions,” Composites Part A, 28(8), pp. 757–765. [CrossRef]
Pardoen, G. C., 1989, “Effect of Delamination on the Natural Frequencies of Composite Laminates,” J. Compos. Mater., 23(12), pp. 1200–1215. [CrossRef]
Xie, D., and Waas, A. M., 2006, “Discrete Cohesive Zone Model for Mixed-Mode Fracture Using Finite Element Analysis,” Eng. Fract. Mech., 73(13), pp. 1783–1796. [CrossRef]
Amaro, A. M., Reis, P. N. B., and De Moura, M. F. S. F., 2011, “Delamination Effect on Bending Behaviour in Carbon–Epoxy Composites,” Strain, 47(2), pp. 203–208. [CrossRef]
Lammerant, L., and Verpoest, I., 1996, “Modelling of the Interaction Between Matrix Cracks and Delaminations During Impact of Composite Plates,” Compos. Sci. Technol., 56(10), pp. 1171–1178. [CrossRef]
Prabhakar, P., and Waas, A. M., 2013, “Micromechanical Modeling to Determine the Compressive Strength and Failure Mode Interaction of Multidirectional Laminates,” Composites, Part A, 50, pp. 11–21. [CrossRef]


Grahic Jump Location
Fig. 1

Two-dimensional illustration of the predelaminated composite panel. The panel is assumed in a state of plane strain in the xz plane and simply supported at its left and right ends.

Grahic Jump Location
Fig. 2

The piecewise linear spring model, which is a slightly modified version of Ref. [16]

Grahic Jump Location
Fig. 5

Free body diagrams of modified CLT

Grahic Jump Location
Fig. 3

Three-point bend test setup

Grahic Jump Location
Fig. 4

The four-section partition of a predelaminated composite panel in modified CLT

Grahic Jump Location
Fig. 6

Simple fracture models. The model only considers the potential energy change of the center delaminated region form pristine to delaminated configuration.

Grahic Jump Location
Fig. 13

The load–displacement response for 12-ply specimen with delamination length of 0.5L. The delamination is located at: (a) upper-quarter-plane; (b) lower-quarter-plane; and (c) mid-plane.

Grahic Jump Location
Fig. 7

Transverse stress distributions at delaminated interface of specimens with midplane delamination length 0.3L. The dashed line represents the location of crack tip.

Grahic Jump Location
Fig. 8

Transverse shear stress distribution near the crack tip at delaminated interface of specimens with midplane delamination length of 0.3L. The dashed line represents the location of crack tip.

Grahic Jump Location
Fig. 9

The convergence study on average transverse stresses of 32-ply specimen with midplane delamination length. d/L is the relative delamination length.

Grahic Jump Location
Fig. 10

The variation of delamination threshold load versus existing center delamination length of 12-ply orthotropic laminated panel. The delamination is located at: (a) upper-quarter-plane; (b) mid-plane; and (c) lower-quarter-plane. Two-dimensional elastic solutions use the combination of three contact models (M1: general model, M2: friction-free model, and M3: constrained model) and two delamination criteria (stress: quadratic stress criterion and energy: fracture mechanics based criterion).

Grahic Jump Location
Fig. 11

The variation of delamination threshold load versus existing center delamination length of 32-ply cross-ply laminated panel

Grahic Jump Location
Fig. 12

The errors between simple fracture models and 2D elasticity theory using general contact model and fracture mechanics based failure criterion

Grahic Jump Location
Fig. 14

The load–displacement response variation with delamination length for 12-ply and 32-ply specimen with midplane delamination



Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In