移动接触下求解二维闭合裂纹的双重迭代边界元法

使用子域边界元法对受移动接触弹性体作用下的二维闭合裂纹问题进行了数值计算。由于两弹性体的接触界面和裂纹表面的接触范围的大小和接触状态事先是未知的 ,对此 ,在两个接触表面同时采用迭代的方法进行了求解。在裂纹的每个裂尖上都采用了四分之一的奇异单元以保证裂尖位移场和应力场奇异性的满足。用我们编制的二维裂纹问题程序对一些中心裂纹问题进行了计算 ,计算结果与经典断裂力学的理论值比较吻合。在无摩擦的条件下 ,对一些具有不同角度且受移动接触弹性体作用下的闭合裂纹问题进行了数值计算 ,得到了一些耦合作用下的应力强度因子的计算结果     

正交各向异性摩擦接触分析的一个二阶锥线性互补法

建立了用于各向异性摩擦接触分析的一个新型互补模型:二阶锥线性互补模型.通过与弹塑性力学的类比,并引入"线性二阶锥互补"的概念,将正交各向异性摩擦接触定律表达为一个二阶锥线性互补模型.在此基础上,在参变量变分原理的框架下,建立了正交各向异性摩擦接触问题的有限元列式.由于避免了对摩擦锥的显式线性化,使得该文模型的变量个数较线性互补模型小很多,而解的精度又高很多.文末通过算例,说明了方法对于正交各向异性摩擦接触分析的有效性.     

考虑摩擦作用闭合裂纹应力强度因子计算的边界元分区算法

本文采用边界元分区算法研究考虑摩擦的闭合裂纹问题,简单裂纹系问题及非均匀介质中裂纹问题.在裂纹尖端采用了1/4面力奇异单元,并对相应的奇异积分给出了数值处理.对于复杂载荷下裂纹面计算给出了增量迭代算法,并采用方程减缩技术使迭代仅在裂纹上进行.计算实例表明方法是可行的.     

基于互补理论和扩展有限元法的摩擦接触裂纹应力强度因子计算

从虚功方程出发,结合扩展有限元离散技术与接触条件的非线性互补表述,建立了摩擦接触裂纹问题的扩展有限元非线性互补模型,将不等式接触条件转化为非线性互补类的非光滑方程组,并采用基于广义导数的非光滑阻尼牛顿法求解方程组,无需引入任何额外人工变量以及迭代求解。以含中心倾斜裂纹平板和边裂纹平板为例,运用相互作用积分法计算摩擦接触裂纹的应力强度因子,将其结果与理论解进行对比分析,该方法都能给出精确的计算结果;基于扩展有限元方法对单轴压缩作用下倾斜裂纹扩展过程进行了数值模拟,计算结果表明,受压裂纹数值结果与实验结果比较吻合,从而验证了本文方法的有效性与正确性。     

缺陷孔对平面裂纹应力强度因子的影响

带有裂纹和缺陷孔洞的板的问题是一个多连域的边值问题，这类问题适合用边界元法所具有的高精度特性来求解．采用子域边界元法，在平面应变的条件下对存在中心裂纹的平面板受远处拉伸和剪切裁荷的作用进行了数值分析．研究了圆形孔洞对Ⅰ型和Ⅱ型应力强度因子的影响，与有限元法进行了对比，求解结果更加精确．计算了椭圆形孔对Ⅰ型应力强度因子的影响，得到了一些有意义的结果，并对移动接触弹性体作用下的带裂纹板进行了钻孔研究．     

任意多椭圆孔多裂纹无限大各向异性板应力强度因子求解的一种新方法

采用各向异性体平面弹性理论中的复势方法,应用保角变换技术,以F aber级数为工具,导出含任意多椭圆孔和裂纹群无限大各向异性板在远场载荷作用下其应力场和位移场的级数解,并在此基础上利用断裂力学方法确定裂纹尖端的应力强度因子,通过算例讨论了材料参数及裂纹、孔的尺寸等对应力强度因子的影响规律,得出了一些有益的结论。数值结果表明本文方法具有计算精度高、收敛速度快、方便快捷等优点,有利于全面系统地研究各参数对结构断裂性能的影响。     

比例边界有限元法求解裂纹面接触问题

