Mode-I Metal-Composite Interface Fracture Testing for Fibre Metal Laminates

-e main contribution of the present paper is the determination of the mode-I fracture of metal-composite interface region for fibre metal laminates (FMLs). A hybrid DCB configuration is proposed to investigate the mode-I fracture between metalcomposite interface using experimental and numerical approaches. A computationally efficient and reliable finite element model was developed to account for the influence of metal plasticity on the measured fracture energy.-e results of the experimental and numerical studies showed that metal plasticity increases the fracture energy of the metal-composite interface as the fracture event progresses. -e applied energy truly utilized to propagate metal-composite interface fracture was predicted numerically by extracting the elastic strain energy data. -e predicted true fracture energy was found to be approximately 50% smaller than the experimentally measured average propagation energy. -e study concluded that metal plasticity in hybrid DCB configuration overpredicted the experimentally measured fracture energy, and this can be alleviated through numerical methodology such as the finite element approach as presented in this paper.


Introduction
Hybridization of material not only makes it structurally efficient but also makes it multifunctional.With sustained research and advancements in adhesive bonding technology, designers are able to successfully develop a novel metalbonded hybrid composite structure such as the fibre metal laminates (FMLs) in the aviation industry.Currently, FML structures are the third largest primary load-carrying aviation material by volume housing nearly 470 sq•m of the upper fuselage section of Airbus A380 passenger air carrier [1].FML is a versatile and high performance hybrid material made up of alternating layers of thin metal sheets and thin fibre reinforce polymeric (FRP) composite layers.Compounding metal and composite makes FML stronger, lighter, better damage tolerant, good fatigue resistant, better corrosion, and impact resistant material which is not feasible through single material alone [2,3].e effectiveness of FML as a good load-carrying material relies heavily on the level of the bonding integrity between the metal and the FRP composite layers [4].e bonding interface must possess sufficient interfacial strength to be able to distribute the load between two characteristically different materials.ere are two different bonding interfaces in a traditional FML: (i) adhesive-metal and (ii) adhesive-composite.e level of bonding integrity is strictly dependent on the appropriate surface pretreatment of bonded metal to composite interfaces and the chemical compatibility of adhesive introduced.Apart from those, the potential bonding strength is also dependent on the stiffness imbalance and the thermal mismatch of dissimilar materials during testing and curing, respectively [5].For the case of FML, the problem of thermal mismatch is a challenging one where the difference in thermal expansion coefficients between adhesive, metal, and composite layers needs to be resolved.A study [6] had shown that poststretching after curing helps to enhance delamination resistance in the FML panels.
Adhesive bonding between monolithic material substrates like metal-metal and composite-composite is fairly well understood, and the testing procedures to evaluate their fracture properties are fairly standardized.However, this is not the case for evaluating dissimilar material joints like metal-composite interface.Furthermore, plastic yielding is predominately sensitive in the DCB sample having thin metal adherent and may attenuate the crack propagation in fracture test thereafter [7].Such plasticity can result in abnormally high values of fracture toughness, and it is necessary to account for its influence in order to obtain the true magnitude of the fracture energy.So far in literature, researchers had employed three different approaches to account for the plasticity effects: (i) experimental approach using thick metal doublers to attenuate plastic flow [6,8], (ii) theoretical approach [9], and (iii) finite element approach [7].In the last two approaches, efforts were made to estimate the fraction of plastically absorbed energy from the overall internal strain energy.
Reeder et al. [8] portrayed the method of bonding doublers to thin substrates in a manner to avoid plastic deformation or premature adherent failure due to bending before the interface delamination initiate.
e standard Irwin-Kies fracture energy equation was revised for the samples having bonded doublers for all three modes of fracture.Vlot and Van Ingen [6] utilized a similar concept to evaluate the delamination fracture energy of metal-composite interface of FML.To achieve stable delamination propagation in the mode-I test using the double cantilever beam (DCB) geometry, the transverse displacement of loading platen was stopped at every interval when there exists a delamination growth and stayed until the load magnitude leveled off.Lawcock et al. [10] followed the standard testing procedure of metallic joints to evaluate mode-I fracture energy of FML assuming that the volume fraction of the composite in the considered DCB geometry is negligibly small compared to the metal volume.
Reyes et al. [11] evaluated the mode-I, mode-II, and mixed-mode fracture energy of metal-composite interface of FML using a single cantilever beam (SCB), a double end notched flexure, and a modified mixed-mode flexure geometry, respectively.Using the SCB geometry, stable interface crack propagation without causing plastic yielding of metal sheet would become feasible.However, if the fracture energy of the adhesive material is far larger than the yielding stress of the metal sheet, this method could enforce more challenge to initiate interface crack growth.
e aim of the present paper is to investigate the mode-I fracture energy evaluation of metal-composite joints in order to understand the interface crack growth in FML that involves the influence of the plastic yielding in thin metal sheets.Finite element approach was utilized to estimate the amount of energy dissipated through plastic deformation.
e robustness of the developed numerical model was validated by comparing the fracture test results of monolithic composite-composite sample with those obtained from the numerical study.Such comparison also explores the influence of plasticity effects and allows one to quantitatively estimate how far the fracture energy magnitude has been altered due to such additional irreversible deformation.

