A coupled electromechanical analysis of a piezoelectric layer bonded to an elastic substrate: Part I,development of governing equations

This two-part contribution presents a novel and efficient method to analyze the two-dimensional (2-D) electromechanical fields of a piezoelectric layer bonded to an elastic substrate, which takes into account the fully coupled electromechanical behavior. In Part I, Hellinger–Reissner variational principle for elasticity is extended to electromechanical problems of the bimaterial, and is utilized to obtain the governing equations for the problems concerned. The 2-D electromechanical field quantities in the piezoelectric layer are expanded in the thickness-coordinate with seven one-dimensional (1-D) unknown functions. Such an expansion satisfies exactly the mechanical equilibrium equations, Gauss law, the constitutive equations, two of the three displacement–strain relations as well as one of the two electric field-electric potential relations. For the substrate the fundamental solutions of a half-plane subjected to a vertical or horizontal concentrated force on the surface are used. Two differential equations and two singular integro-differential equations of four unknown functions, the axial force, N, the moment, M, the average and the first moment of electric displacement, D₀ and D₁, as well as the associated boundary conditions have been derived rigorously from the stationary conditions of Hellinger–Reissner variational functional. In contrast to the thin film/substrate theory that ignores the interfacial normal stress the present one can predict both the interfacial shear and normal stresses, the latter one is believed to control the delamination initiation.     

Subinterface crack in an anisotropic piezoelectric bimaterial

In this paper, the problem of a subinterface crack in an anisotropic piezoelectric bimaterial is analyzed. A system of singular integral equations is formulated for general anisotropic piezoelectric bimaterial with kernel functions expressed in complex form. For commonly used transversely isotropic piezoelectric materials, the kernel functions are given in real forms. By considering special properties of one of the bimaterial, various real kernel functions for half-plane problems with mechanical traction-free or displacement-fixed boundary conditions combined with different electric boundary conditions are obtained. Investigations of half-plane piezoelectric solids show that, particularly for the mechanical traction-free problem, the evaluations of the mechanical stress intensity factors (electric displacement intensity factor) under mechanical loadings (electric displacement loading) for coupled mechanical and electric problems may be evaluated directly by considering the corresponding decoupled elastic (electric) problem irrespective of what electric boundary condition is applied on the boundary. However, for the piezoelectric bimaterial problem, purely elastic bimaterial analysis or purely electric bimaterial analysis is inadequate for the determination of the generalized stress intensity factors. Instead, both elastic and electric properties of the bimaterial’s constants should be simultaneously taken into account for better accuracy of the generalized stress intensity factors.     

Electro-elastic green's functions for a piezoelectric half-space and their application

I.IntroductionItiswell'knownthatoneofthemostpowerfultoolsinlinearfieldtheoriesistheGreen'sfunction.Fore1asticity,considerableresearchcanbefoundintheliterature.However,theGreen'sfunctionforpiezoe1ectricityisratherlimitedduetotheanisotropyandelectromechanic…     

Electro-elastic interaction between a piezoelectric screw dislocation and circular interfacial rigid lines

The electro-elastic interaction between a piezoelectric screw dislocation located either outside or inside inhomogeneity and circular interfacial rigid lines under anti-plane mechanical and in-plane electrical loads in linear piezoelectric materials is dealt with in the framework of linear elastic theory. Using Riemann–Schwarz’s symmetry principle integrated with the analysis of singularity of complex functions, the general solution of this problem is presented in this paper. For a special example, the closed form solutions for electro-elastic fields in matrix and inhomogeneity regions are derived explicitly when interface containing single rigid line. Applying perturbation technique, perturbation stress and electric displacement fields are obtained. The image force acting on piezoelectric screw dislocation is calculated by using the generalized Peach–Koehler formula. As a result, numerical analysis and discussion show that soft inhomogeneity can repel screw dislocation in piezoelectric material due to their intrinsic electro-mechanical coupling behavior and the influence of interfacial rigid line upon the image force is profound. When the radian of circular rigid line reaches extensive magnitude, the presence of interfacial rigid line can change the interaction mechanism.     

Interaction between a screw dislocation and a viscoelastic piezoelectric bimaterial interface

The complex variable method is employed to derive analytical solutions for the interaction between a piezoelectric screw dislocation and a Kelvin-type viscoelastic piezoelectric bimaterial interface. Through analytical continuation, the original boundary value problem can be reduced to an inhomogeneous first-order partial differential equation for a single function of location z = x + iy and time t defined in the lower half-plane, which is free of the screw dislocation. Once the initial, steady-state and far-field conditions are known, the solution to the first order differential equation can be obtained. From the solved function, explicit expressions are then derived for the stresses, strains, electric fields and electric displacements induced by the piezoelectric screw dislocation. Also presented is the image force acting on the screw dislocation due to its interaction with the Kelvin-type viscoelastic interface. The derived solutions are verified by comparing with existing solutions for the simplified cases, and various interesting features are observed, particularly for those associated with the image force.     

