An efficient iteration approach for nonlinear boundary value problems in 2D piezoelectric semiconductors

Based on initial nonlinear constitutive equations, we establish the extended displacement and traction boundary integral equations for a piezoelectric medium with a volume electric charge, along with electron and electric current density boundary integral equations for a conductor with a volume electric current. Then, an iterative approach is proposed for investigation of boundary value problems in two-dimensional piezoelectric semiconductors (PSCs). Compared with extended displacements obtained by finite element analysis, this approach is validated via a rectangular PSC under extended external loads. Furthermore, as a numerical example, extended displacements across an elliptical hole in a rectangular PSC are investigated. It is shown that there is a stress concentration near the elliptical hole, which is closely dependent on its shape.     

Dynamic analysis of the different types of elliptic cylindrical inclusions subjected to plane SH-wave scattering

The complex boundary of the elliptical inclusion rendered it difficult to solve the problem of wave scattering. In this study, the steady-state response was analyzed using the wave function expansion method. Subsequently, the Ricker wavelet was employed as the transient disturbance, and Fourier transform was used to determine the distribution of transient dynamic stress concentration around the elliptical inclusion. The effects of wave number, elliptical axial ratio, and difference in material properties on the distribution of the dynamic stress concentration around the elliptical inclusion were evaluated. The numerical results revealed that the dynamic stress concentration always appeared at both ends of the major axis and minor axis of the elliptical inclusion, and the difference in material properties between the inclusion and medium influenced the variations in the dynamic stress concentration factor with the wave number and elliptical axial ratio.     

含椭圆形夹杂的压电体平面应变问题

利用复变函数理论,对在无限远处均匀应力和电位移载荷作用下的含有椭圆形弹性夹杂的横观各向同性压电材料,作了力电分析．在该文有限元结果和前人相关理论解的基础上,提了出一个可接受的认为弹性夹杂体内的应力场为常应力场的假设．在采用了不导通电边界条件之后,获得了以复势形式表示的压电基体的和弹性夹杂体内部的应力场解．     

Dynamic analysis of anisotropic half space containing an elliptical inclusion under SH waves

Based on the methods of complex function, conformal mapping, and multipolar coordinate system, dynamic response of an elliptical inclusion embedded in an anisotropic half space is investigated. In order to find the solution of SH waves, the governing equation is transferred into its normalized form. Then, the scattering wave induced by the inclusion and the standing wave in the inclusion is deduced. Different incident wave angles and the corresponding anisotropy of the half space are considered to obtain the reflected waves. The elliptical inclusion is transferred into a unit circle by conformal mapping method, and then the undetermined coefficients in scattering wave and standing wave are solved by using the continuous condition at the boundary of the inclusion. Subsequently, the dynamic stress concentration factor (DSCF) around the inclusion is calculated and analyzed. Numerical results demonstrate that the distribution of the DSCF is mainly influenced by the incident wave angle and the incident wave number. It is affected by anisotropic parameters as well.     

Three-phase inclusions of arbitrary shape with internal uniform hydrostatic thermal stresses

We investigate the internal thermal stress field of a three-phase inclusion of arbitrary shape which is bonded to an infinite matrix through an interphase layer. The three phases have different thermoelastic constants. It is found that the internal thermal stress field induced by a uniform change in temperature can be uniform and hydrostatic within an inclusion of elliptical or hypotrochoidal shape when the thickness of the interphase layer is properly designed for given material parameters of the three-phase composite. Several examples are presented to demonstrate the solution. The thermal stress analysis of a (Q + 2)-phase inclusion of arbitrary shape with Q ≥ 2 is also carried out under the assumption that all the phases except the internal inclusion share the same elastic constants. It is found that the irregular inclusion shape permitting internal uniform hydrostatic thermal stresses becomes really arbitrary if a sufficiently large number of interphase layers are added between the inclusion and the matrix.     

Elastic fields in two imperfectly bonded half-planes with a thermal inclusion of arbitrary shape

A general method is presented for the rigorous solution of Eshelby’s problem concerned with an arbitrary shaped inclusion embedded within one of two dissimilar elastic half-planes in plane elasticity. The bonding between the half-planes is considered to be imperfect with the assumption that the interface imperfections are uniform. Using analytic continuation, the basic boundary value problem is reduced to a set of two coupled nonhomogeneous first-order differential equations for two analytic functions defined in the lower half-plane which is free of the thermal inclusion. Using diagonalization, the two coupled differential equations are decoupled into two independent nonhomogeneous first-order differential equations for two newly defined analytic functions. The resulting closed-form solutions are given in terms of the constant imperfect interface parameters and the auxiliary function constructed from the conformal mapping which maps the exterior of the inclusion onto the exterior of the unit circle. The method is illustrated using several examples of an imperfect interface. In particular, when the same degree of imperfection is realized in both the normal and tangential directions between the two half-planes, a thermal inclusion of arbitrary shape in the upper half-plane does not cause any mean stress to develop in the lower half-plane. Alternatively, when the imperfect interface parameters are not equal, then a nonzero mean stress will be induced in the lower half-plane by the thermal inclusion of arbitrary shape in the upper half-plane. Detailed results are presented for the mean stress and the interfacial normal and shear stresses caused by a circular and elliptical thermal inclusion, respectively. Results from these calculations reveal that the imperfect bonding condition has a significant effect on the internal stress field induced within the inclusion as well as on the interfacial normal and shear stresses existing between the two half-planes especially when the inclusion is near the imperfect interface.     

