Next Article in Journal
A Review of Cu3BiS3 Thin Films: A Sustainable and Cost-Effective Photovoltaic Material
Previous Article in Journal
Molecular Dynamics Analysis of Collison Cascade in Graphite: Insights from Multiple Irradiation Scenarios at Low Temperature
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Three-Dimensional Singular Stress Fields and Interfacial Crack Path Instability in Bicrystalline Superlattices of Orthorhombic/Tetragonal Symmetries

by
Reaz A. Chaudhuri
Department of Materials Science & Engineering, 122 S. Central Campus Dr. Room 304, University of Utah, Salt Lake City, UT 84112-0560, USA
Retired Faculty.
Crystals 2024, 14(6), 523; https://doi.org/10.3390/cryst14060523
Submission received: 11 March 2024 / Revised: 24 May 2024 / Accepted: 25 May 2024 / Published: 30 May 2024
(This article belongs to the Section Crystal Engineering)

Abstract

:
First, a recently developed eigenfunction expansion technique, based in part on the separation of the thickness variable and partly utilizing a modified Frobenius-type series expansion technique in conjunction with the Eshelby–Stroh formalism, is employed to derive three-dimensional singular stress fields in the vicinity of the front of an interfacial crack weakening an infinite bicrystalline superlattice plate, made of orthorhombic (cubic, hexagonal, and tetragonal serving as special cases) phases of finite thickness and subjected to the far-field extension/bending, in-plane shear/twisting, and anti-plane shear loadings, distributed through the thickness. Crack-face boundary and interface contact conditions as well as those that are prescribed on the top and bottom surfaces of the bicrystalline superlattice plate are exactly satisfied. It also extends a recently developed concept of the lattice crack deflection (LCD) barrier to a superlattice, christened superlattice crack deflection (SCD) energy barrier for studying interfacial crack path instability, which can explain crack deflection from a difficult interface to an easier neighboring cleavage system. Additionally, the relationships of the nature (easy/easy, easy/difficult, or difficult/difficult) interfacial cleavage systems based on the present solutions with the structural chemistry aspects of the component phases (such as orthorhombic, tetragonal, hexagonal, as well as FCC (face-centered cubic) transition metals and perovskites) of the superlattice are also investigated. Finally, results pertaining to the through-thickness variations in mode I/II/III stress intensity factors and energy release rates for symmetric hyperbolic sine-distributed loads and their skew-symmetric counterparts that also satisfy the boundary conditions on the top and bottom surfaces of the bicrystalline superlattice plate under investigation also form an important part of the present investigation.

1. Introduction

Leading technological developments during the last decades relate to advanced materials, which are concerned with the deposition of thin films over substrates through the employment of techniques such as epitaxy, chemical vapor deposition (CVD), and physical vapor deposition (PVD) [1]. Bicrystals made of orthorhombic/cubic crystalline phases are common occurrences in many modern advanced technological applications, such as sensors [2], semiconductors [3], superconductivity [4], and so on. For example, gold nanocrystal superlattices can be formed on silicon nitride substrates with long-range ordering over several microns [2]. Pashley et al. [3] have reported the preparation of mono-crystalline films of gold and silver onto molybdenum disulfide inside an electron microscope that permits the direct observation of the mode of growth [1]. Although considered to be non-conventional for microelectronics, grain boundary junction engineering is frequently employed in metal oxide superconductor (MOS) THz frequency applications [4]. Yin et al. [5] employed the magnetron sputtering technique to deposit superconducting Yba2Cu3O7−δ (Yttrium barium copper oxide, or in short, YBCO) thin films on four polycrystalline metal substrates, with Yttrium-stabilized zirconia (YSZ) and silver serving as buffer layers.
A brief review of the literature reveals the extensive analytical treatment of singular stress fields at the tips of cracks, anticracks (through slit cracks filled with infinitely rigid lamellas), wedges, and junctions weakening/reinforcing homogeneous/biomaterial/trimaterial isotropic as well as anisotropic plates, from a two-dimensional point of view (see Refs. [6,7,8,9,10,11,12,13] and citations therein). In contrast, the mathematical difficulties posed by three-dimensional stress singularity problems have, until recently [14,15], posed an unsurmountable challenge to the researchers in the field. To start with, the governing PDE’s are much more complicated, and the exact satisfaction of the boundary and interface conditions is another matter altogether [16]. Earlier attempts to solve the three-dimensional through-crack problem resulted in controversies that lasted for about a quarter century [15,16]. More specifically, two-dimensional anisotropic fracture mechanics literature is more or less entirely based upon the inplane stress functions-based formulation due to Lekhnitskii [6], which is an extension of its isotropic counterpart due to Kolosov and Mushkhelishvili [17]. In contrast, a three-dimensional anisotropic fracture mechanics formulation cannot by definition employ the Lekhnitskii-pioneered inplane stress functions approach, and must be formulated in terms of displacement functions. More importantly, two-dimensional isotropic/anisotropic fracture mechanics studies, e.g., [8,9,10,11,12,13], are all based on complex variables. Since the three-dimensional space is too small to accommodate the next higher-dimensional analog of complex variables (for which at least a four-dimensional space will be required), these complex variable-based analyses are not by themselves adequate for the analysis of three-dimensional cracked anisotropic solids [18]. In addition, only an analytical solution can detail the structures of singularities related to the sharpness of a crack or anticrack, while the majority of weighted residual-type methods, e.g., the finite elements, finite difference, and boundary element approaches that are often beset by a lack of convergence, and oscillation (in the absence of the knowledge of the order of the stress singularity) resulting in poor accuracy, can hardly be expected to produce the necessary resolution [19,20]. The primary objective of the present investigation is to analytically solve the crack front stress singularity problems of bicrystalline superlattice plates made of orthorhombic, tetragonal, or cubic phases, subjected to mode I/II/III far-field loading, from a three-dimensional perspective.
During the last three decades, three-dimensional crack/anticrack/notch/antinotch/wedge front stress singularity problems have been solved by introducing a novel eigenfunction expansion technique. Various categories of three-dimensional stress singularity problems include: (i) through-thickness crack/anticrack [15,16,21] as well as their biomaterial and trimaterial interface counterparts [14,22,23], (ii) corresponding wedges/notches [24,25,26,27,28], (iii) biomaterial free/fixed straight edge-face [29,30,31], (iv) tri-material junction [32], (v) interfacial bond line of a tapered jointed plate [33], (vi) circumferential junction corner line of an island/substrate [1], (vii) fiber–matrix interfacial debond [34,35,36], (viii) fiber breaks and matrix cracking in composites [37], (v) penny-shaped crack/anticrack [38,39] and their biomaterial interface counterparts [40,41], (vi) through/part-through hole/rigid inclusion [42,43] and their biomaterial counterparts [44,45], as well as elastic inclusion [46,47], among others. Only the penny-shaped crack/anticrack [38] (and their biomaterial counterparts [40]), the hole [42], the biomaterial hole [44], and inclusion problems [46] have previously been adequately addressed in the literature. A unified three-dimensional eigenfunction approach has recently been developed by Chaudhuri and co-workers [1,14,15,16,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,41,43,45,47] to address three-dimensional stress singularity problems covering all the aspects mentioned above. These facts not only lend credence to the validity of the afore-mentioned three-dimensional eigenfunction expansion approach, but also reinforce the afore-mentioned conceptual as well as mathematical similarities of and linkages among the afore-cited classes of three-dimensional stress singularity problems.
The separation of the variables approach above has recently been extended to cubic and orthorhombic/orthotropic and monoclinic/anisotropic materials [48,49,50,51,52,53,54,55]. Cracked/anticracked transversely isotropic (smeared-out composite) [48] as well as cubic/orthorhombic/diamond cubic mono-crystalline plates subjected to mode I/II far-field loadings [49,50,51] and cubic/orthorhombic/monoclinic/diamond cubic mono/tricrystalline plates under mode III loading [49,51,52,53,54] have been solved by a novel three-dimensional eigenfunction expansion technique, based in part on the above-mentioned separation of the thickness variable and partly an affine transformation, which is similar (but not identical) in spirit to that due to Eshelby et al. [56] and Stroh [7]. This eigenfunction expansion approach has also recently been employed to obtain three-dimensional asymptotic stress fields in the vicinity of the front of the kinked carbon fiber–matrix junction [18] (see also Chaudhuri et al. [55] for its two-dimensional counterpart with isotropic glass fibers).
In what follows, the afore-mentioned modified eigenfunction expansion technique [18,48,49,50,51,52,53,54], based in part on the separation of the thickness variable and partly utilizing a modified Frobenius-type series expansion technique in conjunction with the Eshelby–Stroh formalism, is employed to derive heretofore unavailable three-dimensional singular stress fields in the vicinity of the front of an interfacial crack weakening an infinite bicrystalline superlattice plate, made of orthorhombic (cubic, hexagonal, and tetragonal serving as special cases) phases of finite thickness and subjected to far-field extension/bending, in-plane shear/twisting, and anti-plane shear loadings, distributed through the thickness. Crack-face boundary and interface contact conditions as well as those that are prescribed on the top and bottom surfaces of the bicrystalline superlattice plate are exactly satisfied. This has applications in semiconductors, photovoltaics, superconductors, and fiber-reinforced composites [18], among many others.
The second and more Important objective is to extend a recentlyIed concept of a lattice crack deflection (LCD) barrier [49,57] to a superlattice, christened the superlattice crack deflection (SCD) energy barrier, for studying interfacial crack path instability, which can explain crack deflection from a difficult interface to an easier neighboring cleavage system. What follows is an attempt to fill a gap or address an incompleteness regarding the concept of the sufficient condition for fracture in the contexts of three-dimensionality and a bicrystalline superlattice. That idea was previously explored for a cracked mono-crystalline lattice of orthorhombic symmetry only [57]. The second idea is the relationship of superlattice interfacial crack path stability/instability to the elastic parameter, which was not explored earlier. The present study is an attempt to tie up some of those loose and missing strings of ideas, and extend them to a bicrystalline superlattice. In other words, the present study combines two primary novelties: (i) three-dimensionality, and (ii) crack path stability/instability that satisfies the sufficient condition (non-Griffith/Irwin approach).
Third, the hitherto unexplored relationships of the nature (easy/easy, easy/difficult, or difficult/difficult) interfacial cleavage systems based on the present solutions with the structural chemistry aspects of the component phases (such as orthorhombic, tetragonal, hexagonal, as well as FCC (face-centered cubic) transition metals and perovskites) of the superlattice are also investigated. Finally, results pertaining to the through-thickness variations in mode I/II/III stress intensity factors and energy release rates for symmetric hyperbolic sine distributed loads and their skew-symmetric counterparts that also satisfy the boundary conditions on the top and bottom surfaces of the bicrystalline superlattice plate under investigation also form an important part of the present investigation.

2. Statement of the Problem

The Cartesian coordinate system (x, y, z) is convenient to describe the deformation behavior in the vicinity of the front of an interfacial semi-infinite crack weakening an infinite pie-shaped bicrystalline superlattice plate. The thickness of the bicrystalline superlattice plate is 2h (Figure 1). The cylindrical polar coordinate system (r, θ, z) is, however, more convenient to describe the boundary and interfacial contact conditions. The x-y plane serves as the interfacial plane. Here, the z-axis (|z|h) is placed along the straight crack front, while the coordinates r, θ, or x, y are used to define the position of a point in the plane of the plate (see Figure 1). The bicrystalline superlattice interface between materials 1 and 2 is located at θ = 0, i.e., it coincides with the positive x-axis, while the crack-side faces are located at θ = ±π (Figure 1). Components of the displacement vector in the jth crystal (j = 1, 2) in x and y directions are denoted by uj and vj, while the component in the z-direction is denoted by wj, j = 1, 2. The displacement components in the radial and tangential directions are represented by urj and uθj, respectively.
In what follows, σ x j , σ y j , and σ z j represent the normal stresses, while τ y z j , τ x z j , and τ x y j denote the shear stresses. ε x j , ε y j , and ε z j are the normal strains, while γ y z j , γ x z j , and γ x y j represent shear strains. c i l ( j ) , i, l = 1, …, 6, denotes the elastic stiffness constants of the jth orthorhombic crystal, given in the form:
{ σ x j σ y j σ z j τ y z j τ x z j τ x y j } = [ c 11 ( j ) c 12 ( j ) c 13 ( j ) 0 0 0 c 12 ( j ) c 22 ( j ) c 23 ( j ) 0 0 0 c 13 ( j ) c 23 ( j ) c 33 ( j ) 0 0 0 0 0 0 c 44 ( j ) 0 0 0 0 0 0 c 55 ( j ) 0 0 0 0 0 0 c 66 ( j ) ] { ε x j ε y j ε z j γ y z j γ x z j γ x y j }
In the absence of body forces, the coupled partial differential equations for a linear elastic crystalline solid in terms of the displacement functions uj, vj, and wj are given as follows:
c 11 ( j ) 2 u j x 2   +   c 66 ( j ) 2 u j y 2   +   c 55 ( j ) 2 u j z 2   +   ( c 12 ( j ) + c 66 ( j ) ) 2 v j x y +   ( c 13 ( j ) + c 55 ( j ) )   2 w j x z = 0 ,
( c 12 ( j ) + c 66 ( j ) ) 2 u j x y + c 66 ( j ) 2 v j x 2 +   c 22 ( j ) 2 v j y 2   +   c 44 ( j ) 2 v j z 2   +   ( c 23 ( j ) + c 44 ( j ) )   2 w j y z = 0 ,
( c 13 ( j ) + c 55 ( j ) ) 2 u j x z +   ( c 23 ( j ) + c 44 ( j ) )   2 v j y z + c 55 ( j ) 2 w j x 2 +   c 44 ( j ) 2 w j y 2   +   c 33 ( j ) 2 w j z 2   = 0 ,   j = 1 , 2
The boundary conditions on the bicrystalline superlattice plate surfaces, z = ±h, are given as follows:
Traction-free (see Section 6 below):
σ z j   =   τ x z j   =   τ y z j = 0 , j = 1 ,   2 .
The interfacial contact conditions as well as boundary conditions at the crack-side surfaces (valid for r > 0) are more conveniently expressed in local cylindrical polar coordinates (Figure 1). Assuming that the bicrystalline superlattice interface at θ = 0 is always perfectly bonded, it is easy to establish the continuity conditions of the stresses and displacements along this interface as follows:
( i )   θ = 0 :   u θ 1 =   u θ 2 ,   u r 1 =   u r 2 ,   u z 1 =   u z 2 ,
σ θ 1 =   σ θ 2 ,   τ r θ 1 =   τ r θ 2 ,   τ θ z 1 = τ θ z 2 .
The boundary conditions at the crack-side surfaces, at θ = ±π, are presented below:
σ θ 1 π = τ r θ 1 π = τ θ z 1 π = 0 ,
σ θ 2 π = τ r θ 2 π = τ θ z 2 π = 0 ,
where σrj, σθj, and σzj represent the normal stresses, and τrθj, τrzj, and τθzj denote the shear stresses in the jth crystal (j = 1, 2) in the cylindrical polar coordinate system (r, θ, z).

3. Singular Stress Fields in the Vicinity of a Crack Front Weakening a Bicrystalline Superlattice with Orthorhombic Phases under General Loading

The assumed displacement functions for the jth crystal (j = 1, 2) for the three-dimensional interfacial crack problem under consideration are selected on the basis of the separation of z-variables. These are presented below [18,48,49,50,51,52,53,54]:
u j ( x , y , z ) = e i α z U j ( x , y ) ,
v j ( x , y , z ) = e i α z V j ( x , y ) ,
w j ( x , y , z ) =   e i α z W j ( x , y ) ,
where α is a constant, called the wave number, required for a Fourier series expansion in the z-direction. It may be noted that, since the separated z-dependent term and its first partial derivative can either be bounded and integrable at most admitting ordinary discontinuities, or the first partial derivative at worst be a square integrable (in the sense of the Lebesgue integration) in its interval, z ∈ [−h, h], i.e., admitting singularities weaker than a square root (i.e., z (−1/2+ε), ε > 0 being a very small number), it can be best represented by the Fourier series [18,48,49,50,51,52,53,54]. The latter case is justified by the Parseval theorem [58], and its physical implication is that of satisfying the criterion of finiteness of local strain energy and path independence [59]. The substitution of Equation (6) into Equation (2) yields the following system of coupled partial differential equations (PDEs):
c 11 ( j ) 2 U j x 1 2   +   c 66 ( j ) 2 U j y 1 2   +   c 55 ( j ) U j   + ( c 12 ( j ) + c 66 ( j ) ) 2 V j x 1 y 1 +   ( c 13 ( j ) + c 55 ( j ) )   W j x 1 = 0 ,  
( c 12 ( j ) + c 66 ( j ) ) 2 U j x 1 y 1 + c 66 ( j ) 2 V j x 1 2   +   c 22 ( j ) 2 V j y 1 2   +   c 44 ( j ) V j   +   ( c 23 ( j ) + c 44 ( j ) )   W j y 1 = 0 ,
( c 13 ( j ) + c 55 ( j ) ) U j x 1 +   ( c 23 ( j ) + c 44 ( j ) )   V j y 1   + c 55 ( j ) 2 W j x 1 2 +   c 44 ( j ) 2 W j y 1 2   +   c 33 ( j ) W j   = 0 ,
where
x 1 = i α x ,       y 1 = i α y .
The solution to the system of coupled partial differential Equation (7) subjected to the most general loading can now be sought in the form of the following modified Frobenius-type series in terms of the variable x1 + py1 as follows:
U ( x 1 , y 1 ) = n = 0 a s + n ( x 1 + p y 1 ) s + 2 n + 1 + n = 0 a ¯ s + n ( x 1 + p y 1 ) s + 2 n ,
V ( x 1 , y 1 ) = n = 0 b s + n ( x 1 + p y 1 ) s + 2 n + 1 + n = 0 b ¯ s + n ( x 1 + p y 1 ) s + 2 n ,
W ( x 1 , y 1 ) = n = 0 c s + n ( x 1 + p y 1 ) s + 2 n + n = 0 c ¯ s + n ( x 1 + p y 1 ) s + 2 n + 1 .
For an orthorhombic lamina, out of the various combinations presented above, such as ( a s n , b s n , c s n ), ( a ¯ s n , b ¯ s n , c ¯ s n ), ( a s n , b ¯ s n , c ¯ s n ), ( a ¯ s n , b s n , c ¯ s n ), ( a ¯ s n , b ¯ s n , c s n ), ( a s n , b s n , c ¯ s n ), ( a s n , b ¯ s n , c s n ), and ( a ¯ s n , b s n , c s n ), only the first two groupings can produce meaningful solutions for mode III and mode I/II loading cases, respectively. This step permits the separation of mode III from modes I/II for the problem under investigation. The first grouping is described below, while the second one has previously been employed for the anti-plane shear case [49,51,52,53,54].