Electro-mechanical frictionless contact behavior of a functionally graded piezoelectric layered half-plane under a rigid punch

The frictionless contact problem of a functionally graded piezoelectric layered half-plane in-plane strain state under the action of a rigid flat or cylindrical punch is investigated in this paper. It is assumed that the punch is a perfect electrical conductor with a constant potential. The electro-elastic properties of the functionally graded piezoelectric materials (FGPMs) vary exponentially along the thickness direction. The problem is reduced to a pair of coupled Cauchy singular integral equations by using the Fourier integral transform technique and then is numerically solved to determine the contact pressure, surface electric charge distribution, normal stress and electric displacement fields. For a flat punch, the normal stress intensity factor and electric displacement intensity factor are also given to quantitatively characterize the singularity behavior at the punch ends. Numerical results show that both material property gradient of the FGPM layer and punch geometry have a significant influence on the contact performance of the FGPM layered half-plane.     

PLANE STRAIN PROBLEM OF PIEZOELECTRIC SOLID WITH ELLIPTIC INCLUSION

By using the complex variables function theory, a plane strain electro-elastic analysis was performed on a transversely isotropic piezoelectric material containing an elliptic elastic inclusion, which is subjected to a uniform stress field and a uniform electric displacement loads at infinity. Based on the present finite element results and some related theoretical solutions, an acceptable conjecture was found that the stress field is constant inside the elastic inclusion. The stress field solutions in the piezoelectric matrix and the elastic inclusion were obtained in the form of complex potentials based on the impermeable electric boundary conditions.     

The interaction between a screw dislocation and a piezoelectric fiber composite with a wedge crack

The interaction problem between a screw dislocation and a piezoelectric fiber composite with a semi-infinite wedge crack is investigated in this paper. The piezoelectric media are assumed to be transversely isotropic with the poling direction along the x ₃ direction. The screw dislocation considered here involves a Burgers vector parallel to the poling direction with a line force and a line charge being applied at the core of the dislocation. Both cases of the screw dislocation located at the matrix and inclusion are observed. The analytical derivation is based on the complex variable and the conformal mapping methods. The exact solutions are obtained to calculate the forces on the dislocation and the crack-tip stress and electric displacement intensity factors. Based on these results, the anti-shielding and shielding effects for different loadings, material combinations, and geometric configurations are discussed in detail.     

Solutions for transversely isotropic piezoelectric infinite body,semi-infinite body and bimaterial infinite body subjected to uniform ring loading and charge

Based on the fundamental solutions for transversely isotropic piezoelectric materials, the fundamental solutions of axisymmetric problems are derived by integration and explicit expressions for three possible cases of different characteristic roots and multiple roots are all presented. In the case of s₁ ≠ s₂ ≠ s₃ ≠ s₁, based on the Greens functions for semi-infinite piezoelectric body and bimaterial infinite piezoelectric body, the Greens functions for axisymmetric problems of semi-infinite body and bimaterial infinite body are obtained. Taking PZT-4 as an example, numerical computations are conducted by use of the fundamental solutions to axisymmetric problems. Comparison of the calculated results with those of FEM shows good agreement between them.     

Electro-elastic fundamental solutions of anisotropic piezoelectric materials with an elliptical hole

By using Stroh' complex formalism and Cauchy's integral method, the electro-elastic fundamental solutions of an infinite anisotropic piezoelectric solid containing an elliptic hole or a crack subjected to a line force and a line charge are presented in closed form. Particular attention is paid to analyzing the characteristics of the stress and electric displacement intensity factors. When a line force-charge acts on the crack surface, the real form expression of intensity factors is obtained. It is shown through a special example that the present work is correct. The project supported by the Fund of the State Education Commission of China for Excellent Young Teachers     

A conducting arc crack between a circular piezoelectric inclusion and an unbounded matrix

