BOUNDARY INTEGRAL FORMULA OF ELASTIC PROBLEMS IN CIRCLE PLANE

IntroductionConcerningtheelasticplaneprobleminaunitcircle ,ZhengShenzhouandZhengXueliangdevelopedaboundaryintegralformulaofthestressfunction[1]:Φ(r,θ) =-( 1 -r2 ) 24π ∫2π0ν( φ)1 -2rcos(θ-φ) r2 dφ　　 12π∫2π011 -2rcos(θ-ω) r2 dω∫2π0μ( φ)1 -cos(ω-φ) dφ　　 1 -r22π∫2π0μ( φ)1 -2rcos(θ -φ) r2 dφ　　 ( 0 ≤r <1 ) ,( 1 )whereμ(θ) =Φ(r,θ) ｜r=1,ν(θ) = Φ n r=1= Φ r r=1.Intheformula ( 1 )theseconditemisastrongsingularintegral,itshouldbeunderstoodasanintegra…     

ANALYSIS OF SYMMETRIC LAMINATED RECTANGULAR PLATES IN PLANE STRESS

Symmetric laminated plates used usually are anisotropic plates. Based on the fundamental equation for anisotropic rectangular plates in plane stress problem, a general analytical solution is established accurately by method of stress function. Therefore the general formula of stress and displacement in plane is given. The integral constants in general formula can be determined by boundary conditions. This general solution is composed of solutions made by trigonometric function and hyperbolic function, which can satisfy the problem of arbitrary boundary conditions along four edges, and the algebraic polynomial solutions which can satisfy the problem of boundary conditions at four corners. Consequently this general solution can be used to solve the plane stress problem with arbitrary boundary conditions. For example, a symmetric laminated square plate acted with uniform normal load, tangential load and nonuniform normal load on four edges is calculated and analyzed.     

应用边界积分法求圆形夹杂问题的解析解

边界元方法作为一种数值方法, 在各种科学工程问题中得到了广泛的应用.本文参考了边界元法的求解思路, 从Somigliana等式出发, 利用格林函数性质,得到了一种边界积分法, 使之可以用来寻求弹性问题的解析解.此边界积分法也可以从Betti互易定理得到. 应用此新方法, 求解了圆形夹杂问题.首先设定夹杂与基体之间完美连接, 将界面处的位移与应力按照傅里叶级数展开,根据问题的对称性与三角函数的正交性来简化假设, 减少待定系数的个数.其次选择合适的试函数(试函数满足位移单值条件以及无体力的线弹性力学问题的控制方程),应用边界积分法, 求得界面处的位移与应力的值. 然后再求解域内位移与应力.得到了问题的精确解析解, 当夹杂弹性模量为零或趋向于无穷大时,退化为圆孔或刚性夹杂问题的解析解. 求解过程表明,若问题的求解区域包含无穷远处时, 所取的试函数应满足无穷远处的边界条件.若求解区域包含坐标原点, 试函数在原点处位移与应力应是有限的.结果表明了此方法的有效性.      

A boundary-field equation method for a nonlinear exterior elasticity problem in the plane

This paper presents a new method for solving a nonlinear exterior boundary value problem arising in two-dimensional elasto-plasticity. The procedure is based on the introduction of a sufficiently large circle that divides the exterior domain into a bounded region and an unbounded one. This allows us to consider the Dirichlet-Neumann mapping on the circle, which provides an explicit formula for the stress in terms of the displacement by using an appropriate infinite Fourier series. In this way we can reduce the original problem to an equivalent nonlinear boundary value problem on the bounded domain with a natural boundary condition on the circle. Hence, the resulting weak formulation includes boundary and field terms, which yields the so called boundary-field equation method. Next, we employ the finite Fourier series to obtain a sequence of approximating nonlinear problems from which the actual Galerkin schemes are derived. Finally, we apply some tools from monotone operators to prove existence, uniqueness and approximation results, including Cea type error estimates for the corresponding discrete solutions.     

On multiple circular inclusions in plane magnetoelasticity

A general series solution to the magnetoelastic problem of interacting circular inclusions in plane magnetoelasticity is provided in this paper. By the use of complex variable theory and Laurent series expansion method, the general expression of the magnetic and the magnetoelastic complex potentials for the circular inclusion problem is derived. Expanding the definition of the Airy’s stress function of pure elastic field into the magnetoelastic field and applying the superposition method, the general expression then can be reduced to a set of linear algebraic equations and solved in a series form. An approximate closed form solution for the case of two arbitrarily located inclusions is also provided. For illustrating the effect of the pertinent parameters, the numerical results of the interfacial magnetoelastic stresses are displayed in graphic form.     

