# Singularity of Electromechanical Coupling Field in Piezoelectric Composite Junctions

Authors: Li Jiang; Renyu Ge; Jinlun Zhang
Journal: Mathematical Problems in Engineering
 (2020)

---

## Abstract

The double singularities including singular stress field and singular electric displacement field, in the tips of piezoelectric composite junctions, are analyzed by the interpolation matrix method (IMM). The double singularity analysis problem of piezoelectric composite junctions is converted into eigenvalue solution problem of ordinary differential equations with variable coefficients under corresponding boundary conditions. In numerical examples, the first couple of singularity orders and the corresponding characteristic angular functions of displacement and electric potential for the electromechanical coupling field are obtained and comparisons are presented to validate the accuracy of the proposed method. The singularity of the electromechanical coupling field at the tip of piezoelectric composite material junctions is closely related to the bonding angle and fiber direction. According to the numerical results, the best scheme can be configured for the combination of dissimilar materials.


---

## Body

## 1. Introduction

Piezoelectric devices such as brakes and sensors are composed of piezoelectric materials, conductors, and composite materials. The problem of material and geometric discontinuities is common at the junctions and interface edge of piezoelectric materials, conductors, and composite materials, which leads to the singularity of stress and electric displacement. The singular field at the tips of piezoelectric composite junctions directly affects the safety and engineering structures [1, 2]. Therefore, it is of great significance to study the singularity of the notch to ensure the reliability and safety of the structure in service.Since Williams [3] first gave the classical solution of stress field of interface crack in isotropic bimaterials, many research studies have been performed on the singularity in piezoelectric composite material wedges. The failure behavior and failure criteria of deep notch of piezoelectric ceramics were experimentally analyzed by Zhang et al. [4]. Liu and Chue [5] found that for some special notch angles of composite materials, the degree of stress singularity can be reduced or even disappeared by selecting the appropriate fiber direction. However, the experimental method is limited to obtaining at most the first two singular orders. In order to better understand the notch singularity, some numerical methods are applied to the research of notch singularity. Delale and Erdogan [6] and Chen and Harada [7] studied the stress singularity at the tip of cracks perpendicular to the interface of combined anisotropic materials and proposed the singular characteristic equations. Zhou et al. [8, 9] adopted the Hamiltonian transformation method to analytically determine the stress and electrical strength factors of edge cracks for type III piezoelectric materials and investigate the Mode-III interface V-shape notch in a finite-size one-dimensional hexagonal quasicrystalline bimaterial with piezoelectric effect in succession. In addition, Zhou et al. [10, 11] presented a symplectic approach based on the Hamiltonian system to analyze the electroelastic singularities and intensity factors of an interface crack in piezoelectric-elastic bimaterials and also to find the complex stress intensity factors and T-stress at an edge bimaterial interface crack.Wippler and Kuna [12] used the boundary element method to analyze the singular electroelastic field at the tip of cracks of piezoelectric materials and pointed out that the problem of notch singularity is more complicated than that of crack. Stro transform was applied by Shen et al. [13] for the study of notch singularity of three-phase piezoelectric materials. Xu and Rajapakse [14] used the complex potential function method to investigate the singular electroelastic field at the tip of junctions of conductors, composite materials, and piezoelectric materials under plane deformation, and the influencing parameters including material properties, binding angle, and polarization angle were considered. But there are some defects in the analysis of interface crack by the complex potential function method. Chen et al. [15] made use of a modified one-dimensional finite element form to calculate the singular field of piezoelectric composite material junctions and that of the piezoelectric material and conductor junction. Li et al. [16] presented a semianalytical technique based on the scaled boundary finite element method to analyze two-dimensional cracks and notches inside piezoelectric composites. Luangarpa and Koguchi [17] extended a conservative integral based on the Betti reciprocal principle for calculating intensities of singularity at tips of interfaces in three-dimensional piezoelectric bonded joints. Islam et al. [18] evaluated the singularity orders in three-dimensional two-phase transversely isotropic piezoelectric dissimilar joints by the eigenanalysis method. Hrstka et al. [19] used the expanded Lekhnitskii-Eshelby-Stroh formalism to investigate an asymptotic in-plane problem of bimaterial sharp notches with various geometries and interface cracks in several generally monoclinic piezoelectric bimaterials.In view of the above methods, transcendental equations are generally needed to be solved, and some of them cannot be degenerated to crack problems. Therefore, a new numerical method which adopts the interpolation matrix method is proposed to analyze the singularity at the tip of the notch. The displacement field at the tip of piezoelectric composite material junctions was firstly asymptotically expanded by power series, so that the electromechanical boundary/interface conditions of junctions were transformed into a combination expression of singularity-order characteristic equation. Based on the stress equilibrium equation and Maxwell equation, the characteristic equation about the singularity order of piezoelectric composite material junctions was derived. The developed IMM program [20–23] was used to solve the formed characteristic ordinary differential equations, which can calculate each singularity order and the corresponding characteristic angular function of displacement and electric potential for piezoelectric composite material junctions at one time. The present method has the characteristic of strong adaptability, high accuracy, and small calculation.


