Fracture analysis of magnetoelectroelastic bimaterials with imperfect interfaces by symplectic expansion

A Hamiltonian-based analytical method is used to study the mode Ⅲ interface cracks in magnetoelectroelastic bimaterials with an imperfect interface. By introducing an unknown vector, the governing equations are reformulated in sets of first-order ordinary differential equations. Using separation of variables, eigensolutions in the symplectic space are obtained. An exact solution of the unknown vector is obtained and expressed in terms of symplectic eigensolutions. Singularities of mechanical, electric, and magnetic fields are evaluated with the generalized intensity factors. Comparisons are made to verify accuracy and stability of the proposed method. Numerical examples including mixed boundary conditions are given.     

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.     

Symplectic solution system for reissner plate bending

Based on the Hellinger-Reissner variatonal principle for Reissner plate bendingand introducing dual variables, Hamiltonian dual equations for Reissner plate bending werepresented. Therefore Hamiltonian solution system can also be applied to Reissner platebending problem, and the transformation from Euclidian space to symplectic space and fromLagrangian system to Hamiltonian system was realized. So in the symplectic space whichconsists of the original variables and their dual variables, the problem can be solved viaeffective mathematical physics methods such as the method of separation of variables andeigenfunction-vector expansion. All the eigensolutions and Jordan canonical formeigensolutions for zero eigenvalue of the Hamiltonian operator matrix are solved in detail, and their physical meanings are showed clearly. The adjoint symplectic orthonormal relation of the eigenfunction vectors for zero eigenvalue are formed. It is showed that the alleigensolutions for zero eigenvalue are basic solutions of the Saint-Venant problem and theyform a perfect symplectic subspace for zero eigenvalue. And the eigensolutions for nonzeroeigenvalue are covered by the Saint-Venant theorem. The symplectic solution method is notthe same as the classical semi-inverse method and breaks through the limit of the traditional semi-inverse solution. The symplectic solution method will have vast application.     

Penny-shaped interfacial crack between dissimilar magnetoelectroelastic layers

A penny-shaped interfacial crack between dissimilar magnetoelectroelastic layers subjected to magnetoelectromechanical loads is investigated,where the magnetoelectrically impermeable crack surface condition is adopted. By using Hankel transform technique,the mixed boundary value problem is firstly reduced to a system of singular integral equations,which are further reduced to a system of algebraic equations. The field intensity factors and energy release rate are finally derived. Numerical results elucidate the eects of crack configuration,electric and/or magnetic loads,and material parameters of the magnetoelectroelastic layers on crack propagation and growth. This work should be useful for the design of magnetoelectroelastic composite structures.     

Paradox solution on elastic wedge dissimilar materials

According to the Hellinger-Reissner variational principle and introducing proper transformation of variables , the problem on elastic wedge dissimilar materials can be led to Hamiltonian system, so the solution of the problem can be got by employing the separation of variables method and symplectic eigenfunction expansion under symplectic space, which consists of original variables and their dual variables . The eigenvalue - 1 is a special one of all symplectic eigenvalue for Hamiltonian system in polar coordinate . In general, the eigenvalue - 1 is a single eigenvalue, and the classical solution of an elastic wedge dissimilar materials subjected to a unit concentrated couple at the vertex is got directly by solving the eigenfunction vector for eigenvalue - 1. But the eigenvalue - 1 becomes a double eigenvalue when the vertex angles and modulus of the materials satisfy certain definite relationships and the classical solution for the stress distribution becomes infinite at this moment, that is, the para     

ELECTRO-MAGNETO-THERMOELASTIC PLANE WAVES IN MICROPOLAR SOLID INVOLVING TWO TEMPERATURES

<正>The model of equations of micropolar generalized magneto-thermoelasticity is introduced within the context of the theory of two temperatures generalized thermoelasticity and we consider a problem of an isotropic homogeneous micropolar medium taking into account the heat effects and allowing the magnetic field effects.A plane wave analysis is employed to obtain the exact formulas of the two temperatures(conductive and mechanical),displacement components, micro-rotation components,stresses,couple stresses,induced electric current,electric field and magnetic field.Arbitrary application is chosen to enable us to get the complete solution.The considered variables are presented graphically and discussions are made for the results.     

PLANE STRAIN PROBLEM OF PIEZOELECTRIC SOLID WITH ELLIPTIC INCLUSION

By using the complex variables function theory, a plane strain electro-elastic analysis was performed on a transversely isotropic piezoelectric material containing an elliptic elastic inclusion, which is subjected to a uniform stress field and a uniform electric displacement loads at infinity. Based on the present finite element results and some related theoretical solutions, an acceptable conjecture was found that the stress field is constant inside the elastic inclusion. The stress field solutions in the piezoelectric matrix and the elastic inclusion were obtained in the form of complex potentials based on the impermeable electric boundary conditions.     

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.     

