Next Article in Journal
Surface Segregation of Amphiphilic PDMS-Based Films Containing Terpolymers with Siloxane, Fluorinated and Ethoxylated Side Chains
Next Article in Special Issue
Ripple Formation during Oblique Angle Etching
Previous Article in Journal
Synthesis of Solution-Stable PEDOT-Coated Sulfonated Polystyrene Copolymer PEDOT:P(SS-co-St) Particles for All-Organic NIR-Shielding Films
Previous Article in Special Issue
Stress-Affected Lithiation Reactions in Elasto-Viscoplastic Si Particles with Hyperelastic Polymer Coatings: A Nonlinear Chemo-Mechanical Finite-Element Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fracture Mechanics Solutions for Interfacial Cracks between Compressible Thin Layers and Substrates

by
Roberta Massabò
1,*,
Konstantin Ustinov
2,
Luca Barbieri
1 and
Christian Berggreen
3
1
Department of Civil, Chemical and Environmental Engineering, University of Genoa, 16145 Genova, Italy
2
Institute for Problems in Mechanics, Russian Academy of Sciences, 119526 Moscow, Russia
3
Department of Mechanical Engineering, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
*
Author to whom correspondence should be addressed.
Coatings 2019, 9(3), 152; https://doi.org/10.3390/coatings9030152
Submission received: 30 January 2019 / Revised: 17 February 2019 / Accepted: 21 February 2019 / Published: 26 February 2019
(This article belongs to the Special Issue Modelling and Simulation of Coating 2019)

Abstract

:
The decohesion of coatings, thin films, or layers used to protect or strengthen technological and structural components causes the loss of their functions. In this paper, analytical, computational, and semi-analytical 2D solutions are derived for the energy release rate and mode-mixity phase angle of an edge-delamination crack between a thin layer and an infinitely deep substrate. The thin layer is subjected to general edge loading: axial and shear forces and bending moment. The solutions are presented in terms of elementary crack tip loads and apply to a wide range of material combinations, with a large mismatch of the elastic constants (isotropic materials with Dundurs’ parameters 1 α 1 and 0.4 β 0.4 ). Results show that for stiff layers over soft substrates ( α 1 ), the effects of material compressibility are weak, and the assumption of substrate incompressibility is accurate; for other combinations, including soft layers over stiff substrates ( α 1 ), the effects may be relevant and problem specific. The solutions are applicable to edge- and buckling-delamination of thin layers bonded to thick substrates, to mixed-mode fracture characterization test methods, and as benchmark cases.

1. Introduction

The decohesion of thin layers from thicker substrates is a well-known failure mechanism which may occur at very different scales in structural and technological components due to various loading conditions, and is controlled by the interfacial fracture toughness, a material property which depends on the relative magnitude of shear and normal tractions on the interface [1,2,3,4]. Delamination (or debonding) of the face sheets from the core is one of the dominant failure mechanisms of thick sandwich structures used in naval, aeronautical, automotive, and wind energy applications, and can be due to direct loadings, stress waves, or internal pressure cycles [4,5,6,7,8]. In these systems, the face sheets are typically laminate composite materials, and the core is made up of cellular foams, honeycomb or balsa; exemplary applications are naval bulkheads and deck panels, primary aerostructures, e.g., wing leading edges and rudders, and wind turbine blades. The debonding of Fiber Reinforced Polymer (FRP) sheets used to repair or strengthen reinforced concrete members in civil engineering applications prevents load transfer and nullifies the reinforcing action [9,10]. In thermal barrier coatings designed to reduce the temperature of the underlying metallic substrates for aircraft and power-generating turbines, delamination and spalling are common failure modes [11,12,13,14]. The cracking of thin films used in electronic packages or of metal alloy coatings on elastomers, used in soft electronics, is often driven by residual stresses [2,3,15,16,17,18]. The cracking and spalling of brittle layers may be a consequence of edge loadings [19]. The decohesion of polymer coatings used to protect metals for different applications, e.g., in beverage or food cans, biomedical and mechanical devices, or with anti-fouling functionality, leads to the loss of the protective and aesthetic properties of the coatings [20,21]. Many other thin film applications depend on the adhesion of the film to the substrate [22,23,24].
This paper presents fracture mechanics solution methods and results for the fundamental 2D problem in Figure 1, which describes an edge delamination crack between a thin layer and a semi-infinite substrate. The thin layer is subjected to general edge loadings and mode mixity phase angle and energy release rate associated with the collinear propagation of the crackare obtained in terms of crack tip force and moment resultants, Figure 1b. Isotropic bi-material interfaces with a general mismatch of the elastic constants are examined for a wide range of the Dundurs parameters α and β , which measure the mismatch in uniaxial and bulk compliances of the materials (see Section 2 for definitions). Results are presented for 1 α 1 and 0.4 β 0.4 in order to describe also stiff layers on soft substrates ( α 1 ) and soft layers on stiff substrates ( α 1 ), and to account for the mismatch in the volumetric stiffnesses of layer and substrate ( β 0 ). The work builds up and expands the works in [25,26,27,28] to provide analytical solutions for axial, bending, and shear loading with β = 0 , an efficient and accurate numerical method for the derivation of the fracture mechanics parameters in the numerically challenging film/substrate problem, and numerical coefficients for the application of the semi-analytic solutions in [25,26] to bi-materials with a general mismatch of the elastic constants and β 0 .
For its technological importance, the problem in Figure 1 has been studied and used extensively in the literature since the early work by Obreimoff on the splitting strength of mica [1,2,3,29,30,31,32]. Most studies and applications are limited to conventional bi-material systems with 0.8 α 0.8 and neglect the effects of the mismatch of the volumetric stiffness of the layers; this is done by assuming the second Dundurs’ parameter β = 0 “because of its secondary role”. In addition, most studies neglect the effects of crack tip shear, V in Figure 1b, because of its limited importance when dealing with long cracks. However, new material combinations with a larger mismatch of the elastic constants are increasingly being used [2,3,8,11,33,34,35], where the effects of β 0 on the fracture parameters might be important. Table 1 lists examples of thin-layer/substrate material combinations used in electronics, structural composite sandwiches, and thermal barrier coatings. In addition, shear does contribute to the fracture parameters in many applications and thin-layer tests, such as the blister/bulge test and the peel test for thin films [26,31], or the Single Cantilever Beam (SCB) test for sandwich structures [36]. The analytical and semi-analytical expressions presented in this paper can be used for the solution of various film/substrate problems, for instance for interface toughness measurements or to study buckle-driven delamination. They can be used to validate and support the simplifying assumptions. They also provide a useful reference limiting solutions for problems where the substrate is not semi-infinite, e.g., in SCB sandwich tests.
The paper is organized as follows. Analytical solutions are presented in Section 2 and Appendix A for the interface crack in Figure 1 with a , β = 0 and a general mismatch of the effective Young’s moduli, 1 α 1 . A numerical model which implements a crack surface displacement extrapolation method into the finite element code ANSYS for the accurate definition of the fracture parameters is presented in Section 2 and Appendix B; the model is verified in Section 3 using the analytical results. In Section 2, the semi-analytic solutions in References [25,26] are particularized to the problem in Figure 1, and the numerical coefficients are derived and presented in tabular form in Section 3 for bi-materials with a general mismatch of the elastic constants and Dundurs’ parameters 0.99 α 0.998 and 0.4 β 0.4 . Extreme asymptotics and intermediate expansions of the fracture parameters are derived in Section 3. A discussion on the results is presented in Section 4.

2. Materials and Methods

