Singular integral equations and boundary element method of cracks in thermally stressed planar solids

Ｉｔｉｓｋｎｏｗｎｔｈａｔｍｏｓｔａｇｒｉｃｕｌｔｕｒａｌｐｒｏｄｕｃｔｓａｎｄｆｏｏｄｓａｒｅｐｒｏｃｅｓｓｅｄａｎｄｔｒａｎｓｐｏｒｔｅｄｕｎｄｅｒｃｅｒｔａｉｎｔｅｍｐｅｒａｔｕｒｅｃｏｎｄｉｔｉｏｎｓ,ａｎｄｔｈｅｓｔｒｕｃｔｕｒａｌｃｏｍｐｏｎｅｎｔｓａｌｓｏｗｏｒｋｕｎｄｅｒａｔｈｅｒｍａｌｅｎｖｉｒｏｎｍｅｎｔ.Ｔｅｍｐｅｒａｔｕｒｅｉｎｄｕｃｅｄｓｔｒｅｓｓｅｓｕｓｕａｌｌｙｌｅａｄｔｏｄａｍａｇｅｏｆｆｌａｗｅｄｓｏｌｉｄｓ.Ｔｈｕｓ,ｔｈｅｉｎｖｅｓｔｉｇａｔｉｏｎｏｆｔｈｅｃｒ…     

一种处理弹性动力边界元法中奇异积分的方法

利用边界元法求解瞬态弹性动力学问题时,时域基本解函数的分段连续性和奇异性为该问题的求解带来很大的困难。为了解决时域基本解中的奇异性问题,本文依据柯西主值的定义,对经过时间解析积分之后的时域基本解进行奇异值分解,将其分成奇异和正则积分两部分;其中正则部分可通过采用常规高斯积分方法来计算,而奇异部分具有简单的形式,可以利用解析积分计算。经过上述操作之后,就可以达到直接消除时域基本解中奇异积分的目的。和传统方法相比,本文方法并不依赖静力学基本解来消除奇异性,是一种直接求解方法。最后给定两个数值算例来验证本文提出方法的正确性和可行性,结果表明使用本文算法可以解决弹性动力学边界积分方程中的奇异性问题。     

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

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

势问题的重构核粒子边界无单元法

将重构核粒子法和势问题的边界积分方程方法结合,提出了势问题的重构核粒子边界无单元法. 推导了势问题的重构核粒子边界无单元法的公式,研究其数值积分方案,建立了重构核粒子边界无单元法的离散化边界积分方程,并推导了重构核粒子边界无单元法的内点位势的积分公式. 重构核粒子法形成的形函数具有重构核函数的光滑性,且能再现多项式在插值点的精确值,所以该方法具有更高的精度. 最后给出了数值算例,验证了所提方法的有效性和正确性. }      

Helmholtz边界积分方程中奇异积分间接求解方法

提出了间接求解传统Helmholtz边界积分方程CBIE的强奇异积分和自由项系数,以及Burton-Miller边界积分方程BMBIE中的超强奇异积分的特解法。对于声场的内域问题,给出了满足Helmholtz控制方程的特解,间接求出了CBIE中的强奇异积分和自由项系数。对于声场外域对应的BMBIE中的超强奇异积分,按Guiggiani方法计算其柯西主值积分需要进行泰勒级数展开的高阶近似,公式繁复,实施困难。本文给出了满足Helmholtz控制方程和Sommerfeld散射条件的特解,提出了间接求出超强奇异积分的方法。推导了轴对称结构外场问题的强奇异积分中的柯西主值积分表达式,并通过轴对称问题算例证明了本文方法的高效性。数值结果表明,对于内域问题,采用本文特解法的计算结果优于直接求解强奇异积分和自由项系数的结果,且本文的特解法可避免针对具体几何信息计算自由项系数,因而具有更好的适用性。对于外域问题,两者精度相当,但本文的特解法可避免对核函数进行高阶泰勒级数展开,更易于数值实施。     

用样条边界元与奇异积分方程法联合求解多裂纹柱体的扭转

通过间解的分离，本文将径向多裂纹柱体的导曲函两个调和函数表示，使问题归为解一组混混合型积分方程。针对方程的特点，本文联合使用三次样条边界法与奇异积分方程的数值方法对所得方程建立了数值法，并对裂纹相交情形作了特殊处理。最后对工程中感兴趣的一些典型的多裂纹柱体的扭转作了例题计算，结果表明，本文方法具有收敛快，精度高的特点。     

A CLASS OF NONLINEAR BOUNDARY VALUE PROBLEMS FOR THE SECOND-ORDER E2 CLASS ELLIPTIC SYSTEMS IN GENERAL FORM

