Abstract

Large deviations have been observed while analysing composite double cantilever beam (DCB) specimens assuming each cracked half as a simple cantilever beam. This paper examines the effect of rotational spring stiffness on the critical fracture energy considering nonzero slope at the crack-tip of the DCB specimen by modelling each cracked half as the spring-hinged cantilever beam. The critical load estimates of DCB specimens from are found to be in good agreement with in-house and existing test results of different composite material systems.

1. Introduction

The lightweight carbon fiber-reinforced composite (having high stiffness, strength, fatigue and impact resistance, thermal conductivity, and corrosion resistance) is found to be more economical for commercial aircraft [1]. It has become the primary material in many aircraft components (namely, fuselage, wings, and empennage components) [2]. The high strength and stiffness laminated composites in aircraft and spacecraft structures possess different properties in different directions due to their anisotropic nature. The excellent properties of composites can be utilized when the loading direction coincides with the reinforced direction. In general, the loading is multidirectional for many complex structures, which may initiate delamination and matrix cracking causing loss of structural integrity. Testing of thin skin stiffened panels of aircraft fuselage experiences bond failure at the tip of the frame flange [3]. The interlaminar fracture toughness (or the critical strain energy release rate at the onset of crack or damage propagation) is an important fracture property to be evaluated utilizing the fracture mechanics principles. Walker et al. [4, 5] have examined the tension fracture of laminates for transport fuselage, whereas Wang et al. [6] have utilized the R-curve method for the residual strength evaluation of damaged composite fuselage panel. Research and developmental activities are going on the design of the lightweight composite aircraft to sustain the discrete source of damage (see Figure 1).

The mode-I fracture energy evaluation relies heavily on the interpretation of fracture data, namely, the load-displacement records of specimens containing artificially created crack like defects. The double cantilever beam (DCB) test is widely used for evaluation of the mode-I fracture energy originally developed for studying the fracture of adhesively bonded joints [7]. Hashemi et al. [8] have evaluated the fracture energy of carbon/PEEK and carbon/epoxy composites. A good agreement is noticed between the fracture energy values from the area and compliance methods, but there is a poor agreement with those values obtained from the “load” and “displacement” methods based on simple beam theory. Several factors (namely, errors in the measurement of crack length and displacement, shear correction, and large displacement correction) are being considered to account for these discrepancies. None of these possible errors are significant enough to eliminate the discrepancies and correct the analytical method. Rao and Acharya [9, 10] have proposed a procedure for small as well as large deflections considering nonzero slope at the crack-tip of the DCB specimens.

Motivated by the work of the above researchers, this paper examines the effect of rotational spring stiffness on the critical strain energy release rate useful for accurate failure load evaluation of delaminated composites. The critical load estimates of DCB specimens from are found to be in good agreement with in-house and existing test data of different composites.

2. Theoretical Background

Composite laminates are linear elastic materials. The interlaminar fracture of composites can be dealt with the linear elastic fracture mechanics (LEFM). The double cantilever beam (DCB) test specimens are being used for mode-I fracture of unidirectional composites. Fracture energy or the critical strain energy release rate evaluation in the opening mode (mode-I) for a composite depends on the interpretation of fracture data. Normally, the fracture data consists of a load-displacement record for cracked specimens. The two different LEFM approaches are the area and compliance methods. According to the area method, the released energy due to crack extension can be evaluated by using Irwin’s formula for the critical energy release rate [11, 12]: where is the total energy of the structure and and are the width and length of the crack, respectively. For the elastic body the critical energy release rate can be determined by the formula where is the strain energy, and are the load and the displacement at crack length , and and are the respective values at crack length .

Since the strain energy for linear response is and , the energy release rate in (1) can be expressed through the compliance [11, 12] where is the critical load at which crack growth is observed, is the compliance, and is the displacement. In order to evaluate the compliance must be calculated by using the experimental data, that is, the measurements of the load and displacement at different crack length.

For the DCB test specimen the compliance calculated by the linear beam theory is given by [11, 12] where is modulus of elasticity in the direction of beam axis, is the area moment, and is half the thickness of the DCB specimen. Using (4) in (3), the critical energy release rate for DCB test can be calculated by

Equation (5) can be used for mean value of by using the data reduction with the least square technique. The analytical formula (5), which is based on the simple beam theory, can be corrected by using the fracture data. In Berry method, the compliance is approximated by a power-law: where and are determined experimentally. Using (6) in (3), one can find

According to beam theory, the parameter should be 3, but from the fracture data, another value of this parameter can be obtained. It is noted that different data reduction methods give different values, because, in the linear beam theory, some effects are not taken into account: the rotation and deflection at the crack-tip, the large displacements present in the test specimen, and stiffening effects due to the presence of the bonded end blocks. To compensate for the rotation of the adherends at the crack-tip and the flexibility of the adhesive layer, the crack length is increased with correction term . Accordingly is replaced with in (4). The compliance is proportional to . It should be noted that should be positive. Thus, in a plot of versus the crack length , the correction is found as the intercept of the extrapolation of experimental values to . This procedure is being followed in the modified beam analysis. Analysis of experimental results can also be performed by the finite element method. For calculating the energy release rates the modified virtual crack closure integral (MVCCI) can be used [12].

Hashemi et al. [8] have found good agreement between the values of from the area and compliance methods. They noticed poor agreement with those values of from the load and displacement methods based on simple beam theory. Rao and Acharya [9, 10] have considered nonzero slope at the crack-tip of the DCB specimen by modelling each cracked half as the spring-hinged cantilever beam and provided the compliance as

