Solution of the planar Newtonian stick–slip problem with the singular function boundary integral method

A singular function boundary integral method (SFBIM) is proposed for solving biharmonic problems with boundary singularities. The method is applied to the Newtonian stick–slip flow problem. The streamfunction is approximated by the leading terms of the local asymptotic solution expansion which are also used to weight the governing biharmonic equation in the Galerkin sense. By means of the divergence theorem the discretized equations are reduced to boundary integrals. The Dirichlet boundary conditions are weakly enforced by means of Lagrange multipliers, the values of which are calculated together with the singular coefficients. The method converges very fast with the number of singular functions and the number of Lagrange multipliers, and accurate estimates of the leading singular coefficients are obtained. Comparisons with the analytical solution and results obtained with other numerical methods are also made. Copyright © 2005 John Wiley & Sons, Ltd.     

An analytical method for a mixed boundary value problem of circular plates under arbitrary lateral loads

In this paper by using the concept of mixed boundary functions,an analytical method is proposed for a mixed boundary value problem of circular plates.The trial functions are constructed by using the series of particular solutions of the biharmonic equations in the polar coordinate system.Three examples are presented to show the stability and high convergence rate of the method.     

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…     

A new boundary integral equation for plane harmonic functions

In this paper, we propose a new boundary integral equation for plane harmonic functions. As a new approach, the equation is derived from the conservation integrals. Every variable in the integral equation has direct engineering interest. When this integral equation is applied to the Dirichlet problem, one will get an integral equation of the second kind, so that the algebraic equation system in the boundary element method has diagonal dominance. Finally, this equation is applied to elastic torsion problems of cylinders of different sections, and satisfactary numerical results are obtained.     

势问题的重构核粒子边界无单元法

将重构核粒子法和势问题的边界积分方程方法结合,提出了势问题的重构核粒子边界无单元法. 推导了势问题的重构核粒子边界无单元法的公式,研究其数值积分方案,建立了重构核粒子边界无单元法的离散化边界积分方程,并推导了重构核粒子边界无单元法的内点位势的积分公式. 重构核粒子法形成的形函数具有重构核函数的光滑性,且能再现多项式在插值点的精确值,所以该方法具有更高的精度. 最后给出了数值算例,验证了所提方法的有效性和正确性. }      

Necessary and sufficient condition for the correct formulation of boundary integral equations of harmonic functions

A necessary and sufficient condition for the correct formulation of boundary integral equations of harmonic functions is established in the present paper. A super-determined problem of harmonic functions is proposed for the first time. Then a necessary and sufficient condition for the existence of solution of the super-determined problem is proved. At the same time, it is a necessary and sufficient condition for the correct formulation of boundary integral equations with direct unknown quantities. A relation between boundary integral equations and variational principles is discovered for the first time. And a one-to-one correspondence between boundary integral equations with direct and indirect unknown quantities is indicated. The concept and route of the present paper can be applied to other boundary value problems possessing variational principles.     

Boundary integral equations in quasisteady problems of capillary fluid mechanics Part 2: Application of the stress-stream function

In this paper planar viscous flows with a free boundary are further studied using the quasisteady approximation [1]. The introduction of the bianalytical stress-stream function provides an opportunity to adopt the theory of analytical functions. The mode of construction of the Fredholm boundary integral equations is here proposed through the explicit solutions of two Hilbert problems for holomorphic functions with the application of the conformal mappings. The stabiligy of the equilibrium of the annulus liquid layer is investigated by way of example.

---

Sommario Si prosegue lo studio di flussi piani viscosi con frontiera libera applicando l'approssimazione quasistazionaria [1]. L'introduzione della funzione stress-stream bianalitica consente l'uso della teoria delle funzioni olomorfe. La costruzione delle equazioni integrali di Fredholm al contorno proposta qui si basa sulla risoluzione esplicita di due problemi di Hilbert per funzioni analitiche mediante applicazione della tecnica delle trasformazioni conformi. Come esempio si studia la stabilità dell'equilibrio di uno strato liquido anulare.

     

各向异性势问题充要的随机边界积分方程