IntroductionThispaperdiscussesaclassofnonlinearboundaryvalueproblems (BVP)forthesecond_orderE2 classellipticsystemsingeneralform ,whichareofgreatimportanceinmathematicalsciencesandengineering ,andattractmoreandmoreinterestsfrommanyresearchersdomesticandabroad .ThesenonlinearBVPwerefirstintroducedbyRussianresearcher,L .V .Wolfersdorf,whereaclassofnonlinearBVPforholomorphicfunctionswasstudiedandsomeefficientresultsderived[1,2 ].LaterSigridKenchtandAliSeifMashimba[3],R .P .Gilbertandthea…     

Solution of the planar Newtonian stick–slip problem with the singular function boundary integral method

A singular function boundary integral method (SFBIM) is proposed for solving biharmonic problems with boundary singularities. The method is applied to the Newtonian stick–slip flow problem. The streamfunction is approximated by the leading terms of the local asymptotic solution expansion which are also used to weight the governing biharmonic equation in the Galerkin sense. By means of the divergence theorem the discretized equations are reduced to boundary integrals. The Dirichlet boundary conditions are weakly enforced by means of Lagrange multipliers, the values of which are calculated together with the singular coefficients. The method converges very fast with the number of singular functions and the number of Lagrange multipliers, and accurate estimates of the leading singular coefficients are obtained. Comparisons with the analytical solution and results obtained with other numerical methods are also made. Copyright © 2005 John Wiley & Sons, Ltd.     

半平面多边缘裂纹反平面问题的奇异积分方程

利用复变函数和奇异积分方程方法,求解弹性范围内半平面多边缘裂纹的反平面问题．提出了满足半平面边界自由的由分布位错密度表示的单边缘裂纹的基本解,此基本解由主要部分和辅助部分组成．将半平面多边缘裂纹问题看作是许多单边缘裂纹问题的叠加,建立了一组Cauchy型奇异积分方程．然后,利用半开型积分法则求解该奇异积分方程,得到了裂纹端处的应力强度因子．最后,给出了几个数值算例．     

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

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.     

求解位势问题的虚边界元法

本文提出了求解位势问题的虚边界元法，建立了位势问题的虚边界元的离散方程式，推导了离散化求系数的积分解析式。该方法与传统边界元法相比具有不存在奇异积分和边界附近精度较高等优点，可用来计算真空静电场、稳定温度场、流体绕流、介质中的渗流等各类位势问题。大量算例均获得了满意的结果。     

Numerical solution of singular integral equations in stress concentration problems under longitudinal shear loading

Summary The paper deals with numerical solutions of singular integral equations in stress concentration problems for longitudinal shear loading. The body force method is used to formulate the problem as a system of singular integral equations with Cauchy-type singularities, where unknown functions are densities of body forces distributed in the longitudinal direction of an infinite body. First, four kinds of fundamental density functions are introduced to satisfy completely the boundary conditions for an elliptical boundary in the range 0≤φ _k ≤2π. To explain the idea of the fundamental densities, four kinds of equivalent auxiliary body force densities are defined in the range 0≤φ _k ≤π/2, and necessary conditions that the densities must satisfy are described. Then, four kinds of fundamental density functions are explained as sample functions to satisfy the necessary conditions. Next, the unknown functions of the body force densities are approximated by a linear combination of the fundamental density functions and weight functions, which are unknown. Calculations are carried out for several arrangements of elliptical holes. It is found that the present method yields rapidly converging numerical results. The body force densities and stress distributions along the boundaries are shown in figures to demonstrate the accuracy of the present solutions. Received 26 May 1998; accepted for publication 27 November 1998     

各向异性平面含斜裂纹的奇异积分方程方法

本文应用平面弹性复变方法，将无限各向异性平面中的任意斜裂纹问题归结为求解一组解析函数边值问题，通过构造适当的积分变换将边值问题转化为奇异积分方程，进而应用Lobotto-Chebyshev数值求积公式，求出该奇异积分方程的数值解，并得到了应力强度因子的近似表达式，最后，给出了一些实例的数值结果，对特例的数值结果与精确结果进行比较，吻合的很好。     

薄体位势问题边界元法中的解析积分算法

薄体结构的数值分析是边界元法的难点问题之一。该文导出了一种完全解析积分算法，用这种算法计算了薄体平面位势问题边界元法中出现的几乎弱奇异、强奇异和超奇异积分。当边界离散为一系列线性单元，边界积分方程离散计算的积分可归纳为三种形式。对薄体问题，源点与积分单元距离通常相距很近，这些积分产生显著几乎奇异性，直接采用常规高斯积分不能有效计算。为此该文导出了这些几乎奇异积分的全解析计算公式。按源点与单元的距离是否为零，公式分两种情况。新算法采用全解析积分公式处理几乎奇异积分，首先精确计算出薄体问题边界未知位势和法向位势梯度，然后再进一步计算了域内点的物理参量。算例表明该文算法可处理狭长比为1．E-08的薄体问题，显示了边界元法分析薄体问题具有独特的优势。     

