General solution for the coupled equations of transversely isotropic magnetoelectroelastic solids

IntroductionMechanicsandphysicsofmediapossessingsimultaneouslypiezoelectric ,piezomagneticandmagnetoelectriceffects ,namely ,magnetoelectroelasticsolids,haveattractedmoreandmoreattentionduetotheirgreatpotentialapplicationsinthetechnologiesofsmartandadaptivematerialsystem[1] .Sometheoreticalinvestigationsappearedintheliteratureinclude :1)Theexistenceproblemofsurfacewavesinsemi_infiniteanisotropicmagnetoelectroelasticmediawithvariousboundaryconditions[2 ,3 ] ;2 )Green’sfunctions[4～ 7] ;3)Inho…     

磁电弹性圆锥顶端作用集中荷载的解析解

当磁电弹性材料特征根互异时，用5个势函数表示的通解出发，对圆锥顶端作用集中扭矩Mx的扭转、集中力Px和点电荷Q的压缩、集中力Px和集中力矩My的弯曲变形问题，用一些调和函数的线性组合分别构造了势函数，并根据边界条件求出了势函数中的待定系数从而确定势函数，再将势函数代入通解得到磁电弹性圆锥顶端作用集中载荷时的位移、电势、磁势、应力、电位移和磁感应强度的三维解析解。此解形式简单便于应用。当圆锥角2α=π时，可退化得到半空间问题的解。     

The general solution of spherically isotropic magnetoelectroelastic media and its applications

A general solution of the three-dimensional equilibrium problem of spherically isotropic magnetoelectroelastic media is presented. Base on the obtained general solution, exact and compact form solutions are obtained for (1) a spherically isotropic magnetoelectroelastic cone subjected to concentrated force, concentrated couple, a point charge and a point electric current at its apex; (2) a spherically isotropic magnetoelectroelastic space with a concentrated force at the origin; (3) a spherical shell under spherically symmetric deformation; and (4) stress concentration around a spherical cavity subjected to remote uniform tensile force, electric charge and electric current.     

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.     

磁-电-弹性半空间在轴对称热载荷作用下的三维问题研究

考虑力-电-磁-热等多场耦合作用, 基于线性理论给出了磁-电-弹性半空间在表面轴对称温度载荷作用下的热-磁-电-弹性分析, 并得到了问题的解析解. 利用Hankel 积分变换法求解了磁-电-弹性材料中的热传导及控制方程, 讨论了在磁-电-弹性半空间在边界表面上作用局部热载荷时的混合边值问题, 利用积分变换和积分方程技术, 通过在边界表面上施加应力自由及磁-电开路条件, 推导得到了磁-电-弹性半空间中位移、电势及磁势的积分形式的表达式. 获得了磁-电-弹性半空间中温度场的解析表达式并且给出了应力, 电位移和磁通量的解析解. 数值计算结果表明温度载荷对磁-电-弹性场的分布有显著影响. 当温度载荷作用的圆域半径增大时, 最大正应力发生的位置会远离半无限大体的边界; 反之当温度载荷作用的圆域半径减小时, 最大应力发生的位置会靠近半无限大体的边界. 电场和磁场在温度载荷作用的圆域内在边界表面附近有明显的强化, 而磁-电-弹性场强化区域的强化程度跟温度载荷的大小和作用区域大小相关. 本研究的相关结果对智能材料和结构在热载荷作用下的设计和制造具有指导意义.      

Time-harmonic Green’s functions for anisotropic magnetoelectroelasticity

Two-dimensional (2-D) and three-dimensional (3-D) time-harmonic Green’s functions for linear magnetoelectroelastic solids are derived in this paper by means of Radon-transform. Displacement field and electric and magnetic potentials in a fully anisotropic magnetoelectroelastic infinite solid due to a time-harmonic point force, point charge and magnetic monopole are obtained in form of line integrals over a unit circle in 2-D case and surface integrals over a unit sphere in 3-D case. This dynamic fundamental solution is then split into the sum of regular dynamic plus singular terms. The singular terms coincide with the Green’s functions for the static problem and may be further reduced to closed form expressions. The proposed Green’s functions can be used in the corresponding boundary element method (BEM) formulation.     

THREE-DIMENSIONAL ANALYSIS OF A MAGNETOELECTROELASTIC HALF-SPACE UNDER AXISYMMETRIC TEMPERATURE LOADING 1)

考虑力-电-磁-热等多场耦合作用, 基于线性理论给出了磁-电-弹性半空间在表面轴对称温度载荷作用下的热-磁-电-弹性分析, 并得到了问题的解析解. 利用Hankel 积分变换法求解了磁-电-弹性材料中的热传导及控制方程, 讨论了在磁-电-弹性半空间在边界表面上作用局部热载荷时的混合边值问题, 利用积分变换和积分方程技术, 通过在边界表面上施加应力自由及磁-电开路条件, 推导得到了磁-电-弹性半空间中位移、电势及磁势的积分形式的表达式. 获得了磁-电-弹性半空间中温度场的解析表达式并且给出了应力, 电位移和磁通量的解析解. 数值计算结果表明温度载荷对磁-电-弹性场的分布有显著影响. 当温度载荷作用的圆域半径增大时, 最大正应力发生的位置会远离半无限大体的边界; 反之当温度载荷作用的圆域半径减小时, 最大应力发生的位置会靠近半无限大体的边界. 电场和磁场在温度载荷作用的圆域内在边界表面附近有明显的强化, 而磁-电-弹性场强化区域的强化程度跟温度载荷的大小和作用区域大小相关. 本研究的相关结果对智能材料和结构在热载荷作用下的设计和制造具有指导意义.     