On the Green's Functions and Boundary Integral Formulation of Elastic Planes with Cutouts

ABSTRACT

The problem of an infinite elastic plane that contains a hole of arbitrary shape and is subjected to a concentrated unit load is considered. The Green's function (influence function) for the problem is formulated by means of two complex potential functions. This is accomplished by mapping the region that is exterior to the hole onto a unit circle. A class of closed contour hole shapes is analyzed. Green's functions for an elliptical hole and a class of triangular holes are determined. Green's functions for a class of rectangular holes are also discussed. In order to determine stress and displacement fields for the finite plane problem, Green's function is employed and an indirect boundary integral equation is formulated, with the integrand of the integral equation incorporating the effect of the hole. The contour of the hole is no longer considered a part of the boundary and only the contour of the region that is exterior to the hole is subdivided into boundary elements. Examples for elliptical and triangular holes are solved.     

含界面层圆形夹杂中螺旋位错的弹性分析

研究了复合材料含界面层圆形夹杂内部的一个螺旋位错在夹杂，界面层与基体材料中产生的弹性干涉，将复变函数的分区亚纯函数理论与柯西型积分，罗朗级数相结合，求出了各分区复势的解析关系，化为一个关于界面层复势的函数方程，用显示式表达了问题的结果，揭示了界面层参数对位错干涉能与位错力的影响规律。该解析方法较经典级数方法未知量大量减少，表达式更加简洁，结果的特殊情形包含了若干已有成果。     

Interaction of a rigid punch and a circular plate with a fixed side and a stress-free face

The axisymmetric contact problem of a rigid punch indentation into an elastic circular plate with a fixed side and a stress-free face is considered. The problem is solved by a method developed for finite bodies which is based on the properties of a biorthogonal system of vector functions. The problem is reduced to a Volterra integral equation (IE) of the first kind for the contract pressure function and to a system of two Volterra IE of the first kind for functions describing the derivative of the displacement of the plate upper surface outside the punch and the normal (or tangential) stress on the plate lower fixed surface. The last two functions are sought as the sum of a trigonometric series and a power-law function with a root singularity. The obtained ill-conditioned systems of linear algebraic equations are regularized by introducing small parameters and have a stable solution. A method for solving the Volterra IE is given. The contact pressure functions, the normal and tangential stresses on the plate fixed surface, and the dimensionless indentation force are found. Several examples of a plane punch computation are given.     

A Representation Formula for the Velocity Part of 3D Time-Dependent Oseen Flows

We prove a representation formula for the solution of an initial-boundary value problem associated with the 3D time-dependent Oseen system. This formula involves the solution of an integral equation on the lateral boundary of the space-time cylinder. The key point of this article is that the assumptions on the data are chosen in such a way that the formula in question may be used in a theory of the asymptotic behaviour of solutions to the nonstationary 3D Navier–Stokes system with Oseen term (incompressible viscous flow in an exterior domain with nonzero velocity at infinity).     

Dynamic stress intensity factors around two rectangular cracks in an infinite elastic plate under impact load

A dynamic problem for two equal rectangular cracks in an infinite elastic plate is considered. The two cracks are placed perpendicular to the plane surfaces of the plate. An incoming shock tensile stress is returned by the cracks. In the Laplace transform domain, the boundary conditions at the two sides of the plate are satisfied using the Fourier transform technique. The mixed boundary conditions are reduced to dual integral equations. Crack displacement is expanded in a series of functions which are zero outside of the cracks. The unknown coefficients in the series are determined by the Schmidt method. The stress intensity factors are defined in the Laplace transform domain and these are inverted using a numerical method.     

An embedding method for modeling micromechanical behavior and macroscopic properties of composite materials