Isotropic/orthotropic two-layered strips containing cracks terminating at and going through an interface

Crack problems for isotropic/orthotropic two-layered strips have been investigated. A system of two singular integral equations can be derived by using Fourier integral transformation and boundary conditions of crack problems. After stress singularities at crack tips or other special points are determined for internal and edge cracks, and for cracks terminating at and going through the interface, the system of singular integral equations is solved numerically by Gauss-Jacobi or Gauss-Chebyshev integration formulas for stress intensity factors at the tips and other singular points of cracks. Finally, possible crack growth behavior for cracks approaching and going through the interface is discussed.     

A class of nonlinear boundary value problems for the second-order<Emphasis Type="Italic">E</Emphasis><Subscript>2</Subscript> class elliptic systems in general form

A class of nonlinear boundary value problems (BVP) for the second-order E₂ class elliptic systems in general form is discussed. By introducing a kind of transformation, this kind of BVP is reduced to a class of generalized nonlinear Riemann-Hilbert BVP. And then some singular integral operators are introduced to establish the equivalent nonlinear singular integral equations. The solvability is proved under some suitable hypotheses by means of the properties of singular integral operators and the function theoretic methods. Foundation items: the National Natural Science Foundation of China (19671056); Shanghai Municipal Natural Scientific Foundation (99ZA14030, 01ZA14023); Jiangxi Provincial Natural Scientific Foundation (981102, 0211014) Biographies: LI Ming-zhong (1935−); XU Ding-hua (1967−)     

ELASTIC ANALYSIS OF ORTHOTROPIC PLANE PROBLEMS BY THE SPLINE FICTITIOUS BOUNDARY ELEMENT METHOD

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｔｈｅｃｏｎｖｅｎｔｉｏｎａｌｂｏｕｎｄａｒｙｅｌｅｍｅｎｔｍｅｔｈｏｄ (ＢＥＭ) ,ｓｉｎｇｕｌａｒｐｏｉｎｔｓａｒｅｄｉｓｔｒｉｂｕｔｅｄａｌｏｎｇｔｈｅｂｏｕｎｄａｒｙｏｆｔｈｅｄｏｍａｉｎｕｎｄｅｒｓｔｕｄｙ ,ｗｈｉｃｈｌｅａｄｓｔｏｓｉｎｇｕｌａｒｉｎｔｅｇｒａｌｅｑｕａｔｉｏｎｓｄｕｅｔｏｔｈｅｓｉｎｇｕｌａｒｉｔｙｏｆｆｕｎｄａｍｅｎｔａｌｓｏｌｕｔｉｏｎｓ .Ｔｈｅｒｅｆｏｒｅ ,ｓｉｎｇｕｌａｒｉｎｔｅｇｒａｔｉｏｎｍｕｓｔｂｅｈａｎｄｌｅｄｗｈｅｎｓｏｌ…     

三维变系数热传导问题边界元分析中几乎奇异积分计算

在边界积分的数值计算过程中,当源点离积分单元很近时,边界积分就会具有几乎奇异性,此时不能直接用高斯数值积分公式计算几乎奇异积分。本文以三维非均质热传导问题为例,介绍了一种计算几乎奇异边界积分的新方法。首先,采用Newton-Raphson迭代算法确定积分单元上离源点最近的点;然后,将积分单元上任意一点的坐标在最近点处展开成泰勒级数,并计算源点到积分单元任意点的距离;最后,将距离函数代入几乎奇异边界积分中,并运用指数变换方法导出积分单元上几乎奇异积分的计算公式。文中给出了两个非均质热传导问题的算例来验证所述方法的正确性、有效性和稳定性。     

黏弹性体界面裂纹的冲击响应

研究两半无限大黏弹性体界面Griffith裂纹在反平面剪切突出载荷下,裂纹尖端动应力强度因子的时间响应,首先,运用积分变换方法将黏弹性混合黑社会问题化成变换域上的对偶积分方程,通过引入裂纹位错密度函数进一步化成Cauchy型奇异积分方程,运用分片连续函数法数值求解奇异积分方程,得到变换域内的动应力强度因子,再用Laplace积分变换数值反演方法,将变换域的解反演到时间域内,最终求得动应力强度因子的时间响应,并对黏弹性参数的影响进行分析。