A boundary integral equation method for three-dimensional crack problems in elasticity

This paper presents a solution procedure for three-dimensional crack problems via first kind boundary integral equations on the crack surface. The Dirichlet (Neumann) problem is reduced to a system of integral equations for the jump of the traction (of the field) across the crack surface. The calculus of pseudodifferential operators is used to derive existence and regularity of the solutions of the integral equations. With the concept of the principal symbol and the Wiener-Hopf technique we derive the explicit behavior of the densities of the integral equations near the edge of the crack surface. Based on the detailed regularity results we show how to improve the boundary element Galerkin method for our integral equations. Quasi-optimal asymptotic estimates for the Galerkin error are given.     

Analysis of arbitrarily shaped planar cracks in two-dimensional hexagonal quasicrystals with thermal effects. Part I: Theoretical solutions

An extended displacement discontinuity (EDD) boundary integral equation method is proposed for analysis of arbitrarily shaped planar cracks in two-dimensional (2D) hexagonal quasicrystals (QCs) with thermal effects. The EDDs include the phonon and phason displacement discontinuities and the temperature discontinuity on the crack surface. Green's functions for unit point EDDs in an infinite three-dimensional medium of 2D hexagonal QC are derived using the Hankel transform method. Based on the Green's functions and the superposition theorem, the EDD boundary integral equations for an arbitrarily shaped planar crack in an infinite 2D hexagonal QC body are established. Using the EDD boundary integral equation method, the asymptotic behavior along the crack front is studied and the classical singular index of 1/2 is obtained at the crack edge. The extended stress intensity factors are expressed in terms of the EDDs across crack surfaces. Finally, the energy release rate is obtained using the definitions of the stress intensity factors.     

Transient response of dissimilar piezoelectric layers with multiple interacting interface cracks

In this study, we examine the dynamic behavior of two bonded dissimilar piezoelectric layers containing multiple interfacial cracks subjected to electro-mechanical impact loading. The problem was formulated through Fourier transformation into singular integral equations in which the unknown variables are the jumps of displacement and electric potential across the crack surface in the Laplace transform domain. The resulting integral equations together with the corresponding single-valued conditions are solved numerically for the densities of electro-elastic dislocations on a crack surface. The dynamic field intensity factors and dynamic energy release rate (DERR) history are obtained for both permeable and impermeable crack. The stress field is also obtained for the interface crack under impact loads. The results show that the field intensity factors at the crack tips and dynamic energy release rate depend on the interfacial crack geometry, electromechanical coupling and the electric boundary conditions on the crack surface.     

求解平片裂纹问题的有限部积分与边界元法

本文利用位移的Somigliana公式和有限部积分的概念,导出了求解三维弹性力学中的任意形状平片裂纹问题的超奇异积分方程组,进而联合使用有限部积分法与边界元法对所得方程建立了数值法.为验证本文的方法,计算了若干数值例子的裂纹面的位移间断及裂纹前沿的应力强度因子,它们与理论值相比符合很好.     

A subregion DRBEM formulation for the dynamic analysis of two-dimensional cracks

The dual reciprocity boundary element method employing the step by step time integration technique is developed to analyse two-dimensional dynamic crack problems. In this method the equation of motion is expressed in boundary integral form using elastostatic fundamental solutions. In order to transform the domain integral into an equivalent boundary integral, a general radial basis function is used for the derivation of the particular solutions. The dual reciprocity boundary element method is combined with an efficient subregion boundary element method to overcome the difficulty of a singular system of algebraic equations in crack problems. Dynamic stress intensity factors are calculated using the discontinuous quarter-point elements. Several examples are presented to show the formulation details and to demonstrate the computational efficiency of the method.     

Boundary integral equation method in micropolar elasticity

This paper presents the integral representations of the displacement and rotation fields of the micropolar continuum, and the regularized integral representations of the traction and couple vector. Thus the boundary integral equations and the boundary integro-differential equations are derived. In addition, the formulation for crack problems is given by both the boundary integral equations and the boundary integro-differential equations.     

压电材料中两个非对称平行裂纹的基本解

采用Schmidt方法分析压电材料中非对称平行的双可导通裂纹的断裂性能．利用Fourier变换使问题的求解转换为求解两对以裂纹面位移之差为未知变量的对偶积分方程．为了求解对偶积分方程,直接把裂纹面位移差函数展开成Jacobi多项式形式．最终得到了裂纹的应力强度因子与电位移强度因子之间的关系．数值结果表明,应力强度因子和电位移强度因子与裂纹间的距离、裂纹的几何尺寸有关;与不可导通裂纹有关结果相比,可导通裂纹的电位移强度因子远小于相应问题不可导通裂纹的电位移强度因子．同时可以发现裂纹间的“屏蔽”效应也在压电材料中出现．     

