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

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

半平面内多个纳米圆形夹杂的表/界面效应

基于Gurtin-Murdoch表/界面理论,采用边界积分法,讨论了各向同性的弹性半平面中含有任意多个纳米圆形夹杂问题,得到了受表面/界面影响的纳米复合结构的应力和位移的数值解.最后,给出了半平面中含有单个纳米孔洞和纳米夹杂的数值算例,分析了纳米界面存在对整个半平面结构应力场的影响.     

变边界粘弹性轴对称问题的复变函数法

对边界几何形状、位置随时间变化的变边界结构,给出了用复变函数求解粘弹问题的解析方法。文中用拉普拉斯变换结合平面弹性复变方法,对内外边界变化时粘弹性轴对称问题进行求解。引入两个与时间、空间相关的解析函数,给出了变边界情况下应力、位移以及边界条件与解析函数的关系。当解析函数形式部分确定,则可用边界条件求解其中与时间相关的待定函数。求解待定函数的方程一般情况下为一系列积分方程,特殊情况可求得解析解。对轴对称问题中应力边值问题、位移边值问题以及混合边值问题,分别利用边界条件求得相关系数,从而得到了应力与位移的解析表达。当取Boltzmann粘弹模型时,进行不同边值问题的分析。分析显示,应力、位移的形态与大小均与边界变化过程相关,与固定边界粘弹性问题有较大不同。本文解答可用于粘弹性轴对称问题内外边界任意变化及各种边值问题的力学分析。此外,该法可进一步进行荷载非对称、复杂孔型变边界问题的求解。     

直角平面内弹性圆夹杂对入射平面SH波的散射

利用复变函数法、多极坐标移动技术及傅立叶级数展开求解二维直角平面内圆形弹性夹杂对稳态入射平面SH波的散射问题。首先写出直角平面内不含夹杂时的入射波场和反射波场;其次建立直角平面内含夹杂时夹杂外的散射波解和夹杂内的驻波解,并利用叠加原理写出问题的总波场,借助夹杂边界处应力和位移的连续条件建立求解散射波解和驻波解中未知系数的无穷代数方程组并求解,通过算例具体讨论了直角平面水平边界点的位移幅度比和夹杂边界处径向应力集中系数随不同无量纲波数、入射角及圆孔位置的变化情况,结果表明了算法的有效实用性。     

基于边界元法的非均质薄板弯曲问题的解

本文用边界元法研究非均质无限域弹性薄板弯曲问题.在数值实施过程中,对于夹杂和基体分别形成边界积分方程.通过离散边界积分方程,得到相应的方程组,然后结合界面条件,最终获得问题的求解方程组.在界面的相关量求得之后,可以根据需要来求解基体和夹杂中的有关位置的弯矩.数值结果与已有的解做了对比.     

横观各向同性介质中椭圆夹杂受非弹性剪切变形引起的弹性场的确定

本文求解了横观各向同性介质中椭圆夹杂内受非弹性剪切变形引起的弹性场。采用各向异性弹性力学平面问题的复变函数解法,结合保角变换,获得夹杂内应变能和基体内边界的应力分布和相应的应变能的表达式。进一步,根据最小应变能原理,获得表征夹杂平衡边界的两个特征剪切应变,从而得到了弹性场的解析解。通过应力转换关系,验证了应力解满足夹杂边界上法向正应力和剪应力的连续条件,表明了该解的正确性。本文解可用于复合材料断裂强度的分析中。     

直角平面区域内固定圆形刚性夹杂问题的Green函数解

利用复变函数法、多极坐标移动技术研究了直角平面区域内含有固定圆形夹杂时的反平面问题Green函数解.首先构造出不含夹杂的完整直角平面区域内满足边界应力条件的入射位移场;其次,建立直角平面区域内固定圆形夹杂对该入射场产生的满足直角边界应力自由条件的散射波解,并由叠加原理得到介质内的总波场.最后利用夹杂边界处的位移条件确定出散射波解中的未知系数,最终得到问题的Green函数解,还通过算例讨论了夹杂边界处的径向应力和环向应力随不同波数、角度和不同夹杂位置及不同点源位置的变化情况.算例结果表明了该文方法的有效实用性.     

弹性力学位移构造解待定参数的优化算法研究

从弹性力学平面问题位移解析构造通解的基本原理出发,针对含未知参量的位移函数确定问题,分析了应力边界、位移边界、混合边界的离散节点需要满足的函数关系,构建了以位移解析构造解中未知参量为设计变量,以边界离散节点满足的代数关系为目标函数的优化问题,提出了获得任意边界平面问题的位移构造解中未知参量的优化求解算法,编制了任意节点边界条件的未知参量通用求解程序,给定误差计算的判定方法。求解了平面应力问题的具体实例,通过本文算法与有限元计算结果的误差对比,表明所研究算法的正确性,为任意边界的复杂工程问题求解提供依据。     