4. Singular Stress Fields in the Vicinity of a (010)[001] Through-Thickness Crack Front Propagating under Mode I (Extension/Bending) and Mode II (Sliding Shear/Twisting) in the [100] Direction

The solution to the system of coupled partial differential Equation (7), subjected to far-field mode I (extension/bending) and mode II (sliding shear/twisting) loading, can now be sought in the form of the following modified Frobenius-type series in terms of the variable x1 + py1 as follows [48,49,50,51], although unlike in the case of isotropic materials [14,15,16,21,22,23,24,25,26,29,30,31,32,33], the x1 and y1 variables are no longer separable:
U j ( x 1 , y 1 ) = n = 0 a ¯ s j + n ( x 1 + p ( j ) y 1 ) s + 2 n ,
V j ( x 1 , y 1 ) = n = 0 b ¯ s j + n ( x 1 + p ( j ) y 1 ) s + 2 n ,
W j ( x 1 , y 1 ) = n = 0 c ¯ s j + n ( x 1 + p ( j ) y 1 ) s + 2 n + 1 .
Here, s represents a fractional exponent, while n denotes an integer in the power series. The combined variable, x1 + py1, represents an affine transformation in the same spirit as that by Eshelby et al. [56] and Stroh [7], although these latter authors employed completely different techniques. The substitution of Equation (10) into Equation (7) and equating the coefficients of ( x 1 + p ( j ) y 1 ) s + 2 n 2 yields the following two recurrent relations:
( s + 2 n ) ( s + 2 n 1 ) ( c 11 + c 66 p ( j ) 2 ) a ¯ s + n + c 55 a ¯ s + n 1 + ( s + 2 n ) ( s + 2 n 1 ) ( c 12 + c 66 ) p ( j ) b s + n + ( s + 2 n 1 ) ( c 13 + c 55 ) c ¯ s + n 1 = 0 ,
( s + 2 n ) ( s + 2 n 1 ) ( c 12 + c 66 ) p ( j ) a ¯ s + n + ( s + 2 n ) ( s + 2 n 1 ) ( c 22 + c 66 p ( j ) 2 ) b ¯ s + n + c 44 b ¯ s + n 1 + ( s + 2 n 1 ) ( c 23 + c 44 ) c ¯ s + n 1 = 0 ,
which, for n = 0, reduces to
c 11 + c 66 p ( j ) 2 ( c 12 + c 66 ) p ( j ) ( c 12 + c 66 ) p ( j ) c 22 p ( j ) 2 + c 66 a ¯ s b ¯ s = 0 0 ,   for   s 0 ,   1 .
The characteristic equation for the coupled partial differential Equation (2) or Equation (7) can now be written as follows:
p ( j ) 4 + 2 χ ( j ) p ( j ) 2 + c 11 ( j ) c 22 ( j ) = 0 ,   j = 1 , 2
where the normalized elastic parameter, κ ( j ) = 1 / χ ( j ) , is given by
κ ( j ) = 1 χ ( j ) = 2 c 22 ( j ) c 66 ( j ) ( c 11 ( j ) c 22 ( j ) c 12 ( j ) 2 2 c 12 ( j ) c 66 ( j ) ) , j = 1 , 2 .
It is also sometimes convenient to relate κ = 1 / χ , as will be seen later, to A = 1 / λ , with A being the planar anisotropic ratio (in the x [100] − y [010] plane), as shown below.
κ ( j ) , j = 1, 2, can also be expressed in terms of the planar anisotropic ratio (in the x-y plane), A ( j ) , as follows:
κ ( j ) = A ( j ) c 22 ( j ) c 11 ( j ) c 22 ( j ) + c 12 ( j ) 1 A ( j ) , j = 1 , 2 ,
where
A ( j ) = 2 c 66 ( j ) c 11 ( j ) c 22 ( j ) c 12 ( j ) , j = 1 , 2 ,
It can easily be seen from Equation (15) that A is higher when the shear stiffness (modulus) and major Poisson’s ratio in the x [100], y [010] plane assume larger magnitudes. This simple fact assumes great importance as this investigation aims to solve one Holy Grail issue, in the fracture mechanics of anisotropic media, of coming up with a dimensionless parameter akin to the Reynold’s number in fluid flow problems, crossing a critical value that signifies a transition from one regime to another, such as the critical value of the Reynold’s number above which the flow is turbulent and below which it is laminar. It is an attendant issue relating to crack deflection in mono-crystalline orthorhombic laminas [57].
Equation (13) has either (a) four complex or (b) four imaginary roots, depending on whether:
( a )   A ( j ) >   1   or   equivalently ,   κ ( j ) > c 11 ( j ) c 22 ( j ) ,
or   ( b )   A ( j ) <   1   or   equivalently ,   κ ( j ) < c 11 ( j ) c 22 ( j ) ,
A ( j ) = 1 or κ ( j ) = 1 represents the degenerate isotropic material case, for which the solution is available in Chaudhuri and Xie [15,16].
Case (a): Complex Roots
p 1 , 2 = ξ ( j ) ± i η ( j ) , p 3 , 4 = ξ ( j ) ± i η ( j ) ,
where
ξ ( j ) = 1 2 c 11 ( j ) c 22 ( j ) 1 / 2 1 κ ( j ) 1 / 2 ,
η ( j ) = 1 2 c 11 ( j ) c 22 ( j ) 1 / 2 + 1 κ ( j ) 1 / 2 ,
valid for κ ( j ) > c 11 ( j ) c 22 ( j ) , .
Case (b): Imaginary Roots
The four imaginary roots of Equation (13) are given by
p 1 , 2 = ± i ξ ( j ) + i η ( j ) , p 3 , 4 = ± i ξ ( j ) i η ( j ) ,
where
ξ ( j ) = 1 2 c 11 ( j ) c 22 ( j ) 1 / 2 + 1 κ ( j ) 1 / 2 ,
ξ ( j ) = 1 2 c 11 ( j ) c 22 ( j ) 1 / 2 + 1 κ ( j ) 1 / 2 ,
valid for κ ( j ) < c 11 ( j ) c 22 ( j ) ,

5. Satisfaction of Crack-Face Boundary and Interfacial Contact Conditions

5.1. Both Crystal Layers with Complex Roots

For the cases of both crystal layers with complex roots, the substitution of Equations (A2b,c) and (A6a,b) in conjunction with Equation (A7) into the left- and right-hand sides of Equations (4a) and (4b), respectively, would yield four homogeneous equations. In addition, the substitution of Equation (A2b,c), in conjunction with Equation (A7), into the left-hand sides of Equation (5a,b) would yield another set of four homogeneous equations. These can be expressed in the compact form as follows:
Δ s A i j = 0 .
The existence of a nontrivial solution for A i j requires a vanishing of the coefficient determinant
Δ s = 0 .
Δ s is an 8×8 matrix involving s in a transcendental form. The solution that has physical meaning is given by 0 < Re(s) < 1. Here, s = 0.5 ± i ε , where ε is obtained from the following relationships, obtained by equating the real and imaginary parts:
Real Part:
a cosh 4 ε π b cosh 2 ε π + c = 0 ,
Equation (22) yields
cosh ε π = b ± b 2 4 a c 2 a 1 / 2 ,
where
a = a 1 e 1 + b 1 f 1 c 2 2 d 1 f 1 ,
b = c 2 f 1 a 1 c 2 a 1 c 1 + b 1 d 1 c 1 f 1 d 1 e 1 b 2 d 1 + d 2 + a 2 c 2 ,
c = a 1 b 2 c 2 d 1 c 1 f 1 d 1 e 1 .
where
a 1 = c 12 ( 1 ) + c 22 ( 1 ) H 11 ξ ( 1 ) H 21 η ( 1 ) ,
b 1 = c 22 ( 1 ) H 11 ξ ( 1 ) + H 21 η ( 1 ) ,
c 1 = H 11 + ξ ( 1 ) ,
d 1 = H 21 + η ( 1 ) ,
a 2 = c 12 ( 2 ) + c 22 ( 2 ) H 12 ξ ( 2 ) H 22 η ( 2 ) ,
b 2 = c 22 ( 2 ) H 22 ξ ( 2 ) + H 12 η ( 2 ) ,
c 2 = H 12 + ξ ( 2 ) ,
d 2 = H 22 + η ( 2 ) ,
e 1 = H 11 a 2 c 2 + b 2 d 2 c 1 H 12 a 2 H 22 b 2 H 22 c 2 + H 12 d 2 ,
f 1 = H 21 a 2 c 2 + b 2 d 2 d 1 H 12 a 2 H 22 b 2 H 22 c 2 + H 12 d 2 .
Imaginary Part:
a cosh 4 ε π + b cosh 2 ε π + c = 0 ,
where
a = c + d e ,   b = e 2 c .
Equation (28) yields
cosh δ π = e 2 c ± e 2 4 d c 2 d e + c 1 / 2 ,
where
c = a 1 b 2 2 d 1 c 1 f 1 d 1 e 1 ,
d = c 2 f 1 a 1 e 1 + b 1 f 1 b 2 d 2 c 2 a 1 a 2 ,
e = a 1 b 2 c 2 d 1 a 1 e 1 + b 1 f 1 + c 1 f 1 d 1 e 1 a 2 b 2 c 2 f 1 b 2 2 d 2 f 1 a 1 2 b 2 c 2 .

5.2. Both Crystal Layers with Imaginary Roots

Similarly, for the case of both crystal layers with imaginary roots, the substitution of Equations (A9b,c) and (A12a,b) in conjunction with Equation (A7) into the left- and right-hand sides of Equations (4a) and (4b), respectively, would also yield four homogeneous equations. Additionally, the substitution of Equation (A9b,c), in conjunction with Equation (A7), into the left-hand sides of (5a,b) would yield another set of four homogeneous equations. Proceeding in a similar manner, one obtains s = 0.5 ± i ε π , where ε is obtained from the relationship:
a cosh 4 ε π + b cosh 2 ε π + c = 0 .
Equation (32) yields
cosh ε π = b ± b 2 4 a c 2 a 1 / 2 ,
where
a = f 1 f 2 2 a 1 b 1 c 1 d 1 + a 1 2 d 1 2 + e 1 2 + 2 a 1 2 d 1 e 1 e 2 f 1 + a 1 c 1 e 2 a 2 b 2 e 1 2 + b 1 2 c 1 2 + a 1 2 d 1 2 ,
b = f 1 f 2 2 a 1 b 1 c 1 d 1 a 1 2 d 1 2 e 1 2 + a 1 e 1 e 2 f 1 b 1 c 1 2 a 1 d 1 + a 1 c 1 a 1 d 1 b 1 c 1 f 1 g 2 f 2 g 1 + a 1 c 1 a 2 b 2 c 1 e 1 g 2 d 1 e 1 f 2 + b 1 2 c 1 2 e 2 2 e 1 2 e 2 a 1 2 d 1 2 e 2 ,
c = a 1 b 1 c 1 d 1 f 1 f 2 + a 1 c 1 e 1 e 2 a 1 g 1 b 1 f 1 + a 1 c 1 a 1 d 1 + b 1 c 1 f 1 g 2 + f 2 g 1 a 1 2 c 1 2 g 1 g 2 + a 1 c 1 e 1 a 2 b 2 e 1 e 2 c 1 g 2 + d 1 f 2 ,
where
a 1 = c 12 ( 1 ) + c 22 ( 1 ) H 11 ξ 1 ( 1 ) + η 1 ( 1 ) ,
b 1 = c 12 ( 1 ) + c 22 ( 1 ) H 21 ξ 1 ( 1 ) η 1 ( 1 ) ,
c 1 = H 11 ξ 1 ( 1 ) + η 1 ( 1 ) ,
d 1 = H 21 ξ 1 ( 1 ) η 1 ( 1 ) ,
a 2 = c 12 ( 2 ) + c 22 ( 2 ) H 12 ξ 1 ( 2 ) + η 1 ( 2 ) ,
b 2 = c 12 ( 2 ) + c 22 ( 2 ) H 22 ξ 1 ( 2 ) η 1 ( 2 ) ,
c 2 = H 12 ξ 1 ( 2 ) + η 1 ( 2 ) ,
d 2 = H 22 ξ 1 ( 2 ) η 1 ( 2 ) ,
e 1 = b 1 c 1 a 1 d 1 ,
f 1 = c 2 a 1 b 2 d 2 a 1 a 2 ,
g 1 = c 2 b 1 b 2 d 2 b 1 a 2 ,
e 2 = H 22 c 2 H 12 d 2 ,
f 2 = b 2 H 12 c 1 H 11 c 2 H 22 a 1 c 1 + H 11 a 2 d 2 ,
g 2 = H 21 b 2 c 2 a 2 d 2 + d 1 H 22 a 2 H 12 b 2 .

5.3. Top Crystal (Layer 1) with Complex Roots and Bottom Crystal (Layer 2) with Imaginary Roots

For the case of the top crystal (layer 1) with complex roots and bottom crystal (layer 2) with imaginary roots, the substitutions of Equation (A6a,b) as well as Equation (A2b,c) in conjunction with Equation (A7) into the left-, and those of Equation (A12a,b) as well as Equation (A9b,c) in conjunction with Equation (A7) into the right-hand sides of Equation (4a,b) would yield four homogeneous equations. In addition, the substitution of Equation (A2b,c), in conjunction with Equation (A7), into the left-hand sides of Equation (5a) would yield another set of two homogeneous equations. Likewise, the substitution of Equation (A9b,c), in conjunction with Equation (A7), into the left-hand sides of (5b) would yield the remaining set of two homogeneous equations. Proceeding in a similar manner as above, one obtains s = 0.5 ± i ε ¯ π , where ε ¯ is obtained from the relationship:
a ¯ cosh 4 ε ¯ π + b ¯ cosh 2 ε ¯ π + c ¯ = 0 ,
Equation (39) yields
cosh ε ¯ π = b ¯ ± b ¯ 2 4 a ¯ c ¯ 2 a ¯ 1 / 2 ,
where
a ¯ = d ¯ e ¯ g ¯ ,   b ¯ = d ¯ e ¯ h ¯ f ¯ g ¯ ,   c ¯ = f ¯ h ¯ ,
where
d ¯ = e ¯ 1 2 H 11 f ¯ 1 + c 1 e ¯ 2 + b 1 f ¯ 2 a 1 c 2 d 2 d 1 a 2 b 2 f ¯ 1 ,
e ¯ = e ¯ 1 H 21 f ¯ 1 + d 1 e ¯ 2 a 1 f ¯ 2 ,
f ¯ = d 1 e ¯ 1 f ¯ 1 + e ¯ 2 ,
g ¯ = b 1 e ¯ 1 c 2 d 2 ,
h ¯ = b 1 c 1 f ¯ 1 + e ¯ 1 c 2 d 2 ,
where
e ¯ 1 = a 1 c 1 + b 1 d 1 ,
f ¯ 1 = b 2 c 2 a 2 d 2 ,
e ¯ 2 = a 2 H 22 b 2 H 12 ,
f ¯ 2 = c 2 H 22 d 2 H 12 .
It may be noted that, for a homogeneous solution, such as what is presented above, the seven eigenvector components, A 12 z , A 21 z , …, A 42 z , cannot be evaluated in absolute terms, but can only be expressed in terms of the eighth component, A 11 z , which can be related to the “stress intensity factors”, K I j z and K I I j z , which can, in turn, be resolved as follows:
K I j z = K I s j z + K I j a z
K I I j z = K I I s j z + K I I j a z
where K I s j z and K I I s j z , j = 1, 2, represent the symmetric (with respect to z) mode I and mode II stress intensity factors, which correspond to in-plane extension and sliding shear, respectively. K I a j z and K I I a j z , j = 1, 2, denote the antisymmetric (with respect to z) mode I and mode II stress intensity factors, which correspond to bending and twisting, respectively.
The stress intensity factors and the energy release rates are found to be as follows:
K I / I I j z = K I / I I j , 2 D D b j z , G I / I I j z = G I / I I j , 2 D D b j z 2 ,
where K I / I I j , 2 D and G I / I I j , 2 D are available in the published literature; see e.g., Wu [10] and Wang et al. [12]. In what follows, the primary focus will be on the three-dimensionality and lattice (more specifically super-lattice) aspects of bicrystalline interfacial fracture.

6. Boundary Conditions on the Bicrystalline Superlattice Plate Surfaces and Through-Thickness Distribution of Singular Stress Fields

6.1. General Distributed Far-Field Loading

By using the traction-free boundary conditions at the top and bottom surfaces, z* = z/h = ±1, of the plate, given by Equation (3), the general form of D b j z , defined in Equation (A3), can be obtained. The stress field that satisfies traction-free boundary conditions at z* = ±1, in the vicinity of the front of a bicrystalline interface crack under extension/sliding shear, can be recovered if, in Equation (A3):
D b j z   =   D b s   j z = D 2 j cos α z , j = 1 , 2
By substituting Equation (A2d) or Equation (A9d) into the boundary condition on the singular stress-free surface of a cracked bicrystalline superlattice plate, given by Equation (3), the general form of D b s   j z can be obtained as follows:
  D b s j z = m = 0 ± D 1 m j cos ( 2 m + 1 ) π z * 2 , j = 1 , 2 .
Hence, KI = KIs and KII = KIIs represent symmetric “stress intensity factors”. If the odd functions are selected from D b j z that satisfies traction-free boundary condition at z* = ± 1, it can yield the out-of-plane bending/twisting case given by
D b j z =   D b a j z = D 1 j sin α z , j = 1 , 2
D b a   j z that satisfies the singular stress-free conditions on the cracked bicrystalline superlattice plate surfaces is given by
D b a j z = m = 1 ± D 2 m j sin m π z * , j = 1 , 2 ,
Here, KI = KIa and KII = KIIa are anti-symmetric stress intensity factors. Equations (48) and (49) are valid, provided the loading function vanishes at z* = 0, thus eliminating the possibility of discontinuity of the function at z* = 0. In the presence of the discontinuity of the function at z* = 0, D b a j z can be written as follows:
D b a j z = m = 0 ± D 1 m j cos ( 2 m + 1 ) π z * 2 , j = 1 ,   2
If the singular stress-free boundary conditions, given by Equation (3), on the cracked bicrystalline superlattice plate surfaces are satisfied, all the displacements also vanish on these surfaces in the vicinity of the front of a bicrystalline interface crack.

6.2. Hyperbolic Sine Distributed Far-Field Loading