The problem in Figure 1a, representing an edge-delamination crack of length a between a layer of thickness h 1 and a semi-infinite plane h 1 / h 2 0 , subjected to general edge loadings per unit width, P , M , V with M = M V a , is examined. The crack tip loads are defined by the elementary, self-equilibrated systems in Figure 1b: pure bending moments, M , axial forces, P (plus the compensating moment, M * = M + P h 2 / 2 , since h 1 / h 2 0 ) and shear forces, V [25,26].
The layer, material #1, and the half plane, material #2, are linearly elastic and isotropic, with Young’s moduli and Poisson ratios E 1 , ν 1 and E 2 , ν 2 ; the interface between the two materials is abrupt, smooth, and adherent for x 0 . The mismatch of the elastic constants is described by the Dundurs’ parameters [37]:
α = E ¯ 1 E ¯ 2 E ¯ 1 + E ¯ 2 β = G 1 ( κ 2 1 ) G 2 ( κ 1 1 ) G 1 ( κ 2 + 1 ) + G 2 ( κ 1 + 1 ) ε = 1 2 π ln ( 1 β 1 + β )
where E ¯ i = E i and κ i = ( 3 ν i ) / ( 1 + ν i ) for plane-stress, and E ¯ i = E i / ( 1 ν i 2 ) and κ i = 3 4 ν i for plane-strain; G i = E i / ( 2 + 2 ν i ) is the shear modulus. The ratio of the effective Young’s moduli is Σ = E ¯ 1 / E ¯ 2 and is related to α through α = ( Σ 1 ) / ( Σ + 1 ) . The parameters α , β define measures of the mismatch in uniaxial and bulk compliances of the two materials.
The crack length a is sufficiently long to ensure that the crack tip fields are unaffected by the actual distribution of the end tractions; since the layer is homogeneous and isotropic, this condition is satisfied if a is a few times the layer thickness, h 1 (Saint Venant’s principle). Details on the applicability of these solutions to systems where h 1 / h 2 0 , namely with finite substrates or finite thickness layers, are given in Section 4.
The stress and displacement fields in the singularity dominated zone at the crack tip are described by the stress and displacement components, σ y y , σ x y and u y , u x . For general interfaces, they are oscillating [38,39], and the oscillations are described by the oscillatory index ε given in Equation (1):
( σ y y + i σ x y ) θ = 0 = K 2 π r r i ε ( u y + i u x ) θ = π ( u y + i u x ) θ = π = 8 E * K r i ε + 1 / 2 2 π ( 1 + 2 i ε ) cosh ( π ε )
with r is the distance from the crack tip (ahead for θ = 0 or behind for θ = ± π ), Figure 1b, E * = ( 1 α ) E ¯ 1 , and K = K 1 + i K 2 is the complex stress intensity factor [1,40]. The oscillations imply crack surface interpenetration at the crack tip, which can be discounted if contact occurs in a small region near the crack tip (small-scale contact) [40]. For β = 0 , the index ε = 0 , the fields are similar to those in a homogeneous material and K I = K 1 and K I I = K 2 define the mode I and mode II components of the stress intensity factor.
The energy release rate for the collinear propagation of the crack is related to the modulus of the complex stress intensity factor, | K | = ( K 1 2 + K 2 2 ) 1 / 2 , by:
G = 1 β 2 E * | K | 2
[41]. The crack tip conditions are generally mixed mode and defined by the mode-mixity phase angle. For bi-material interfaces with β = 0 , K = K I + i K I I = K e i ψ , and the mode mixity phase angle is uniquely defined by:
ψ = tan 1 ( K I I / K I ) ,   for   ε = β = 0
For bi-material interfaces with β 0 , the mode-mixity phase angle is not uniquely defined due to the oscillating fields. Here, following [40], the mode-mixity angle is related to the ratio between the shear and normal stresses at a reference distance, r ^ , ahead of the crack tip:
ψ r ^ = tan 1 ( Im [ K r ^ i ε ] / Re [ K r ^ i ε ] ) ,   for   ε β 0
Results will be presented in this paper assuming r ^ = h 1 . The choice of r ^ is arbitrary, since the phase angle at a different reference length is given by ψ r ^ = ψ h 1 + ε ln ( r ^ / h 1 ) [40]; as a consequence, the interfacial fracture toughness defined through fracture experiments will depend on the mode mixity and on the value of r ^ used to define ψ r ^ [35,42].
For the geometry in Figure 1a, with h 1 / h 2 0 , the energy release rates and phase angles associated with the elementary crack tip loads, P , M , V , assume the forms [26]:
E ¯ 1 G M = f M M h 1 3 / 2 , E ¯ 1 G P = f P P h 1 1 / 2 , E ¯ 1 G V = f V ( α , β ) V h 1 1 / 2 ψ P = ω ( α , β ) , ψ M = ω π / 2 , ψ V = ψ V ( α , β )
where, f M = 6 and f P = 1 / 2 [19,25]. The phase angle ψ P = ω has been derived semi-analytically in Reference [25] for 0.4 β 0.4 and 0.8 α 0.8 and analytically in [28] for 1 α 1 and β = 0 ; f V and ψ V can be obtained from the analytical results in [28] for β = 0 (see below), and for 0.8 α 0.8 and β = 0 , they have been defined numerically in [26]. Results for orthotropic materials with α = β = 0 have been presented in References [43,44] for axial forces and bending moments and for shear in [27,45]. In this paper, solutions for f V , ψ V and ψ P = ω will be derived for isotropic materials with 0.99 α 0.998 and 0.4 β 0.4 .
For most applications of interface mechanics in the literature, the effect of the second Dundurs’ parameter on the solution is neglected by assuming β = 0 . This is done to simplify the treatment of the problem and avoid dealing with the oscillatory fields. This assumption is valid for many material combinations, especially when α is small, e.g., 0.4 α 0.4 , since β would then vary in the approximate range 0.1 < ˜ β < ˜ 0.1 , which is exact for ν 1 = ν 2 = 1 / 3 under plane-strain conditions [34]. However, the assumption is questionable for other cases and larger values of α (see Table 1).

2.1. Analytical Model: Materials with β = 0

An analytical solution of the problem in Figure 1 for β = 0 has been derived in [28] using the Laplace transform and reducing the problem to a homogeneous matrix Riemann-Hilbert problem following [46,47]; see Appendix A for more details. The solution has been obtained under the assumption that the loads are applied at an infinite distance from the crack tip. The crack is in mixed-mode conditions with the Mode I and mode II stress intensity factors given by:
{ K I K I I } = k M M h 1 3 / 2 + k P P h 1 1 / 2 + k V V h 1 1 / 2
where
k M = 6 ( 1 α ) ( sin ω cos ω ) , k P = ( 1 α ) 2 ( cos ω sin ω ) , k V = δ k M + 3 k P
and ω and δ are given in Equations (A1) and (A2); they depend on Σ = E ¯ 1 / E ¯ 2 (or α ) and have been tabulated for ease of application in Table 2.
The mode-mixity phase angles associated with the elementary loads in Figure 1b are derived from Equations (4), (7) and (8):
ψ P = ω , ψ M = ω π / 2 , ψ V = ω tan 1 ( 2 δ )
and the dimensionless functions which define the elementary energy release rates in Equation (6) are obtained as:
f M = 6 , f P = 1 / 2 , f V = 6 δ 2 + 3 2
where f M and f P coincide with those derived in [25]. The analytical results for ψ P = ω , ψ V , and f V will be presented in Section 3.3 for selected values of α ; extreme asymptotics and intermediate expansions, which describe very soft and stiff substrates and near homogeneous problems, will be derived in Section 3.3.
The global mode mixity angle when P , M , V 0 in Figure 1, is obtained through Equations (4) and (7), and the global energy release rate in Equation (3) becomes:
G = 1 2 E ¯ 1 h [ 12 ( M h 1 + δ V ) 2 + ( P + 3 V ) 2 ]
which coincides with results in [25] when V = 0 .

2.2. Computational Model

The finite element model used to analyze the problem in Figure 1 is presented in Appendix B. The model implements the Crack Surface Displacement Extrapolation (CSDE) technique presented in [48] to define the energy release rate and mode-mixity phase angle. The CSDE is based on a linear extrapolation technique, where the energy release rate and mode-mixity phase angle calculated from relative nodal pair displacements at locations behind the crack tip are extrapolated linearly to the crack tip location. The energy release rate and mode-mixity phase angle are defined in terms of relative crack opening, Δ u y , and shearing, Δ u x , displacements at a distance r behind the crack tip [40]:
G = ( Δ u y 2 + Δ u x 2 ) r ( 1 + 4 ε 2 ) π E * 32 r ψ = ψ h 1 = tan 1 ( Δ u x / Δ u y ) r ε ln ( r / h 1 ) + tan 1 ( 2 ε )
with Δ u i = ( u i ) θ = π ( u i ) θ = π . The displacements are calculated using the finite element model presented in Appendix B.
The CSDE technique was originally developed to solve convergence problems for high stiffness mismatch bi-material interfaces, e.g., sandwich face/core interfaces, encountered with other mode separation methods, such as the Virtual Crack Closure Technique (VCCT), where nodal forces at or close to the crack tip are utilized to calculate the mode-mixity. Beuth [49] showed that the traditional VCCT by Rybicki and Kanninen [50] is not able to generate converged results for the mode-mixity due to the near-tip stress and displacement oscillations. Therefore, a modified version of the VCCT was presented in [49] where the near-tip oscillations are neutralized analytically during the post-processing of the near tip nodal force and displacement results to calculate the mode-mixity, and in this way achieve converged mode-mixity results. However, in most numerical finite element representations for applied bi-material crack problems, nodal results at or close to the crack tip are often associated with high numerical errors. This is both due to high stiffness mismatches over the interface, and thus the inability of the elements to accurately capture the complex stress field near the crack tip, but also to the inability of most models to capture the oscillations close to the crack tip. As a result, even the modified VCCT presented in [49] is unable to accurately predict the mode-mixity for high stiffness mismatch bi-material problems, like sandwich face/core interfaces [51]. The CSDE method, due to its linear extrapolation technique, bypasses the numerical error-zone close to the crack tip, and therefore constitutes a very stable and robust mode separation technique which is able to calculate the mode-mixity phase angle for even very challenging extreme stiffness mismatch bi-material crack problems [51].

