Wave propagation in magneto-electro-elastic multilayered plates

An analytical treatment is presented for the propagation of harmonic waves in magneto-electro-elastic multilayered plates, where the general anisotropic and three-phase coupled constitutive equations are used. The state-vector approach is employed to derive the propagator matrix which connects the field variables at the upper interface to those at the lower interface of each layer. The global propagator matrix is obtained by propagating the solution in each layer from the bottom of the layered plate to the top using the continuity conditions of the field variables across the interfaces. From the global propagator matrix, we finally obtain the dispersion relation by imposing the traction-free boundary condition on the top and bottom surfaces of the layered plate. Dispersion curves, modal shapes, and natural frequencies are presented for layered plates made of orthotropic elastic (graphite–epoxy), transversely isotropic PZT-5A, piezoelectric BaTiO₃ and magnetostrictive CoFe₂O₄ materials. While the numerical results show clearly the influence of different stacking sequences as well as material properties on the field response, the general methodology presented in the paper could be useful to the analysis and design of layered composites made of smart piezoelectric and piezomagnetic materials.     

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.     

Magnetic and electric poling effects associated with crack growth in BaTiO3–CoFe2O4 composite

Magnetoelectroelastic composite possesses the dual feature that the application of magnetic field induces electric polarization and electric field induces magnetization. The poling directions introduced magnetically and electrically can be different in addition to those for the applied magnetic and electric field. Their choices can influence the character of crack growth which could be enhanced or retarded. The details of how the directions of poling and applied field would affect crack initiation and growth are discussed in relation to the volume fraction of inclusions for a BaTiO₃–CoFe₂O₄ two phase composite. The multi-functional aspects of magnetoelectroelastic materials are involved since they entail multi-scaling features. Failure criteria that applies to isotropic elastic materials may not hold for composites exhibiting piezomagnetic and piezoelectric properties. For instance, a negative energy release rate has been obtained for cracks in piezoelectric materials.In view of what has been said with reference to the energy release rate approach, it is desirable to use the strain energy density function as a failure criterion, even if it is only for its positive definiteness character. Physically speaking, it is attractive to have a function that could rank the proportion of energy related to volume and shape change. They determine the proportion of the hard and soft phase of the composite and hence the volume fraction of the constituent. Strength and toughness parameters used for ranking isotropic and homogeneous materials will not apply for anisotropic and/or nonhomogeneous materials if these microstructure effects could not be suppressed to a lower scale and represented as an average at the macroscopic scale. Too much emphases cannot be placed on the need to clarify the multi-scaling aspects of piezoelectric and piezomagnetic materials. Their behavior as affected by the presence of crack-like defects should be understood prior to deciding whether the material characterization approach would be suitable. That is whether simplicity could justify at the expense of conceptual rigor. Much of this would depend on scaling the time and size related to loading and material structure interaction. The magnetoelectroelastic crack model selected in the work to follow perhaps will provide an insight into the complexicity of the state of affairs for treating the finer details of material behavior with rigor.The proposed test model shows that crack growth in the magnetoelectroelastic materials can be suppressed by increasing the magnitude of the piezomagnetic constants in relation to those for piezoelectricity. A more rational means of evaluating the resistance of materials against fracture is thus proposed, particularly when anisotropy and inhomogeneity might be present.     

Reflection and transmission of plane waves at an imperfectly bonded interface between piezoelectric and piezomagnetic media

