摘要
边界元法(BEM)是一种应用广泛的求解偏微分方程的方法,它具有精度高,降维等特点。无网格局部Petrov–Galerkin (MLPG)法是一种很受关注的数值方法,适合于求解非齐次,非线性,各向异性等问题。
本文首先将无网格局部Petrov-Galerkin法和改进的移动最小二乘近似相结合,形成了改进的无网格局部Petrov-Galerkin法,并求解了二维类Helmholtz方程。在无网格局部Petrov-Galerkin法中,移动最小二乘近似被用来构造近似函数,在移动最小二乘近似中的代数方程组有时是病态的。因此改进的移动最小二乘近似被提出,改进的移动最小二乘近似采用加权正交函数系作为基函数,与传统的移动最小二乘近似相比,改进的移动最小二乘近似中的系数矩阵变成了非奇异的对角矩阵,因而无需计算系数矩阵的逆。数值算例的研究结果均表明改进的无网格局部Petrov-Galerkin法精度高,收敛速度快。
无网格局部Petrov-Galerkin法是一个真正的无网格方法,它不需要单元或网格,但是它的计算量比有限元和边界元都大。因此本文基于边界元法和无网格局部Petrov-Galerkin法提出了一种直接耦合法,该方法将问题区域分解为不相重叠的边界元子域和无网格子域,连续性条件在两子域的公共边界上得到满足。然后将边界元方程、无网格方程以及连续性条件耦合成最终的方程组。在不同的子区域划分模式下讨论了该方法,一些数值算例被给出,证明了该方法的有效性。
耦合法需要将边界元方程和无网格方程联立在一起,形成一个统一的大型方程组,因此本文又研究了基于边界元法和无网格局部Petrov-Galerkin法的区域分解法,该方法也将问题区域分解为不相重叠的边界元子域和无网格子域,连续性条件要在两子域的公共边界上得到满足,必须通过迭代程序。为了加速收敛,引进了固定松弛因子和动态松弛因子。然后通过丰富的数值算例详细讨论了两种松弛因子对迭代次数的影响以及公共边界上的初始值对迭代次数的影响。
The boundary element method (BEM) is a widely used method because of its accuracy and dimensionality reduction. The meshless local Petrov–Galerkin (MLPG) method is a promising meshless method for solving partial differential equations. It is well-suited to problems involving non-homogeneous, anisotropic and non-linear problems.
Firstly, the MLPG method and improved MLS (IMLS) approximation are combined, and an improved MLPG (IMLPG) method for two-dimensional analogical Helmholtz equation is discussed in this paper. In the MLPG method, the MLS approximation is used to obtain the approximation function. The algebra equation system in MLS approximation is sometimes ill conditioned. So a new method to establish the approximation function, IMLS approximation, was presented. In IMLS approximation, an orthogonal function system with a weight function is used as the basis function. In comparison with MLS approximation, the algebra equation system in IMLS approximation is not ill conditioned, and can be solved without having to obtain the inverse matrix. In our numerical tests, the numerical convergence of IMLPG method is studied and accurate results are obtained.
Secondly, The MLPG method is a genuine meshless method which does not need‘‘elements’’or‘‘mesh,’’but uses a distributed set of nodes for both field interpolation and background integration. Although the MLPG method has many advantages, the method is more computationally expensive than the finite element method (FEM) and the boundary element method (BEM). So a coupled BEM and MLPG method for analyzing two-dimensional potential problems is presented in this paper. In this method, the analysis domain is divided into two non-overlapping regions (BEM region and MLPG region). The proposed coupling method directly couples the two methods without transition region, the continuity conditions are satisfied on the interface of the two sub-regions. The final system of algebra equation is composed of the BEM equation, the MLPG equation and the equation formed by the continuity conditions. This coupling method is numerically treated under several subdomain models. Some numerical examples for the potential problems governed by the Laplace and Poisson equations are presented to evaluate the accuracy and efficiency of the proposed technique.
Thirdly, the coupling method must employs an entire unified equation for the whole domain, by combining the discretized equations for the BEM and MLPG sub-domains, so a non-overlapping domain decomposition algorithm based on the BEM and MLPG method is presented in this paper. This algorithm is iterative in nature. It essentially involves subdivision of the problem domain into subregions being respectively modeled by the two methods, as well as restoration of the original problem with continuity and equilibrium being satisfied along the interface. To speed up the rate at which the algorithm converges, static and dynamic relaxation parameters are employed. The validity of the algorithm is verified by solving some potential problems.
引文
[1] C. A. Brebbia. The boundary element method for engineers[M].London: pentech press,1978.
[2]杨德全.赵忠生.边界元理论及应用[M].
[3]余德浩.自然边界元方法的数学理论[M].科学出版社.1993.
[4]祝家麟.椭圆边值问题的边界元分析[M].科学出版社.1991.
[5]严重.边界单元法基础[M],重庆,重庆大学出版社,1986.
[6]王元淳.边界元法基础[M].上海交通大学出版社.1988.9.
[7]申光宪着.边界元方法[M].机械工业出版社, 1998.
[8]董海云,祝家麟,张守贵.二维Laplace方程Dirichlet问题直接边界积分方程的Galerkin解法[J].重庆大学学报. 2006,29(4):122-125.
[9]袁政强,黄剑,祝家麟.三维Laplace方程边界元中线性单元的精确积分法[J].应用力学学报. 2004,21(3):117-121.
[10] W.T.Ang. A boundary integral equation method for the two-dimensional diffusion equation subject to a non-local condition[J]. Engineering Analysis with Boundary Elements.2005, 25:1-6.
[11]李开海,祝家麟.弹性动力学问题的边界元区域分解算法[J].数值计算与计算机应用. 2000,4:276-286.
[12] T.S.A.Ribeiro, G.Beer, C.Duenser. Efficient elastoplastic analysis with the boundary element method[J]. Comput Mech.2008,41:715-732.
[13] Yumin CHENG, Miaojuan PENG. Boundary element-free method for elastodynamics[J]. Science in China Ser. G Physics, Mechanics & Astronomy. 2005, 48 (6 ):641-657.
[14] C. A. Brebbia and L. C. Wrobel. Boundary element method for fluid flow[J]. Advances in Water Resources, 1979, 2:83-89.
[15] T. A. Cruse and S. T. Raveendra. A comparison of long and short crack elastoplastic response using the boundary element method Engineering Fracture Mechanics[J], 1988,30(1):59-75.
[16] G. C. Hsiao, B. N. Khoromskij and W. L. Wendland. Preconditioning for boundary element methods in domain decomposition[J]. Engineering Analysis with Boundary Elements, 2001,25(4-5): 323-338.
[17] G. C. Hsiao, O. Steinbach and W. L. Wendland. Domain decomposition methods via boundary integral equations[J]. Journal of Computational and Applied Mathematics, 2000,125(1-2): 521-537.
[18] W.L. Wendland and J. Zhu. The boundary element method for three-dimensional stokesflows exterior to an open surface[J]. Mathematical and Computer Modelling, 1991,15(6):19-41.
[19]杜庆华,姚振汉.边界积分方程-边界元法的基本理论及其在弹性力学方面的若干工程应用[J].固体力学学报, 1982,1:1-22.
[20]杜庆华,嵇醍.弹性与弹塑性边界元应力分析的若干最近成果[J].应用力学学报, 1985,2(2):1-15.
[21]姚振汉,杜庆华.边界元法应用的若干近期研究及国际新进展[J].清华大学学报, 2001, 1(4): 89-93.
[22]张雄,刘岩编着.无网格法[M].清华大学出版社,2004.
[23]刘桂荣,顾元通着.无网格法理论及程序设计[M].山东大学出版社,2007.
[24] R A Gingold, J J Monaghan. Smoothed particle hydrodynamics: theory and application to non-spherical stars[J]. Mon. Not. Roy. Astrou. Soc., 1977, 18:375-389.
[25] B Nayroles, G Touzot, P Villon. Generalizing the finite method: diffuse approximation and diffuse elements[J]. Comput. Mech., 1998, 21:91-100.
[26] T Belytschko, Y Y Lu, L Gu. Element free Galerkin methods[J]. Int. J. Num. Meth. Engng., 1994,37:229-256.
[27] WK Liu, S Jun, YF Zhang. Reproducing Kernel Particle methods[J]. Int.J.Numer.meth. Fluids., 1995,20:1081-1106.
[28] WK Liu, SF Li, T Belytsehko. Moving least-square reproducing kernel methods .1. Methodology and convergence[J] . Comput. Meth.Appl.Mech Eng., 1997, 143:113-154.
[29] S Li, WK Liu. Reproducing Kernel Hierarchical Partition of Unity Part l: Theoretical Formulation[J]. Int.J.Numer.methods Eng.,1999,45:251-288.
[30] WKLiu, WM Han, HS Lu,etal. Reproducing kernel Particle method Part l: Theoretical formulation[J]. Comput. Meth.Appl.Mech Eng.,2004,193:933-951.
[31] JS Chen, C Pan, CT Wu. Large Deformation Analysis of Rubber Based on a Reproducing Kernel Particle Method[J]. Comput.meth.,1997,19:211-217.
[32] JS vChen, CMOL Roque, CH pan, etal. Analysis of Metal Forming Process Based on Meshless Method[J]. J.Mater.Proeess.Technol.,1998.80-81:642-646.
[33] RA Uras, CT Chang, Y Chen, WK Liu. Multi resolution Reproducing Kernel Particle Methods in Acoustics[J]. J. Comput. Acoust., 1997,5(l):71-94.
[34] SF Li, WK Liu, AJ Rosakis,etal. Mesh-free Galerkin simulations of dynamic shear band Propagation and failure mode transition[J]. Int. J. Solids Struct.,2002,39(5):1213-1240.
[35] I Babuska, JM Melenk. The partition of unity method[J]. Int. J. Numer. Methods Eng.,1997,
[36] JM Melenk, I Babuska, The partition of unity finite element method: Basic theory and applications[J]. Comput. Meth.Appl.Mech Eng.,1996,139:289-314.
[37] T Strouboulis, I Babuska, KL Copps. The design and analysis of the Generalized Finite Element Method[J]. Comput. Meth.Appl.Mech Eng.,2000,181:43-69.
[38] CA Duarte, JT Oden. An h-P adaptive method using louds[J]. Comput. Meth.Appl.Mech Eng., 1996,139(l-4):237-262.
[39]刘欣,朱德憋,陆明万,张雄.平面裂纹问题的h, p, hp型自适应无网格方法的研究[J]..力学学报. 2000,32(3):308-318.
[40] E Onate, S Idesohn, O C Zienkiewicz. A finite point method in computational mechanics: Applications to convective transport and fluid flow[J]. Int. J. Num. Meth. Engng., 1996, 39: 3839-3866.
[41] E Onate, S Idelsohn, OC Zienkiewiez,etal. A stabilized finite point method for analysis of fluidMechanics problems[J]. Comput. Meth.Appl.Mech Eng., 1996,139(l-4):315-346.
[42] X Zhang, XH Liu, KZ Song,etal. Least-squares collocation meshless method[J]. Int. J. Num. Methods Eng.,2001,51(9):1089-1100.
[43] XF Pan, X Zhang, MW Lu. Meshless Galerkin least-squares method[J].Comput.meth.,2005, 35(3):182-189.
[44] S N Atluri, T Zhu. A new Meshless Local Petrov-Galerkin (MLPG) approach in computational mechanics[J]. Comput. Mech., 1998,22:117-127.
[45] S N Atluri, J Y Cho, H G Kim. Analysis of thin beams, using the meshless local Petrov-Galerkin method, with generalized moving least squares interpolations[J]. Comput. Mech., 1999,24:334-347.
[46] T Zhu, J D Zhang, S N Atluri. A local boundary integral equation (LBIE) method in computational mechanics, and a meshless discretization approach[J]. Comput. Mech., 1998,21:223-235.
[47]蔡永昌,朱合华,王建华.基于Voronoi结构的无网格局部Petrov-Galerkin方法[J].力学学报,2003,35(2):187-193.
[48] SA Viana, D Rodger, HC Lai. Meshless local Petrov-Galerkin method with radial basis functions applied to electromagnetics[J]. IEE Proc.-Sci. Meas.Technol.,2004,151(6):449-451.
[49] J Braun. A numerical method for solving Partial differential equation s on highly irregular revolving grids[J]. Nature,1995,76:655.
[50] N Sukumar, B Moran, T Belytsehko. The natural element method in solid mechanics[J]. Int. J Numer. Methods Eng., 1998,43(5):839-887.
[51] GR Liu, YT Gu. A Point interpolation method for two-dimensional solids[J]. Int. J Numer.Methods Eng., 2001,50(4):937-951.
[52] YX Mukherjee, S Mukherjee. The boundary node method for Potential Problems, Int. J Numer[J]. Methods Eng., 1997,40(5):797-815.
[53] MK Chati, S Mukherjee.The boundary node method for three dimensional Problems in Potential Theory[J]. Int. J Numer. Methods Eng., 2000,47(9):1523-1547.
[54] G Li, NR Aluru. A boundary cloud method with a cloud-by-cloud Polynomial basis[J]. Eng. Anal. Bound.Elem.,2003,27(1):57-71.
[55] LC Nieolazzi, CS Bareellos, EA Faneello,etal. Generalized boundary element method for galerkin Boundary integrals[J]. Eng. Anal. Bound.Elem.,2005,29(5):494-510.
[56] J M Zhang, ZH Yao, H Li. A hybrid boundary node method[J]. Int. J Numer. Methods Eng., 2002, 5(4):751-763.
[57] JM Zhang, ZH Yao. Meshless regular hybrid Boundary Node Method[J]. CMES-Comp. Model, Eng.Sci.,2001,2(3):307-318.
[58] Zhang JM, Tanaka M, Endo M. The hybrid boundary node method accelerated by fast multipole expansion technique for3D potential problems[J]. Int. J Numer. Methods Eng., 2005,63(5):660-65.
[59] JM Zhang, ZH Yao, M Tanaka. The meshless regular hybrid boundary node method for 2D line elasticity[J]. Eng. Anal. Bound.Elem.,2003,27(3):259-268.
[60] W Chen, YC Hon. Numerical investigation on convergence of boundary knot method in the Analysis of homogeneous Helmholtz,modified Helmholtz,and convection-diffusion Problems[J].Comput. Meth.Appl.Mech Eng., 2003,192(15):1859,1875.
[61] W Chen. Meshfree boundary Particle method applied to Helmholtz Problems[J]. Eng. Anal. Bound.Elem.,2002,26(7):577,581.
[62] O.C. Zienkiewicz, D.W. Kelly and P. Bettes, The coupling of the finite element method and boundary solution procedures[J]. Int J Numer Methods Engng, 1977,11:355–375.
[63] C.A. Brebbia and P. Georgion, Combination of boundary and finite elements in elastostatics[J]. Appl Math Model 1979,3: 212–220.
[64] G. Beer and J.L. Meek, The coupling of boundary and finite element methods for infinite domain problems in elasticity[J]. In: C.A. Brebbia, Editor, Boundary element methods, Springer, Berlin 1981 :575–591.
[65] G. Swoboda, W. Mertz and G. Beer, Rheological analysis of tunnel excavations by means of coupled finite element (FEM)–boundary element (BEM) analysis[J]. Int J Numer Anal Methods Geomech ,1987,11:115–129.
[66] H.B. Li, G.M. Han, H.A. Mang and P. Torzicky, A new method for the coupling of finiteelement and boundary element discretized sub-domains of elastic bodies[J]. Comput Methods Appl Mech Engng ,1986,54:161–185.
[67] K.L. Leung, P.B. Zavareh and D.E. Beskos, 2-D elastostatic analysis by a symmetric BEM/FEM scheme[J]. Engng Anal Boundary Elem ,1995,15: 67–78.
[68] A.S. Chen, G. Hofstetter, Z.K. Li, H.A. Mang and P. Torzicky, Coupling of FE- and BE-discretizations for 3D-stress analysis of tunnels in layered anisotropic rock. In: Discretization methods in structural mechanics[J]. Springer, Berlin ,1989:427–436.
[69] N. Kamiya, H. Iwase and E. Kita, Parallel computing for the combination method of BEM and FEM[J]. Engng Anal Boundary Elem ,1996,18:221–229.
[70] N. Kamiya and H. Iwase, BEM and FEM combination parallel analysis using conjugate gradient and condensation[J]. Engng Anal Boundary Elem ,1997,20:319–326.
[71] W M Elleithy, H J Al-Gahtani, M El-Gebeily. Iterative coupling of BE and FE methods in elastostatics[J] . Engineering Analysis with Boundary Elements 2001,25 (8):685-695.
[72] Jr D Soares, WJ Mansur,O Von Estorff. An efficient time-domain FEM/BEM coupling approach based on FEM implicit Green's functions and truncation of BEM time convolution process[J]. Computer Methods in Applied Mechanics and Engineering 2007,196 (9-12): 1816-1826.
[73] T. Belytschko, D Organ. A coupled finite element–element-free Galerkin method[J]. Comput Mech 1995,17:186–195.
[74] D. Hegen, Element-free Galerkin methods in combination with finite element approaches[J]. Comput Meth Appl Mech. Engng . 1996,135:143–166.
[75] G.R. Liu and Y.T. Gu, Meshfree local Petrov–Galerkin (MLPG) method in combination with finite element and boundary element approaches[J]. Comput Mech 2000,26 (6):536–546.
[76] Q Z Xiao , M D hanasekar. Coupling of FE and EFG using collocation approach[J]. Advances in Engineering Software, 2002,33:507-515.
[77] L Jia, Q Huang, Z Yin. New and better EFGM-FEM direct coupling method[J]. Journal of Northwestern Polytechnical University 2007,25 (3): 337-341.
[78] S Li, B Sui, W Zhu. Coupling meshless method and finite element method for seepage of jointed rock mass[J]. Chinese Journal of Rock Mechanics and Engineering 2007,26 (1):75-80.
[79] Gu, Y.T., Zhang, L.C. Coupling of the meshfree and finite element methods for determination of the crack tip fields[J]. Engineering Fracture Mechanics 2008,75 (5):986-1004.
[80] Y T Gu, G R Liu. A coupled element-free Galerkin/boundary element method for stress analysis of two-dimension so1id[J]. Computer Methods in Applied Mechanics andEngineering, 2001, 190:4405-4419.
[81] Y T Gu, G R Liu. Meshless Methods Coupled with Other Numerical Methods[J]. Tsinghua Science and technology, 2005,10(1):8-15.
[82] Zan Zhang, K M Liew, Yumin Cheng. Coupling of the improved element-free Galerkin and boundary element methods for two-dimensional elasticity problems[J]. Engineering Analysis with Boundary Elements. 2008,32:100-107.
[83] Y T Gu, G R Liu. Coupling of element free Galerkin and hybrid boundary element methods using modified variational formulation[J]. Computational Mechanics, 2000,26:166-173.
[84] C K Lee, S T Lie, Y Y Shuai. On coupling of reproducing kernel particle method and boundary element method[J]. Computational Mechanics, 2004,34:282-297.
[85] Gang Li, Glaucio H Paulino, N R Aluru. Coupling of the mesh-free finite cloud method with the boundary element method a collocation approach[J]. Comput. Methods Appl. Mech. Engng., 2003,192:2355-2375.
[86] J Euripides. Sellountos, Demothenes Polyzos. A MLPG(LBIE) approach in combination with BEM[J]. Comput. Methods Appl. Mech. Engrg. 2005,194: 859–875.
[87]吕涛.石济民.林振宝.区域分解算法---偏微分方程数值解新技术[M].科学出版社,1992.
[88] NS Bakhvalov.and GM Kobel’kov. 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,NO.190,1988.
[89] Kuznetsov,YU.A.,Multilevel domain decomposition methods[J].Applied Numer. Math., 1990,6: 303--314.
[90] JH Bramble, JE Pasciak, XU J. Parallel multilevel Preconditioners[J]. Muth.Comp., 1990, 55:1-22.
[91] WM Elleithy and HJ Al-Gahtani, An overlapping domain decomposition approach for coupling the finite and boundary element methods[J]. Engrg. Anal. Boundary Elem. 2000,24(5) :391–398.
[92] WM Elleithy, HJ Al-Gahtani, M El-Gebeily, Convergence of the iterative coupling of BEM and FEM, in: 21st World Conference on the Boundary Element Method, BEM21, Oxford University, UK, 1999:281–290.
[93] WM Elleithy, M Tanaka. Interface relaxation algorithms for BEM-BEM coupling and FEM-BEM coupling[J]. Computer Methods in Applied Mechanics and Engineering, 2003, 192(26-27): 2977-2992.
[94] W M Elleithy, H J Al-Gahtani and M El-Gebeily, Iterative coupling of BE and FE methods in elastostatics[J]. Engrg. Anal. Boundary Elem. 2001,258: 685–695.
[95] N. Kamiya, H. Iwase and E. Kita, Parallel implementation of boundary element method with domain decomposition[J]. Engrg. Anal. Boundary Elem. 1997,18:209–216.
[96]张太平,祝家麟.抛物型问题的边界元重叠型区域分解法[J].重庆大学学报,2002,2(25): 75-78.
[97]袁政强,谭宏,祝家麟. Stokes方程的边界元区域分裂算法及其应用[J].重庆大学学报,2001,4(24):71-73.
[98]田玲玲.地下水渗流的直接边界元非重叠性区域分解算法[D].重庆大学硕士学位论文,2006.
[99] F Ihlenburg, I Babuska. Finite element solution of the Helmholtz equation with high wave number Part I: The h-version of the FEM[J]. Comput. Math. Appl., 1995, 30(9):9–37.
[100] I Singer, E Turkel. High-order finite difference methods for the Helmholtz equation[J]. Computer methods in Applied Mechanics and Engineering, 1998, 163(1-4):343-358.
[101] T.Chen I.S.Raju. A coupled finite element method and meshless local Petrov-Galerkin method for two-dimensional potential problems[J]. Comput. Methods Appl. Mech. Engrg. 2005,194: 859–875.
[102] Chia-Ching Lin, E.C.Lawton, J.A.Caliendo and L.R.Anderson. An iterative finite element-boundary element algorithm[J]. Computers&Structures,1996,59(5):899-909.