2.3. Semi-Analytical Model—General Material Mismatch β 0

Following the formulations in [25,26], the energy release rate and mode mixity phase angle for the bi-material system in Figure 1 are defined in terms of crack tip force and moment resultants, P , M , V , by operating in the complex plane of the stress intensity factors and noting that for each elementary load, e.g., V, the following relations hold: K V h 1 i ε = | K V | e i ψ V and | K V | = ( 1 α ) ( 1 β 2 ) 1 f V V h 1 1 / 2 . This and Equation (6) yield:
G = 6 E ¯ 1 h 1 3 M 2 + 1 2 E ¯ 1 h 1 P 2 + f V 2 E ¯ 1 h 1 V 2 + 2 6 f V E ¯ 1 h 1 2 cos ( ψ M ψ V ) M V   + 2 f V E ¯ 1 h 1 cos ( ω ψ V ) P V ψ = tan 1 ( Im [ K h 1 i ε ] Re [ K h 1 i ε ] ) = tan 1 ( 2 3 M sin ψ M + P h 1 sin ω + 2 f V V h 1 sin ψ V 2 3 M cos ψ M + P h 1 cos ω + 2 f V V h 1 cos ψ V )
where ψ M = ω π / 2 . For β = 0 , f V , ω and ψ V , are given in Equations (9) and (10), and G coincides with Equation (11); the expressions in Equation (13) are then fully analytical. For β 0 , ω ( α , β ) , f V ( α , β ) and ψ V ( α , β ) have to be defined numerically; this will be done in Section 3.3. A discussion on the applicability of the analytical and semi-analytical solutions to problems where substrate and crack have finite sizes is presented in Section 4.

3. Results

In this section, the accuracy of the finite element model and the CSDE technique for the challenging problem in Figure 1 is verified using the analytical results for ω , ψ V and f V in materials with β = 0 . The finite element model is then applied to derive the dimensionless coefficients of the semi-analytic solutions in Equation (13) for β 0 . Extreme asymptotics and intermediate expansions for the coefficients of the analytical solutions are derived to facilitate the application of the analytical results to systems with no mismatch in the volumetric stiffnesses.

3.1. Finite Element Results for Materials with β = 0

The finite element model, Appendix B, and the CSDE technique implemented in ANSYS are first applied to define the energy release rate, stress intensity factors, and mode mixity angles associated with the elementary loads in Figure 1b for a bi-material interface with Dundurs’ parameter β = 0 . The phase angle ψ P = ω associated to the elementary axial load P (plus the compensating bending moment) are presented in Table 3 and Figure 2 on varying α and compared with the analytical solutions in Equations (7)–(9). Phase angle, energy release rate, and stress intensity factors for the elementary load of pure crack tip shear, V , are presented in Table 4 and Table 5 and Figure 3.
The absolute error between the finite element based CSDE results and the analytical results is always below 0.04° on the phase angle, and below 0.001 on the dimensionless stress intensity factors for 0.8 α 0.8 ; errors are slightly larger for 0.99 α 0.9 and 0.9 α 0.99 , but always below 0.2° and 0.01 . The results have also been successfully compared with solutions in the literature [25,26] (a discrepancy is observed in the solutions for pure shear loading and α 0 in [26]; given the agreement with the analytical solution, the present results appear to be more accurate). The results in Figure 2 and Figure 3 show large variations for α > 0.8 , which corresponds to stiff layers on soft material substrates. These cases will be examined later in Section 3.3, where extreme asymptotics of the analytical solutions will be derived.

3.2. Tabulated Coefficients for the Semi-Analytic Solutions in Equation (13)—Materials with β 0

The finite element model and the CSDE technique have been applied to the geometry in Figure 1 to define energy release rate and mode mixity phase angle for the elementary loads in bi-material interfaces with general elastic mismatch and β 0 .
The dimensionless function f V = E ¯ 1 h G V / V and the mode-mixity phase angle ψ V are presented in Table 4 and Table 5 for 0.99 α 0.998 and 0.4 β 0.4 . The phase angle ψ P = ω for 0.99 α 0.9 and 0.9 α 0.998 and 0.4 β 0.4 , which has not been presented before in the literature, has also been derived, and results are listed in Table 3; these results can be combined with those tabulated in [25] for 0.8 α 0.8 to cover all material combinations (a verification performed on the results in [25] using the analytical model, Equation (7), and the finite element model formulated here indicates an uncertainty of ± 0.5 ° ). The uncertainties are given in the tables and have been defined through comparison with the analytical results and results in [25,26]. They are mainly determined by cases with α 1 , for which the numerical results are less accurate. However, as it evident also from Figure 3b, Figure 4a, Figure 5 and Figure 6, the coefficients approach a plateau for large negative values of α and the numerical derivation for α 1 is not necessary.

3.3. Explicit Solutions for Soft/Stiff Material Substrates and Nearly Homogeneous Systems— β = 0

To facilitate the application of the analytical results presented in Section 2, extreme asymptotics and intermediate expansions of ω , ψ V and f V are derived for stiff layers over soft substrates, where Σ = E ¯ 1 / E ¯ 2 > > 1 is very large and α 1 , for soft layers over stiff substrates, where Σ = E ¯ 1 / E ¯ 2 < < 1 is very small and α 1 , and near homogeneous cases, where Σ = E ¯ 1 / E ¯ 2 1 and α 0 . They are listed below.
  • Phase angle associated with the axial force, ψ P = ω (in radians) [28]:
    Large   Σ = E ¯ 1 / E ¯ 2 : ω ( Σ ) = π 2 ( 4 3 ) 1 / 6 Σ 1 / 3 + o ( Σ 1 / 3 )
    Small   Σ = E ¯ 1 / E ¯ 2 : ω ( Σ ) = 0.794 + 0.171 Σ 0.0934 Σ 2 + O ( Σ 3 )
    Intermediate   Σ = E ¯ 1 / E ¯ 2 1 : ω ( α ) = 0.909 + 0.159 α + 0.0697 α 2 + 0.0421 α 3 + O ( α 4 )
    The uncertainties are below 0.2 for E ¯ 1 / E ¯ 2 330 ( α 0.993 ) in Equation (14), for E ¯ 1 / E ¯ 2 0.4 ( α 0.43 ) in Equation (15), and for 0.4 E ¯ 1 / E ¯ 2 3.2 ( 0.65 α 0.51 ) in Equation (16). These uncertainties imply a relative percent error always lower than 0.4%. The equations can be used with confidence for all values of α outside the interval 0.51 < α < 0.993 where the exact analytical solutions in Equation (A1) or interpolation of the results in Table 3 and Figure 2 should be applied.
  • Phase angle associated with shear, ψ V (in radians):
    Large   Σ = E ¯ 1 / E ¯ 2 : ψ V ( Σ ) = 2 5 / 3 3 1 / 6 Σ 1 / 3 + o ( Σ 1 / 3 )
    Small   Σ = E ¯ 1 / E ¯ 2 : ψ V ( Σ ) = 0.099 0.146 Σ + 0.113 Σ 2 + O ( Σ 3 )
    Intermediate   Σ = E ¯ 1 / E ¯ 2 1 : ψ V ( Σ ) = 0.0168 0.0917 α 0.00913 α 2 + O ( α 3 )
    The uncertainties are below 0.2 for E ¯ 1 / E ¯ 2 330 ( α 0.994 ) in Equation (17), for E ¯ 1 / E ¯ 2 0.37 ( α 0.46 ) in Equation (18), and for 0 E ¯ 1 / E ¯ 2 6.2 ( 1 α 0.72 ) in Equation (19). These uncertainties imply a relative percent error higher than that on ω , since ψ V is much smaller, which can be as high as 9% (Figure 2 and Figure 3). However, this occurs when ψ V is closed to zero, and the problem is essentially mode I, and the error is therefore irrelevant for engineering applications. The equations can be used with confidence for all values of α outside the interval 0.72 < α < 0.994 , where the exact analytical solutions in Equations (9) and (A2), or interpolation of the results in Table 4, should be applied.
  • The dimensionless function which defines the energy release rate associated with shear, f V :
    Large   Σ = E ¯ 1 / E ¯ 2 : f V ( Σ ) = 6 δ 2 + 3 2 with δ = 2 2 / 3 3 5 / 6 Σ 1 / 3 + o ( Σ 1 / 3 )
    Small Σ = E ¯ 1 / E ¯ 2 : f V ( Σ ) = 6 δ 2 + 3 2 ,   with δ = 0.417 + 0.269 Σ 0.104 Σ 2 + o ( Σ 3 )
    Intermediate Σ = E ¯ 1 / E ¯ 2 1 : f V ( Σ ) = 6 δ 2 + 3 2 ,   with δ = 0.62 + α π + 5 8 π α 2 + o ( α 3 )