本文针对各向异性势问题提出了一类充分必要的随机边界积分方程。数值计算结果表明在退化尺度附近，充要的随机边界积分方程较习用的随机边界积分方程有较大的优越性。     

On the uniqueness of boundary integral equation for the exterior Helmholtz problem

From the point of view of energy analysis,the cause that the uniqueness of theboundary integral equation induced from the exterior Helmholtz problem does not hold isinvestigated in this paper.It is proved that the Sommerfeld’s condition at the infinity ischanged so that it is suitable not only for the radiative wave but also for the absorptive wavewhen we use the boundary integral equation to describe the exterior Helmholtz problem.Therefore,the total energy of the system is conservative.The mathematical dealings toguarantee the uniqueness are discussed based upon this explanation     

三维位势问题的梯度边界积分方程的新解法

三维位势问题的边界元分析中,关于坐标变量的边界位势梯度的计算是一个困难的问题. 已有一些方法着手解决这个问题,然而,这些方法需要复杂的理论推导和大量的数值计算. 本文提出求解一般边界位势梯度边界积分方程的辅助边值问题法. 该方法构造了与原边界值问题具有相同解域的辅助边值问题,该辅助边值问题具有已知解,因此通过求解此辅助边值问题,可获得梯度边界积分方程对应的系统矩阵,然后将此系统矩阵应用于求解原边值问题,求解过程非常简单,只需求解一个线性系统即可获得原边值问题的解. 值得注意的是,在求解原边值问题时,不再需要重新计算系统矩阵,因此辅助边值问题法的效率并不很差. 辅助边值问题法避免了强奇异积分的计算,具有数学理论简单、程序设计容易、计算精度高等优点,为坐标变量梯度边界积分方程的求解提供了一个新的途径. 3个标准的数值算例验证了方法的有效性.      

弹性板边界元中计算边界高阶导数的渐近收敛积分方程

本文针对板弯曲边界元方法中计算边界曲率等高阶导数项时边界积分方程中出现的高阶奇异积分项,通过对未知挠曲函数作渐近展开并加以适当摄动,获得了渐近收敛的边界积分方程。采用这一方法计算板边界上的曲率分布,获得了满意的数值结果。     

EQUIVALENT BOUNDARY INTEGRAL EQUATIONS WITH INDIRECT VARIABLES FOR PLANE ELASTICITY PROBLEMS

IntroductionIthasbeenratheralonghistorythattheBoundaryElementMethod (BEM )isappliedtosolvetheplaneelasticityproblems[1～2 ].However,theEBIE ,whichisequivalenttotheoriginalboundaryvalueproblem ,hasnotbeenfullyappreciatedandsolvedinBEMcommunity .TheconventionalboundaryintegralequationswithindirectvariablesarediscussedthoroughlyanditisshownthatthepreviousresultsarenotEBIE ,i.e .,sometimes,thereexistsnosolutionormorethanonesolutiontothem .Themainkeyliesintheexactformoftheexteriorproblems.I…     

平面问题等价边界积分方程的三次边界轮廓法

基于弹性力学平面问题等的边界积分方程,给出了三次单元的边界轮廓法。根据平面问题解的复变函数表示,构造了三次形函数。给出了对于混合边值问题求解系统方程确定的边界轮廓方程配置和三次单元界轮廓法的实施。     

Bianalytic functions, biharmonic functions and elastic problems in the plane

1　ＳｔｒｅｓｓＦｕｎｃｔｉｏｎＥｑｕａｔｉｏｎｓｏｆＰｌａｎａｒＥｌａｓｔｉｃＢｏｄｙａｎｄＣｏｎｄｉｔｉｏｎｓｏｆＤｅｆｉｎｉｔｅＳｏｌｕｔｉｏｎｓ　　Ａｎａｌｙｔｉｃｆｕｎｃｔｉｏｎｓｈａｖｅｍａｎｙｉｍｐｏｒｔａｎｔａｐｐｌｉｃａｔｉｏｎｓｔｏｔｈｅｐｒｏｂｌｅｍｓｏｆｐｌａｎａｒｅｌａｓｔｉｃｍｅｃｈａｎｉｃｓａｎｄｆｌｕｉｄｍｅｃｈａｎｉｃｓ[1～ 3].Ｉｎｏｒｄｅｒｔｏｍａｋｅａｃａｒｅｆｕｌｓｔｕｄｙｏｆｔｈｅｖｅｃｔｏｒｆｉｅｌｄｓｗｉｔｈｓｏｕｒｃｅｓａｎｄｃｕｒｌｓ,ｔｈｅｂｉａ…     