## 2. Singularity in Wedge Notch of Piezoelectric Materials

For the planar wedge notch of piezoelectric materials as shown in Figure1, two rectangular coordinate systems (oxz and ox′z′) and one polar coordinate system orθ are defined and the angle between the z′ axis and z axis is set as β.Figure 1
A piezoelectric wedge notch.In the rectangular coordinate systemoxz, the constitutive relations for piezoelectric material domain (PZT) are expressed as follows:(1)σxxσzzσxzDxDz=c^11c^1300−e^31c^13c^3300−e^3300c^44−e^15000e^15k^110e^31e^3300k^33εxxεzzεxzExEz,where c^11, c^13, c^33, and c^44 are the elastic modulus, e^15, e^31, and e^33 are the piezoelectric constant, and k^11 and k^33 are the dielectric constant. εij and σij ij=xx,zz,xz are the strain and stress, respectively. Ei and Di are the electric field and electric displacement, respectively.In the polar coordinate systemorθ, the constitutive relations can be expressed as(2)σrrσθθσrθDrDθ=Q11Q12Q13Q14Q15Q21Q22Q23Q24Q25Q31Q32Q33Q34Q35Q41Q42Q43Q44Q45Q51Q52Q53Q54Q55εrrεθθεrθErEθ,where εij and σij ij=xx,zz,xz are the strain and stress, respectively. Ei and Di are the electric field and electric displacement, respectively. The matrix elements Qiji,j=1,…,5 are the product of the elastic modulus (c^11, c^13, c^33, and c^44), the piezoelectric constant (e^15, e^31, and e^33) and the dielectric constant (k^11 and k^33), and the θ−β trigonometric function [12], in which(3)εrrεθθεrθErEθ=∂urr,θ∂rurr,θr+∂uθr,θ/r∂θ∂urr,θr∂θ+∂uθr,θ/∂r−uθr,θ/r−∂φr,θ∂r−∂φr,θr∂θ.The displacement fieldurr,θ and potential field φr,θ near the tip of the notch are expressed in the form of asymptotic expansion as follows [14]:(4)urr,θ=∑k=1nAkrλku˜rkθ,uθr,θ=∑k=1nAkrλku˜θkθ,φr,θ=∑k=1nAkrλkφ˜kθ,where r is the radial distance from the tip of notch, Ak is the combined amplitude coefficient, λk is the singularity order for the tip of notch, n is the number of series terms intercepted, and u˜ikθi=r,θ and φ˜kθ are the characteristic angular function of displacement and characteristic angular function of electric potential, respectively.In the absence of body force and free charge density, the equilibrium equations of the piezoelectric material are given by(5)∇rσ˜r−σθθ∇rσ˜θ+σrθ∇rD=0,where ∇=∂/∂r,1/r∂/∂θ, σr=σrrσrθ, σθ=σrθσθθ, and D=DrDθ.The typical termsAkrλku˜rkθ, Akrλku˜θkθ, and Akrλkφ˜kθ in equation (4) are substituted into equations (2) and (3) and then substituted into equilibrium equation (5) of the piezoelectric material (for brevity, u˜ikθi=r,θ, φ˜kθ, and λk are abbreviated as u˜ii=r,θ, φ˜, and λ, respectively), and the equations can be obtained as follows:(6)A˜u˜r″θu˜θ″θφ˜″θ+B˜u˜r′θu˜θ′θφ˜′θ+C˜u˜rθu˜θθφ˜θ+λE˜u˜r′θu˜θ′θφ˜′θ+λF˜u˜rθu˜θθφ˜θ+λ2G˜u˜rθu˜θθφ˜θ=000,where(7)A˜=A˜11θA˜12θA˜13θA˜21θA˜22θA˜23θA˜31θA˜32θA˜33θB˜=B˜11θB˜12θB˜13θB˜21θB˜22θB˜23θB˜31θB˜32θB˜33θC˜=C˜11θC˜12θC˜13θC˜21θC˜22θC˜23θC˜31θC˜32θC˜33θ,E˜=E˜11θE˜12θE˜13θE˜21θE˜22θE˜23θE˜31θE˜32θE˜33θF˜=F˜11θF˜12θF˜13θF˜21θF˜22θF˜23θF˜31θF˜32θF˜33θG˜=G˜11θG˜12θG˜13θG˜21θG˜22θG˜23θG˜31θG˜32θG˜33θ,where matrix elements A˜ijθ, B˜ijθ, C˜ijθ, E˜ijθ, F˜ijθ, and G˜ijθ are the linear combinations of the matrix element Qij and its first derivative Qij′ with respect to the coordinate θ. ⋯′ is the first derivative with respect to the coordinate θ, and ⋯″ is the second derivative with respect to the coordinate θ. The following is the same.