Materials and Test Methodology
In this section, the mode-I fracture test geometry and the methodology used in the study are described.Subsequently, the investigation of mode-I fracture behavior of compositecomposite and metal-composite joints are presented, and some conclusive remarks are narrated in the final section.
2.1.Materials.Aluminum alloy Al 2024-0, L-530 8-harness satin weave 7781 glass epoxy prepreg, and Redux 335K adhesive film are the materials used to fabricate the adhesively bonded fibre-metal laminated samples investigated in this study.Before bonding, the aluminum alloy bonded surfaces were abraded using 320 grit size sand paper, cleaned with acetone solvent, and appropriately surface treated in a motive to promote good chemical linkage between aluminum-adhesive-composite interface.
Metal-composite layers of designated stacking were laminated and bonded with the single layer of adhesive film.Stacked laminates were cured in accordance with the composite manufacturer data sheet using an autoclave chamber under 3 bar and 120 °C in the presence of vacuum.
e cured laminates were cut to the required dimension without disturbing the adhesively bonded interface using a high-speed water jet cutter, and quality checks were carefully carried out under the optical microscope.

Experimental Tests.
In total, three different configurations of DCB tests were experimented to investigate the mode-I fracture behavior of metal-composite adhesive interface fracture growth.ey are as follows: (i) ick metal doublers bonded with thin composite adherent, M-C-M configuration (ii) Monolithic composite adherents with in situ epoxy matrix adhesive interface, C-C configuration (iii) Hybrid geometry with firmly bonded metalcomposite adherent as one arm and monolithic composite adherent as another arm, C-M-C configuration Schematic representation of above listed configurations and relevant dimensions are shown in Figures 1(a)-1(c), respectively.Samples were tested at a rate of 3 mm/min using a universal testing machine.Magnitude of load P, cross head displacement δ, and crack length a were recorded during both crack initiation and instant propagation.
e fracture energy for M-C-M configuration is evaluated using the following guidelines of ASTM D3833 assuming that the composite layer is very thin compared to thick aluminum doublers [10].
e equation of mode-I fracture energy, G I , is given by where E, B, and t 1 are Young's modulus, width, and thickness of the aluminum adherends.

Advances in Materials Science and Engineering
While for the C-C and C-M-C con gurations, the value of fracture energy is ascertained using the Irwin-Kies equation [12] given by where C δ/P is the specimen compliance and a is the corresponding crack length for a given displacement of δ.

