三维有限体平片裂纹的超奇异积分方程与边界元法

利用Ｓｏｍｉｇｌｉａｎａ公式及有限部积分的概念，导出了含任意平片裂纹三维有限体问题的超奇异积分方程组，并联合使用有限部积分与边界元法，建立了数值求解方法．在裂纹前沿附近单元，采用与理论分析一致的平方根位移模型，以提高数值结果的精度．最后计算了若干典型例子的应力强度因子．     

用退化变换法消除三维弹塑性裂纹分析边界元模式中的奇异性

本文提出了一套坐标变换关系式,运用于带1/r,1/r~2,1/r~3乃至1/r~4型奇异性的积分运算,并可有效地消去积分核函数的奇异性,显著提高数值积分的精度。所述方法应用于三维弹塑性裂纹分析,对相应边界元表达式中所有奇异面积分(弹性区)和体积分(局部塑佳区)都作了恰当处理,并通过实例计算和有关理论结果进行了比较,表明吻合良好。     

多裂纹扩展分析的边界元方法

采用边界元数值模拟和即时等效材料常数计算相结合的办法,只需模拟一个裂纹的扩展情况便可预测裂纹体的整体响应.针对准脆性材料的特点,采用粘性裂纹模型模拟裂纹开裂行为;采用二次裂纹扩展量作为增量控制变量,避免了软化及失稳分析中用力或位移作控制变量时遇到的困难.针对二次裂纹扩展路径未知情况,给出了预测和修正裂纹开裂界面的迭代技术.分析计算了含多个规则分布裂纹石膏板受压时的响应,并与实验结果进行了比较.结果表明该文方法的可行性与有效性.     

各向异性柱体扭转的充分必要的边界积分方程

本文列出各向异性柱体扭转问题无量纲化后的有关方程的边界条件推导和验证了基本解，并指出一些书中基本解列式有误＾［９－１１］，列出了充分必要的边界积分方 程，进行了数据计算，并与习用的边界积分方程所得结果进行了比较，表明在退化值附近，习用的边界积分方程所得的边界剪应力会出现巨大的误差，扭转刚度的误差则要小得多，而充要的边界积分方程计算的结果则如终保持良好的精度，再次显示了它的优点。     

二维边界元奇异积分和多域缩聚法分析

基于基本解的一种新的表达式，对二维边界元分析中奇异积分的精确求解进行了讨论，从几何方面对基本解的奇异性进行了分析，给出了超参非连续元离散位势和弹性力学问题边界积分方程时奇异积分计算的精确式，从而为判断各种近似方法的优劣和间接方法的精度提供了依据，也为精确地分析了大规模问题提供了一条有效的途径。     

三维裂纹扩展轨迹的边界元数值模拟

提出了一种对三维裂纹扩展轨迹进行数值模拟的新方法。采用一种新的具有C^1连续性、高精度的单节点二次边界单元，使边界元(BEM)的分析效率和裂纹张开位移(COD)、应力强度因子(SIF)的精度大大提高。采用裂纹张开位移全场拟合法(GCDFP)求出裂纹面前缘的SIF，所得到的SIF达到与所用的COD资料同样的精度。使用Paris公式求出裂纹前缘各点的裂纹扩展增量，并用三次B样条函数对这些增量进行拟合，得到新的光滑裂纹前缘。根据以上思想方法，开发了具有较高的计算效率和精度的数值模拟软件。此软件可以自动跟踪裂纹扩展，得到裂纹扩展的轨迹。运用该软件对椭圆和矩形裂纹的扩展轨迹进行了数值模拟。其结果与理论上的预言完全一致，裂纹最后都趋于一个圆裂纹，具有实际指导意义。     

弹性力学平面问题的无奇异边界积分方程

将弹性力学平面问题归化成无奇异边界积分方程,避免了传统的边界元法中的柯西主值（CPV）积分和Hadamard-Finite-Parts(HFP)积分的计算,建立完整的数值求解体系。     

确定裂纹体等效弹性模量的边界元方法

采用边界元方法计算含有序分布裂纹的裂纹体在压缩载荷作用下的等效弹性模量，利用一种能适当考虑裂纹有间相作用的自洽理论，建立了相应的迭代格式，通过算例研究了裂纹方向，裂纹面间摩擦系数对裂纹体等效弹性模量的影响。     

Torsion problems for a cylinder with a rectangular hole and a rectangular cylinder with a crack

Using the method of singular integral equation and the crack-cutting technique, the rigorous solutions are obtained for a cylinder with a rectangular hole and a rectangular cylinder with a crack, which exactly satisfy the boundary conditions and the conditions at the corner points. After that the torsional rigidities and the stress intensity factors at the crack tip are determined. Next, for the doubly connected circular cylinder with a rectangular hole the expressions for the singular stresses around the concave corner points are derived and the generalized stress intensity factors are then defined. Since the crack-cutting technique is used in this paper, the solution of the matching rectangular cylinder is also obtained and its numerical results coincide with those in references. Thus the method proposed here is verified. The project supported by National Natural Science Foundation of China     

