正交各向异性平面问题应力强度因子的边界元分析

本文给出了正交各向异性平面弹性问题的边界元方程，导出了常单元离散化时求系数的解析式。作为数值算例，计算了正交各向异性板的应力强度因子。结果表明，本文所导出公式的正确性。     

裂纹面受荷载作用的应力强度因子的计算

基于比例边界有限元法计算了裂纹面有荷载作用情况下裂纹尖端的应力强度因子,给出了有限介质裂纹面作用荷载的比例边界有限元方程的基本求解过程.对于随径向坐标任意变化的一类面荷载的积分能够显式计算,不需要引入额外的近似;并将计算结果与解析解和数值结果进行对比,结果表明比例边界有限元法在计算裂纹面作用荷载时的应力强度因子是有效且精确的.此外,该方法可方便地处理各向异性材料裂纹问题,本文给出了正交各向异性矩形盘裂纹面受均布荷载情况的应力强度因子.     

正交各向异性平面热弹性问题的边界元分析

给出了正交各向异性平面热弹性问题的边界积分方程、内点应力公式和常单元离散化时计算奇异积分的解析式，计算了正交各向异性板的热应力强度因子，结果表明了文中所导公式的正确性     

正交异性光弹性应力分离的边界元法

本文对平面正交各向异性复合材料模型引入正应力线性和及边界上正应力线性和流的概念,提出从应力相容方程出发.用边界元法计算正交异性光弹性模型内任一点的正应力线性和位的方法,再与正交异性光弹性法中所给出的应力同的关系结合,即可进行正交异性光弹性应力的分离.最后,对边界元方法的精度进行了讨论.     

各向异性介质波动问题的时域边界元法及实验验证

对前已建立的各向异性介质波动问题的时域边界元计算模型进行了实验验证．用单向纤维增强光弹性复合材料模拟正交异性介质，用冲击加载，加载方向与纤维方向分别成0°，90°及45°角度，进行了正交异性动态光弹性实验及动态应变测量，并同时对该模型进行了时城边界元计算．将时域边界元方法计算出的应力分量代入正交异性动态光弹性的动态应力-光性定律，得到双折射条纹级数随时间的变化曲线，将其与动态光弹性实验的结果进行比较；此外，由动态电测获得的应变响应曲线推算出应力时程；与时城边界元计算出的应为响应曲线也进行了比较．两种情况下，时域边界元的计算成果均与实验成果吻合较好，从而证明该各向异性介质波动问题的时域边界元计算模型具有较好的精度及稳定性，可用于各向异性介质的动态问题特别是波传播问题的分析研究．     

含中心裂纹正交各向异性有限大板的反平面剪切问题

本文采用边界配置方法计算了含中心裂纹正交各向异性有限大板问题的应力强度因子。在其特例——材料各向同性情形,本文结果与[2]采用Fourier变换方法所得结果一致;同时,本文给出了不同材料参数和裂纹长度情形的计算结果。     

一种新的计算各向异性材料裂纹尖端应力强度因子杂交元法

利用有限元特征分析法研究了平面各向异性材料裂纹端部的奇性应力指数以及应力场和位移场的角分布函数,以此构造了一个新的裂纹尖端单元。文中利用该单元建立了研究裂纹尖端奇性场的杂交应力模型,并结合Hellinger-Reissner变分原理导出应力杂交元方程,建立了求解平面各向异性材料裂纹尖端问题的杂交元计算模型。与四节点单元相结合,由此提出了一种新的求解应力强度因子的杂交元法。最后给出了在平面应力和平面应变下求解裂纹尖端奇性场的算例。算例表明,本文所述方法不仅精度高,而且适应性强。     

用边界元方法分析复合材料中的裂纹问题

利用层状材料的广义Klevin基本解，建立了计算三维层状材料中的裂纹边界元方法。采用边界元方法中的多区域方法和能反映均匀介质中裂纹尖端应力场和位移场特征的面力奇异单元。裂纹的应力强度因子由裂纹面上的位移经插值计算得到。算例分析表明，本文建议的方法可以获得较高的计算精度。     

表面钝裂纹的计算模型及其边界元法模拟

研究了表面钝裂纹问题的边界元模拟方法。文中通过平面应变比拟，建立了三维钝型纹的计算模型和局部场结构；并由三维边界元程序计算了表面钝裂纹前沿附近的位移场和应力场；进而利用裂纹面前沿的“张开位移”推算应力强度因子的分布文中的钝裂纹模型的有效性和离散格式的收敛性进行了考核，应力强度因子计算针对含表面裂纹的平板进行。     

含裂纹粘接结构的边界元应力分析

本文采用子域法和直接应力奇异单元法求解二维粘接结构中的裂纹问题。子域法把粘接结构划分为三个子域,根据每个子域的边界积分方程和子域间的界面条件,可以建立粘接结构的边界积分方程组。直接应力奇异单元能够在整个单元长度上反映裂纹端部的1/r~(1/2)奇异性,在计算时可以通过坐标变换消除奇异单元积分中的奇异性,直接计算出应力强度因子。含裂纹多层结构的数值示例和粘接补强单边裂纹板的应力测试和疲劳试验结果证实了本文方法的有效性。     

Application of an improved boundary element method on the problem of elastic torsion

