Application of the energy release rate approach for delamination growth in Glare

Fatigue crack growth in a fibre metal laminate such as Glare is accompanied by delamination growth at the interface between the aluminium and glass fibre/adhesive layers. To incorporate this delamination growth in crack growth prediction methods, the energy release rate approach is applied to describe the delamination growth rate. Tests were performed to determine the relationship between the delamination growth rate and the calculated energy release rate.     

基于虚拟裂纹闭合技术的应变能释放率分析

基于虚拟裂纹闭合技术(VCCT),建立了复合材料层合板层间裂纹尖端的应变能释放率(SERR)三维有限元计算模型。该模型考虑了裂纹尖端大转动和离散单元形状变化对应变能释放率计算的影响,修正了裂纹尖端应变能释放率的计算方法。利用该模型计算了裂纹长度为15 mm和35 mm时纯Ⅰ型和纯Ⅱ型的应变能释放率,纯Ⅰ型应变能释放率分别为 207 J/m2和 253 J/m2;纯Ⅱ型应变能释放率分别为 758 J / m 2和 1040 J / m2;计算值与试验值吻合得很好。同时,该模型计算了混合型不同比值 R=(G_Ⅱ/G_Ⅰ+G_Ⅱ)的长裂纹层合板层间断裂过程的应变能释放率,其中Ⅰ型和Ⅱ型应变能释放率计算值与试验平均值的最大误差为 11.4%,最小误差为 0.4%。该模型能有效计算裂纹尖端的应变能释放率。     

Prediction of energy release rates for crack growth using FEM-based energy derivative technique

This paper presents a new method, named energy derivative technique, to calculate energy release rate for a variety of crack growth scenarios. The new method is based on energy conservation principle for crack growth, and is applicable to crack development in any quasi-static condition in which dynamic energy for crack growth is negligible. The method has the advantage over existing finite element-based methods in that the former does not require an elaborate fine mesh in the vicinity of a crack tip, and is not limited to linear deformation behaviour. Several case studies are presented to demonstrate validity of the method, which are (i) growth of penny-shaped crack for linear elastic fracture behaviour, (ii) crack growth in rubber sheet under tension for nonlinear elastic fracture behaviour, (iii) delamination in end-notched flexure specimen with friction, and (iv) crack growth with plastic deformation in double-edge-notched plate under tension. Results from these case studies show excellent agreement with data available in the literature, which were determined using either analytical or other FEM-based techniques.     

Off-axis fatigue crack growth and the associated energy release rate in composite laminates

Stable matrix crack growth behaviour under mechanical fatigue loading has been studied in a quasi-isotropic (0/90/-45/+45)_s GFRP laminate. Detailed experimental observations were made on the accumulation of cracks and on the growth of individual cracks in +45° as well as 90° plies. A generalised plain strain finite element model of the damaged laminate has been constructed. This model has been used to relate the energy release rate of growing cracks to the crack growth rate via a Paris relation.     

A new conservation integral for circular arc crack under multiple loads

In this paper, a path independent integral that represents the rate of energy flux at the tip during crack extension in a homogeneous and isotropic material has been derived from the principle of virtual work for a two-dimensional stationary circular arc crack subjected to multiple loads. This integral is an extension of the two-dimensional version of F-integral and includes the presence of the effects of thermal strains, initial strains and body forces, hitherto, unavailable in open literature, to the best of the authors’ knowledge. It has been further demonstrated that Rice’s J-integral, the -integral derived by Kishimoto et al. and the F-integral proposed by Lorentzon et al. are special cases of the generalized integral . The integral has been implemented into a finite element post-processing program for examining the path independence behavior under elastic and elastic-plastic deformation subjected to mechanical loads and thermo-elastic analyses under pure thermal loads. Within the limits of numerical accuracy, the application demonstrates that the solutions for the energy release rate on different contours preserve nearly identical values over the computational range.     

Crack-tip force method for computing energy release rate in delaminated plates

A new method called the crack-tip force method (CTFM) is derived for computing the energy release rate in delaminated beams and plates. In this method the delaminated plate is divided into two laminates on either side of the plane of delamination. The interaction forces, called crack-tip forces, between the sub-laminates at the crack-tip are computed. The energy release rate is expressed as a quadratic function of the crack-tip forces and the plate compliance coefficients. The CTFM is compared to the virtual crack closure technique (VCCT) as well as to a previously derived method called the strain energy density method using double cantilevered beam specimens as examples. The CTFM is found to be very efficient as the crack-tip forces are part of the solution of finite element analysis of delaminated plates, and they can be readily used to compute the point-wise energy release rate along the delamination front.     

