边界元法计算切口多重应力奇性指数

提出采用边界元法直接计算V形切口的多重应力奇性指数。首先在切口尖端挖出一微小扇形域,在该域边界列常规边界积分方程,后将扇形域内的位移场和应力场表示成关于切口尖端距离ρ的渐近级数展开式,回代入切口边界积分方程,离散后得到关于切口奇性指数的代数特征方程,从而求解获得V形切口的应力奇性指数。该法避免了常规边界元法和有限元法在切口尖端附近布置细密单元的缺陷,并可同时求得多阶应力奇性指数。     

双材料基本解弹塑性边界元法

提出并研究采用双材料基本解的弹塑性边界元法,得到了内点应力公式中有关奇点塑性应变自由项的完整表达式,并利用非连续边界单元和非连续区域单元解决了当奇点位于界面上时该自由项难于确定,以及计算区域Cauchy主值积分的常塑性应变场法在与界面相连的奇异区域单元上无法实施的困难．采用双材料基本解的弹塑性边界无法针对双材料的结构特点,特别适于分析有关弹塑性双材料界面及界面裂纹问题．     

Solution of elastic-plastic shallow shell problems by the boundary element method

The boundary integral equations for elasto-plastic problems of shallow shells are established by using the fundamental solutions of shallow shells derived previously. The strains and stress-resultants in the plastic region are used as unknown variables. The simultaneous nonlinear equations of these variables and unknown boundary values are established and solved by direct iteration method.     

Analytical solutions for finite cylindrical dynamic cavity expansion in compressible elastic-plastic materials

Analytical solutions for the dynamic cylindrical cavity expansion in a com-pressible elastic-plastic cylinder with a finite radius are developed by taking into account of the effect of lateral free boundary, which are different from the traditional cavity expan-sion models for targets with infinite dimensions. The finite cylindrical cavity expansion process begins with an elastic-plastic stage followed by a plastic stage. The elastic-plastic stage ends and the plastic stage starts when the plastic wave front reaches the lateral free boundary. Approximate solutions of radial stress on cavity wall are derived by using the Von-Mise yield criterion and Forrestal’s similarity transformation method. The effects of the lateral free boundary and finite radius on the radial stress on the cavity wall are discussed, and comparisons are also conducted with the finite cylindrical cavity expansion in incompressible elastic-plastic materials. Numerical results show that the lateral free boundary has significant influence on the cavity expansion process and the radial stress on the cavity wall of metal cylinder with a finite radius.     

近场波动模拟的一种应力人工边界

采用平面波和远场散射波混合透射，引入无限介质线弹性本构关系建立了一种应力人工 边界条件. 其优点在于边界结点反应与内部有限元结点反应采用相同的积分格式计算，有限 元积分方法稳定时不存在人工边界失稳问题. 数值算例表明：边界精度高于现有黏性边 界、黏弹性人工边界，以及一、二阶透射人工边界.     

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

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

STRESS INTENSITY FACTOR OF AN ANTI-PLANE CRACK PARALLEL TO THE WEAK/MICRO-DISCONTINUOUS INTERFACE IN A BI-FGM COMPOSITE

The problem considered is a mode Ⅲ crack lying parallel to the interface of an exponential-type functional graded material （FGM） strip bonded to a linear-type FGM substrate with infinite thickness. By applying the Fourier integral transform, the problem was reduced as a Cauchy singular integral equation with an unknown dislocation density function. The collocation method based on Chebyshev polynomials proposed by Erdogan and Gupta was used to solve the singular integral equation numerically. With the numerical solution, the effects of the geometrical and physical parameters on the stress intensity factor （SIF） were analyzed and the following conclusions were drawn： （a） The region affected by the interface or free surface varies with the material rigidity, and higher material rigidity will lead to bigger affected region. （b） The SIF of the crack in the affected region and parallel to the micro-discontinuous interface is lower than those of the weak discontinuous cases. Reducing the weak-discontinuity of the interface will be beneficial to decrease the SIF of the interface-parallel crack in the region affected by the interface. （c） The effect of the free surface on SIF is more remarkable than that of the interface, and the latter is still more notable than that of the material rigidity. When the effects of the interface and free surface are fixed, increase of the material rigidity will enhance the value of SIF.     

New type of dual bem and Green's-function-library strategy for fracture analysis in complex structures

A new type of dual boundary integral equations (DBIE) is presented first, through which, a smaller system of equations needs to be solved in fracture analysis. Then a non-conforming crack tip element in two-dimensional problems is proposed. The exact formula for the hypersingular integral over the non-conforming crack tip element is given next. By virtue of Green's-function-library strategy, a series of stress intensity factors (SIF) of different crack orientations, locations and/or sizes in a complicated structure can be obtained easily and efficiently. Finally, several examples of fracture analysis in two dimensions are given to demonstrate the accuracy and efficiency of the method proposed. Partially supported by the Aeronautical Science Foundation of China (No. 99C53026)     

BOUNDARY INTEGRAL FORMULA OF ELASTIC PROBLEMS IN CIRCLE PLANE

IntroductionConcerningtheelasticplaneprobleminaunitcircle ,ZhengShenzhouandZhengXueliangdevelopedaboundaryintegralformulaofthestressfunction[1]:Φ(r,θ) =-( 1 -r2 ) 24π ∫2π0ν( φ)1 -2rcos(θ-φ) r2 dφ　　 12π∫2π011 -2rcos(θ-ω) r2 dω∫2π0μ( φ)1 -cos(ω-φ) dφ　　 1 -r22π∫2π0μ( φ)1 -2rcos(θ -φ) r2 dφ　　 ( 0 ≤r <1 ) ,( 1 )whereμ(θ) =Φ(r,θ) ｜r=1,ν(θ) = Φ n r=1= Φ r r=1.Intheformula ( 1 )theseconditemisastrongsingularintegral,itshouldbeunderstoodasanintegra…     

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）.     