The present paper investigates the problem of a conducting arc crack between a circular piezoelectric inclusion and an unbounded piezoelectric matrix. The original boundary value problem is reduced to a standard Riemann–Hilbert problem of vector form by means of analytical continuation. Explicit solutions for the stress singularities δ=−(1/2)±iε are obtained, closed form solutions for the field potentials are then derived through adopting a decoupling procedure. In addition, explicit expressions for the field component distributions in the whole field and along the circular interface are also obtained. Different from the interface insulating crack, stresses, strains, electric displacements and electric fields at the crack tips all exhibit oscillatory singularities. We also define a complex electro-elastic field concentration vector to characterize the singular fields near the crack tips and derive a simple expression for the energy release rate, which is always positive, in terms of the field concentration vector. The condition for the disappearance of the index ε is also discussed. When the index ε is zero, we obtain conventionally defined electro-elastic intensity factors. The examples demonstrate the physical behavior and the correctness of the obtained solution.     

A mode III crack crossing the magnetoelectroelastic bimaterial interface under concentrated magnetoelectromechanical loads

A mode III crack cutting perpendicularly across the interface between two dissimilar semi-infinite magnetoelectroelastic solid is studied under the combined loads of a line force, a line electric charge and a line magnetic charge at an arbitrary location. The impermeable conditions are implied on the crack faces. The technique developed in literature for the elastic bimaterial with a crack cutting interface is exploited to treat the magnetoelectroelastic bimaterial. The Riemann-Hilbert problem can be formulated and solved based on complex variable method. Analytical solutions can be obtained for the entire plane. The intensity factors around crack tips can be defined for the elastic, electric and magnetic fields. It shows that, no matter where the load position is, the electric displacement intensity factors (EDIFs), as well as the magnetic induction intensity factors (MIIFs), are identical in magnitude but opposite in sign for both crack tips, on condition that a line force is solely applied. Alternatively, if only a line electric charge is considered, then the stress intensity factors (SIFs) and the MIIFs exhibit the behavior. Likewise, if only a line magnetic charge is applied, it turns to the SIFs and the EDIFs instead. In addition, the dependence of the intensity factors is graphically shown with respect to the location of a line force. It is found that the SIF for a crack tip tends to be infinite if the applied force is approaching the tip itself, but the EDIF, with the complete opposite trend, tends to be vanishing. Finally, focusing on the more practical case of piezoelectric/piezomagnetic bimaterial, variation of the SIF along with the moduli as well as the piezo constitutive coefficients is explored. These analyses may provide some guidance for material selection by minimizing the SIF. It is also believed that the results obtained in this paper can serve as the Green’s function for the dissimilar magnetoelectroelastic semi-infinite bimaterial with a crack cutting the interface under general magnetoelectromechanical loads.     

A piezoelectric screw dislocation near a wedge-shaped bi-material interface

The electro-elastic stress field due to a piezoelectric screw dislocation near the tip of a wedge-shaped bi-material interface is derived. The screw dislocation is subjected to a line charge and a line force at the core. The explicit closed-form analytical solutions for the stress field are derived by means of the complex variable and conformal mapping methods. The stress and electric intensity factors of the wedge tip induced by the dislocation and the image force acting on the dislocation are also formulated and calculated. The influence of the wedge angle and the different bi-material constant combinations on the image force is discussed. Numerical results for three particular wedge angles are calculated and compared.     

Elastic fields of dislocation loops in three-dimensional anisotropic bimaterials

By applying semi-analytical point-force Green's functions obtained via the Stroh formulism, we derive simple line integrals to calculate the elastic displacement and stress fields for a three-dimensional dislocation loop in an anisotropic bimaterial system. The solutions for the case of anisotropy are more convenient for treating an arbitrary dislocation loop compared with traditional area integration. With this new formulation, we numerically examine the displacement, stress, and energy due to the interaction between a dislocation loop and the bimaterial interface in an Al–Cu system. The interactive image energy due to the elastic moduli mismatch across the interface is then numerically evaluated. The result shows that a dislocation loop is subjected to an attractive force by the interface when it lies in the stiff material, and a repulsive force when it lies in the soft material. Moreover, the dependence of the interactive image energy of a dislocation loop on the position and size of the dislocation loop are also demonstrated and discussed. Significantly, it is found that the interactive image energy for a dislocation loop depends only on the ratio d/a, where a is the loop diameter and d is its distance to the interface. The examples studied provide benchmark solutions for anisotropic bimaterial dislocation problems.     

Numerical analysis for a crack in piezoelectric material under impact