Hyperbolic sine distributed far-field loading, which is proportional to sinh ( z * ) , | z * | < 1, is applied. The applied antisymmetric loading function and the corresponding stress intensity factors (valid for | z ˜ * | ≤ 1) are proportional to
D b a z * = sinh ( z * ) = exp ( z * ) exp ( z * ) 2 .
The corresponding Fourier series can be derived as follows:
D b a z * = m = 0 e 1 1 m e 1 m 1 + m 2 π 2 m π sin m π z * .
The applied symmetric loading function (valid for | z * | < 1) and the corresponding stress intensity factors (valid for | z * | ≤ 1) are proportional to
D b s z * = sinh ( z ) = 1 2 exp ( z * ) exp ( z * ) .
The corresponding Fourier series can be obtained as given below:
D b s z * = m = 0 e 1 1 m e 1 m 1 + m 2 π 2 m π sin m π z * .
Alternatively, since D b s z * has no discontinuity at z * = 0, it can also be derived as given below:
D b s z * = m = 0 e 1 1 m e 1 m 1 + m + 1 2 2 π 2 m + 1 2 π cos 2 m + 1 2 π z * .

7. Singular Stress Fields in the Vicinity of a Through-Crack Front Propagating under Mode III (Anti-Plane Shear) in the [100] Direction

Yoon and Chaudhuri [54] have already addressed the problem of a cracked tricrystal plate subjected to mode III (anti-plane shear loading). The present problem is a special case of that, and can be obtained by substituting θ1 = π (refer to Figure 1 of Yoon and Chaudhuri [54]).

8. Crack Path Stability/Instability Criteria

8.1. Necessary Condition—Griffith–Irwin Theory-Based Crack Deflection Criterion

The important issue of a cleavage plane being deemed easy or difficult can be related to a crack deflection criterion, which is based on the relative fracture energy (or the energy release rate) available for possible fracture pathways [13]. The deflection or kinking of a crack from cleavage system 1 to cleavage system 2 is favored if (however, not iff, i.e., if and only if):
G 1 2 Γ 1 < 1 < G 2 2 Γ 2 G 2 G 1 > Γ 2 Γ 1 ,
where Gj and Γj, j = 1, 2, are the energy release rate (fracture energy) and surface energy, respectively, of the jth cleavage system.
It may be noted here that an experimental determination of the critical surface energy, Γj, j = 1, 2, of the component phase or the corresponding interfacial energy, Γint, of a bicrystalline superlattice can sometimes be notoriously challenging, due to the presence of micro-to-nanoscale defects, such as porosity, dislocation, twin boundaries, misalignment of bonds with respect to the loading axis, and the like. In contrast, the bond shear strain at superlattice crack deflection, γ b d s , and superlattice crack deflection (SCD) barrier, Δ K s * , to be discussed below, are, relatively speaking, much easier to determine in comparison to the critical surface or interfacial energy.

8.2. Sufficient Condition

The above Griffith–Irwin theory-based crack deflection criterion (albeit being still very useful and widely employed) is not accepted as a sufficient condition for a cleavage system deemed to be easy or difficult for crack propagation in single crystals [60], the most enigmatic being the {100} cleavage of BCC transition metals, especially macroscopically isotropic tungsten (W). Although the close-packed {110} surface has a lower surface energy, the preferred cleavage plane has experimentally been observed by Hull et al. [61] to be {100} (see also Riedle et al. [60]). For W, {100}〈001〉 ({crack plane}<crack front>) is deemed to be the preferred cleavage system for the reason of having d2sp3 hybridized orbitals, while {−110}<001> is considered to be difficult for crack propagation.
Atomistic scale modeling of cracks requires the consideration of both long-range elastic interactions and short-range chemical reactions. The Griffith theory does not take the latter into account [49]. Secondly, and more importantly, fracture criteria derived from equilibrium theories, such as the Griffith (thermodynamics-based) energy balance criterion, are not equipped to meet the sufficiency condition because of prevailing non-equilibrium conditions, such as physico-chemical reactions during crack propagation. Hence, such criteria can only be regarded as necessary conditions for fractures, but not as sufficient [60,62]. The effect of short-range chemical reactions can obviously be encapsulated by atomic-scale simulations, such as the investigation of low-speed propagation instabilities in silicon using quantum–mechanical hybrid multi-scale modeling due to Kermode et al. [63], which, however, entails extensive computational and other resources. Alternatively, and more importantly, such short-range interactions can also be captured by elastic property-based parameters (with a few exceptions), such as the planar anisotropic ratio, Aj, or equivalently, the normalized elastic parameter, κj, j = 1, 2 [49]. The general theory behind these characteristics pertaining to the structural chemistry of crystals is available in well-known treatises (see, e.g., [64,65,66]). More specifically, the elastic properties of superconducting YBa2Cu3O7−δ are strongly influenced by oxygen non-stoichiometry (as well as various structural defects).
A single dimensionless parameter, such as the planar anisotropic ratio, Aj, or equivalently, the normalized elastic parameter, κj, j = 1, 2, can serve as the Holy Grail quantity for an a priori determination of the status of a cleavage system to be easy or difficult, very much akin to the Reynold’s number for fluid flow problems, crossing a critical value that signifies the transition from one regime to another. Here, the planar anisotropic ratio, Aj, or equivalently, the normalized elastic parameter, κj, j = 1, 2, for a (010)[001]×[100] cleavage system, crossing the critical value of 1 or c 22 ( j ) / c 11 ( j ) , j = 1, 2, respectively, signifies transitioning from self-similar crack growth or propagation to crack deflection or turning from a difficult cleavage system onto a nearby easy one. This is a significant qualitative as well as quantitative improvement over a two-parameter-based model, introduced by Wang et al. [12], in the context of two-dimensional anisotropic fracture mechanics.
Finally, just as the introduction of the Reynold’s number facilitated the design and setting up of experiments in addition to the experimental verification of analytical and computational solutions in fluid dynamics, the accuracy and efficacy of the available experimental results on elastic constants of mono-crystalline superconducting YBa2Cu3O7−δ, measured by modern experimental techniques with resolutions at the atomic scale or nearly so, such as X-ray diffraction [67], ultrasound technique [68,69,70], neutron diffraction [71]/scattering [72], Brillouin spectroscopy [73,74]/scattering [75], resonant ultrasound spectroscopy [76,77], and the like, in a way best suited to preserve the characteristics associated with short-range interactions [57], are assessed with a powerful theoretical analysis of crack path stability/instability, in part based on a single dimensionless parameter, such as the planar anisotropic ratio, Aj, or equivalently, the normalized elastic parameter, κj, j = 1, 2.

9. Numerical Results and Discussions (Table 1, Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8, Table 9, Table 10, Table 11, Table 12, Table 13, Table 14, Table 15, Table 16, Table 17, Table 18, Table 19, Table 20, Table 21, Table 22, Table 23, Table 24, Table 25, Table 26, Table 27, Table 28, Table 29, Table 30, Table 31, Table 32, Table 33, Table 34, Table 35, Table 36, Table 37, Table 38, Table 39, Table 40, Table 41, Table 42, Table 43, Table 44, Table 45, Table 46, Table 47, Table 48 and Table 49)

9.1. Structure–Fracture Property Relations for Certain Model Bicrystalline Superlattices