Finite Element Simulation.
e simulation of interface crack initiation and propagation was performed using cohesive element in the commercial nite element platform, ABAQUS/Standard.e constitutive formulation of cohesive element is modeled based on the traction (t)separation (δ) law that relates the crack opening displacement in the process zone to the resisting tractions [13].e material model of cohesive element is de ned by an initial elastic sti ness, K, peak traction which causes crack initiation, t 0 α , and critical energy release rate equal to area under tractionseparation curve responsible to cause crack-propagation rate, G c α .e schematic of the bilinear traction-separation law used to simulate the progressive interface crack of DCB con guration in this paper is shown in Figure 2.
e input parameters required to model the cohesive elements are determined based on Turon's methodology [14].e methodology provides a systematic procedure to calculate the length of the cohesive zone, l cz , size of the cohesive element, l e , and number of cohesive elements, N e , required to model the fracture growth accurately.e methodology also allows us to relax the requirement of extremely ne meshes through arti cially increasing the length of cohesive zone by reducing the peak traction, t 0 α , magnitude.e mathematical derivations to obtain adjusted traction strength, t a α , and the optimum mesh parameters using Turon's methodology are summarized in the Appendix.Both metal and composite adherents are modeled using isotropic and orthotropic elastic material properties, respectively.
ese layers are discretized using an incompatible-mode eight-node brick element, C3D8I, and the interface region using three-dimensional cohesive element, COH3D8.
e material properties used in the nite element model for the substrates and interface region are provided in Table 1.
e FE model was validated by comparing the predicted load-displacement results with the experimental results.Using this curve, fracture energy data are calculated by employing the area method given by [15] where ΔU is the area under numerical P-δ curve between consecutive crack growth intervals and Δa is the corresponding crack length increment as shown in Figure 3.

Selection of DCB Con guration.
e initial emphasis to determine the metal-composite fracture behavior was analyzed through the experimental study of the M-C-M conguration.Enhanced resistance to bending because of thick aluminum adherends resulting in multiple interfacial and interlaminar crack propagation was observed between M-C and C-C regions, respectively, as shown in Figure 4(a).is behavior indicates that the presence of thick aluminum adherends desensitized these multiple surface nucleation.Only the initiation fracture energy is feasible from this experiment using (1) which amounted to about 0.656 N/mm.e DCB geometry with dissimilar adherends having bimaterial metal-composite adherent as one arm and monolithic composite as another arm (C-M-C con guration) was found to be a feasible design that allows the  Advances in Materials Science and Engineering evaluation of the fracture energy not only in the crack initiation phase but also throughout crack propagation path. in adherends admit the arm to bend freely and facilitate the preimitated crack front propagation along metalcomposite interface in stable manner even though there exist a significant amount of plastic deformation in the thin metal layer.
For the case of the hybrid configuration of DCB arms, the bending centroids of the two arms are generally different, and this difference results in unsymmetric flexure and inherently induces the mixed-mode loading state on the precrack front instead of pure mode-I loading state.Such mixed-mode inference was attenuated by selecting appropriate thickness for two DCB arms in a manner to have equal bending rigidity, EI, for both adherent arms.e derived geometric dimensions of the C-M-C configured DCB design is given in Figure 1(c).

DCB Results of C-C Configuration.
e DCB test using the C-C configuration is performed to obtain an accurate fracture energy of adhesive interface where the adherents undergo only elastic deformation.
is study primarily serves as a benchmark that allows one to understand the metal sheet plastic deformation influence on G I magnitude in the C-M-C configuration test.Moreover, it also provides opportunity to check the reliability of developed finite element model by ensuring that the internal energy utilized to delete the cohesive interface elements is the same as that of the experimentally predicted G I .
e experimental and predicted mode-I fracture crack interface between composite adherends is shown in Figures 5(a) and 5(b), respectively.No visible fibre bridging patterns were found along the cracked interface.is is because the weave pattern of the woven composite adherends completely hindered the fibre pullout failure.Figure 6 illustrates the comparison plots of the load P versus deflection δ and the fracture energy G I versus crack extension Δa. e smooth experimental load variation with no discrete load drops and change in fracture energy attributed in propagation phase confirms that the crack growth surfaces are pure cohesive in nature as shown in Figure 5(a).e average experimental fracture energy to cause crack initiation is found to be about 0.2 N/mm, and it increased to about 0.65 N/mm and maintained constant throughout the crack propagation phase.e latter propagation phase G I is used as a damage evolution material data input for cohesive elements in finite element simulation.

