# Antiplane Fracture Problem of Three Nanocracks Emanating from an Electrically Permeable Hexagonal Nanohole in One-Dimensional Hexagonal Piezoelectric Quasicrystals

Authors: Dongsheng Yang; Guanting Liu
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

Based on the Gurtin-Murdoch surface/interface model and complex potential theory, by constructing a new conformal mapping, the electrically permeable boundary condition with surface effect is established, and the antiplane fracture problem of three nanocracks emanating from a hexagonal nanohole in one-dimensional hexagonal piezoelectric quasicrystals with surface effect is studied. The exact solutions of the stress intensity factor of the phonon field and the phason field, the electric displacement intensity factor, and the energy release rate are obtained under the two electrically permeable and the electrically impermeable boundary conditions. The numerical examples show the influence of surface effect on the stress intensity factors of the phonon field and the phason field, the electric displacement intensity factor, and the energy release rate under the two boundary conditions. It can be seen that the surface effect leads to the coupling of the phonon field, phason field, and electric field, and with the decrease of cavity size, the influence of surface effect is more obvious.


---

## Body

## 1. Introduction

Quasicrystals (QCs) [1] were a new kind of structural material which was discovered in 1984. It has unique properties [2–4] such as low conductivity and good toughness. The application of this material is broad. Because defects are closely related to the safety of materials, it has attracted many scholars’ attention. Therefore, the research on QCs fracture mechanics is of great significance. The presence of the phason field makes quasicrystals different from traditional crystals. In reference [5], though the researchers studied elasticity from different descriptions, the unit-cell description based on the Penrose tiling is adopted too, but the density wave description based on the Laudau [6] phenomenological theory on symmetry-breaking of condensed matter has played the central role and been widely acknowledged. This means there are two elementary excitations of low energy, phonon u and phason w for quasicrystals, in which vector u is in the parallel space and vector w is in the perpendicular space, respectively. While w can be understood as the phason degrees of freedom in quasicrystals, which described the local rearrangement of unit-cell description based on the Penrose tiling. Li et al. [7] obtained a set of 3D general solutions to static problems of 1D hexagonal piezoelectric quasicrystals. Li et al. [8] discussed the problem of an infinite one-dimensional hexagonal quasicrystal medium weakened by a penny-shaped crack and subjected to a pair of antisymmetric and identical uniform heat fluxes. Reference [9] shows that the corresponding parameters have a significant influence on the size of the plastic zone and the distribution of normal stress. Li [10] points out that when the lips of crack with special configurations (penny-shaped, external circular, and half-infinite plane) are subjected to nonaxisymmetric external concentrated forces, the resulting phonon and phason elastic fields are expressed completely and exactly in terms of some elementary functions. It is of great significance to study various problems of quasicrystals with defects. The materials studied in this paper are quasicrystals with piezoelectric effect. At present, the research on piezoelectric materials has been more profound. Reference [11] studied new exact solutions for modified polarization saturation models for arbitrary polarized and semipermeable 2D piezoelectric materials. Sharma et al. [12] analyzed the numerical study of finite-specimen-size effects on the cracked 2D piezoelectric body. Sharma et al. [13] studied new numerical results accounted for an array of equidistant inclined cracks in a 2D piezoelectric strip by the distributed dislocation method. Bui [14] obtained numerical solutions of static and transient dynamic generalized fracture parameters of impermeable cracks in a 2D linear piezoelectric body. Bui and Zhang [15] studied transient stationary dynamic crack analysis in 2D homogeneous and linear piezoelectric materials. Reference [16] applied the XFEM and the level set method to numerical analysis of subinterface crack problems in piezoelectric bimaterials. For one-dimensional (1D) QCs antiplane fracture mechanics, the complex potential method is an important research method. A great deal of research has been done on the antiplane fracture of different defects in 1D hexagonal QCs. In references [17–21], the exact solutions of asymmetric cracks emanating from circular and elliptical hole in 1D hexagonal QCs have been obtained. In reference [22], the antiplane fracture problem of a circular hole with 2k periodic radial cracks in 1D hexagonal QCs is studied by using the complex method. In reference [23], the numerical conformal mapping is introduced to study the antiplane problem of asymmetric cracks emanating from a hexagonal hole in piezoelectric materials.In the study of fracture mechanics, if the size of hole and cracks approaches nanoscale, the ratio of surface area to volume will be very large. At this time, the influence of surface effect on the field intensity factor and energy release rate is very significant [24, 25]. By introducing surface stress, the Gurtin-Murdoch surface/interface model reasonably characterizes the phenomenon of surface effect caused by too small defect sizes. Based on this model, many scholars have done a lot of research on this kind of problem. Xioa et al. [26–28] studied the antiplane fracture problem of nanohole in piezoelectric materials and discussed the size-dependent effect of the stress intensity factor (SIF). References [29–31] increased consideration of electric field and study the antiplane problem of nanohole in the piezoelectric body under the condition of electrically impermeable boundary and the analytical solutions of the SIF, and the energy release rate at the crack tip is obtained. In reference [32], the problems of elliptical nanohole or nanocrack in piezoelectric materials under the electrically permeable boundary condition are studied. In references [33, 34], the antiplane problems of a nanocrack or elliptical nanohole in magnetoelastic materials are studied. In the past research, the complex potential method has been used to solve many problems of cracks emanating from an elliptical hole or circular hole. There are few researches on the SIFs of cracks emanating from a regular hexagonal hole, but the problem of cracks emanating from a regular hexagonal hole widely exists in engineering which has important practical and theoretical significance for this kind of research. In reference [35], the influence of the maximum stress and the stress distribution on the stress concentration of the hexagonal hole in the boat hull structure is studied. Because of the excellent seismic performance of the honeycomb members, the hexagonal honeycomb beams are widely used at present. In references [36, 37], the stress distribution characteristics of asymmetric cracks emanating from a regular hexagonal honeycomb beam are studied.This paper, like the published studies, is based on the Gurtin-Murdoch surface/interface model and complex potential theory. However, for the study of the antiplane problem of nanoscale defect in one-dimensional hexagonal piezoelectric quasicrystals, the analytical solutions under the electrically permeable boundary condition is also considered in this paper; this new study was not available in previous studies. In the published studies on antiplane fracture problems, the research on cracks emanating from a polygonal hole focuses on one crack emanating from a regular triangle hole, two or four cracks emanating from a square hole, two cracks emanating from a regular hexagonal hole. However, a new configuration of three cracks emanating from a regular hexagonal hole is studied in this paper. Inspired by references [20, 21], this paper constructs the conformal mapping from the exterior of the three cracks emanating from a regular hexagonal hole to the inside of a circle (radius R) hole. Based on the Gurtin-Murdoch surface/interface model theory and the complex potential theory [38–40], the electrically permeable boundary conditions of 1D hexagonal piezoelectric QCs with surface effect are established. By constructing a new conformal mapping and using the analytical function theory, the fracture problem of three cracks emanating from a hexagonal hole in one-dimensional hexagonal piezoelectric quasicrystals is studied for the first time at nanoscale. Under the electrically permeable boundary condition, the exact analytical solutions of the stress intensity factor and energy release rate are obtained. It has a good reference value for analyzing the change of stress intensity factor and energy release rate with material defect size.