Fracture analysis for torsion problems of a composite cylinder with curvilinear cracks

The Salnt-Venant torsion problems of a composite cylinder with curvilinear cracks were investigated. By considering the bimaterial interface as a boundary of the outer bar or inner one, the problem was reduced to the solution of boundary integral equations on the crack, external boundary and interface. Using the new boundary element method, some typical torsion problems of a composite cylinder involving a straight or kinked crack were calculated. The obtained results were compared with data in the literature to show validity and applicability of the present method.     

The solution of integral equations with strongly singular kernels applied to the torsion of cracked circular cylinder

In this paper,the functions of warping displacement interruption defined on the crack lines are taken for the fundamental unknown functions.The torsion problem of cracked circular cylinder is reduced to solving a system of integral equations with strongly singular kernels.Using the numerical method of these equations,the torsional rigidities and the stress intensity factors are calculated to solve the torsion problem of circular cylinder with star-type and other different types of cracks.The numerical results are satisfactory.     

Solving the Three-Dimensional Elastic Problem for a Cylinder of Finite Length Subject to Bending

A new technique is proposed for solving the three-dimensional antisymmetric elastic problem for a thick-walled shell of finite length. The boundary-value problem is reduced to a system of one-dimensional singular integral equations. Characteristic stresses are calculated     

Plane Problem of Three-Dimensional Stability for a Hinged Plate with Two Symmetric End Cracks

The plane stability problem for a rectangular plate with two symmetric end cracks is solved in three-dimensional formulation. The three-dimensional linearized theory of stability and the finite-difference method are used. The effect of the crack parameter on the critical load is examined__________Translated from Prikladnaya Mekhanika, Vol. 41, No. 4, pp. 47–52, April 2005.     

On the Torsion of a Class of Nonlinear Fiber Composite Cylinders

This paper presents an analysis of the torsion of a solid or annular circular cylinder consisting of nonlinear material in the form of an elastic matrix with embedded unidirectional elastic fibers parallel to the cylinder axis. The specific class of composite considered is one for which nonlinear fiber-matrix interface slip is captured by uniform cohesive zones of vanishing thickness. Previous work on the effective antiplane shear response of this material leads to a stress–strain relation depending on the interface slip together with an integral equation governing its evolution. Here, we obtain an approximate single mode solution to the integral equation and utilize it to solve the torsion problem. Equations governing the radial distributions of shear stress and interface slip are obtained and formulae for torque–twist rate are presented. The existence of singular surfaces, i.e., surfaces across which the slip and the shear stress experience jump discontinuities are analyzed in detail. Specific results are presented for an interface force law that allows for interface failure in shear.     

横观各向同性圆柱体端部问题的三维弹性理论解

从位移的通解出发，用分离变量法得到横观各向同性圆柱体的位移和应力的特征函数展开式，并把位移势函数的解用付里叶积分的形式表示。利用留数运算，该积分解可以转换成类似于特征函数的展开式。通过混合端部边界问题，得到与特征函数解成双正交关系的另一组函数。利用这种双正交关系，可以处理不同的端部边界问题。     

Stress intensity factors of the eccentric and edge cracked cylinders

Using the single crack solution and the regular solution of harmonic function,thetorsion problem of a cracked cylinder is reduced to solving a set of mixed-type integralequations which can be solved by combining the numerical method of singular integralequation with the boundary element method.Several numerical examples arecalculated and the stress intensity factors are obtained.     

一维六方准晶中星形静态裂纹和运动裂纹的解析解

利用复变函数方法研究了一维六方准晶中星形静态裂纹和运动裂纹的反平面剪切问题,得到了星形裂纹尖端处应力强度因子和动应力强度因子的解析解.当裂纹条数给定时,由此可得到直线裂纹,Griffith裂纹,共点均匀分布三裂纹,对称十字形裂纹,米字型裂纹(对称八裂纹)静力学和动力学问题的解析解.当k=4时,用数值算例讨论了声子场-相位子场耦合系数和裂纹运动速度对动应力强度因子的影响.当速度趋于0时,运动裂纹的解可以退化为静态裂纹的解.     

非平面应变状态下的叠层厚壁筒

抛弃任何有关位移或应力模式的人为假设，在轴对称情况下，导出正交异性厚壁筒的状态方程。在沿筒轴方向任意分布的轴对称荷载下，给出叠层厚壁筒静力问题的精确解。此解满足所有弹性力学基本方程，包含了全部弹性常数，可满足任意精度要求。数值结果和ＳＡＰ５有限元解进行了对比。