Accuracy assessment of a three-dimensional,crack tip element based approach for predicting delamination growth in stiffened-skin geometries

Experimental observations of delamination growth in two stiffened-skin geometries are compared to predictions made using a three-dimensional crack tip element based approach. Each geometry consists of a six-ply graphite/epoxy skin co-cured to a six-ply, hat-shaped stiffener containing a preimplanted teflon delamination between the skin and stiffener at the stiffener termination point. One stiffened-skin geometry was loaded in three-point bending and the other had in-plane tension loads applied to the skin. To predict delamination growth, a three-dimensional crack tip element analysis was first performed on each geometry in order to determine the total energy release rate, G, as well as its mode I, II and III components, G_I, G_II and G_III, respectively. These results were used to define a mode mix at each point along the delamination front, G_s/G, where G_s=G_II + G_III. To obtain the delamination toughness, G_c, it was assumed that G_c exhibits the same dependence on G_s/G as on G_II/G, where the results for G_c versus G_II/G were taken from an earlier experimental study. Next, a comparison of the energy release rate to the toughness at each position along the delamination front was performed, and these results were scaled appropriately in order to predict the sequence of loads and corresponding locations at which the delamination will advance. The predictions were then compared to experimental results that included c-scan images of the test specimens taken at each increment of observed growth, and very good quantitative and qualitative correlations were obtained for both geometries. These results indicate the practicality of, and considerable computational savings that may be achieved by, employing crack tip element analyses for delamination growth predictions in realistic structural geometries.     

Computations of energy release rate under monotonic and cyclic loading conditions

In case of an elastic–plastic fracture mechanics analysis, the determination of the energy release rate distribution is a crucial point. In the present paper, three numerical techniques: the virtual crack closure technique (VCCT), J-integral and energy derivative technique (EDT), are used to compute the energy release rate in a middle-crack tension specimen with the combined isotropic/kinematic hardening model. The results obtained by these methods are compared with each other under monotonic and cyclic loading conditions. Finally, it comes out that the difference of the VCCT method to the J-integral is rather insensitive to load increasing, especially when the traction >40% of yield stress, however, the deviation of the VCCT and J-integral results are within 10%, suggesting that one may use the VCCT for plastic cracked specimen analysis. The computations show that the EDT provides the same values for the monotonic as the J-integral if the plastic deformations are not large, but for high plastic loading the EDT overestimates the fracture energy. For cyclic loading case, VCCT method offers closer results as the elastic analytical results, also suggesting that the whole plastic dissipated energy in the loading process should be integrated. While EDT method gives the smaller results than the J-integral because of the energy dissipated in the unloading phase is considered in the loading process.     

The surface effect on the strain energy release rate of buckling delamination in thin film-substrate systems

Gurtin-Murdoch continuum surface elasticity model is employed to study the buckling delamination of ultra thin film-substrate system. The effects of surface deformation and residual stress on the large deflection of ultra thin film are considered in analysis. A concept of effective bending rigidity (EBR) for ultra thin plate is proposed on the basis of Gurtin-Murdoch continuum theory and the principle of minimum potential energy. The governing equations with EBR are formally consistent with the classical plate theory, including both small deflection and large deflection. A surface effect factor is introduced to decide whether there is need to consider the surface effect or not. Combining the buckling theory and interface fracture mechanics, we obtain analytical solutions of the critical buckling load and the energy release rate of the interface crack in the film-substrate system. It is seen that the surface deformation and residual stress have significant effects on the buckling delamination of ultra thin film-substrate system.     

Strain energy release rate calculation for a moving delamination front of arbitrary shape based on the virtual crack closure technique. Part I: Formulation and validation

This paper proposes a simple, efficient algorithm to trace a moving delamination front with an arbitrary and changing shape so that delamination growth can be analyzed by using stationary meshes. Based on the algorithm, a delamination front can be defined by two vectors that pass through any point on the front. The normal vector and the tangent vector for the local coordinate system can then be obtained based on the two delamination front vectors. An important feature of this approach is that it does not require the use of meshes that are orthogonal to the delaminations front. Therefore, the approach avoids adaptive re-meshing techniques that may create a large computational burden in delamination growth analysis. An interface element that can trace the instantaneous delamination front, determine the local coordinate system, approximate strain energy release rate components and apply fracture mechanics criteria has been developed and implemented into ABAQUS^® with its user-defined element (UEL) feature. In this Part I of a two-part paper, the approach and its implementation are described and validated by comparison to results from existing cases having analytical solutions or other established FEA predictions.     

