三维无单元伽辽金法及其在拱坝的分析中的应用.pdfVIP

水 利 学 报 2005 年6 月 SHUILI XUEBAO 第36 卷 第6 期 文章编号: 0559-9350(2005)06-0644-06 三维无单元伽辽金法及其在拱坝分析中的应用 周维垣，黄岩松，林鹏 (清华大学 水利水电工程系，北京 100084) 摘要：本文介绍了三维无单元伽辽金法的基本理论，并将其应用于拱坝结构分析中，基于变分原理推导出三维线 弹性无单元伽辽金法的整体平衡方程。为了处理实际水利工程中断层、宏观裂纹等不连续情况，采用可视准则处 理不连续面对积分点影响域的影响，同时考虑不连续面两表面相互作用对整体平衡方程的贡献。编制了三维无单 元法计算程序并用以分析拉西瓦水电工程拱坝-地基系统模型，结果表明，其反映了拱坝在水荷载作用下的位移分 布规律。 关键词：三维无单元伽辽金法；滑动最小二乘法；三维不连续面；拱坝分析 中图分类号：TV642.2 文献标识码：A 无网格法摆脱单元的束缚，只采用节点信息及其紧支域上的权函数实行局部精确逼近，通过配点法或 是伽辽金法得到微分方程弱形式，再选用合适的积分方案聚合整体平衡方程并求解，从而实现对问题的求 解。无网格法能够有效地克服有限元法的一些问题，例如：(1)前处理困难，尤其对三维问题；(2)计算结 果通常是位移连续，而应力、应变不连续，需要进行后处理的修匀；(3)在材料接近不可压缩时，计算收 敛缓慢，计算结果失真，产生“体积自锁”现象；(4)在处理如裂纹扩展或是结构大变形等网格发生改变 的情况时存在困难。无网格法最早由Nayroles[1]等提出，称之为模糊单元法(Diffuse Element Method)， 此后各种无网格法不断涌现，至今已有十多种，如模糊单元法(DEM)、无单元伽辽金法(EFGM)、无网格局 部伽辽金法(MLPGM)、小波伽辽金法(WGM)、再生核质点法(RKPM)、单位分解法(PUM)、自然单元法(NEM)和 无单元流形法(MM)等，其中由Belytschko等[2] [3] 、Lu Y Y等 对Nayroles的模糊单元法进行了改进而提出的 无单元伽辽金法(Element-free Galerkin Method)的影响最大、应用最广。无单元伽辽金法采用滑动最小 二乘法拟合来近似场函数，采用伽辽金法对偏微分方程求解，并采用背景网格高斯积分聚合整体平衡方程。 由于拟合函数的自由度数远远低于插值函数的基点数，因而无单元法既能提高计算精度，又能保证形函数 导数(即应力场)的连续性。国内对无网格法的研究相对较晚，最早见于周维垣、寇小东[4]引入无单元伽辽 金法并将其应用在岩土工程中，之后还应用无单元法追踪二维裂纹扩展[5]。后来文献[6,7]对无网格法作了 系统的研究，此后国内很多学者在不同领域中应用无网格法以解决二维问题[8。9] 。本文将无单元伽辽金法 推广应用到三维问题，基于变分原理推导三维线弹性无单元伽辽金法的整体平衡方程，并编制计算程序， 通过算例验证方法及程序的有效性。 1 三维无单元伽辽金法 1.1 滑动最小二乘法 设场函数为u(x)，其中x=(x,y,z)T为场点，给定u(x)在n个节点上的值 收稿日期：2004-10-12 基金项目：国家基础科学研究(973)项目(2002-CB412708) 作者简介：周维垣(1928-)，男，河南人，教授，主要从事岩石力学及高坝结构方面研究。 E-mail:zhouwy@tsinghua.edu.cn 1 水 利 学 报 2005 年6 月 SHUILI XUEBAO 第36 卷 第6 期 u(x )=u (i=1,2