Interface crack between two interface deformable piezoelectric layers

Based on an interface deformable piezoelectric bi-layer beam model, a bonded piezoelectric bi-material beam with an interface crack perpendicular to the poling axis is analyzed within the framework of the theory of linear piezoelectricity. The layer-wise approximations of both the elastic displacements and electric potential are employed, and each sub-layer is modeled as a single linearly elastic Timoshenko beam perfectly bonded together through a deformable interface. Using the impermeable crack assumption, the closed form solutions for the energy release rate (ERR) and crack energy density (CED) are derived for the layered piezoelectric beam subjected to combined uniformly distributed electromechanical loading. Based on superposition principle, both the ERR and CED and their components are all reduced to the functions of the crack tip loading parameters. Loading dependence of the total CED with respect to the applied electric field is manifested with the analytical results, showing that there is a transformation from an even dependence to an odd dependence for the normalized CED when the applied mechanical loading increases. Compared with the commonly used equivalent single layer model, the proposed analysis augments the crack driving force by alleviating the stress concentration along the interface and thus increases the loading parameters at the crack tip. The proposed model provides improved solutions for fracture analysis of piezoelectric layered structures and sheds light on the loading dependence of the fracture parameters (i.e., the ERR and CED) with respect to the applied electromechanical loadings.     

Delamination of layered structures on elastic foundation

This study presents an interface fracture mechanics analysis of delamination of a layered beam resting on a Winkler elastic foundation subject to general mechanical loads. A crack tip element on elastic foundation model is established first, through which, two concentrated forces existing at the crack tip are determined in closed-form. Then total energy release rate of the crack can be expressed in term of these two forces. By using available numerical results in the literature, the phase angle of the total energy release rate is also obtained. To verify the validity and accuracy of the solutions, debonding of a bonded overlay from the base structure resting on a Winkler elastic foundation is analyzed using the present solution. Comparisons with the baseline results obtained by finite element analysis suggest that the present analytical solution provides an excellent estimation of the total energy release rate and its phase angle for interface cracks in layered structure on elastic foundation. This study provides an approximated analysis of the debonding of a thin overlay debonding from the concrete pavement, where the effect of the base structure is simplified by a Winkler elastic foundation. This solution can also be used to analyze other similar delamination problems, such as local delamination in laminated composites, and face sheet delimitation in sandwich beams.     

Interface crack in sandwich specimen

Fracture of a sandwich specimen loaded with axial forces and bending moments is analyzed in the context of linear elastic fracture mechanics. A closed form expression for the energy release rate for interface cracking of a sandwich specimen with isotropic face sheets is found from analytical evaluation of the J-integral. An approach is applied, whereby the mode mixity for any combination of the loads can be calculated analytically when a load-independent phase angle has been determined. This load-independent phase angle is determined for a broad range of sandwich configurations of practical interest. The load-independent phase angle is determined using a novel finite element based method called the crack surface displacement extrapolation method. The expression for the energy release rate is based on the J-integral and certain stress distributions along the ends of the sandwich specimen. When the stresses from the crack tip interacts with the stresses at the ends, the present analytical calculation of the J-integral becomes inaccurate. The results show that for the analytically J-integral to be accurate the crack tip must be a certain distance away from the uncracked end of the specimen. For a sandwich specimen with face sheet/core stiffness ratio of 100, this distance is in the order 10 times the face sheet thickness. For sandwich structures with face sheet/core stiffness ratio of 1,000, the distance is 30 times the face sheet thickness.     

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.     

Delamination tolerance studies in laminated composite panels

Determination of levels of tolerance in delaminated composite panels is an important issue in composite structures technology. The primary intention is to analyse delaminated composite panels and estimate Strain Energy Release Rate (SERR) parameters at the delamination front to feed into acceptability criteria. Large deformation analysis is necessary to cater for excessive rotational deformations in the delaminated sublaminate. Modified Virtual Crack Closure Integral (MVCCI) is used to estimate all the three SERR components at the delamination front from the finite element output containing displacements, strains and stresses. The applied loading conditions are particularly critical and compressive loading on the panel could lead to buckling of the delaminated sublaminate and consequent growth of delamination. Numerical results are presented for circular delamination of various sizes and delamination at various interfaces (varying depth-wise location) between the base- and the sub-laminates. Numerical data are also presented on the effect of bi-axial loading and in particular on compressive loading in both directions. The results can be used to estimate delamination tolerance at various depths (or at various interfaces) in the laminate.     

In-plane elastic constants and strengths of circular cell honeycombs