Micromechanics predictions of the effective moduli of magnetoelectroelastic composite materials

In this paper, the self-consistent, generalized Mori–Tanaka and dilute micromechanics theories are extended to study the coupled magnetoelectroelastic composite materials. The heterogeneous inclusion problem of magnetoelectroelastic behavior is formulated in terms of five interaction tensors related to the Green's functions for an infinite three-dimensional transversely isotropic magnetoelectroelastic solid. These tensors are then used to predict the effective moduli of the magnetoelectroelastic solid based on the self-consistent, Mori–Tanaka and the dilute approaches. Numerical results are obtained for various types of inclusions. These results are employed to study the effects of the inclusion properties, such as moduli, volume fractions, shapes, etc., on the effective moduli of magnetoelectroelastic composites, in particular, the related magnetic properties. The results obtained using the self-consistent model, the generalized Mori–Tanaka's model and the dilute approach are compared with the existing experimental and theoretical results.     

Analysis method of planar cracks of arbitrary shape in the isotropic plane of a three-dimensional transversely isotropic magnetoelectroelastic medium

The hyper-singular boundary integral equation method of crack analysis in three-dimensional transversely isotropic magnetoelectroelastic media is proposed. Based on the fundamental solutions or Green’s functions of three-dimensional transversely isotropic magnetoelectroelastic media and the corresponding Somigliana identity, the boundary integral equations for a planar crack of arbitrary shape in the plane of isotropy are obtained in terms of the extended displacement discontinuities across crack faces. The extended displacement discontinuities include the displacement discontinuities, the electric potential discontinuity and the magnetic potential discontinuity, and correspondingly the extended tractions on crack face represent the conventional tractions, the electric displacement and the magnetic induction boundary values. The near crack tip fields and the intensity factors in terms of the extended displacement discontinuities are derived by boundary integral equation approach. A solution method is proposed by use of the analogy between the boundary integral equations of the magnetoelectroelastic media and the purely elastic materials. The influence of different electric and magnetic boundary conditions, i.e., electrically and magnetically impermeable and permeable conditions, electrically impermeable and magnetically permeable condition, and electrically permeable and magnetically impermeable condition, on the solutions is studied. The crack opening model is proposed to consider the real crack opening and the electric and magnetic fields in the crack cavity under combined mechanical-electric-magnetic loadings. An iteration approach is presented for the solution of the non-linear model. The exact solution is obtained for the case of uniformly applied loadings on the crack faces. Numerical results for a square crack under different electric and magnetic boundary conditions are displayed to demonstrate the proposed method.     

反平面剪切下电磁弹性纤维增强复合材料界面相效应研究

研究具有界面相电磁弹性纤维增强复合材料的反平面剪切问题,利用复变函数方法,获得了无穷域中带界面相纤维问题在远场力、电、磁多场作用下的闭合解,得到了复合材料内部各区域电磁弹性物理量的精确表达式.利用所得结果,考虑纤维和基体间的界面相效应,研究了界面相厚度及弹性模量对复合材料内部应力场、电场强度和磁场强度的影响,数值结果给出了复合材料电磁弹性物理量随界面相参数变化的规律,为该类复合材料的设计与计算提供了有价值的参考.     

Bending analysis of magnetoelectroelastic nanoplates resting on Pasternak elastic foundation based on nonlocal theory

Based on the nonlocal theory and Mindlin plate theory, the governing equations (i.e., a system of partial differential equations (PDEs) for bending problem) of magnetoelectroelastic (MEE) nanoplates resting on the Pasternak elastic foundation are first derived by the variational principle. The polynomial particular solutions corresponding to the established model are then obtained and further employed as basis functions with the method of particular solutions (MPS) to solve the governing equations numerically. It is confirmed that for the present bending model, the new solution strategy possesses more general applicability and superior flexibility in the selection of collocation points. The effects of different boundary conditions, applied loads, and geometrical shapes on the bending properties of MEE nanoplates are evaluated by using the developed method. Some important conclusions are drawn, which should be helpful for the design and applications of electromagnetic nanoplate structures.     

Modified H-R mixed variational principle for magnetoelectroelastic bodies and state-vector equation

Based upon the Hellinger-Reissner (H-R) mixed variational principle for three-dimensional elastic bodies, the modified H-R mixed variational theorem for magnetoelectroelastic bodies was established. The state-vector equation of magnetoelectroelastic plates was derived from the proposed theorem by performing the variational operations. To lay a theoretical basis of the semi-analytical solution applied with the magnetoelectroelastic plates, the state-vector equation for the discrete element in plane was proposed through the use of the proposed principle. Finally, it is pointed out that the modified H-R mixed variational principle for pure elastic, single piezoelectric or single piezomagnetic bodies are the special cases of the present variational theorem.     