Selection of Optimum Mesh Parameters.
e optimum mesh parameters required to simulate DCB mode-I fracture were found using the equations presented in the Appendix.
e minimum cohesive zone length required to model mode-I fracture based on the experimental fracture energy for the C-C configuration is estimated about 0.29 mm using (A.1), and the material properties used are listed in Table 1.Based on that magnitude, parametric studies for three levels of mesh l e ( ≥ l cz ) were conducted for four different N e in the cohesive zone length.
e corresponding magnitude of adjusted traction strengths and artificially increased cohesive zone length for crack propagation was found using (A.3) and (A.1), respectively, and these are tabulated in Table 2. e predicted load-deflection response for different levels of mesh refinement is shown in Figure 7. Irrespective of l e , the response predicted by the adjusted interfacial strength t a α with N e � 1 did not converge and no numerical solution exists.For all other cases, the predicted softening crack propagation response is almost similar compared to the experimental measurement with some minor variations.
For the mesh size of l e � 1 mm, the softening response shows spurious oscillations for smaller N e which means the mesh size is too coarse to accurately predict the crack propagation.Take note that the amplitude of the spurious oscillation is found to be attenuated with increase in N e .It can be seen that the responses obtained with the finer mesh size of l e � 0.5 mm and l e � 0.3 mm overpredicted the peak load.Meanwhile for this mesh size, the softening response exhibited smooth crack propagation.On comparing the number of iterations required to complete the solution, mesh sizes with more N e in its artificially increased cohesive zone length l cz have shown to be computationally efficient.
e above comparative study confirms that the interface strength t 0 α required to initiate crack propagation does not have a significant influence on the overall response.us, the interface strength modification and the artificial cohesive zone length increment strategy initially proposed by Turon [14] can be used to model accurate crack propagation using cohesive elements at coarse mesh level.For the present DCB simulation, results obtained using lowest adjusted interfacial strength t a α with larger N e (l e �1 mm and N e � 8) yield computationally efficient solution.All finite element simulations presented in the current paper are based on the abovementioned optimum mesh size.
It can be seen in Figure 6 that the numerical solution obtained for the C-C configuration model having artificially increased cohesive length and adjusted traction strength shows good agreement with the experimental results.Cohesive element deletion replicating mode-I crack growth happened approximately around the stated fracture energy material input of 0.65 N/mm.

DCB Results of C-M-C Configuration.
Having validated the reliability of the developed finite element model to simulate mode-I crack growth, the fracture of adhesive interface between plastically deformed metal adherent and composite adherent was investigated using numerical approach.Experimental observation confirmed that no stiffness imbalance exists between the two dissimilar DCB arms, and complete crack growth is caused strictly by the opening mode of fracture.Figure 8 shows the load-displacement plots and the final permanently deformed shape of the fractured C-M-C configuration sample obtained from the experiments and the finite element model.With an increase in crack opening displacement, the nature of crack growth in experiments was found to change from stable to unstable manner.e continuous cumulative effect of plastic flow in metal layers makes the crack growth in a stick-slip manner in 6 Advances in Materials Science and Engineering the form of periodic short bursts and crack arrest.is scenario leaves a combined cohesive and adhesive kind of fracture surface along the bonded adherent surfaces.Figure 9(a) shows the variation of the measured fracture energy G I using (2) against crack extension Δa; it can be seen that the plasticity e ect causes G I to increase with increase in crack length.It seems that the peak propagation fracture energy increases nearly by 3.5 times more than the initiation fracture energy (i.e., G I at Δa 0). e absence of a single value of G I from the experimental study makes the selection of an appropriate G I input data required for numerical simulation a very di cult one.e fact that the magnitude of the fracture energy of the adhesives is constant and independent of the adherent material and geometry when their failure surface along the adherent interface is completely cohesive [16].Assuming that the simulated fracture surface of the C-M-C model is presumably cohesive, a propagation G I magnitude of  Advances in Materials Science and Engineering 0.65 N/mm obtained in the monolithic C-C DCB sample is used instead.Two di erent numerical cases were simulated with the distinct material model de nition for thin metal sheets: (i) pure elastic (El) and (ii) combined elastic-plastic (El-Pl) to investigate the in uence of plastic e ects.e material data required to model plastic response of metal layer were obtained from the standard tensile test as given in [17].
From Figure 8(a), it appears that the predicted initial elastic responses of both cases are correlated well with the experimental load-displacement curve.After reaching the plastic state, only the case with metal plasticity inclusion coincides with the experimental crack propagation phase until the stable crack region.Numerically, the plasticity e ects were found to increase the propagation fracture energy nearly by an average of 1.5 times compared to the pure elastic model as shown in Figure 9(a).e true fracture energy value was calculated by extracting the elastic strain energy magnitude from the fractured elastic-plastic numerical model and substituted into (4) which is similar to (3): where U R is the recoverable (or elastic) strain energy obtained using ABAQUS output variable ALLSE in postprocessing.By plotting U R versus a as shown in Figure 9(b) and taking the derivative, the average true propagation fracture energy (G I ) Tr can be found using ( 4).
e values of the fracture energy obtained from the di erent cases are summarized in Table 3.
e obtained (G I ) Tr values are approximately equal to the initiation fracture energy (G I ) In predicted from experiment while it is 50% smaller than the average propagation fracture energy (G I ) Pr .Such huge di erences are ultimately caused by the in uence of the plasticity e ect in the thin metal layer and its induced unstable nature of crack propagation.