In what follows, 22 model bicrystal (superlattice) cleavage systems, each comprising a nano-film deposited on a substrate, are investigated. Table 1 lists the structures and elastic stiffness constants (with respect to <100> axes) of these mono-crystalline materials. Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8, Table 9, Table 10, Table 11, Table 12, Table 13, Table 14, Table 15, Table 16, Table 17, Table 18, Table 19, Table 20, Table 21, Table 22 and Table 23 list the cleavage systems, (crack plane)[crack front]×[initial propagation direction], and elastic stiffness constants (with respect to suitably rotated coordinates) of the component phases of bicrystalline superlattice systems 1–22. For example, the hexagonal substrate or material 2 is rotated such that {11 2 ¯ 0} (a prism plane) is parallel to the crack face. It may be noted that the rotated crystal displays tetragonal-type symmetry.
The bicrystalline superlattice systems investigated here are namely as follows: (i) Au (gold), nano-layer/film or material 1 deposited on Si3N4 (silicon nitride), substrate or material 2 (Table 2, Table 3, Table 4 and Table 5); (ii) Au (nano-layer/film) deposited on substrate MgO (magnesium oxide) (Table 6, Table 7, Table 8 and Table 9); (iii) YBa2C3O7 (tetragonal/fully oxidized or non-superconducting YBa2Cu3O7, in short YBCOT) nano-layer/film, deposited on substrate Si3N4 (Table 10); (iv) YBa2C3O7 (nano-layer/film) deposited on substrate SrTiO3 (Table 11, Table 12, Table 13 and Table 14); (v) YBa2C3O7−δ (superconducting YBCO, in short YBCO) nano-layer/film, deposited on substrate Si3N4 (Table 15); (vi) YBa2C3O7−δ, nano-layer/film deposited on substrate MgO (Table 16 and Table 17); and (vii) YBa2C3O7−δ (nano-layer/film) deposited on substrate SrTiO3 (strontium titanate) (Table 18, Table 19, Table 20, Table 21, Table 22 and Table 23). For the bicrystalline superlattice systems under investigation, the computed mode I, II, or mixed-mode I/II order-of-stress singularity, λ i = 1 − s i , i = 1, 2, is found to be 0.5 ± iε or 0.5. In contrast, the computed mode III order-of-stress singularity, λ 3 = 1 – s3, is always equal to 0.5.
Table 1. Structures and elastic properties of various single crystals.
Table 1. Structures and elastic properties of various single crystals.
Single CrystalBravais LatticeStructure c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
Au [78]FCCFCC192.9192.9192.9163.8163.8163.841.541.541.5
MgO [78]FCCRock Salt289.3289.3289.387.7087.7087.70154.77154.77154.77
SrTiO3 [78]Simple CubicPerovskite348.17348.17348.17100.64100.64100.64454.55454.55454.55
Si3N4HCPHCP343.0343.0600.0136.0120.0120.0124.0124.0103.5
YBa2C3O7 [72] *TetragonalPerovskite230.0230.0150.0100.0100.0100.050.050.085.0
YBa2C3O7−δ [57,70,76] OrthorhombicPerovskite231.0268.0186.066.071.095.049.037.082.0
* All values measured by neutron scattering by Reichardt et al. [72]; all values measured by resonant ultrasound spectroscopy by Lei et al. [76], except c 12 and c 66 measured by ultrasound by Saint-Paul and Henry [70].
Table 2. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 1.
Table 2. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 1.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1Au (FCC)(010)[001]×[100]192.9192.9192.9163.8163.8163.841.541.541.5
2 *Si3N4 (HCP)(001)[0 1 ¯ 0]×[ 1 ¯ 00]343.0600.0343.0120.0136.0120.0124.0103.5124.0
* Rotated about the z-axis by 90°, then rotated about the new x-axis by −90°.
Table 3. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 2.
Table 3. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 2.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1Au (FCC)(010)[001]×[100]192.9192.9192.9163.8163.8163.841.541.541.5
2 *Si3N4 (HCP)( 1 ¯ 00)[0 1 ¯ 0]×[001]600.0343.0343.0120.0120.0136.0103.5124.0124.0
* Rotated first about the y-axis by −90°, then rotated about the new x-axis by 90°.
Table 4. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 3.
Table 4. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 3.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1Au (FCC)( 1 ¯ 10)[001]×[110]219.85219.85192.9136.85163.8163.841.541.514.5
2 *Si3N4 (HCP)(001)[0 1 ¯ 0]×[ 1 ¯ 00]343.0600.0343.0120.0136.0120.0124.0103.5124.0
* Rotated about the z-axis by 90°, then rotated about the new x-axis by −90°.
Table 5. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 4.
Table 5. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 4.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1Au (FCC)( 1 ¯ 10)[001]×[110]219.85219.85192.9136.85163.8163.841.541.514.5
2 *Si3N4 (HCP)( 1 ¯ 00)[0 1 ¯ 0]×[001]600.0343.0343.0120.0120.0136.0103.5124.0124.0
* Rotated first about the y-axis by −90°, then rotated about the new x-axis by 90°.
Table 6. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 5.
Table 6. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 5.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1Au (FCC)(010)[001]×[100]192.9192.9192.9163.8163.8163.841.541.541.5
2MgO(FCC)(010)[001]×[100]289.3289.3289.387.787.787.7154.77154.77154.77
Table 7. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 6.
Table 7. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 6.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1Au (FCC)(010)[001]×[100]192.9192.9192.9163.8163.8163.841.541.541.5
2MgO (FCC)( 1 ¯ 10)[001]×[110]343.27343.27289.333.7387.787.7154.77154.77100.8
Table 8. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 7.
Table 8. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 7.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1Au (FCC)( 1 ¯ 10)[001]×[110]219.85219.85192.9136.85163.8163.841.541.514.5
2MgO (FCC)(010)[001]×[100]289.3289.3289.387.787.787.7154.77154.77154.77
Table 9. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 8.
Table 9. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 8.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1Au (FCC)( 1 ¯ 10)[001]×[110]219.85219.85192.9136.85163.8163.841.541.514.5
2MgO (FCC)( 1 ¯ 10)[001]×[110]343.27343.27289.333.7387.787.7154.77154.77100.8
In the FCC metal nano-film, Au, listed in Table 2, Table 3, Table 4, Table 5, Table 6, Table 7, Table 8, Table 9, Table 24, Table 25, Table 26, Table 27, Table 28, Table 29, Table 30 and Table 31, all the orbitals belong to the d-block with partially filled d-shells [66]. In a mono-crystalline FCC metal, the bonds are oriented along the face diagonals, <110>. Such a metal contains linear chains of near-neighbor bonds in these directions, resulting in higher elastic stiffness constants along them. As shown in Table 20, Table 25, Table 28, and Table 29, A = 2.8522 > 1 and κ = 4.9777 > c 22 / c 11 = 1 for Au, giving rise to complex roots for a {010}<001>×<100> through-crack. Similar calculations yield A = 0.3494 < 1 and κ = 0.2487 < c 22 / c 11 = 1 for Au, giving rise to imaginary roots for the { 1 ¯ 10 }<001>×<110> through-crack [49], as can be seen in Table 26, Table 27, Table 30, and Table 31. It can then be inferred that 1 ¯ 10 <001>×<110> would constitute an easy cleavage system, while {010}<001>×<100> would be deemed difficult.
Table 10. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 9.
Table 10. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 9.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1YBa2C3O7 (Tetragonal)(010)[001]×[100]230.0230.0150.0100.0100.0100.050.050.085.0
2 *Si3N4 (HCP)( 1 ¯ 00)[0 1 ¯ 0]×[001]600.0343.0343.0120.0120.0136.0103.5124.0124.0
* Rotated first about the y-axis by −90°, then rotated about the new x-axis by 90°.
Table 11. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 10.
Table 11. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 10.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1YBa2C3O7 (Tetragonal)(010)[001]×[100]230.0230.0150.0100.0100.0100.050.050.085.0
2SrTiO3 (Simple Cubic) (010)[001]×[100]348.17348.17348.17100.64100.64100.64454.55454.55454.55
Table 12. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 11.
Table 12. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 11.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1YBa2C3O7 (Tetragonal)(010)[001]×[100]230.0230.0150.0100.0100.0100.050.050.085.0
2SrTiO3 (Simple Cubic) ( 1 ¯ 10)[001]×[110]678.96678.96348.17−230.15100.64100.64454.55454.55123.77
Table 13. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 12.
Table 13. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 12.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1YBa2C3O7 (Tetragonal)( 1 ¯ 10)[001]×[110]250.0250.0150.080.0100.0100.050.050.065.0
2SrTiO3 (Simple Cubic) ( 1 ¯ 10)[001]×[110]678.96678.96348.17−230.15100.64100.64454.55454.55123.77
Table 14. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 13.
Table 14. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 13.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1YBa2C3O7 (Tetragonal)( 1 ¯ 10)[001]×[110]250.0250.0150.080.0100.0100.050.050.065.0
2SrTiO3 (Simple Cubic)(010)[001]×[100]348.17348.17348.17100.64100.64100.64454.55454.55454.55
The ionic crystal MgO (an alkaline earth metal oxide), used as a substrate, which is listed in Table 6, Table 7, Table 8, Table 9, Table 16 and Table 17, in addition to in Table 28, Table 29, Table 30, Table 31, Table 38 and Table 39, is structurally of the rock salt type, but is an exception to the general rule for ionic crystals, such as alkali halides (e.g., NaCl and KCl) with a rock salt structure [49,64]. The reason is, as explained by Newnham [64], due to Mg2+ (and also Li+) having small cations, which permit the anions to be in contact with one another and, consequently, restrict bending actions. As a result, elastic stiffness coefficients in the <110> and <111> directions become larger than their <100> counterparts. This is in contrast to NaCl and KCl, wherein Cl− anions are not in contact. Additionally, the importance of anion–anion forces were pointed out by Weidner and Simmons [79]. These researchers found, in connection with the computation of elastic properties of several alkali halides from a two-body central force model, the necessity to include anion–anion interactions in addition to cation–anion forces. As shown in Table 29, Table 31, and Table 39, A = 0.6513 and κ = 0.6297 are both less than unity ( A < 1, κ < c 22 / c 11 = 1), giving rise to imaginary roots for the ( 1 ¯ 10)[001]×[110] through-crack. This is in contrast to A = 1.5354 and κ = 1.8329 being both larger than unity (A > 1, κ > c 22 / c 11 = 1), giving rise to complex roots for the (010)[001]×[100] through-crack, as shown in Table 28, Table 30, and Table 38. It can then be inferred that 1 ¯ 10 <001>×<110> would constitute an easy cleavage system, while {010}<001>x<100> would be deemed difficult.
Table 15. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 14.
Table 15. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 14.
Mater-ial (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1YBa2C3O7−δ (Orthorhombic)(001)[100]×[010]268.0186.0231.095.066.071.037.082.049.0
2 *Si3N4 (HCP)(001)[0 1 ¯ 0]×[ 1 ¯ 00]343.0600.0343.0120.0136.0120.0124.0103.5124.0
* Rotated about the z-axis by 90°, then rotated about the new x-axis by −90°.
Table 16. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 15.
Table 16. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 15.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1YBa2C3O7-δ (Orthorhombic)(001)[100]×[010]268.0186.0231.095.066.071.037.082.049.0
2MgO(FCC)(010)[001]×[100]289.3289.3289.387.787.787.7154.77154.77154.77
Table 17. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 16.
Table 17. Cleavage system: (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 16.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1YBa2C3O7−δ (Orthorhombic)(001)[100]×[010]268.0186.0231.095.066.071.037.082.049.0
2MgO (FCC)( 1 ¯ 10)[001]×[110]343.27343.27289.333.7387.787.7154.77154.77100.8
Table 18. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 17.
Table 18. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 17.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1YBa2C3O7−δ (Orthorhombic)(001)[100]×[010]268.0186.0231.095.066.071.037.082.049.0
2SrTiO3 (Simple Cubic) (010)[001]×[100]348.17348.17348.17100.64100.64100.64454.55454.55454.55
Table 19. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 18.
Table 19. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 18.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1YBa2C3O7−δ (Orthorhombic)(001)[100]×[010]268.0186.0231.095.066.071.037.082.049.0
2SrTiO3 (Simple Cubic) ( 1 ¯ 10)[001]×[110]678.96678.96348.17−230.15100.64100.64454.55454.55123.77
Table 20. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 19.
Table 20. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 19.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1YBa2C3O7−δ (Orthorhombic)(010)[001]×[100]231.0268.0186.066.071.095.049.037.082.0
2SrTiO3 (Simple Cubic) (010)[001]×[100]348.17348.17348.17100.64100.64100.64454.55454.55454.55
Table 21. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 20.
Table 21. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 20.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1YBa2C3O7−δ (Orthorhombic)(010)[001]×[100]231.0268.0186.066.071.095.049.037.082.0
2SrTiO3 (Simple Cubic) ( 1 ¯ 10)[001]×[110]678.96678.96348.17−230.15100.64100.64454.55454.55123.77
Table 22. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 21.
Table 22. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 21.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1YBa2C3O7−δ (Orthorhombic)( 1 ¯ 00)[001]×[010]268.0231.0186.066.095.071.037.049.082.0
2SrTiO3 (Simple Cubic) (010)[001]×[100]348.17348.17348.17100.64100.64100.64454.55454.55454.55
Table 23. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 22.
Table 23. Cleavage system (crack plane)[crack front]×[initial propagation direction] and elastic stiffness constants of the component phases of bicrystalline superlattice system 22.
Material (j) #Single Crystal PhaseCleavage System c 11 (GPa) c 22 (GPa) c 33 (GPa) c 12 (GPa) c 13 (GPa) c 23 (GPa) c 44 (GPa) c 55 (GPa) c 66 (GPa)
1YBa2C3O7−δ (Orthorhombic)( 1 ¯ 00)[001]×[010]268.0231.0186.066.095.071.037.049.082.0
2SrTiO3 (Simple Cubic) ( 1 ¯ 10)[001]×[110]678.96678.96348.17−230.15100.64100.64454.55454.55123.77
The next substrate, Strontium titanate (SrTiO3), which is a perovskite, is listed in Table 11, Table 12, Table 13 and Table 14 and Table 18, Table 19, Table 20, Table 21, Table 22, Table 23 as well as in Table 33, Table 34, Table 35 and Table 36 and Table 40, Table 41, Table 42, Table 43, Table 44, Table 45. SrTiO3 has at room temperature, an ideal cubic perovskite structure with TiO6 octahedra being connected by straight chains [63]. Table 34, Table 36, and Table 41 show that A = 0.2723 and κ = 0.3614 are both less than unity ( A < 1, κ < c 11 / c 11 = 1), giving rise to imaginary roots for the ( 1 ¯ 10)[001]×[110] through-crack. This is in contrast to A = 3.6727 and κ = 16.1473 being both larger than unity (A > 1, κ > c 22 / c 11 = 1), giving rise to complex roots for the (010)[001]×[100] through-crack, as shown in Table 33, Table 35, and Table 40. It can then be inferred that 1 ¯ 10 <001>×<110> would constitute an easy cleavage system, while {010}<001>×<100> would be deemed difficult.
The third and last substrate studied here is hexagonal close-packed (HCP) Si3N4, which is listed in Table 2, Table 3, Table 4 and Table 5, Table 10, and Table 15 as well as in Table 24, Table 25, Table 26 and Table 27, Table 32, and Table 37. Table 24, Table 26, and Table 32 show that A = 0.7433 < 1 and κ = 0.9206 < c 22 / c 11 = 1.3226, giving rise to imaginary roots for the (001)[0 1 ¯ 0]×[ 1 ¯ 00] through-crack. In a similar vein, A = 0.7433 < 1 and κ = 0.5263 < c 22 / c 11 = 0.5717, giving rise to imaginary roots for the ( 1 ¯ 00)[0 1 ¯ 0]×[001] through-crack, as shown in Table 25, Table 27, and Table 37. It can then be inferred that {001}<0 1 ¯ 0>×< 1 ¯ 00> and { 1 ¯ 00}<0 1 ¯ 0>×<001> would both constitute easy cleavage systems.
Fully oxidized (non-superconducting) tetragonal YBa2C3O7, the second nano-film investigated here, is listed in Table 10, Table 11, Table 12, Table 13 and Table 14 and also in Table 32, Table 33, Table 34, Table 35 and Table 36. Granozio and di Uccio [80] have also presented approximate theoretical results of fully oxidized YBCO’s (δ = 0, 1), and concluded that the three lowest surface energies follow the inequality: γ (001) < γ (100) < γ (010). Furthermore, based on the experimental results from transmission electron microscopy [81], X-ray photo-emission microscopy [82], low-energy ion scattering spectroscopy [83], and surface polarity [84] analyses performed on fully oxidized YBa2C3O7 crystals, these authors [80] have shown that the low energy cut is between the Ba=O and Cu=O planes. Table 32, Table 33 and Table 34 show that A = 1.9077 and κ = 3.1514 are both larger than unity (A > 1, κ > c 22 / c 11 = 1), giving rise to complex roots for the (010)[001]×[100] through-crack. This is in contrast to A = 0.7647 and κ = 0.7112 being both less than unity ( A < 1, κ < c 11 / c 11 = 1), giving rise to imaginary roots for the ( 1 ¯ 10)[001]×[110] through-crack, as shown in Table 35 and Table 36. It can then be inferred that 1 ¯ 10 <001>×<110> would constitute an easy cleavage system, while {010}<001>×<100> would be deemed difficult.
Finally, the third nano-film, YBa2C3O7−δ (orthorhombic), investigated here is a high TC superconductor, and is listed in Table 15, Table 16, Table 17, Table 18, Table 19, Table 20, Table 21, Table 22 and Table 23 in addition to Table 37, Table 38, Table 39, Table 40, Table 41, Table 42, Table 43, Table 44 and Table 45. As can be seen in Table 2, Table 3, Table 4, Table 5, Table 6 and Table 7 of Chaudhuri [57], all the cleavage systems are predicted to be easy, which are in agreement with the experimentally observed fracture characteristics of YBa2C3O7−δ due to Cook et al. [85], Raynes et al. [86], and Goyal et al. [87], among others; see also Granozio and di Uccio [80] for a summary of the available experimental results. Here, A = 0.764 < 1 and κ = 0.5784 < c 22 / c 11 = 0.8331, giving rise to imaginary roots for the (001)[100]×[010] through-crack, as shown in Table 37, Table 38, Table 39, Table 40 and Table 41; see also Table 6 of Chaudhuri [57]. Similarly, A = 0.8971 < 1 and κ = 0.9406 < c 22 / c 11 = 1.0771, giving rise to imaginary roots for the (010)[001]×[100] through-crack, as shown in Table 42 and Table 43; see also Table 2 of Chaudhuri [57]. Likewise, A = 0.8971 < 1 and κ = 0.817 < c 22 / c 11 = 0.9284, giving rise to imaginary roots for the ( 1 ¯ 00)[001]×[010] through-crack, as shown in Table 44 and Table 45; see also Table 4 of Chaudhuri [57]. It can then be inferred that {001}<100>×<010> would constitute an easy cleavage system.
Table 24. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 1.
Table 24. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 1.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
Au (FCC)(010)[001]×[100]2.85221.04.9777ComplexDifficult
Si3N4 (HCP)(001)[0 1 ¯ 0]×[ 1 ¯ 00]0.74331.32260.9206ImaginaryEasy
Table 25. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 2.
Table 25. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 2.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
Au (FCC)(010)[001]×[100]2.85221.04.9777ComplexDifficult
Si3N4 (HCP)( 1 ¯ 00)[0 1 ¯ 0]×[001]0.74330.57170.5263ImaginaryEasy
Table 26. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 3.
Table 26. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 3.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
Au (FCC)( 1 ¯ 10)[001]×[110]0.34941.00.2487ImaginaryEasy
Si3N4 (HCP)(001)[0 1 ¯ 0]×[ 1 ¯ 00]0.74331.32260.9206ImaginaryEasy
Table 27. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 4.
Table 27. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 4.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
Au (FCC)( 1 ¯ 10)[001]×[110]0.34941.00.2487ImaginaryEasy
Si3N4 (HCP)( 1 ¯ 00)[0 1 ¯ 0]×[001]0.74330.57170.5263ImaginaryEasy
Table 28. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 5.
Table 28. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 5.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
Au (FCC)(010)[001]×[100]2.85221.04.9777ComplexDifficult
MgO (FCC)(010)[001]×[100]1.53541.01.8329ComplexDifficult
Table 29. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 6.
Table 29. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 6.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
Au (FCC)(010)[001]×[100]2.85221.04.9777ComplexDifficult
MgO (FCC)( 1 ¯ 10)[001]×[110]0.65131.00.6297ImaginaryEasy
Table 30. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 7.
Table 30. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 7.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
Au (FCC)( 1 ¯ 10)[001]×[110]0.34941.00.2487ImaginaryEasy
MgO (FCC)(010)[001]×[100]1.53541.01.8329ComplexDifficult
Table 31. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 8.
Table 31. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 8.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
Au (FCC)( 1 ¯ 10)[001]×[110]0.34941.00.2487ImaginaryEasy
MgO (FCC)( 1 ¯ 10)[001]×[110]0.65131.00.6297ImaginaryEasy
Table 32. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 9.
Table 32. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 9.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
YBa2C3O7 (Tetragonal)(010)[001]×[100]1.90771.03.1514ComplexDifficult
Si3N4 (HCP)( 1 ¯ 00)[0 1 ¯ 0]×[001]0.74330.57170.9206ImaginaryEasy
Table 33. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 10.
Table 33. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 10.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
YBa2C3O7 (Tetragonal)(010)[001]×[100]1.90771.03.1514ComplexDifficult
SrTiO3 (Simple Cubic)(010)[001]×[100]3.67271.016.1473ComplexDifficult
Table 34. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 11.
Table 34. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 11.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
YBa2C3O7 (Tetragonal)(010)[001]×[100]1.90771.03.1514ComplexDifficult
SrTiO3 (Simple Cubic) ( 1 ¯ 10)[001]×[110]0.27231.00.3614ImaginaryEasy
Table 35. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 12.
Table 35. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 12.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
YBa2C3O7 (Tetragonal)( 1 ¯ 10)[001]×[110]0.76471.00.7112ImaginaryEasy
SrTiO3 (Simple Cubic)(010)[001]×[100]3.67271.016.1473ComplexDifficult
Table 36. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 13.
Table 36. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 13.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
YBa2C3O7 (Tetragonal)( 1 ¯ 10)[001]×[110]0.76471.00.7112ImaginaryEasy
SrTiO3 (Simple Cubic) ( 1 ¯ 10)[001]×[110]0.27231.00.3614ImaginaryEasy
Table 37. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 14.
Table 37. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 14.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
YBa2C3O7−δ (Orthorhombic)(001)[100]×[010]0.7640.83310.5784ImaginaryEasy
Si3N4 (HCP)(001)[0 1 ¯ 0]×[ 1 ¯ 00]0.74331.32260.9206ImaginaryEasy
Table 38. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 15.
Table 38. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 15.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
YBa2C3O7−δ (Orthorhombic)(001)[100]×[010]0.7640.83310.5784ImaginaryEasy
MgO (FCC)(010)[001]×[100]1.53541.01.8329ComplexDifficult
Table 39. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 16.
Table 39. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 16.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
YBa2C3O7−δ (Orthorhombic)(001)[100]×[010]0.7640.83310.5784ImaginaryEasy
MgO (FCC)( 1 ¯ 10)[001]×[110]0.65131.00.6297ImaginaryEasy
Table 40. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 17.
Table 40. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 17.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
YBa2C3O7−δ (Orthorhombic)(001)[100]×[010]0.7640.83310.5784ImaginaryEasy
SrTiO3 (Simple Cubic) (010)[001]×[100]3.67271.016.1473ComplexDifficult
Table 41. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 18.
Table 41. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 18.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
YBa2C3O7−δ (Orthorhombic)(001)[100]×[010]0.7640.83310.5784ImaginaryEasy
SrTiO3 (Simple Cubic) ( 1 ¯ 10)[001]×[110]0.27231.00.3614ImaginaryEasy
Table 42. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 19.
Table 42. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 19.
Material (j) #Cleavage SystemA c 22 c 11 κ RootsCleavage System: Easy or Difficult
YBa2C3O7−δ (Orthorhombic)(010)[001]×[100]0.89711.07710.9406ImaginaryEasy
SrTiO3 (Simple Cubic) (010)[001]×[100]3.67271.016.1473ComplexDifficult
Table 43. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 20.
Table 43. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 20.
Material (j) #Cleavage SystemA c 22 c 11 κ RootsCleavage System: Easy or Difficult
YBa2C3O7−δ (Orthorhombic)(010)[001]×[100]0.89711.07710.9406ImaginaryEasy
SrTiO3 (Simple Cubic) ( 1 ¯ 10)[001]×[110]0.27231.00.3614ImaginaryEasy
Table 44. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 21.
Table 44. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 21.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
YBa2C3O7−δ (Orthorhombic)( 1 ¯ 00)[001]×[010]0.89710.92840.817ImaginaryEasy
SrTiO3 (Simple Cubic) (010)[001]×[100]3.67271.016.1473ComplexDifficult
Table 45. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 22.
Table 45. Normalized elastic parameter, roots of characteristic equation, and the nature (easy or difficult) of the through-thickness cleavage system of bicrystalline superlattice system 22.
Material (j) #Cleavage System A c 22 c 11 κ RootsCleavage System: Easy or Difficult
YBa2C3O7−δ (Orthorhombic)( 1 ¯ 00)[001]×[010]0.89710.92840.817ImaginaryEasy
SrTiO3 (Simple Cubic) ( 1 ¯ 10)[0 1 ¯ 01]×[110]0.27231.00.3614ImaginaryEasy
Table 46 summarizes the nature (easy or difficult) of the cleavage system in component phases of the afore-mentioned bicrystalline superlattice systems. It also lists the real or complex eigenvalues, s = 0.5 ± iε, of these bicrystalline superlattices. These results suggest that the interfacial cracks would propagate in the mixed (I/II) mode, primarily when both the component phases are characterized by difficult cleavage systems (complex roots), the exception being perovskite SrTiO3 serving as the substrate (for Yba2C3O7−δ or Yba2C3O7 nano-films). A plausible reason for this exceptional behavior of SrTiO3 may lie in its unusually high shear stiffness, c66, which is substantially greater than its longitudinal stiffness, c11, in combination with easiest cleavage system of Yba2C3O7−δ (or Yba2C3O7). This is in contrast with other cubic mono-crystals, such as FCC rock salt MgO. This fact results in a negative Poisson’s ratio effect, when rotated about the [001] axis by 45°, which is not generally encountered in cubic crystal elasticity. Additionally, this behavior is also in contrast with other easy cleavage systems of YBa2C3O7−δ deposited on the same 45° rotated SrTiO3, as shown in Table 43, Table 45, and Table 46.
Table 46. Real or complex eigenvalues of the bicrystalline superlattice systems with through interfacial cracks.
Table 46. Real or complex eigenvalues of the bicrystalline superlattice systems with through interfacial cracks.
Bicrystal System #Nano-Film/SubstrateCleavage SystemsRootsCleavage System: Easy or Difficults = 0.5 or
s = 0.5 ± iε
1Au/Si3N4(010)[001]×[100]/(001)[0 1 ¯ 0]×[ 1 ¯ 00]Complex/ImaginaryDifficult/Easy0.5
2Au/Si3N4(010)[001]×[100]/( 1 ¯ 00)[0 1 ¯ 0]×[001]Complex/ImaginaryDifficult/Easy0.5
3Au/Si3N4( 1 ¯ 10)[001]×[110]/(001)[00]×[ 1 ¯ 00]Imaginary/ImaginaryEasy/Easy0.5
4Au/Si3N4( 1 ¯ 10)[001]×[110]/( 1 ¯ 00)[0 1 ¯ 0]×[001]Imaginary/ImaginaryEasy/Easy0.5
5Au/MgO(010)[001]×[100]/(010)[001]×[100]Complex/ComplexDifficult/Difficult0.5 ± 0.3814i
0.5 ± 0.2108i
6Au/MgO(010)[001]×[100]/(10)[001]×[110]Complex/ImaginaryDifficult/Easy0.5
7Au/MgO( 1 ¯ 10)[001]×[110]/(010)[001]×[100]Imaginary/ComplexEasy/Easy0.5
8Au/MgO( 1 ¯ 10)[001]×[110]/( 1 ¯ 10)[001]×[110]Imaginary/ImaginaryEasy/Easy0.5
9YBa2C3O7/Si3N4(010)[001]×[100]/( 1 ¯ 00)[0 1 ¯ 0]×[001]Complex/ImaginaryDifficult/Easy0.5
10YBa2C3O7/SrTiO3(010)[001]×[100]/(010)[001]×[100]Complex/ComplexDifficult/Difficult0.5 ± 0.7636i
11YBa2C3O7/SrTiO3(010)[001]×[100]/( 1 ¯ 10)[001]×[110]Complex/ImaginaryDifficult/Easy0.5
12YBa2C3O7/SrTiO3( 1 ¯ 10)[001]×[110]/(010)[001]×[100]Imaginary/ComplexEasy/Difficult0.5
13YBa2C3O7/SrTiO3( 1 ¯ 10)[001]×[110]/( 1 ¯ 10)[001]×[110]Imaginary/ImaginaryEasy/Easy0.5 ± 0.0580i
14YBa2C3O7−δ/Si3N4(001)[100]×[010]/(001)[0 1 ¯ 0]×[ 1 ¯ 00]Imaginary/ImaginaryEasy/Easy0.5
15YBa2C3O7−δ/MgO(001)[100]×[010]/(010)[001]×[100]Imaginary/ComplexEasy/Difficult0.5
16YBa2C3O7−δ/MgO(001)[100]×[010]/( 1 ¯ 10)[001]×[110]Imaginary/ImaginaryEasy/Easy0.5
17YBa2C3O7−δ/SrTiO3(001)[100]×[010]/(010)[001]×[100]Imaginary/ComplexEasy/Difficult0.5
18YBa2C3O7−δ/SrTiO3(001)[100]×[010]/(10)[001]×[110]Imaginary/ImaginaryEasy/Easy0.5 ± 0.1757i
19YBa2C3O7−δ/SrTiO3(010)[001]×[100]/(010)[001]×[100]Imaginary/ImaginaryEasy/Easy0
20YBa2C3O7−δ/SrTiO3(010)[001]×[100]/( 1 ¯ 10)[001]×[110]Imaginary/ImaginaryEasy/Easy0
21YBa2C3O7−δ/SrTiO3( 1 ¯ 00)[001]×[010]/(010)[001]×[100]Imaginary/ImaginaryEasy/Easy0
22YBa2C3O7−δ/SrTiO3( 1 ¯ 00)[001]×[010]/( 1 ¯ 10)[001]×[110]Imaginary/ImaginaryEasy/Easy0

9.2. Superlattice Trapping and Superlattice Crack Deflection (SCD)

The theory of lattice crack deflection (LCD) is discussed in Chaudhuri [49]. Table 47 displays the structures and elastic compliance constants of mono-crystalline FCC transition metal Au, FCC rock salt MgO, cubic perovskite SrTiO3, HCP ceramic Si3N4, fully oxidized tetragonal YBa2C3O7, and orthorhombic (superconducting) YBa2C3O7−δ [49,57,78]. Table 48 shows the results for computed lattice crack deflection (LCD) parameters (energy barrier) and associated bond shear strains at crack deflection from a difficult cleavage system to an easy one, and their correlations with the anisotropic ratios relating to the difficult cleavage system along with Bravais lattice and structure. Only two crack systems are considered: {010}〈001〉×<100]> and 1 ¯ 10 <001>×<110]>.
Table 47. Structures and elastic compliance constants of selected single crystals [49,78].
Table 47. Structures and elastic compliance constants of selected single crystals [49,78].
Single CrystalBravais LatticeStructure s 11 (10−2 GPa−1) s 22 (10−2 GPa−1) s 33 (10−2 GPa−1) s 12 (10−2 GPa−1) s 13 (10−2 GPa−1) s 23 (10−2 GPa−1) s 44 (10−2 GPa−1) s 55 (10−2 GPa−1) s 66 (10−2 GPa−1)
Au [79,88]FCCFCC2.3552.3552.355−1.081−1.081−1.0812.40962.40962.4096
MgO [79]FCCRock Salt0.40240.40240.4024−0.0936−0.0936−0.09360.64610.64610.6461
SrTiO3 [79]Simple CubicPerovskite0.330.330.33−0.074−0.074−0.0740.220.220.22
Si3N4HCPHCP0.35760.35760.1852−0.1255−0.0464−0.04640.80650.80650.9662
YBa2C3O7TetragonalPerovskite0.63890.63891.1186−0.1304−0.339−0.3392.02.01.1765
Y Ba2C3O7−δOrthorhombicPerovskite0.50030.46480.7054−0.0678−0.1564−0.21152.04082.70271.2195
For mono-crystalline FCC transition metals, {010}〈001〉×<100> is deemed to be a difficult cleavage system for the reasons explained above, while 1 ¯ 10 <001>×<110]> is considered to be the preferred one for crack propagation. This is illustrated in Figure 26a in Chaudhuri [49]. A nonvanishing lattice crack deflection (LCD) energy barrier implies that a {010}〈001〉×<100]> through-crack in such single crystals would not deflect right at the appropriate Griffith/Irwin critical stress intensity factor (Kc) for mixed-mode propagation because of the lattice effect, but would require additional bond shear strains for Au (Table 48). In the case of the nonvanishing lattice crack deflection (LCD) barrier, e.g., in Au with a moderately high anisotropic ratio, A = 2.8481 > 1, the difficult {010}〈001〉×<100]> crack may initially become lattice-trapped and/or propagate in a “difficult” manner till an applied load somewhat higher than its Griffith mixed-mode counterpart is reached, and then only deflect into the easy cleavage system, 1 ¯ 10 <001>×<110]>. In addition, bond breaking would not be continuous, but abrupt. In contrast, for the same crystal with a very low modified anisotropic ratio, A = 0.3494 < 1, the lattice crack deflection (LCD) barrier vanishes and the easy 1 ¯ 10 <001>×<110]> crack would begin to propagate right at the Griffith/Irwin critical stress intensity factor. There would be no crack turning.
Table 48. Easy/difficult cleavage system, lattice crack deflection (LCD) barrier, and associated bond shear strains in selected single crystals.
Table 48. Easy/difficult cleavage system, lattice crack deflection (LCD) barrier, and associated bond shear strains in selected single crystals.
CrystalEasy Cleavage System (ECS)Difficult Cleavage System (DCS)Aniso. Ratio, A, at DCSBond Shear Strain at Lattice Crack Deflection γ b d Lattice Crack Deviation LCD Parameter in DCS Δ K *
Au [88] 1 ¯ 10 <001>×<110]>{010}〈001〉×<110]>2.8522 0.64380.8364
MgO 1 ¯ 10 <001>×<110]>{010}〈001〉×<100]>1.53540.53530.6414
SrTiO3 1 ¯ 10 <001>×<110]>{010}〈001〉×<100]>3.67270.60750.5114
Si3N4( 1 ¯ 00)[0 1 ¯ 0]×[001] *-------
YBa2C3O7 1 ¯ 10 <001>×<110]>{010}〈001〉×<100]>1.90770.50550.6071
YBa2C3O7−δ{001}<100>×<010>-------
* Rotated about the z-axis by 90°, then rotated about the new x-axis by −90°.
Bicrystals form superlattices have not been discussed, to the author’s knowledge, in anisotropic fracture mechanics literature [88]. However, the geometric mean of the two constituent phases would serve as a reasonably accurate procedure for the computation of the bond shear strain at the superlattice crack deviation, γ b d s , and superlattice crack deflection (SCD) barrier, Δ K s * ; see also Chaudhuri [88]. Rigorous proof of this is currently being worked out and will be presented in the near future. The numerical results are shown in Table 49. For a bicrystalline superlattice, e.g., Au/MgO (respectively, YBa2C3O7/SrTiO3), with both difficult cleavage systems, (010)[001]×[100]/(010)[001]×[100], serving as the interface, with the SCD barrier, Δ K s * , value of 0.7324 (resp., 0.6240), the interfacial crack would encounter a tough interface and would initially be superlattice-trapped and/or experience a mixed-mode propagation in a “difficult” manner till an applied load somewhat higher than its Griffith/Irwin mixed-mode interfacial fracture toughness counterpart — quantified by Δ K s * — is reached, and thence deflect into the available easier cleavage system, 1 ¯ 10 <001>×<110]>, of the component phase with the lower LCD barrier, Δ K * = 0.6414 for MgO (resp. 0.5114 for SrTiO3). In addition, bond breaking would not be continuous, but abrupt. In contrast, for the same bicrystalline superlattice, Au/MgO, with both easy ( 1 ¯ 10)[001]×[110]/( 1 ¯ 10)[001]×[110] cleavage systems serving as the interface, the SCD barrier, Δ K s * , vanishes, and the easy interfacial crack would begin to propagate (in the absence of mode mixity) in a self-similar manner right at the Griffith/Irwin critical stress intensity factor. The bond breaking would be smooth and continuous. Interestingly, for the Au/MgO or YBa2C3O7/SrTiO3 superlattice, with one easy and the second one difficult, either (010)[001]×[100]/( 1 ¯ 10)[001]×[110] or ( 1 ¯ 10)[001]×[110]/(010)[001]×[100] cleavage systems serving as the interface, the SCD barrier, Δ K s * , also vanishes, and the interfacial crack would begin to propagate (in the absence of mode mixity) on the easier side of and parallel to the interface at the Griffith/Irwin critical stress intensity factor. Bond breaking would be smooth and continuous on the easier side, but discontinuous and abrupt on the tougher side of the interface.
Table 49. Cleavage system, lattice crack deflection barrier, and associated bond shear strains in selected bicrystalline superlattice systems.
Table 49. Cleavage system, lattice crack deflection barrier, and associated bond shear strains in selected bicrystalline superlattice systems.
Bicrystal System #Nano-film/SubstrateCleavage SystemsCleavage System: Easy or DifficultBond Shear Strain at Superlattice Crack Deflection γ b d s Superlattice Crack Deviation (SCD) Parameter in DCS Δ K s *
1Au/Si3N4(010)[001]×[100]/(001)[0 1 ¯ 0]×[ 1 ¯ 00]Difficult/Easy00
2Au/Si3N4(010)[001]×[100]/( 1 ¯ 00)[0 1 ¯ 0]×[001]Difficult/Easy00
3Au/Si3N4( 1 ¯ 10)[001]×[110]/(001)[00]×[ 1 ¯ 00]Easy/Easy00
4Au/Si3N4( 1 ¯ 10)[001]×[110]/( 1 ¯ 00)[0 1 ¯ 0]×[001]Easy/Easy00
5Au/MgO(010)[001]×[100]/(010)[001]×[100]Difficult/Difficult0.47100.7324
6Au/MgO(010)[001]×[100]/(10)[001]×[110]Difficult/Easy00
7Au/MgO( 1 ¯ 10)[001]×[110]/(010)[001]×[100]Easy/Difficult00
8Au/MgO( 1 ¯ 10)[001]×[110]/( 1 ¯ 10)[001]×[110]Easy/Easy00
9YBa2C3O7/Si3N4(010)[001]×[100]/( 1 ¯ 00)[0 1 ¯ 0]×[001]Difficult/Easy00
10YBa2C3O7/SrTiO3(010)[001]×[100]/(010)[001]×[100]Difficult/Difficult0.52020.6240
11YBa2C3O7/SrTiO3(010)[001]×[100]/( 1 ¯ 10)[001]×[110]Difficult/Easy00
12YBa2C3O7/SrTiO3( 1 ¯ 10)[001]×[110]/(010)[001]×[100]Easy/Difficult00
13YBa2C3O7/SrTiO3( 1 ¯ 10)[001]×[110]/( 1 ¯ 10)[001]×[110]Easy/Easy00
14YBa2C3O7−δ/Si3N4(001)[100]×[010]/(001)[0 1 ¯ 0]×[ 1 ¯ 00]Easy/Easy00
15YBa2C3O7−δ/MgO(001)[100]×[010]/(010)[001]×[100]Easy/Difficult00
16YBa2C3O7−δ/MgO(001)[100]×[010]/( 1 ¯ 10)[001]×[110]Easy/Easy00
17YBa2C3O7−δ/SrTiO3(001)[100]×[010]/(010)[001]×[100]Easy/Difficult00
18YBa2C3O7−δ/SrTiO3(001)[100]×[010]/(10)[001]×[110]Easy/Easy00
19YBa2C3O7−δ/SrTiO3(010)[001]×[100]/(010)[001]×[100]Easy/Difficult00
20YBa2C3O7−δ/SrTiO3(010)[001]×[100]/( 1 ¯ 10)[001]×[110]Easy/Easy00
21YBa2C3O7−δ/SrTiO3( 1 ¯ 00)[001]×[010]/(010)[001]×[100]Easy/Difficult00
22YBa2C3O7−δ/SrTiO3( 1 ¯ 00)[001]×[010]/( 1 ¯ 10)[001]×[110]Easy/Easy00
For the Au/Si3N4 superlattice, with both easy cleavage systems, either ( 1 ¯ 10)[001]×[110]/(001)[0 1 ¯ 0]×[ 1 ¯ 00] or ( 1 ¯ 10)[001]×[110]/( 1 ¯ 00)[0 1 ¯ 0]×[001], serving as the interface, the SCD barrier, Δ K s * , again vanishes, and the easy interfacial crack would begin to propagate (in the absence of mode mixity) in a self-similar manner right at the Griffith/Irwin critical stress intensity factor. Bond breaking would be smooth and continuous. As before, the same superlattice, with one easy and the second one difficult, either (010)[001]×[100]/(001)[0 1 ¯ 0]×[ 1 ¯ 00] or (010)[001]×[100]/( 1 ¯ 00)[0 1 ¯ 0]×[001] cleavage systems serving as the interface, Δ K s * , also vanishes, and the interfacial crack would begin to propagate (in the absence of mode mixity) on the easier side of and parallel to the interface at the Griffith/Irwin critical stress intensity factor. Bond breaking would be smooth and continuous on the easier side, but discontinuous and abrupt on the tougher side of the interface. A similar situation prevails for the YBa2C3O7/Si3N4 superlattice, with one difficult and the second one easy, (010)[001]×[100]/( 1 ¯ 00)[0 1 ¯ 0]×[001], serving as the interface.
For the orthorhombic perovskite/HCP YBa2C3O7−δ/Si3N4 superlattice, with both easy (001)[100]×[010]/(001)[0 1 ¯ 0]×[ 1 ¯ 00] cleavage systems serving as the interface, the SCD barrier, Δ K s * , again vanishes, and the easy interfacial crack would begin to propagate (in the absence of mode mixity) in a self-similar manner right at the Griffith/Irwin critical stress intensity factor. Similar results follow for the orthorhombic perovskite/FCC rock salt YBa2C3O7−δ/MgO superlattice, with both easy (001)[100]×[010]/( 1 ¯ 10)[001]×[110] cleavage systems serving as the interface.
For the perovskite orthorhombic/FCC rock salt MgO bicrystalline superlattice, with both easy (001)[100]×[010]/( 1 ¯ 10)[001]×[110] cleavage systems serving as the interface, Δ K s * would again vanish, and the resulting easy interfacial crack would begin to propagate (in the absence of mode mixity) in a self-similar manner right at the Griffith/Irwin critical stress intensity factor. Bond breaking would be smooth and continuous. As before for Δ K s * for the same superlattice, with one easy and the second one difficult, (001)[100]×[010]/(010)[001]×[100] would also vanish, and the interfacial crack would begin to propagate (in the absence of mode mixity) on the easier side of and parallel to the interface at the Griffith/Irwin critical stress intensity factor. Bond breaking would be smooth and continuous on the easier side, but discontinuous and abrupt on the tougher side of the interface.
Finally, for the perovskite orthorhombic/cubic perovskite YBa2C3O7−δ/SrTiO3 (respectively, YBa2C3O7/SrTiO3) bicrystalline superlattice, with both easy (001)[100]×[010]/( 1 ¯ 10)[001]×[110] (respectively, ( 1 ¯ 10)[001]×[110]/( 1 ¯ 10)[001]×[110]) cleavage systems serving as the interface, Δ K s * would again vanish, and the resulting easy interfacial crack would experience a mixed-mode propagation/growth right at the Griffith/Irwin critical complex stress intensity factor. Rice [9] has discussed the computation and interpretation of the resulting complex stress intensity factor (S.I.F.) for an isotropic bimaterial interface crack; see Section 9.3 below. Bond breaking is expected to be smooth and continuous. The reason for this exceptional behavior of SrTiO3 lies, as has been explained above, in its unusually high shear stiffness, c66, which is substantially greater than its longitudinal stiffness, c11. However, Δ K s * for the same superlattices, with one easy and the second one difficult, (001)[100]×[010]/(010)[001]×[100] for YBa2C3O7−δ/SrTiO3, or ( 1 ¯ 10)[001]×[110]/(010)[001]×[100] and (010)[001]×[100]/( 1 ¯ 10)[001]×[110] for YBa2C3O7/SrTiO3 also vanish, and the interfacial crack would begin to propagate (in the absence of mode mixity) on the easier side of and parallel to the interface at the Griffith/Irwin critical stress intensity factor. Bond breaking would be smooth and continuous on the easier side, but discontinuous and abrupt on the tougher side of the interface.

