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.     

Interfacial fracture analysis of a graded piezoelectric layer on a substrate with finite dimension

In fracture analysis of piezoelectric devices, the structural dimension is often assumed to be infinite at least in one direction. However, all practical piezoelectric structures are finite and their dimensions in different directions are often comparable and cannot be simplified as infinite. The assumption of infinite dimension may lead to inexact theoretical results. The present work aims at studying the interfacial fracture behavior of a functionally graded piezoelectric layer on a dielectric substrate with finite dimension. The crack problem is solved by the methods of Fourier series and Cauchy singular integral equation. Parametric studies on the stress intensity factor (SIF) reveal the following: (a) when a crack tip is near to an interface end, its SIF is mainly governed by the end effect; (b) when a crack is far from the interface ends and the piezoelectric layer is thin, its SIF is principally affected by the thickness of the piezoelectric layer, and (c) only when a crack is far from the interface ends and meanwhile the piezoelectric layer is thick will its SIF be dominated by the non-homogeneity parameter, and in this case, the SIF increases with the increasing non-homogeneity parameter.     

Exact solutions for magneto-electro-elastic laminates in cylindrical bending

Analytical solutions are derived for the cylindrical bending of multilayered, linear, and anisotropic magneto-electro-elastic plates under simple-supported edge conditions. We construct the general solution in terms of a simple formalism for any homogeneous layer, from which any physical quantities can be solved for the given boundary conditions. For multilayered plates, we derive the solution in terms of the propagator matrices. A special feature of cylindrical bending, which distinguishes itself from the three-dimensional plate problem, is that the associated eigenvalues for any homogeneous layer are independent of the sinusoidal mode, and thus need to be solved only once. Typical numerical examples are also presented for a piezomagnetic plate, a two-layered piezoelectric/piezomagnetic plate, and a four layered piezoelectric/piezomagnetic plate, with different span-to-thickness ratios. In particular, the piezoelectric and piezomagnetic fields show certain interesting features, which give guidance on the development of piezoelectric/piezomagnetic thin-plate theories. Furthermore, it is shown that the variations of the elastic, electric, and magnetic quantities with thickness depend strongly upon the material property and layering, which could be useful in the analysis and design of smart composite structures with sensors/actuators.     

Propagation behavior of SH waves in a piezomagnetic substrate with an orthorhombic piezoelectric layer

The dispersion behavior of the shear horizontal (SH) waves in the coupled structure consisting of a piezomagnetic substrate and an orthorhombic piezoelectric layer is investigated with different cut orientations. The surface of the piezoelectric layer is mechanically free, electrically shorted, or open, while the surface of the piezomagnetic substrate is mechanically free, magnetically open, or shorted. The dispersion relations are derived for four electromagnetic boundary conditions. The dispersion characteristics are graphically illustrated for the layered structure with the PMN-PT layer perfectly bonded on the CoFe₂O₄ substrate. The effects of the PMN-PT cut orientations, the electromagnetic boundary conditions, and the thickness ratio of the layer to the substrate on the dispersion behavior are analyzed and discussed in detail. The results show that, (i) the effect of the cut orientation on the dispersion curves is very obvious, (ii) the electrical boundary conditions of the PMN-PT layer dominate the propagation feature of the SH waves, and (iii) the thickness ratio has a significant effect on the phase velocity when the wave number is small. The results of the present paper can provide valuable theoretical references to the applications of piezoelectric/piezomagnectic structure in acoustic wave devices.     

一维磁电弹准周期结构带隙特性的研究

采用传递矩阵方法,研究了横波(SV波)垂直入射时压电/(弹性/压磁)和(压电/弹性)/压磁两种Fibonacci准周期结构的频带特性,通过计算局部化因子和位移透射系数,数值揭示了此两种Fibonacci准周期结构频带特性的差异以及与相应理想周期结构频带特性的不同,而且表明(压电/弹性)/压磁Fibonacci准周期结构的频带特性与纯弹性材料Fibonacci准周期结构的频带特性是相似的.     