弹性椭圆夹杂纵向剪切问题

获得纵向剪切下弹性椭圆夹杂问题的精确解。将复变函数的分区全纯函数理论，Cauchy型积分和Riemann边值问题相结合，求得各复势函数之间的解析关系，从而得到问题的封闭形式解，并给出了界面应力的解析表达式。本文解答与已有文献结果一致。本文发展的分析方法，为求解复杂多连通域的平面弹性问题提供了一条有效途径。     

改进型扩展比例边界有限元法

结合了扩展有限元法(extended finite elementmethods,XFEM)和比例边界有限元法(scaled boundary finite elementmethods,SBFEM)的主要优点,提出了一种改进型扩展比例边界有限元法(improvedextended scaled boundary finite elementmethods,$i$XSBFEM),为断裂问题模拟提供了一条新的途径.类似XFEM,采用两个正交的水平集函数表征材料内部裂纹面,并基于水平集函数判断单元切割类型;将被裂纹切割的单元作为SBFE的子域处理,采用SBFEM求解单元刚度矩阵,从而避免了XFEM中求解不连续单元刚度矩阵需要进一步进行单元子划分的缺陷;同时,借助XFEM的主要思想,将裂纹与单元边界交点的真实位移作为单元结点的附加自由度考虑,赋予了单元结点附加自由度明确的物理意义,可以直接根据位移求解结果得出裂纹与单元边界交点的位移;对于含有裂尖的单元,选取围绕裂尖单元一圈的若干层单元作为超级单元,并将此超级单元作为SBFE的一个子域求解刚度矩阵,超级单元内部的结点位移可通过SBFE的位移模式求解得到,应力强度因子可基于裂尖处的奇异位移(应力)直接获得,无需借助其他的数值方法.最后,通过若干数值算例验证了建议的$i$XSBFEM的有效性,相比于常规XFEM,$i$XSBFEM的基于位移范数的相对误差收敛性较好;采用$i$XSBFEM通过应力法和位移法直接计算得到的裂尖应力强度因子均与解析解吻合\较好.      

The effect of an elastic triangular inclusion on a crack

IntroductionWiththedevelopmentofparticleandfiberreinforcedcomposites,theinclusion_crackinteractionproblemisbecominganimportantfieldbeingstudied .Andasamodel,itisalsousedtostudytheeffectsofmaterialdefectsonthestrengthandfractureofengineeringstructure.TheinterationbetweencircularinclusionandcrackwasstudiedinRefs.[1 -6 ] ;InRefs.[7-1 2 ] ,theinterationbetweenlineinclusionandcrackswasdiscussed ;TheinterationbetweenellipticalinclusionandcrackwasstudiedinRefs.[1 3,1 4] .However,withthedevelopmento…     

波浪与外圆弧开孔壁双圆筒柱的相互作用

应用比例边界有限元法(SBFEM)研究了短峰波与圆筒外接圆弧开孔结构物的相互作用. 求解时将外接圆弧延伸构建一个虚拟圆, 该圆的孔隙影响系数可由矩阵G_0统一进行表达. 整个流场可划分成一个有限域和一个无限域. SBFEM只需对虚拟圆边界进行离散, 使空间维数降低一阶, 在圆的半径方向保持解析, 并且无限域处的辐射边界条件能自动满足. 利用变分原理推导SBFEM方程, 有限域和无限域分别采用贝塞尔函数和汉克尔函数作为基函数来求得对应域的解. 将计算结果与解析解和其他数值方法进行了比较, 验证了该方法是一种用很少单元便能得到精确结果的高效算法. 进一步研究了诸如短峰波波向、结构的几何、材料参数等因素对结构所受波浪载荷及绕射波轮廓的影响, 并进行了分析.      

Stress-strain state of an anisotropic plate with an elliptic hole and thin rigid inclusions

The stress-strain state of an anisotropic plate containing an elliptic hole and thin, absolutely rigid, curvilinear inclusions is studied. General integral representations of the solution of the problem are constructed that satisfy automatically the boundary conditions on the elliptic-hole contour and at infinity. The unknown density functions appearing in the potential representations of the solution are determined from the boundary conditions at the rigid inclusion contours. The problem is reduced to a system of singular integral equations which is solved by a numerical method. The effects of the material anisotropy, the degree of ellipticity of the elliptic hole, and the geometry of the rigid inclusions on the stress concentration in the plate are studied. The numerical results obtained are compared with existing analytical solutions. __________ Translated from Prikladnaya Mekhanika i Tekhnicheskaya Fizika, Vol. 48, No. 4, pp. 173–180, July–August, 2007.     

Boundary element analysis of interaction between an elastic rectangular inclusion and a crack