## 2. Basic Equation

As shown in Figure1, the shape of the defect is three equal length nanocracks emanating from a hexagonal nanohole. The hexagonal side length is a, and the crack length is L, the defect penetrates along the quasiperiodic direction. It is affected by the phonon field mechanical load τzy∞, the phason field mechanical load Hzy∞, and the electric load Dy∞ at infinity. For 1D hexagonal piezoelectric QCs, we take the z-axis as the quasiperiodic direction and the electric poling direction, and the xoy plane perpendicular to the z-axis as the isotropic plane.Figure 1
Regular hexagonal nanohole with three nanocracks in 1D hexagonal piezoelectric QCs.Constitutive equations:(1)σzx=2c44εzx+R3ωzx−e15Ex,σzy=2c44εzy+R3ωzy−e15Ey,Hzx=2R3εzx+K2ωzx−d15Ex,Hzy=2R3εzy+K2ωzy−d15Ey,Dx=2e15εzx+d15ωzx+λ11Ex,Dy=2e15εzy+d15ωzy+λ11Ey.Equilibrium equations:(2)σzj,j=0,Hzj,j=0,Dj,j=0.Gradient equations:(3)εzj=12uz,j,ωzj=wz,j,Ej=−φ,j,where j=x,y in Cartesian coordinates; σzj, εzj, and uz are the stress, strain, and displacement of the phonon field, respectively; Hzj, ωzj, and wz are the stress, strain, and displacement of the phason field; Dj, Ej, and φ are the electric displacement, the electric field, and the electric potential; c44 is the elastic constant of the phonon field; K2 is the elastic constant of the phason field; R3 is the coupling constant of the phonon field and the phason field; e15 is the piezoelectric coefficient of the phonon field; d15 is the piezoelectric coefficient of the phason field; and λ11 is the dielectric permittivity.The material matrixA and the generalized displacement u are defined as follows:(4)A=c44R3e15R3K2d15e15d15−λ11,u=uzωzφT,where A≠0 and the superscript T denotes the transpose of the vector or matrix.According to equations (1)–(3), we have(5)∇2u=0,where ∇2 is a Laplace operator.According to the method of complex potential, the solutionu3, w3, and φ of equation (5) can be expressed as the real part or imaginary part of three analytical functions F1z, F2z, and F3z. It can be assumed that(6)u=uzωzφT=ReF1zF2zF3zT=ReF,where z=x+iy and ℜe represents the real part of the analytical function.SinceFizi=1,2,3 are analytical functions, we have(7)∂F∂x=F′,∂F∂y=iF′,where F′=dF/dz.According to the above results, we have(8)σzx−iσzyHzx−iHzyDx−iDyT=AF′.The above equation can be expressed in polar coordinates:(9)σzθ−iσzrHzθ−iHzrDθ−iDrT=eiθAF′.