## 3. Singularity in Wedge Notch of Composite Materials

For the planar wedge notch of orthotropic composite materials as shown in Figure2, the material principal axes are identified as axes 1 and 3. γ is the angle between the principal axis system (1, 3) and the global coordinate system x,y of the orthotropic material, that is, the direction angle of the composite fiber. θ is the angle between the polar coordinate system r,θ and global coordinate system x,y, which is positive in the counterclockwise direction.Figure 2
A wedge notch of orthotropic composite materials.According to the basic theory of elastic mechanics, in the case of plane stress, the stress-strain relationships of orthotropic materials in the principal axis system (1, 3) are as follows:(8)σ11σ33σ13=C11C130C13C33000C66ε11ε33ε13,where Ciji,j=1,2,…,6 is the stiffness coefficient matrix of composite materials in the principal axis system (1, 3), in which(9)C11=E11−υ13υ31,C13=υ13E31−υ13υ31=υ31E11−υ13υ31,C33=E31−υ13υ31,C66=G13.According to the coordinate transformation, the constitutive equation of orthotropic materials in the polar coordinate systemorθ can be expressed as(10)ΠrrΠθθΠrθ=D11D12D16D12D22D26D16D26D66ΣrrΣθθΣrθ,where the matrix element Diji,j=1,…,6 is the product of material elastic moduli E1 and E3, shear modulus G13, Poisson’s ratios υ13 and υ31, and θ−γ trigonometric functions in the principal axis system (1, 3).The strain geometry equation of orthotropic materials in polar coordinate systemorθ can be expressed as(11)ΣrrΣθθΣrθ=∂ψrr,θ∂rψrr,θr+1r∂ψθr,θ∂θ1r∂ψrr,θ∂θ+∂ψθr,θ∂r−ψθr,θr,where ψrr,θ and ψθr,θ are the displacement fields at the tip of wedge notch and Σrr, Σθθ, and Σrθ are the stress field at the tip of the wedge notch.The origin of the polar coordinate system is set at the tip of the notch and the displacement field of elastic deformation near the tip of the notch can be asymptotically expanded to a series of forms about radial direction [14]:(12)ψrr,θ=∑k=1nAkrλkψ⌢rkθ,ψθr,θ=∑k=1nAkrλkψ⌢θkθ,where λk is the stress singularity order, ψ⌢rkθ and ψ⌢θkθ are the characteristic angular function of displacement, Ak is the corresponding displacement amplitude coefficient k=1,2,……n, and n is the number of series terms intercepted.The equilibrium equations of the orthotropic composite materials without considering the body force are(13)∇rΠr−Πθθ∇rΠθ+Πrθ=0,where ∇=∂/∂r,1/r∂/∂θ, Πr=ΠrrΠrθ, and Πθ=ΠrθΠθθ.The typical termsAkrλkψ⌢rkθ and Akrλkψ⌢θkθ in equation (12) are substituted into equations (10) and (11) and then substituted into the equilibrium equation (13). For brevity, ψ⌢ikθi=r,θ and λk are abbreviated as ψ⌢iθ and λ, respectively, and the equations can be obtained as follows:(14)A⌢ψ⌢r″θψ⌢θ″θ+B⌢ψ⌢r′θψ⌢θ′θ+C⌢ψ⌢rθψ⌢θθ+λE⌢ψ⌢r′θψ⌢θ′θ+λF⌢ψ⌢rθψ⌢θθ+λ2G⌢ψ⌢rθψ⌢θθ=00,where(15)A⌢=A⌢11θA⌢12θA⌢21θA⌢22θ,B⌢=B⌢11θB⌢12θB⌢21θB⌢22θ,C⌢=C⌢11θC⌢12θC⌢21θC⌢22θ,E⌢=E⌢11θE⌢12θE⌢21θE⌢22θ,F⌢=F⌢11θF⌢12θF⌢21θF⌢22θ,G⌢=G⌢11θG⌢12θG⌢21θG⌢22θ,where matrix elements A⌢ijθ, B⌢ijθ, C⌢ijθ, E⌢ijθ, F⌢ijθ, and G⌢ijθ are the linear combinations of the matrix element Dij and its first derivative Dij′ with respect to the coordinate θ.


