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

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

Pre-kinking of a moving crack in a magnetoelectroelastic material under in-plane loading

A constant moving crack in a magnetoelectroelastic material under in-plane mechanical, electric and magnetic loading is studied for impermeable crack surface boundary conditions. Fourier transform is employed to reduce the mixed boundary value problem of the crack to dual integral equations, which are solved exactly. Steady-state asymptotic fields near the crack tip are obtained in closed form and the corresponding field intensity factors are expressed explicitly. The crack speed influences the singular field distribution around the crack tip and the effects of electric and magnetic loading on the crack tip fields are discussed. The crack kinking phenomena is investigated using the maximum hoop stress intensity factor criterion. The magnitude of the maximum hoop stress intensity factor tends to increase as the crack speed increases.     

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.     

<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.     

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.     

Two-dimensional generalized thermoelasticity problem for a half-space subjected to ramp-type heating

In this paper, we will consider a half-space filled with an elastic material, which has constant elastic parameters. The governing equations are taken in a unified system from which the field equations for coupled thermoelasticity as well as for generalized thermoelasticity can be easily obtained as particular cases. A linear temperature ramping function is used to more realistically model thermal loading of the half-space surface. The medium is assumed initially quiescent. Laplace and Fourier transform techniques are used to obtain the general solution for any set of boundary conditions. The general solution obtained is applied to a specific problem of a half-space subjected to ramp-type heating. The inverse Fourier transforms are obtained analytically while the inverse Laplace transforms are computed numerically using a method based on Fourier expansion techniques. Some comparisons have been shown in figures to estimate the effect of the ramping parameter of heating with different theories of thermoelasticity.     

Dynamic fracture analysis of a penny-shaped crack in a magnetoelectroelastic layer

This paper analyzes the dynamic magnetoelectroelastic behavior induced by a penny-shaped crack in a magnetoelectroelastic layer subjected to prescribed stress or prescribed displacement at the layer surfaces. Two kinds of crack surface conditions, i.e., magnetoelectrically impermeable and permeable cracks, are adopted. The Laplace and Hankel transform techniques are employed to reduce the problem to Fredholm integral equations. Field intensity factors are obtained and discussed. Numerical results of the crack opening displacement (COD) intensity factors are presented and the effects of magnetoelectromechanical loadings, crack surface conditions and crack configuration on crack propagation and growth are examined. The results indicate that among others, the fracture behaviors of magnetoelectroelastic materials are affected by the sizes and directions of the prescribed magnetic and/or electric fields, and the effects are strongly dependent on the elastic boundary conditions.     

Semi-permeable crack analysis in a 2D magnetoelectroelastic medium based on the EMPS model

For a crack in a magnetoelectroelastic plane under the electrically and magnetically semi-permeable boundary condition, we derive the non-linear analytical solution of the strip electric–magnetic polarization saturation (EMPS) model. Using the extended dislocation theory and integral equation method, we obtain the electric and magnetic yielding zones, as well as the field intensity factor and local J-integral. Adapting an iterative method, numerical examples were performed to analyze the effect of different boundary conditions and the electric–magnetic saturated properties on the electric displacement and magnetic induction in the crack cavity, electric and magnetic yielding zones, stress intensity factor and local J-integral.     

Theoretical model of crack branching in magnetoelectric thermoelastic materials

Thermomagnetoelectroelastic crack branching of magnetoelectro thermoelastic materials is theoretically investigated based on Stroh formalism and continuous distribution of dislocation approach. The crack face boundary condition is assumed to be fully thermally, electrically and magnetically impermeable. Explicit Green’s functions for the interaction of a crack and a thermomagnetoelectroelastic dislocation (i.e., a thermal dislocation, a mechanical dislocation, an electric dipole and a magnetic dipole located at a same point) are presented. The problem is reduced to two sets of coupled singular integral equations with the thermal dislocation and magnetoelectroelastic dislocation densities along the branched crack line as the unknown variables. As a result, the formulations for the stress, electric displacement and magnetic induction intensity factors and energy release rate at the branched crack tip are expressed in terms of the dislocation density functions and the branch angle. Numerical results are presented to study the effect of applied thermal flux, electric field and magnetic field on the crack propagation path by using the maximum energy release rate criterion.     

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 analysis for an opening crack in a piezoelectromagnetic solid

Magnetoelectroelastic analysis of a cracked piezoelectromagnetic solid is made within the framework of the theory of linear magnetoelectroelasticity. The associated mixed boundary-value problem is solved by the Fourier integral transform. For general electromagnetic crack-face boundary conditions, a full magnetoelectroelastic field in the entire plane induced by a crack is obtained explicitly, and field intensity factors and energy release rate are given. The influences of applied electric and magnetic loadings on the energy release rate, the strain intensity factor, and the stress distribution are presented graphically.     

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

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

电磁弹性复合材料双圆柱夹杂问题

研究了双圆柱压电夹杂嵌于无限大压磁基体中的力学问题，获得了基体受到无究远处电，磁，力载荷作用下的解析解。利用复变函数呆角变换和解析延拓，以及圆环域内的洛朗级数展开和Cauchy积分公式等得到了基体和两个夹杂中的复位势。由所获得的复位势继而得到了应力，电位移和磁通密度等物理量的解析表达式。并给出了算例分析，以此来表明夹杂对系统电磁弹性耦合行为的影响，以及验证所得解的正确性和实用性。     

Dynamic analysis of a cracked magnetoelectroelastic medium under antiplane mechanical and inplane electric and magnetic impacts