Here is the rotational spring constant. For the case of infinitely large support stiffness (i.e., ), (8) reduces to (4) for the case of a simple beam. In reality, infinitely large support stiffness at the crack-tip of the DCB specimen will not exist. For a realistic analysis, the rotational stiffness of the support has to be taken. The rotational spring constant in (8) has to be obtained from (8) by substituting the measured value of load and displacement for the measured crack size of the DCB specimen.

Eliminating in (8), the fracture toughness can be evaluated from the fracture data as [13, 14]

The critical load is evaluated from

The rotational spring constant in (10) is obtained from the fracture data as

The mode-I delamination fracture toughness is evaluated considering the fracture data of glass/epoxy [15]. Table 1 gives comparison of values obtained from the modified beam analysis, compliance method, finite element method, Berry method, and the present approach based on the spring-hinged cantilever beam theory. The critical energy release rates obtained from the different data reduction methods are in good agreement with each other.

3. Experimentation

Fracture toughness tests have been performed on the DCB specimens made of carbon/epoxy and carbon/PEEK composites. To create the initial crack in the composite specimens, a thin Teflon film (having 13 μm thickness) is inserted during stacking. As per the ASTM standard D5528 [16], the optimum size of the DCB specimen should have 20 to 25 mm width, above 125 mm length, and 3 to 5 mm thickness. The DCB specimens of 130 mm length, 20 mm width, and 3 mm thickness are made from carbon/epoxy and 3.25 mm thick carbon/PEEK composites. Fracture toughness tests are valid only when the initial delamination length and the thickness of DCB specimens satisfy the following criteria [14]:

The specimen surfaces are scrubbed with sand paper and cleaned thoroughly with acetone. To maintain good bonding, the base of aluminium piano hinge is scratched with file and cleaned with acetone. A thin layer of araldite adhesive is used for fixing of the piano hinge to the specimen, which is meant for application of the load. Proper care is taken while applying the araldite without covering the sides of Teflon insert. The maximum load anticipated during a DCB test is [17]

Instron3367 universal testing machine is used for testing of specimens subjected to wedge loading under displacement control. The crosshead speed is set at 1 mm/min to ensure steady crack propagation and ease of recording. The load-displacement history had been recorded by the computer. Specimen edges are coated with typewriter correction fluid and marking on the specimen started from the end of the insert to monitor the crack propagation. As per ASTM standards, the first five markings are made at an interval of 1 mm and the rest at an interval of 5 mm. The crack growth from the starter insert is observed through magnification lens and history.

Fracture data of carbon/epoxy and carbon/PEEK composites generated from DCB specimens and the values of are presented in Tables 2 and 3. The relative error in critical load estimates is found to be less than 5%, whereas it is less than 10% on displacement estimates. This demonstrates the accuracy of evaluation utilizing (9). The rotational spring stiffness value for carbon/PEEK DCB specimens is found to be very large. Hence the critical load obtained from the present analysis is close to the value obtained from (13). The crack growth resistance curve (R-curve) for the two composites is generated from the fracture data and represented by the power-law: . It is observed from Figures 2 and 3 that the critical fracture energy of carbon/PEEK is much higher than that of carbon/epoxy.

4. Critical Load Estimates of DCB Specimens

In order to examine the adequacy of the fracture toughness evaluation, fracture data [1827] of the double cantilever beam (DCB) specimens on carbon-fiber/PEEK, carbon/polyether-sulphone, carbon/epoxy, and carbon/PEK-C composites are considered. Fracture toughness of the materials and estimates of the critical load for the DCB specimens are presented.

The critical load at the onset of delamination from the critical strain energy release rate is compared with the fracture data of the materials in Figures 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, and 18. The critical load estimates at the onset of delamination are found to be in good agreement with existing test results. For the case of large rotational stiffness, the difference in critical load estimates with and without considering will be insignificant.

It is observed from Figure 19 that the crack propagation in unidirectional laminate is stable along the length of specimen, whereas in the case of multidirectional laminates (with +45°/+45° and +45°/−45° interface) the crack growth is found to be zigzag because of fiber bridging.

5. Conclusions

This paper presents a simple and reliable relation for direct evaluation of delamination fracture toughness from the fracture data of composite DCB specimens. The critical load estimates of DCB specimens from the evaluated are found to be in good agreement with in-house and existing test results of different composites. This study clearly demonstrates the effect of root rotation at the crack-tip of the DCB specimen on the critical fracture toughness evaluation. However, if the rotational spring stiffness is very large (i.e., ), the crack-tip behaves like built-in support and the critical load of (9) coincides with the results of simple cantilever beam theory (see Figures 10 to 15).

It is noted from Figures 4 to 9 and Figures 16 to 18 that cantilever beam theory overestimates the critical load and underestimates deflection of the test results. The multidirectional laminates (with +45°/+45° and +45°/−45° interfaces) show higher fracture energy than the unidirectional one.

Nomenclature

:Crack length
:Initial crack length
, :Constants in power-law
:Width of the DCB specimen
:Compliance of the specimen
DCB:Double cantilever beam
:Longitudinal tensile modulus
:Fracture toughness or critical strain energy release rate
:Strain energy release rate at any point
2:Specimen thickness
:Rotational spring stiffness at the crack-tip
:Number of fracture data
:Applied load on both sides of the specimen
:Crack mouth opening displacement.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.