## 3. Analytical Solutions under the Electrically Permeable Boundary Condition

The material matrix of the surface layerAS is(10)AS=c44SR3Se15SR3SK2Sd15Se15Sd15S−λ11S.Based on the Gurtin-Murdoch surface/interface model theory [24, 25], the boundary conditions of 1D hexagonal piezoelectric QCs with defect at nanoscale are as follows:(11)uzct=uzmt,wzct=wzmt,φct=φmt,t∈S,(12)−σzrmtHzrmtDrmt−DrctT=ASρεzθ,θSωzθ,θS−Eθ,θST,t∈S,where the subscripts c, S, and m refer to quantities for the hexagonal nanohole, surface layer, and matrix, respectively.If we consider the dielectric permeability of the medium inside the hexagonal nanohole, the constitutive and equilibrium equations can be reduced to, respectively,(13)Dic=−λ11c∂φ,ic,i=x,y,(14)∇2φc=0,where λ11c denotes the dielectric permittivity of the medium such as air inside the nanohole. The solution φc of equation (14) can be expressed as the real part of the analytical function F3cz.According to reference [38], the boundary condition in the hole is(15)σzycHzycDycT=00Dy∞−d0T,where d0=e15K2−d15R3τzy∞+d15c44−e15R3Hzy∞/c44K2−R32.In order to solve this boundary valued problem, inspired by reference [20], the following conformal mapping is constructed. The exterior of the three nanocracks emanating from a regular hexagonal nanohole in the z plane is mapped to the interior of the circular nanohole (radius R) in the ζ plane. The detailed construction process is shown in Appendix section.(16)z=ωζ=Rμζ+115μζ−5+199μζ−11+11377μζ−17,μζ=141/3ζRζ3−R32+lζ3+R32+1+lζ3+R31+lζ6+2l−3ζ3R3+1+lR61/3,where R≈0.9258a, l=1+c3+1+c−3/2. The positive real parameter c is determined by the following formula:(17)L+a=R1+c+1151+c−5+1991+c−11+113771+c−17.By expanding the analytical functionsFizi=1,2,3 in the z plane,(18)Fiζ=ai∗lnζ+∑−∞+∞aikζk,i=1,2,3,where ai∗ and aik are the undetermined constants.The solutions of the problem can be obtained by selecting the finite term.(19)F3cζ=A1ζ,(20)Fmζ=F1mζF2mζF3mζ=B1C1D1ζ+B−1C−1D−1ζ−1,where A1, B1, C1, D1, B−1, C−1, and D−1 are complex constants.According to equations (8), (13), and (20), we have(21)B1C1D1=−i1+l1/321/3A−1τzy∞Hzy∞dy∞.According to the boundary condition equations (15) and (19), the following result can be obtained by the complex method [38]:(22)A1=i21+l1/321/3λ11cDy∞−d0.According to boundary condition equation (12), we have(23)ASRB1C1D1−1R2B−1C−1D−1=AB1C1D1+1R2B−1C−1D−1−00−λ11cA1.According to equations (21)–(23),(24)B−1D−1F−1=−i1+l1/321/3R2ASR+A−1ASR−AA−1τzy∞Hzy∞Dy∞+200Dy∞−d0.Substituting equations (20), (21), and (24) into equation (8),(25)σzyHzyDy=1+l1/321/3ω′ζE−R2ζ2AASR+A−1ASR−AA−1τzy∞Hzy∞Dy∞−R2ζ2f,where E is a third-order unit matrix and f=2AAS/R+A−100Dy∞−d0.


## 4. Field Intensity Factor and Energy Release Rate

### 4.1. Field Intensity Factor

Define the SIFs of the phonon field and phason field and EDIF in thez plane as [5](26)KIIIσKIIIHKIIID=limz⟶z02πz−z0σzyHzyDy,where z0=L+a.


#### 4.1.1. Electrically Permeable Problem