The paper analyzes the reflection and transmission of plane waves at an imperfect interface between piezoelectric (PE) and piezomagnetic (PM) media. The materials are assumed to be transversely isotropic. The linear spring model is used to describe the imperfection of bonding behavior at the interface. According to this model, the properties of the interface can be characterized by the normal and tangential interfacial stiffnesses. Numerial examples are performed for BaTiO₃/CoFe₂O₄ material combination. Four cases, a perfect, slip, normal weak bonding and unbonding interfaces for the coupled quasi-pressure (QP) wave incidence from BaTiO₃ medium are compared in detail. Numerical results of the reflection and transmission coefficients (RTCs) varying with incident angle for different interfacial stiffnesses are presented. Results show that the transmitted QP/reflected QSV waves are the strongest in the perfect/unbonded cases, respectively. The scattered waves in the slip and normal weak bonding cases are between those of the perfect and unbonded cases generally. Critial angles have noticeable effect on the RTCs and energy coefficients for the coupled scattered waves of the perfect and slip cases but have a little/no effect on those of the normal weak bonding/unbonded cases. It is found that the sum of the energy carried by the transmitted/reflected QP and QSV waves is less than unit for the imperfect bonding between BaTiO₃/CoFe₂O₄ solids. These results may provide some useful reference datum for the imperfection measurment at the PE and PM interface.     

Wave band gaps in two-dimensional piezoelectric/piezomagnetic phononic crystals

In this paper, the elastic wave propagation in phononic crystals with piezoelectric and piezomagnetic inclusions is investigated taking the magneto-electro-elastic coupling into account. The electric and magnetic fields are approximated as quasi-static. The band structures of three kinds of piezoelectric/piezomagnetic phononic crystals—CoFe₂O₄/quartz, BaTiO₃/CoFe₂O₄ and BaTiO₃–CoFe₂O₄/polymer periodic composites are calculated using the plane-wave expansion method. The piezoelectric and piezomagnetic effects on the band structures are analyzed. The numerical results show that in CoFe₂O₄/quartz structures, only one narrow band gap exists along the Γ–X direction for the coupling of xy-mode and z-mode for the filling fraction f being 0.4; while in BaTiO₃/CoFe₂O₄ composites, only one narrow band gap exists along the Γ–X direction forxy-mode and no band gap exists for z-mode as the filling friction f is 0.5. Moreover, for the new type of magneto-electro-elastic phononic crystal—BaTiO₃–CoFe₂O₄/polymer periodic composite, the band gap characteristics are more superior in the whole considered frequency regions due to the big contrast of the material properties in the two constituents and the effects of the piezoelectricity and piezomagneticity on the band gap structures are remarkable.     

Piezomagnetic and piezoelectric poling effects on mode I and II crack initiation behavior of magnetoelectroelastic materials

The ferrite and ferroelectric phase of magnetoelectroelastic (MEE) material can be selected and processed to control the macroscopic behavior of electron devices using continuum mechanics models. Once macro- and/or microdefects appear, the highly intensified magnetic and electric energy localization could alter the response significantly to change the design performance. Alignment of poling directions of piezomagnetic and piezoelectric materials can add to the complexity of the MEE material behavior to which this study will be concerned with.Appropriate balance of distortional and dilatational energy density is no longer obvious when a material possesses anisotropy and/or nonhomogeneity. An excess of the former could result in unwanted geometric change while the latter may lead to unexpected fracture initiation. Such information can be evaluated quantitatively from the stationary values of the energy density function dW/dV. The maxima and minima have been known to coincide, respectively, with possible locations of permanent shape change and crack initiation regardless of material and loading type. The direction of poling with respect to a line crack and the material microstructure described by the constitutive coefficients will be specified explicitly with reference to the applied magnetic field, electric field and mechanical stress, both normal and shear. The crack initiation load and direction could be predicted by finding the direction for which the volume change is the largest. In contrast to intuition, change in poling directions can influence the cracking behavior of MEE dramatically. This will be demonstrated by the numerical results for the BaTiO₃–CoFe₂O₄ composite having different volume fractions where BaTiO₃ and CoFe₂O₄ are, respectively, the inclusion and matrix.To be emphasized is that mode I and II crack behavior will not have the same definition as that in classical fracture mechanics where load and crack extension symmetry would coincide. A striking result is found for a mode II crack. By keeping the magnetic poling fixed, a reversal of electric poling changed the crack initiation angle from θ₀=+80° to θ₀=−80° using the line extending ahead of the crack as the reference. This effect is also sensitive to the distance from the crack tip. Displayed and discussed are results for r/a=10⁻⁴ and 10⁻¹. Because the theory of magnetoelectroelasticity used in the analysis is based on the assumption of equilibrium where the influence of material microstructure is homogenized, the local space and temporal effects must be interpreted accordingly. Among them are the maximum values of (dW/dV)_max and (dW/dV)_min which refer to as possible sites of yielding and fracture. Since time and size are homogenized, it is implicitly understood that there is more time for yielding as compared to fracture being a more sudden process. This renders a higher dW/dV in contrast to that for fracture. Put it differently, a lower dW/dV with a shorter time for release could be more detrimental.     

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.     