The relative error is below 1% for E ¯ 1 / E ¯ 2 200 ( α 0.99 ) in Equation (20), for E ¯ 1 / E ¯ 2 0.6 ( α 0.25 ) in Equation 21, and for 0.4 E ¯ 1 / E ¯ 2 2.1 ( 0.44 α 0.36 ) in Equation (22). The equations can be used with confidence for all values of α outside the interval 0.36 < α < 0.99 where the exact analytical solutions in Equations (10) and (A2), or the results in Table 5, should be applied. If a maximum relative error of 5% is acceptable, the validity of Equation (20) can be extended to all E ¯ 1 / E ¯ 2 0.36 ( α 0.47 ).
The asymptotics for small Σ = E ¯ 1 / E ¯ 2 and for Σ 1 have been obtained by developing the integrands of the corresponding values into series followed by integration. The asymptotics for large Σ = E ¯ 1 / E ¯ 2 have been obtained by first substituting s = s Σ 1 / 3 into Equations (A1) and (A2) and then developing them into series. The asymptotics (20) which define the energy release rate associated with shear in the limit of a soft substrate coincides with that obtained by approximating the thin layer as a beam (plate) [52]. An approximate expression for ω , valid in the range 0.8 α 0.8 , and obtained through a numerical fit of the results in [25] is presented in [3].
Figure 4 shows results for ψ P = ω as a function of E ¯ 1 / E ¯ 2 in a logarithmic scale to highlight the behaviors. The analytical results are compared with the numerical results obtained here and presented in Table 5 for 0.99 α 0.8 and 0.8 α 0.999 , and with those [25] for 0.8 α 0.8. The figure shows the two asymptotics in Equations (14) and (15) for small and large values of E ¯ 1 / E ¯ 2 and the asymptotics presented in [2] for large E ¯ 1 / E ¯ 2 . The latter has been obtained using the scaling, K I / K I I ( E ¯ 1 / E ¯ 2 ) 1 / 3 , derived from the model of a beam on a semi-infinite substrate, and a numerical fit for α = 0.95 ; comparison with the analytical results and the analytical asymptotic in Equation (14) shows that the approximation is quite accurate for large values of Σ = E ¯ 1 / E ¯ 2 .
Figure 5 shows results for ψ V as function of E ¯ 1 / E ¯ 2 on a logarithmic scale to highlight behavior in the extreme cases. The analytical results are compared with the numerical results in Table 4, for 0.99 α 0.998 . The figure also shows the two asymptotics in Equations (17) and (18), for small and large values of E ¯ 1 / E ¯ 2 , and the intermediate expansion in Equation (19).

4. Discussion

Analytical expressions for the energy release rate and mode-mixity phase angle associated with the elementary loads in Figure 1b, Equations (9) and (10), have been derived from the results in Reference [28] for materials with β = 0 , and depend on Σ = E ¯ 1 / E ¯ 2 or α = ( E ¯ 1 E ¯ 2 ) / ( E ¯ 1 + E ¯ 2 ) . Useful asymptotics solutions have been derived in Equations (14)–(22) for stiff layers over soft substrates, where α 1 or E ¯ 1 / E ¯ 2 , e.g., metal coatings on elastomers, and for soft layers over stiff substrates, where α 1 or E ¯ 1 / E ¯ 2 0 , e.g., polymer coatings over metals; furthermore, intermediate expansions have been derived for near-homogeneous problems.
The effects of β 0 on the fracture parameters are presented in Table 3, Table 4 and Table 5. In the absence of shear forces, for V = 0 ,   M , P 0 , the energy release rate, Equation (13), does not depend on β and the effects on the mode mixity phase angle ψ P = ω are synthesized in Figure 6. For stiff layers over soft substrates, where α 1 or E ¯ 1 / E ¯ 2 , the second Dundurs parameter β = 1 / 2 ( 1 2 ν 2 ) / ( 1 ν 2 ) does not depend on the Poisson ratio of the layer, which can be approximated as a beam [2,52,53]. The effects of β 0 in this limit are not very pronounced, and this implies that (i) the assumption of incompressible substrate β = 0 is accurate and (ii) the fracture mechanics parameters could be obtained using the approximation of a beam on a Winkler foundation, for M , or the membrane approximation, for P . For soft layers over stiff substrates, where α 1 or E ¯ 1 / E ¯ 2 0 , the second Dundurs parameter β = 1 / 2 ( 1 2 ν 1 ) / ( 1 ν 1 ) depends on the Poisson ratio of the layer only, and the beam/membrane approximations are not acceptable. The effects of β 0 are important and controlled by the near-tip deformations in the layer ahead of the crack tip. These effects should not be discounted a priori.
Similar conclusions apply to the mode mixity phase angle associated with pure shear loading V 0 ,   M , P = 0 . The energy release rate associated with pure shear, Table 4, is also affected by β 0 , even if the effect is weaker. The effects are limited for intermediate values of α and become negligible in the limit α 1 where the solution does not depend on the Poisson ratios of the substrate and layer and the approximation of a beam over a Winkler foundation is accurate. In the limit α 1 , neglecting the compressibility of the layer could lead to errors as high as 10%. The assumption β = 0 should therefore be used with care. This conclusion is supported by experimental results on bi-material systems which show an important variation of the interfacial fracture toughness on varying the mode mixity, especially when crack tip conditions are mode II dominated [2,35,42].
Material systems with β 0 can be studied using the semi-analytical solutions in Equation (13) and the numerical coefficients in Table 3, Table 4 and Table 5. They can be applied also to problems with large deformations, e.g., buckle-driven delamination or peeling problems if the crack tip forces are previously calculated using a model which accounts for the geometrical nonlinearities [16,26]. Both analytical and semi-analytical solutions presented in this paper account for the presence of crack tip shear. This contribution is often neglected under the assumption that the effect of the shear becomes negligible for long cracks. This is indeed correct. However, a simple example may clarify the importance of shear in problems where the layer is subjected to a transverse load, e.g., in a blister test, peel test, or SCB face/core interface characterization sandwich testing. Consider the problem in Figure 1 and assume that only a transverse load V 0 acts at a distance a from the crack tip (similar to a shaft-loaded blister test but with free rotations at the layer end). The transverse load will induce a bending moment, M = V a , and shear, V , at the crack tip, Figure 1b. The energy release rate of Equation (13) then simplifies:
G = 6 E ¯ 1 h 1 3 ( V a ) 2 + f V 2 E ¯ 1 h 1 V 2 + 2 6 f V E ¯ 1 h 1 2 cos ( ψ M ψ V ) V 2 a  
The dominant term in G is related to the crack tip bending moment and depends quadratically on the crack length. The term which is directly related to the crack tip shear (the second on the righthand side) is independent of the crack length, and its effect is negligible apart from the case of very short cracks. The effects of shear are important through the third term, which depends linearly on the crack length and has a non-negligible effect also for intermediate crack lengths. The mechanical significance of this term has been first explained in [27]. It describes the effect of crack tip shear on the near-tip deformations generated by the bending moment. These deformations are typically modeled as a crack tip root-rotation, which is the rotation the cross-section experiences at the crack tip (here assumed as positive if clockwise) due to the elastic deformability. If the layer is sufficiently long, the crack tip root rotations are related to the dimensionless crack tip bending moment, ( V a ) / ( E ¯ 1 h 1 2 ) , through a compliance coefficient, a 1 M . For the problem in Figure 1, the compliance coefficient is a 1 M = 2 6 f V cos ( ψ M ψ V ) and the third term on the righthand side of Equation (23) defines the work done by the crack tip shear force on the crack tip root rotation. Equation (23) then shows that, in the presence of shear, the crack tip root rotation due to the crack tip bending moment does contribute to the energy release rate. The effect is important for stiff layers over soft substrates, where the function f V becomes very large (Figure 3b and Table 3). In the absence of shear, the root rotations due to pure bending would not have any effects on the fracture parameters [25], but only on the deformability of the system, e.g., in buckle-driven delamination [15,16,54]. These concepts have recently been discussed in [55] for sandwich beams and in [56] for bi-material DCB specimens.
The solutions derived here have been obtained assuming the crack faces to be traction free and the loads to act at a sufficient distance from the crack tip. However, they can be used to obtain approximate solutions for problems where the loads act as distributed tractions applied along the crack faces, e.g., the blister test or pressure-loaded debonds, using the crack tip force and moment resultants produced by the loads at the crack tip [44,45,57].
The results obtained in this paper can be applied to problems where the geometry has finite dimensions h 2 × 2 c and h 1 / h 2 0 (with h 2 and 2 c depth and width, see Figure A1). This is possible if a , h 2 , and 2 c are above minimum values, which must ensure that (i) the crack tip fields are unaffected by the actual distribution of the end loads applied to the thin layer and by the presence of boundaries (free or clamped) and (ii) the stresses generated at the crack tip do not interact with the boundaries. For Saint Venant’s principle, the crack length should be just a few times the thickness of the layer h 1 . The analyses in [2] using the model of a beam on a semi-infinite half-plane indicate that the width of the domain should be c / h 1 > 3 E ¯ 1 / E ¯ 2 to ensure that the interfacial tractions ahead of the crack tip drop to almost zero. This would imply c / h > 3000 , 3 × 10 4 for E ¯ 1 / E ¯ 2 = 1 × 10 3 , 1 × 10 4 or α = 0.998 , 0.999 . The calculations performed here and the results in Table 3, Table 4 and Table 5 show that good engineering accuracy is achieved on the fracture parameters, for 0.99 α 0.999 , with much smaller widths ( c / h = 1000 ). This conclusion is supported by analyses performed in [8,55,58] on three-layers structures. According to [2], the depth of the substrate should be h 2 / h 1 > 10 E ¯ 1 / E ¯ 2 ; again, the convergence analyses performed here indicate that a smaller depth is necessary to obtain engineering accuracy for large E ¯ 1 / E ¯ 2 ( h 2 / h 1 = 1000 ). The problem of a corner interface crack, where the crack emerges at the left substrate edge, imposes different restrictions on the crack length a , which are similar to those for c [12,59].
The results and methods in this paper can be extended to orthotropic materials using an orthotropic rescaling technique [43,44,45,60]. For the problem in Figure 1, this extension might be important when dealing with debonds in sandwich structures where the core is often orthotropic. The models applied here to the thin-layer/substrate system can be extended to general bi-material interfaces in beam-like structures. The extension allows for accurate and efficient calculation of the fracture parameters, avoiding case specific numerical calculations, in a large number of fracture tests used for material characterization, e.g., Double-Cantilever Beam tests, mixed-mode bending tests, end-notched flexural tests, four-point bending, and inverted four-point bending tests. Analytical solutions for bi-material beams with h 1 / h 2 = 1 and β = 0 have been recently derived in [61], and semi-analytical solutions for beams with h 1 / h 2 1 and β 0 are presented in [62].