## 4. Singularity in Piezoelectric Composite Junctions

As shown in Figure3, the electromechanical boundary conditions for the two free boundaries of junctions are free stress and free electric displacement, i.e.,Figure 3
Configuration of piezoelectric-composite wedge junctions.Whenθ=θ1, σθθ=σrθ=0, and Dθ=0,(16)R1u˜r′θu˜θ′θφ˜′θ+R2u˜rθu˜θθφ˜θ+λR3u˜rθu˜θθφ˜θ=000,where R1=Q23Q22−Q25Q33Q32−Q35Q53Q52−Q55, R2=Q22−Q230Q32−Q330Q52−Q530, and R3=Q21Q23−Q24Q31Q33−Q34Q51Q53−Q54.Whenθ=θ3, the stress is free, that is,(17)S1ψ⌢r′θψ⌢θ′θ+S2ψ⌢rθψ⌢θθ+λS3ψ⌢rθψ⌢θθ=00,where S1=D26D22D66D26, S2=D22−D26D26−D66, and S3=D12D26D16D66.Whenθ=θ2, the interface conditions between piezoelectric and composite materials are σθθ=Πθθ, σrθ=Πrθ, Dθ=0, ur=ψr, and uθ=ψθ, that is,(18)T1u˜r′θu˜θ′θφ˜′θψ⌢r′θψ⌢θ′θ+T2u˜rθu˜θθφ˜θψ⌢rθψ⌢θθ+λT3u˜rθu˜θθφ˜θψ⌢rθψ⌢θθ=00000,where(19)T1=Q23Q22−Q25−D26−D22Q33Q32−Q35−D66−D26Q53Q52−Q55000000000000,T2=Q22−Q230−D22D26Q32−Q330−D26D66Q52−Q53000100−100100−1,T3=Q21Q23−Q24−D12−D26Q31Q33−Q34−D16−D66Q51Q53−Q54000000000000.Thus, the problem of calculating the singularity orders of the junctions bonded by piezoelectric and composite wedges is transformed into the problem of solving the eigenvalues of ordinary differential equations (6) and (14) under boundary conditions equations (16) and (17) and interface conditions equation (18).


## 5. Numerical Method