On the crack-tip stress singularity of interfacial cracks in transversely isotropic piezoelectric bimaterials

In this paper, characteristics of the interface crack-tip stress and electric displacement fields in transversely isotropic piezoelectric bimaterials are studied. The authors have proven, within the framework of the generalized Stroh formalism for piezoelectric bimaterials, that there is no coexistence of the parameters (oscillating) and κ (non-oscillating) in the interface crack-tip generalized stress field for all transversely isotropic piezoelectric bimaterials. This leads to the classification of piezoelectric bimaterials into one group that exhibits the oscillating property in the interface crack-tip generalized stress field and the other that does not. Fifteen (15) pair-combinations of six (6) piezoelectric materials PZT-4, PZT-5H, PZT-6B, PZT-7A, P-7, and BaTiO₃, which are commonly used in practice, are numerically analyzed in this study, and the results backup the above theoretical conclusions. Moreover, the associated eigenvectors for such material systems (with either =0 or κ=0) are also obtained numerically, and the result show that there still exist four linear independent associate eigenvectors for each bimaterial.     

EXACT SOLUTION TO THE PROBLEM OF A HALF-PLANE CRACK IN A TRANSVERSELY ISOTROPIC PIEZOELECTRIC BODY SUBJECTED TO ANTISYMMETRIC TANGENTIAL POINT FORCES

An exact and complete solution of the problem of a half-plane crack in an infinite transversely isotropic piezoelectric body is presented. The upper and lower crack faces are assumed to be loaded antisymmetrically by a couple of tangential point forces in opposite directions. The solution is derived through a limiting procedure from that of a penny-shaped crack. The expressions for the electroelastic field are given in terms of elementary functions. Finally, the numerical results of the second and third mode stress intensity factorsk ₂ andk ₃ of piezoelectric materials and elastic materials are compared in figures. Project supported by the National Natural Science Foundation of China (No. 19872060 and 69982009) and the Postdoctoral Foundation of China.     

Large deflection of a rectangular magnetoelectroelastic thin plate

Based on the von Karman plate theory of large deflection, we derive the nonlinear partial differential equation for a rectangular magnetoelectroelastic thin plate under the action of a transverse static mechanical load. By employing the Bubnov-Galerkin method, the nonlinear partial differential equation is transformed to a third-order nonlinear algebraic equation for the maximum deflection where a coupling factor is introduced for determining the coupling effect on the deflection. Numerical results are carried out for the thin plate made of piezoelectric BaTiO₃ and piezomagnetic CoFe₂O₄ materials. Some interesting results are obtained which could be useful to future analysis and design of multiphase composite plates.     

Analysis of a permeable interface crack in elastic dielectric/piezoelectric bimaterials