5. Conclusions

Analytical and semi-analytical solutions have been derived for the energy release rate and mode mixity phase angle of an edge-delamination between a thin layer and an infinitely deep substrate. The 2D solutions are in terms of three elementary crack tip load systems, axial forces, shear forces, and bending moment, and can be used to define the fracture parameters for general loading conditions. Isotropic materials with an exceptionally large mismatch of the elastic constants have been analyzed, and solutions derived for Dundurs’ parameters 1 α 1 and 0.4 β 0.4 . The solutions account for the effects of shear, which can be important also for intermediate crack lengths, especially when the near-tip deformations are large as in the presence of soft substrates.
Extreme asymptotics and intermediate expansions have been presented for materials with β = 0 , which highlight the scaling of the fracture parameters on varying E ¯ 1 / E ¯ 2 and provide useful closed-form solutions applicable to stiff layers over soft substrates, to soft layers over stiff substrates and to near homogeneous problems. Results show that for stiff layers over soft substrates ( α 1 ), the effects of material compressibility, β 0 , are weak; the assumption of substrate incompressibility is accurate as well as the approximation of the layer as a beam on an elastic foundation and also as a beam on a Winkler foundation. For other combinations, including soft layers over stiff substrates ( α 1 ), the effects may be relevant and problem specific.
A computational model, previously formulated and verified for the derivation of the fracture parameters in bi-material beam-like structures, which uses a crack surface displacement extrapolation technique for mode partitioning, has been verified, and its accuracy proven also for this numerically challenging problem (long slender layer + thick substrate + bi-materials with large mismatch of the elastic constants).

Author Contributions

Conceptualization, R.M., C.B. and K.U.; Methodology and Formal Analysis, C.B., R.M. and K.U.; Investigation, R.M., C.B. and K.U.; Software, L.B., C.B.; Validation, R.M., K.U. and L.B.; Writing, R.M., C.B. and K.U.

Funding

This research was funded by the U.S. Navy, Office of Naval Research, ONR (Nos. N00014-17-1-2914 (Roberta Massabò) and N00014-16-1-2977 (Christian Berggreen)) monitored by Yapa D.S. Rajapakse; by Italian Department for University and Scientific and Technological Research, MIUR Prin15 project (No. 2015LYYXA8); by the Russian State Assignment under contract (No. AAAA- А17-117021310371-9 (Konstantin Ustinov)).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The analytical solution [28] has been obtained by relating the Laplace transforms of the discontinuous (normal and tangential) displacement derivatives in the strip and in the half plane along the delaminated part with the (normal and shear) stresses acting along the intact part of the interface by means of a 2 × 2 matrix, and reducing therefore the problem under consideration to a matrix Riemann-Hilbert problem. For the particular restriction β = 0 , the matrix equation have been solved using the technique in [30] which yields the expressions for the stress intensity factors in Equations (7) and (8). The quantities ω and δ in the equations depend on Σ = E ¯ 1 / E ¯ 2 and are given by:
ω ( Σ ) = π 4 1 π 0 ( 1 + s 2 ) 1 / 2 log [ A B s coth s 1 + s 2 + 1 ] d s
δ ( Σ ) = 4 log 2 π + 1 π 0 d d s { log A B + ( 1 + s 2 ) 1 / 2 log [ A B s coth s 1 + s 2 + 1 ] } d s s
with
2 A = Σ tanh 3 s + tanh s ( tanh s + s 1 + s 2 / cosh 2 s ) / ( 1 s 2 / sinh 2 s ) 2 B = Σ tanh s + ( 1 s 1 + s 2 / ( sinh s cosh s ) ) / ( 1 s 2 / sinh 2 s )
Here ω is given in radians. The solution above covers the general case α = ( Σ 1 ) / ( Σ + 1 ) 0 and coincides, when the half plane and the strip are of the same material and α = 0 , with the solution previously derived in [30]. Insurmountable difficulties arise when the technique is applied to the most general case α 0 , β 0 or even for α = 0 , β 0 , for which cases a numerical solution is required (see Section 2).

Appendix B. Finite Element Model

A sketch of the finite element model is shown in Figure A1. The model is defined in the finite element code ANSYS and the model definitions are shown in Figure A1. The dimensions of the domain used to approximate the semi-infinite plane and mesh refinement have been defined through a convergence study.
Figure A1. (a) Dimensions and domains of the finite element model in ANSYS. (b) Zooms of the Finite element mesh showing internal, intermediate and external domains (parts of).
Figure A1. (a) Dimensions and domains of the finite element model in ANSYS. (b) Zooms of the Finite element mesh showing internal, intermediate and external domains (parts of).
Coatings 09 00152 g0a1
The mesh is subdivided into three domains. The external domain uses 8-noded isoparametric elements (PLANE183) with a fine mesh near the layer (size no larger than 5 / 8 h 1 ) and a coarser mesh at the lower boundary (size up to 80 h 1 ). The intermediate domain uses 8-noded isoparametric elements with size no larger than 5 / 8 h 1 . The internal domain is composed of a succession of approximately 200 rings with elements of decreasing size towards the crack tip; the 8 innermost rings surrounding the crack tip are composed of 256, 4-noded, isoparametric elements (PLANE 182) of size h 1 / 4000 ; the 75 surrounding rings are composed of 8-noded isoparametric elements of size h 1 / 1500 . The remaining 117 rings are used to smoothly transition to the intermediate mesh.

