基于积分方程方法的弹性波多重散射计算

运用积分方程方法计算了含多个随机分布椭圆柱型孔洞的随机非均匀介质中相干波的速度和衰减系数,分析了这种介质的频散特性。首先,建立了散射位移场满足的积分方程,推导了单个椭圆柱孔洞的散射截面计算公式。接着分析了在含多个随机分布椭圆柱型孔洞的随机非均匀介质中弹性波的多重散射,给出在统计平均意义下的相干波的波速和衰减系数计算公式。然后用Matlab进行了编程,给出了一个数值算例,并将计算结果与波函数展开法进行了比较,分析了随机空隙介质的频散特征及其孔洞椭圆偏心率和材料空隙率的影响。     

刚性斜桩顶部受任意力的位移的线载荷积分方程法的分析

刚性斜桩顶部受任意力作用的位移分析可以分解为在倾斜平面xoz及其法平面yoz内受力的位移分析.xoz(或yoz)平面内的位移分析,可以用集度为未知函数X(t)(或Y(t))和Z(t)的Mindlin水平点力(平行x轴(或y轴)),垂直点力,在xoz(或yoz)平面内沿桩轴[0,L]内分布,根据边界条件,可将问题归结为Fredholm第一种积分方程.用离散的方法可获数值解.文中给出数值计算的例子.计算的精度用功的互等定理来检查,并将直桩的结果与别人的直桩结果作比较.     

The parameter perturbation method on elastic wave equation in inhomogeneous medium

I.IntroductionTheelasticwaveininhomogeneousmediumiscomplicatedbecauseofthediffracting,scattering,andtransmutingofthewavetapes.Exceptforsomesimpleandregularmediummode1s,thesolutionofelasticwavehasnotbeengotyet.Nowadays,theresearchoftheelasticwavescattering…     

Using fredholm integral equation of the second kind to solve the vertical vibration of elastic plate on an elastic half space

I.IntroductionThedynamicanalysisofelasticplateonanelasticplateisveryimportanttotheresearchandengineering.Article[l]--131studiedtheverticalvibratiollof'rigidplatcl']alldelastic1llatcwitharigidcorely]orarigidedgel'joilallelastichalfspace.BecauseloadillgactatrigidPltltcof'therigidregionofelasticplate,onlytheresultantofh7rcesiscollsidcrcdillldwitll11oI.cliltlolltothedistributionofloading.Thevibrationofplateiscollsidcredastileproblelllof-sillglcdegreeoffi.eedolnsystem..Thepurposeofthepresentpaper…     

基体裂纹与界面刚性线的弹性干涉

研究两种材料界面上的刚性线与其它任意位置处直线裂纹弹性干涉的反平面问题。基于界面上刚性线与任意位置处螺型位错干涉的基本解,运用连续位错密度模型法将问题转化为奇异积分方程。用半开型积分法求解奇异积分方程,得到位错密度函数的离散值,计算裂纹尖端处的应力强度因子。算例说明该方法可用于工程实际问题。     

一类非线性奇异积分方程及其数值方法研究

探讨了一类非线性奇异积分方程的数理性质以及在颗粒雷诺数Rep＜1时此类方程解的存在条件，然后详细研究了该方程的数值计算方法并构造称之为P（EC）^k多步法的差分格式，分析了该格式的收敛性和代数精度，得到时域离散步长的约束关系。运用该格式计算了静止流场和均匀振荡流中球形小颗粒的非恒定运动，将计算结果与其解析解及有关实验数据的比较表明，上述数值方法具有良好的计算精度。     

Improved non-singular local boundary integral equation method

When the source nodes are on the global boundary in the implementation of local boundary integral equation method (LBIEM),singularities in the local boundary integrals need to be treated specially. In the current paper,local integral equations are adopted for the nodes inside the domain trod moving least square approximation (MLSA) for the nodes on the global boundary,thus singularities will not occur in the new al- gorithm.At the same time,approximation errors of boundary integrals are reduced significantly.As applications and numerical tests,Laplace equation and Helmholtz equa- tion problems are considered and excellent numerical results are obtained.Furthermore, when solving the Hehnholtz problems,the modified basis functions with wave solutions are adapted to replace the usually-used monomial basis functions.Numerical results show that this treatment is simple and effective and its application is promising in solutions for the wave propagation problem with high wave number.     

Isotropicalized spline integral equation method for the analysis of anisotropic plates

