完备型二次边界轮廓法的研究及求解断裂力学Ja,L,M积分

对线弹性平面问题的边界轮廓法，选用完备的二次位移形函数，使求问题的维数降低两维，给出了求解边界位移和面力以及内点应力的求解方法。证明平面弹怀断鲜明力学Ｊａ积分、Ｍ积分、Ｌ积分方程的被积函数的散度均等于零，将它们分别转化为边界点的位移和面力的线性迭加，无需计算数值积分，算例表明，本文方法具有较高的精度。     

完备型二次边界轮廓法的研究及求解断裂力学J_α、L、M积分

对线弹性平面问题的边界轮廓法,选用完备的二次位移形函数,使求解问题的维数降低两维,给出了求解边界位移和面力以及内点应力的求解方法。证明平面弹性断裂力学J_a积分、M积分、L积分方程的被积函数的散度均等于零,将它们分别转化为边界点的位移和面力的线性迭加,无需计算数值积分。算例表明,本文方法具有较高的精度。     

用面力边界积分方程求解断裂力学J积分

证明面力边界积分方程被积函数的散度等于零，应用Stokes公式，对平面线弹性问题，将面力边界积分的求解转化为边界点的位移势函数的点值计算。应用边界积分方程的求解结果，推导出J积分亦可表示为边界点的积分势函数的点值计算，无需进行数值积分，实例计算说明该方法的有效性。     

完备型二次边界轮廓法的研究及J积分和应力强度因子的计算

选择二次完全多项多作为位移形函数,对边界轮廓法作了进一步的发展,证明二维弹性断裂问题的Ｊ积分方程的被积分函数的散度等于零,将Ｊ积分化为边界点的势函数数值的计算,无需计算数值积分,算例表明,该方法较传统边界元法求得的结果精度更好。     

一种基于直接计算高阶奇异积分的断裂力学双边界积分方程分析法

相对于有限元法,边界单元法在求解断裂问题上有着独特的优势,现有的边界单元法中主要有子区域法和双边界积分方程法.采用一种改进的双边界积分方程法求解二维、三维断裂问题的应力强度因子,对非裂纹边界采用传统的位移边界积分方程,只需对裂纹面中的一面采用面力边界积分方程,并以裂纹间断位移为未知量直接用于计算应力强度因子.采用一种高阶奇异积分的直接法计算面力边界积分方程中的超强奇异积分;对于裂纹尖端单元,提供了三种不同形式的间断位移插值函数,采用两点公式计算应力强度因子.给出了多个具体的算例,与现存的精确解或参考解对比,可得到高精度的计算结果.      

导数场边界积分方程分析近边界应力分布

研究二维弹性力学问题边界积分方程,通过分部积分变换消除了常规导数边界积分方程中的超奇异积分,获得仅含强奇异积分的应力自然边界积分方程.对于近边界应力的计算,进一步运用正则化算法解析计算其中的几乎强奇异积分.较常规边界元法相比,应力自然边界积分方程可以求解离边界更加接近的内点应力值.算例证明了文中方法的可应用性和有效性.     

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

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

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

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

三维Laplace方程边界元中线性单元的精确积分法

边界元中的边界积分计算影响计算精度和计算速度。非奇异积分一般采用数值积分，当配置点接近积分单元时，计算精度降低。未知函数线性插值得到的解是连续解，但计算难度增大。本文采用积分区域变换，将三维Laplace问题的二维积分化为一维积分，这样奇异积分和非奇异积分能采用精确积分的方法计算，使求解精度，计算速度都得到提高。     

边界元法计算近边界点参量的一个通用算法

针对边界元法存在近力界点参量计算的困难,给出了一个通用性方法,将近边界点到边界单元的距离参数通过分部积分变换到积分式之外,从而计算出二维问题近边界点参量的几乎强奇异和超奇异积分,由此,对任何近边界点参量,提出了整套计算方案,算例证明了本法的有效性。     

应力边界元法解平面弹性问题

本文提出了求解平面弹性问题的应力边界元法。简述了边界积分方程的建立,给出了常单元离散化时求系数的解析式。这种方法适用于应力边界值问题。边界积分方程中的一个边界函数就是边界点法向应力和切向应力之和,因此计算孔边应力非常方便。作为数值算例,计算了有孔无限板的孔边应力。应力边界元法也可应用于平面热弹性问题和平板弯曲问题。     

一种新型的求解扭转问题的边界元方法

证明了柱体自由扭转的边界积分方程被积函数的散度等于零,将翘曲函数表示为翘曲势函数在边界点的数值计算,避免求解奇异的数值积分。实例计算表明,该表度较高。     

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

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

边界元法中计算几乎奇异积分的一种无奇异算法

边界元法中存在几乎奇异积分的计算困难。引起边界单元上几乎奇异积分的因素是源点到其邻近单元的最小距离δ。本文拓展文[1]的思想,进一步采用分部积分将δ移出奇异积分式中积分核之外,转换后积分核是δ的正则函数。所以几乎强奇异和超奇异积分被化为无奇异的规则积分与解析积分的和,可由通常的Gauss数值积分解。文中应用此正则化技术求解了弹性力学平面问题的近边界点位移和应力。     

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

结合了扩展有限元法(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通过应力法和位移法直接计算得到的裂尖应力强度因子均与解析解吻合\较好.      

边界变换法求边界近傍的应力

本文介绍的边界变换法,可以避免边界层效应,解决边界近傍的应力(或位移)的准确计算问题.     

压电介质二维边界积分方程中的基本解

由于压电介质的变形－电场耦合效应及压电响应的各向异性，使解析求解压电介质问题的工作变量十分复杂，若采用边界元数值方法求解，必须具备积分方程中的基本解，本文根据电磁场方程及连续介质力学的耦合性质论层出了二维无限域中分别在单位力及单位电荷载作用下的位移场，电势场、应力场和电位移场的解，从而确立了边界积分方程中所必需的八个基本解。     

一类新变量的Reissner板弯曲边界元法

文［４］导出了二维弹性力学平面问题的一类新型边界积分方程，本文将该理论和方法推广到三变量的Ｒｅｉｓｓｎｅｒ板弯曲中，给出边界场变量含广义位移和新型广义力的边界积分方程。从而边界弯矩应力张量可直接由离散边界积分方程求出。     

二维直角平面内固定圆形夹杂对稳态入射反平面剪切波的散射

利用复变函数法、多极坐标及傅立叶级数展开技术求解了二维直角平面内固定圆形夹杂对稳态入射反平面剪切（shearing horizontal, SH）波的散射问题。首先构造出介质内不存在夹杂时的入射波场和反射波场,然后建立介质内存在夹杂时由夹杂边界产生的能够自动满足直角边应力自由条件的散射波解，从而利用叠加原理写出介质内的总波场。利用夹杂边界处位移条件和傅立叶级数展开方法列出求解散射波中未知系数的无穷代数方程组，在满足计算精度的前提下通过有限项截断，得到相应有限代数方程组的解，最后通过算例具体讨论了二维直角平面水平边界点的位移幅度比和相位随量纲一波数、入射波入射角及夹杂位置的不同而变化的情况，结果表明了算法的有效实用性。     

物体边界层域内位移及应力的边界元法新公式

1.前言当分析塑性加工工具表层剥落破坏原因和求解断裂力学强度因子或其他工程问题时,仅计算物体表面应力是不够的,还需要计算边界层域内的位移及应力场。用边界元法计算上述工程问题时,当域内点和边界积分点距离足够小时,基本解中