Fracture mechanics approach to facesheet delamination in honeycomb: measurement of energy release rate of the adhesive bond

Two types of experiments were designed and performed to evaluate the adhesive bond in honeycomb sandwich panels. The tensile bond strength between the facesheet and the core was determined through the flatwise tension test. The fracture toughness of the bond line was measured through the double cantilever beam test. Fracture toughness values varied for different facesheet thicknesses and core materials. Toughness was also different for the bag and tool sides of the panels for all specimen types.     

Energy release rate and bifurcation angles of piezoelectric materials with antiplane moving crack

The dynamic fracture problems of the piezoelectric materials with antiplane moving crack are analysed by using function of complex variable in the paper. The results show that the coupled elastic and electric fields inside piezoelectric media depend on the speed of the crack propagation, and have singularity at the crack tip. The stress intensity factor is independent of the speed of the crack propagation, which is identical to the conclusion of purely elasticity. Moreover, independent of the electric loading, the dynamic energy release rate can be expressed by the stress intensity factor and enlarge with the increase of crack speed. High speed of the crack moving could impede the crack growth. At the same time, the crack can be propagated into either curve or bifurcation if the crack speed is higher than the critical speed.     

Strain energy release rate calculation for a moving delamination front of arbitrary shape based on the virtual crack closure technique. Part II: Sensitivity study on modeling details

The interface element and VCCT process described in Part I of this two-part paper, developed to compute strain energy release rates of an arbitrary delamination front using non-orthogonal finite element meshes, are further investigated in this paper for robustness and ease of use in tracking delamination growth. Standard 3-D elements are used in conjunction with the interface elements. No special singularity elements are required. Stationary meshes that are independent of the shape of the delamination front can be used. Three cases having different initial delamination shapes are examined. The process is shown to be insensitive to the values used for the interfacial spring stiffness, the orientation of the interface element, or even the mesh pattern if the mesh has a reasonable degree of refinement. Therefore, the method can be used with ease and confidence in general-purpose delamination growth analysis for engineering applications.     

Numerical evaluation of the effects of friction and geometric nonlinearities on the energy release rate in three- and four-point bend end-notched flexure tests

Nonlinear finite element analyses are used to examine the effects of friction and geometric nonlinearities on the energy release rate in three- and four-point bend end-notched flexure tests. Energy release rates are first determined by a recently developed direct energy balance approach. It is shown that the finite diameter loading rollers that are typically used in practical test set-ups cause both tests to be inherently nonlinear. The effect of these nonlinearities on the energy release rate is shown to be larger in the four point than the three point test and to increase with increasing roller diameter, increasing coefficient of friction along the crack plane, and decreasing supporting span length. For the four point test, the effect of these nonlinearities is also shown to increase with increasing ratio of inner to outer span length. Next, energy release rates at the onset of crack advance are determined by a simulated compliance calibration technique. This “perceived toughness” is compared with predictions of the “true toughness” given by the direct energy balance approach at the same load. It is shown that perceived toughnesses from this simulated compliance calibration procedure are larger than previously reported results that were obtained in a similar fashion using linear theory. In addition, the perceived toughness is shown to strongly depend upon the range used for fitting the load versus deflection data to obtain compliance. These findings are used to make some general recommendations regarding use of the two test methods and their associated data reduction techniques.     

Calculation of transient strain energy release rates under impact loading based on the virtual crack closure technique

This paper describes an interface element to calculate the strain energy release rates based on the virtual crack closure technique (VCCT) in conjunction with finite element analysis (FEA). A very stiff spring is placed between the node pair at the crack tip to calculate the nodal forces. Dummy nodes are introduced to extract information for displacement openings behind the crack tip and the virtual crack jump ahead of the crack tip. This interface element leads to a direct calculation of the strain energy release rate (both components G_I and G_II) within a finite element analysis without extra post-processing. Several examples of stationary cracks under impact loading were examined. Dynamic stress intensity factors were converted from the calculated transient strain energy release rate for comparison with the available solutions by the others from numerical and experimental methods. The accuracy of the element is validated by the excellent agreement with these solutions. No convergence difficulty has been encountered for all the cases studied. Neither special singular elements nor the collapsed element technique is used at the crack tip. Therefore, the fracture interface element for VCCT is shown to be simple, efficient and robust in analyzing crack response to the dynamic loading. This element has been implemented into commercial FEA software ABAQUS^® with the user defined element (UEL) and should be very useful in performing fracture analysis at a structural level by engineers using ABAQUS^®.     