9.3. Complex Stress Intensity Factor (S.I.F.) and Raman Spectroscopic Surface Measurement

9.3.1. Complex Stress Intensity Factor (S.I.F.)

Rice [9] has discussed the validity of the complex stress intensity factor, K = KI + iKII, as a characterizing parameter in the context of a two-dimensional bimaterial interface crack for material systems (of isotropic phases), for small-scale yielding and small-scale contact zones at the crack tip. The proportion of opening and shearing modes at the interfacial crack tip is characterized by means of the phase angle, given by ψ = Arctan K I I / K I . Wang et al. [12] and Wu [10] have extended Rice’s [9] analysis, and presented corresponding relations for interface cracks in monoclinic composites.

9.3.2. Raman Spectroscopic Surface Measurement of Carbon/Graphite Fiber–Epoxy Interfacial Bond

It has been argued by Chaudhuri et al. [36] that, in connection with a Raman spectroscopic study of a carbon fiber–epoxy interfacial debond, in real-life composite materials, the fiber–matrix interface (a sharp material discontinuity) is replaced by an interphase region (of the order of 0.5 μm thickness), which permits us to ignore the contact of the Comninou [89] type (ε = 0). Various factors that contribute to the formation of this interphase region are molecular entanglement following interdiffusion, electrostatic attraction, cationic groups at the ends of molecules attracted to an anionic surface that result in polymer matrix orientation at the fiber surface, chemical reaction, mechanical keying, etc. [90]. Additionally, the residual stress and coupling agents (e.g., silanes on glass fibers) applied to the fiber’s surface have pronounced effects on the interfacial bond strength. Under these circumstances, the macroscopic (mean field, e.g., linear elasticity) analysis presented above may be considered to be only a first approximation to the detailed microscopic (fluctuating) state of stress at the carbon fiber–epoxy debond tip. It may further be noted that the oscillatory behavior characterizing the interpenetration of the component phases is not that unrealistic after all, but is rather a first-order approximation of the interdiffusion followed by molecular entanglement, chemical bonding, and other similar microscopic (kinetic) phenomena studied by materials scientists.
A close scrutiny of two asymptotic solutions for bimaterials involving (i) the complex eigenvalues (ε > 0) and (ii) their real counterpart (ε = 0), in light of the Raman spectroscopic measurements, reveals that they represent two extreme cases of the interfacial (interphase) region of the carbon/graphite fibers with a polymeric matrix, such as epoxy [36]. While case (i) corresponds to an idealized carbon fiber with the surface layer(s) being comprised of completely disordered (amorphous) carbon resembling activated charcoal (R = I A 1 g / I E 2 g = 1), case (ii) represents an idealized version of graphite fiber with the surface layer(s) consisting of a completely ordered (crystalline) carbon lattice resembling stress-annealed pyrolite graphite (R = 0). Actual carbon/graphite fibers, both commercial and experimental (experimentally investigated by Chaudhuri et al. [36]), fall somewhere in between the two idealized extremes with 0 < R < 1. Therefore, the interfacial debond nucleation and propagation in an actual carbon/graphite fiber/epoxy matrix composite will be governed by the solution with ε replaced by [36]:
ε ¯ = ε R = ε I A 1 g / I E 2 g ,
where I A 1 g and I E 2 g represent the relative (Raman) intensities of the peaks corresponding to the A1g (weak) and E2g (strong) modes, respectively. R represents the degree of chemical functionality (potential for strong or covalent bond formation). R ranges from 0.22 (Morganite I) to 0.85 (Thornel 10), with the industry standard AS4 being 0.811, as reported by Tuinstra and Koenig [91]. For the experimental fibers (grown in a carefully controlled environment to mimic their commercial counterparts) investigated by Chaudhuri et al. [36], R ranges from 0.279 to 0.785.

