无单元法解拉普拉斯方程和泊松方程(hanshu(.pdfVIP

第 21 卷第 1 期 青岛建筑工程学院学报 Vol. 21 №1 2000 Journal of Qingdao Institute 　 of Architecture and Engineering 无单元法解拉普拉斯方程和泊松方程① 庞　辉 张伟星 (青岛建筑工程学院土木工程系 ,青岛 266033) (青岛建筑工程学院教务处) 摘 　要 　将无单元法理论与微分方程等价性结合 ,用来求解拉普拉斯方程和泊松方程 ,并编制相 应的计算程序进行计算. 关键词 　无单元法 ,边值问题 ,微分方程等价性 中图分类号 　O241 有限元发展至今 , 已经成为工程数值分析的有力工具. 有限元法的成功之处在于场函数的 局部化近似 ,可以处理相当复杂的问题. 有限元法在如下问题上仍值得改进 : ( ) 1 在有限元分析中计算网络的产生 ,单元信息输入一直是个难题 , 目前虽有许多前后处 理软件 ,但其功能有限. (2) 通常有限元都采用 C0 阶连续性单元, 虽然可构造出 C1 、C2 阶连续性的近似函数, 但 需增加节点参数, 增加计算复杂性. 无单元法采用滑动最小二乘法所产生的光滑函数近似场函数, 保留了有限元的一些特点, 但克服了有限元的上述不足. 无单元法的思想最早由 Nayroles ,etal [1 ] 于 1992 年提出 ,Belytschko ,et al [2 ] 于 1994 年对 ( ) Nayroles 的方法作了改进 ,提出无单元伽辽金法 Element - free Galerkin Method ,之后 Lu Y , Y et al [3 ]又将无单元伽辽金法作了进一步的改进. 1998 年 ,周维垣[4 ] 对其在平面弹性体中的 应用作了研究. 笔者将无单元法与微分方程的等价性[5 ] 相结合 ,用来求解拉普拉斯和泊松方程 ,计算表 明 ,此法不仅减少了人工处理工作量 ,而且结果的精度也是令人满意的. 1 　无单元法的插值函数 根据滑动最小二乘法 ,场函数 U ( X) 可被近似为 : m h ( ) ( ) ( ) T ( ) ( ) ( ) U X = P X a X = P X a X 1 ∑j j j = 1 其中 : X = ( x , y) T 为场点 ; a ( X) 为 m 维系数向量; P ( X) 为 m 维基向量. 对于二维情况, P ( X) 可作如下选取 ① 收稿日期 :1999 —04 —21 © 1994-2006 China Academic Journal Electronic Publishing House. All rights reserved. 100 　 　 　 　 　 　 　 　青 岛 建 筑 工 程 学 院 学 报 　　　　　　　　第 21 卷 ( ) ( ) T ( ) P X = 1 , x , y 　　　　 线形基 m