Electroelastic analysis of an interface anti-plane shear crack in a layered piezoelectric plate

The problem of an anti-plane interface crack in a layered piezoelectric plate composed of two bonded dissimilar piezoelectric ceramic layers subjected to applied voltage is considered. It is assumed that the crack is either impermeable or permeable. An integral transform technique is employed to reduce the problem considered to dual integral equations, then to a Fredholm integral equation by introducing an auxiliary function. Field intensity factors and energy release rate are obtained in explicit form in terms of the auxiliary function. In particular, by solving analytically a resulting singular integral equation, they are determined explicitly in terms of given electromechanical loadings for the case of two bonded layers of equal thickness. Some numerical results are presented graphically to show the influence of the geometric parameters on the field intensity factors and the energy release rate.     

Effect of electric boundary conditions on crack energy density and its derivatives for piezoelectric material

This paper deals with crack energy density (hereafter CED) as a possible fracture parameter in piezoelectricity under arbitrary electric boundary conditions on a notch-like crack surface. The definitions of CED and its derivatives are given first under exact boundary condition. Next, their path independent integrals are also derived and their approximate expressions are discussed under some restrictions on the crack surfaces. It is found that electrical terms along the notch-like crack surface do not vanish unlike in the case of impermeable crack. Then, we introduce evaluation methods of CED, and, with the help of the results of finite element analyses (FEA), we closely examine how electric boundary conditions along the notch surface and initial notch width influence CED and its derivatives. It is shown from the FEA results that because of the difficulties of computing path integral terms along the notch-like crack tip in the path independent expressions, the evaluation by the definitions of CED and its derivatives is preferable and more convenient than the evaluation of their path independent expressions. It is also found that all the parameters are significantly affected by both permittivity inside the electric inclusion and root radius of the notch. Finally, the possibility of mechanical CED as a governing fracture parameter is discussed.     

Dynamic mixed-mode I/II delamination fracture and energy release rate of unidirectional graphite/epoxy composites

Mixed-mode open-notch flexure (MONF), anti-symmetric loaded end-notched flexure (MENF) and center-notched flexure (MCNF) specimens were used to investigate dynamic mixed I/II mode delamination fracture using a fracturing split Hopkinson pressure bar (F-SHPB). An expression for dynamic energy release rate G_d is formulated and evaluated. The experimental results show that dynamic delamination increases linearly with mode mixing. At low input energy E_i ? 4.0 J, the dynamic (G_d) and total (G_T) energy rates are independent of mixed-mode ratio. At higher impact energy of 4.0 ? E_i ? 9.3 J, G_d decreases slowly with mixed I/II mode ratio while G_T is observed to increase more rapidly. In general, G_d increases more rapidly with increasing delamination than with increasing energy absorbed. The results show that for the impact energy of 9.3 J before fragmentation of the plate, the effect of kinetic energy is not significant and should be neglected. For the same energy-absorption level, the delamination is greatest at low mixed-mode ratios corresponding to highest Mode II contribution. The results of energy release rates from MONF were compared with mixed-mode bending (MMB) formulation and show some agreement in Mode II but differences in prediction for Mode I. Hackle (Mode II) features on SEM photographs decrease as the impact energy is increased but increase as the Mode I/II ratio decreases. For the same loading conditions, more pure Mode II features are generated on the MCNF specimen fractured surfaces than the MENF and MONF specimens.     

Fracture analysis of a surface through-thickness crack in PZT thin film under a continuous laser irradiation

In this paper, the surface crack problem in PZT thin films under a continuous laser irradiation has been investigated by the superposition principle. Using commercial (FEM) software ANSYS 9.0, the piezoelectric fields near the crack tip were solved for surface crack in the finite PZT thin film. The SIFs for crack-tip fields were obtained by using the limited stress extrapolation technique (LSET) and then the energy release rates were calculated by the relation to the intensity factors. When the irradiation time and crack location were changed, the energy release rates G, G_I, and G_e for total, mechanical terms (mode I) and electric contribution were investigated. The results show that the mechanical opening mode I is the main mode for the surface crack under a continuous laser irradiation. However, electric mode IV has inhibiting effect on crack growth. At the beginning of laser irradiation, the surface tiny crack which is close to the centre of film will propagate more easily. During the laser irradiation, the crack which is far from the centre of film will propagate more easily.