The transient analysis of a magnetoelectroelastic medium containing a crack is made under antiplane mechanical and inplane electric and magnetic impacts. The crack is assumed to penetrate through the solid along the poling direction. By using the Fourier and Laplace transforms, the associated mixed boundary value problem is reduced to a Fredholm integral equation of the second kind, which is solved numerically. By means of a numerical inversion of the Laplace transform, dynamic field intensity factors are obtained in the time domain. Numerical results are presented graphically to show the effects of the material properties and applied electric and magnetic impacts on the dynamic intensity factors of COD and stress, and dynamic energy density factors. The results indicate that except for the intensity factors of electric displacement and magnetic induction, other field intensity factors exhibit apparent transient feature. Moreover, they depend strongly on mechanical input as well as electric and magnetic impacts.     

A Hybrid Boundary/Finite Element Method for Simulating Viscous Flows and Shapes of Droplets in Electric Fields

A mixed boundary element and finite element numerical algorithm for the simultaneous prediction of the electric fields, viscous flow fields, thermal fields and surface deformation of electrically conducting droplets in an electrostatic field is described in this paper. The boundary element method is used for the computation of the electric potential distribution. This allows the boundary conditions at infinity to be directly incorporated into the boundary integral formulation, thereby obviating the need for discretization at infinity. The surface deformation is determined by solving the normal stress balance equation using the weighted residuals method. The fluid flow and thermal fields are calculated using the mixed finite element method. The computational algorithm for the simultaneous prediction of surface deformation and fluid flow involves two iterative loops, one for the electric field and surface deformation and the other for the surface tension driven viscous flows. The two loops are coupled through the droplet surface shapes for viscous fluid flow calculations and viscous stresses for updating the droplet shapes. Computing the surface deformation in a separate loop permits the freedom of applying different types of elements without complicating procedures for the internal flow and thermal calculations. Tests indicate that the quadratic, cubic spline and spectral boundary elements all give approximately the same accuracy for free surface calculations; however, the quadratic elements are preferred as they are easier to implement and also require less computing time. Linear elements, however, are less accurate. Numerical simulations are carried out for the simultaneous solution of free surface shapes and internal fluid flow and temperature distributions in droplets in electric fields under both microgravity and earthbound conditions. Results show that laser heating may induce a non-uniform temperature distribution in the droplets. This non-uniform thermal field results in a variation of surface tension along the surface of the droplet, which in turn produces a recirculating fluid flow in the droplet. The viscous stresses cause additional surface deformation by squeezing the surface areas above and below the equator plane.     

A periodic array of cracks in a transversely isotropic magnetoelectroelastic material

Magnetoelectroelastic materials are inherently brittle and prone to fracture. Therefore, it is important to evaluate the fracture behavior of these advanced materials. In this paper, a periodic array of cracks in a transversely isotropic magnetoelectroelastic material is investigated. Hankel transform is applied to solve elastic displacements, electric potential and magnetic potential. The problem is reduced into a system of integral equations. Both impermeable and permeable crack-face electromagnetic boundary conditions assumptions are investigated. Quantities of the stress, electric displacement and magnetic induction and their intensity factor are obtained. Effect of the crack spacing on these quantities is investigated in details.     

Magneto-thermo-elastic stresses induced by a transient magnetic field in a conducting solid circular cylinder

In the present paper, dynamic and quasi-static behaviors of magneto-thermo-elastic stresses induced by a transient magnetic field in a conducting solid circular cylinder are investigated. It is assumed that a transient magnetic field which is defined by an arbitrary function of time acts on the surface of the solid cylinder in the direction parallel to its surface. Fundamental equations of plane axisymmetrical electromagnetic, temperature, and elastic fields are formulated. Then, solutions of magnetic field, eddy current, temperature change and both dynamic solutions and quasi-static ones of stresses and deformations are analytically derived in the forms including the arbitrary function. The solutions of stresses are determined to be sums of thermal stress caused by eddy current loss and magnetic stress caused by Lorentz force. For this case that the arbitrary function is given by the smoothed ramp function with sine function, the dynamic and quasi-static behaviors of the stresses are examined by numerical calculations.     

刚性圆形压头作用下热电材料半平面的无摩擦接触问题

热电材料可以将热能转化为电能,反之亦然,这一优良的性质将有助于研发更具成本效益的设备和器件。本文研究了刚性圆形压头作用在热电材料半平面的无摩擦接触问题。假定压头为电导体、热导体,且压头压入深度及与材料的接触区域宽度未知。首先求解电场和温度场,利用傅里叶变换得到了电势函数、温度、电流密度和能量通量的解析表达式。然后求解弹性场,利用积分变换和边界条件,将该热弹性接触问题转化为第一类奇异积分方程并数值求解。数值结果讨论了压头半径和热电载荷对法向接触应力、电流强度因子和能量通量强度因子的影响。结果表明,对于圆压头,热电材料的法向电流密度、法向能量通量在接触边缘表现出奇异性,而表面法向接触应力在接触边缘为零。本文建立的研究模型有助于更深层次的了解热电材料的接触行为。     

Dynamic fracture analysis of an annular interfacial crack between dissimilar magnetoelectroelastic layers

Transient response of an annular interfacial crack between dissimilar magnetoelectroelastic layers under impacts is investigated. On the crack surface, magnetoelectrically impermeable boundary condition is adopted. Using Laplace and Hankel transform techniques, the mixed boundary value problem is reduced to a system of singular integral equations. The integral equations are further reduced to a system of algebraic equations with the aid of Jacobi polynomials. The dynamic field intensity factor and dynamic energy release rate are determined. Numerical results reveal the effects of electric or magnetic loadings and material parameters of composite on crack propagation and growth.     

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

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