A permeable interface crack between elastic dielectric material and piezoelectric material is studied based on the extended Stroh’s formalism. Motivated by strong engineering demands to design new composite materials, the authors perform numerical analysis of interface crack tip singularities and the crack tip energy release rates for 35 types of dissimilar bimaterials, respectively, which are constructed by five kinds of elastic dielectric materials: Epoxy, Polymer, Al₂O₃, SiC, and Si₃N₄ and seven kinds of practical piezoelectric ceramics: PZT-4, BaTiO₃, PZT-5H, PZT-6B, PZT-7A, P-7, and PZT-PIC 151, respectively. The elastic dielectric material with much smaller permittivity than commercial piezoelectric ceramics is treated as a special transversely isotropic piezoelectric material with extremely small piezoelectricity. The present investigation shows that the structure of the singular field near the permeable interface crack tip consists of three singularities: and , which is quite different from that in the impermeable interface crack. It can be concluded that different far field loading cases have significant influence on the near-tip fracture behaviors of the permeable interface crack. Based on the present theoretical treatment and numerical analysis, the electric field induced crack growth is well explained, which provides a better understanding of the failure mechanism induced from interface crack growth in elastic dielectric/piezoelectric bimaterials. The project supported by the National Natural Science Foundation of China (10572110), Doctor Foundation of the Chinese Education Ministry and Doctorate Foundation of Xi’an Jiaotong University. The English text was polished by Yunming Chen.     

压电界面裂纹裂尖强度因子的显式计算公式

对常见横观各向同性压电材料(TIP)中界面裂纹的裂纹面与压电材料的极化方向成任意夹角的一般情况进行了研究,通过推导得到了计算裂尖强度因子的显式外推公式,同时给出了裂纹面与极化方向垂直的典型情况下的外推公式.这些显式计算公式为常见数值方法如有限元法及边界元法在压电材料断裂力学中的应用提供了便利.     

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.     

A piezoelectric screw dislocation interacting with a half-plane trimaterial composite

The electro-elastic stress investigation on the interaction between a screw dislocation and a half-plane trimaterial composite composed of three bonded dissimilar transversely isotropic piezoelectric materials is analyzed in the framework of linear piezoelectricity. Each layer is assumed to have the same material orientation with x ₃ in the poling direction. The dislocations are characterized by a discontinuous displacement and electric potential across the slip plane and are subjected to a line force and a line charge at the core. Based on the complex variable and the method of alternating technique, the solution of electric field and displacement field is expressed in terms of explicit series form. The solutions derived here can be applied to a variety of problems, for example, a half-plane bimaterial, a quarter-plane bimaterial, a quarter-plane material and a rectangular strip etc. Numerical results are provided to show the influences of the material combinations and geometric configurations on the electro-elastic fields and image force calculated through the generalized Peach-Koehler formula. The solutions proposed here can be served as Green??s functions for the analyses corresponding piezoelectric cracking problems.     

A moving dislocation in a magneto–electro-elastic solid

This work considers the generalized plane problem of a moving dislocation in an anisotropic elastic medium with piezoelectric, piezomagnetic and magnetoelectric effects. The closed-form expressions for the elastic, electric and magnetic fields are obtained using the extended Stroh formalism for steady-state motion. The radial components, E_rand H_r, of the electric and magnetic fields as well as the hoop components, D_θ and B_θ, of electric displacement and magnetic flux density are found to be independent of θ in a polar coordinate system. This interesting phenomenon is proven to be is a consequence of the electric and magnetic fields, electric displacement and magnetic flux density that exhibit the singularity r⁻¹ near the dislocation core. As an illustrative example, the more explicit results for a moving dislocation in a transversely isotropic magneto–electro-elastic medium are provided and the behavior of the coupled fields is analyzed in detail.     

Inclined circular flat punch on a transversely isotropic piezoelectric half-space

Summary Utilizing the general solution of transversely isotropic piezoelectricity, the paper analyzes the problem of an inclined rigid circular flat punch indenting a transversely isotropic piezoelectric half-space. The potential theory method is employed and generalized to take into account the effect of the electric field in piezoelectric materials. Assuming that the punch is maintained at a constant electric potential, exact expressions for the elastoelectric field are derived in terms of elementary functions. It is noted that the solution corresponding to a flat circular punch centrally loaded by a concentrated force can be obtained as a special case. Received 15 December 1998; accepted for publication 9 March 1999     