非对称载荷作用的外部圆形裂纹问题

使用边界积分方程方法,研究了三维无限弹性体中受非对称载荷作用的外部圆形裂纹问题。通过使用Fourier级数和超几何函数,将问题的二维边界奇异积分方程简化为Abel型方程,获得了一般非对称载荷作用的外部圆形裂纹问题的应力强度因子精确解,比用Hankel变换法得到的结果更为一般。结果表明:边界积分方程法在解析分析方面还有很大的潜力。     

An application of a boundary element method to natural convection

The direct boundary element method is applied to the numerical modelling of thermal fluid flow in a transient state. The Navier-Stokes equations are considered under the Boussinesq approximation and the viscous thermal flow equations are expressed in terms of stream function, vorticity, and temperature in two dimensions. Boundary integral equations are derived using logarithmic potential and time-dependent heat potential as fundamental solutions. Boundary unknowns are discretized by linear boundary elements and flow domains are divided into a series of triangular cells. Charged points are translated upstream in the numerical evaluation of convective terms. Unknown stream function, vorticity, and temperature are staggered in the computational scheme.

Simple iteration is found to converge to the quasi steady-state flow. Boundary solutions for two-dimensional examples at a Reynolds number 100 and Grashoff number 10⁷ are obtained.     

Analysis solution method for 3D planar crack problems of two-dimensional hexagonal quasicrystals with thermal effects

An analysis solution method (ASM) is proposed for analyzing arbitrarily shaped planar cracks in two-dimensional (2D) hexagonal quasicrystal (QC) media. The extended displacement discontinuity (EDD) boundary integral equations governing three-dimensional (3D) crack problems are transferred to simplified integral-differential forms by introducing some complex quantities. The proposed ASM is based on the analogy between these EDD boundary equations for 3D planar cracks problems of 2D hexagonal QCs and those in isotropic thermoelastic materials. Mixed model crack problems under combined normal, tangential and thermal loadings are considered in 2D hexagonal QC media. By virtue of ASM, the solutions to 3D planar crack problems under various types of loadings for 2D hexagonal QCs are formulated through comparison to the corresponding solutions of isotropic thermoelastic materials which have been studied intensively and extensively. As an application, analytical solutions of a penny-shaped crack subjected uniform distributed combined loadings are obtained. Especially, the analytical solutions to a penny-shaped crack subjected to the anti-symmetric uniform thermal loading are first derived for 2D hexagonal QCs. Numerical solutions obtained by EDD boundary element method provide a way to verify the validity of the presented formulation. The influences of phonon-phason coupling effect on fracture parameters of 2D hexagonal QCs are assessed.     

Effect of shearing forces on the stressed state of a half space with crack

Using the method of boundary integral equations, we study the stressed state in the neighborhood of a plane crack perpendicular to the boundary of a half space. The crack surfaces are subjected to the action of shearing forces. The problem is reduced to two-dimensional hypersingular integral equations, and their regular kernels, taking into account interaction between the crack and boundary of the half space, are written in explicit form. The dependences of stress intensity factors on the angular coordinate are presented for different loads of the crack. Translated from Matematychni Metody ta Fizyko-Mekhanichni Polya, Vol. 51, No. 1, pp. 112–120, January–March, 2008.     

An extended thermal-medium crack model

It is important to investigate the effects of heat conduction of crack interior on thermoelastic fields of a cracked material. In this paper, an extended thermal-medium crack model is proposed to address the influences of the thermal conductivity inside an opening crack on the induced thermoelastic fields. Then the problem of a penny-shaped crack in a transversely isotropic material is investigated under applied mechanical and uniform heat flow loadings. Based on the Hankel transform technique, the governing partial differential equations are transformed to ordinary differential equations, then to a system of coupled dual integral equations. The thermoelastic fields around the penny-shaped crack are obtained explicitly by solving the derived dual integral equations. Numerical results are reported to show the influences of the thermal conductivity of crack interior on partial insulation coefficient, temperature change across crack and thermal stress intensity factor. As compared to the known thermal-medium crack model, the proposed one exhibits more applicability.     

A NUMERICAL BOUNDARY EIGENVALUE PROBLEM FOR ELASTIC CRACKS IN FREE AND HALF SPACE

We present in this paper a numerical method for hypersingular boundary integral equations.This method was developed for planar crack problems:additional edge si...     

Elastodynamic Problem of an Interface Linear Crack under Harmonic Loading

The present work is devoted to application of boundary integral equations to the 2D problem for a linear crack located on the bimaterial interface under harmonic loading. The system of linear algebraic equations is derived to solve the problem numerically. The distribution of the displacements and tractions at the bonding interface and the surface of the crack are obtained for the case of the tension–compression wave which propagates normally to the interface. The results are compared with those obtained for the static case. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

An interface crack moving between magnetoelectroelastic and functionally graded elastic layers

A constant crack moving along the interface of magnetoelectroelastic and functionally graded elastic layers under anti-plane shear and in-plane electric and magnetic loading is investigated by the integral transform method. Fourier transforms are applied to reduce the mixed boundary value problem of the crack to dual integral equations, which are expressed in terms of Fredholm integral equations of the second kind. The singular stress, electric displacement and magnetic induction near the crack tip are obtained asymptotically and the corresponding field intensity factors are defined. Numerical results show that the stress intensity factors are influenced by the crack moving velocity, the material properties, the functionally graded parameter and the geometric size ratios. The propagation of the moving crack may bring about crack kinking, depending on the crack moving velocity and the material properties across the interface.     

三维横观各向同性介质界面裂纹的边界积分方程方法

基于两相三维横观各向同性介质的基本解和Somigliana恒等式,对三维横观各向同性介质中的任意形状的平片界面裂纹,以裂纹面上的不连续位移为待求参量建立了超奇异积分_微分方程,界面平行于横观各向同性面.根据发散积分的有限部积分理论,应用积分方程方法研究得到裂纹前沿的位移和应力场的表达式、奇性指数以及应力强度因子的不连续位移表达式.在非震荡情形下,超奇异积分_微分方程退化为超奇异积分方程,与均匀介质的超奇异积分方程形式完全相同.     

Special Green’s function boundary element approach for steady-state axisymmetric heat conduction across low and high conducting planar interfaces

The problem of determining the steady-state axisymmetric temperature distribution in a bimaterial with a planar interface is considered here. The interface is either low or high conducting. Special Green’s functions satisfying the thermal conditions on the interface are derived and employed to obtain boundary integral equations whose path of integration does not include the interface. Boundary element procedures that do not require the interface to be discretized into elements are proposed for solving the problem under consideration.     

热载荷作用下平面裂纹问题的奇异积分方程与边界无法

从边界积分方程出发,导出了二维裂纹体热传导问题及热弹性问题的积分方程组,继而使用奇异积分方程与边界元相结合的方法,为其建立了相应的数值求解方法。此外,利用奇异积分方程的主部分析法,严格地证明了裂纹尖端温度梯度场的1/√r 奇异性,并且给出了奇性温度梯度场的精确解。最后。对一些典型例子,做了数值计算。     

An augmented galerkin procedure for the boundary integral method applied to two-dimensional screen and crack problems

Here we present a new solution procedure for Helm-holtz and Laplacian Dirichlet screen and crack problems in IR² via boundary integral equations of the first kind having as an unknown the jump of the normal derivative across the screen or a crack curve T. Under the assumption of local finite energy we show the equivalence of the integral equations and the original boundary value problem. Via the method of local Mellin transform in [5]-[lo] and the calculus of pseudodifferential operators we derive existence, uniqueness and regularity results for the solution of our boundary integral equations together with its explicit behaviour near the screen or crack tips.With our integral equations we set up a Galerkin scheme on T and obtain high quasi-optimal convergence rates by using special singular elements besides regular splines as test and trial functions.     

压电压磁复合材料中界面裂纹对弹性波的散射

利用Schmidt方法分析了压电压磁复合材料中可导通界面裂纹对反平面简谐波的散射问题．经过富里叶变换得到了以裂纹面上的间断位移为未知变量的对偶积分方程A·D2在求解对偶积分方程的过程中,裂纹面上的间断位移被展开成雅可比多项式的形式．数值模拟分析了裂纹长度、波速和入射波频率对应力强度因子、电位移强度因子、磁通量强度因子的影响A·D2从结果中可以看出,压电压磁复合材料中可导通界面裂纹的反平面问题的应力奇异性形式与一般弹性材料中的反平面问题应力奇异性形式相同．