含椭圆形夹杂的压电材料平面问题

应用复变函数的Ｆａｂｅｒ级数展开方法，本文研究了含椭圆形夹杂的压电材料平面问题，给出了问题的封闭解·解答表明，椭圆夹杂内的应力、应变、电场强度和电位移均为常量·通过算例，还讨论了正、逆压电效应在基体孔周处的机电行为·     

轴对称金属模具电磁热裂纹止裂中热应力场的分析

为求解金属模具脉冲放电止裂瞬间裂纹尖端附近的热应力场,选择具有半埋藏环形裂纹的金属凹模为研究对象,采用复变函数方法求解了凹模内外环面均匀通入强脉冲电流放电止裂时的热应力场．理论分析结果证实:由于放电瞬间脉冲电流的绕流集中效应,使金属凹模内部环形裂纹尖端附近金属迅速升温,金属熔化形成堆焊,并由于瞬间温升产生热压应力场．研究结果表明:应用电磁热效应止裂技术可以减小裂纹尖端的应力集中,形成的热压应力场有效地阻止金属模具中干线裂纹源的开裂趋势,达到了裂纹止裂目的．     

SH波对界面含有半圆形脱胶的圆柱形弹性夹杂的散射的近似分析

采用Green函数法、复变函数法研究了SH波对界面附近含有半圆形脱胶的圆柱形弹性夹杂的散射,并给出了动应力集中系数的数值结果．首先,界面将整个空间分成上下两部分．在下半空间,给出在含有半圆形凸起的圆柱形弹性夹杂的弹性半空间中,水平表面上任意一点承受时间谐和的出平面线源荷载作用时的位移函数．其次,取该位移函数作为Green函数．上下空间连接时在界面处满足连续性条件,构造出半圆形脱胶裂纹,进而求出应力和位移的表达式．最后作为算例,给出了动应力集中系数的数值结果,分析了介质参数和入射波参数对动应力集中的影响情况．     

On the uniform stress state inside an inclusion of arbitrary shape in a three-phase composite

The stress field inside a two-dimensional arbitrary-shape elastic inclusion bonded through an interphase layer to an infinite elastic matrix subjected to uniform stresses at infinity is analytically studied using the complex variable method in elasticity. Both in-plane and anti-plane shear loading cases are considered. It is shown that the stress field within the inclusion can be uniform and hydrostatic under remote constant in-plane stresses and can be uniform under remote constant anti-plane shear stresses. Both of these uniform stress states can be achieved when the shape of the inclusion, the elastic properties of each phase, and the thickness of the interphase layer are properly designed. Possible non-elliptical shapes of inclusions with uniform hydrostatic stresses induced by in-plane loading are identified and divided into three groups. For each group, two conditions that ensure a uniform hydrostatic stress state are obtained. One condition relates the thickness of the interphase layer to elastic properties of the composite phases, while the other links the remote stresses to geometrical and material parameters of the three-phase composite. Similar conditions are analytically obtained for enabling a uniform stress state inside an arbitrary-shape inclusion in a three-phase composite loaded by remote uniform anti-plane shear stresses.     

Detecting inclusions with a generalized impedance condition from electrostatic data via sampling

In this paper, we derive a sampling method to solve the inverse shape problem of recovering an inclusion with a generalized impedance condition from electrostatic Cauchy data. The generalized impedance condition is a second order differential operator applied to the boundary of the inclusion. We assume that the Dirichlet‐to‐Neumann mapping is given from measuring the current on the outer boundary from an imposed voltage. A simple numerical example is given to show the effectiveness of the proposed inversion method for recovering the inclusion. We also consider the inverse impedance problem of determining the impedance parameters for a known material from the Dirichlet‐to‐Neumann mapping assuming the inclusion has been reconstructed where uniqueness for the reconstruction of the coefficients is proven.     

A lattice Boltzmann model for Maxwell’s equations

In this paper, a lattice Boltzmann model for the Maxwell’s equations is proposed by taking separate sets of distribution functions for the electric and magnetic fields, and a lattice Boltzmann model for the Maxwell vorticity equations with third order truncation error is proposed by using the higher-order moment method. At the same time, the expressions of the equilibrium distribution function and the stability conditions for this model are given. As numerical examples, some classical electromagnetic phenomena, such as the electric and magnetic fields around a line current source, the electric field and equipotential lines around an electrostatic dipole, the electric and magnetic fields around oscillating dipoles are given. These numerical results agree well with classical ones.     

Sharp estimate of electric field from a conductive rod and application