9.4. Through-Thickness Distribution of Stress Intensity Factors (Fracture Toughness) and Energy Release Rates (Fracture Energy)

Figure 2a,b show variations in the normalized stress intensity factor, K i j * z = K i j z / K i j , max ) , i = I, II, III, j = 1, 2, through the thickness of a bicrystalline superlattice plate, weakened by the through-crack investigated here. Figure 2a shows the through-thickness variation in the stress intensity factor for a far-field symmetrically distributed hyperbolic sine (modes I, II, and III) load, while its antisymmetric counterpart, not encountered in two-dimensional analyses can be deemed associated with the singular residual stress field [92] and is displayed in Figure 2b. Figure 3 shows the corresponding variation in the energy release rate, G * , through the top half of the plate thickness. For through-thickness symmetric far-field loading, the crack is expected to grow through thickness in a stable manner till the stress intensity factor or the energy release rate reaches its critical value at the maximum location (just below the top and bottom surfaces). With a further increase in the magnitude of the far-field loading, unstable crack growth is expected to progressively spread throughout the plate thickness. For skew-symmetric loading, as reported on earlier occasions [16], the bottom half will experience crack closure. Such types of results describing the three-dimensional distribution of stress intensity factors and energy release rates have only recently become available in fracture mechanics literature.

10. Summary and Conclusions

A recently developed eigenfunction expansion method, based in part on the separation of the thickness variable and partly utilizing a modified Frobenius-type series expansion in terms of affine-transformed x-y coordinate variables of the Eshelby–Stroh type, is employed for obtaining three-dimensional asymptotic displacement and stress fields in the vicinity of the front of an interfacial crack weakening an infinite pie-shaped bicrystalline superlattice plate, of finite thickness, formed as a result of a mono-crystalline metal or superconductor film deposited over a substrate. The bicrystalline superlattice is made of orthorhombic (tetragonal, hexagonal, and cubic as special cases) phases, and is subjected to the far-field extension/bending, in-plane shear/twisting, and anti-plane shear loadings distributed through the thickness. Crack-face boundary and interface contact conditions as well as those that are prescribed on the top and bottom surfaces of the bicrystalline superlattice plate are exactly satisfied.
It also extends a recently developed concept of the lattice crack deflection (LCD) barrier to its superlattice counterpart, christened the superlattice crack deflection (SCD) energy barrier, for studying interfacial crack path instability, which can explain crack deflection from a difficult interface to an easier neighboring cleavage system. Additionally, the relationships of the nature (easy/easy, easy/difficult, or difficult/difficult) of interfacial cleavage systems based on the present solutions with the structural chemistry aspects of the component phases (such as orthorhombic, tetragonal, hexagonal, as well as FCC (face-centered cubic) transition metals and perovskites) of the superlattice are also investigated.
Important conclusions drawn from this study can be listed as follows:
(i) Atomistic-scale modeling of interfacial cracks requires a consideration of both long-range elastic interactions and short-range chemical reactions. The Griffith thermodynamic-based theory does not take the latter into account, and hence must be regarded as only a necessary condition (albeit being still very useful and widely employed), but not as sufficient.
(ii) The effect of short-range chemical reactions can be adequately captured by elastic properties-based parameters, such as the planar anisotropic ratio, Aj, or equivalently, the normalized elastic parameter, κj, j = 1, 2. This is because the elastic properties are controlled by various aspects of the underlying structural chemistry of single crystals, such as the Bravais lattice type, bonding (covalent, ionic, and metallic), bonding (including hybridized) orbitals, electro-negativity of constituent atoms in a compound, polarity, etc. More specifically, the elastic properties of superconducting YBa2Cu3O7−δ are strongly influenced by oxygen non-stoichiometry (as well as various structural defects).
(iii) A single dimensionless parameter, such as the planar anisotropic ratio, Aj, or equivalently, the normalized elastic parameter, κj, j = 1, 2, can serve as the Holy Grail quantity for an a priori determination of the status of a cleavage system to be easy or difficult, very much akin to the Reynold’s number for fluid flow problems, crossing a critical value that signifies transition from one regime to another. Here, the planar anisotropic ratio, Aj, or equivalently, the normalized elastic parameter, κj, j = 1, 2, for a (010)[001]×[100] cleavage system, crossing the critical value of 1 or c 22 ( j ) / c 11 ( j ) , j = 1, 2, respectively, signifies transitioning from self-similar crack growth or propagation to crack deflection or turning from a difficult cleavage system onto a nearby easy one. This is a significant qualitative as well as quantitative improvement over two-parameter-based models, suggested by previous researchers, e.g., [13], in the context of two-dimensional anisotropic fracture mechanics.
(iv) Just as the introduction of the Reynold’s number facilitated the design and setting up of experiments in addition to experimental verification of analytical and computational solutions in fluid dynamics, the accuracy and efficacy of the available test results on elastic constants of YBa2Cu3O7−δ single crystals, measured by modern experimental techniques with resolutions at the atomic scale or nearly so, such as X-ray diffraction, ultrasound technique, neutron diffraction/scattering, Brillouin spectroscopy/scattering, resonant ultrasound spectroscopy, and the like, are assessed with a powerful theoretical analysis of crack path stability/instability, in part based on a single dimensionless parameter, such as the planar anisotropic ratio, Aj, j = 1, 2.
(v) Experimental determination of the surface energy, Γj, j = 1, 2, of the component phases or the corresponding interfacial energy, Γint, of a bicrystalline superlattice can sometimes be notoriously challenging, due to the presence of micro-to-nanoscale defects, such as porosity, dislocation, twin boundaries, misalignment of bonds with respect to the loading axis, and the like. In contrast, the above-derived bond shear strain at superlattice crack deflection, γ b d s , and superlattice crack deflection (SCD) barrier, Δ K s * , are, relatively speaking, much easier in comparison to the determination of surface or interfacial energy.
(vi) Computed complex eigenvalues, s = 0.5 ± iε, for Au/MgO, and YBa2C3O7/SrTiO3, bicrystalline superlattices, with (010)[001]×[100]/(010)[001]×[100] serving as the interface suggest that the corresponding interfacial cracks would propagate in a mixed (I/II) mode. Likewise, for the bicrystalline superlattice, YBa2C3O7−δ/SrTiO3 (respectively, YBa2C3O7/SrTiO3) with (001)[100]×[010]/( 1 ¯ 10)[001]×[110] (respectively, ( 1 ¯ 10)[001]×[110]/( 1 ¯ 10)[001]×[110]) cleavage systems serving as the interface, the computed eigenvalues are also complex, resulting in a mixed (I/II)-mode interfacial crack growth.
(vii) For the bicrystalline superlattice, Au/MgO (respectively, YBa2C3O7/SrTiO3), with both difficult cleavage systems, (010)[001]×[100]/(010)[001]×[100], serving as the interface, with the bond shear strain at superlattice crack deflection, γ b d s , value of 0.4710 (respectively, 0.5202) and superlattice crack deflection (SCD) barrier, Δ K s * , value of 0.7324 (resp., 0.6240), the interfacial crack would encounter a tough interface, and would initially be superlattice-trapped and/or experience a mixed-mode propagation in a “difficult” manner till an applied load somewhat higher than its Griffith/Irwin mixed-mode interfacial fracture toughness counterpart—quantified by Δ K s * —is reached, and thence deflect onto the available easier cleavage system, 1 ¯ 10 <001>×<110]>, of the component phase with the lower LCD barrier, Δ K * = 0.6414 for MgO (0.5114 for SrTiO3). In addition, bond breaking would not be continuous, but abrupt.
(viii) In contrast, for the perovskite orthorhombic/cubic perovskite YBa2C3O7−δ/SrTiO3 (respectively, YBa2C3O7/SrTiO3) bicrystalline superlattice, with both easy (001)[100]×[010]/( 1 ¯ 10)[001]×[110] (respectively, ( 1 ¯ 10)[001]×[110]/( 1 ¯ 10)[001]×[110]) cleavage systems serving as the interface, both bond shear strain at superlattice crack deflection, γ b d s , and superlattice crack deflection (SCD) barrier, Δ K s * , vanish, and the resulting easy interfacial crack would experience a mixed-mode propagation/growth right at the Griffith/Irwin-based critical complex stress intensity factor (S.I.F.), the computation and interpretation of which is expounded by Rice’s [9] extension from a two-dimensional isotropic bimaterial interface crack to its to anisotropic counterpart. Bond breaking is expected to be smooth and continuous.
(ix) For the Au/Si3N4 or YBa2C3O7/Si3N4 superlattice, with both easy cleavage systems, either ( 1 ¯ 10)[001]×[110]/(001)[0 1 ¯ 0]×[ 1 ¯ 00] or ( 1 ¯ 10)[001]×[110]/( 1 ¯ 00)[0 1 ¯ 0]×[001], serving as the interface, both bond shear strain at superlattice crack deflection, γ b d s , and superlattice crack deflection (SCD) barrier, Δ K s * , vanish, and the easy interfacial crack would begin to propagate (in the absence of mode mixity) in a self-similar manner right at the Griffith/Irwin critical stress intensity factor. Likewise, a YBa2C3O7−δ/Si3N4 superlattice, with both easy cleavage systems, either {001}<100>×<010>/(001)[0 1 ¯ 0]×[ 1 ¯ 00] or {001}<100>×<010>/( 1 ¯ 00)[0 1 ¯ 0]×[001], serving as the interface, would elicit a similar behavior. Bond breaking would be smooth and continuous.
(x) Other examples of superlattices, with both easy cleavage systems, include the Au/MgO (respectively, YBa2C3O7−δ/MgO) bicrystalline superlattice with ( 1 ¯ 10)[001]×[110/( 1 ¯ 10)[001]×[110] (respectively, {001}<100>×<010>/( 1 ¯ 10)[001]×[110]) cleavage systems serving as the interface. A similar response also ensues for the YBa2C3O7−δ/SrTiO3 bicrystalline superlattice with (010)[001]×[100]/( 1 ¯ 10)[001]×[110] or ( 1 ¯ 00)[001]×[010]/( 1 ¯ 10)[001]×[110] cleavage systems serving as the interface. Again, bond breaking would be smooth and continuous.
(xi) Interestingly, for the Au/MgO or YBa2C3O7/SrTiO3 superlattice, with one easy and the second one difficult, either (010)[001]×[100]/( 1 ¯ 10)[001]×[110] or ( 1 ¯ 10)[001]×[110]/(010)[001]×[100] cleavage systems serving as the interface, the SCD barrier, Δ K s * , also vanishes, and the interfacial crack would begin to propagate (in the absence of mode mixity) on the easier side of and parallel to the interface at the Griffith/Irwin critical stress intensity factor. Bond breaking is expected to be smooth and continuous on the easier side, but discontinuous and abrupt on the tougher side of the interface.
(xii) Response similar to item (xi) ensues for the YBa2C3O7−δ/SrTiO3 bicrystalline superlattice with {001}<100>×<010>/(010)[001]×[100] cleavage systems serving as the interface. Likewise, the Au/Si3N4 or YBa2C3O7/Si3N4 superlattice with one easy and the second one difficult, either (010)[001]×[100]/(001)[0 1 ¯ 0]×[ 1 ¯ 00] or (010)[001]×[100]/( 1 ¯ 00)[0 1 ¯ 0]×[001] cleavage systems serving as the interface, produces the same outcome. Again, bond breaking is expected to be smooth and continuous on the easier side, but discontinuous and abrupt on the tougher side of the interface.
(xiii) Finally, the hitherto unavailable results, pertaining to the through-thickness variations in normalized stress intensity factors for a symmetrically distributed hyperbolic sine load and its skew-symmetric counterpart that also satisfy the boundary conditions on the top and bottom surfaces of the bicrystalline superlattice plate, in the vicinity of an interfacial crack front under investigation, bridge a longstanding gap in interfacial stress singularity/fracture mechanics literature.

Funding

This research received no external funding.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Details of the Derivation of the Solution Involving Complex Roots for a (010)[001]×[100] Crack (Mode I/II Loading)

Some of the details of the mathematical derivation of the solution, involving complex roots, for an orthorhombic crystal plate, reinforced by (010) [001] through-anticrack and subjected to mode I/II loading (Section 4, case (a)), are presented here. The components of displacement that satisfy the equilibrium Equation (1) can be expressed in the following form:
u j ( r , θ , z ) = r s D b j ( z ) ( i k ) s ( cos ( θ ) + ξ ( j ) sin ( θ ) ) 2 + η ( j ) 2 sin 2 ( θ ) s / 2 { A 1 j   cos ( s ψ j )   +   A 2 j   sin ( s ψ j ) } + { ( cos ( θ ) ξ ( j ) sin ( θ ) ) 2 + η ( j ) 2 sin 2 ( θ ) } s / 2 { A 3 j   cos ( s ψ j ) + A 4 j   sin ( s ψ j ) } + O ( r s + 2 ) ,
v j ( r , θ , z ) = r s D b j z i k s cos 2 ( θ ) + ξ ( j ) + η ( j ) 2 sin 2 ( θ ) s / 2 ( H 1 j A 1 j + H 2 j A 2 j ) cos ( s ψ j ) + ( H 1 j A 2 j H 2 j A 1 j )   sin ( s ψ j ) } + { ( cos ( θ ) ξ ( j ) sin ( θ ) ) 2 + η ( j ) 2 sin 2 ( θ ) } s / 2 { ( H 1 j A 3 j + H 2 j A 4 j )   cos ( s ψ j )   ( H 1 j A 4 j H 2 j A 3 j ) sin ( s ψ j ) } + O ( r s + 2 ) ,
w j ( r , θ , z ) = O ( r s + 2 ) ,
and
σ x j ( r , θ , z ) = r s 1 D b j ( z ) ( i k ) s s { ( cos ( θ ) + ξ ( j ) sin ( θ ) ) 2 + η ( j ) 2 sin 2 ( θ ) } ( s 1 ) / 2 [ ( A 1 j { c 11 ( j ) + ( ξ ( j ) H 1 j η ( j ) H 2 j ) c 12 ( j ) } + A 2 j { ( η ( j ) H 1 j + ξ ( j ) H 2 j ) c 12 ( j ) } ) cos ( ( s 1 ) ψ j ) + ( A 1 j { ( η ( j ) H 1 j + ξ ( j ) H 2 j ) c 12 ( j ) } + A 2 j { c 11 ( j ) + ( ξ ( j ) H 1 j η ( j ) H 2 j ) c 12 ( j ) } ) sin ( ( s 1 ) ψ j ) ] + { ( cos ( θ ) ξ ( j ) sin ( θ ) ) 2 + η ( j ) 2 sin 2 ( θ ) } ( s 1 ) / 2 [ ( A 3 j { c 11 ( j ) + ( ξ ( j ) H 1 j η ( j ) H 2 j ) c 12 ( j ) } + A 4 j { ( η ( j ) H 1 j + ξ ( j ) H 2 j ) c 12 ( j ) } ) cos ( ( s 1 ) ψ j ) + ( A 3 j { ( η ( j ) H 1 j + ξ ( j ) H 2 j ) c 12 ( j ) } + A 4 j { c 11 ( j ) + ( ξ ( j ) H 1 j η ( j ) H 2 j ) c 12 ( j ) } ) sin ( ( s 1 ) ψ j ) ]
σ y j r , θ , z = r s 1 D b j z i k s s cos ( θ ) + ξ ( j ) sin ( θ ) 2 + η ( j ) 2 sin 2 ( θ ) ( s 1 ) / 2 A 1 j c 12 ( j ) + ξ ( j ) H 1 j η ( j ) H 2 j c 22 ( j ) + A 2 j η ( j ) H 1 j + ξ ( j ) H 2 j c 22 ( j ) cos s 1 ψ j + A 1 j η ( j ) H 1 j + ξ ( j ) H 2 j c 22 ( j ) + A 2 j c 12 ( j ) + ξ ( j ) H 1 j η ( j ) H 2 j c 22 ( j ) sin s 1 ψ j + cos ( θ ) ξ ( j ) sin ( θ ) 2 + η 2 sin 2 ( θ ) ( s 1 ) / 2 A 3 j c 12 ( j ) + ξ ( j ) H 1 j η ( j ) H 2 j c 22 ( j ) + A 4 j η ( j ) H 1 j + ξ ( j ) H 2 j c 22 ( j ) cos s 1 ψ j + A 3 j η ( j ) H 1 j + ξ ( j ) H 2 j c 22 ( j ) + A 4 j c 12 ( j ) + ξ ( j ) H 1 j η ( j ) H 2 j c 22 ( j ) sin s 1 ψ j
τ x y j r , θ , z = r s 1 D b 1 z i k s s c 66 ( j ) cos ( θ ) + ξ ( j ) sin ( θ ) 2 + η ( j ) 2 sin 2 ( θ ) ( s 1 ) / 2 A 1 j ξ ( j ) + H 1 j + A 2 j η ( j ) + H 2 j cos s 1 ψ j + A 1 j η ( j ) + H 2 j + A 2 j ξ ( j ) + H 1 j sin s 1 ψ j + cos ( θ ) ξ ( j ) sin ( θ ) 2 + η ( j ) 2 sin 2 ( θ ) ( s 1 ) / 2 A 3 j ξ ( j ) + H 1 j A 4 j η ( j ) + H 2 j cos s 1 ψ j + A 3 j η ( j ) + H 2 j A 4 j ξ ( j ) + H 1 j sin s 1 ψ j + O r s + 1 ,
σ z j r , θ , z = r s 1 D b j z i k s s cos ( θ ) + ξ ( j ) sin ( θ ) 2 + η ( j ) 2 sin 2 ( θ ) ( s 1 ) / 2 A 1 j c 12 ( j ) + ξ ( j ) H 1 j η ( j ) H 2 j c 23 ( j ) + A 2 j η ( j ) H 1 j + ξ ( j ) H 2 j c 23 ( j ) cos s 1 ψ j + A 1 j η ( j ) H 1 j + ξ ( j ) H 2 j c 23 ( j ) + A 2 j c 12 ( j ) + ξ ( j ) H 1 j η ( j ) H 2 j c 23 ( j ) sin s 1 ψ j + cos ( θ ) ξ ( j ) sin ( θ ) 2 + η ( j ) 2 sin 2 ( θ ) ( s 1 ) / 2 A 3 j c 12 ( j ) + ξ ( j ) H 1 j η ( j ) H 2 j c 23 ( j ) + A 4 j η ( j ) H 1 j + ξ ( j ) H 2 j c 23 ( j ) cos s 1 ψ j + A 3 j η ( j ) H 1 j + ξ ( j ) H 2 j c 23 + A 4 j c 13 ( j ) + ξ ( j ) H 1 j η ( j ) H 2 j c 23 ( j ) sin s 1 ψ j
τ x z j r , θ , z = O r s , τ y z j r , θ , z = O r s .
τ x z j r , θ , z = O r s , τ y z j r , θ , z = O r s .
where s = 1 – λ, λ being the strength or order of the stress singularity, and
D b j z =   D 1 j sin α z + D 2 j cos α z , j = 1 , 2
H 1 j = ξ ( j ) ( c 11 ( j ) c 22 ( j ) + c 66 ( j ) ) ( c 12 ( j ) + c 66 ( j ) ) , H 2 j = η ( j ) ( c 11 ( j ) c 22 ( j ) c 66 ( j ) ) ( c 12 ( j ) + c 66 ( j ) ) .
cos ( ψ j ( θ ) ) = cos ( θ ) + ξ ( j ) sin ( θ ) { ( cos ( θ ) + ξ ( j ) sin ( θ ) ) 2 + η ( j ) 2 sin 2 ( θ ) } 1 / 2 ,
sin ( ψ j ( θ ) ) = η ( j ) sin ( θ ) { ( cos ( θ ) + ξ ( j ) sin ( θ ) ) 2 + η ( j ) 2 sin 2 ( θ ) } 1 / 2 ,
cos ( ψ j ( θ ) ) = cos ( θ ) ξ ( j ) sin ( θ ) { ( cos ( θ ) ξ ( j ) sin ( θ ) ) 2 + η ( j ) 2 sin 2 ( θ ) } 1 / 2 ,
sin ( ψ ( θ ¯ ) ) = η sin ( θ ¯ ) { ( cos ( θ ¯ ) ξ sin ( θ ¯ ) ) 2 + η 2 sin 2 ( θ ¯ ) } 1 / 2 ,
A i j , i = 1, …,4, j = 1, 2, are undetermined coefficients (i.e., eigenvectors related to “stress intensity factors”). It may be noted that, since s or Re s (when s is complex) is positive, all the higher-order terms in Equation (A2) vanish as r → 0. The components of displacement can now be expressed in the cylindrical polar coordinate system as follows:
u r j ( r , θ , z ) = r s D b j z i k s cos ( θ ) + ξ ( j ) sin ( θ ) 2 + η ( j ) 2 sin 2 ( θ ) s / 2 A 1 j   cos θ   + H 1 j A 1 j + H 2 j A 2 j sin ( θ ) cos s ψ j   +   A 2 j   cos θ   + H 1 j A 2 j H 2 j A 1 j sin ( θ )   sin s ψ j + cos ( θ ) ξ ( j ) sin ( θ ) 2 + η ( j ) 2 sin 2 ( θ ) s / 2 A 3 j   cos θ   + H 1 j A 3 j + H 2 j A 4 j sin ( θ ) cos s ψ j +   A 4 j   cos θ   H 1 j A 4 j H 2 j A 3 j sin ( θ )   sin s ψ j + O r s + 2 ,
u θ j ( r , θ , z ) = r s D b j z i k s cos ( θ ) + ξ ( j ) sin ( θ ) 2 + η ( j ) 2 sin 2 ( θ ) s / 2 A 1 j   sin θ   + H 1 j A 1 j + H 2 j A 2 j cos ( θ ) cos s ψ j   +   A 2 j   sin θ   + H 1 j A 2 j H 2 j A 1 j cos ( θ )   sin s ψ j + cos ( θ ) ξ ( j ) sin ( θ ) 2 + η ( j ) 2 sin 2 ( θ ) s / 2 A 3 j   sin θ   + H 1 j A 3 j H 2 j A 4 j cos ( θ ) cos s ψ j     A 4 j   sin θ   + H 1 j A 4 j + H 2 j A 3 j cos ( θ )   sin s ψ j + O r s + 2 ,
w j ( r , θ , z ) = O r s + 1 .
Similarly, the components of the asymptotic stress field can be conveniently expressed by using the standard transformation rule:
σ r j σ θ j τ r θ j = cos 2 ( θ ) sin 2 ( θ ) sin ( 2 θ ) sin 2 ( θ ) cos 2 ( θ ) sin ( 2 θ ) 1 2 sin ( 2 θ ) 1 2 sin ( 2 θ ) cos ( 2 θ ) σ x j σ y j τ x y j ,
τ r z j τ θ z j = cos ( θ ) sin ( θ ) sin ( θ ) cos ( θ ) τ x z j τ y z j .
The stress component, σzj, is as presented in (A2d).