The interaction between an elastic rectangular inclusion and a kinked crack inan infinite elastic body was considered by using boundary element method. The new complexboundary integral equations were derived. By introducing a complex unknown function H(t)related to the interface displacement density and traction and applying integration by parts,the traction continuous condition was satisfied automatically. Only one complex boundaryintegral equation was obtained on interface and involves only singularity of order l/ r. Toverify the validity and effectiveness of the present boundary element method, some typicalexamples were calculated. The obtained results show that the crack stress intensity factorsdecrease as the shear modulus of inclusion increases. Thus, the crack propagation is easiernear a softer inclusion and the harder inclusion is helpful for crack arrest.     

An analytical solution of an axisymmetric problem of deforming an isotropic half-space with an elastic fixed boundary under the action of a distributed load

An analytical solution to the axisymmetric problem on the action of a distributed load on an isotropic half-space when the load is given by a function dependent on the radial coordinate is obtained. The surface of the half-space is elastically fixed outside the circular domain of load application, the shear stresses are absent along the entire boundary, and the stresses vanish at infinity. At the boundary and inside the elastic half-space, the solutions are represented by the formulas for the stress tensor components and for the displacement vector components.     

Computational model for short-fiber composites with eigenstrain formulation of boundary integral equations

A computational model is proposed for short-fiber reinforced materials with the eigenstrain formulation of the boundary integral equations(BIE)and solved with the newly developed boundary point method(BPM).The model is closely derived from the concept of the equivalent inclusion of Eshelby tensors.Eigenstrains are iteratively determined for each short.fiber embedded in the matrix with various properties via the Eshelby tensors,which can be readily obtained beforehand either through analytical or numerical means.As unknown variables appear only on the boundary of the solution domain,the solution scale of the inhomogeneity problem with the model is greatly reduced.This feature is considered significant because such a traditionally time-consuming problem with inhomogeneity can be solved most cost-effectively compared with existing numerical models of the FEM or the BEM.The numerical examples are presented to compute the overall elastic properties for various short-fiber reinforced composites over a representative volume element(RVE),showing the validity and the effectiveness of the proposed computational modal and the solution procedure.     

Computational model for short-fiber composites with eigenstrain formulation of boundary integral equations

A computational model is proposed for short-fiber reinforced materials with the eigenstrain formulation of the boundary integral equations （BIE） and solved with the newly developed boundary point method （BPM）. The model is closely derived from the concept of the equivalent inclusion Of Eshelby tensors. Eigenstrains are iteratively determined for each short-fiber embedded in the matrix with various properties via the Eshelby tensors, which can be readily obtained beforehand either through analytical or numerical means. As unknown variables appear only on the boundary of the solution domain, the solution scale of the inhomogeneity problem with the model is greatly reduced. This feature is considered significant because such a traditionally time-consuming problem with inhomogeneity can be solved most cost-effectively compared with existing numerical models of the FEM or the BEM. The numerical examples are presented to compute the overall elastic properties for various short-fiber reinforced composites over a representative volume element （RVE）, showing the validity and the effectiveness of the proposed computational modal and the solution procedure.     

双周期圆柱形夹杂纵向剪切问题的精确解

研究无限介质中矩形排列双周期圆柱形夹杂的纵向剪切问题．利用Eshelby等效夹杂理论并结合双周期与双准周期解析函数工具，为这类考虑夹杂相互影响的问题提供了一个严格又实用的分析方法，求得了问题的全场级数解．作为退化情形得到单夹杂问题的经典解答，双周期孔洞、双周期刚性夹杂及单行(列)周期弹性夹杂等问题也可作为特殊情况被解决．数值结果揭示了这类非均匀材料力学性质随微结构参数变化的规律．     

The null field approach to two-dimensional elastostatics

The vector basis functions, necessary for solving two-dimensional inclusion problems in an elastic solid under time independent conditions by means of the null field approach (T-matrix method), are obtained as a zero frequency limit of the corresponding basis functions commonly used in elastodynamics. The expansion of the fields appearing in the surface integral representation of the static displacement can thus be achieved, leading to the T-matrix equations of 2d-elastostatics. We specialize the problem to the simple boundary condition case of a single cavity and develop the analytical expressions as much as possible before numerical implementation. A numerical test for the ellipse and some examples for the superellipse, with applied static pressure or shear stress at infinity, are given.     

Axisymmetric contact problem for an elastic half-space and a circular cover plate

The contact interaction problem for a thin circular rigid cover plate and an elastic half-space loaded at infinity by a tensile force directed in parallel to the boundary of the half-space is considered. It is assumed that the cover plate is not resistant to bending deformations. The problem can be reduced to an integral equation of the first kind whose kernel has a logarithmic singularity. The equation is solved approximately by the Multhopp-Kalandia method. The resulting approximate solution is compared with the previously obtained asymptotic solution.