As shown in Figure3, the dimensionless coordinate ξ is introduced into the wedge domain of piezoelectric material θ1≤θ≤θ2:(20)ξ=θ−θ1θ2−θ1,0≤ξ≤1.The characteristic angular functions of displacement, characteristic angular functions of electric potential, matrix elementsQij, and their derivatives of piezoelectric materials can be expressed as(21)u˜r″θu˜θ″θφ˜″θ=1θ2−θ12u˜r″ξu˜θ″ξφ˜″ξ,u˜r′θu˜θ′θφ˜′θ=1θ2−θ1u˜r′ξu˜θ′ξφ˜′ξ,u˜rθu˜θθφ˜θ=u˜rξu˜θξφ˜ξ,(22)Qij′θ=1θ2−θ1Qij′ξ,Qijθ=Qijξ.By substituting equations (21) and (22) into equation (6) for dimensionless coordinate, we can obtain(23)1θ2−θ12A˜u˜r″ξu˜θ″ξφ˜″ξ+1θ2−θ1B˜u˜r′ξu˜θ′ξφ˜′ξ+C˜u˜rξu˜θξφ˜ξ+λ1θ2−θ1E˜u˜r′ξu˜θ′ξφ˜′ξ+λF˜u˜rξu˜θθφ˜θ+λ2G˜u˜rθu˜θθφ˜θ=000.Then, the analysis of the electroelastic singularity field at the apex of the wedge of piezoelectric materials is transformed into the eigenvalue problem of ordinary differential equations with dimensionless coordinateξ.The dimensionless coordinateξ is introduced into the wedge domain of orthotropic composite material θ2≤θ≤θ3:(24)ξ=θ−θ2θ3−θ2,0≤ξ≤1.The characteristic angular functions of displacement, matrix elementsDij, and their derivatives of composite materials can be expressed as(25)ψ⌢r″θψ⌢θ″θ=1θ3−θ22ψ⌢r″ξψ⌢θ″ξ,ψ⌢r′θψ⌢θ′θ=1θ3−θ2ψ⌢r′ξψ⌢θ′ξ,ψ⌢rθψ⌢θθ=ψ⌢rξψ⌢θξ,(26)Dij′θ=1θ3−θ2Dij′ξ,Dijθ=Dijξ.By substituting equations (25) and (26) into equation (14) for dimensionless coordinate, we can obtain(27)1θ3−θ22A⌢ψ⌢r″ξψ⌢θ″ξ+1θ3−θ2B⌢ψ⌢r′ξψ⌢′θξ+C⌢ψ⌢rξψ⌢θξ+λ1θ3−θ2E⌢ψ⌢r′ξψ⌢θ′ξ+λF⌢ψ⌢rξψ⌢θξ+λ2G⌢ψ⌢rξψ⌢θξ=00.Then, the analysis of the singular stress field in wedge notch of orthotropic composite materials is transformed into the eigenvalue problem of ordinary differential equations with dimensionless coordinateξ.By introducing dimensionless coordinateξ, the boundary and interface conditions are dimensionless and expressed as(28)1θ2−θ1R1u˜r′ξu˜θ′ξφ˜′ξ+R2u˜rξu˜θξφ˜ξ+λR3u˜rξu˜θξφ˜ξ=000,(29)1θ3−θ2S1ψ⌢r′ξψ⌢θ′ξ+S2ψ⌢rξψ⌢θξ+λS3ψ⌢rξψ⌢θξ=00,(30)T11θ2−θ1u˜r′ξ1θ2−θ1u˜θ′ξ1θ2−θ1φ˜′ξ1θ3−θ2ψ⌢r′ξ1θ3−θ2ψ⌢θ′ξ+T2u˜rξu˜θξφ˜ξψ⌢rξψ⌢θξ+λT3u˜rξu˜θξφ˜ξψ⌢rξψ⌢θξ=00000.The wedge domainsθ1≤θ≤θ2 and θ2≤θ≤θ3 are transformed into dimensionless region 0≤ξ≤1, and the interval 0,1 is divided into N segments, 0=ξ0,ξ1,ξ2,…,ξn−1,ξn=1, and Δli=ξi−ξi−1=1/N. By using the difference method, the derivative values of the five functions u˜rξ, u˜θξ, φ˜ξ, ψ⌢rξ, and ψ⌢θξ in the system of ordinary differential equations with variable coefficients are expressed by the function values on the interval partition points:(31)u˜r′ξj−u˜r′ξ0=∫ξ0ξju˜r″ξdξ,u˜θ′ξj−u˜θ′ξ0=∫ξ0ξju˜θ″ξdξ,φ˜′ξj−φ˜′ξ0=∫ξ0ξjφ˜″ξdξ,ψ⌢r′ξj−ψ⌢r′ξ0=∫ξ0ξjψ⌢r″ξdξ,ψ⌢θ′ξj−ψ⌢θ′ξ0=∫ξ0ξjψ⌢θ″ξdξ,j=0,1,2,…,n.In the above formula,u˜r″ξ, u˜θ″ξ, φ˜″ξ, ψ⌢r″ξ, and ψ⌢θ″ξ are approximated by interpolation functions without loss of generality. Take function u˜r″ξ as an example to illustrate(32)u˜r″ξ=∑i=1Nu˜r″ξiLiξ,where Liξ is the Lagrange interpolation basis function, so(33)u˜r′ξj−u˜r′ξ0=∑i=0nu˜r″ξjJji,i=0,1,2,…,N;j=0,1,2,…,N,where Jji=∫ξ0ξjLiξdξ, and vectors and matrix notations are introduced(34)τ=00⋮0,σ=11⋮1,J=000⋯0J10J11J12⋯J1nJ20J21J22⋯J2n⋯⋯⋯⋯⋯Jn0Jn1Jn2⋯Jnn,I=100⋯0010⋯0001⋯0⋯⋯⋯⋯⋯000⋯1,where the matrix J is called as the integral matrix, which only depends on the interpolation basis function Liξ. Quadratic parabola interpolation is adopted, and the basis function is expressed as(35)Liξ=ξ−ξiξ−ξi+12Δli2ξ−ξi−1ξ−ξi+1−Δli2ξ−ξi−1ξ−ξi2Δli2,ξ∈ξi−1,ξi+1.The functionu˜rξ and its derivative functions of each order in equation (31) are written in the vector form as follows:(36)u˜r′ξi=τu˜rξ0+σu˜r′ξ0+Ju˜r″ξi=τ,σ,JN+1×N+3u˜rξ0u˜r′ξ0u˜r″ξiN+3×1.Let(37)H1=τ,σ,Jn+1×n+3,Y˜r=u˜rξ0u˜r″ξ0u˜r″ξiN+3×1,such that equation (36) can be expressed as(38)u˜r′ξi=H1Y˜r.Then, the low-order derivative function is replaced by the higher-order derivative function in sequence, and the following formula can be obtained by stepwise recursion:(39)u˜rξi=σu˜rξ0+Jσu˜r′ξ0+J2u˜r″ξi=σ,Jσ,J2N+1×N+3u˜rξ0u˜r′ξ0u˜r″ξiN+3×1.Let(40)H0=σ,Jσ,J2N+1×N+3.such that equation (39) can be expressed as(41)u˜rξi=H0Y˜r,u˜r″ξi=τu˜rξ0+τu˜r′ξ0+Iu˜r″ξi=τ,τ,IN+1×N+3u˜rξ0u˜r′ξ0u˜r″ξiN+3×1.Equation (41) can be expressed as(42)H2=τ,τ,IN+1×N+3,u˜r″ξi=H2Y˜r.The other function vectors are handled in the same way and given by(43)u˜θξi=H0Y˜θ,u˜θ′ξi=H1Y˜θ,u˜θ″ξi=H2Y˜θ,φ˜ξi=H0Y˜φ,φ˜′ξi=H1Y˜φ,φ˜″ξi=H2Y˜φ,ψ⌢rξi=H0Y⌢r,ψ⌢r′ξi=H1Y⌢r,ψ⌢r″ξi=H2Y⌢r.ψ⌢θξi=H0Y⌢θ,ψ⌢θ′ξi=H1Y⌢θ,ψ⌢θ″ξi=H2Y⌢θ,where(44)Y˜θ=u˜θξ0u˜θ′ξ0u˜θ″ξiN+3×1,Y˜φ=φ˜ξ0φ˜′ξ0φ˜″ξiN+3×1,Y⌢r=ψ⌢rξ0ψ⌢r′ξ0ψ⌢r′ξiN+3×1,Y⌢θ=ψ⌢θξ0ψ⌢θ′ξ0ψ⌢θ″ξiN+3×1.Each coefficient in equations (23), (27), and (28)–(30) is written as a diagonal matrix of N+1×N+1. For example, A˜ijξ and A⌢ijξ are transformed into the diagonal matrix after discretization, and others are similar:(45)A˜ij=A˜ijξ000⋯00A˜ijξ10⋯000A˜ijξ3⋯0⋯⋯⋯⋯⋯000⋯A˜ijξN,A⌢ij=A⌢ijξ000⋯00A⌢ijξ10⋯000A⌢ijξ3⋯0⋯⋯⋯⋯⋯000⋯A⌢ijξN.By substituting equations (43)–(45) into equations (23), (27), and (28)–(30), the calculation of the singularity order of the electromechanical coupling field at the tip of piezoelectric composite material junctions is transformed into a set of algebraic equation eigenvalue problems:(46)LY+λMY=0,where Y=Y˜rY˜θY˜φY⌢rY⌢θT, the matrices L and M are the linear combinations of the diagonal matrices A˜ij, B˜ij, C˜ij, … and A⌢ij, B⌢ij, C⌢ij, …. Equation (46) is a general characteristic equation set of 5n+3 unknown vectors Y=Y˜rY˜θY˜φY⌢rY⌢θT. By using the QR method to solve the equation, many eigenvalues λ and corresponding eigenvectors Y=Y˜rY˜θY˜φY⌢rY⌢θT can be obtained.


