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.     

Line-integral representations for extended displacements,stresses, and interaction energy of arbitrary dislocation loops in transversely isotropic magneto-electro-elastic bimaterials

In addition to the hexagonal crystals of class 6 mm, many piezoelectric materials （e.g., BaTiO3）, piezomagnetic materials （e.g., CoFe2O4）, and multiferroic com-posite materials （e.g., BaTiO3-CoFe2O4 composites） also exhibit symmetry of transverse isotropy after poling, with the isotropic plane perpendicular to the poling direction. In this paper, simple and elegant line-integral expressions are derived for extended displace-ments, extended stresses, self-energy, and interaction energy of arbitrarily shaped, three-dimensional （3D） dislocation loops with a constant extended Burgers vector in trans-versely isotropic magneto-electro-elastic （MEE） bimaterials （i.e., joined half-spaces）. The derived solutions can also be simply reduced to those expressions for piezoelectric, piezo-magnetic, or purely elastic materials. Several numerical examples are given to show both the multi-field coupling effect and the interface/surface effect in transversely isotropic MEE materials.     

Crack tip field in piezoelectric/piezomagnetic media

This paper considers the magnetoelectroelastic problem of a crack in a medium possessing coupled piezoelectric, piezomagnetic and magnetoelectric effects. Based on the extended Stroh formalism, the general two-dimensional solutions to the magnetoelectroelastic problem are obtained, involving five analytic functions of different variables. The magnetoelectroelastic field around the crack tip is given. It contains five modes of square root singularities. Expressions of the stresses, electric displacements and magnetic inductions in the vicinity of the crack tip are derived and the field intensity factors are provided. The path-independent conservative integral is derived. The energy release rate is written in terms of those field intensity factors. The explicit algebraic results are given for a special case of an anti-plane crack in a magnetoelectroelastic medium.     

Magnetoelectroelastic composite with poling parallel to plane of line crack under out-of-plane deformation

The three-dimensional field equations can in general be regarded as the sum of in-plane and out-of-plane deformation. The method for the general solution is the same for both although the boundary conditions could make a difference. If a particular solution in exact form may be found for the out-of-plane case, the same may not hold for the in-plane case. Hence, there may be a good reason for discussing the out-of-plane crack problem in certain situations that should be emphasized. Otherwise, the reason may lie in the exploration of possible application to the in-plane problem, a direct solution of which would have required a considerable effort. The contribution of this work rests on the new findings for the case of poling parallel to the crack in a magnetoelectroelastic composite made of BaTiO₃–CoFe₂O₄. The inclusions are BaTiO₃ and the matrix is CoFe₂O₄. Several new features of the solution were not expected before hand.Unlike in-plane deformation with poling normal to the crack plane, maximum crack growth enhancement is found to occur in the BaTiO₃–CoFe₂O₄ composite for a volume fraction of about 50%. Crack retardation increases as the volume fraction of the inclusions either increase or decrease. The occurrence of this same phenomenon in Mode I and II remain to be investigated. Poling direction of magnetic and electric field for line defects can have a significant effect on crack growth for magnetoelectroelastic materials. The foregoing conclusions are based on predictions made from the strain energy density criterion.     

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.     

Antiplane internal crack normal to the edge of a functionally graded piezoelectric/piezomagnetic half plane

This paper deals with the antiplane magnetoelectroelastic problem of an internal crack normal to the edge of a functionally graded piezoelectric/piezomagnetic half plane. The properties of the material such as elastic modulus, piezoelectric constant, dielectric constant, piezomagnetic coefficient, magnetoelectric coefficient and magnetic permeability are assumed in exponential forms and vary along the crack direction. Fourier transforms are used to reduce the impermeable and permeable crack problems to a system of singular integral equations, which is solved numerically by using the Gauss-Chebyshev integration technique. The stress, electric displacement and magnetic induction intensity factors at the crack tips are determined numerically. The energy density theory is applied to study the effects of nonhomogeneous material parameter β, edge conditions, location of the crack and load ratios on the fracture behavior of the internal crack.     

THE BASIC SOLUTION OF A 3-D RECTANGULAR PERMEABLE CRACK IN A PIEZOELECTRIC/PIEZOMAGNETIC COMPOSITE MATERIAL

The solution of a 3-D rectangular permeable crack in a piezoelectric/piezomagnetic composite material was investigated by using the generalized Almansi’s theorem and the Schmidt method.The problem was formulated through Fourier transform into three pairs of dual integral equations,in which the unknown variables are the displacement jumps across the crack surfaces.To solve the dual integral equations,the displacement jumps across the crack surfaces were directly expanded as a series of Jacobi polynomials.Finally,the relations between the electric filed,the magnetic flux field and the stress field near the crack edges were obtained and the efects of the shape of the rectangular crack on the stress,the electric displacement and magnetic flux intensity factors in a piezoelectric/piezomagnetic composite material were analyzed.     