比例边界有限元侧面上有任意荷载时,将侧面载荷分解成关于径向方向局部坐标的多项式函数的和,推导给出了考虑侧面载荷存在的新型形函数,并基于该形函数推导了刚度矩阵和等效节点载荷列阵.首次对比例边界有限元法求解裂纹面接触问题进行了研究,运用Lagrange乘子引入接触界面约束条件,推导给出了比例边界有限元求解裂纹面接触问题的控制方程.将裂纹面单元分为非裂尖单元和含有侧面的裂尖单元.在非裂尖单元中的裂纹面,裂纹面作为多边形单元的边界,边界上的接触力可等效到节点上,通过在节点上构造Lagrange乘子,采用点对点接触约束进行处理.对于含有侧面的裂尖单元,在整个侧面上构造Lagrange乘子的插值场,采用边对边接触约束进行处理.对三个不同的接触约束状态下的算例进行了数值计算,通过与解析解及有限元软件ABAQUS计算结果的对比,验证了本文提出的比例边界有限元点对点和边对边接触求解裂纹面接触问题的精确性与有效性.     

用子域边界元法研究各向异性材料中的界面裂纹

用子域边界元法研究各向异性材料中的界面裂纹，在边界元公式中，采用了带特征根的基本解，以增量形式的边界积分方程为基础，通过二次等参元及国分之一面力奇异远离散化处理，可以得到各子域的代数方程组，依据凝集技术，可得到仅含有子域公共边界及裂纹边界未知量的求解方程组，通过迭代法，可以寻求到每种载荷作用下的裂纹所处的真实状态，然后，由文献「２」中的方法求解界面裂纹的应力强度因子。结果表明，子域边界方法是正确的     

界面裂纹尖端的动态奇异特性

本文研究了界面裂纹尖端的动态应力场的奇异特性.引入尖端无摩擦接触的界面裂纹模型并采用具有运动边界的控制积分方程.证明了在动态界面裂纹尖端仅存在平方根奇异的应力场.数值结果表明接触区中的正应力确保持为压应力.为表现界面裂纹的动态特性,给出了应力强度因子和裂纹面接触区尺寸的数值结果.     

基于近场动力学微分算子的含裂纹正交各向异性板热传导分析方法

采用近场动力学微分算子(Peridynamic Differential Operator, PDDO)理论建立正交各向异性板热传导的非局部模型。通过构造近场动力学函数,将边界条件和热传导方程由局部微分形式转化为非局部积分形式,引入Lagrange乘子,采用变分分析对含裂纹正交各向异性板温度及裂纹尖端的热通量分布进行求解。通过对比算例,验证了该模型具有较好的收敛性和有效性。分析了正交各向因子、材料铺设角、裂纹倾角及间距对裂纹尖端热通量的影响。结果表明,基于PDDO建立的含裂纹正交各向异性板热传导模型,考虑了热传导问题中的非局部性,能有效提高计算精度,预测含裂纹板中裂纹尖端出现的奇异性。     

On interface crack models with contact zones situated in an anisotropic bimaterial