8
Advances in Materials Science and Engineering bre metal laminates (FMLs) using experimental and nite element study.A robust numerical approach was proposed to calculate the true fracture energy of the interface crack propagation between the plastically deformed metal and elastic composite surfaces.Di erent hybrid DCB congurations based on the type of adherends were used in the experimental study to observe metal-composite interface fracture not only in initiation phase but also in the complete propagation phase.e M-C-M con guration containing thick metal adherent resulted in highly unstable crack growth via nucleating multiple interfacial and interlaminar cracks, and no reliable fracture data can be extracted.

Advances in Materials Science and Engineering
An accurate fracture energy G I of the adhesive used was found using monolithic woven composite C-C configuration.e crack propagation surfaces observed in the test were predominantly cohesive in nature, and no sign of fibre bridging was apparent.e obtained G I was used as a material input data in the developed finite element model to simulate damage evolution of interface cohesive element.Good correlation was found with the developed finite element model, and this ensured the reliability in modeling mode-I fracture in the interface of dissimilar layers.
Finally, the influence of the plastic deformation on G I was studied employing the C-M-C configuration where one arm of DCB adherent was designed to be metal-composite hybrid and the other with only composite material.is experimental study indicated that the plastic flow in thin metal layer continuously increases the measured G I by up to 3.5 times as the metal-composite interface fracture progresses.Similar response was also predicted in finite element study where the inclusion of the plastic constitutive model for metal layer increases the average propagation energy by 1.5 times compared to the pure elastic numerical model.A true fracture energy (G I ) Tr was obtained using the numerically predicted internal elastic energy in the finite element model.e predicted result indicated that (G I ) Tr was nearly equal to experimentally measured initiation fracture energy (G I ) In and 50% smaller than the average propagation fracture energy (G I ) Pr . is difference indicated the influence of metal plasticity effects on the measured G I magnitude through experiments.

Figure 2 :
Figure 2: Bilinear traction-separation law for cohesive elements in the DCB model.

Figure 7 :
Figure 7: Load-displacement curve of C-C con guration for di erent cohesive element sizes: (a) l e 1 mm, (b) l e 0.5 mm, and (c) l e 0.3 mm.

Figure 8 :
Figure 8: Mode-I fracture of the C-M-C con guration: (a) load-displacement curve and (b) permanent plastic deformation.

Figure 9 :
Figure 9: Mode-I fracture of C-M-C con guration: (a) fracture energy-crack extension and (b) recoverable elastic strain energy-crack length variation.

Table 1 :
Material properties used in the FE model.
0 Figure 3: Area method of fracture energy determination.

Table 2 :
Adjusted normal interfacial strength with respect to di erent mesh sizes.