References

  1. Hutchinson, J.W.; Suo, Z. Mixed Mode Cracking in Layered Materials. Adv. Appl. Mech. 1991, 29, 63–191. [Google Scholar]
  2. Begley, M.R.; Hutchinson, J.W. The Mechanics and Reliability of Films, Multilayers and Coatings; Cambridge University Press: Cambridge, UK, 2017; ISBN 9781316443606. [Google Scholar]
  3. Freund, L.B.; Suresh, S. Thin Film Materials: Stress, Defect Formation, and Surface Evolution; Cambridge University Press: Cambridge, UK, 2003; ISBN 9780511165658. [Google Scholar]
  4. Carlsson, L.A.; Kardomateas, G.A. Structural and Failure Mechanics of Sandwich Composites; Springer: Dordrecht, The Netherlands, 2011; ISBN 9781402032240. [Google Scholar]
  5. Rinker, M.; Krueger, R.; Ratcliffe, J. Analysis of an Aircraft Honeycomb Sandwich Panel with Circular Face Sheet/Core Disbond Subjected to Ground-Air Pressurization; NASA/CR-2013-217974; NASA Center for AeroSpace Information: Hanover, MD, USA, 2013. [Google Scholar]
  6. Chen, Z.M.; Krueger, R.; Rinker, M. Face Sheet/Core Disbond Growth in Honeycomb Sandwich Panels Subjected to Ground-Air-Ground Pressurization and In-Plane Loading. In Proceedings of the NAFEMS World Congress 2015, San Diego, CA, USA, 21–24 June 2015. NF1676L-19962. [Google Scholar]
  7. Farshidi, A.; Berggreen, C. Analysis of disbonded aircraft sandwich panels with unvented honeycomb cores. Int. J. Rev. 2018. under review. [Google Scholar]
  8. Østergaard, R.C.; Sørensen, B.F. Interface crack in sandwich specimen. Int. J. Fract. 2007, 143, 301–316. [Google Scholar] [CrossRef]
  9. Balaguru, P.N.; Nanni, A.; Giancaspro, J. FRP Composites for Reinforced and Prestressed Concrete Structures: A Guide to Fundamentals and Design for Repair and Retrofit; Taylor & Francis: Abingdon, UK, 2009; ISBN 9780415448543. [Google Scholar]
  10. Yu, H.; Bai, Y.-L.; Dai, J.-G.; Gao, W.-Y. Finite element modeling for debonding of FRP-to-concrete interfaces subjected to mixed-mode loading. Polymers 2017, 9, 438. [Google Scholar] [CrossRef]
  11. Vaunois, J.-R.; Poulain, M.; Kanouté, P.; Chaboche, J.-L. Development of bending tests for near shear mode interfacial toughness measurement of EB-PVD thermal barrier coatings. Eng. Fract. Mech. 2017, 171, 110–134. [Google Scholar] [CrossRef]
  12. Hutchinson, R.G.; Hutchinson, J.W. Lifetime assessment for thermal barrier coatings: Tests for measuring mixed mode delamination toughness. J. Am. Ceram. Soc. 2011, 94, S85–S98. [Google Scholar] [CrossRef]
  13. Evans, A.G.; Mumm, D.R.; Hutchinson, J.W.; Meier, G.H.; Pettit, F.S. Mechanisms controlling the durability of thermal barrier coatings. Prog. Mater. Sci. 2001, 46, 505–553. [Google Scholar] [CrossRef]
  14. Di Leo, C.V.; Luk-Cyr, J.; Liu, H.; Loeffel, K.; Al-Athel, K.; Anand, L. A new methodology for characterizing traction-separation relations for interfacial delamination of thermal barrier coatings. Acta Mater. 2014, 71, 306–318. [Google Scholar] [CrossRef] [Green Version]
  15. Cotterell, B.; Chen, Z. Buckling and cracking of thin films on compliant substrates under compression. Int. J. Fract. 2000, 104, 169–179. [Google Scholar] [CrossRef]
  16. Yu, H.-H.; Hutchinson, J.W. Influence of substrate compliance on buckling delamination of thin films. Int. J. Fract. 2002, 113, 39–55. [Google Scholar] [CrossRef]
  17. Volinsky, A.A.; Moody, N.R.; Gerberich, W.W. Interfacial toughness measurements for thin films on substrates. Acta Mater. 2002, 50, 441–466. [Google Scholar] [CrossRef]
  18. Pelassa, M.; Massabò, R. Explicit solutions for multi-layered wide plates and beams with perfect and imperfect bonding and delaminations under thermo-mechanical loading. Meccanica 2015, 50, 2497–2524. [Google Scholar] [CrossRef] [Green Version]
  19. Thouless, M.D.; Evans, A.G.; Ashby, M.F.; Hutchinson, J.W. The edge cracking and spalling of brittle plates. Acta Metall. 1987, 35, 1333–1341. [Google Scholar] [CrossRef]
  20. van den Bosch, M.J.; Schreurs, P.J.G.; Geers, M.G.D. Identification and characterization of delamination in polymer coated metal sheet. J. Mech. Phys. Solids 2008, 56, 3259–3276. [Google Scholar] [CrossRef]
  21. Jiménez-Pardo, I.; van der Ven, L.; van Benthem, R.; de With, G.; Esteves, A. Hydrophilic self-replenishing coatings with long-term water stability for anti-fouling applications. Coatings 2018, 8, 184. [Google Scholar] [CrossRef]
  22. Xiao, J.; Gong, C.; Qi, M.; Jiag, A.; Wang, Z.; Li, M.; Ma, J. Effect of TiN/C microstructure composite layer on the adhesion of FDLC film onto silicon substrate. Coatings 2018, 8, 18. [Google Scholar] [CrossRef]
  23. Rao, J.; Sharma, A.; Rose, T. Titanium aluminium nitride and titanium boride multilayer coatings designed to combat tool wear. Coatings 2018, 8, 12. [Google Scholar] [CrossRef]
  24. Collino, R.R.; Philips, N.R.; Rossol, M.N.; McMeeking, R.M.; Begley, M.R. Detachment of compliant films adhered to stiff substrates via van der Waals interactions: Role of frictional sliding during peeling. J. R. Soc. Interface 2014, 11. [Google Scholar] [CrossRef] [PubMed]
  25. Suo, Z.; Hutchinson, J.W. Interface crack between two elastic layers. Int. J. Fract. 1990, 43, 1–18. [Google Scholar] [CrossRef]
  26. Li, S.; Wang, J.; Thouless, M.D. The effects of shear on delamination in layered materials. J. Mech. Phys. Solids 2004, 52, 193–214. [Google Scholar] [CrossRef]
  27. Andrews, M.G.; Massabò, R. The effects of shear and near tip deformations on energy release rate and mode mixity of edge-cracked orthotropic layers. Eng. Fract. Mech. 2007, 74, 2700–2720. [Google Scholar] [CrossRef]
  28. Ustinov, K. On separation of a layer from the half-plane: Elastic fixation conditions for a plate equivalent to the layer. Mech. Solids 2015, 50, 62–80. [Google Scholar] [CrossRef]
  29. Obreimoff, J.W. The splitting strength of mica. Proc. R. Soc. A Math. Phys. Eng. Sci. 1930. [Google Scholar] [CrossRef]
  30. Zlatin, A.N.; Khrapkov, A.A. Semi-infinite crack parallel to the boundary of an elastic half-plane. Sov. Phys. Dokl. 1986, 31, 1009. [Google Scholar]
  31. Jensen, H.M.; Thouless, M.D. Effects of residual stresses in the blister test. Int. J. Solids Struct. 1993, 30, 779–795. [Google Scholar] [CrossRef]
  32. Sundaram, S.; Lipkin, D.M.; Johnson, C.A.; Hutchinson, J.W. The influence of transient thermal gradients and substrate constraint on delamination of thermal barrier coatings. J. Appl. Mech. Trans. ASME 2013, 80. [Google Scholar] [CrossRef]
  33. Parry, G.; Colin, J.; Coupeau, C.; Foucher, F.; Cimetière, A.; Grilhé, J. Effect of substrate compliance on the global unilateral post-buckling of coatings: AFM observations and finite element calculations. Acta Mater. 2005, 53, 441–447. [Google Scholar] [CrossRef]
  34. Suga, T.; Elssner, G.; Schmauder, S. Composite parameters and mechanical compatibility of material joints. J. Compos. Mater. 1988, 22, 917–934. [Google Scholar] [CrossRef]
  35. Liechti, K.M.; Chai, Y.S. Asymmetric shielding in interfacial fracture under in-plane shear. J. Appl. Mech. Trans. ASME 1992, 59, 295–304. [Google Scholar] [CrossRef]
  36. Ratcliffe, J.G.; Reeder, J.R. Sizing a single cantilever beam specimen for characterizing facesheet-core debonding in sandwich structure. J. Compos. Mater. 2011, 45, 2669–2684. [Google Scholar] [CrossRef]
  37. Dundurs, J. Discussion: Edge-bonded dissimilar orthogonal elastic wedges under normal and shear loading. J. Appl. Mech. Trans. ASME 1969, 36, 650–651. [Google Scholar] [CrossRef]
  38. Williams, M.L. The stresses around a fault or crack in dissimilar media. Bull. Seismol. Soc. Am. 1959, 49, 199–204. [Google Scholar]
  39. Salganik, R.L. The brittle fracture of cemented bodies. J. Appl. Math. Mech. 1963, 27, 1468–1478. [Google Scholar] [CrossRef]
  40. Rice, J.R. Elastic fracture mechanics concepts for interfacial cracks. J. Appl. Mech. 1988, 55, 98–103. [Google Scholar] [CrossRef]
  41. Malyshev, B.M.; Salganik, R.L. The strength of adhesive joints using the theory of cracks. Int. J. Fract. Mech. 1965, 1, 114–128. [Google Scholar] [CrossRef]
  42. Banks-Sills, L. Interface fracture mechanics: Theory and experiment. Int. J. Fract. 2015, 191, 131–146. [Google Scholar] [CrossRef]
  43. Suo, Z. Delamination specimens for orthotropic materials. J. Appl. Mech. 1990, 57, 627–634. [Google Scholar] [CrossRef]
  44. Brandinelli, L.; Massabò, R. Mode II weight functions for isotropic and orthotropic double cantilever beams. Int. J. Fract. 2006, 139, 1–25. [Google Scholar] [CrossRef]
  45. Massabò, R.; Brandinelli, L.; Cox, B.N. Mode I weight functions for an orthotropic double cantilever beam. Int. J. Eng. Sci. 2003, 41, 1497–1518. [Google Scholar] [CrossRef]
  46. Gakhov, F.D. Boundary Value Problems; Pergamon Press: Oxford, UK, 1966; Volume 85. [Google Scholar]
  47. Rogosin, S.; Mishuris, G. Constructive methods for factorization of matrix-functions. IMA J. Appl. Math. 2016, 81, 365–391. [Google Scholar] [CrossRef]
  48. Berggreen, C.; Simonsen, B.C.; Borum, K.K. Experimental and numerical study of interface crack propagation in foam-cored sandwich beams. J. Compos. Mater. 2007, 41. [Google Scholar] [CrossRef]
  49. Beuth, J.L. Separation of crack extension modes in orthotropic delamination models. Int. J. Fract. 1996, 77, 305–321. [Google Scholar] [CrossRef] [Green Version]
  50. Rybicki, E.F.; Kanninen, M.F. A finite element calculation of stress intensity factors by a modified crack closure integral. Eng. Fract. Mech. 1977, 9, 931–938. [Google Scholar] [CrossRef]
  51. Berggreen, C. Damage Tolerance of Debonded Sandwich Structures. Ph.D. Thesis, Technical University of Denmark, Lyngby, Denmark, 2004. [Google Scholar]
  52. Salganik, R.L.; Ustinov, K.B. Deformation problem for an elastically fixed plate modeling a coating partially delaminated from the substrate (Plane Strain). Mech. Solids 2012, 47, 415–425. [Google Scholar] [CrossRef]
  53. Shield, T.W.; Kim, K.S. Beam theory models for thin film segments cohesively bonded to an elastic half space. Int. J. Solids Struct. 1992, 29, 1085–1103. [Google Scholar] [CrossRef]
  54. Ustinov, K.B. On influence of substrate compliance on delamination and buckling of coatings. Eng. Fail. Anal. 2015, 47, 338–344. [Google Scholar] [CrossRef] [Green Version]
  55. Barbieri, L.; Massabò, R.; Berggreen, C. The effects of shear and near tip deformations on interface fracture of symmetric sandwich beams. Eng. Fract. Mech. 2018, 201, 298–321. [Google Scholar] [CrossRef]
  56. Thouless, M.D. Shear forces, root rotations, phase angles and delamination of layered materials. Eng. Fract. Mech. 2018, 191, 153–167. [Google Scholar] [CrossRef]
  57. Andrews, M.G.; Massabò, R. Delamination in flat sheet geometries with material imperfections and thickness variations. Compos. Part B Eng. 2008, 39, 139–150. [Google Scholar] [CrossRef]
  58. Kardomateas, G.A.; Berggreen, C.; Carlsson, L.A. Energy-release rate and mode mixity of face/core debonds in sandwich beams. AIAA J. 2013, 51, 885–892. [Google Scholar] [CrossRef]
  59. Yu, H.-H.; He, M.Y.; Hutchinson, J.W. Edge effects in thin film delamination. Acta Mater. 2001, 49, 93–107. [Google Scholar] [CrossRef]
  60. Suo, Z.; Bao, G.; Fan, B.; Wang, T.C. Orthotropy rescaling and implications for fracture in composites. Int. J. Solids Struct. 1991, 28, 235–248. [Google Scholar] [CrossRef] [Green Version]
  61. Ustinov, K. On semi-infinite interface crack in bi-material elastic layer. Eur. J. Mech. A Solids 2019, 75, 56–69. [Google Scholar] [CrossRef]
  62. Massabó, R.; University of Genoa, Genova, Italy; Barbieri, L.; University of Genoa, Genova, Italy; Berggreen, C.; Technical University of Denmark, Lyngby, Denmark. The influence of shear on an interface crack between two elastic layers: Large mismatch of the elastic constants. Personal communication, 2019. [Google Scholar]