The in-plane elastic modulus, Poisson’s ratio, brittle crushing strength and plastic yielding strength of honeycombs with hexagonally packed circular cells are analyzed theoretically. The resulting theoretical expressions are compared with the numerical results obtained from a series of finite element analyses for circular cell honeycombs with various relative densities, leading to a good correlation. It is also found that the in-plane mechanical properties of circular cell honeycombs are significantly affected by the ratio of cell-wall thickness to cell radius. Though the elastic constants along the two principal directions of circular cell honeycombs are the same, the brittle crushing strength and plastic yielding strength along the two principal directions are not identical. Furthermore, the in-plane mechanical properties of circular cell honeycombs are compared to those of regular hexagonal honeycombs with straight and uniform-thickness cell walls to evaluate their microstructural efficiency.     

Nonlinear response of a shear deformable S-FGM shallow spherical shell with ceramic-metal-ceramic layers resting on an elastic foundation in a thermal environment

This article presents an analytical approach to investigate the nonlinear stability of thick, functionally graded material (FGM) shallow spherical shells resting on elastic foundations, subjected to uniform external pressure and exposed to thermal environments. Material properties are assumed to be temperature dependent and graded in the thickness direction according to a Sigmoid power law distribution (S-FGM) in terms of the volume fractions of constituents. Using the first-order shear deformation theory and the Galerkin method, the effects of materials, geometry, elastic foundation parameters, and temperature on the nonlinear response of the thick S-FGM shells are analyzed and discussed in detail.     

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.     

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.     

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.     

The role of flaw geometry in thin film delamination from two-dimensional interface flaws along free edges

Delamination from planar interface edge flaws between a thin film and a semi-infinite substrate is examined to determine the roles of flaw width and depth relative to the film thickness. The flaws have curved and straight sections, and the crack front intersects the free edge at a right angle. Three-dimensional finite element models are used to extract local energy release rates and mode-mixity angles along the entire crack front. This paper focuses the crack dimensions required to reach steady state, wherein energy release rates are independent of flaw dimensions along the entire crack front. Results indicate that moderate elastic mismatch, although affecting mode mixity, plays a small role in determining the crack aspect ratios required to reach steady state. For wide cracks, the energy release rate for crack advance into the film interior approaches the plane-strain steady-state value when the half-width of the crack is approximately four times its depth (for cracks whose depths is several times the film thickness). For narrow cracks, the energy release rate near the free edge is significantly greater than the plane-strain steady-state result, and reaches a steady state when the depth approximately 10 times its width (for widths several time the film thickness). The results imply that delamination from wide cracks is reasonably accurately predicted via plane-strain analyses. Conversely, two-dimensional models are incapable of accurately predicting delamination from narrow cracks, which have a tendency to widen into flaws with more balanced aspect ratios (i.e. without growth in the depth direction).     

Delamination growth in Fibre Metal Laminates under variable amplitude loading

This paper presents the experimental and analytical investigation of the effect of variable amplitude (VA) load sequences on delamination behavior in Fibre Metal Laminates (FMLs). Delamination tests were performed and results are compared with linear damage accumulation predictions. Scanning Electronic Microscopy (SEM) was used to analyse the delaminated surfaces to study the delamination growth rate under VA loading in more detail. The correlation between test results and predictions highlighted the absence of load sequence and interaction effects in delamination growth rate under VA loading. This correlation is supported by the SEM observations.     

Effect of damage levels and curing stresses on delamination growth behaviour emanating from circular holes in laminated FRP composites

This paper deals with the study of the effect of drilling induced delamination damage levels and residual thermal stresses (developed during manufacturing process of cooling the laminate form curing temperature to room temperature) on delamination growth behaviour emanating form circular holes in graphite/epoxy laminated FRP composites. Two sets of full three dimensional finite element analyses (one with the residual thermal stresses developed while curing the laminate and the other without residual thermal stresses i.e. with mechanical loading only) have been performed to calculate the displacements and interlaminar stresses along the delaminated interfaces responsible for the delamination onset and propagation. Modified crack closure integral (MCCI) techniques based on the concepts of linear elastic fracture mechanics (LEFM) have been used to calculate the distribution of individual modes of strain energy release rates (SERR) to investigate the interlaminar delamination initiation and propagation characteristics. Asymmetric variations of SERR obtained along the delamination front are caused by the overlapping stress fields due to the coupling effect of thermal and mechanical loadings. It is found that parameters such as ply orientation, drilling induced damage levels and material heterogeneity at the delaminated interface dictate the interlaminar fracture behaviour of laminated FRP composites.     