The phonon field stress, phason field stress, and electric displacement at the crack tip are as follows:(27)σzyHzyDy=1+l1/321/3ω′ζE−AASR+A−1ASR−AA−1τzy∞Hzy∞Dy∞−f.Substituting equation (27) into equation (26),(28)KIIIσKIIIHKIIID=π1+l1/321/3ω″RE−AASa+A−1ASa−AA−1τzy∞Hzy∞Dy∞−f.Suppose that equation (28) can be expressed as(29)KIIIσKIIIHKIIID=π1+l1/321/3ω″RKτ∗KH∗KD∗,where(30)Kτ∗KH∗KD∗=E−AASa+A−1ASa−AA−1τzy∞Hzy∞Dy∞−f.When the surface effect is neglected, i.e.,AS=0, equation (29) degenerates to the classical field intensity factor of 1D hexagonal piezoelectric QCs under the electrically permeable boundary condition:(31)KIIIσKIIIHKIIID=2π1+l1/321/3ω″Rτzy∞Hzy∞e15K2−q15R3τzy∞+q15c44−e15R3Hzy∞c44K2−R32.Except for the difference of conformal mapping, these results are consistent with reference [38] (under the electrically permeable boundary condition).


#### 4.1.2. Electrically Impermeable Problem

Due toD∞−d0=0 under the electrically impermeable boundary condition, equation (28) degenerates to(32)KIIIσKIIIHKIIID=π1+l1/321/3ω″RE−AASa+A−1ASa−AA−1τzy∞Hzy∞Dy∞.WhenAS=0, equation (32) degenerates to the classical field intensity factor of 1D hexagonal piezoelectric QCs under the electrically impermeable boundary condition,(33)KIIIσKIIIHKIIID=2π1+l1/321/3ω″Rτzy∞Hzy∞Dy∞.Equation (33) can be rewritten as(34)KIIIσKIIIHKIIID=πL′Kτzy∞Hzy∞Dy∞.Except for the difference of conformal mapping, this result is consistent with reference [38] (under the electrically impermeable boundary condition), where L′ is the equivalent crack length, and L′=L+R. K is the dimensionless field intensity factor,(35)K=2π1+l1/321/3L′ω″R.


### 4.2. Energy Release Rate

According to reference [5, 40], it is assumed that the energy release rate can be expressed as follows:(36)G=KIIIσKIIIε+KIIIHKIIIω+KIIIDKIIIE2,where KIIIε is the phonon field strain intensity factor, KIIIω is the phason field strain intensity factor , and KIIIE is the electric field intensity factor determined by the following formula:(37)KIIIεKIIIωKIIIE=A−1KIIIσKIIIHKIIID.According to equations (36) and (37), we have(38)G=12KIIIσKIIIHKIIIDA−1KIIIσKIIIHKIIID.


#### 4.2.1. Electrically Permeable Problem

Substituting equation (28) into equation (38), the energy release rate of three nanocracks emanating from an electrically permeable regular hexagonal nanohole in 1D hexagonal piezoelectric QCs with surface effect can be expressed as(39)G=π1+l2/32×22/3ω″RBτzy∞Hzy∞Dy∞−fTA−1Bτzy∞Hzy∞Dy∞−f,where B=E−AAS/R+A−1AS/R−AA−1.In particular, whenAS=0, the energy release rate of three nanocracks emanating from an electrically permeable regular hexagonal nanohole in 1D hexagonal piezoelectric QCs without surface effect can be obtained by substituting equation (31) into equation (38)(40)G=2π1+l2/322/3ω″Rc44K2−R32Γ.Except for the difference of conformal mapping, these results are consistent with reference [40] (under the electrically permeable boundary condition), where(41)Γ=−K2e152K2+c44d152+K2λ11−R32d15e15+λ11R3τzy∞2−c44e152K2+c44d152+K2λ11−R32d15e15+λ11R3Hzy∞2+2R3e152K2+c44d152+K2λ11−R32d15e15+λ11R3τzy∞Hzy∞.


#### 4.2.2. Electrically Impermeable Problem

Substituting equation (32) into equation (38), the energy release rate of three nanocracks emanating from an electrically impermeable regular hexagonal nanohole in 1D hexagonal piezoelectric QCs with surface effect can be expressed as(42)G=π1+l2/32×22/3ω″Rτzy∞Hzy∞Dy∞TBTA−1Bτzy∞Hzy∞Dy∞.In particular, whenAS=0, the energy release rate of three nanocracks emanating from an electrically impermeable regular hexagonal nanohole in 1D hexagonal piezoelectric QCs without surface effect can be obtained by substituting equation (33) into equation (38)(43)G=2π1+l2/322/3ω″RdetAΛ.Except for the difference of conformal mapping, these results are consistent with reference [40] (under the electrically impermeable boundary condition), where(44)Λ=d152+K2λ11τzy∞2+e152+c44λ11Hzy∞2+R32−K2c44Dy∞2+2e15K2−d15R3τzy∞Hzy∞+2c44d15−e15R3Hzy∞Dy∞−d15e15+λ11R3τzy∞Hzy∞.


## 4.1. Field Intensity Factor