弱界面压电／压磁层状结构中的SH波

研究了弱界面压电层/压磁半空间结构中SH波的频散特性, 界面性能由 ``弹簧'模型表征, 压电层的表面考虑电学短路和电学开路两种边界条件. 推导了显函形式的 频散方程, 并结合算例分析了界面性能和电学边界条件对SH波传播特性的影响. 数值结果表 明: 弱界面降低了SH波传播的相速度; 对于短路电边界条件, 弱界面可以使一阶模态的相速 度小于B-G表面波的速度; 对于开路电边界条件, 一阶模态的相速度总是大于剪切体波的波 速.     

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.     

压电/压磁层合纳米梁屈曲及自由振动的研究 1)

针对压电/压磁层合纳米梁屈曲、自由振动问题,基于非局部理论与正弦剪切型变形梁理论,建立了力学模型;利用哈密顿原理推导出层合梁运动方程与边界条件;通过数值解法求得层合梁临界屈曲载荷与自由振动频率。对数值结果分析可知：磁电弹夹层对压电/压磁层合纳米梁屈曲和自由振动的影响不能忽略;磁电弹夹层中压电或压磁材料的体积分数和夹层厚度为主要影响因素;分析得到的影响规律可为此类材料在工程中的应用提供理论参考。     

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.     

Solution of two-dimensional scattering problem in piezoelectric/piezomagnetic media using a polarization method

Using a polarization method, the scattering problem for a two-dimensional inclusion embedded in infinite piezoelectric/piezomagnetic matrices is investigated. To achieve the purpose, the polarization method for a two-dimensional piezoelectric/piezo-magnetic "comparison body" is formulated. For simple harmonic motion, kernel of the polarization method reduces to a 2-D time-harmonic Green's function, which is ob-tained using the Radon transform. The expression is further simplified under condi-tions of low frequency of the incident wave and small diameter of the inclusion. Some analytical expressions are obtained. The analytical solutions for generalized piezoelec-tric/piezomagnetic anisotropic composites are given followed by simplified results for piezoelectric composites. Based on the latter results, two numerical results are provided for an elliptical cylindrical inclusion in a PZT-5H-matrix, showing the effect of different factors including size, shape, material properties, and piezoelectricity on the scattering cross-section.     

Why traction-free? Piezoelectric crack and Coulombic traction

The non-zero traction condition is introduced in piezoelectric crack problems with the unknown Coulombic traction acting on the crack surfaces. An analytical solution under this condition is obtained by means of the generalized Stroh formalism and by accounting for the permittivity of medium inside the crack gap. As the crack in such materials can be thought of as a low-capacitance medium carrying a potential drop, the Coulombic traction always pulls the two opposite surfaces of the crack together. It is proved that under relatively larger mechanical loading and relatively smaller electrical field, the Coulombic traction may be negligible and the previous investigations under the traction-free crack condition may be accepted in a tolerant way, otherwise the Coulombic traction may lead to some erroneous results with over 10% relative errors. It is also shown that, unlike the traction-free crack condition, the applied electric field does change the Mode I stress intensity factor (SIF) for a central crack in an infinite plane piezoelectric material, and in this way may significantly influence piezoelectric fracture. It is also concluded that the variable tendencies of the normalized SIF and the ERR against the applied electric field depend on the mechanical loading levels. This load-dependence feature may lead to a transformation of the normalized SIF and the ERR from an even functional dependence to an odd functional dependence on the applied electric field.     

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.     

A domain-independent interaction integral for fracture analysis of nonhomogeneous piezoelectric materials

Piezoelectric materials and structures contain more or less electromechanical interfaces in engineering applications. It is difficult to obtain the fracture parameters efficiently of the piezoelectric materials with complex interfaces. This paper presents a domain-independent interaction integral for material nonhomogeneity and discontinuity which can be used for solving the stress intensity factors (SIFs) and the electric displacement intensity factor (EDIF) of piezoelectric materials with complex interfaces efficiently. The interaction integral is based on the J-integral by superimposition of two admissible states and the present formulation does not involve any derivatives of mechanical and electric properties. Moreover, it is proved that the interface in the integral domain does not affect the value of the interaction integral and thus, the present method does not require electromechanical parameters of piezoelectric materials to be continuous. The interaction integral method combined with the extended finite element method (XFEM) is used to investigate the influences of material continuity on the SIF and the EDIF and the results show that the material parameters and their first-order derivatives affect both the SIF and the EDIF greatly, while the higher-order derivatives affect both of them slightly.     

