基于双层位势的虚边界元计算方法
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
用边界元法来求解位势问题有效而简单,但通常需要求解奇异积分,特别是当公式中有双层位势的法向导数时,会遇到超强奇异积分。若采用虚边界元法[23]就可以避开这些弱点,通过在所研究的区域之外的延拓区域中选择一个虚拟的边界,根据原问题的边界条件确定出在虚边界上分布的虚拟密度,以此来求解原区域上的解。在数值计算时,是对虚边界进行离散,由于场点和源点在不同的边界上,从而避免了传统边界元中关于奇异与超强奇异积分的处理。
     通常所谓的虚边界元方法都是基于单层位势的延拓,利用在虚拟边界上的密度函数作为间接未知量,令此虚密度函数对真实边界所产生的位势及通量与边界条件一致,从而达到求解虚拟密度函数的目的[58]。
     本文则利用在虚边界上分布的矩密度,得出基于双层位势的一种虚边界积分方程。作为一种间接边界积分公式,由于积分点和场点分别位于虚边界和实边界上,可以避免强奇异积分和超强奇异积分的计算,并将其应用到各种椭圆边值问题的求解上。
     首先,本文针对Laplace方程的三类内边值问题建立出相应的基于双层位势的虚边界积分方程,并采用常单元来求解。其次,对广义线性Poisson类方程,本文考虑用将虚边界元和径向基函数耦合的方式来求解,对齐次方程的解采用虚边界元求解,对特解则用径向基函数逼近,从而避免区域积分项处理的繁琐。再者,进一步考虑,对双调和方程,通过引入中间变量,将其分解为两个Poisson方程,再用径向基函数同虚边界元耦合的方法求解,并应用于非线性双调和方程的求解。最后,利用自己编写的程序,通过数值算例验证了该方法的可行性和正确性。
The boundary element method is an effective numerical method for solving potential problem, but it needs to calculate singular integrals. Especially, we will meet the hyper-singular integral when we calculate the normal derivative of a double layer potential. When virtual boundary method is used,the singularity of integral can be avoided. We set a closed virtual curve, so called virtual boundary outside the domain under consideration. So the unknown virtual density function distributed on the virtual boundary can be determined by the boundary conditions given on the physical boundary via nonsingular integral equation. In calculation, since the source points are on the physical boundary, and the integral points on virtual boundary respectively, so that the integral is no longer singular or hyper-singular. The formulations used in virtual boundary are the indirect boundary integrals.
     The virtual boundary method used before is based on single layer potential, in which the indirect unknown function is the density function distributed on the virtual boundary. And then the potential and its derivative on the physical boundary can be determined by virtual density function, which can be obtained according to the boundary conditions of the original problem.
     In this paper, we use another kind of virtual boundary elements method, which is based on double layer potential by using of the moment density distributed on the virtual boundary. Since there is a distance between the virtual boundary and physical one, the integral is no longer singular. It can be applied into various elliptic boundary condition problems.
     Firstly, we write the virtual boundary integral equation according to three boundary conditions for Laplace problem and solve it by constant elements. Secondly, for generalized Poisson equation, we combine virtual boundary method with RBFs to solve it. We use RBF approximation is employed to construct particular solution and VBCM for homogenous solution instead of domain discretization. Thirdly, for harmonic equation,we divide it into two Poisson equations by introducing intermediate function. And then, we can solve it by use of virtual boundary method and RBF, this approach can be used to solve nonlinear harmonic equation also. Finally, we write a program in Matlab and did some numerical tests, the numerical results of examples demonstrate that the scheme presented is effective and accurate.