The Problem of an External Circular Crack Under Asymmetric Loadings

ＩｎｔｒｏｄｕｃｔｉｏｎＷｈｅｎｔｈｅｂｏｕｎｄａｒｙｉｎｔｅｇｒａｌｅｑｕａｔｉｏｎｍｅｔｈｏｄｉｓａｐｐｌｉｅｄｔｏｃｒａｃｋａｎａｌｙｓｉｓ,ｏｎｌｙｎｕｍｅｒｉｃａｌｓｏｌｕｔｉｏｎｓｃａｎｂｅｏｂｔａｉｎｅｄ ,ｓｕｃｈａｓ:ｔｈｅｔｙｐｉｃａｌｗｏｒｋｓｏｆＳｎｙｄｅｒａｎｄＣｒｕｓｅ[1],Ｃｒｏｕｃｈ[2 ],Ｂｌａｎｄｆｏｒｄｅｔａｌ.[3],Ｐｏｒｔｅｌａｅｔａｌ.[4 ],Ｂｕｉ[5 ],Ｗｅａｖｅｒ[6 ]ａｎｄＷＡＮＧｅｔａｌ.[7- 9].Ｉｔｉｓｖｅｒｙｄｉｆｆｉｃｕｌｔｔｏａｐｐｌｙｔｈｅｂｏｕｎｄａｒ…     

奇性校正特解场法计算任意点应力和位移

本文将边界应力积分方程中的强奇性边界积分化为Ｃａｕｃｈｙ主值积分，证明了该主值积分的存在性，并给出它在变量替换中的比尺附加项显式；运用刚体位移法和单位应力场法将强奇性主子块和应力奇性系数矩阵表为同行副子块的线性组合，再配用极坐标变换和适当的数值求积计算各副子块。本法无论对平面问题还是空间问题，是光滑边界点、侧棱点还是角点，有限域或无限域，受载还是变温，都一概适用，而且可以扩展用来计算近边界点，公式统一，程序通用，数值效果良好。     

穿过反平面圆形夹杂界面的曲线型刚性线

利用复变函数方法和叠加原理建立了求解刚性线夹杂问题的弱奇积分方程,利用Ｃａｕｃｈｙ型奇异积分方程主部分方法,研究了穿过反平面圆夹杂界面的曲线型刚性线夹杂在界面交点处点处的奇性应力指数以及交点处角形域内的奇性应力,并定义了交点处的应力奇性因子。利用所得的奇性应力指数,通过对弱奇异积分方程的数值求解,得出了刚性线端点和交点处的应力奇性因子。     

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.     

数值流形方法的粘性边界问题初探

在实际工程数值流形方法分析中,采用固定约束边界的方法处理无限域或者半无限域的情况,边界处应力波的反射造成模拟结果与实际情况不符.本文基于Lysmer等人提出的粘性边界理论,在边界上设置阻尼器,推导相应粘性边界条件下流形单元刚度矩阵的数值计算格式,经岩石长条中弹性波传播算例,并与有限元结果对比,验证了该粘性边界的有效性,有利于数值流形方法的工程中推广应用.     

Interaction of elastic waves with a periodic array of collinear inplane cracks

A boundary integral equation method is applied to the study of the interaction of plane elastic waves with a periodic array of collinear inplane cracks. Numerical results are presented for the dynamic stress intensity factors. The effects of the wave type, wave frequency, wave incidence angle, and crack spacing on the dynamic stress intensity factors are analyzed in detail. The project supported by the Committee of Science and Technology of Shanghai and Tongji University     

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

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

NEW BOUNDARY ELEMENT METHOD FOR TORSION PROBLEMS OF CYLINDER WITH CURVILINEAR CRACKS

The Saint-Venant torsion problems of a cylinder with curvilinear cracks were considered and reduced to solving the boundary integral equations only on cracks. Using the interpolation models for both singular crack tip elements and other crack linear elements, the boundary element formulas of the torsion rigidity and stress intensity factors were given. Some typical torsion problems of a cylinder involving a straight, kinked or curvilinear crack were calculated. The obtained results for the case of straight crack agree well with those given by using the Gauss-Chebyshev integration formulas, which demonstrates the validity and applicability of the present boundary element method.     

三维裂纹分析的主值型面力边界积分方程法

给出了一组只包含Ｃａｕｃｈｙ主值积分、不含有强奇异积分的三维静动力边界积分方程及其应用于裂纹问题的具体列式，并给出了几何轴对称问题的相应半解析边界元求解方法，将三维问题降阶为一维数值问题．文中分析了无限、半无限介质中圆裂纹、平行圆裂纹系、球面裂纹等在静载及应力波作用下的静力或瞬态动力响应问题，求得了相应的应力强度因子．     

Elastic-plastic analysis of antiplane crack near crack surface region

The elastic-plastic stress distribution and the elastic-plastic boundary con- figuration near a crack surface region are significant but hard to obtain by means of the conventional analysis. A crack line analysis method is developed in this paper by consid- ering the crack surface as an extension of the crack line. The stresses in the plastic zone, the length, and the unit normal vector of the elastic-plastic boundary near a crack surface region are obtained for an antiplane crack in an elastic-perfectly plastic solid. The usual small scale yielding assumptions are not needed in the analysis.