Numerical solutions of singular integral equations for planar rectangular interfacial crack in three dimensional bimaterials

Stress intensity factors for a three dimensional rectangular interfacial crack were considered using the body force method. In the numerical calculations, unknown body force densities were approximated by the products of the fundamental densities and power series; here the fundamental densities are chosen to express singular stress fields due to an interface crack exactly. The calculation shows that the numerical results are satisfied. The stress intensity factors for a rectangular interface crack were indicated accurately with the varying aspect ratio, and bimaterial parameter.     

Numerical methods for solving singular integral equations obtained by fracture mechanical analysis of cracked wedge简

The numerical solutions to the singular integral equations obtained by the fracture mechanical analyses of a cracked wedge under three different conditions are considered. The three considered conditions are: (i) a radial crack on a wedge with a non-finite radius under the traction-traction boundary condition, (ii) a radial crack on a wedge with a finite radius under the traction-traction boundary condition, and (iii) a radial crack on a finite radius wedge under the traction-displacement boundary condition. According to the boundary conditions, the extracted singular integral equations have different forms. Numerical methods are used to solve the obtained coupled singular integral equations, where the Gauss-Legendre and the Gauss-Chebyshev polynomials are used to approximate the responses of the singular integral equations. The results are presented in figures and compared with those obtained by the analytical response. The results show that the obtained Gauss-Chebyshev polynomial response is closer to the analytical response.     

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.     

Boundary integral equations of unique solutions in elasticity

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｂｏｕｎｄａｒｙｅｌｅｍｅｎｔｍｅｔｈｏｄ（ＢＥＭ）ｐｒｏｖｉｄｅｓａｎａｔｔｒａｃｔｉｖｅａｌｔｅｒｎａｔｉｖｅｆｏｒｔｈｅａｎａｌｙｓｉｓｏｆｅｎｇｉｎｅｅｒｉｎｇｐｒｏｂｌｅｍｓ．Ｉｔｓｍａｉｎａｄｖａｎｔａｇｅｓａｒｅｅｃｏｎｏｍｉｃａｌａｎｄｐａｒｔｉｃｕｌａｒｌｙｃｏｎｖｅｎｉｅｎｔｆｏｒｕｎｂｏｕｎｄｅｄｄｏｍａｉｎａｎｄｓｔｒｅｓｓｃｏｎｃｅｎｔｒａｔｉｏｎｐｒｏｂｌｅｍｓ．Ｔｈｅｂｏｕｎｄａｒｙｉｎｔｅｇｒａｌｅｑｕａｔｉｏｎ（ＢＩＥ）ｉｓｔｈｅ…     

Interaction between elliptical and ellipsoidal inclusions under bending stress fields

Summary  This paper deals with interaction problems of elliptical and ellipsoidal inclusions under bending, using singular integral equations of the body force method. The problems are formulated as a system of singular integral equations with Cauchy-type or logarithmic-type singularities, where unknown functions are densities of body forces distributed in the x,y and r,θ,z directions in infinite bodies having the same elastic constants as those of the matrix and inclusions. In order to satisfy the boundary conditions along the elliptical and the ellipsoidal boundaries, the unknown functions are approximated by a linear combination of fundamental density functions and polynomials. The present method is found to yield the exact solutions for a single elliptical or spherical inclusion under a bending stress field. It yields rapidly converging numerical results for interface stresses in the interaction of inclusions. Received 9 September 1999; accepted for publication 15 January 2000     

Numerical solution of multiple hole problem by using boundary integral equation

This paper studies a numerical solution of multiple hole problem by using a boundary integral equation. The studied problem can be considered as a supposition of many single hole problems. After considering the interaction among holes, an algebraic equation is formulated, which is then solved by using an iteration technique. The hoop stress around holes can be finally determined. One numerical example is provided to check its accuracy.     

The dynamic stress intensity factor analysis of adhesively bonded material interface crack with damage under shear loading

This paper studies the dynamic stress intensity factor (DSIF) at the interface in an adhesive joint under shear loading. Material damage is considered. By introducing the dislocation density function and using the integral transform, the problem is reduced to algebraic equations and can be solved with the collocation dots method in the Laplace domain. Time response of DSIF is calculated with the inverse Laplace integral transform. The results show that the mode Ⅱ DSIF increases with the shear relaxation parameter, shear module and Poisson ratio, while decreases with the swell relaxation parameter. Damage shielding only occurs at the initial stage of crack propagation. The singular index of crack tip is -0.5 and independent on the material parameters, damage conditions of materials, and time. The oscillatory index is controlled by viscoelastic material parameters.     

Uniquenes of solution of field point of singular source outside-region-distribution method

The uniqueness of solution of field point, inside a convex region due to singular source(s) with kernel function decreasing with distance increasing, outsidergeion-distribution(s) such that the boundary condition expressed by the response of the source(s) is satisfied, is proved by using the condition of kernel function decreasing with distance increasing and an integral inequality. Examples of part ofthese singular sources such as Kelvin’s point force, Point-Ring-Couple (PRC) etc.are given. The proof of uniqueness of solution of field point in a twisted shaft of revolution due to PRC distribution is given as an example of application.     