横观各向同性材料三维裂纹问题的数值分析

严格从三维横观各向同性材料弹性空间问题的Green函数出发,采用Hadamard有限部积分概念,导出了三维状态下单位位移间断(位错)集度的基本解.在此基础上,将三维任意形状的片状裂纹问题归结为求解-组以未知位移间断表示的超奇异积分方程;并给出了边界元离散形式.对方程中出现的超奇异积分,采用了Had-alnard定义的有限部积分来处理.论文最后给出了若干典型片状裂纹问题的数值算例,数值结果表明了本文方法是非常有效的.     

Multiple reciprocity method with two series of sequences of high-order fundamental solution for thin plate bending

IntroductionIt’swell_knownthatthecomplicatedfundamentalsolution[1,2 ]forHelmholtzequationΔu(x) +k2 u(x) =0　　(x∈Ω:boundedopenregioninR2 )isu (x,y) =-iH(2 )0 (k x-y ) 4,thusit’snotconvenientfornumericalcomputation .IfapplyingthesimplefundamentalsolutionofLaplaceequationu 0 (x ,y) =-ln｜x-y｜(2π) ,theexpressionforthesolutionofequationintheclosedregion Ωisc(y)u(y) + ∫Γu(x) u 0 (x,y) nx -u 0 (x ,y) u(x) n dsx =-k2∫Ωu(x)u 0 (x,y)dΩx.Astherightsideappearstheregionalintegrationinclu…     

ELASTIC ANALYSIS OF ORTHOTROPIC PLANE PROBLEMS BY THE SPLINE FICTITIOUS BOUNDARY ELEMENT METHOD

ＩｎｔｒｏｄｕｃｔｉｏｎＩｎｔｈｅｃｏｎｖｅｎｔｉｏｎａｌｂｏｕｎｄａｒｙｅｌｅｍｅｎｔｍｅｔｈｏｄ (ＢＥＭ) ,ｓｉｎｇｕｌａｒｐｏｉｎｔｓａｒｅｄｉｓｔｒｉｂｕｔｅｄａｌｏｎｇｔｈｅｂｏｕｎｄａｒｙｏｆｔｈｅｄｏｍａｉｎｕｎｄｅｒｓｔｕｄｙ ,ｗｈｉｃｈｌｅａｄｓｔｏｓｉｎｇｕｌａｒｉｎｔｅｇｒａｌｅｑｕａｔｉｏｎｓｄｕｅｔｏｔｈｅｓｉｎｇｕｌａｒｉｔｙｏｆｆｕｎｄａｍｅｎｔａｌｓｏｌｕｔｉｏｎｓ .Ｔｈｅｒｅｆｏｒｅ ,ｓｉｎｇｕｌａｒｉｎｔｅｇｒａｔｉｏｎｍｕｓｔｂｅｈａｎｄｌｅｄｗｈｅｎｓｏｌ…     

关于加权残值法的充要性问题及其权函数的选择

本文论证了偏微分方程初-边值混合问题的加权残值方程的充要性条件和有关定理。关于充要性条件的讨论也为正确选择解析函数作为权函数,建立充要的加权残值方程提供了理论指导。重点讨论了边界积分方程的充要性。指出了两种满足充要条件的权函数和相应的边界积分方程。     

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

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

弹性力学中一种新的边界轮廓法

利用基本解的特性,将面力积分方程化成仅含有Ｃａｕｃｈｙ主值积分的形式,基于这种边界积分方程,提出了一种新的边界轮廓法,对于三维问题,该方法只须计算沿边界单元界线的线积分,对二维问题,则只需计算边界单元两点的热函数之差,无须进行数值积分计算,实例计算说明该方法是有效的。