Appendix B. Details of the Derivation of the Solution Involving Imaginary Roots for a (010)[001]×[100] Crack (Mode I/II Loading)

This appendix provides some of the details of the mathematical derivation of the solution, involving imaginary roots, for an orthorhombic crystal plate, weakened by (010) [001] through-crack and subjected to mode I/II loading (Section 4, case (b)). The components of displacement that satisfy the equilibrium Equation (1) can be expressed in the following form:
u j ( r , θ , z ) = r s D b j z i k s cos 2 ( θ ) + ξ ( j ) + η ( j ) 2 sin 2 ( θ ) s / 2 A 1 j   cos s ψ 1 j   +   A 2 j   sin s ψ 1 j + cos 2 ( θ ¯ ) + ξ ( j ) η ( j ) 2 sin 2 ( θ ) s / 2 A 3 j   cos s ψ 1 j   +   A 4 j   sin s ψ 1 j + O r s + 2 ,
v j ( r , θ , z ) = r s D b j z i k s cos 2 ( θ ) + ξ ( j ) + η ( j ) 2 sin 2 ( θ ) s / 2 H 1 j A 2 j   cos s ψ 1 j   +   H 1 j A 1 j   sin s ψ 1 j + cos 2 ( θ ) + ξ ( j ) η ( j ) 2 sin 2 ( θ ) s / 2 H 2 j A 41   cos s ψ 1 j   +   H 2 j A 3 j   sin s ψ 1 j + O r s + 2 ,
w j ( r , θ , z ) = O r s + 1 .
and
σ x j r , θ , z = r s 1 D b j z i k s s cos 2 ( θ ) + ξ ( j ) + η ( j ) 2 sin 2 ( θ ) ( s 1 ) / 2 c 11 ( j ) + H 1 j ξ + η c 12 ( j ) A 1 j cos s 1 ψ 1 j + A 2 j sin s 1 ψ 1 j + cos 2 ( θ ) + ξ ( j ) η ( j ) 2 sin 2 ( θ ) } ( s 1 ) / 2 c 11 ( j ) + H 2 j ξ ( j ) η ( j ) c 12 ( j ) A 3 j cos s 1 ψ 1 j + A 4 j sin s 1 ψ 1 j
σ y j r , θ , z = r s 1 D b j z i k s s cos 2 ( θ ) + ξ + η 2 sin 2 ( θ ) ( s 1 ) / 2 c 12 ( j ) + H 1 j ξ ( j ) + η ( j ) c 22 ( j ) A 1 j cos s 1 ψ 1 j + A 2 j sin s 1 ψ 1 j + cos 2 ( θ ) + ξ ( j ) η ( j ) 2 sin 2 ( θ ) } ( s 1 ) / 2 c 12 ( j ) + H 2 j ξ ( j ) η ( j ) c 22 ( j ) A 3 j cos s 1 ψ 1 j + A 4 j sin s 1 ψ 1 j
τ x y j r , θ , z = r s 1 D b j z i k s s c 66 ( j ) cos 2 ( θ ) + ξ ( j ) + η ( j ) 2 sin 2 ( θ ) ( s 1 ) / 2 H 1 j ξ ( j ) + η ( j ) A 2 j   cos ( s 1 ) ψ 1 j   +   A 1 j   sin ( s 1 ) ψ 1 j + cos 2 ( θ ) + ξ ( j ) η ( j ) 2 sin 2 ( θ ) ( s 1 ) / 2 H 2 j ξ ( j ) η ( j ) A 4 j   cos ( s 1 ) ψ 1 j   +   A 3 j   sin ( s 1 ) ψ 1 j
σ z j r , θ , z = r s 1 D b j z i k s s cos 2 ( θ ) + ξ ( j ) + η ( j ) 2 sin 2 ( θ ) ( s 1 ) / 2 c 12 ( j ) + H 1 j ξ ( j ) + η ( j ) c 23 ( j ) A 1 j cos s 1 ψ 1 j + A 2 j sin s 1 ψ 1 j + cos 2 ( θ ) + ξ ( j ) η ( j ) 2 sin 2 ( θ ) } ( s 1 ) / 2 c 12 ( j ) + H 2 j ξ ( j ) η ( j ) c 23 ( j ) A 3 j cos s 1 ψ 1 j + A 4 j sin s 1 ψ 1 j
τ x z j r , θ , z = O r s , τ y z j r , θ , z = O r s .
where
H 1 j = c 11 ( j ) c 22 ( j ) c 66 ( j ) ξ ( j ) + η ( j ) 2 c 12 ( j ) + c 66 ( j ) ξ ( j ) + η ( j ) , H 2 j = c 11 ( j ) c 22 ( j ) c 66 ( j ) ξ ( j ) η ( j ) 2 c 12 ( j ) + c 66 ( j ) ξ ( j ) η ( j ) .
and
cos ψ 1 j ( θ ) = cos ( θ ) cos 2 ( θ ) + ( ξ ( j ) + η ( j ) ) 2 sin 2 ( θ ) 1 / 2 ,
sin ψ 1 j ( θ ) = ( ξ ( j ) + η ( j ) ) sin ( θ ) cos 2 ( θ ) + ( ξ ( j ) + η ( j ) ) 2 sin 2 ( θ ) 1 / 2 ,
cos ψ 1 j ( θ ) = cos ( θ ) cos 2 ( θ ) + ( ξ ( j ) η ( j ) ) 2 sin 2 ( θ ) 1 / 2 ,
sin ψ 1 j ( θ ) = ( ξ ( j ) η ( j ) ) sin ( θ ) cos 2 ( θ ) + ( ξ ( j ) η ( j ) ) 2 sin 2 ( θ ) 1 / 2 ,
It may be noted that, since s or Re s (when s is complex) is positive, all the higher-order terms in Equation (A9) vanish as r → 0. The components of displacement can now be expressed in the cylindrical polar coordinate system as follows:
u r j ( r , θ , z ) = r s D b j z i k s cos 2 ( θ ) + ξ ( j ) + η ( j ) 2 sin 2 ( θ ) s / 2 A 1 j cos ( θ ) H 1 j A 2 j sin ( θ ) cos s ψ 1 j + A 2 j cos ( θ ) + H 1 j A 1 j sin ( θ ) sin s ψ 1 j + cos 2 ( θ ) + ξ ( j ) η ( j ) 2 sin 2 ( θ ) s / 2 A 3 j cos ( θ ) H 2 j A 4 j sin ( θ ) cos s ψ 1 j + A 4 j cos ( θ ) + H 2 j A 3 j sin ( θ ) sin s ψ 1 j + O r s + 2 ,
u θ j ( r , θ , z ) = r s D b j z i k s cos 2 ( θ ) + ξ ( j ) + η ( j ) 2 sin 2 ( θ ) s / 2 A 1 j sin ( θ ) + H 1 j A 2 j cos ( θ ) cos s ψ 1 j + A 2 j sin ( θ ) + H 1 j A 1 j cos ( θ ) sin s ψ 1 j + cos 2 ( θ ) + ξ ( j ) η ( j ) 2 sin 2 ( θ ) s / 2 A 3 j sin ( θ ) + H 2 j A 4 j cos ( θ ) cos s ψ 1 j + A 4 j sin ( θ ) + H 2 j A 3 j cos ( θ ) + O r s + 2 ,
w j ( r , θ , z ) = O r s + 1 .
Similarly, the components of the asymptotic stress field can be conveniently expressed by using the standard transformation rule, provided earlier in Equation (A7). The stress component, σz, is presented in Equation (A9d).