Summary A boundary value problem for two semi-infinite anisotropic spaces with mixed boundary conditions at the interface is considered. Assuming that the displacements are independent of the coordinate x ₃, stresses and derivatives of displacement jumps are expressed via a sectionally holomorphic vector function. By means of these relations the problem for an interface crack with an artificial contact zone in an orthotropic bimaterial is reduced to a combined Dirichlet-Riemann problem which is solved analytically. As a particular case of this solution, the contact zone model (in Comninou's sense) is derived. A simple transcendental equation and an asymptotic formula for the determination of the real contact zone length are obtained. The classical interface crack model with oscillating singularities at the crack tips is derived from the obtained solution as well. Analytical relations between fracture mechanical parameters of different models are found, and recommendations concerning their implementation are given. The dependencies of the contact zone lengths on material properties and external load coefficients are illustrated in graphical form. The practical applicability of the obtained results is demonstrated by means of a FEM analysis of a finite-sized orthotropic bimaterial with an interface crack. Received 19 October 1998; accepted for publication 13 November 1998     

任意多孔多裂纹有限大板的应力强度因子分析

采用各向异性体平面弹性理论中的复势方法,以Faber级数为工具,应用保角映射技术和最小二乘边界配点法,导出内边界条件精确满足,外边界条件近似满足的含多椭圆孔及裂纹群有限大板在任意载荷作用下的应力场、位移场的级数解,建立了任意多椭圆孔及裂纹群有限大板应力强度因子的有效分析方法,讨论了各参数对裂尖应力强度因子及孔边应力集中的影响．数值结果表明,该方法具有计算精度高、收敛速度快、方便快捷等优点,有利于全面系统地研究各参数对结构断裂性能的影响．     

Dynamic stress intensity factor KⅢ and dynamic crack propagation characteristics of anisotropic materials

Based on the mechanics of anisotropic materials,the dynamic propagation problem of a mode Ⅲ crack in an infinite anisotropic body is investigated.Stress,strain and displacement around the crack tip are expressed as an analytical complex function,which can be represented in power series.Constant coefficients of series are determined by boundary conditions.Expressions of dynamic stress intensity factors for a mode Ⅲ crack are obtained.Components of dynamic stress,dynamic strain and dynamic displacement around the crack tip are derived.Crack propagation characteristics are represented by the mechanical properties of the anisotropic materials,i.e.,crack propagation velocity M and the parameter α.The faster the crack velocity is,the greater the maximums of stress components and dynamic displacement components around the crack tip are.In particular,the parameter α affects stress and dynamic displacement around the crack tip.     

Dynamic stress intensity factor K III and dynamic crack propagation characteristics of anisotropic materials

Based on the mechanics of anisotropic materials, the dynamic propagation problem of a mode Ⅲ crack in an infinite anisotropic body is investigated. Stress, strain and displacement around the crack tip are expressed as an analytical complex function, which can be represented in power series. Constant coefficients of series are determined by boundary conditions. Expressions of dynamic stress intensity factors for a mode Ⅲ crack are obtained. Components of dynamic stress, dynamic strain and dynamic displacement around the crack tip are derived. Crack propagation characteristics are represented by the mechanical properties of the anisotropic materials, i.e., crack propagation velocity M and the parameter ~. The faster the crack velocity is, the greater the maximums of stress components and dynamic displacement components around the crack tip are. In particular, the parameter α affects stress and dynamic displacement around the crack tip.     

Hypersingular integral equation method for a three-dimensional crack in anisotropic electro-magneto-elastic bimaterials

Using Green’s functions, the extended general displacement solutions of a three-dimensional crack problem in anisotropic electro-magneto-elastic (EME) bimaterials under extended loads are analyzed by the boundary element method. Then, the crack problem is reduced to solving a set of hypersingular integral equations (HIE) coupled with boundary integral equations. The singularity of the extended displacement discontinuities around the crack front terminating at the interface is analyzed by the main-part analysis method of HIE, and the exact analytical solutions of the extended singular stresses and extended stress intensity factors (SIFs) near the crack front in anisotropic EME bimaterials are given. Also, the numerical method of the HIE for a rectangular crack subjected to extended loads is put forward with the extended crack opening dislocation approximated by the product of basic density functions and polynomials. At last, numerical solutions of the extended SIFs of some examples are obtained.     

复合材料承受剪切时双裂纹周围的应力强度因子

本文研究了由各向同性和各向异性半无限接合而成的复合材料中的应力强度因子问题，在复合材料的接合面附近处具有与接合面平行且共线的两个Ｇｒｉｆｆｉｔｈ裂纹，裂纹面上作用有剪应力，本文利用付利叶变换将混合边值问题归毕为求解奇异积分方程问题，为求解这些方程，将裂纹面上，下的位移差展成级数，并满足理解纹面外侧边界条件，级数中的待定系数利用裂纹面内的边界条件和施密特方法求得，本文对硼纤维塑料和铝板接合的复合材料     

Study on failure process of delaminated stiffened composite plates under compression

Failure behavior of the delaminated stiffened composite plates under compression is studied by the finite element method, based on a Global-Local variational model. A virtual crack closure technique and a self-adaptive grid moving scheme are proposed to predict the delamination growth process. The contact effect along the delamination front is considered. The numerical results show that the influences of the distribution and location of the stiffeners, the configuration and size of the delamination, the boundary condition and the contact upon the failure behavior of the plates are significant. The project supported by the National Natural Science Foundation of China (59975014)     

Analysis of bonded anisotropic wedges with interface crack under anti-plane shear loading

The antiplane stress analysis of two anisotropic finite wedges with arbitrary radii and apex angles that are bonded together along a common edge is investigated. The wedge radial boundaries can be subjected to displacement-displacement boundary condi- tions, and the circular boundary of the wedge is free from any traction. The new finite complex transforms are employed to solve the problem. These finite complex transforms have complex analogies to both kinds of standard finite Mellin transforms. The traction free condition on the crack faces is expressed as a singular integral equation by using the exact analytical method. The explicit terms for the strength of singularity are extracted, showing the dependence of the order of the stress singularity on the wedge angle, material constants, and boundary conditions. A numerical method is used for solving the resul- tant singular integral equations. The displacement boundary condition may be a general term of the Taylor series expansion for the displacement prescribed on the radial edge of the wedge. Thus, the analysis of every kind of displacement boundary conditions can be obtained by the achieved results from the foregoing general displacement boundary condition. The obtained stress intensity factors （SIFs） at the crack tips are plotted and compared with those obtained by the finite element analysis （FEA）.