Basic solutions of multiple parallel symmetric mode-III cracks in functionally graded piezoelectric/piezomagnetic material plane

The Schmidt method is adopted to investigate the fracture problem of multiple parallel symmetric and permeable finite length mode-III cracks in a functionally graded piezoelectric/piezomagnetic material plane. This problem is formulated into dual integral equations, in which the unknown variables are the displacement jumps across the crack surfaces. In order to obtain the dual integral equations, the displacement jumps across the crack surfaces are directly expanded as a series of Jacobi polynomials. The results show that the stress, the electric displacement, and the magnetic flux intensity factors of cracks depend on the crack length, the functionally graded parameter, and the distance among the multiple parallel cracks. The crack shielding effect is also obviously presented in a functionally graded piezoelectric/piezomagnetic material plane with mul- tiple parallel symmetric mode-III cracks.     

FREE VIBRATION OF FUNCTIONALLY GRADED, MAGNETO-ELECTRO-ELASTIC, AND MULTILAYERED PLATES

The state-space method is employed to evaluate the modal parameters of functionally graded, magneto-electro-elastic, and multilayered plates. Based on the assumption that the properties of the functionally graded material are exponential, the state equation of structural vibration which takes the displacement and stress of the structure as state variables is derived. The natural frequencies and modal shapes are calculated based on the general solutions of the state equation and boundary conditions given in this paper. The influence of the functionally graded exponential factor on the elastic displacement, electric, and magnetic fields of the structure are discussed by assuming a sandwich plate model with different stacking sequences.     

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

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

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.     

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

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

电磁弹性动力学初边值问题12类变量广义变分原理

基于Yao建立的电磁弹性固体广义变分原理,运用关于非传统Hamilton型广义变分原理的方法,建立了电磁弹性动力学初边值问题的12类变量广义变分原理,可反映该问题的全部特征,其独立变分变量为该问题的全部变量,即位移、速度、动量、应变、应力、电位移、磁感应强度、电场强度、磁场强度、电标量势、磁标量势和磁矢量势.本文建立的...     

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.     

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.     

Dielectric crack problem for a magnetoelectroelastic strip with functionally graded properties

The paper presents a fracture analysis for an electromagnetically dielectric crack in a functionally graded magnetoelectroelastic strip. It is considered that the material properties are varying exponentially along the width direction. Under the assumption of the in-plane magneto-electro-mechanical loadings, the dielectric crack is simulated by using the semi-permeable crack-face boundary conditions. The Fourier transform technique is applied to solve the boundary-value problem and four coupling singular integral equations are determined. A nonlinear system of algebraic equations is further derived and solved numerically to determine the electromagnetic field inside the crack. Then the field intensity factors of stress, electric displacement, and magnetic induction are given. Through the numerical computations, the effects of the material non-homogeneity and the permeability of crack interior on the electric displacement and the magnetic induction at the crack faces are studied. The variations of the intensity factors of stress, electric displacement, and magnetic induction versus the geometry of the crack, the strip width, and the material non-homogeneity are presented in graphics respectively.     

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.     

Nonlinear fracture of 2D magnetoelectroelastic media: Analytical and numerical solutions

A strip electric–magnetic polarization saturation (SEMPS) model is developed to study the electric and magnetic yielding effects on a crack in magnetoelectroelastic (MEE) media. In this model, the MEE medium is treated as being mechanically brittle, and electrically and magnetically ductile. Analogously to the classic Dugdale model, the electric and magnetic yielding zones in front of the crack are represented for simplicity by two strips. In the electric yielding strip the electric displacement equals the electric displacement saturation and meanwhile in the magnetic yielding zone the magnetic induction equals the magnetic induction saturation. The nonlinear analytical solution of this SEMPS model of crack in an infinite MEE medium is obtained using an integral equation approach. The equivalence between the proposed SEMPS model and the existing strip electric–magnetic breakdown (SEMB) model is demonstrated.To analyze the nonlinear fracture problem in the corresponding finite MEE media, the non-linear hybrid extended displacement discontinuity-fundamental solution (NLHEDD-FS) method is modified, and a multiple iteration approach is adapted to determine the electric and magnetic yielding zones. Comparing with the analytical solution, the applicability and effectiveness of the NLHEDD-FS method is verified. Numerical results based on the SEMPS model for a center crack in infinite and finite MEE strip are presented.     

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

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