In the view of Reissner’s and Kirchhoff’s theories,respectively,we formulate theisotropicalized governing equations for the anisotropic plates,and give the proof of theequivalence relation between these two plate-models for the simply-supported rectangularorthotropic plates.The well-known fundamental solutions of the isotropic plates are utlizedfor the spline integral equation analysis of anisotropic plates.Even with sparse meshes thesatisfactory results can be obtained.The analysis of plates on two-parameter elasticfoundation is so simple as the common case that only a few terms should be added to theformulas of fictitious loads.     

Dynamic interaction of a pile with a transversely isotropic elastic half-space under transverse excitations

This investigation is concerned with a mathematical analysis of an elastic circular cylindrical pile embedded in a transversely isotropic half-space under lateral dynamic excitations. A combination of time-harmonic horizontal shear force and moment are applied at the top end of the pile. The boundary value problem is formulated by decomposing the pile-medium system into a fictitious pile and an extended transversely isotropic half-space. A Fredholm integral equation of the second kind governs the interaction problem, whose solution is then computed numerically. Selected results for dynamic compliance bending moment, displacement and slope profiles are presented for different transversely isotropic half-spaces to portray the influence of degree of anisotropy of the medium on various aspects of the solution.     

Solution of A 2-D weak singular integral equation with constraint

In this paper, the solution of a 2-D weak singular integral equation of the first kind

SOLVING VIBRATION PROBLEM OF THIN PLATESUSING INTEGRAL EQUATION METHOD

ＳＯＬＶＩＮＧＶＩＢＲＡＴＩＯＮＰＲＯＢＬＥＭＯＦＴＨＩＮＰＬＡＴＥＳＵＳＩＮＧＩＮＴＥＧＲＡＬＥＱＵＡＴＩＯＮＭＥＴＨＯＤ￥(许明田，程德林)ＸｕＭｉｎｇｔｉａｎ；ＣｈｅｎｇＤｅｌｉｎ（Ｄｅｐａｒｔｍｅｎｔ．ｏｆＭａｔｈｅｍａｔｉｃｓａｎｄＰｈｙｓｉ...     

Perturbation finite volume method for convective-diffusion integral equation

A perturbation finite volume (PFV) method for the convective-diffusion integral equation is developed in this paper. The PFV scheme is an upwind and mixed scheme using any higher-order interpolation and second-order integration approximations, with the least nodes similar to the standard three-point schemes, that is, the number of the nodes needed is equal to unity plus the face-number of the control volume. For instance, in the two-dimensional (2-D) case, only four nodes for the triangle grids and five nodes for the Cartesian grids are utilized, respectively. The PFV scheme is applied on a number of 1-D linear and nonlinear problems, 2-D and 3-D flow model equations. Comparing with other standard three-point schemes, the PFV scheme has much smaller numerical diffusion than the first-order upwind scheme (UDS). Its numerical accuracies are also higher than the second-order central scheme (CDS), the power-law scheme (PLS) and QUICK scheme. The project supported by the National Natural Science Foundation of China (10272106, 10372106)     

Solution of nonlinear wave equation of elastic rod

The longitudinal oscillation of a nonlinear elastic rod with lateral inertia are studied. A nonlinear wave equation is derived. The equation is solved by the method of full approximation.     

A method using the reciprocal theorem to solve the bending of thin elastic semicircular plates

The bending of the thin elastic semicircular plates, because of its complicated boundary conditions, brings some difficulties for us to obtain its solution. This paper applies the reciprocal theorem to propose a general simple convenient method to obtain the transverse defiectional equations of the plates.     

BOUNDARY INTEGRAL FORMULA OF ELASTIC PROBLEMS IN CIRCLE PLANE

IntroductionConcerningtheelasticplaneprobleminaunitcircle ,ZhengShenzhouandZhengXueliangdevelopedaboundaryintegralformulaofthestressfunction[1]:Φ(r,θ) =-( 1 -r2 ) 24π ∫2π0ν( φ)1 -2rcos(θ-φ) r2 dφ　　 12π∫2π011 -2rcos(θ-ω) r2 dω∫2π0μ( φ)1 -cos(ω-φ) dφ　　 1 -r22π∫2π0μ( φ)1 -2rcos(θ -φ) r2 dφ　　 ( 0 ≤r <1 ) ,( 1 )whereμ(θ) =Φ(r,θ) ｜r=1,ν(θ) = Φ n r=1= Φ r r=1.Intheformula ( 1 )theseconditemisastrongsingularintegral,itshouldbeunderstoodasanintegra…     