This paper presents a numerical method for modeling the micromechanical behavior and macroscopic properties of fiber-reinforced composites and perforated materials. The material is modeled by a finite rectangular domain containing multiple circular holes and elastic inclusions. The rectangular domain is assumed to be embedded within a larger circular domain with fictitious boundary loading represented by truncated Fourier series. The analytical solution for the complementary problem of a circular domain containing holes and inclusions is obtained by using a combination of the series expansion technique with a direct boundary integral method. The boundary conditions on the physical external boundary are satisfied by adopting an overspecification technique based on a least squares approximation. All of the integrals arising in the method can be evaluated analytically. As a result, the elastic fields and effective properties can be expressed explicitly in terms of the coefficients in the series expansions. Several numerical experiments are conducted to verify the accuracy and efficiency of the numerical method and to demonstrate its application in determination of the macroscopic properties of composite materials.     

Helmholtz边界积分方程中奇异积分间接求解方法

提出了间接求解传统Helmholtz边界积分方程CBIE的强奇异积分和自由项系数,以及Burton-Miller边界积分方程BMBIE中的超强奇异积分的特解法。对于声场的内域问题,给出了满足Helmholtz控制方程的特解,间接求出了CBIE中的强奇异积分和自由项系数。对于声场外域对应的BMBIE中的超强奇异积分,按Guiggiani方法计算其柯西主值积分需要进行泰勒级数展开的高阶近似,公式繁复,实施困难。本文给出了满足Helmholtz控制方程和Sommerfeld散射条件的特解,提出了间接求出超强奇异积分的方法。推导了轴对称结构外场问题的强奇异积分中的柯西主值积分表达式,并通过轴对称问题算例证明了本文方法的高效性。数值结果表明,对于内域问题,采用本文特解法的计算结果优于直接求解强奇异积分和自由项系数的结果,且本文的特解法可避免针对具体几何信息计算自由项系数,因而具有更好的适用性。对于外域问题,两者精度相当,但本文的特解法可避免对核函数进行高阶泰勒级数展开,更易于数值实施。     

The Problem of an External Circular Crack Under Asymmetric Loadings

ＩｎｔｒｏｄｕｃｔｉｏｎＷｈｅｎｔｈｅｂｏｕｎｄａｒｙｉｎｔｅｇｒａｌｅｑｕａｔｉｏｎｍｅｔｈｏｄｉｓａｐｐｌｉｅｄｔｏｃｒａｃｋａｎａｌｙｓｉｓ,ｏｎｌｙｎｕｍｅｒｉｃａｌｓｏｌｕｔｉｏｎｓｃａｎｂｅｏｂｔａｉｎｅｄ ,ｓｕｃｈａｓ:ｔｈｅｔｙｐｉｃａｌｗｏｒｋｓｏｆＳｎｙｄｅｒａｎｄＣｒｕｓｅ[1],Ｃｒｏｕｃｈ[2 ],Ｂｌａｎｄｆｏｒｄｅｔａｌ.[3],Ｐｏｒｔｅｌａｅｔａｌ.[4 ],Ｂｕｉ[5 ],Ｗｅａｖｅｒ[6 ]ａｎｄＷＡＮＧｅｔａｌ.[7- 9].Ｉｔｉｓｖｅｒｙｄｉｆｆｉｃｕｌｔｔｏａｐｐｌｙｔｈｅｂｏｕｎｄａｒ…     

New Green’s function for stress field and a note of its application in quantum-wire structures

The elastic field caused by the lattice mismatch between the quantum wires and the host matrix can be modeled by a corresponding two-dimensional hydrostatic inclusion subjected to plane strain conditions. The stresses in such a hydrostatic inclusion can be effectively calculated by employing the Green’s functions developed by Downes and Faux, which tend to be more efficient than the conventional method based on the Green’s function for the displacement field. In this study, Downes and Faux’s paper is extended to plane inclusions subjected to arbitrarily distributed eigenstrains: an explicit Green’s function solution, which evaluates the stress field due to the excitation of a point eigenstrain source in an infinite plane directly, is obtained in a closed-form. Here it is demonstrated that both the interior and exterior stress fields to an inclusion of any shape and with arbitrarily distributed eigenstrains are represented in a unified area integral form by employing the derived Green’s functions. In the case of uniform eigenstrain, the formulae may be simplified to contour integrals by Green’s theorem. However, special care is required when Green’s theorem is applied for the interior field. The proposed Green’s function is particularly advantageous in dealing numerically or analytically with the exterior stress field and the non-uniform eigenstrain. Two examples concerning circular inclusions are investigated. A linearly distributed eigenstrain is attempted in the first example, resulting in a linear interior stress field. The second example solves a circular thermal inclusion, where both the interior and exterior stress fields are obtained simultaneously.     

