COMSOL 6.4 - Mixed-Mode Debonding of a Laminated Composite

Mixed-Mode Debonding of a Laminated Composite

Interfacial failure by delamination or debonding is one of the main failure modes of laminate structures. Such failures can be simulated with a cohesive zone model (CZM). A key ingredient of a cohesive zone model is a traction-separation law that describes the softening in the cohesive zone near the delamination tip. This example shows how to model debonding using the decohesion model in COMSOL Multiphysics. The capabilities of the CZM to predict mixed-mode softening and delamination propagation are demonstrated in a model of the mixed-mode bending of a composite beam.

Model Definition

Cohesive Zone Model (CZM)

The CZM is in this example defined using the displacement-based damage model available in the node under . The model is used to predict crack propagation at the interface of a laminated composite beam under mixed-mode loading. The material properties needed for this constitutive model are summarized in Table 1.

Table 1: Summary of material properties of the CZM interface. The values are for AS4/PEEK.
Property	Symbol	Value
Normal tensile strength	σt	80 MPa
Shear strength	σs	100 MPa
Penalty stiffness	pn	106 N/mm3
Critical energy release rate, tension	Gct	969 J/m2
Critical energy release rate, shear	Gcs	1719 J/m2
Exponent of the Benzeggagh and Kenane (B-K) criterion	α	2.284

The CZM is defined using a bilinear traction-separation law. Traction increases linearly with a stiffness pn until the opening crack reaches a damage initiation displacement u0. When the crack opens beyond u0, the material softens irreversibly and the stiffness decreases as a function of increasing damage d. The material fails once the stiffness has decreased to zero, that is, when d = 1. This happens at the ultimate displacement uf.

The values of u0 and uf depend on whether the separation displacement is normal (mode I) or tangential (mode II and III) to an interface. For the mixed mode, a combination is used. For the displacement based damage model, two different criteria are available to define this combination. Here the model by Benzeggagh and Kenane is used.

Mixed-mode bending of a laminated composite beam

A commonly used method to measure the delamination resistance of composite materials is the mixed-mode bending (MMB) test, see Ref. 1 and Ref. 2. This experimental procedure is here modeled to demonstrate the capabilities of the CZM.

The geometry of the test specimen is illustrated in Figure 1. It consists of a beam cracked along a ply interface halfway through its thickness. The initial crack length is cl. The beam is supported at the outermost bottom edges. A mixed-mode bending load is produced as the result of forces applied to the top edges at the cracked end and at the center of the beam.

Because of the symmetry, only half of the beam is modeled and a boundary condition is applied.

Figure 1: The geometry of the test specimen.

The material properties are those of AS4/PEEK unidirectional laminates and are listed in Table 2. The transverse isotropic linear elastic properties assume that the longitudinal direction is aligned with the global X direction. The AS4/PEEK unidirectional laminate is a built-in material in the material library.

Table 2: Laminated composite material properties.
Property	Symbol	Value
Young’s modulus, along fibers	EX	122.7 GPa
Young’s modulus, across fibers	EY=EZ	10.1 GPa
Poisson’s ratio	νYZ	0.45
Poisson’s ratio	νXY=νXZ	0.25
Shear modulus	GXY=GXZ	5.5 GPa

The beam is supported on the bottom at its outer edges. A lever that sits on top of the beam applies a load. The lever is also attached to the cracked end and swivels around a contact area at the center of the beam. The lever is pushed down at the opposite free end, thereby simultaneously applying mode I and mode II loads on the test specimen. Arbitrary ratios of mixed-mode loading can be adjusted by varying the length of the lever ll.

In this example, the lever is omitted. Instead, the forces that the lever transmits to the beam are applied directly. A pulling force Fe is acting on the cracked side of the beam. At the center, a force Fm pushes down. The desired mixed-mode ratio mm regulates the ratio of their magnitudes lr via

Further details on the background of the equation above can be found in Ref. 1 and Ref. 2.

Results and Discussion

The model is analyzed for a mixed-mode ratio of 50%. The von Mises stress distribution of the last computed parameter step is shown in Figure 2. At this step the initial crack has propagated along the interface as shown in Figure 3. The crack now extends beyond the center of the beam.

Figure 2: The von Mises stress distribution at the last computation step.

Figure 3: Plot showing the health of the laminate interface. The debonded part is shown in red, the intact part in green.

If, however, you zoom to the crack front region in the plot showing the interface health, you will notice that the color surfaces do not look perfectly regular. The reason is that the damage function has sharp or even discontinuous gradients within some elements. If you want examine the damage state in detail, then it is more relevant to look at the results in the individual integration points, since it is those values that are actually used during the solution. In Figure 4, such a plot is shown. Gauss point plots are suited for examining local details in the solution, but not for providing a general overview.

Figure 4: Interface health visualized in the integration points.

One of the outputs of the MMB test is a load-displacement curve. Both the load and displacement are measured at the endpoint of the lever that is used to apply the load to the test specimen. Since the lever is not explicitly modeled, the load-displacement data has to be deduced from the simulation results. Details of the analysis are contained in Ref. 1 and Ref. 2, with the following result.

The force Flp at the load point of the lever can be determined from the load applied to the cracked edge in the model Fe and the lengths of the test beam lb and load lever ll:

The length of the load lever above depends on the desired mode mixture mm:

ll measures the length from the center of the test specimen to the free end of the load lever.

The displacement at the load point ulp is computed from the mode I opening at the cracked edge uIe and the z-displacement at the center of the beam wc according to

The resulting load-displacement curve is shown in Figure 5. The curve confirms what Figure 3 displayed. The maximum load that the beam with the initial crack can carry is exceeded and delamination occurs. After exceeding a peak load, the load decreases until the displacement reaches around 7 mm. This point approximately corresponds to when the crack reaches the center of the specimen. Thereafter, the load starts to increase again, but with a much lower stiffness than before delamination.

Figure 5: Load-displacement curve of the MMB test at 50% mixed-mode loading.

Notes About the COMSOL Implementation

Since no large relative deformations of the laminates are expected, the node is utilized to model the behavior between laminates. This is an alternative way to model contact without requiring an assembly.

To implement a cohesive zone model in COMSOL Multiphysics, add an subnode with a subnode to the node.

To trace the solution after the peak load, the simulation must be displacement controlled. This is obtained by using a global equation, in which the distributed load is controlled by the sweep parameter, which is the displacement at the free edge.

To accurately resolve the decohesion at the interface, a dense mesh is used in the zone where delamination is expected. If the mesh elements were too large, the crack could jump several elements in one parameter step, resulting in an unstable solution.

Reference

1. P.P. Camanho, C.G. Davila, and M.F. De Moura, “Numerical Simulation of Mixed-mode Progressive Delamination in Composite Materials,” Journal of composite materials 37.16 (2003): 1415–1438.

2. J.R. Reeder, and J.R. Crews Jr., “Mixed-mode bending method for delamination testing,” AiAA Journal 28.7 (1990): 1270–1276.

Structural_Mechanics_Module/Contact_and_Friction/cohesive_zone_debonding

Modeling Instructions

From the menu, choose .

New

In the window, click