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

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

Transient elastodynamic three-dimensional problems in cracked bodies

The general formulation of the transient elastodynamic second boundary value problem in an isotropic linear elastic body with a crack of arbitrary shape by combining the boundary integral equation method and the Laplace transform with respect to time is presented in this paper. Both finite and infinite elastic bodies are considered. A numerical solution of the transformed boundary integral equations is proposed.     

Linear integral equations for conformal mapping of bounded multiply connected regions onto a disk with circular slits

In this paper we present a boundary integral equation method for the numerical conformal mapping of bounded multiply connected region Ω onto a disk with circular slits. The method is based on some uniquely solvable boundary integral equations with classical adjoint and generalized Neumann kernels. These boundary integral equations are constructed from a boundary relationship satisfied by a function analytic on a multiply connected region. Some numerical examples are presented to illustrate the efficiency of the presented method.     

A wavelets method for plane elasticity problem

A Neumann boundary value problem of plane elasticity problem in the exterior circular domain is reduced into an equivalent natural boundary integral equation and a Poisson integral formula with the DtN method. Using the trigonometric wavelets and Galerkin method, we obtain a fast numerical method for the natural boundary integral equation which has an unique solution in the quotient space. We decompose the stiffness matrix in our numerical method into four circulant and symmetrical or antisymmetrical submatrices, and hence the solution of the associated linear algebraic system can be solved with the fast Fourier transform (FFT) and the inverse fast Fourier transform (IFFT) instead of the inverse matrix. Examples are given for demonstrating our method has good accuracy of our method even though the exact solution is almost singular.     

含圆形孔口和水平直裂纹的平面问题的应力分析

利用复变函数方法和积分方程理论研究了既含有圆形孔口又含有水平裂纹的无限大平面的平面弹性问题,将复杂的解析函数的边值问题化成了求解只在裂纹上的奇异积分方程的问题.此外,还给出了裂纹尖端附近的应力场和应力强度因子的公式.     

Dynamic stresses in a compound body with circular crack under sliding contact on an interface

We have considered the three-dimensional problem of harmonic loading of a circular crack in an elastic composite consisting of two dissimilar half-spaces under sliding contact on the surface of their bonding. The defect is situated in one of the half-spaces perpendicularly to the interface of materials. Using the representations of solutions in the form of Helmholtz potentials, we have reduced the problem to a boundary integral equation for the function of dynamic defect opening. Based on the numerical solution of this equation, we have obtained the frequency dependences of mode I stress intensity factor near the crack for different relations between the elastic moduli of components of the composite and the depths of crack location with respect to the interface.     

新的三维力学GELD正演和反演算法

在本文中 ,我们提出了新的整体积分和局部微分GILD的力学正演和反演方法 .我们建立了弹性和塑性力学的体积分微分方程 .我们证明了这个体积分方程和伽辽金虚功原理等价 .新的GILD方法是基于这个体积分微分方程 .GL方法是进一步的发展 ,GL是一种整体场和局部场相互作用的全新方法 .在这个方法中 ,仅仅需要解 3× 3或者 6 × 6的局部小矩阵 .特别是 ,用GL方法求解无限域的偏微分方程时 ,不需要任何人工边界 ,不需要任何吸收边界条件和不需要任何边界积分方程 .新的三维力学GILD正演和反演算法已被应用研究奈米材料的力学性质的模拟计算 .我们获得非常好的奈米材料的力学变形的超拉力的力学性质 .我们提出了新的奈米地球物理新概念和发现了GILD数值量子     

层状横观各向同性饱和地基上圆板的非轴对称振动

研究了层状横观各向同性饱和地基上弹性圆板的非轴对称振动问题．首先,通过方位角的Fourier变换,将圆柱坐标系下横观各向同性饱和土的三维动力方程转化为一阶常微分方程组,基于径向Hankel变换,建立问题的状态方程,求解状态方程后得到传递矩阵;其次,利用传递矩阵,结合层状饱和地基的边界条件、排水条件及层间接触和连续条件,给出了任意简谐激振力作用下层状横观各向同性饱和地基动力响应的通解;然后,按混合边值问题建立层状饱和地基上弹性圆板非轴对称振动的对偶积分方程,并将对偶积分方程化为易于数值计算的第二类Fredholm积分方程,并给出了算例．     

