任意拉格朗日2欧拉描述法研究进展.PDF

第 14 卷第 1 期 计 算 力 学 学 报 V o l. 14 N o. 1 1997 年 2 月 CH IN ESE JOU RNAL O F COM PU TA TLONAL M ECHAN ICS February 1997 任意拉格朗日欧拉描述法研究进展 张　雄　　陆明万　　王建军 (清华大学工程力学系, 北京, 100084) 摘　要　任意拉格朗日欧拉( ) 描述综合了纯拉格朗日和纯欧拉描述的优点, 克服了各 AL E 自的缺点, 成为非线性连续介质力学中大变形分析的非常有效的方法。本文论述了AL E 法 的研究进展及其在流体动力学、流体结构相互作用、加工成型、碰撞、接触等大变形问题中的 应用。 关键词　非线性连续介质力学; 任意拉格朗日欧拉描述; 非线性有限元 分类号　 33; 35; 242 21 O O O 1　引　　言 非线性连续介质力学的有限元分析方法已取得了很多令人瞩目的进展, 其中在大变形问 ( ) ( 题中绝大多数研究工作都采用拉格朗日 主要用在固体力学中 或欧拉 主要用在流体力学 ) 中 描述方法。 ( ) 在拉格朗日描述 即物质描述或L 描述 中, 计算网格固定在物体上随物体一起运动, 即 网格点与物质点在物体的变形过程中始终保持重合, 因此物质点与网格点之间不存在相对运 ( ) 动 即迁移运动, 也称对流运动 。这大大地简化了控制方程的求解过程, 而且能准确描述物体 的移动界面, 并可跟踪质点的运动轨迹。但在涉及到特大变形的问题中, 物质的扭曲将导致计 算网格的畸形而使得计算失败。 ( ) 在欧拉描述 即空间描述或E 描述 中, 网格固定在空间中, 即计算网格在物体的变形过 程中保持不变, 因此可很容易处理物质的扭曲。但对运动界面需要引入非常复杂的数学映射, 将可能导致较大的误差。另外当使用普通的伽辽金离散时, 由于迁移项的影响, 有限元方程中 的系数矩阵是非对称的, 而且还可能得到振荡解。 纯拉格朗日和纯欧拉描述都存在严重的缺陷, 但也具有各自的优势。如果能将二者有机地 结合, 充分吸收各自的优势, 克服各自的缺点, 则可解决一大批只用纯拉格朗日和纯欧拉描述 所解决不了的问题。任意拉格朗日欧拉( ) 方法就是基于此目的最早由 ( 1964) [ 1 ] 以耦 AL E N oh 合欧拉拉格朗日的术语提出的, 并用有限差分法求解带有移动边界的二维流体动力学问题。 在N oh 的研究工作中, 网格点可以随物质点一起运动, 但也可以在空间中固定不动, 甚至网格 点可以在一个方向上固定, 而在另一个方向上随物体一起运动, 因此AL E 描述也被称为耦合 欧拉拉格朗日描述。例如在液体表面波的传播问题中, 网格点在垂向随物质点一起运动, 而在 核工业科学基金资助。 本文于 1996 年 5 月 27 日收到。 © 1995-2004 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 92 计 算 力 学 学 报 　　14 卷 水平方向上固定不动。这样可很容易描述液体表面的运动, 而且网格不会发生扭曲。