In this paper, a numerical analysis of impact interfacial fracture for a piezoelectric bimaterial is provided. Starting from the basic equilibrium equation, a dynamic electro-mechanical FEM formulation is briefly presented. Then, the path-independent separated dynamic J integral is extended to piezoelectric bimaterials. Based on the relationship of the path-independent dynamic J integral and the stress and electric displacement intensity factors, the component separation method is used to calculate the stress and electric displacement intensity factors for piezoelectric bimaterials in this finite-element analysis. The response curves of the dynamic J integral, the stress and electric displacement intensity factors are obtained for both homogeneous material (PZT-4 and CdSe) and CdSe/PZT-4 bimaterial. The influences of the piezoelectricity and the electro-mechanical coupling factor on these responses are discussed. The effects of an applied electric field are also discussed.     

Magneto-electro-elastic antiplane analysis of a bimaterial BaTiO3–CoFe2O4 composite wedge with an interface crack

The antiplane analysis is made for a bimaterial BaTiO₃–CoFe₂O₄ composite wedge containing an interface crack. The coupled magneto-electro-elastic field is induced by the piezoelectric/piezomagnetic BaTiO₃–CoFe₂O₄ composite materials. For the crack problems, the intensity factors of stress, strain, electric displacement, electric field, magnetic induction and magnetic field at crack tips are derived analytically. Also, the energy density criterion is applied to predict the fracture behavior of the interface crack. The numerical results also show that the energy release rate for a crack in a single wedge is negative.     

Fracture analysis of circular-arc interface cracks in piezoelectric materials

The anti-plane problem of N arc-shaped interfacial cracks between a circular piezoelectric inhomogeneity and an infinite piezoelectric matrix is investigated by means of the complex variable method. Cracks are assumed to be permeable and then explicit expressions are presented, respectively, for the electric field on the crack faces, the complex potentials in media and the intensity factors near the crack-tips. As examples, the corresponding solutions are obtained for a piezoelectric bimaterial system with one or two permeable arc-shaped interfacial cracks, respectively. Additionally, the solutions for the cases of impermeable cracks also are given by treating an impermeable crack as a particular case of a permeable crack. It is shown that for the case of permeable interfacial cracks, the electric field is jumpy ahead of the crack tips, and its intensity factor is always dependent on that of stress. Moreover all the field singularities are dependent not only on the applied mechanical load, but also on the applied electric load. However, for the case of a homogeneous material with permeable cracks, all the singular factors are related only to the applied stresses and material constants.     

Green’s functions of two-dimensional piezoelectric quasicrystal in half-space and bimaterials

In this paper, we obtain Green’s functions of two-dimensional (2D) piezoelectric quasicrystal (PQC) in half-space and bimaterials. Based on the elastic theory of QCs, the Stroh formalism is used to derive the general solutions of displacements and stresses. Then, we obtain the analytical solutions of half-space and bimaterial Green’s functions. Besides, the interfacial Green’s function for bimaterials is also obtained in the analytical form. Before numerical studies, a comparative study is carried out to validate the present solutions. Typical numerical examples are performed to investigate the effects of multi-physics loadings such as the line force, the line dislocation, the line charge, and the phason line force. As a result, the coupling effect among the phonon field, the phason field, and the electric field is prominent, and the butterfly-shaped contours are characteristic in 2D PQCs. In addition, the changes of material parameters cause variations in physical quantities to a certain degree.     

THREE-DIMENSIONAL INTERACTIONS OF A HALF-PLANE CRACK IN A TRANSVERSELY ISOTROPIC PIEZOELECTRIC SPACE WITH RESULTANT SOURCES

I. INTRODUCTIONBecause of the widespread application and intrinsic brittleness of piezoelectric ceramics, significantattention is being paid to the crack problems of piezoelectric ceramics. The last decade has seen a lotof three-dimensional studies of crack in piezoelectric ceramics[1??9]. In addition, the electroelastic fieldin a transversely isotropic piezoelectric space with a half-plane crack subjected to symmetric normalpoint forces, antisymmetric tangential point forces and point ch…     

Elastic field due to an edge dislocation in a multilayered composite

A numerical procedure is presented for the analysis of the elastic field due to an edge dislocation in a multilayered composite. The multilayered composite consists of n perfectly bonded layers having different material properties and thickness, and two half-planes adhere to the top and bottom layers. The stiffness matrices for each layer and the half-planes are first derived in the Fourier transform domain, then a set of global stiffness equations is assembled to solve for the transformed components of the elastic field. Since the singular part of the elastic field corresponding to the dislocation in the full-plane has been extracted from the transformed components, regular numerical integration is needed only to evaluate the inverse Fourier transform. Numerical results for the elastic field due to an edge dislocation in a bimaterial medium are shown in fairly good agreement with analytical solutions. The elastic field and the Peach–Kohler image force are also presented for an edge dislocation in a single layered half-plane, a two-layered half-plane and a multilayered composite made of alternating layers of two different materials.