## 6. Example Analysis and Discussion

The piezoelectric composite material wedge junction is shown in Figure4. The junction angle of PZT-4/PZT-5 piezoelectric material wedge notch is ϕ2=180∘. The junction angle of composite material wedge notch is within the range 0∘<ϕ1≤180∘. When ϕ1=180∘, the bimaterial junction structure degenerates into interface crack structure. The first four-order values of the stress and electric displacement singular characteristic order in the interface crack tip region are shown in Table 1 by using the interpolation matrix method (IMM). In Table 1, N refers to the number of segments at the tip of each wedge notch. The unknown function in the subinterval is simulated by quadratic piecewise polynomial function. The more segments N is, the higher the calculation accuracy is, and Reλ1−1 converges to -0.500000. Here, −1<Reλ1−1<0 is the singular term and Reλ2−1>0, Reλ3−1>0, and Reλ4−1>0 are the nonsingular terms. These nonsingular parameters are indispensable to describe the complete stress field at the tip of piezoelectric composite material junctions [24]. At present, there is still a lack of these singular order values for the tip of piezoelectric composite material junctions, especially for the region far from the tip of junctions. When ϕ1=90∘, the first four-order values of the stress singular characteristic order in the bimaterial junctions are shown in Table 2. The results obtained in this paper are in good agreement with the calculated results in the existing literature [15].Figure 4
Geometry for piezoelectric-composite wedge junctions.Table 1
The first four singularity orders of the interface crack with interval segmentation numberN.
N
λ
1
−
1
λ
2
−
1
λ
3
−
1
λ
4
−
1
10
−0.510769–0.0406356i
−0.510769 + 0.0406356i
0.504289–0.0455896i
0.504289 + 0.0455896i
20
−0.502625–0.0448070i
−0.502625 + 0.0448070i
0.500793–0.0440447i
0.500793 + 0.0440447i
30
−0.500487–0.0444804i
−0.500487 + 0.044480a4i
0.500149–0.0440330i
0.500149 + 0.0440330i
40
−0.500063–0.0442202i
−0.500063 + 0.0442202i
0.500035–0.0440385i
0.500035 + 0.0440385i
50
−0.499986–0.0441066i
−0.499986 + 0.0441066i
0.500010–0.0440370i
0.500010 + 0.0440370i
60
−0.499978–0.0440604i
−0.499978 + 0.0440604i
0.500003–0.0440345i
0.500003 + 0.0440345i
70
−0.499982–0.0440414i
−0.499982 + 0.0440414i
0.500001–0.0440325i
0.500001 + 0.0440325i
80
−0.499987–0.0440334i
−0.499987 + 0.0440334i
0.500000–0.0440311i
0.500000 + 0.0440311i
Ref. [15]
−0.5–0.044i
−0.5 + 0.044i
—
—
Note:i=−1 is the imaginary number.Table 2
The first four singularity orders of piezoelectric composite material wedge junction with interval segmentation numberN, when ϕ1=90∘.
N
λ
1
−
1
λ
2
−
1
λ
3
−
1
λ
4
−
1
10
−0.34958611
−0.022406747
0.62279008–0.38780841i
0.6227900 + 0.387808418i
20
−0.34989329
−0.022582134
0.62231249–0.387750664i
0.62231249 + 0.387750664i
40
−0.34997618
−0.022601267
0.62225081–0.38760694i
0.62225081 + 0.38760694i
80
−0.34998674
−0.022600726
0.62224740–0.38759404i
0.62224740 + 0.38759404i
Ref. [15]
−0.34998
−0.0226004
—
—
Note:i=−1 is the imaginary number.As shown in Figure4, the junction angle is ϕ2=180∘ and the polarization angle of piezoelectric material PZT-4/PZT-5 is β=0∘. When the junction angle of composite wedge notch is 0∘<ϕ1≤180∘ and the number of segments N=40, the order values are calculated by the interpolation matrix method. Figure 5 shows the variations of singularity orders λk−1 with junction angle ϕ1 for the electromechanical coupling field of the junctions between piezoelectric material PZT-4/PZT-5 and composite material Gr/Ep. It can be seen that the calculated values are consistent with the solution in the literature [15], indicating the validity and accuracy of the calculated method in this paper. At the same time, Figure 5 also shows that with the increase of the junction angle ϕ1, the first two singularity orders λk−1 of the coupling field are negative and decrease gradually, that is, the singularity of the electromechanical coupling field becomes stronger. When ϕ1=180∘, the piezoelectric composite material junctions shown in Figure 4 degenerate into the interface crack model, the singularity of the electromechanical coupling field at the tip of junctions is the strongest, and the real part of the singularity order is Reλ1−1=−0.5000638.Figure 5
The variations of singularity orders withϕ1 for piezoelectric-composite wedge junctions.The method in this paper can not only obtain the useful singularity orders of stress and electric displacement at the tip of piezoelectric and composite material junctions but can also obtain the corresponding characteristic angular functions of displacement and electric potential and their derivative functions. Figure6 shows the curves of characteristic angular function corresponding to the first stress singularity order λ1−1=−0.2828998 of PZT-4-Gr/Ep junction at ϕ1=60∘, ϕ2=180∘, and N=40. Figure 7 shows the curves of the first derivative function of the characteristic angular function corresponding to the first stress singularity order λ1−1=−0.2828998 of PZT-4—Gr/Ep junction at ϕ1=60∘, ϕ2=180∘, and N=40. It can be seen from Figures 6 and 7 that the characteristic angular functions of displacement and electric potential are continuous at the bonding position of dissimilar materials, and their first derivative functions change suddenly. These indicate that the displacement field and potential field in the bonding position of piezoelectric composite material junctions are continuous, but the stress and electric displacement are not equal, which is one of the reasons for the failure and damage of piezoelectric and composite material junctions.Figure 6
The characteristic angular function of displacement and electric potential.Figure 7
The first derivative characteristic angular function of displacement and electric potential.Figure8(a) shows the geometry for piezoelectric composite material junctions at ϕ1=90∘, ϕ2=90∘, and β=0∘ and fiber direction 0∘≤γ≤180∘. The singularity orders of electromechanical coupling field at the tip of PZT-4 and Gr/Ep junctions are calculated by the interpolation matrix method at N=40. The variation of singularity orders λ1−1 with the fiber direction γ is plotted in Figure 9. It is found that there is only first stress and electric displacement singularity order at the tip of bimaterial junctions at ϕ1=ϕ2=90∘ and 0∘≤γ≤180∘. Besides, the singularity at the tip of junctions is the strongest at γ=135∘, in which λ1−1=−0.0964289.Figure 8
Geometry for piezoelectric-composite wedge junctions: (a)ϕ1=90∘, ϕ2=90∘, and β=0∘ and (b) ϕ1=90∘ and ϕ2=180∘.
(a)
(b)Figure 9
The variations of singularity orders with the fiber directionγ of graphite/epoxy at ϕ1=ϕ2=90∘.Figure8(b) shows the geometry for piezoelectric composite material junctions at ϕ1=90∘ and ϕ2=180∘ and fiber direction 0∘≤γ≤180∘. The variation of singularity orders λ−1 with the fiber direction γ is plotted in Figure 10. Figure 10 illustrates that there are second stress and electric displacement singularity order at the tip of bimaterial junctions at ϕ1=90∘ and ϕ2=180∘ and 0∘≤γ≤180∘. The fiber direction γ has a great influence on the first singularity order λ1−1, but a little influence on the second singularity order λ2−1. Besides, the singularity at the tip of junctions is the strongest at γ=45∘, where λ1−1=−0.4575995.Figure 10
The variations of singularity orders with the fiber directionγ of graphite/epoxy at ϕ1=90∘ and ϕ2=180∘.Figures11 and 12 show the characteristic angular function of displacement and electric potential and their first derivative function at ϕ1=ϕ2=90∘ and γ=135∘, where the first stress singularity order at the tip of PZT−4 and Gr/Ep junctions λ1−1=−0.0964289. Figures 13 and 14 show the characteristic angular function of displacement and electric potential and their first derivative function at ϕ1=90∘, ϕ2=180∘, and γ=45∘, where the first stress singularity order at the tip of PZT−4 and Gr/Ep junctions λ1−1=−0.4575995.Figure 11
The characteristic angular function of displacement and electric potential atλ1−1=−0.0964289.Figure 12
The first derivative characteristic angular function of displacement and electric potential atλ1−1=−0.0964289.Figure 13
The characteristic angular function of displacement and electric potential atλ1−1=−0.4575995.Figure 14
The first derivative characteristic angular function of displacement and electric potential atλ1−1=−0.4575995.


## 7. Conclusions

The numerical method based on the interpolation matrix method was presented to obtain the first couple of singularity orders of electromechanical coupling field at the tip of piezoelectric composite material junctions. The corresponding characteristic angular functions of displacement and electric potential and their derivative functions of each order were obtained at the same time. All these calculated values have the same accuracy, which is a significant advantage. Moreover, the desired calculating accuracy can be achieved by subdividing the interval mesh. The main conclusions are summarized as follows:(1)
The singularity orders at the tip of piezoelectric composite material junctions obtained in this paper are in good agreement with the calculated results in the existing literature. The effectiveness and high accuracy of the presented method are proved in calculating the singularity orders of the stress field and electric field at the tip of piezoelectric composite material junctions.(2)
The singularity of the electromechanical coupling field at the tip of piezoelectric composite material junctions is closely related to the bonding angle and fiber direction. The numerical results can be used to guide the optimal configuration for the combination of dissimilar materials.(3)
Based on the results of this paper, the complete stress field and electric displacement field can be obtained by using the finite element or boundary element method for the region around the tip of piezoelectric composite material junctions.



---

*Source: 5616138-2020-02-25.xml*