SH波对直角平面区域内圆形衬砌的散射与地震动

研究了以任意方向入射的平面SH 波对直角平面区域内圆形衬砌的散射与地震动问题．利用复变函数法、多极坐标移动技术及镜像法,将直角平面区域内的波场延拓于半无限空间,根据边界条件将该问题的解答可归结为对一组无穷代数方程组的求解问题,通过算例考虑衬砌的动应力集中和地表位移．结果表明,衬砌的动应力集中和地表位移幅值取决于入射波频率、角度、衬砌的位置和厚度比．     

Artificial boundary method for the exterior Stokes flow in three dimensions

In this paper, the artificial boundary method is considered for the numerical simulation of the exterior Stokes flow in three dimensions. First, an exact relation between the normal stress and the velocity field is obtained on a spherical artificial boundary. With the relation specified on the artificial boundary, the original problem is reduced to a new one only defined on a finite domain. After that, an variational problem equivalent to the reduced problem is derived. By truncating the series term in the formulation, a sequence of approximate variational problems are obtained, which can then be solved with a suitable finite‐element scheme. Finally, a numerical example is presented to show the performance of the method. Copyright © 2003 John Wiley & Sons, Ltd.     

Evaluation of strongly singular domain integrals for internal stresses in functionally graded materials analyses using RIBEM

An accurate evaluation of strongly singular domain integral appearing in the stress representation formula is a crucial problem in the stress analysis of functionally graded materials using boundary element method.To solve this problem,a singularity separation technique is presented in the paper to split the singular integral into regular and singular parts by subtracting and adding a singular term.The singular domain integral is transformed into a boundary integral using the radial integration method.Analytical expressions of the radial integrals are obtained for two commonly used shear moduli varying with spatial coordinates.The regular domain integral,after expressing the displacements in terms of the radial basis functions,is also transformed to the boundary using the radial integration method.Finally,a boundary element method without internal cells is established for computing the stresses at internal nodes of the functionally graded materials with varying shear modulus.     

Analysis of an Anisotropic Finite Wedge under Antiplane Deformation

The antiplane deformation of an anisotropic wedge with finite radius is considered in this paper within the classical linear theory of elasticity. The traction-free condition is imposed on the circular segment of the wedge. Three different cases of boundary conditions on the radial edges are considered, which are: traction-displacement, displacement-displacement and traction-traction. The solution to the governing differential equation of the problem is accomplished in the complex plane by relating the displacement field to a complex function. Several complex transformations are defined on this complex function and its first and second derivatives to formulate the problem in each of the three cases of the problem corresponding to the radial boundary conditions, separately. These transformations are then related to integral transforms which are complex analogies to the standard finite Mellin transforms of the first and second kinds. Closed form expressions are obtained for the displacement and stress fields in the entire domain. In all cases, explicit expressions for the strength of singularity are derived. These expressions show the dependence of the order of stress singularity on the wedge angle and material constants. In the displacement-displacement case, depending upon the applied displacement, a new type of stress singularity has been observed at the wedge apex. This revised version was published online in August 2006 with corrections to the Cover Date.     

圆外Stokes问题的边界积分公式

针对圆外区域的Stokes问题,利用Fourier展开法,通过自然边界归化得到了一个只与边界速度有关的Stokes问题的边界积分公式.根据此公式及边界速度值,求得区域内速度及压力分布的解析表达式,并通过数值积分的方法进行求解,计算量小,所得速度及压力分布图曲线光滑.最后借助流体软件进行数值计算,结果验证了边界积分公式的正确性、可行性.     

EQUIVALENT BOUNDARY INTEGRAL EQUATIONS WITH INDIRECT VARIABLES FOR PLANE ELASTICITY PROBLEMS

IntroductionIthasbeenratheralonghistorythattheBoundaryElementMethod (BEM )isappliedtosolvetheplaneelasticityproblems[1～2 ].However,theEBIE ,whichisequivalenttotheoriginalboundaryvalueproblem ,hasnotbeenfullyappreciatedandsolvedinBEMcommunity .TheconventionalboundaryintegralequationswithindirectvariablesarediscussedthoroughlyanditisshownthatthepreviousresultsarenotEBIE ,i.e .,sometimes,thereexistsnosolutionormorethanonesolutiontothem .Themainkeyliesintheexactformoftheexteriorproblems.I…