Effects of substrate compliance on circular buckle delamination of thin films

The present work analyzes circular delamination buckling in a film/substrate system based on the Von Karman nonlinear plate theory with the consideration of elastic deformation of the substrate. Due to the axisymmetry of circular buckling, the substrate deformation is modeled by coupled springs and the spring compliances are determined from the dimension analysis and finite element calculations. The numerical shooting method is used to solve the nonlinear post-buckling problem. The stress intensity factors, the energy release rate, and the phase angle are given here for a variety of the elastic mismatch between the film and the substrate. The results show that in some cases, the energy release rate can be several times larger than that derived from the widely used clamped edge condition.     

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.     

The effect of facesheet/core delamination in sandwich structures under transverse loading

The present work is concerned with the problem of a delamination crack along the facesheet/core interface of a sandwich structure which is submitted to transverse loading. In contrast to a loading by compressive inplane forces or a bending loading the presumed transverse loading does not lead to buckling of the delaminated facesheet but it may provoke further delamination crack growth. As a kind of crack driving force the energy release rate is studied for a virtual crack growth by means of Irwin's crack closure integral within a finite element modelling. The resultant energy release rate is dependent on various geometrical and material parameters which is investigated in some detail.     

Thermal stress intensities at an interface crack between two elastic layers

The thermal stress intensities (energy release rate and stress intensity factors) due to temperature changes are derived in closed-form for an interface crack between two elastic layers of dissimilar materials. The solutions are two-dimensional and tabulated over a wide range of material and layer thickness combinations. The tables serve as rapid evaluations of the thermal stress intensities for given temperature changes. A strain gauge technique is given for determining constraint coefficients which reflect the constraint conditions during the temperature changes. The solutions are compared with results from the literature. The stress intensities due to thermal and mechanical loads are generally superimposed. As an example of application, the solutions are utilized to obtain the complete thermal and mechanical stress intensities for a four-point bend specimen.     

A size-dependent third-order shear deformable plate model incorporating strain gradient effects for mechanical analysis of functionally graded circular/annular microplates

In this paper, we develop a novel size-dependent plate model for the axisymmetric bending, buckling and free vibration analysis of functionally graded circular/annular microplates based on the strain gradient elasticity theory. The displacement field is chosen by using a refined third-order shear deformation theory which assumes that the in-plane and transverse displacements are partitioned into bending and shear components and satisfies the zero traction boundary conditions on the top and bottom surfaces of the microplate. Besides, the present model contains three material length scale parameters to capture the size effect. The material properties of the microplate are assumed to vary in the thickness direction and estimated through the classical rule of mixture. By using Hamilton's principle, the equations of motion and boundary conditions are obtained. Afterward, the differential quadrature method is adopted to discretise the governing differential equations along with various types of edge supports and therefore the deflection, critical buckling load and natural frequency can be determined. Convergence and comparison studies are carried out to establish the reliability and accuracy of the numerical results. Finally, a parametric study is conducted to investigate the influences of material length scale parameters, gradient index, thickness-to-outer radius ratio, outer-to-inner radius ratio and boundary conditions on the mechanical characteristics of the microplate.     

Mechanics of tunnelling cracks in trilayer elastic materials in tension

In this work, the crack driving force for a tunnelling crack in a thin brittle layer confined by dissimilar thick, and more compliant, elastic layers is considered at tensile loading. The steady-state energy release rate is evaluated using distributed dislocation technique and series representation of the complex potentials for an isotropic trimaterial. Evolution of the energy release rate with the crack length is studied by means of FEM. The 3D FEM simulations for tunnel cracks suggest that the ERR can represented by a universal relation (mastercurve) in suitably normalised co-ordinates. An analytical approximation of the ERR mastercurve is obtained as a function of crack length, cracking layer thickness, and a non-dimensional steady-state ERR.     

Influence of bending rotations on three and four-point bend end notched flexure tests

It has been experimentally observed that mode II critical energy release rate (G_IIC) values determined by four-point end notched flexure test and three-point end notched flexure test are different for the same material. At the present work correction factors related to bending rotations are introduced to explain the differences between values of G_IIC obtained by three point and four point end notched flexure tests. The bending angle leads to the contact zones between specimen and supports and specimen and load rollers changing in both test configurations. The present analysis has been carried out by the classical beam theory, neglecting shear effects and assuming the hypothesis of small rotated angles. Results show that the relative differences between corrected and uncorrected values of G_IIC are greater in the case of four-point end notched flexure than in the case of three-point end notched flexure test.