Define the SIFs of the phonon field and phason field and EDIF in thez plane as [5](26)KIIIσKIIIHKIIID=limz⟶z02πz−z0σzyHzyDy,where z0=L+a.


### 4.1.1. Electrically Permeable Problem

The phonon field stress, phason field stress, and electric displacement at the crack tip are as follows:(27)σzyHzyDy=1+l1/321/3ω′ζE−AASR+A−1ASR−AA−1τzy∞Hzy∞Dy∞−f.Substituting equation (27) into equation (26),(28)KIIIσKIIIHKIIID=π1+l1/321/3ω″RE−AASa+A−1ASa−AA−1τzy∞Hzy∞Dy∞−f.Suppose that equation (28) can be expressed as(29)KIIIσKIIIHKIIID=π1+l1/321/3ω″RKτ∗KH∗KD∗,where(30)Kτ∗KH∗KD∗=E−AASa+A−1ASa−AA−1τzy∞Hzy∞Dy∞−f.When the surface effect is neglected, i.e.,AS=0, equation (29) degenerates to the classical field intensity factor of 1D hexagonal piezoelectric QCs under the electrically permeable boundary condition:(31)KIIIσKIIIHKIIID=2π1+l1/321/3ω″Rτzy∞Hzy∞e15K2−q15R3τzy∞+q15c44−e15R3Hzy∞c44K2−R32.Except for the difference of conformal mapping, these results are consistent with reference [38] (under the electrically permeable boundary condition).


### 4.1.2. Electrically Impermeable Problem

Due toD∞−d0=0 under the electrically impermeable boundary condition, equation (28) degenerates to(32)KIIIσKIIIHKIIID=π1+l1/321/3ω″RE−AASa+A−1ASa−AA−1τzy∞Hzy∞Dy∞.WhenAS=0, equation (32) degenerates to the classical field intensity factor of 1D hexagonal piezoelectric QCs under the electrically impermeable boundary condition,(33)KIIIσKIIIHKIIID=2π1+l1/321/3ω″Rτzy∞Hzy∞Dy∞.Equation (33) can be rewritten as(34)KIIIσKIIIHKIIID=πL′Kτzy∞Hzy∞Dy∞.Except for the difference of conformal mapping, this result is consistent with reference [38] (under the electrically impermeable boundary condition), where L′ is the equivalent crack length, and L′=L+R. K is the dimensionless field intensity factor,(35)K=2π1+l1/321/3L′ω″R.


## 4.1.1. Electrically Permeable Problem

The phonon field stress, phason field stress, and electric displacement at the crack tip are as follows:(27)σzyHzyDy=1+l1/321/3ω′ζE−AASR+A−1ASR−AA−1τzy∞Hzy∞Dy∞−f.Substituting equation (27) into equation (26),(28)KIIIσKIIIHKIIID=π1+l1/321/3ω″RE−AASa+A−1ASa−AA−1τzy∞Hzy∞Dy∞−f.Suppose that equation (28) can be expressed as(29)KIIIσKIIIHKIIID=π1+l1/321/3ω″RKτ∗KH∗KD∗,where(30)Kτ∗KH∗KD∗=E−AASa+A−1ASa−AA−1τzy∞Hzy∞Dy∞−f.When the surface effect is neglected, i.e.,AS=0, equation (29) degenerates to the classical field intensity factor of 1D hexagonal piezoelectric QCs under the electrically permeable boundary condition:(31)KIIIσKIIIHKIIID=2π1+l1/321/3ω″Rτzy∞Hzy∞e15K2−q15R3τzy∞+q15c44−e15R3Hzy∞c44K2−R32.Except for the difference of conformal mapping, these results are consistent with reference [38] (under the electrically permeable boundary condition).


## 4.1.2. Electrically Impermeable Problem

Due toD∞−d0=0 under the electrically impermeable boundary condition, equation (28) degenerates to(32)KIIIσKIIIHKIIID=π1+l1/321/3ω″RE−AASa+A−1ASa−AA−1τzy∞Hzy∞Dy∞.WhenAS=0, equation (32) degenerates to the classical field intensity factor of 1D hexagonal piezoelectric QCs under the electrically impermeable boundary condition,(33)KIIIσKIIIHKIIID=2π1+l1/321/3ω″Rτzy∞Hzy∞Dy∞.Equation (33) can be rewritten as(34)KIIIσKIIIHKIIID=πL′Kτzy∞Hzy∞Dy∞.Except for the difference of conformal mapping, this result is consistent with reference [38] (under the electrically impermeable boundary condition), where L′ is the equivalent crack length, and L′=L+R. K is the dimensionless field intensity factor,(35)K=2π1+l1/321/3L′ω″R.


## 4.2. Energy Release Rate