引文
[1]张有天,王镭,陈平著.边界元方法及其在工程中的应用[M].北京:水利电力出版社,1989.
    [2]余德浩,汤华中著.微分方程数值解法[M].北京:科学出版社,2003.
    [3]杨德全,赵忠生著.边界元理论及应用[M].北京:北京理工大学出版社,2002.
    [4]姜礼尚,庞之垣著.有限元方法及其理论基础[M].北京:人民教育出版社,1979.
    [5]张鸿庆,王鸣著.有限元的数学理论[M].北京:科学出版社,1991,12.
    [6]祝家麟.边界元方法—一个老而新的发展中的数值方法[J],数学的实践与认识,1983,(3):74-79.
    [7]祝家麟.边界元方法――它的产生、内容和意义,第一届工程中的边界元法会议论文集[C],杜庆华,1985.
    [8]姚寿广著,边界元数值方法及其工程应用[M].北京:国防工业出版社,1991.
    [9]申光宪等著,边界元法[M].北京:机械工业出版社,1998.
    [10]稽醒,臧跃龙,程玉民著,边界元法进展及通用程序[M].上海:同济大学出版社,1997.
    [11]孙炳楠等著,工程中边界元法及其应用[M].浙江:浙江大学出版社,1991.
    [12] Jaswon,M.A.,Symm,G.T.,Integral Equation Methods in potential theory and elastostatics[M]. London :Academic Press,1977.
    [13]余德浩著.自然边界元方法的数学理论[M].北京:科学出版社,1993.
    [14]祝家麟著.椭圆边值问题的边界元分析.[M].北京:科学出版社,1991.
    [15]严重著.边界单元法基础[M].重庆:重庆大学出版社,1986.
    [16]王元淳著.边界元法基础.[M].上海:上海交通大学出版社,1988.
    [17]杜庆华著.边界积分方程方法-边界元法[M].北京:高等教育出版社,1989.
    [18]卢盛松主编.边界元理论及应用[M].北京:高等教育出版社,1990.
    [19]张凯院,高维数值积分的一种计算方法[J].工程数学学报(14)3:91-95.
    [20] BREBBIA C.A. The Boundary Element Method for Engineers.[M].London: Pentech Press. 1978.
    [21] C.A.Brebbia.J.C.F.Telles.L.C.Wrobel,Boudary Element techniques–Theory and Applications in Engineering[M].Springer-Verlag,Berlin Heidelberg New York Tokyo,1984 .
    [22]冯康著.数值计算方法[M].北京:国防工业出版社,1978.
    [23]翟景春,Laplace方程边值问题的一种数值解法[J],山东科学,1995(3):21-29.
    [24]陈浩,二维Laplace方程解的性质[J],工程数学,1995(12)2:95-99.
    [25]祝家麟,边界元方法中的奇异性[J].重庆建筑工程学院学报1985(13)2:90—102.
    [26]谷超豪等,数学物理方程[M].北京:人民教育出版社,1979.
    [27]胡建伟,汤怀民著.微分数值方法[M].北京:科学出版社.1999.1.
    [28]李岳生,黄友谦著.数值逼近[M].北京:人民教育出版社,1978.
    [29]李荣华.解边值问题的Galerkin方法[M].上海:上海科学技术出版社.1988.12.
    [30]韩厚德,关治,俞崇庆.一个电镀模型的正则边界元分析[J],高校应用数学学报. 1988.3(1):101-111.
    [31]张凯,祝家麟,张洁.椭圆单侧问题的边界元计算方法[J],重庆大学学报,2003.26(10):27-30.
    [32]张太平.祝家麟.抛物型方程的边界元重叠性区域分解法[J].重庆大学学报. 2002(10):27-31.
    [33]祝家麟.用边界积分方程法解平面双调和方程的Dirichlet问题[J].计算数学,1984(3):278-286.
    [34]张太平,抛物型方程的边界元区域分解法[D],重庆大学硕士论文集,2001.
    [35]张守贵,祝家麟,用双层位势求解二维Laplace方程的Neumann问题的Galerkin边界元解法[D],重庆大学硕士论文集,2005.
    [36]董海云,祝家麟,二维Laplace方程的直接边界积分方程的Galerkin边界元解法[D],重庆大学硕士论文集,2005.
    [37]张瑞丰等著,精通MATLAB6.5[M].北京:中国水利水电出版社,2003.
    [38]赵翔龙著,Fortron90学习教程[M].北京:北京大学出版社,2000.
    [39]徐明编,Fortran power station 4.0基础教程[M].北京:清华大学出版社,2000.
    [40]导向科技编著,MATLAB6.0程序设计与实例应用[M],北京:中国铁道出版社,2001.
    [41] AITHISON J M, POOL M W. A numerical algorithm for the solution of Signorini problems[J],Journal of Computational and Applied mathematics1998(94):55-67.
    [42] R.Glowinski,J.L.lions,andR.Tremolieres, Numerical analysis of variational inequalities[M], North-Holland:Amsterdam,1981:82-125.
    [43] J.M. Aitchison, A.A. Lacey, M. Shillor, A model for an electropaint process[J], IMA J. Appl. Math. 1984(33): 17-31.
    [44] Bakhvalov,N.S.and Kobel’kov,G.M.,An iterative method for solving elliptic problems with a rate of convergence that does not depend on the range of coefficients[J],Otdel.Vychisl.mat. AN.SSSR, 1988(190):128-79.
    [45] Kuznetsov,YU.A.,Multilevel domain decomposition methods[J].Applied Numer.Math., 1990(6):303—314.
    [46] Glowinsk,R.,Golub,G.H.,Meurant,G.A.and Periaux,J.,eds,Proceeding of First International Symposium on Domain Decomposition Methods For Partial Differential Equation[M]. SIAM:Philadedphia,PA,1988.
    [47]第二届工程中边界元法会议文集[C],1988,南宁.
    [48]第六届工程中边界元法会议文集[C],2000,重庆.
    [49]张雄,刘岩著.无网格法[M].北京:清华大学出版社.2004.
    [50] Von Karman,T.,calculation of Pressure Distribution on Airship Hulls[J],N.A.C.A. TM1930: 574.
    [51] Cerkirge,H.M. and W.J.Lehr,Auxiliary source point method: the Laplace equation[J], Eng.Analysis,1987(4):186-189.
    [52] Burges G J,Mahajerin E A. A comparison of the boundary element and superposition methods[J].Comput.Structures,1984(19):697-712.
    [53]李晓川,电磁弹性固体辛对偶体系及虚边界元数值法[D],大连理工大学大学博士学位论文. 2007.6
    [54]许永林,唐锦春.域外奇点法及格林公式法解薄板弯曲问题[J].计算结构力学及其应用,1986.3(2):49-55.
    [55]卢西林,杜庆华无奇异边界元法薄板弯曲问题[J].清华大学学报,1989,6(3): 32-37.
    [56]王有成等,边界元技术中的全特解场法[J].力学学报,1995,27(4):33-37.
    [57]王有成.边界元技术新探--边界点方法[J].力学与实践,1997,19(4):42-46.
    [58]孙焕纯,张立洲,许强,张耀明著,无奇异边界元法.[M]大连:大连理工大学出版社.1999:1-25.
    [59]孙焕纯,杨海天.虚边界元法的应用及其求解方法[J].应用力学学报,1994,11(1):44-60.
    [60] Golberg MA, Chen CS. The method of fundamental solutions for potential,Helmholtz and diffusion problems[C].In:Golberg MA,editor.Boundary integral methods-numerical and mathematical aspects.Computational Mechanics Publication;1998.103-76.
    [61]许强,孙焕纯.虚边界元-最小二乘配点法[J].计算力学学报,1997,14(2):55-61.
    [62]邹广德,孙焕纯,柴山等,求解位势问题的虚边界元法[J],计算结构力学及其应用.1994.11(3):271-277.
    [63]陈传璋,侯宗义,李明忠.积分方程论及其应用[M].上海:上海科学技术出版社.1987.
    [64] Onate E, Idelsohn S,Zienkiewicz O C et al. A finite point method in computational mechanics:Applications to convective transport and fluid flow[J].Int. J. Num. Meth Engng., 1996(39):3839~3866.
    [65] Onate E, Idelsohn S, Zienkiewicz O C et al. A stabilized finite point method for analysis of fluid mechanics problems[J]. Comput. Methods Appl.Mech.Engrg.,1996(139):315~346.
    [66]张雄,刘小虎,宋康祖.Least-square collocation meshless method[J],Int J.Num.Meth. Engng.,2001,51(9):1089~1100.
    [67]张雄,胡炜,潘小飞等.加权最小二乘无网格法[J].力学学报,2003,35(4):425~431.
    [68]刘岩,张雄,陆明万.Meshless Weighted Least-Square (MWLS) Method for Steady-state Heat Conduction Problems[C].见:袁明武,陈璞编。中国计算力学大会’2003论文集工程与科学中的计算力学.北京,2003.
    [69] Franke C, Schaback R. Solving partial differential equations by collocation using radial basis functions[J]. Appl. Math. Comput.1998(93):73~82.
    [70] Coleman C J.On the use of radial basis functions in the solution of elliptic boundary value problems[J].Comput.Mech.,1996,(17):418~422.
    [71] Wu Z M. Solving PDE with radial basis function and the error estimation[C]. In:Chen Z,Li Y,Micchelli CA et al.(eds.). Advances in Computational Mathematics, Lecture Notes on Pure and Applied Mathematics.Guangzhou,1998.202.
    [72] Wu Z M.Schaback R. Local error estimates for radial basis function interpolation of scattered data[J].IMA Journal of Numerical Analysis,1993(13):13~27.
    [73] Chen W. Symmetric boundary knot method[J]. Engng. Anal. Bound. Elem.,2002(6):489~494.
    [74] Chen W. Meshfree boundary particle method applied to Helmholtz problems[J].Engng.Anal. Bound.Elem.,2002,26(7):577~581.
    [75] Chen W,Tanaka M. A meshless,integration-free and boundary–only RBF technique[J]. Comput Math Appl 2002(43):379-91.
    [76] Song K, Zhang X, Lu M.Collocation with Modified Compactly Supported Radial basis functions[C] .Proceedings of international conference on computational Engineering & Sciences,August 19-25,2001,Puerto Vallarta.
    [77] Nardini D,Brebbia C A. A new approach to free vibration analysis using boundary element[C] , in Boundary Element Methods in Engineering(ed. Brebbia C A). Springer-Verlag.Berlin and New York:1982.
    [78] Nardini D,Brebbia C A. Boundary integral formulation of mass matrices for dynamics analysis[C], In Topics in Boundary Element Research.Vol.2(ed.Brebbia CA).Berlin and New York:Springer-Verlag,1985.
    [79] Brebbia C A, Nardini.D Dynamic Analysis in Solid Mechanics by Alternative Boundary Element Approach[J] .Int J.Soid Dynamics and Earthquake Engineering,1983,2(4):228-233.
    [80] Katsikadelis JT. The analog equation method–a powerful BEM-based solution technique for solving linear and nonlinear engineering problems[C]. In:Brebbia CA, editor.Boundary element method XVI. Southampton:CLM Publications;1994.p.167.
    [81] WANG HUI,QIN QINGHUA,A meshless method for generalized linear or nonlinear Possion-type prob.[J],Engineering Analysis with Boundary Elements.2006(30):515-521.
    [82] Zhang Y, Zhu S.On the choice of interpolation functions used in the dual-reciprocityboundary-element method[J]. Eng Anal Bound Elem 1994(13):387-96.
    [83] Zhu J.Integral equation solutions for finite and infinite plate[C].C. A. Brebbia,G. Maier.Boundary Elements VII.Berlin:Springer-Verlage,1985.
    [84] Beskos D E.Boundary element analysis of plates and shells [M].Berlin:Springer- Verlage, 1991.
    [85] Chantasiriwan S. Cartesian grid methods using radial basis functions for solving Poisson,Helmholtz,and diffusion-convection equations [J]. Engng Anal Bound Elem,2004(28):1417-1425.
    [86] Yamada PW,Wrobel LC ,Power H. On the convergence of the dual reciprocity method[J]. Eng Anal Bound Elem,1994(13):291-8.
    [87] Mukherjee Y X, Mukherjee S. The boundary node method for potential problems[J].Int.J. Num.Meth.Engng, 1997(40):797~815.
    [88] Xie H, Nogami T, Wang J.A radial boundary node method for two-dimensional elastic analysis [J].Engng Anal Bound Elem,2003(27):853-862.
    [90] Mau-Duy N,Tran-Cong T,Tanner R.I. A domain-type boundary–integral-equation method for 2-D biharmonic Dirichlet problem[J].Engng Anal Bound Elem,2006(30):809-817.

© 2004-2018 中国地质图书馆版权所有 京ICP备05064691号 京公网安备11010802017129号

地址:北京市海淀区学院路29号 邮编:100083

电话:办公室:(+86 10)66554848;文献借阅、咨询服务、科技查新:66554700