Three-dimensional singularities of elastic fields near vertices

For the computation of the singular behavior of an elastic field near a three-dimensional vertex subject to displacement boundary conditions we use a boundary integral equation of the first kind whose unknown is the boundary stress. Localization at the vertex and Mellin transformation yield a one-dimensional integral equation on a piecewise circular curve γ in IR³ depending holomorphically on the complex Mellin parameter. The corresponding spectral points and packets of generalized eigenvectors characterize the desired stress field and are computed by a spline-Galerkin method with graded meshes at the corner points of the curve γ. © 1993 John Wiley & Sons, Inc.     

抛物型初边值问题的自然积分方程及其数值解法

１．引言数值求解无界区域的偏微分方程,自然的处理方式是削去区域的无界部分,即引入一条适当的人工边界ｒ。,将原问题的求解限制在一个适当的有界区域Ｄ内,这样必须在人工边界上引入适当边界的条件．于是很自然地导致这样一个问题：'是否存在一个人工边界条件,使得在这边界条件下,原问题在区域Ｄ内所求得的数值解与原无界区域的解在Ｄ上的限制是完全一致的？"这里我们的着眼点是寻求与原无界区域问题等价的数学形式,以便于数值求解．因为边界元方法可以将区域内的问题转化到区域的边界上去处理,经典的边界元方法常被应用．七十年代…     

Stress field around an arbitrary thin inclusion in a transversely isotropic elastic half-space

We consider a thin flat inclusion of arbitrary shape located inside a transversely isotropic elastic half-space in the plane parallel to its boundary z = 0. An arbitrary tangential displacement is prescribed on the inclusion. The boundary of the half-space is stress-free. We need to find the complete field of stresses and displacements in this half-space. A governing integral equation is derived by the generalized method of images, introduced by the author. The case of circular inclusion is considered as an example. Two methods of solution of the governing integral equation are derived. A detailed solution is presented for the particular cases of radial expansion, torsion and lateral displacement of the inclusion. The solution is also valid for the case of isotropy. The governing integral equation for the case of isotropy is derived.     

无穷凹角区域各向异性问题的自然边界元与有限元耦合法

本文研究无穷凹角区域上一类各向异性问题的自然边界元与有限元耦合法.利用自然边界归化原理,获得圆弧或椭圆弧人工边界上的自然积分方程,给出了耦合的变分形式及其数值方法,以及逼近解的收敛性和误差估计,最后给出了数值例子,以示方法的可行性和有效性.     

Iterative solution of the generalized Love's problem in anisotropic media

An integral equation method is used to show the well posedness of the generalized Love's problem for an elastic anisotropic layer superimposed to a homogeneous substrate. The solution is derived using an iterative technique and the dispersion equation is obtained by imposing the pertinent boundary conditions. Inhomogeneous layers with different profiles of the material parameters are considered as examples and numerical results are given. The corresponding generalized problem is discussed for Lamb's modes.     

Dirichlet-Neumann alternating algorithm for an exterior anisotropic quasilinear elliptic problem

In this paper, by the Kirchhoff transformation, a Dirichlet-Neumann (D-N) alternating algorithm which is a non-overlapping domain decomposition method based on natural boundary reduction is discussed for solving exterior anisotropic quasilinear problems with circular artificial boundary. By the principle of the natural boundary reduction, we obtain natural integral equation for the anisotropic quasilinear problems on circular artificial boundaries and construct the algorithm and analyze its convergence. Moreover, the convergence rate is obtained in detail for a typical domain. Finally, some numerical examples are presented to illustrate the feasibility of the method.     

A fast numerical method for a natural boundary integral equation for the Helmholtz equation