Figure 1. (a) An edge-delamination crack between a thin layer and a semi-infinite half-plane subject to edge loading per unit width, P , M , V , with M = M V a . (b) The crack tip element with elementary loads per unit width P , M , V . #1 and #2 indicate the materials of layer and substrate with α 0 and β 0 .
Figure 1. (a) An edge-delamination crack between a thin layer and a semi-infinite half-plane subject to edge loading per unit width, P , M , V , with M = M V a . (b) The crack tip element with elementary loads per unit width P , M , V . #1 and #2 indicate the materials of layer and substrate with α 0 and β 0 .
Coatings 09 00152 g001
Figure 2. The mode-mixity phase angle ψ P = ω in the edge-delamination crack in Figure 1 with P 0 ,   M = V = 0 , materials with β = 0 . See also Figure 4 for results plotted on a logarithmic scale.
Figure 2. The mode-mixity phase angle ψ P = ω in the edge-delamination crack in Figure 1 with P 0 ,   M = V = 0 , materials with β = 0 . See also Figure 4 for results plotted on a logarithmic scale.
Coatings 09 00152 g002
Figure 3. Fracture parameters in the edge-delamination crack in Figure 1 with V 0 ,   M = P = 0 , and materials with β = 0 ; (a) the mode-mixity phase angle ψ V ; (b) the function f V defining energy release rate, logarithmic scale; (c) the mode I stress intensity factor; (d) the mode II stress intensity factor.
Figure 3. Fracture parameters in the edge-delamination crack in Figure 1 with V 0 ,   M = P = 0 , and materials with β = 0 ; (a) the mode-mixity phase angle ψ V ; (b) the function f V defining energy release rate, logarithmic scale; (c) the mode I stress intensity factor; (d) the mode II stress intensity factor.
Coatings 09 00152 g003
Figure 4. (a) Mode mixity phase angle ψ P = ω as a function of Σ = E ¯ 1 / E ¯ 2 . The analytical results are compared with numerical results and with the asymptotics derived here and in [2]. (b) The zoom for large values of Σ = E ¯ 1 / E ¯ 2 . Range A in (a) is the interval where the asymptotics and expansions in Equations (14)–(16) do not ensure uncertainties below 0.2°.
Figure 4. (a) Mode mixity phase angle ψ P = ω as a function of Σ = E ¯ 1 / E ¯ 2 . The analytical results are compared with numerical results and with the asymptotics derived here and in [2]. (b) The zoom for large values of Σ = E ¯ 1 / E ¯ 2 . Range A in (a) is the interval where the asymptotics and expansions in Equations (14)–(16) do not ensure uncertainties below 0.2°.
Coatings 09 00152 g004
Figure 5. The mode mixity phase angle ψ V as a function of Σ = E ¯ 1 / E ¯ 2 . The analytical results are compared with numerical results and with the asymptotics and intermediate expansion in Equations (17)–(19). Range A in (a) is the interval where the asymptotics and expansions do not ensure uncertainties below 0.2°.
Figure 5. The mode mixity phase angle ψ V as a function of Σ = E ¯ 1 / E ¯ 2 . The analytical results are compared with numerical results and with the asymptotics and intermediate expansion in Equations (17)–(19). Range A in (a) is the interval where the asymptotics and expansions do not ensure uncertainties below 0.2°.
Coatings 09 00152 g005
Figure 6. The mode mixity phase angle ψ P = ω as a function of Σ = E ¯ 1 / E ¯ 2 . The analytical results for β = 0 are compared with numerical results for β 0 after Table 3 and [25]. The angle has been calculated at a reference distance r ^ = h 1 .
Figure 6. The mode mixity phase angle ψ P = ω as a function of Σ = E ¯ 1 / E ¯ 2 . The analytical results for β = 0 are compared with numerical results for β 0 after Table 3 and [25]. The angle has been calculated at a reference distance r ^ = h 1 .
Coatings 09 00152 g006
Table 1. Dundurs’ parameters for bi-material interfaces (Mat.#1/Mat.#2) with a large mismatch of the elastic constants (plane strain conditions) [2,3,8,11,33,34,35]; β is expressed here as a function of α and the Poisson ratios.
Table 1. Dundurs’ parameters for bi-material interfaces (Mat.#1/Mat.#2) with a large mismatch of the elastic constants (plane strain conditions) [2,3,8,11,33,34,35]; β is expressed here as a function of α and the Poisson ratios.
Dundurs’ Parameters α = E ¯ 1 E ¯ 2 E ¯ 1 + E ¯ 2 β = 1 4 [ ( ( 1 2 ν 2 ) ( 1 ν 2 ) + ( 1 2 ν 1 ) ( 1 ν 1 ) ) α + ( ( 1 2 ν 2 ) ( 1 ν 2 ) ( 1 2 ν 1 ) ( 1 ν 1 ) ) ]
Glass-epoxy/polymer foam0.990.28
Glass/epoxy0.940.19
Al/PMMA.0.910.23
ITO/PET.0.970.2
Ni/Polycarbonate0.9750.2
Al/Alfoam.0.820.23
Y2O3-ZrO2/β-(Ni,Pt)Al.0.630.31
Au/PDMS≈1≈0
SiO2/Epoxy0.890.21
SiC/SiO20.720.38
Table 2. Selected values of ω (in radians) and δ on varying α for materials with β = 0 .
Table 2. Selected values of ω (in radians) and δ on varying α for materials with β = 0 .
α −1−0.999−0.99−0.9−0.8−0.6−0.4−0.20
ω 0.7941830.7942680.7950420.8029550.8121570.8320390.8542950.8795920.908934
δ 0.4169840.4171190.4183350.4308850.4457190.4786930.5172080.5632410.619981
α 0.20.40.60.80.90.990.9980.9991
ω 0.9439430.9875181.0457551.1361371.2145481.3957601.4669321.48809π/2
δ 0.6930450.7936300.9490541.2586121.6370333.6765216.332027.98856
Table 3. Mode-mixity phase angle associated with the axial force, Figure 1a, with P 0 ,   M = V = 0 .
Table 3. Mode-mixity phase angle associated with the axial force, Figure 1a, with P 0 ,   M = V = 0 .
ψ p = ω = tan 1 ( Im [ K P h 1 i ε ] / Re [ K P h 1 i ε ] )   *
β = −0.4β = −0.3β = −0.2β = −0.1β = 0β = 0.1β = 0.2β = 0.3β = 0.4
α = −0.9964.0°59.0°54.3°49.8°45.2° (45.553°)
α = −0.964.7°59.8°55.2°50.6°46° (46.006°)
α = −0.864.4°60.2°55.5°51.0°46.7° (46.533°)
α = 0.8 65.7° (65.096°)62.0°57.0°54.9°50.3°
α = 0.9 69.6° (69.588°)66.1°62.6°58.9°55.1°
α = 0.99 80° (79.971°)77.7°75.4°73.1°70.7°
α = 0.998 84.1° (84.049°)82.7°81.4°80.1°78.8°
α = 0.999 85.3° (85.261°)84.3°83.4°82.6°81.8°
* The uncertainties are below ± 0.5 . Analytical results are shown in brackets. Results for α = 0.8 , 0.8 after [25].
Table 4. The mode-mixity phase angle associated with the elementary shear load, Figure 1a, with V 0 ,   P = M = 0 .
Table 4. The mode-mixity phase angle associated with the elementary shear load, Figure 1a, with V 0 ,   P = M = 0 .
ψ V = tan 1 ( Im [ K V h 1 i ε ] / Re [ K V h 1 i ε ] )   *
β = −0.4β = −0.3β = −0.2β = −0.1β = 0β = 0.1β = 0.2β = 0.3β = 0.4
α = −0.9918.4°15.0°11.7°8.6°5.4° (5.634°)
α = −0.917.9°14.6°11.4°8.3°5.2° (5.252°)
α = −0.817.2°14.0°10.9°7.8°4.8° (4.818°)
α = −0.6 12.7°9.7°6.8°3.9° (3.920°)
α = −0.4 11.3°8.5°5.7°3.0° (2.978°)0.2°
α = −0.2 7.2°4.6°2.0° (1.993°)−0.7°
α = 0 5.8°3.4°1.0° (0.963°)−1.5°−4.1°
α = 0.2 2.1°−0.1° (−0.107°)−2.4°−4.8°
α = 0.4 0.7°−1.2° (−1.208°)−3.2°−5.4°−7.6°
α = 0.6 −2.3° (−2.301°)−4.0°−5.8°−7.7°
α = 0.8 −3.2° (−3.238°)−4.4°−5.7°−7.0°−8.5°
α = 0.9 −3.4° (−3.427°)−4.1°−4.8°−5.7°−6.6°
α = 0.99 −2.3° (−2.284°)−1.5°−0.7°0.2°1.0°
α = 0.998 −1.3° (−1.436°)0.4°2.2°4.1°6.0°
* The uncertainties are below ± 0.5 . Analytical results are shown in brackets.
Table 5. The dimensionless function defining the energy release rate associated with the elementary shear load, Figure 1, with V 0 ,   P = M = 0 .
Table 5. The dimensionless function defining the energy release rate associated with the elementary shear load, Figure 1, with V 0 ,   P = M = 0 .
f V = E ¯ 1 h G V / V *
β = −0.4β = −0.3β = −0.2β = −0.1β = 0β = 0.1β = 0.2β = 0.3β = 0.4
α = −0.991.4441.4891.5281.5631.593 (1.597)
α = −0.91.4651.5101.5501.5851.616 (1.617)
α = −0.81.4871.5331.5731.6091.640 (1.641)
α = −0.6 1.5861.6271.6631.695 (1.696)
α = −0.4 1.6501.6921.7291.762 (1.762)1.791
α = −0.2 1.7731.8111.844 (1.845)1.874
α = 0 1.8781.9161.951 (1.951)1.9812.007
α = 0.2 2.0582.093 (2.093)2.1232.150
α = 0.4 2.2612.297 (2.298)2.3282.3552.378
α = 0.6 2.627 (2.628)2.6582.6852.708
α = 0.8 3.316 (3.317)3.3473.3743.3953.411
α = 0.9 4.190 (4.193)4.2224.2474.2664.279
α = 0.99 9.1 (9.089)9.19.19.19.1
α = 0.998 15.3 (15.559)15.615.615.515.5
* The uncertainties are below ± 0.007 ; ± 0.3 for α = 0.998 . Analytical results are shown in brackets.

Share and Cite

MDPI and ACS Style

Massabò, R.; Ustinov, K.; Barbieri, L.; Berggreen, C. Fracture Mechanics Solutions for Interfacial Cracks between Compressible Thin Layers and Substrates. Coatings 2019, 9, 152. https://doi.org/10.3390/coatings9030152

AMA Style

Massabò R, Ustinov K, Barbieri L, Berggreen C. Fracture Mechanics Solutions for Interfacial Cracks between Compressible Thin Layers and Substrates. Coatings. 2019; 9(3):152. https://doi.org/10.3390/coatings9030152

Chicago/Turabian Style

Massabò, Roberta, Konstantin Ustinov, Luca Barbieri, and Christian Berggreen. 2019. "Fracture Mechanics Solutions for Interfacial Cracks between Compressible Thin Layers and Substrates" Coatings 9, no. 3: 152. https://doi.org/10.3390/coatings9030152

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop