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

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

平面定常Stokes问题的无奇异第一类边界积分方程

对无奇异边界积分方程归化法的研究,已有的结果都是针对直接变量的,其核心思想是利用刚体位移(包括刚体的转动和平移)或均匀场．然而,对第一类边界积分方程的无奇异边界归化法的研究,至今还未涉足．本文提交一种新方法,归化出平面定常Stokes问题的第一类无奇异边界积分方程,并建立完整的数值求解体系．一个简单的算例表明本文方法可获得理想的数值结果,特别是边界量的数值结果。     

含曲线裂纹圆柱扭转问题的新边界元法

研究含曲线裂纹圆柱的Saint-Venant扭转,将问题化归为裂纹上边界积分方程的求解．利用裂纹尖端的奇异元和线性元插值模型,给出了扭转刚度和应力强度因子的边界元计算公式．对圆弧裂纹、曲折裂纹以及直线裂纹的典型问题进行了数值计算,并与用Gauss-Chebyshev求积法计算的直裂纹情形结果进行了比较,证明了方法的有效性和正确性．     

改进的无奇异局部边界积分方程方法

在局部边界积分方程方法中,当源节点位于分析域的整体边界上时,局部边界积分将出现奇异积分问题,这些奇异积分需要做特别的处理．为此,提出了对域内节点采用局部积分方程,而对边界节点直接采用移动最小二乘近似函数引入边界条件来解决奇异积分问题,这同时也解决了对积分边界进行插值引入近似误差的问题．作为应用和数值实验,对Laplace方程和Helmholtz方程问题进行了分析,取得了很好的数值结果．进而,在Helmholtz方程求解中,采用了含波解信息的修正基函数来代替单项式基函数进行近似．数值结果显示,这样处理是简单高效的,在高波数声传播问题的求解中非常具有前景．     

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

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

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

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

三维瞬态弹性动力场的基本解及边界积分方程

三维瞬态动力场的求解,一直为世人所瞩目,本文拟采用边界积分方程——边界单元法进行求解。与其它数值方法相比较,本方法由于采用了动力场区域的基本解,因此在处理无限域或半无限域的动力场问题时可以毋须人为地去划定边界,附加相应的边界约束条件。同时在求解过程中,采用了权函数与奇异格林函数,使得域积分可以转化为边界积     

位势平面问题的新的规则化边界积分方程

广泛实践集中在直接变量边界积分方程的规则化研究,其本质是利用简单解消除边界积分的奇异性．然而,至今关于平面位势问题的第一类边界积分方程的规则化研究尚未涉足．致力于间接变量边界积分方程的规则化方法研究,基于一种新的思想和观点,确立平面位势问题的间接变量规则边界积分方程,它不包含CPV强奇异积分和HFP超奇异积分．数值算例表明现在的方法可取得很好的精度和效率,特别是边界量的计算．     

阻尼边界条件散射问题的数值解法

该文研究了光滑区域上二维Helmholtz方程阻尼边界条件外问题的数值解法, 应用单双层位势组合来逼近散射场, 因此积分方程中含有超奇异算子. 给出了超奇异算子的离散化方法, 在Holder空间中给出了误差估计和解析边界的收敛性分析. 最后针对该方法给出数值实例, 以表明该方法的有效性.     

多裂纹问题计算分析的本征COD边界积分方程方法

针对多裂纹问题，若采用常规的数值求解技术，计算效率较低.为实现多裂纹问题的大规模数值模拟，建立了本征裂纹张开位移（crack opening displacement, COD）边界积分方程及其迭代算法，并引入Eshelby矩阵的定义，将多裂纹分为近场裂纹和远场裂纹来处理裂纹间的相互影响.以采用常单元作为离散单元的快速多极边界元法为参照，对提出的计算模型和迭代算法进行了数值验证.结果表明，本征COD边界积分方程方法在处理多裂纹问题时取得较大的改进，其计算效率显著高于传统的边界元法和快速多极边界元法.     

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.     

3D modeling of crack growth in electro-magneto-thermo-elastic coupled viscoplastic multiphase composites

This work presents a time-domain hypersingular integral equation (TD-HIE) method for modeling 3D crack growth in electro-magneto-thermo-elastic coupled viscoplastic multiphase composites (EMTE-CVP-MCs) under extended incremental loads rate through intricate theoretical analysis and numerical simulations. Using Green’s functions, the extended general incremental displacement rate solutions are obtained by time-domain boundary element method. Three-dimensional arbitrary crack growth problem in EMTE-CVP-MCs is reduced to solving a set of TD-HIEs coupled with boundary integral equations, in which the unknown functions are the extended incremental displacement discontinuities gradient. Then, the behavior of the extended incremental displacement discontinuities gradient around the crack front terminating at the interface is analyzed by the time-domain main-part analysis method of TD-HIE. Also, analytical solutions of the extended singular incremental stresses gradient and extended incremental integral near the crack fronts in EMTE-CVP-MCs are provided. In addition, a numerical method of the TD-HIE for a 3D crack subjected to extended incremental loads rate is put forward with the extended incremental displacement discontinuities gradient approximated by the product of time-domain basic density functions and polynomials. Finally, examples are presented to demonstrate the application of the proposed method.     

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.     

Hermite三次样条多小波自然边界元法

针对应用自然边界元方法解上半平面的Laplace方程的Neumann边值问题时存在奇异积分的困难,本文提出了Hermite三次样条多小波自然边界元法.Hermite三次样条多小波具有较短的紧支集、很好的稳定性和显式表达式,而且它们在不同层上的导数还是相互正交的.因此,本文将它与自然边界元法相结合,利用小波伽辽金法离散自然边界积分方程,使自然边界积分方程中的强奇异积分化为弱奇异积分,从而降低了问题的复杂性.文中给出的算例表明了该方法的可行性和有效性.     

Temperature Stresses in a Piecewise-Homogeneous Cylindrical Shell Caused by the Presence of a Longitudinal Crack

By using the method of generalized conjugation problems, we propose a numerical scheme for investigation of the redistribution of temperature stresses in a piecewise-homogeneous cylindrical shell caused by a longitudinal crack. This scheme is based on systems of integral equations (some of them are singular) to determine the unknown jumps of integral characteristics of the disturbed temperature field and displacements and their derivatives on the line of a crack and on the interface as well as the derivatives of these integral characteristics with respect to the longitudinal coordinate at the interface.     

Some integral relations of Hankel transform type and applications to elasticity theory

The dominant part of an integral equation arising in connection with boundary value problems for the circular disc is evaluated in terms of orthogonal polynomials. This relation leads to an efficient method for numerical solution of the complete integral equation even in the presence of a complicated bounded kernel. The static problem of a circular crack in an infinite elastic body under general loads is used to illustrate vector boundary conditions leading to two coupled integral equations, while the problem of a vibrating flexible circular plate in frictionless contact with an elastic half space is solved by use of the associated numerical method.     

A discrete collocation method for Fredholm integro-differential equations with weakly singular kernels

A fully discrete version of a piecewise polynomial collocation method is constructed to solve initial or boundary value problems of linear Fredholm integro-differential equations with weakly singular kernels. Using an integral equation reformulation and special graded grids, optimal global convergence estimates are derived. For special values of parameters an improvement of the convergence rate of elaborated numerical schemes is established. Some of our theoretical results are illustrated by numerical experiments.     

A new boundary integral equation formulation for plane orthotropic elastic media

A new boundary integral equation formulation for solving plane elasticity problems involving orthotropic media is presented in this paper. Based on the real variable fundamental solutions of the considered problems, a limit theorem for the transformation from domain integral equations into boundary integral equations (BIEs) and a novel decomposition technique to the fundamental solutions, the regularized BIEs with indirect unknowns, which do not involve the direct calculation of CPV and HFP integrals, are established. The limiting process is done in global coordinates and no separate numerical treatment for strong and weak singular integrals was necessary. The current method does not need to transform the considered problems into isotropic ones as is normally done in the existing literature, so no inverse transform is required. The numerical implementation is carried out using both discontinuous quadratic elements and exact elements, which is developed to model its boundary with negligible error. The validity of the proposed scheme is demonstrated by three numerical examples. Excellent agreement between the numerical results and exact solutions was obtained even with using small amounts of element.