An improved boundary clement method has been used in analyzing and calculating the problems of the torsion of a prismatic bar with elliptical cross-section. In this paper the calculated results correspond with the values of boundary element method. However, the quantity of data required by the improved boundary element method is much less than that required by boundary element method, and the calculating time will be greatly reduced. Therefore, the procedure of this paper is an economical and efficient numerical computational way for solving Poisson equation problem.     

An effective boundary element method for analysis of crack problems in a plane elastic plate

A simple and effective boundary element method for stress intensity factor calculation for crack problems in a plane elastic plate is presented. The boundary element method consists of the constant displacement discontinuity element presented by Crouch and Starfield and the crack-tip displacement discontinuity elements proposed by YAN Xiangqiao. In the boundary element implementation the left or the right crack-tip displacement discontinuity element was placed locally at the corresponding left or right each crack tip on top of the constant displacement discontinuity elements that cover the entire crack surface and the other boundaries. Test examples (i. e. , a center crack in an infinite plate under tension, a circular hole and a crack in an infinite plate under tension) are included to illustrate that the numerical approach is very simple and accurate for stress intensity factor calculation of plane elasticity crack problems. In addition, specifically, the stress intensity factors of branching cracks emanating from a square hole in a rectangular plate under biaxial loads were analysed. These numerical results indicate the present numerical approach is very effective for calculating stress intensity factors of complex cracks in a 2-D finite body, and are used to reveal the effect of the biaxial loads and the cracked body geometry on stress intensity factors.     

V形切口应力强度因子的一种边界元分析方法

将V形切口结构分成围绕切口尖端的小扇形和剩余结构两部分. 尖端处扇形域应力场表示成关于尖端距离$\rho$的渐近级数展开式,从线弹性理论方程推导出了一组分析平面V形切口奇异性的常微分方程特征值问题,通过求解特征方程,得到前若干个奇性指数和相应的特征向量. 再将切口尖端的位移和应力表示为有限个奇性阶和特征向量的组合. 然后用边界元法分析挖去小扇形后的剩余结构. 将位移和应力的线性组合与边界积分方程联立,求解获得切口根部区域的应力场、应力幅值系数和整体结构的位移和应力. 从而准确计算出平面V形切口的奇异应力场和应力强度因子.      

双轴载荷作用下源于椭圆孔的分支裂纹的一种边界元分析

利用一种边界元方法来研究双轴载荷作用下无限大板中源于椭圆孔的分支裂纹．该边界元方法由Crouch与Starfied建立的常位移不连续单元和笔者提出的裂尖位移不连续单元构成．在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界,文中算例说明本数值方法对计算平面弹性裂纹的应力强度因子是非常有效的。该文对双轴载荷作用下无限大板中源于椭圆孔的分支裂纹的数值结果进一步证实本数值方法对计算复杂裂纹的应力强度因子的有效性,同时该数值结果可以揭示双轴载荷及裂纹体几何对应力强度因子的影响。     

双材料界面裂纹应力强度因子的边界元分析

采用双材料基本解建立边界元法基本方程,计算双材料界面裂纹尖端附近的应用力和位移场。不离散界面,并设置面力奇异四分之一点裂尖单元以提高计算精度。数值结果表明,本文的方法具有较高的精度和效率。     

复杂载荷三维裂纹分析双重边界元法

提出可用于高温、高转速状态下的热动力机械三维含裂构件热弹性分析方法——双重边界元法.首先建立了考虑温度及离心载荷的双重边界积分方程组,并对边界积分方程组的选取及适用范围进行了讨论。然后提出角非快调元模型离散技术。接着提出超奇异积分方程分析去除奇异性方法及数值积分技术.数值算例表明计算结果与有关权函数解十分吻合,说明了用双重边界元法计算复杂载荷条件下三维应力强度因子的有效性.还讨论了有关热应力强度因子权函数解的适用范围.     

若干综合应用递归二次规划法和边界元法的平面弹性结构形状优化问题

结构的边界表示为若干设计变量的函数,结构形状优化问题表示为数学规划问题。本文采用递归二次规划法求解数学规划问题,采用边界元法做结构分析,求解了受拉多边形板、受弯悬臂梁和空腹重力坝的形状优化问题。结果表明本文的求解方案非常有效。     

Two dimensional stress intensity factor computations by multi-domain boundary element method for crack closure with friction

A multi-domain boundary element method is used to compute the stress intensity factor of plane stress/plane strain crack problems with friction. The analysis is performed by using traction-singular quarter-point boundary elements on each side of the crack tips. The increment iteration is given. The technique is applied to some specific examples in order to show that the results will be with good accuracy.The Project 13 supported by National Natural Science Foundation of China.     

内部压力作用下矩形板中源于椭圆孔的分支裂纹应力强度因子的一种数值分析

应用一种边界元方法来研究内部压力作用下矩形板中源于椭圆孔的分支裂纹。该边界元方法由Crouch与Starfied建立的常位移不连续单元和笔者最近提出的裂尖位移不连续单元构成。在该边界元方法的实施过程中,左、右裂尖位移不连续单元分别置于裂纹的左、右裂尖处,而常位移不连续单元则分布于除了裂尖位移不连续单元占据的位置之外的整个裂纹面及其它边界。本数值结果进一步证实这种数值方法对计算有限大板中复杂裂纹的应力强度因子的有效性,同时该数值结果可以揭示裂纹体几何对应力强度因子的影响。