应用边界积分法求圆形夹杂问题的解析解

边界元方法作为一种数值方法,在各种科学工程问题中得到了广泛的应用.本文参考了边界元法的求解思路,从Somigliana等式出发,利用格林函数性质,得到了一种边界积分法,使之可以用来寻求弹性问题的解析解.此边界积分法也可以从Betti互易定理得到.应用此新方法,求解了圆形夹杂问题.首先设定夹杂与基体之间完美连接,将界面处的位移与应力按照傅里叶级数展开,根据问题的对称性与三角函数的正交性来简化假设,减少待定系数的个数.其次选择合适的试函数(试函数满足位移单值条件以及无体力的线弹性力学问题的控制方程),应用边界积分法,求得界面处的位移与应力的值.然后再求解域内位移与应力.得到了问题的精确解析解,当夹杂弹性模量为零或趋向于无穷大时,退化为圆孔或刚性夹杂问题的解析解.求解过程表明,若问题的求解区域包含无穷远处时,所取的试函数应满足无穷远处的边界条件.若求解区域包含坐标原点,试函数在原点处位移与应力应是有限的.结果表明了此方法的有效性.     

Analysis of a three-dimensional dissimilar material joint with one real singularity using a conservative integral

In the present study, a conservative integral based on the Betti reciprocal principle is formulated to determine the intensity of singularity at a vertex of the interface in three-dimensional dissimilar material joints with one real singularity. Eigenanalysis formulated using a three-dimensional finite element method (FEM) is used to calculate the order of stress singularity, angular functions of displacements and stresses. Models with various element sizes and various integral areas are used to investigate the effect of the integration area on the accuracy of the results. The results are compared with those obtained from the boundary element method (BEM) using a curve-fitting technique to calculate the intensity of singularity. In addition, models of various lengths and various material combinations are used to investigate the stress singularity characteristics in three-dimensional dissimilar material joints. The results of the present study indicate that the conservative integral can be used to determine the intensity of singularity in three-dimensional bi-material joints. The accuracy of the results can be improved by mesh refinement. Finally, the relationships among the intensity of singularity, the order of stress singularity and the model geometry are discussed.     

On some basic principles in dynamic theory of elastic materials with voids

According to the basic idea of dual-complementarity, in a simple and unified way proposed by the author^[1], some basic principles in dynamic theory of elastic materials with voids can be established systematically. In this paper, an important integral relation in terms of convolutions is given, which can be considered as the generalized principle of virtual work in mechanics. Based on this relation, it is possible not only to obtain the principle of virtual work and the reciprocal theorem in dynamic theory of elastic materials with voids, but also to derive systematically the complementary functionals for the eight-field, six-field, four-field and two-field simplified Gurtin-type variational principles. Furthermore, with this approach, the intrinsic relationship among various principles can be explained clearly. The project supported by the Foundation of Zhongshan University Advanced Research Center     

A boundary integral equation method for Navier-Stokes equations—application to flow in annulus of eccentric cylinders

A novel Navier-Stokes solver based on the boundary integral equation method is presented. The solver can be used to obtain flow solutions in arbitrary 2D geometries with modest computational effort. The vorticity transport equation is modelled as a modified Helmholtz equation with the wave number dependent on the flow Reynolds number. The non-linear inertial terms partly manifest themselves as volume vorticity sources which are computed iteratively by tracking flow trajectories. The integral equation representations of the Helmholtz equation for vorticity and Poisson equation for streamfunction are solved directly for the unknown vorticity boundary conditions. Rapid computation of the flow and vorticity field in the volume at each iteration level is achieved by precomputing the influence coefficient matrices. The pressure field can be extracted from the converged streamfunction and vorticity fields. The solver is validated by considering flow in a converging channel (Hamel flow). The solver is then applied to flow in the annulus of eccentric cylinders. Results are presented for various Reynolds numbers and compared with the literature.     

On the regularization method of the first kind of Fredholm integral equation with a complex kernel and its application

The regularized integrodifferential equation for the first kind of Fredholm integral equation with a complex kernel is derived by generalizing the Tikhonov regularization method and the convergence of approximate regularized solutions is discussed. As an application of the method, an inverse problem in the two-demensional wave-making problem of a flat plate is solved numerically, and a practical approach of choosing optimal regularization parameter is given. Project supported by the National Natural Science Foundation, of China