A Neumann boundary value problem of the Helmholtz equation in the exterior circular domain is reduced into an equivalent natural boundary integral equation. Using our trigonometric wavelets and the Galerkin method, the obtained stiffness matrix is symmetrical and circulant, which lead us to a fast numerical method based on fast Fourier transform. Furthermore, we do not need to compute the entries of the stiffness matrix. Especially, our method is also efficient when the wave number k in the Helmholtz equation is very large.     

无界区域抛物方程自然边界元方法

本文应用自然边界元方法求解无界区域抛物型初边值问题。首先将控制方程对时间进行离散化,得到关于时间步长离散化的椭圆型问题。通过Fourier展开,导出相应问题的自然积分方程和Poisson积分公式。研究了自然积分算子的性质,并讨论了自然积分方程的数值解法,最后给出数值例子。从而解决了抛物型问题的自然边界归化和自然边界元方法。     

Pressure of a plane stamp of nearly circular cross section on an elastic half-space : PMM vol. 40, n 6, 1976, pp. 1143–1145

An approximate method of solving the contact problem of impressing a plane stamp of nearly circular cross section into an elastic half-space is suggested. The friction of the contact surface is neglected. A numerical algorithm for the method is produced. An elliptical and rectangular stamps are considered as examples.There is no general method of solving the problems for stamps of nearly circular cross section. Apart from the classical problem of a plane elliptical stamp, the literature gives solutions for the problems of polygonal stamps, with each problem however requiring a different approach. An approximate solution for the problem of impressing a stamp of nearly circular cross section into an elastic half-space is given in [1]. The method makes it possible to use the same approach to solve the contact problem for an arbitrary region of contact, and to construct an universal numerical algorithm. The program can be adapted to each particular case by making the corresponding changes in the procedure of computing the Fourier coefficients of the equation of the boundary of the area of contact. Below a numerical algorithm for the approximate method in question is given. A more effective formulation of the solution is given for the case of the elliptical stamp.     

含孔薄板弯曲波动的双互易边界元法

采用双互易边界元法对开孔无限大薄板弹性波的散射与动应力集中问题进行理论分析和数值计算．基于功的互等定理,采用静力基本解建立了薄板弯曲波动问题的双互易边界积分方程．作为数值算例,计算了圆孔附近的动应力集中系数,通过与已有结果进行比较,表明该方法简单有效并能够保证计算精度．     

Free vibration analysis of discrete-continuous functionally graded circular plate via the Neumann series method

The Neumann series method has been used for the first time to solve the boundary value problem of free axisymmetric and nonaxisymmetric vibrations of continuous and discrete-continuous functionally graded circular plate on the basis of the classical plate theory. The equation of motion and the general solution for a functionally graded circular plate with a very complex system of a discrete elements attached, such as concentric ring masses, elastic supports, rotational springs, and damping elements are presented for the first time. The particular continuous solutions to the defined differential equations are obtained as the Neumann power series rapidly, absolutely, and uniformly convergent to the exact eigenfrequencies for any physically justified values of the plate's parameters on the basis of the properties of the obtained closed-form kernels of the Volterra integral equations. The multiparametric nonlinear characteristic equations for plate with classical and nonclassical boundary conditions are defined and numerically solved to obtain the full spectrum of eigenfrequencies in a simple way. The effects of the position and stiffness of ring supports and of singularities as the radii of supports shrink to the center of the plate on the dimensionless eigenfrequencies of homogeneous and functionally graded circular plate with sliding support and elastic constraints are comprehensively studied and presented for the first time. The accuracy of the proposed low-computational-cost method is demonstrated by comparison of the numerical results with those available in the literature.     

Contact with break-off under bending of an elastic strip with a rigid disk

We consider the problem of the theory of elasticity of the contact interaction of a rigid circular disk and an elastic strip, which rests upon two supports with disturbance of contact in the middle part of the contact region. On the basis of the Wiener–Hopf method, an integral equation of the problem is reduced to an infinite system of algebraic equations. The size of the zone of break-off of the boundary of the strip from the disk and the distribution of contact stresses are determined.