According to reference [5, 40], it is assumed that the energy release rate can be expressed as follows:(36)G=KIIIσKIIIε+KIIIHKIIIω+KIIIDKIIIE2,where KIIIε is the phonon field strain intensity factor, KIIIω is the phason field strain intensity factor , and KIIIE is the electric field intensity factor determined by the following formula:(37)KIIIεKIIIωKIIIE=A−1KIIIσKIIIHKIIID.According to equations (36) and (37), we have(38)G=12KIIIσKIIIHKIIIDA−1KIIIσKIIIHKIIID.


### 4.2.1. Electrically Permeable Problem

Substituting equation (28) into equation (38), the energy release rate of three nanocracks emanating from an electrically permeable regular hexagonal nanohole in 1D hexagonal piezoelectric QCs with surface effect can be expressed as(39)G=π1+l2/32×22/3ω″RBτzy∞Hzy∞Dy∞−fTA−1Bτzy∞Hzy∞Dy∞−f,where B=E−AAS/R+A−1AS/R−AA−1.In particular, whenAS=0, the energy release rate of three nanocracks emanating from an electrically permeable regular hexagonal nanohole in 1D hexagonal piezoelectric QCs without surface effect can be obtained by substituting equation (31) into equation (38)(40)G=2π1+l2/322/3ω″Rc44K2−R32Γ.Except for the difference of conformal mapping, these results are consistent with reference [40] (under the electrically permeable boundary condition), where(41)Γ=−K2e152K2+c44d152+K2λ11−R32d15e15+λ11R3τzy∞2−c44e152K2+c44d152+K2λ11−R32d15e15+λ11R3Hzy∞2+2R3e152K2+c44d152+K2λ11−R32d15e15+λ11R3τzy∞Hzy∞.


### 4.2.2. Electrically Impermeable Problem

Substituting equation (32) into equation (38), the energy release rate of three nanocracks emanating from an electrically impermeable regular hexagonal nanohole in 1D hexagonal piezoelectric QCs with surface effect can be expressed as(42)G=π1+l2/32×22/3ω″Rτzy∞Hzy∞Dy∞TBTA−1Bτzy∞Hzy∞Dy∞.In particular, whenAS=0, the energy release rate of three nanocracks emanating from an electrically impermeable regular hexagonal nanohole in 1D hexagonal piezoelectric QCs without surface effect can be obtained by substituting equation (33) into equation (38)(43)G=2π1+l2/322/3ω″RdetAΛ.Except for the difference of conformal mapping, these results are consistent with reference [40] (under the electrically impermeable boundary condition), where(44)Λ=d152+K2λ11τzy∞2+e152+c44λ11Hzy∞2+R32−K2c44Dy∞2+2e15K2−d15R3τzy∞Hzy∞+2c44d15−e15R3Hzy∞Dy∞−d15e15+λ11R3τzy∞Hzy∞.


## 4.2.1. Electrically Permeable Problem

Substituting equation (28) into equation (38), the energy release rate of three nanocracks emanating from an electrically permeable regular hexagonal nanohole in 1D hexagonal piezoelectric QCs with surface effect can be expressed as(39)G=π1+l2/32×22/3ω″RBτzy∞Hzy∞Dy∞−fTA−1Bτzy∞Hzy∞Dy∞−f,where B=E−AAS/R+A−1AS/R−AA−1.In particular, whenAS=0, the energy release rate of three nanocracks emanating from an electrically permeable regular hexagonal nanohole in 1D hexagonal piezoelectric QCs without surface effect can be obtained by substituting equation (31) into equation (38)(40)G=2π1+l2/322/3ω″Rc44K2−R32Γ.Except for the difference of conformal mapping, these results are consistent with reference [40] (under the electrically permeable boundary condition), where(41)Γ=−K2e152K2+c44d152+K2λ11−R32d15e15+λ11R3τzy∞2−c44e152K2+c44d152+K2λ11−R32d15e15+λ11R3Hzy∞2+2R3e152K2+c44d152+K2λ11−R32d15e15+λ11R3τzy∞Hzy∞.


## 4.2.2. Electrically Impermeable Problem

Substituting equation (32) into equation (38), the energy release rate of three nanocracks emanating from an electrically impermeable regular hexagonal nanohole in 1D hexagonal piezoelectric QCs with surface effect can be expressed as(42)G=π1+l2/32×22/3ω″Rτzy∞Hzy∞Dy∞TBTA−1Bτzy∞Hzy∞Dy∞.In particular, whenAS=0, the energy release rate of three nanocracks emanating from an electrically impermeable regular hexagonal nanohole in 1D hexagonal piezoelectric QCs without surface effect can be obtained by substituting equation (33) into equation (38)(43)G=2π1+l2/322/3ω″RdetAΛ.Except for the difference of conformal mapping, these results are consistent with reference [40] (under the electrically impermeable boundary condition), where(44)Λ=d152+K2λ11τzy∞2+e152+c44λ11Hzy∞2+R32−K2c44Dy∞2+2e15K2−d15R3τzy∞Hzy∞+2c44d15−e15R3Hzy∞Dy∞−d15e15+λ11R3τzy∞Hzy∞.


