用户名: 密码: 验证码:
开边界外调和方程的迦辽金边界元解法
详细信息    本馆镜像全文|  推荐本文 |  |   获取CNKI官网全文
摘要
最典型也最简单的椭圆型偏微分方程是调和方程,又称Laplace方程。力学和物理学研究中的许多问题都归结为Laplace方程的边值问题。例如:弹性膜的平衡问题,稳定状态的热传导问题,不可压缩势流问题,静电场问题以及静磁场问题。这些问题虽然有完全不同的物理背景,却往往导致完全相同的数学表达式。本文利用Galerkin边界元方法对该问题加以研究,特别是针对边界是直线段或开弧段情况下无限域上的调和方程进行了研究。该问题属奇异边界问题,对应于裂纹、屏障等问题的数值模拟在实际应用的背景。因为把二维调和方程的边值问题转化为等价的边界积分方程时带有约束条件,用Galerkin边界元方法求解带约束条件的变分方程,在数值离散时,我们采用Lagrange乘子法处理约束条件。
    Galerkin方法是基于变分原理基础上的一种把微分方程或积分方程转化为等价的变分方程,通过离散变分方程求原方程数值解的方法。本文把Laplace方程的边值问题转化为边界积分方程后,通过与边界积分方程等价的变分形式,采用线性单元,利用Galerkin边界元方法求解。在计算单元刚度矩阵时,对二重积分的第一重使用精确积分,第二重使用数值积分,从而有效的克服了奇异积分的计算。本文还把此方法推广到求解区域边界是直线段或者开弧段的情况,采用引入奇异边界单元的方法来模拟解在开边界端点附近的奇异性。最后利用Fortran Power Station4.0程序语言分别编写了用Galerkin方法求解区域边界为闭曲线及区域边界为直线段或开弧段Laplace方程的计算机程序,通过几个算例证明了该方法是有效的、是可行。而且通过数值实验也验证了Galerkin方法误差的理论结果。
The harmonic equation,or Laplace equation is a typical and simple elliptic partial differential equation. The mathematical model of many problems of mechanics and physics are reduced to boundary problems of Laplace equation. Such as: the balance of elastic membrane, the heat conduction of stationary state, the incompressible potential flow, the problem of electrostatic field and magnetostatics field. These problems have different physical nature, but can be reduced to the same mathematical expression. In this paper, we study the numerical solutions of these problems using Galerkin boundary element, especially for the problem with open boundary such as the problems exterior to an open segment or an open curve in the plane. Since the equivalent boundary integral equation for two dimensional Laplace equation has constraint condition. Lagrange multiplicator method is introduced in the numerical computation to release the constraint.
     Galerkin method based on the variation principle is used to solve differential and integral equations. In this paper, the boundary problem of Laplace equation is changed into the variational equation which is equivalent to the boundary integral equation. Using linear element, it is solved by Galerkin boundary element method. In computation of stiffness matrix, the calculation of double singular integration is needed. The exactly integral formula is used in the inner integral expression, and the numerical integral formula is used in the outer integral expression. We extend this method to the problem of a domain which boundary is an open arc or an open segment. In this paper, we simulate the singularity of solution in extreme point of open arc or open segment by singular element. With Fortran Power Station 4.0, we make Galerkin boundary element method program for solving Laplace equation on region which boundary is a closed curve or an open arc, and the numerical experiments also prove this method is reliable. Last we test the error of Galerkin boundary element by numerical experimentation.