FRACTURE ASSESSMENT OF AN INTERFACE CRACK BETWEEN TWO DISSIMILAR MAGNETOELECTROELASTIC MATERIALS UNDER HEAT FLOW AND MAGNETOELECTROMECHANICAL LOADINGS

A magnetoelectrically permeable interface crack between two semi-infinite magnetoelectroelastic planes under the action of a heat flow and remote magnetoelectromechanical loadings is considered, where the assumption of frictionless contact between two dissimilar halfplanes is adopted. Not only the solutions of the interface crack problem are presented in an explicit form, but also the general condition for the transition from a perfect thermal contact of two magnetoelectroelastic bodies to their separation is given.     

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.     

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.     

Analysis of cracked magnetoelectroelastic composites under time-harmonic loading

This paper presents a numerical model for the analysis of cracked magnetoelectroelastic materials subjected to in-plane mechanical, electric and magnetic dynamic time-harmonic loading. A traction boundary integral equation formulation is applied to solve the problem in combination with recently obtained time-harmonic Green’s functions (Rojas-Diaz et al., 2008). The hypersingular boundary integral equations appearing in the formulation are first regularized via a simple change of variables that permits to isolate the singularities. Relevant fracture parameters, namely stress intensity factors, electric displacement intensity factor and magnetic induction intensity factor are directly evaluated as functions of the computed nodal opening displacements and the electric and magnetic potentials jumps across the crack faces. The method is checked by comparing numerical results against existing solutions for piezoelectric solids. Finally, numerical results for scattering of plane waves in a magnetoelectroelastic material by different crack configurations are presented for the first time. The obtained results are analyzed to evaluate the dependence of the fracture parameters on the coupled magnetoelectromechanical load, the crack geometry and the characteristics of the incident wave motion.     

Exact solution for elliptical inclusion in magnetoelectroelastic materials

This paper considers the magneto-electro-mechanical coupling between an inclusion and matrix, which are both of magnetoelectroelastic materials. The general cases including the mode I, mode II and mode III are studied. Analytical solutions for an elliptical cylinder inclusion inside an infinite magnetoelectroelastic medium under combined mechanical–electrical–magnetic loadings are formulated via the Stroh formalism. Crack problem is also investigated and the stress, electric and magnetic fields in the vicinity of the crack tip are determined by a complex vector of intensity factors. Various special cases, including an impermeable inclusion, a permeable crack, a rigid and permeable inclusion, a rigid and permeable line inclusion, a permeable cavity, and an impermeable cavity are discussed.     

A PENNY-SHAPED CRACK IN A MAGNETOELECTROELASTIC LAYER UNDER RADIAL SHEAR IMPACT LOADING

This paper analyzes the dynamic magnetoelectroelastic behavior induced by a penny- shaped crack in a magnetoelectroelastic layer.The crack surfaces are subjected to only radial shear impact loading.The Laplace and Hankel transform techniques are employed to reduce the prob- lem to solving a Fredholm integral equation.The dynamic stress intensity factor is obtained and numerically calculated for different layer heights.And the corresponding static solution is given by simple analysis.It is seen that the dynamic stress intensity factor for cracks in a magnetoelec- troelastic layer has the same expression as that in a purely elastic material.And the influences of layer height on both the dynamic and static stress intensity factors are insignificant as h/a>2.     

Stress state of a transversely isotropic magnetoelectroelastic body with a plane crack under antisymmetric loads

The static equilibrium of a transversely isotropic magnetoelectroelastic body with a plane crack of arbitrary shape in the isotropy plane under antisymmetric mechanical loading is studied. The relationships between the stress intensity factors (SIFs) for an infinite magnetoelectroelastic body and the SIFs for a purely elastic body with the same crack and under the same antisymmetric loading are established. This enables the SIFs for a magnetoelectroelastic body to be found directly from the analogous problem of elasticity. As an example of using this result, the SIFs for penny-shaped, elliptic, and parabolic cracks in a magnetoelectroelastic body under antisymmetric mechanical loading are found Translated from Prikladnaya Mekhanika, Vol. 44, No. 10, pp. 37–51, October 2008.     

<Emphasis Type="Italic">T</Emphasis>-stress analysis for a Griffith crack in a magnetoelectroelastic solid

T-stress as an important parameter characterizing the stress field around a cracked tip has attracted much attention. This paper concerns the T-stress near a cracked tip in a magnetoelectroelastic solid. By applying the Fourier transform, we solve the associated mixed boundary-value problem. Adopting crack-faces electromagnetic boundary conditions nonlinearly dependent on the crack opening displacement, coupled dual integral equations are derived. Then, the closed-form solution for the T-stress is obtained. A comparison of the T stresses for a cracked magnetoelectroelastic solid and for a cracked purely elastic material is made. Obtained results reveal that in addition to applied mechanical loading, the T-stress is dependent on electric and magnetic loadings for a vacuum crack.