Bending of multiferroic laminated rectangular plates with imperfect interlaminar bonding

The bending problem of a multiferroic rectangular plate with magnetoelectric coupling and imperfect interfaces is investigated via three-dimensional exact theory. A generalized spring layer model is proposed to characterize the imperfection of the bonding behavior at interfaces. In particular, the linear relation between the electric displacement and the jump of electric potential, the corresponding one for the magnetic field as well as linear relations among different physical fields are adopted. State space formulations are established, which, compared to the analysis for perfect laminates, only introduces a so-called interfacial transfer matrix. The present analysis can be readily used for the piezoelectric, piezomagnetic and elastic laminates by setting the proper material constants as zero. Numerical results are presented and discussed.     

反平面剪切状态下压电-压磁夹层结构的圣维南端部效应

研究了反平面状态下压电-压磁夹层结构的圣维南端部效应的衰减特征.夹层结构端部的力-电-磁载荷是自平衡的;外层的边界是应力自由的,但承受四种不同的电磁边界条件,即电学短路/磁学短路、电学短路/磁学开路、电学开路/磁学开路和电学开路/磁学短路.基于磁电弹全耦合的控制方程,推导了问题的精确解和关于衰减率的特征方程.利用得到的特征方程,计算了压电-压磁-压电和压磁-压电-压磁两种夹层结构的衰减率,结果表明电磁边界条件、压电层的材料性能和体积分数对衰减特征有显著的影响.作为特例,分析了压电-弹性-压电夹层结构的衰减特征.     

恒磁场作用下压磁/压电半导体复合圆柱壳的耦合响应分析

本文针对由压磁材料和压电半导体材料组成的无限长复合圆柱壳结构,理论研究了其在径向恒磁场作用下结构内的多场耦合力学响应问题。为解析求解方便,文中分别采用单向耦合法和线性全耦合法,导出复合圆柱壳内位移场、电场、载流子等物理量的解析表达式。利用导出的解析表达式,数值分析了磁场大小和半导体层厚度比对结构内电势、电场和载流子的影响。计算结果表明：磁场和厚度比均可用来有效调控半导体层内的电学量,复合圆柱壳结构的厚度比有一个最优区间。     

恒磁场作用下压磁/压电半导体复合圆柱壳的耦合响应分析

论文针对由压磁材料和压电半导体材料组成的无限长复合圆柱壳结构,理论研究了其在径向恒磁场作用下结构内的多场耦合力学响应问题.为解析求解方便,文中分别采用单向耦合法和线性全耦合法,导出复合圆柱壳内位移场、电场、载流子等物理量的解析表达式.利用导出的解析表达式,数值分析了磁场大小和半导体层厚度比对结构内电势、电场和载流子的影响.计算结果表明:磁场和厚度比均可用来有效调控半导体层内的电学量,复合圆柱壳结构的厚度比有一个最优区间.     

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.     

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.     

压电压磁复合材料中裂纹对反平面简谐弹性波的散射问题

分析了压电压磁复合材料中裂纹对反平面简谐弹性波的散射问题。利用傅立叶变换,使问题的求解转换为对一对以裂纹表面上的位移差为未知变量的对偶积分方程的求解。为了求解对偶积分方程,把裂纹面上的位移差展开为雅可比多项式形式,进而得到了裂纹长度、入射波波速及入射波频率对裂纹应力强度因子的影响。从数值结果可以看出,压电压磁复合材料中可导通裂纹的反平面问题的动应力奇异性与一般弹性材料中的反平面断裂问题动应力奇异性相同。