热弹性平面问题的规则化边界积分方程

对于热弹性平面问题,过去广泛集中在直接变量边界元法研究,本文研究间接变量规则化边界元法,建立了间接变量规则化边界积分方程。和直接边界元法相比,间接法具有降低密度函数的连续性要求、位移梯度方程中的热载荷体积分具有较弱奇异性等优点。数值实施中,用精确单元描述边界几何,不连续插值函数逼近边界量。算例表明,本文方法效率高,所得数值结果与精确解相当吻合。     

Numerical solution of singular perturbation problems for the fourth-order elliptic differential equations

In this paper we construct the finite-difference scheme for the singularly perturbedboundary value problem for the fourth-order elliptic differential equation on the basis ofpaper[1],and prove the uniform convergence of this scheme with respect to the smallparameterεin the discrete energy norm,Finally,we give a numerical example.     

Numerical solution of integro-differential and singular integral equations for plate bending problems

Two integral equation formulations for the determination of the vertical displacement and the bending moment around holes in an elastic plate are presented. Each formulation consists of two equations, the first one an integral equation and an integro-differential equation and the second one two singular integral equations. The equations are solved using B-splines as approximations to the unknowns and the method is applied to the case of one elliptic hole in a twisted plate.

---

Zusammenfassung Zwei verschiedene Integralgleichungssysteme für die Bestimmung von die Durchbiegung und die Biegemoment in einer gelochten elastischen Platte werden entwickelt. Die eine Systeme besteht von einer Integralgleichung und einer Integro-Differentialgleichung und die andere von zwei singulären Integralgleichungen von Cauchy'schen Typus. Bei der Auflösung der Systeme werden die unbekannten mit Hilfe B-splines ausgedrückt. Beide Systeme werden benutzt in dem Fall von einer elliptischen Loch in einer uendlichen, tordierten Platte.

---

---

An abbreviated version of this paper was included in a paper Integral equation solutions to mechanical problems. A review and an application to plate theory which will appear in the Proceedings from second national congress of theoretical and applied mechanics, Druzba, Bulgaria, 1973.     

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

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

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

An iteration method for integral equations arising from axisymmetric loading problems

Let the concentrated forces and the centers of pressure with unknown density functions x(ξ) and y(ξ) respectively be distributed along the axis z outside the solid, then one can reduce an axismmetric loading problem of solids of revolution to two simultaneous Fredholm integral equations. An iteration method for solving such equations is duscussed. A lemma equivalent to E. Rakotch’s contractive mapping theorem and a theorem concerning the convergent proof of the iteration method are presented.     

A numerical method for Cauchy principal value of singular integral under parametric transformation in BEM

A numerical method for the Cauchy principal value of the singular integral in BEM is developed using the concept of finite part integration under integral variable transformation. It is applied to the numerical integration on isoparametric element successfully, as shown in the examples in this paper.The project supported by National Natural Foundation of China.     

A unified approach to problems of stress concentration near V-shaped notches with sharp and rounded tip

The paper proposes a unified approach to problems of stress concentration near notches with sharp and rounded tip based on the method of singular integral equations. A solution for an elastic region having a V-shaped notch with rounded tip of large curvature is first found. Then, the stress intensity factor at the tip of a sharp-tipped notch is calculated by passing to the limit. Numerical results are obtained for a slit and a square hole in an elastic plane and an edge notch in a half-plane __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 2, pp. 70–87, February 2007. For the centenary of the birth of G. N. Savin.     

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

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

Physically nonlinear problems of stress concentration in members,structures, and buildings

Results from studies of stress concentration in physically nonlinear elastic structural members are reviewed __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 1, pp. 93–99, January 2007.     

The effect of transverse shear deformation on stress concentration factors for shallow shells with a small circular hole

Simplified equations are derived for the analysis of stress concentration for shear-deformable shallow shells with a small hole.General solutions of the equations are obtained,in terms of series,for shallow spherical shells and shallow circular cylindrical shells with asmall circular hole.Approximate explicit solutions and numerical results are obtianed forthe stress concentration factors of shallow circular cylindrical shells with a small hole onwhich uniform pressure is acting.     

辛解析奇异元在动荷载作用下V型切口问题中的应用

利用辛解析奇异单元,结合时域精细算法,研究了动荷载作用下的含平面V型切口问题。时域上,采用时域精细算法,并结合自适应算法控制展开项数,保证了计算精度。空间域上,切口尖端附近采用辛解析奇异单元,其余区域采用常规有限单元,避免了局部网格加密。本文使用的辛解析奇异单元不需要过渡单元和局部网格加密,且能够通过奇异单元内部的参数关系直接给出切口尖端的应力强度因子,不需要复杂的后处理过程。数值结果表明,本文方法具有良好的精度和稳定性,可以准确地计算动态应力强度因子。