Fracture of piezoelectromagnetic materials

The crack problem in a medium possessing coupled piezoelectric, piezomagnetic and magnetoelectric effects is considered. A conservative integral is derived based on the governing equations for magnetoelectroelastic media. Closed-form solution is obtained for an anti-plane crack in an infinite medium. The conservative integral is used to obtain the path-independent integral near the crack tip. Expressions for stresses, electric displacements and magnetic inductions in the vicinity of a crack tip are derived. It is found that the path-independent integral around the crack tip equals the energy release rate. In the absence of applied mechanical loads, the energy release rate is always negative.     

Strip electric-magnetic breakdown model in a magnetoelectroelastic medium

Extending the polarization saturation model [Gao et al., 1997. Local and global energy release rates for an electrically yielded crack in a piezoelectric ceramic. J. Mech. Phys. Solids 45, 491-510] and the dielectric breakdown (DB) model [Zhang et al., 2005. The strip dielectric breakdown model. Int. J. Fract. 132, 311-327] in piezoelectric materials, the Strip Electric-Magnetic Breakdown (SEMB) model is proposed for electrically and magnetically impermeable crack in a magnetoelectroelastic medium to study the effect of the nonlinear character of electric field and magnetic field on fracture of magnetoelectroelastic materials. In the SEMB model, the electric field in the strip of the electric breakdown zone ahead of the crack tip is equal to the electric breakdown strength, while the magnetic filed in the strip of the magnetic breakdown zone is equal to the magnetic breakdown strength. By using the extended Stroh formalism and the extended dislocation modeling of a crack, the Griffith crack problem under the electrically and magnetically elastic-plastic condition in a magnetoelectroelastic medium is reduced to a set of dual integral equations. The sizes of the electric breakdown zone and the magnetic breakdown zone, the extended intensity factors and the local J-integral are obtained. The effect of the combined mechanical-electric-magnetic loadings on the local J-integral is studied.     

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.     

A generalized screw dislocation interacting with interfacial cracks along a circular inhomogeneity in magnetoelectroelastic solids

The interaction of a generalized screw dislocation with circular arc interfacial cracks under remote antiplane shear stresses, in-plane electric and magnetic loads in transversely isotropic magnetoelectroelastic solids is dealt with. By using the complex variable method, the general solutions to the problem are presented. The closed-form expressions of complex potentials in both the inhomogeneity and the matrix are derived for a single circular-arc interfacial crack. The intensity factors of stress, electric displacement and magnetic induction are provided explicitly. The image forces acting on the dislocation are also calculated by using the generalized Peach–Koehler formula. For the case of piezoelectric matrix and piezomagnetic inclusion, the shielding and anti-shielding effect of the dislocation upon the stress intensity factors is evaluated in detail. The results indicate that if the distance between the dislocation and the crack tip remains constant, the dislocation in the interface will have a largest shielding effect which retards the crack propagation. In addition, the influence of the interfacial crack geometry and materials magnetoelectroelastic mismatch upon the image force is discussed. Numerical computations show that the perturbation effect of the above parameters upon the image force is significant. The main result shows that a stable or unstable equilibrium point may be found when a screw dislocation approaches the surface of the crack from infinity which differs from the perfect bonded case under the same conditions. The present solutions contain a number of previously known results which can be shown to be special cases.     

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.     

Magnetostrictive/electrostrictive fracture of the piezomagnetic and piezoelectric layers in a multiferroic composite: Anti-plane case

The main purpose of the present work is to study the influences of magnetostriction, electrostriction and piezomagnetic/piezoelectric stiffening on the fracture behavior of a layered multiferroic composite. For comparison, it is assumed that there is a crack, parallel to the interface, in each layer. Methods of cosine transform and Cauchy singular integral equations are used to solve the crack problem. Numerical results of the stress intensity factor (SIF) are provided and the computational accuracy is demonstrated. Discussion on the numerical results indicates that the multiferroic composite consisting of cobalt ferrite and barium titanate layers are more prone to fracture under electric loading than under magnetic loading. In the case of magnetostriction, to increase the shear modulus of the piezomagnetic layer would raise the SIF; but to increase that of the piezoelectric layer would reduce the SIF; in the case of electrostriction, inverse results are obtained. Piezomagnetic stiffening can affect the SIF when the composite is under electrostriction; piezoelectric stiffening can influence the SIF if the composite is under magnetostriction. In addition, it is also revealed that two parallel equal cracks may shield each other even if an interface exists between them.     