## 5. Numerical Analysis

A particular 1D hexagonal QCs with piezoelectric effects are considered in reference [7], and the material constants are c44=50GPa, R3=1.2GPa, K2=0.3GPa, e15=−0.138C/m2, d15=−0.160C/m2, and λ11=82.6×10−12C2N/m2. The literature also points out that although the material constants are not verified by experiment, however, the positive definiteness has been checked. So, we use these material constants to calculate numerical results. In addition, the surface property in general can be obtained by atomistic simulations [40]; we choose c44S=62.5N/m, R3S=1.5N/m, K2S=0.5N/m, e15S=1.25×10−8C/m, d15S=0.62×10−8C/m, and λ11S=0 as an approximation according to reasonable estimates for some metal surfaces.Figures2–4 show the variations of the Kτ∗ of electrically permeable and impermeable nanohole with a under the phonon field and phason field mechanical and electrical loads, respectively. It is important that the result of the electrically permeable nanohole is obviously different from that of the electrically impermeable nanohole. It can be observed from Figure 2 that the phonon field SIF of the electrically permeable nanohole under phonon field mechanical loads is always larger than that of the electrically impermeable one. In addition, the phonon field SIFs both of electrically permeable and impermeable nanoholes always increase with increasing size of a. From Figures 3 and 4, it is found that the value of Kτ∗ induced obviously by Hzy∞ and Dy∞ with surface effect. When a=0nm, the influence of surface effect is very small (Kτ∗/Dy∞ tends to be 0). When the side length a increases gradually at nanoscale, the ratio of surface area to volume of defect becomes very large and then decreases gradually. At this time, Kτ∗/Dy∞ increases rapidly with the increase of a and then decreases gradually with the increase of a. According to the analysis of the analytical solution, it can be seen that the electrical loads also have influence on the SIFs of phonon field with surface effect under the lectrically permeable boundary condition, but the influence is very small. This means that the surface effect can result in the coupling among electric, magnetic, and mechanical fields.Figure 2
Variations ofKτ∗/τzy∞ with a.Figure 3
Variations ofKτ∗/Hzy∞ with a.Figure 4
Variations ofKτ∗/Dzy∞ with a.Figures5–7 show the surface effect on the phason field SIF induced by three kinds of loads in 1D hexagonal QCs. When a=0nm, the influence of surface effect is very small (KH∗/Dy∞ tends to be 0). The ratio of the surface area to the volume of the defect becomes very large when the a increases gradually at nanoscale; KH∗/Dy∞ increases rapidly with the increase of a. Then, as the ratio of area to volume decreases gradually, the value of KH∗/Dy∞ decreases with the increase of a. It can be observed that the phason field SIF under the electrically permeable and electrically impermeable boundary conditions with surface effect depend on the τzy∞,Hzy∞, and Dy∞ and are size-dependent which is different from the classical elasticity result. It can be seen from the figures that with the increase of a, the influence of the surface effect under the electric permeable and electric impermeable boundary conditions tends to the result of classical elastic theory. Obviously different from the classical elasticity theory, Dy∞ has an effect on KH∗, but it is very small.Figure 5
Variations ofKH∗/τzy∞ with a.Figure 6
Variations ofKH∗/Hzy∞ with a.Figure 7
Variations ofKH∗/Dy∞ with a.Figures8–10 show the surface effect on the EDIF induced by three kinds of loads in 1D hexagonal QCs. It can be observed that the EDIF under the electrically permeable and electrically impermeable boundary conditions with surface effect are also affected by the τzy∞,Hzy∞, and Dy∞. According to equation (30), if the surface effect is ignored, the KD∗/τzy∞ under the two boundary conditions is equal to 2e15K2−q15R3/c44K2−R32 and 2, respectively, and the KD∗/τzy∞ under the two boundary conditions is equal to 2q15c44−e15R3/c44K2−R32 and 2, respectively. It can be seen from the figures that with the continuous increase of a, the influence of surface effect becomes smaller and smaller. Dy∞ has an effect on KD∗, but it is very small.Figure 8
Variations ofKD∗/τzy∞ with a.Figure 9
Variations ofKD∗/Hzy∞ with a.Figure 10
Variations ofKD∗/Dy∞ with a.The presence of the phason field makes quasicrystals different from traditional crystals, Figures11–13 analyze the effect of phason field mechanical load on KIIIσ, KIIIH, and KIIID. Figure 11 shows that under two boundary conditions, KIIIσ decreases with the increase of Hzy∞. In other words, the increase of the mechanical load of the phason field can restrain the increase of the stress intensity factor of the phonon field, but the restrain degree is relatively low under the electrically permeable boundary condition. Figure 12 shows that under two boundary conditions, KIIIH decreases with the increase of Hzy∞. Figure 13 shows that KIIID decreases with the increase of Hzy∞ under the electrically impermeable boundary condition (at this time, the increase of electric displacement intensity factor will be restrained) and increases with the increase of Hzy∞ under the electrically permeable boundary condition (at this time, the increase of electric displacement intensity factor will be promoted.Figure 11
Variations ofKIIIσ with Hzy.Figure 12
Variations ofKIIIH with Hzy.Figure 13
Variations ofKIIID with Hzy.It can be observed from Figure14, the influence of crack length on the dimensional field intensity factor without surface effect. It can be seen from this figure that the increase of the crack length at nanoscale will promote the crack growth, and the value of K at R=2nm is greater than the value of K at R=4nm, and the value of K at R=4nm is greater than the value of K at R=6nm.Figure 14
Variations ofK with L without surface effect.According to the analysis of the classical 1D hexagonal piezoelectric QCs fracture theory [38], the value of G/G0 under the two boundary conditions without surface effect is almost the same, where G0 represents the energy release rate under the boundary condition without surface effect. So, this paper only analyzes the change of G/G0 with R under the electrically permeable boundary condition without surface effect and the change curve of dimensionless energy release rate with R under the electrically permeable and impermeable conditions with surface effect in Figure 15. It can be seen from the figure that when R increases from 0 to 300, the value of G/G0 under the electrically permeable and impermeable conditions with surface effect gradually increases from 0 to 1, respectively. In other words, with the increase of R, the influence of the surface effect on the dimensionless energy release rate decreases.Figure 15
Variations ofG/G0 with R under two electrical boundary conditions.Figure16 shows the variations of G/G0 with the phonon field mechanical loads τzy∞ under the electrically permeable boundary condition when R=5nm. With the increase of the mechanical loads of the phonon field, G/G0 first decreases and then increases and finally tends to be stable. At the same time, when the phonon field mechanical loads are small (less than 1 MPa) or large, the value of G/G0 will be very close when the phason field mechanical loads are 0.2 MPa, 0.4 MPa, and 0.6 MPa. Figure 17 shows the variations of G/G0 with the phason field mechanical loads Hzy∞ under the electrically permeable boundary conditions when R=5nm. With the increase of the phason field mechanical loads, G/G0 first increases and then decreases and finally tends to be stable. It is also shown that smaller (less than 0.05 MPa) or larger (more than 1 MPa) phason field mechanical loads will shield the influence of the electric loads (Figure 18).Figure 16
Variations ofG/G0 with τzy∞ under the electrically permeable boundary condition with surface effect.Figure 17
Variations ofG/G0 with Hzy∞ under the electrically permeable boundary condition with surface effect.Figure 18
Conformal mapping from the outside region of the regular hexagonal hole and cracks to the interior of a circle (radiusR). (a) z plane. (b) z1 plane. (c) z1 plane 0<θ2π/3. (d) z2 plane. (e) z3 plane. (f) z4 plane. (g) z5 plane. (h) z6 plane. (i) ζ1 plane. (j) ζ plane.
(a)(b)(c)(d)(e)(f)(g)(h)(i)(j)


## 6. Conclusion

Based on the Gurtin-Murdoch surface/interface model and the theory of complex potential, the antiplane fracture problem of three nanocracks emanating from a electrically permeable regular hexagonal nanohole in hexagonal piezoelectric QCs is studied. By constructing a new conformal mapping, the electrically permeable and impermeable boundary conditions with surface effect are established. The exact solutions of the SIFs of the phonon field and phason field, the EDIF, and the energy release rate with or without surface effect are obtained. The conclusions are as follows:(1)
At nanoscale, the SIFs of the phonon field and phason field and EDIF are all affected by the coupling of far-field phonon field mechanical loads, phason field mechanical loads, and electric loads, which is obviously different from the results of 1D hexagonal piezoelectric QCs without surface effect.(2)
At nanoscale, the effect of the surface effect is more obvious when the side length of the regular hexagon is smaller. When the size of defect increases to a certain extent, the final result tends to classical fracture theory.(3)
With increase of crack length, the dimensionless field intensity factor increases, and the value ofK at R=4nmis greater than the value of K at R=6nm.(4)
At nanoscale, the dimensionless energy release rate decreases first and then increases with the increase of the phonon field mechanical loads with surface effect, and the final one tends to be stable. At the same time, the effect of the phonon field mechanical loads on the dimensionless energy release rate will be shielded when the phason field mechanical loads are small or large.



---

*Source: 1372474-2020-10-12.xml*