基于边界元法的疲劳裂纹扩展模拟

用数值方法模拟疲劳裂纹扩展对预测构件残余寿命具有重要意义。该文以基于子域法边界元分析的一点面力法解决了裂纹扩展全过程的应力强度因子精度问题．用坐标变换法导出了Ⅰ－Ⅱ型复合裂纹最大周向应力准则的扩展偏角和有效应力强度因子；建立了复合型裂纹疲劳扩展数值模拟算法，并在ＳＧＩ工作站上研制了具有裂纹构形自动更新和动态图形显示功能的模拟演示软件，该软件对试件裂纹扩展模拟的结果与试验结果一致，说明该文方法具有高     

两个共面矩形裂纹的边界元分析

利用三维裂纹分析的边界积分方程方法，研究了三维无限弹性体中两个相等共面矩形裂纹的相互影响，将问题归结为求解三个二维的边界奇异积分方程并用二次元方法进行了数值计算。获得了一些典型问题的Ⅰ、Ⅱ、Ⅲ型应力强度因子。所得结果可供工程实际应用。     

裂纹问题的一种边界元解

根据线弹性断裂力学分析中常用的叠加方法,把一个有限域内的裂纹问题分解成无裂纹部分和裂纹部分的叠加来求解.其法只需较少的离散单元就可达到理想的精度,是求解裂纹问题的一种简单有效的方法.     

子域边界元法在双材料界面裂纹问题中的应用

采用子域边界元法对双材料界面边缘裂纹进行了研究，通过对两种材料分别沿材料边界进行单元划分，并利用在结合面上的两种材料的面力和位移的关系，得到边界上的所有未知分量，进而得到裂尖附近的位移场和应力场。应用外推法得到应力强度因子及裂尖周围的应力三维度。     

反平面裂纹问题的边界元解法

在传统边界积分公式的基础上运用分步积分等技巧,得到一个适用于求解反平面裂纹问题的新的边界积分方程,积分核只具有1／r阶的奇异性,在裂纹面上以位错密度为未知量,应力强度因子可由裂纹面上的位错密度求出,新的边界积分方程适用于任意形状的裂纹问题,两个数值算例证明了本文边界元法的正确性。     

非匀速动态裂纹应力强度因子的边界元算法

为了精确计算动态裂纹扩展时的裂纹应力强度因子,在Ｆｅｄｅｌｉｎｓｋｉ提出的时域对偶边界元法的基础上,将函数变换引入到弱奇异积分中,提高了弱奇异积分的精度,并将时域对偶边界元法应用到非匀速的裂纹扩展中,建立了裂纹非匀速扩展的边界元格式,进行了计算。将计算结果与根据文献数据对裂纹作的弹塑性有限元分析,及在此基础上模拟焦散法（ｃａｕｓｔｉｃｓ）所得的结果作了比较,两者符合得较好,说明边界元法在分析动态裂纹扩展中有很好的应用前景。     

多裂纹相互作用下断裂行为的边界元分析

基于边界元法,数值计算含多裂纹有限板在拉伸载荷下的断裂参数。利用裂纹面边界单元节点的张开位移,并结合线性拟合外推的方法计算应力强度因子。综合考虑计算精度和计算效率,讨论并选定合适的边界单元划分密度和类型。按照最大周向应力准则,通过断裂参数等效应力强度因子和扩展方向,分析裂纹I和II型复合断裂行为。裂纹间的相互作用会抑制或增强断裂行为,通过多个算例考查多裂纹的角度、长度、分布位置和间距等因素对抑制和增强作用的影响。     

一种新边界积分方程在三维裂纹分析中的应用

在推导出一种新的边界积分方程的基础上，分别对三维单个裂纹、多个周期裂纹进行了数值计算。单个裂纹的计算表明，最大裂纹张开位移能达到误差低于１．０％的精确结果；多裂纹的计算结果表明，在退化为二维问题时与解析解相吻合，该方法可用于断裂力学对多裂纹的研究。     

边界元法计算浅表面裂纹应力强度因子

文章在常规边界元法中引入几乎奇异积分的解析算法,来计算浅表面裂纹的应力强度因子,并给出了一个二维平板中有浅表面裂纹的算例.处理了几乎奇异积分的边界元法可以有效分析离表面仅为0.5 μm的浅表面裂纹.结果显示,随着距表面距离的减小,裂纹的应力强度因子收敛于常数,但其大小与无限大板中裂纹的应力强度因子有很大差别.     

半平面裂纹问题的边界积分方程

研究弹性半平面上的裂纹问题，得到一个适宜于求解各向同性半平面断裂力学问题的新边界积分方程，在裂纹面上以位错密度为未知量，以此求解应力强度因子．新的边界积分方程只具有1／r的奇异性，且适用于求解半平面上任意形状的裂纹问题．     

物性值随温度变化的热弹性边界元分析

提出求解物性值随温度变化的热弹性问题的边界元法，简述了基本积分方程的建立，导出有关计算公式．算例表明了本方法的有效性．     

Analysis of 3-D Frictional Contact Mechanics Problems by a Boundary Element Method ＊

The development of two boundary element algorithms for solving 3-D, frictional, and linear elastostatic contact problems is reported in this paper. The algorithms employ nonconforming discretizations for solving 3-D boundary element models, which provide much needed flexibility in the boundary element modeling for 3-D contact problems. These algorithms are implemented in a new 3-D boundary element code and verified using several examples. For the numerical examples studied, the results using the new boundary element algorithms match very well with the results using a commercial finite element code, and clearly demonstrate the feasibility of the new boundary element approach for 3-D contact analysis     

近似边界元方法及其收敛性分析

以Ｌａｐｌａｃｅ方程Ｄｉｒｉｃｈｌｅｔ问题为例,为椭圆边值问题近似边界元法的建立及其收敛性分析提供了一种框架性的工作。文中给出了近似基本解的确定方法,具近似基本解的离散边界变分方程解的存在惟一性定理以及近似解的误差估计,特别给出了近似基本解中截断数和离散网格宽度应保持的匹配关系,文末给出了数值算例。     

功能梯度材料中裂纹问题的边界元分析

采用基于边界元方法的广义Kelvin解对功能梯度材料中的裂纹问题进行了研究，主要在对裂纹的评价中采用了多域法和八节点面力奇异边界单元，并采用层离散化方法用来近似功能梯度材料，计算出了不同条件下功能梯度材料中币形裂纹的应力强度因子。     

基于边界元方法的边值问题数值解的外推

利用边界元方法求解椭圆边值问题,并通过Poisson积分方程的Galerkin 解讨论了这种方程的外推算法,进而对边值问题的数值解获得了Ｏ（ｈ３）精度的外推结果.     

Signorini接触问题的边界元方法及收敛性分析

以Signorini接触问题为例，讨论了接触问题边界变分不等式的边界元方法，得到了离散边界变分不等式近似解的存在唯一性定理，给出了近似解与精确解的误差估计表达式。     

拉氏变换域中三维广义热弹性力学的边界积分方程及其基本解

将拉氏变换域中热传导方程的耦合项“线性化”处理后，精确地敢出了拉氏变换域中的基本解，并用位移互等定理导出了问题的边界积分方程。     

推力轴承的边界元分析

目前边界元法还很少用于润滑问题的研究,原因在于雷诺方程的基本解很难求得.本文通过拟合推力轴承的油膜厚度,将雷诺方程转化为拉普拉斯方程的形式,用边界元法很快就得到了推力轴承的压力分布.计算结果表明,该方法的计算精度在实际工程中是能够被接受的.