Fracture of magnetoelectroelastic composite materials using boundary element method (BEM)

The behavior of cracked linear magnetoelectroelastic solids is analysed by means of the dual Boundary Element Method (BEM) approach. Media possessing fully coupled piezoelectric, piezomagnetic and magnetoelectric effects are considered. An explicit 2-D Green’s function in terms of the extended Stroh formalism for magnetoelectroelastic full-plane under static loading is implemented. Hypersingular integrals arising in the traction boundary integral equations are computed through a regularization technique. Evaluation of fracture parameters directly from computed nodal values is discussed. The stress intensity factors (SIF), the electric displacement intensity factor (EDIF), the magnetic induction intensity factor (MIIF) as well as the mechanical strain energy release rate (MSERR) are evaluated for different crack configurations in both finite and infinite solids subjected to in-plane combined magnetic–electric–mechanical loading conditions. The accuracy of the boundary element solution is confirmed by comparison with selected analytical solutions in the literature. The new results that can be of interest in the design and maintenance of novel magnetoelectroelastic devices are also discussed.     

压电材料中的Bueckner功共轭积分及和J积分M积分的关系

研究横观各向同性压电材料中裂纹问题,提出了Bueckner功共轭积分在这类材料中的表达式：并通过引出两类辅助的应力-位移-电位移-电势场,证明功共轭积分和这类材料中的J积分和M积分仍然存在简单的两倍关系由此,各类在脆性材料断裂问题中已广泛应用的权函数方法可顺理成章地推广到压电材料的研究中来．这对独立地确定电位移强度因子和经典的I、II型应力强度因子提供了有力的数学上的工具．进而通过计算机械应变能释放率对压电材料中裂纹的稳定做出判断．     

Applicability of the crack-face electromagnetic boundary conditions for fracture of magnetoelectroelastic materials

This paper discusses the different electromagnetic boundary conditions on the crack-faces in magnetoelectroelastic materials, which possess coupled piezoelectric, piezomagnetic and magnetoelectric effects. A notch of finite thickness in these materials containing air (or vacuum) is also addressed. Four ideal crack-face electromagnetic boundary condition assumptions, that is, (a) electrically and magnetically impermeable crack, (b) electrically impermeable and magnetically permeable crack, (c) electrically permeable and magnetically impermeable crack and (d) electrically and magnetically permeable crack, are investigated separately. The influence of notch thickness on the field intensity factors at notch tips and the electromagnetic field inside the notch are obtained in closed-form. The results are compared with the ideal crack solutions. Applicability of crack-face electromagnetic boundary condition assumptions is discussed.     

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.     

MULTIPLE PARALLEL SYMMETRIC PERMEABLE MODEL-Ⅲ CRACKS IN A PIEZOELECTRIC/PIEZOMAGNETIC COMPOSITE MATERIAL PLANE

In this paper, the interactions of multiple parallel symmetric and permeable finite length cracks in a piezoelectric/piezomagnetic material plane subjected to anti-plane shear stress loading are studied by the Schmidt method.The problem is formulated through Fourier transform into dual integral equations, in which the unknown variables are the displacement jumps across the crack surfaces.To solve the dual integral equations, the displacement jumps across the crack surfaces are directly expanded as a series of Jacobi polynomials.Finally, the relation between the electric field, the magnetic flux field and the stress field near the crack tips is obtained.The results show that the stress, the electric displacement and the magnetic flux intensity factors at the crack tips depend on the length and spacing of the cracks.It is also revealed that the crack shielding effect presents in piezoelectric/piezomagnetic materials.     

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.     

Coupled field state around three parallel non-symmetric cracks in a piezoelectric/piezomagnetic material plane

In this paper, the behavior of three parallel non-symmetric permeable cracks in a piezoelectric/piezomagnetic material plane subjected to anti-plane shear stress loading was studied by the Schmidt method. The problem was formulated through Fourier transform into three pairs of dual integral equations, in which unknown variables are jumps of displacements across the crack surfaces. To solve the dual integral equations, the jumps of displacements across the crack surfaces were directly expanded as a series of Jacobi polynomials. Finally, the relations among the electric displacement, the magnetic flux and the stress fields near the crack tips can be obtained. The results show that the stress, the electric displacement and the magnetic flux intensity factors at the crack tips depend on the lengths and spacing of cracks. It was also revealed that the crack shielding effect is present in piezoelectric/piezomagnetic materials.