引文
[1] C.A. 布瑞比亚,J.C.F 泰勒斯,L.C. 诺贝尔著, 边界元法的理论和工程应用,龙述兴,刘腾喜,蔡松柏译.国防工业出版社,1988.
    [2] C.A.Brebbia著,边界单元法进展,陈祥福,王家林译,中国展望出版社,1986.
    [3] 杜庆华等合著,边界积分方程法--边界元法,高等教育出版社,1989.
    [4] 王元淳著,边界元法基础,上海,上海交通大学出版社,1988.
    [5] 余德浩著,自然边界元法的数学理论,科学出版社,1993.
    [6] 严重著,边界单元法基础,重庆,重庆大学出版社,1986.
    [7] 姚寿广编著,边界元数值方法及其工程应用,国防工业出版社,1995.
    [8] 祝家麟著,椭圆边值问题的边界元分析,科学出版社,1991.
    [9] 祝家麟,边界元方法--一个老而新的发展中的数值方法.数学的实践与认识,1983,第三期.
    [10] 祝家麟,边界元方法--它的产生,内容和意义,第一届工程中的边界元法会议论文集,杜庆华,1985.
    [11] 刘希云,赵润祥编著,流体力学中的有限元与边界元方法,上海交通大学出版社,1993.
    [12] 李瑞遐,边界是光滑开弧Helmholtz方程的边界积分法,华东理工大学学报,25:4(1999),410-412.
    [13] 李瑞遐,外边值问题的边界元法与有限元法组合及奇性处理,应用数学与计算数学学报,6:1(1992).34-41.
    [14] 刘光廷,王宗敏,周鸿钧,缝端奇异边界单元和界面裂缝的应力强度因子计算,清华大学学报,36:1(1996),34-40.
    [15] 韩宇,单辉祖,正交各向异性平面断裂问题应力强度因子的边界元解法,计算结构力学及其应用,8:4(1991),373-377.
    [16] 祝家麟,边界元方法中的奇异性,重庆建筑工程学院学报,13:2(1991),90-102.
    [17] 李瑞遐,Helmholtz方程外边值问题的自然边界元法,高校应用数学学报(A辑),12:3(1997),369-373.
    [18] 陆志良,杨柞生,不可压粘流N-S方程的边界积分解法,力学学报,28:2(1996),225-231.
    [19] 刘希云,杨柞生,广义Stokes方程的完全边界积分表示式及其在求解N-S方程中的应用,力学学报,24:6(1992),645-651.
    [20] 祝家麟,用边界积分方程法解平面双调和方程的Dirichlet问题,计算数学,3(1984),278-288.
    [21] 祝家麟,三维定常流Stokes问题的边界积分方程法,计算数学,1(1985),40-49.
    [22] 祝家麟,金朝嵩,带滑动边界条件的Stokes方程的边界元逼近,应用数学学报,14:3(1991),296-303.
    
    
    [23] 祝家麟,定常Stokes问题的边界积分方程法,计算数学,3(1986),281-289.
    [24] 李瑞遐,曲线积分的误差分析,华东理工大学学报,21:2(1995),267-272.
    [25] 蒋劲,翁晓红,边界积分方程中强奇异积分的解法,武汉水利电力大学学报,29:6(1996).
    [26] 蒋礼尚,庞之垣著,有限元方法及其理论基础,人民教育出版社,1979.
    [27] 李开海,弹性动力学方程的边界元重叠型区域分解法,重庆建筑大学硕士论文,1999.
    [28] 张太平,抛物型方程的边界元区域分解法,重庆大学硕士论文,2001.
    [29] 李瑞遐编著,有限元法与边界元法,上海科技教育出版社,1993.
    [30] 孙炳楠,项玉寅,张永元编,工程中边界元法及其应用,浙江大学出版社,1991.
    [31] 卢盛松主编,边界元理论及应用,高等教育出版社,1990.
    [32] 冈村弘之(日)著,线性断裂力学入门,李顺林译,江苏科学技术出版社,1981.
    [33] 嵇醒,臧跃龙,程玉民编著,边界元法进展及通用程序(M),同济大学出版社,1997.
    [34] C. A. 布瑞比亚, S. 沃克著,边界元法的工程应用,张治强,周天孝,陈瀚译,陕西科学出版社,1985.
    [35] 沃德·切尼,戴维·金凯德著,数值数学和计算,薛密译,复旦大学出版社,1991.
    [36] Germund, Dahlqist, Ake Bjorck(瑞典)著,数值方法,包雪松译,高等教育出版社,1990.
    [37] 蔡大用,白峰杉,高等数值分析,清华大学出版社,1997.
    [38] 沐定夷,胡鸿钊编著,数值分析,上海交通大学出版社,1994.
    [39] 申光宪,肖宏,陈一鸣编著,边界元法,机械工业出版社,1998.
    [40] 田宗若著,复合材料中的边界元法,西北工业大学出版社,1992.
    [41] 臧跃龙,嵇醒,关于边界元法中奇异积分的处理,固体力学学报,15:2(1994),161-165.
    [42] 汪鸿振,郭芃,用边界元法计算声辐射时高次奇异积分的处理方法,声学技术,15:3(1996),97-100.
    [43] 吴凤林,任家骏,边界元法中奇异积分核的一种有效数值方法,太原理工大学学报,29:4(1998),373-375.
    [44] 张效松,叶天麒,二维边界元奇异积分和多域缩聚法分析,计算力学学报,14:2(1997),196-203.
    [45] P. K. 班努杰,R. 白脱费尔德编著,工程科学中的边界元法,冯振兴,李正秀,陶维本译,国防工业出版社,1988.
    [46] 余德浩,计算数学与科学工程计算及其在中国的若干发展,数学进展,31:1(2002),1-6.
    [47] 张耀明,平面Poisson外边值问题,应用力学学报,19:1(2002),39-43.
    [48] 王小贞,臧跃龙,粘性液体小幅晃动常单元及线性单元的比较分析,重庆建筑大学学报,22:6(2000),67-69.
    [49] 唐少武,冯振兴,李正秀,高阶奇异边界积分的新型求积公式,重庆建筑大学学报,22:6(2000),54-57.
    
    
    [50] 徐明编,Fortran PowerStation 4.0基础教程,清华大学出版社,2000.
    [51] 高福成,姜玉泉,梁静毅编著,最新FORTRAN程序设计实用教程,天津科学技术出版社,1993.
    [52] 张有天,王镭,陈平,边界元方法及其在工程中的应用,水利电力出版社,1989.
    [53] 黎在良,王元汉,李廷芥著,断裂力学中的边界数值方法,地震出版社,1996.
    [54] 朱谷君主编,工程传热传质学,航空工业出版社,1989.
    [55] 《现代数学手册》数学编纂委员会,现代数学手册(计算机数学卷),华中科技大学出版社,2001.
    [56] 《数学手册》编写组,数学手册,高等教育出版社,1979.
    [57] J. P. 霍尔曼著,姚维信,陈宏模译,热传递学问题详解(上、下册),晓园出版社,1993.
    [58] 邝国能,熊振南,宋振熊编著,工程实用边界单元法,中国铁道出版社,1989.
    [59] 饶寿期编,有限元法和边界元法基础,北京航空航天大学出版社,1990.
    [60] 厉家尚,陆大有,王光麟,于广经,陆寿娥编译,传热基础600题解,宇航出版社,1990.
    [61] D. R. 匹茨,L. E. 西逊姆编著,夏雅君译,传热学的理论与习题,机械工业出版社,1983.
    [62] 田中正隆,田中喜九昭著,郎德宏译,边界元法的基础与应用,煤炭工业出版社,1987.
    [63] 陆宁编著,MATLAB语言即学即会,机械工业出版社,2000.
    [64] 导向科技编著,MATLAB6.0程序设计与实例应用,中国铁道出版社,2001.
    [65] 清源计算机工作室编著,MATLAB高级应用——图形及影像处理,机械工业出版社,2000.
    [66] C. F. 杰拉尔德著,颜仁鸿译,应用数值分析,晓图出版社,世界图书出版公司,1993.
    [67] 赵亚楠,张晓峰,刘应华,岑章志,迦辽金边界元法在平面问题极限分析中的应用,重庆建筑大学学报,22:6(2000),112-116.
    [68] 蒋豪贤,梁峰,何志伟,陈学辉,二维边界元法的新算法,华南理工大学学报,24:1(1996),131-137.
    [69] 邢彭龄,木德,温度场的边界元分析,内蒙古民族师院学报,9:2(1994),32-35.
    [70] 丁立,袁修干,人体温度场的边界元分析,北京航空航天大学学报,26:6(2000),673-675.
    [71] A. R. Kukreti, M. M. Zaman, A. Issa, Analysis of fluid storage tanks including foundation-superstructure interaction, Appl. Math. Modelling, 17: (1993), 618-631.
    [72] Andreas Karageorghis, Graeme Fairweather, The simple layer potential method of fundamental solutions for certain biharmonic problems, International Journal for Numerical methods in Fluids, 9: (1989), 1221-1234.
    [73] C.A.Brebbia, the Boundary Element Method for Engineers , Pentech Press ,London, 1978.
    [74] C. A. Brebbia, Boundary element methods in engineering, Proceeding of the fourth international seminar, Southampton, England, 1982.
    [75] Cho Lik Chan, Abhijit Chandra, An algorithm for handling corners in the boundary element
    
    method: Application to conduction-convection equations, Appl. Math. Modelling, 15: (1991), 244-255.
    [76] E. S. Shoukralla, Numerical solution of Helmholtz equation for an open boundary in space, Appl. Math. Modelling, 21: (1997), 231-232.
    [77] George E.Blandford, Anthony R.Ingraffea, James A.Liggett, Two-dimensional stress intensity factor computations using the boundary element method, International Journal for Numerical Methods in Engineering,17:(1981), 387-404.
    [78] G. S. Padhi, R. A. Shenoi, S. S. J. Moy, M. A. McCarthy, Analytic integration of kernel shape function product integrals in the boundary element method, Computers and Structures, 79(2001),1325-1333.
    [79] Hajime Igarashi, Toshihisa Honma, A boundary element method for potential fields with corner singularities, Appl. Math. Modelling, 20: (1996), 843-852.
    [80] Henry Power, Guillermo Mirand, Second kind integral formulation of stokes flows past a particle of arbitrary shape, SIAM J. Appl. Math. ,47:(1987),689-698.
    [81] I. H. Sloan, A. Spence, The Galerkin method for integral equations of the first kind with Logarithmic Kernel: Theory, IMA Journal of Numerical Analysis, 8: (1988), 105-122.
    [82] I. H. Sloan, A. Spence, The Galerkin method for integral equations of the first kind with Logarithmic Kernel: Application, IMA Journal of Numerical Analysis, 8: (1988), 123-140.
    [83] I-Lin, Ming-Te Liang, Shyh-Rong Kuo, Jeng-Tzong Chen, Dual boundary integral equations for helmholtz equation at a corner using contour approach around singularity, Journal of marine science and technology, 1: (2001), 53-63.
    [84] Jean-Claude Nedelec, Finite elements for exterior problems using integral equations, International Journal for Numerical methods in Fluids, 7: (1987), 1229-1234.
    [85] Jurgen Friedrich, A linear analytical boundary element method(BEM) for 2D homogeneous potential problems, Computers & Geosciences, 28(2002),679-692.
    [86] KENDALL E. ATKINSION and IAN. H. SLOAN, Numerical solution of first-kind Logarithmic-kernel integral equations on smooth open arcs, mathematics of computation 193:(1991), 119-139.
    [87] Luigi Morino, Giovanni Bernardini, Singularities in BIEs for the Laplace equation; Joukowski trailing-edge conjecture revisited, Engineering Analysis with Boundary Elements, 25(2001), 805-818.
    [88] Lei Gu, A compatible boundary element method for plane elasticity and fracture mechanics, Appl. Math. Modelling, 17: (1993), 394-405.
    [89] Li Kaihai, Zhu Jialin ,Boundary element methods with domain decomposition method for
    
    elastodynamics, Chinese J. Num. Method & Comp. Appl. 23:2(2001), 43-55.
    [90] Martin Costable , Ernst Stephan, Curvature terms in asymptotic expansion for solutions of boundary integral equations on curved polygons, Journal of Integral Equations, 5(1983), 353-371.
    [91] Masataka Tanaka, New crack elements for boundary element analysis of elastostatics considering arbitrary stress singularities, Appl. Math. Modelling, 11: (1987), 357-363.
    [92] M. Costable, V. J. Ervin, E. P. Stephan, On the convergence of collocation methods for symm's integral equation on open curves, Mathematics of Computation, 51: 183(1988), 167-179.
    [93] M. Bourlard, S. Nicaise, L. Paquet, An adapted boundary element method for the Dirichlet problem in polygonal domains, SIMA J. Numer. Anal. 28: 3(1991), 728-743.
    [94] Nengquan Liu, Nicholas J. Altiero, A new boundary element method for the solution of plane steady-state thermoelastic fracture mechanics problems, Appl. Math. Modelling, 16: (1992), 618-629.
    [95] N. K. Mukhopadhyay, S. K. Maiti, A. Kakodkar, Variable singularity boundary element and its applications in computation of SIFs, Computers and Structures, 77(2000), 141-154.
    [96] Mirela Kohr, An indirect boundary integral method for a Stokes flow problem, Computer methods in applied mechanics and engineering, 190(2000), 487-497
    [97] P. Kunkel, Efficient computation of singular points, IMA Journal of Numerical Analysis, 9: (1989), 421-433.
    [98] Qing Zhang, Subrata Mukherjee, Design sensitivity coefficients for linear elastic bodies with zones and corners by the derivative boundary element method, Int. J. Solids Structure, 27: 8(1991), 983-998.
    [99] R.D. Henshell, K.G. Shaw, Crack tip finite elements are unnecessary, International Journal for Numerical Methods in Engineering, 9:(1975), 495-507.
    [100] Roshdy S.Barsoum, On the use of isoparametric finite elements in linear fracture mechanics, International Journal for Numerical Methods in Engineering, 10:(1976), 25-37
    [101] Ruixia Li, BEM and estimate of actual error for exterior Dirichlet problem for two-dimensional Helmholtz equation, Advances in Engineering Software , 15(1992), 181-183.
    [101] S.Kobayashi, N.Nishimura and T.Kawakami, Simple layer potential for domains having external corners, Appl.Math.Modelling, 8(1984), 61-65.
    [102] Stephan E.P, Wendland W.L, An augmented Galerkin procedure for the boundary integral method applied to two-dimensional screen and crack problem. Applicable Analysis, 18(1984)183-219.
    [103] S. Mukherjee, Finite parts of singular and hypersingular integrals with irregular boundary
    
    source points, Engineering Analysis with Boundary Elements, 24(2000), 767-776.
    [104] Trevor W. Dawson, On the singularity of the axially symmetric Helmholtz Green's function ,with application to BEM, Appl. Math. Modelling, 19: (1995), 590-600.
    [105] W.L.Wendland, J.Zhu, The boundary element method for three-dimensional stokes flows exterior to an open surface, Mathl.Comput.Modelling, 15:6(!991), 19-41.
    [106] Z. C. Li, T. T. Lu, Singularity and Treatments of Elliptic Boundary Value Problem, Mathematical and Computer Modelling, 31(2000), 97-145.
    [107] Zhenhan Yao, M. H. Aliabadi, Boundary Element Techniques, Tsinghua University Press, 2002.

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

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

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