SINGULAR BEHAVIORS OF INTERFACIAL CRACKS IN 2D MAGNETOELECTROELASTIC BIMATERIALS

The singular characteristics of stress, electric displacement and magnetic induction fields near the tip of impermeable interracial cracks in two-dimensional magnetoelectroelastic bimaterials are studied using the generalized Stroh formalism. Two types of singularities are obtained： one is the oscillating singularity 1/2±iε, the other is the non-oscillating singularity 1/2±κ. It is found that the non-zero parameters ε and κ cannot coexist for one transversely isotropic MEE bimaterial, a similar result is obtained for transversely isotropic piezoelectric bimaterials.     

On the generalized Barnett–Lothe tensors for monoclinic piezoelectric materials

The three generalized Barnett–Lothe tensors L, S and H, appearing frequently in the investigations of the two-dimensional deformations of anisotropic piezoelectric materials, may be expressed in terms of the material constants. In this paper, the eigenvalues and eigenvectors for monoclinic piezoelectric materials of class m, with the symmetry plane at x₃ = 0 are constructed based on the extended Stroh formalism. Then the three generalized Barnett–Lothe tensors are calculated from these eigenvectors and are expressed explicitly in terms of the elastic stiffness instead of the reduced elastic compliance. The special case of transversely isotropic piezoelectric materials is also presented.     

Rigid biological composite materials: Structural examples for biomimetic design

Biological hard tissues are composites of inorganics and biopolymers, and, therefore, represent hybrid systems. The inorganic components may be oxides (e.g., SiO₂, Fe₃O₄), carbonates (e.g., CaCO₃) sulfides (e.g., FeS, CdS), or others, mostly in crystalline forms but also occasionally in glassy forms. The biopolymer is often proteinaceous, but can also involve lipids and especially polysaccharides (e.g., chitin). These hybrid materials can be found in single celled organisms (such as bacteria and protozoa), invertebrates (such as mollusks), insects (such as beetles), and vertebrates (such as mammals). A common denominator of all hard tissues is that they are hierarchically structured from the nanometer scale to the microscale and the macroscale. It is these controlled structures that give biological hard tissues their unique and highly evolved functional properties. The engineering properties include mechanical, piezoelectric, optical, and magnetic. The hard tissues can be in the form of nanoparticles, spines, spicules, skeletons, and shells. The objective of this paper is to demonstrate mechanical aspects of some of these hard tissues, to discuss their structure-function relationships (with examples from the literature as well as from our research), and to reveal their potential utility in materials science and engineering applications.     

A domain-independent interaction integral for magneto-electro-elastic materials

Magneto-electro-elastic (MEE) materials usually consist of piezoelectric (PE) and piezomagnetic (PM) phases. Between different constituent phases, there exist lots of interfaces with discontinuous MEE properties. Complex interface distribution brings a great difficulty to the fracture analysis of MEE materials since the present fracture mechanics methods can hardly solve the fracture parameters efficiently of a crack surrounded by complex interfaces. This paper develops a new domain formulation of the interaction integral for the computation of the fracture parameters including stress intensity factors (SIFs), electric displacement intensity factor (EDIF) and magnetic induction intensity factor (MIIF) for linear MEE materials. The formulation derived here does not involve any derivatives of material properties and moreover, it can be proved that an arbitrary interface in the integral domain does not affect the validity and the value of the interaction integral. Namely, the interaction integral is domain-independent for material interfaces and thus, its application does not require material parameters to be continuous. Due to this advantage, the interaction integral becomes an effective approach for extracting the fracture parameters of MEE materials with complex interfaces. Combined with the extended finite element method (XFEM), the interaction integral is employed to solve several representative problems to verify its accuracy and domain-independence. Good results show the effectiveness of the present method in the fracture analysis of MEE materials with continuous and discontinuous properties. Finally, the particulate MEE composites composed of PE and PM phases are considered and four schemes of different property-homogenization level are proposed for comparing their effectiveness.