References

  1. Chaudhuri, R.A. Three-dimensional singular stress fields near the circumferential junction corner line of an island/substrate system either free-standing or fully/partially bonded to a rigid block. Eng. Fract. Mech. 2013, 107, 80–97. [Google Scholar] [CrossRef]
  2. Lin, X.M.; Jaeger, H.M.; Sorensen, C.M.; Klabunde, K.J. Formation of long-range-ordered nanocrystal superlattices on silicon nitride substrates. J. Phys. Chem. B 2001, 105, 3353–3357. [Google Scholar] [CrossRef]
  3. Ovsyannikov, G.A.; Borisenko, I.V.; Constantinian, K.Y.; Kislinski, Y.V.; Hakhoumian, A.; Pogosyan, N.G.; Zakaryan, T.V. Josephson bicrystal junctions on sapphire substrates for THz frequency application. J. Supercond. Nov. Magn. 2006, 19, 669–673. [Google Scholar] [CrossRef]
  4. Pashley, D.W.; Stowell, M.J.; Jacobs, M.H.; Law, T.J. The growth and structure of gold and silver deposits formed by evaporation inside an electron microscope. Philos. Mag. 1964, 10, 127–158. [Google Scholar] [CrossRef]
  5. Yin, E.; Rubin, M.D.; Dixon, M.A. Sputtered YBCO Films on Metal Substrates. J. Mater. Res. 1992, 7, 1636–1640. [Google Scholar] [CrossRef]
  6. Lekhnitskii, S.G. Anisotropic Plates; Gordon and Breach: New York, NY, USA, 1968. [Google Scholar]
  7. Stroh, A.N. Dislocations and cracks in anisotropic elasticity. Philos. Mag. 1958, 7, 625–646. [Google Scholar] [CrossRef]
  8. Nazarov, S.A. Stress intensity factors and crack deviation conditions in a brittle anisotropic solid. J. Appl. Mech. Tech. Phys. 2005, 36, 386–394. [Google Scholar] [CrossRef]
  9. Rice, J.R. Elastic fracture mechanics concepts for interfacial cracks. ASME J. Appl. Mech. 1988, 55, 98–103. [Google Scholar] [CrossRef]
  10. Wu, K.-C. Stress intensity factors and energy release rate for interfacial cracks between dissimilar anisotropic materials. ASME J. Appl. Mech. 1990, 57, 882–886. [Google Scholar] [CrossRef]
  11. Qu, J.; Bassani, J.L. Cracks on bimaterial and bicrystal interfaces. J. Mech. Phys. Solids 1989, 37, 417–433. [Google Scholar] [CrossRef]
  12. Wang, T.-C.W.; Shih, C.F.; Suo, Z. Crack extension and kinking in laminates and bicrystals. Int. J. Solids Struct. 1992, 29, 327–344. [Google Scholar]
  13. 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]
  14. Xie, M.; Chaudhuri, R.A. Three-dimensional stress singularity at a bimaterial interface crack front. Compos. Struct. 1997, 40, 137–147. [Google Scholar] [CrossRef]
  15. Chaudhuri, R.A.; Xie, M. A novel eigenfunction expansion solution for three-dimensional crack problems. Compos. Sci. Technol. 2000, 60, 2565–2580. [Google Scholar] [CrossRef]
  16. Chaudhuri, R.A. On through-thickness distribution of stress intensity factors and energy release rates in the vicinity of crack fronts. Eng. Fract. Mech. 2019, 216, 106478. [Google Scholar] [CrossRef]
  17. Muskhelishvili, N.I. Some Basic Problems of the Mathematical Theory of Elasticity; Academy of Sciences of USSR: Moscow, Russia, 1954. [Google Scholar]
  18. Chaudhuri, R.A. Comparison of stress singularities of kinked carbon and glass fibers weakening compressed unidirectional composites: A three-dimensional trimaterial junction stress singularity analysis. Philos. Mag. 2014, 94, 625–667. [Google Scholar] [CrossRef]
  19. Stenger, F.; Chaudhuri, R.; Chiu, J. Sinc solution of boundary integral form for two-dimensional bi-material elasticity problems. Compos. Sci. Technol. 2000, 60, 2197–2211. [Google Scholar] [CrossRef]
  20. Dal Corso, F.; Bigoni, D. The interactions between shear bands and rigid lamellar inclusions in a ductile metal matrix. Proc. R. Soc. A 2009, 465, 143–163. [Google Scholar] [CrossRef]
  21. Chaudhuri, R.A. Eigenfunction expansion solutions for three-dimensional rigid planar inclusion problem. Int. J. Fract. 2003, 121, 95–110. [Google Scholar] [CrossRef]
  22. Yoon, J.; Chaudhuri, R.A. Three-dimensional asymptotic antiplanet shear stress fields at the front of interfacial crack/anticrack type discontinuities in trimaterial bonded plates. Compos. Struct. 2011, 93, 1505–1515. [Google Scholar] [CrossRef]
  23. Chaudhuri, R.A.; Yoon, J. Three-dimensional asymptotic mode I/II stress fields at the front of interfacial crack/anticrack discontinuities in trimaterial bonded plates. Compos. Struct. 2012, 94, 351–362. [Google Scholar] [CrossRef]
  24. Chaudhuri, R.A.; Xie, M. On three-dimensional asymptotic solution, and applicability of Saint–Venant’s principle to pie-shaped wedge and end face (of a semi-infinite plate) boundary value problems. Eng. Fract. Mech. 2015, 142, 93–107, Erratum in Eng. Fract. Mech. 2019, 217, 106506. [Google Scholar] [CrossRef]
  25. Chaudhuri, R.A. On applicability and uniqueness of the correspondence principle to pie-shaped wedge (“wedge paradox”) with various boundary conditions. Eng. Fract. Mech. 2020, 231, 106991. [Google Scholar] [CrossRef]
  26. Chaudhuri, R.A.; Xie, M. A tale of two saints: St. Venant and “St. Nick”—Does St. Venant’s principle apply to bimaterial straight edge and wedge singularity problems? Compos. Sci. Technol. 2000, 60, 2503–2515. [Google Scholar] [CrossRef]
  27. Xie, M.; Chaudhuri, R.A. Three-dimensional asymptotic stress field at the front of a bimaterial wedge of symmetric geometry under antiplane shear loading. Compos. Struct. 2001, 54, 509–514. [Google Scholar] [CrossRef]
  28. Chiu, J.S.H.; Chaudhuri, R.A. Three-dimensional asymptotic stress field at the front of an unsymmetric bimaterial pie-shaped wedge under antiplane shear loading. Compos. Struct. 2002, 58, 129–137. [Google Scholar] [CrossRef]
  29. Chaudhuri, R.A.; Xie, M. Free-edge stress singularity in a bimaterial laminate. Compos. Struct. 1997, 40, 129–136. [Google Scholar] [CrossRef]
  30. Chaudhuri, R.A.; Chiu, S.H.J. Three-Dimensional Asymptotic Stress Field in the Vicinity of an Adhesively Bonded Scarf Joint Interface. Compos. Struct. 2009, 89, 475–483. [Google Scholar] [CrossRef]
  31. Chiu, S.H.J.; Chaudhuri, R.A. A three-dimensional eigenfunction expansion approach for singular stress field near an adhesively-bonded scarf joint interface in a rigidly-encased plate. Eng. Fract. Mech. 2011, 78, 2220–2234. [Google Scholar] [CrossRef]
  32. Yoon, J.; Chaudhuri, R.A. Three-dimensional asymptotic stress fields at the front of a trimaterial junction. Compos. Struct. 2012, 94, 337–350. [Google Scholar] [CrossRef]
  33. Chaudhuri, R.A.; Chiu, S.H.J. Three-dimensional singular stress field near the interfacial bond line of a tapered jointed plate either free-standing (notch) or (fully/partially) attached to a super-rigid inclusion (antinotch). Eng. Fract. Mech. 2012, 91, 87–102. [Google Scholar] [CrossRef]
  34. Chaudhuri, R.A. Three-dimensional asymptotic stress field in the vicinity of the circumferential tip of a fiber-matrix interfacial debond. Int. J. Eng. Sci. 2004, 42, 1707–1727. [Google Scholar] [CrossRef]
  35. Chaudhuri, R.A. Three-dimensional singular stress field near a partially debonded cylindrical rigid fiber. Compos. Struct. 2006, 72, 141–150. [Google Scholar] [CrossRef]
  36. Chaudhuri, S.N.; Chaudhuri, R.A.; Benner, R.E.; Penugonda, M. Raman spectroscopy for characterization of interfacial debonds between carbon fibers and polymer matrices. Compos. Struct. 2006, 76, 375–387. [Google Scholar] [CrossRef]
  37. Chaudhuri, R.A.; Chiu, S.H.J. Three-dimensional asymptotic stress field at the front of an unsymmetric bimaterial wedge associated with matrix cracking or fiber break. Compos. Struct. 2007, 78, 254–263. [Google Scholar] [CrossRef]
  38. Chaudhuri, R.A. Three-dimensional asymptotic stress field in the vicinity of the circumference of a penny shaped discontinuity. Int. J. Solids Struct. 2003, 40, 3787–3805. [Google Scholar] [CrossRef]
  39. Kaczynski, A.; Kozlowski, W. Thermal stresses in an elastic space with a perfectly rigid flat inclusion under perpendicular heat flow. Int. J. Solids Struct. 2009, 46, 1772–1777. [Google Scholar] [CrossRef]
  40. Willis, J.R. The penny shaped crack on an interface. Q. J. Mech. Appl. Math. 1972, 25, 367–385. [Google Scholar] [CrossRef]
  41. Chaudhuri, R.A. Three-dimensional asymptotic stress field in the vicinity of the circumference of a bimaterial penny shaped interface discontinuity. Int. J. Fract. 2006, 141, 207–221. [Google Scholar] [CrossRef]
  42. Folias, E.S. The 3D stress field at the intersection of a hole and a free surface. Int. J. Fract. 1987, 35, 187–194. [Google Scholar] [CrossRef]
  43. Chaudhuri, R.A. Three-dimensional asymptotic stress field in the vicinity of the line of intersection of a circular cylindrical through/part-through open/rigidly plugged hole and a plate. Int. J. Fract. 2003, 122, 65–88. [Google Scholar] [CrossRef]
  44. Folias, E.S. On interlaminar stresses of a composite plate around the neighborhood of a hole. Int. J. Solids Struct. 1989, 25, 1193–1200. [Google Scholar] [CrossRef]
  45. Chaudhuri, R.A. An eigenfunction expansion solution for three-dimensional stress field in the vicinity of the circumferential line of intersection of a bimaterial interface and a hole. Int. J. Fract. 2004, 129, 361–384. [Google Scholar] [CrossRef]
  46. Folias, E.S. On the stress singularities at the intersection of a cylindrical inclusion with the free surface of a plate. Int. J. Fract. 1989, 39, 25–34. [Google Scholar] [CrossRef]
  47. Chaudhuri, R.A. Three-dimensional asymptotic stress field in the vicinity of the line of intersection of an inclusion and plate surface. Int. J. Fract. 2002, 117, 207–233. [Google Scholar] [CrossRef]
  48. Chaudhuri, R.A. Three-dimensional singular stress field at the front of a crack weakening a unidirectional fiber reinforced composite plate. Compos. Struct. 2011, 93, 513–527. [Google Scholar] [CrossRef]
  49. Chaudhuri, R.A. Three-dimensional singular stress field at the front of a crack and lattice crack deviation (LCD) in a cubic single crystal plate. Philos. Mag. 2010, 90, 2049–2113. [Google Scholar] [CrossRef]
  50. Chaudhuri, R.A. On three-dimensional singular stress field at the front of a planar rigid inclusion (anticrack) in an orthorhombic mono-crystalline plate. Int. J. Fract. 2012, 174, 103–126. [Google Scholar] [CrossRef]
  51. Chaudhuri, R.A. Three-dimensional mixed mode I+II+III singular stress field at the front of a (111)[2]×[10] crack weakening a diamond cubic mono-crystalline plate with crack turning and step/ridge formation. Int. J. Fract. 2014, 187, 15–49. [Google Scholar] [CrossRef]
  52. Chaudhuri, R.A. On three-dimensional singular stress/residual stress fields at the front of a crack/anticrack in an orthotropic/orthorhombic plate under anti-plane shear loading. Compos. Struct. 2010, 92, 1977–1984, Erratum in Compos. Struct. 2011, 93, 1058. [Google Scholar] [CrossRef]
  53. Chaudhuri, R.A. Three-dimensional stress/residual stress fields at crack/anticrack fronts in monoclinic plates under antiplane shear loading. Eng. Fract. Mech. 2012, 87, 16–35. [Google Scholar] [CrossRef]
  54. Yoon, J.; Chaudhuri, R.A. Three-dimensional singular antiplane shear stress fields at the fronts of interfacial crack/anticrack/contact type discontinuities in tricrystal anisotropic plates. Eng. Fract. Mech. 2013, 102, 15–31. [Google Scholar] [CrossRef]
  55. Chaudhuri, R.A.; Xie, M.; Garala, H.J. Stress singularity due to kink band weakening a unidirectional composite under compression. J. Compos. Mater. 1996, 30, 672–691. [Google Scholar] [CrossRef]
  56. Eshelby, J.D.; Read, W.T.; Shockley, W. Anisotropic elasticity with application to dislocation theory. Acta Metall. 1953, 1, 251–259. [Google Scholar] [CrossRef]
  57. Chaudhuri, R.A. Employment of fracture mechanics criteria for accurate assessment of the full set of elastic constants of orthorhombic/tetragonal mono-crystalline YBCO. Appl. Mech. 2023, 4, 585–643. [Google Scholar] [CrossRef]
  58. Carslaw, H.S. Introduction to the Theory of Fourier Series and Integrals, 3rd ed.; Dover: New York, NY, USA, 1930. [Google Scholar]
  59. Wilcox, C.H. Uniqueness theorems for displacement fields with locally finite energy in linear elastostatics. J. Elast. 1979, 9, 221–243. [Google Scholar] [CrossRef]
  60. Riddle, J.; Gumbsch, P.; Fischmeister, H.F. Cleavage Anisotropy in Tungsten Single Crystals. Phys. Rev. Lett. 1996, 76, 3594–3597. [Google Scholar] [CrossRef]
  61. Hull, D.; Beardmore, P. Velocity of propagation of cleavage cracks in tungsten. Int. J. Fract. Mech. 1966, 2, 468–487. [Google Scholar] [CrossRef]
  62. Perez, R.; Gumbsch, P. Directional anisotropy in the cleavage fracture of silicon. Phys. Rev. Lett. 2000, 84, 5347–5350. [Google Scholar] [CrossRef]
  63. Kermode, J.R.; Albaret, T.; Sherman, D.; Bernstein, N.; Gumbsch, P.; Payne, M.C.; Csányi, G.; De Vita, A. Low-speed fracture instabilities in a brittle crystal. Nature 2008, 455, 1224–1227. [Google Scholar] [CrossRef]
  64. Newnham, R.E. Structure-Property Relations; Springer: New York, NY, USA, 1975. [Google Scholar]
  65. Pauling, L. The Chemical Bond; Cornell Univ. Press: Ithaca, NY, USA, 1967. [Google Scholar]
  66. Cotton, F.A.; Wilkinson, G. Advanced Inorganic Chemistry, 4th ed.; John Wiley & Sons: New York, NY, USA, 1980. [Google Scholar]
  67. Alexandrov, I.V.; Goncharov, A.F.; Stishov, S.M. State equation and compressibility of YBa2Cu3Ox high temperature superconductor monocrystals under pressure to 20 GPa. Pis’ Ma V Zhurnal Ehksperimental’noj I Teor. Fiz. 1988, 47, 357–360. [Google Scholar]
  68. Golding, B.; Haemmerle, W.H.; Schneemeyer, L.F.; Waszczak, J.V. Gigahertz ultrasound in single crystal superconducting YBa2Cu3O7. In Proceedings of the IEEE 1988 Ultrasonics Symposium Proceedings, Chicago, IL, USA, 2–5 October 1988; Volume 2, pp. 1079–1083. [Google Scholar] [CrossRef]
  69. Saint-Paul, M.; Tholence, J.L.; Noel, H.; Levet, J.C.; Potel, M.; Gougeon, P. Ultrasound study on YBa2Cu3O7−δ and GdBa2Cu3O7−δ single crystals. Solid State Comm. 1989, 69, 1161–1163. [Google Scholar] [CrossRef]
  70. Saint-Paul, M.; Henry, J.Y. Elastic anomalies in YBa2Cu3O7−δ single crystals. Solid State Commun. 1989, 72, 685–687. [Google Scholar] [CrossRef]
  71. Jorgensen, J.D.; Pei, S.; Lightfoot, P.; Hinks, D.G.; Veal, B.W.; Dabrowski, B.; Paulikas, A.P.; Kleb, R.; Brown, I.D. Pressure-induced charge transfer and dTc/dP in YBa2Cu3O7−x” and was published in the journal. Phys. C Superconduct. 1990, 171, 93. [Google Scholar] [CrossRef]
  72. Reichardt, W.; Pintschovius, L.; Hennion, B.; Collin, F. Inelastic neutron scattering study of YBa2Cu3O7-x. Supercond. Sci. Technol. 1988, 1, 173–176. [Google Scholar] [CrossRef]
  73. Baumgart, P.; Blumenröder, S.; Erle, A.; Hillebrands, B.; Güntherodt, G.; Schmidt, H. Sound velocities of YBa2Cu3O7−δ single crystals measured by Brillouin spectroscopy. Solid State Commun. 1989, 69, 1135–1137. [Google Scholar] [CrossRef]
  74. Baumgart, P.; Blumenröder, S.; Erle, A.; Hillebrands, B.; Splittgerber, P.; Güntherodt, G.; Schmidt, H. Sound velocities of YBa2Cu3O7−δ and Bi2Sr2CaCu2Ox single crystals measured by Brillouin spectroscopy. Phys. C Supercond. Appl. 1989, 162–164, 1073–1074. [Google Scholar] [CrossRef]
  75. Zouboulis, E.; Kumar, S.; Chen, C.H.; Chan, S.K.; Grimsditch, M.; Downey, J.; McNeil, L. Surface waves on the a, b and c faces of untwinned single crystals of YBa2C3O7−δ. Phys. C Supercond. 1992, 190, 329–332. [Google Scholar] [CrossRef]
  76. Lei, M.; Sarrao, J.L.; Visscher, W.M.; Bell, T.M.; Thompson, J.D.; Migliori, A.; Welp, U.W.; Veal, B.M. Elastic constants of a monocrystal of superconducting YBa2Cu3O7-δ. Phys. Rev. B 1993, 10, 6154–6156. [Google Scholar] [CrossRef]
  77. Migliori, A.; Visscher, W.M.; Brown, S.E.; Fisk, Z.; Cheong, S.W.; Alten, B.; Ahrens, E.T.; Kubat-Martin, K.A.; Maynard, J.D.; Huang, Y.; et al. Elastic constants and specific-heat measurements on single crystals of La2CuO4. Phys. Rev. B 1990, 41, 2098–2102. [Google Scholar] [CrossRef]
  78. Nye, J.F. Physical Properties of Crystals; Oxford University Press: Oxford, UK, 1979. [Google Scholar]
  79. Weidner, D.J.; Simmons, G.J. The effect of many-body forces on the elastic properties of simple oxides and olivine explores the influence of non-central or many-body forces on the elastic properties of simple oxides and more complex structures like olivine. Geophys. Res. 1972, 77, 826–847. [Google Scholar] [CrossRef]
  80. Granozio, F.M.; di Uccio, U.S. Gibbs energy and growth habits of YBCO. J. Alloys Compd. 1997, 251, 56–64. [Google Scholar] [CrossRef]
  81. Streiffer, S.K.; Lairson, B.M.; Eom, C.B.; Clemens, B.M.; Bravman, J.C.; Geballe, T.H. Microstructure of ultrathin films of YBa2Cu3O7− δ on MgO. Phys. Rev. B 1991, 43, 13007. [Google Scholar] [CrossRef] [PubMed]
  82. Fowler, D.E.; Brundle, C.R.; Lerczac, J.; Holtzberg, F.J. Core and valence XPS spectra of clean, cleaved single crystals of YBa2Cu3O7. Electron Spectrosc. Relat. Phenom. 1990, 52, 323–339. [Google Scholar] [CrossRef]
  83. Tanaka, S.; Nakamura, T.; Tokuda, H.; Iiyama, M. All in situ deposition and characterization of YBa2Cu3O7−x thin films by low-energy electron diffraction and low-energy ion scattering spectroscopy. Appl. Phys. Lett. 1993, 62, 3040. [Google Scholar] [CrossRef]
  84. Goyal, A.; Funkenbusch, P.D.; Kroeger, D.M.; Burns, S.J. Anisotropic hardness and fracture toughness of highly aligned YBa2Cu3O7−δ. J. Appl. Phys. 1992, 71, 2363–2366. [Google Scholar] [CrossRef]
  85. Lin, C.T.; Liang, W.Y. Etch defects in YBa2Cu3O7− δ single crystals grown from flux. Phys. C 1994, 225, 275–286. [Google Scholar] [CrossRef]
  86. Cook, R.F.; Dinger, T.R.; Clarke, D.R. Fracture toughness measurements of YBa2Cu3Ox single crystals. Appl. Phys. Lett. 1987, 61, 454–456. [Google Scholar] [CrossRef]
  87. Raynes, A.S.; Freiman, S.W.; Gayle, F.W.; Kaiser, D.L. Fracture toughness of YBa2Cu3O6+δ single crystals: Anisotropy and twinning effects. J. Appl. Phys. 1991, 70, 5254–5257. [Google Scholar] [CrossRef]
  88. Campilho, R.D.S.G.; Madani, K.; Belhouari, M. (Eds.) Interfacial Crack Turning in Bi-Crystals of Cubic Symmetry: A Three-Dimensional Stress Singularity Analysis. In Fracture Mechanics: Advances in Research and Applications; Campilho, Nova Science Publishers, Inc.: New York, NY, USA, 2024; pp. 17–56, Chapter 2. [Google Scholar]
  89. Comninou, M. The interface crack. ASME J. Appl. Mech. 1977, 44, 631–636. [Google Scholar] [CrossRef]
  90. Hull, D.; Clyne, T.W. An Introduction to Composite Materials, 2nd ed.; Cambridge University Press: Cambridge, UK, 1996. [Google Scholar]
  91. Tuinstra, F.; Koenig, J.L. Characterization of graphite fiber surfaces with Raman spectroscopy. J. Compos. Mater. 1970, 4, 492–499. [Google Scholar] [CrossRef]
  92. Chaudhuri, R.A.; El-Gerwi, L.A. Estimation of residual stresses in a pressurized thin cylindrical vessel weakened by a finite line crack. In Fracture Mechanics: Advances in Research and Applications; Campilho, R.D.S.G., Madani, K., Belhouari, M., Eds.; Nova Science Publishers, Inc.: New York, NY, USA, 2024; pp. 259–288, Chapter 11. [Google Scholar]
Figure 1. Schematic of a bicrystalline superlattice plate with an interfacial crack.
Figure 1. Schematic of a bicrystalline superlattice plate with an interfacial crack.
Crystals 14 00523 g001
Figure 2. Variation in stress intensity factors through thickness for sine hyperbolic load: (a) symmetric and (b) skew-symmetric.
Figure 2. Variation in stress intensity factors through thickness for sine hyperbolic load: (a) symmetric and (b) skew-symmetric.
Crystals 14 00523 g002
Figure 3. Variation in (modes I, II, or III) energy release rate through thickness due to far-field sine hyperbolic load.
Figure 3. Variation in (modes I, II, or III) energy release rate through thickness due to far-field sine hyperbolic load.
Crystals 14 00523 g003
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Chaudhuri, R.A. Three-Dimensional Singular Stress Fields and Interfacial Crack Path Instability in Bicrystalline Superlattices of Orthorhombic/Tetragonal Symmetries. Crystals 2024, 14, 523. https://doi.org/10.3390/cryst14060523

AMA Style

Chaudhuri RA. Three-Dimensional Singular Stress Fields and Interfacial Crack Path Instability in Bicrystalline Superlattices of Orthorhombic/Tetragonal Symmetries. Crystals. 2024; 14(6):523. https://doi.org/10.3390/cryst14060523

Chicago/Turabian Style

Chaudhuri, Reaz A. 2024. "Three-Dimensional Singular Stress Fields and Interfacial Crack Path Instability in Bicrystalline Superlattices of Orthorhombic/Tetragonal Symmetries" Crystals 14, no. 6: 523. https://doi.org/10.3390/cryst14060523

APA Style

Chaudhuri, R. A. (2024). Three-Dimensional Singular Stress Fields and Interfacial Crack Path Instability in Bicrystalline Superlattices of Orthorhombic/Tetragonal Symmetries. Crystals, 14(6), 523. https://doi.org/10.3390/cryst14060523

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