We are concerned with the quantitative study of the electric field perturbation due to the presence of an inhomogeneous conductive rod embedded in a homogenous conductivity. We sharply quantify the dependence of the perturbed electric field on the geometry of the conductive rod. In particular, we accurately characterize the localization of the gradient field (i.e., the electric current) near the boundary of the rod where the curvature is sufficiently large. We develop layer‐potential techniques in deriving the quantitative estimates and the major difficulty comes from the anisotropic geometry of the rod. The result complements and sharpens several existing studies in the literature. It also generates an interesting application in EIT (electrical impedance tomography) in determining the conductive rod by a single measurement, which is also known as the Calderón's inverse inclusion problem in the literature.     

The formulation of linearized boundary integral equations of the anisotropic theory of elasticity and their application in geometrical inverse problems

An elastic bounded anisotropic solid with an elastic inclusion is considered. An oscillating source acts on part of the boundary of the solid and excites oscillations in it. Zero displacements are specified on the other part of the solid and zero forces on the remaining part. A variation in the shape of the surface of the solid and of the inclusion of continuous curvature is introduced and the problem of the theory of elasticity with respect to this variation is linearized. An algorithm for constructing integral representations for such linearized problems is described. The limiting properties of the linearized operators are investigated and special boundary integral equations of the anisotropic theory of elasticity are formulated, which relate the variations of the boundary strain and stress fields with the variations in the shape of the boundary surface. Examples are given of applications of these equations in geometrical inverse problems in which it is required to establish the unknown part of the body boundary or the shape of an elastic inclusion on the basis of information on the wave field on the part of the body surface accessible for observation.     

两压电介质之间的界面夹杂问题

应用Stroh理论,研究了两压电介质之间的刚性介电线夹杂问题。首先该问题被化为Hilbert问题,然后分别给出了压电介质内的复势函数解、夹杂内的电场解和夹杂尖端场的解析表达式。结果表明,在夹杂尖端附近,所有的场变量均呈现奇异性和振荡性,且其强度取决于介质的材料常数和无限场远处的应变。此外,结果还表明,当从夹杂内部趋近夹杂尖端时,夹杂内的电场也呈现奇异性和振荡性。     

An inclusion in one of two joined isotropic elastic half-spaces

Two dissimilar, homogeneous and istropic, elastic half-spacesare bonded together over thier infinite plane of contract. Anarbitrarily shaped finite part of one of them (an inclusion)tends spontaneously to undergo a unifrom infinitesimal strain,but, as it remains attached to and restrained by the surroundingmaterial, an equilibrated state of stress and strain is establishedeverywhere instead. By adopting a convenient expression forthe fundamental field of a point force, we transformed inclusion.For a general shape of the inclussion and for particular sphericaland finite cylindrical shapes in detail, we consider the evaluationof the elastic strain energy, especially of the interactionterm which depends on the location of the inclusion and bothpairs of elastic moduli, and which is of great significancein physical applications.     

应力偶对孔洞附近应力集中的影响

将求解无限弹性平面中孔洞附近应力集中问题的复变函数方法,推广到微极弹性介质的应力集中问题上去,在复平面上给出了二维微极弹性理论应力集中问题的一般解,它可由解析函数与“域函数”构造出来,并利用保角映射的方法来满足非圆孔洞的边界条件。在此基础上建立了求解微极弹性理论中应力集中问题的一般求解方法。最后,对圆形孔洞附近的应力集中系数作了数值计算,并给出了具体结果。     

含正三角形孔口热电材料的断裂力学分析

采用复变函数方法研究了含正三角形孔口热电材料中受到无穷远处均匀电流密度和能量流作用下的断裂问题.在电绝缘和热绝缘边界条件下,获得了温度场和应力场的解析表达式,分析了三角形尺寸、能量流载荷对热电材料性能的影响.结果表明,能量流载荷和三角形尺寸的变化对环向能量流、环向应力、环向热流有较明显的影响.     

Special finite elements including stress concentration effects of a hole

Special finite elements are developed for efficient evaluation of stress concentration around a hole in complex structures. The complex variable formulation is used to derive a special set of stress functions which embody the stress concentration effects of a hole. The stress functions in combination with an independent displacement field assumed along the element boundary are used to construct the special elements with the hybrid displacement finite element method. Several numerical examples are presented to show that the used of special finite elements to model critical regions around a hole, together with conventional finite elements to model other regions away from the hole, is not only very convenient but also highly accurate.     

Basic solutions of a 3D rectangular limited-permeable crack or two 3D rectangular limited-permeable cracks in the piezoelectric/piezomagnetic composite materials

The solutions of a three-dimensional rectangular limited-permeable crack or two three-dimensional rectangular limited-permeable cracks in the piezoelectric/piezomagnetic composite materials were investigated by using the generalized Almansi’s theorem and the Schmidt method. Finally, the relations among the electric field, the magnetic flux field and the stress field near the crack tips were obtained and the effects of the electric permittivity, the magnetic permeability of the air inside the crack, the shape of the rectangular crack on the stress, the electric displacement and magnetic flux intensity factors in the piezoelectric/piezomagnetic composite materials were analyzed.