任意拉格朗日-欧拉描述法研究进展

理形状并准确描述物体的移动界面∀所以网格运动算法在 ( )∞ 描述中占很重要的地位∀
- ( )∞ 描述的控制方程
为了建立在 ( )∞ 描述下的控制方程"利用非线性连续介质力学的方法以及高斯定理可 以推出物理量
的物质导数为Y
Θ Θ Γ !τ' Z γ !Ν"τ'_ςΝ 和 Η !τ' Z η!ξ "τ'_ςξ
物体一起运动的∀在空间描述方法中的有限单元剖分是对空间进行的"网格点就是空间点"因
此网格是固定在空间中不随物体运动的∀而在 ( )∞ 描述方法中的有限单元剖分是对参考构形
进行的"网格点就是参考点"网格是独立于物体和空间运动的"可以根据需要自由选择∀
以上映射应该是一对一映射"雅可比行列式 ϑ!描述从初始构形到现时构形的映射关系' 和混合雅可比行列式 ϑδ!描述从参考构形到现时构形的映射关系' 都应不等于 $"即Y
式中Y
Θ Θ Θ Μ Z
Θδ_ςΝ Z
Θ_ςξ Z
Θ$_ς Ξ
8Ν
8ξ
8Ξ
Θδ!Ν"τ' Z ϑδΘ!ξ "τ'" Θ$!Ξ "τ' Z ϑΘ!ξ "τ'
!#/'
质量守恒定律表明"质量的整体变化率!即物质导数'应等于零∀利用式!#/'和式!##1J'可
以得到在参考坐标系下的质量守恒方程为Y
55Θτδ" Ν ]
网格运动规律可以将 ( )∞ 描述退化为 Λ 描述和 Ε 描述∀ # ϖδ Z $"即计算网格在空间中固定不动"退化为欧拉描述∀ # ϖδ Z ϖ"即计算网格随同物体一起运动"退化为拉格朗日描述∀ # ϖδ Ξ ϖΞ $"即计算网格在空间中独立运动"对应于一般的 ( )∞ 描述∀
在 ( )∞ 描述中参考构形是已知的"各物理量用参考坐标 Ν描述比较方便"即 Φ Z Φ!Ν"O'∀
中的位置矢量 Ν确定"因此
ΝZ Ν!Ξ "τ'
描述了物质点 Ξ 在参考坐标系中的运动规律"
!#1χ'
#期
张 雄等Y任意拉格朗日2欧拉描述法研究进展
3-
ξ Z ξ !Ν"τ'
!#1δ '
则描述参考点 Ν在空间中的运动规律∀
在物质描述方法中的有限单元剖分是对物体进行剖分的"网格点就是物质点"即网格是随
Ξ 核工业科学基金资助∀ 本文于 #334 年 . 月 /: 日收到∀
但对运动界面需要引入非常复杂的数学映射"
3/
τ
计 算 力 而学 学 报
#& 卷
水平方向上固定不动∀ 这样可很容易描述液体表面的运动"而且网格不会发生扭曲∀ 在 ( )∞ 描述中计算网格可以在空间中以任意的形式运动"即可以独立于物质坐标系和
!:'
式中Yχι Z ϖι ^ ϖδι 为物质点相对于网格点的运动速度"即迁移速度!也称对流速度'∀从式!:' 可
以看出"迁移速度 χ 实质上是将质点在参考坐标系中的运动速度 ω 变换到现时坐标系中而得
到的"也就是说 χ 是在现时坐标系中物质点相对于参考构形的运动速度∀式!:' 也可化为Y
χι Z 5ξ ι!Ν!Ξ5τ"τ'"τ'" Ξ
格点'Ν在空间中的位置矢量 ξ !Ν"τ' 对时间的导数"即
ϖδ Z 5ξ !5Ντ"τ'" Ν
!&'
物质点 Ξ 在参考坐标系中的位置矢量 Ν!Ξ "τ' 对时间的导数
ω Z 5Ν!Ξ5τ"τ'" Ξ
!.'
则为物质点 Ξ 在参考坐标系中的速度∀
在 ( )∞ 描述中"参考构形!也就是计算网格'的运动规律可以是任意给定的"指定特殊的
在欧拉描述!即空间描述或 Ε 描述' 中"网格固定在空间中"即计算网格在物体的变形过 程中保持不变"因此可很容易处理物质的扭曲∀ 将可能导致较大的误差∀另外当使用普通的伽辽金离散时"由于迁移项的影响"有限元方程中 的系数矩阵是非对称的"而且还可能得到振荡解∀
纯拉格朗日和纯欧拉描述都存在严重的缺陷"但也具有各自的优势∀如果能将二者有机地 结合"充分吸收各自的优势"克服各自的缺点"则可解决一大批只用纯拉格朗日和纯欧拉描述 所解决不了的问题∀任意拉格朗日2欧拉!( )∞'方法就是基于此目的最早由 012!#34&'≈#6以耦 合欧拉2拉格朗日的术语提出的"并用有限差分法求解带有移动边界的二维流体动力学问题∀ 在 012 的研究工作中"网格点可以随物质点一起运动"但也可以在空间中固定不动"甚至网格 点可以在一个方向上固定"而在另一个方向上随物体一起运动"因此 ( )∞ 描述也被称为耦合 欧拉2拉格朗日描述∀例如在液体表面波的传播问题中"网格点在垂向随物质点一起运动"而在
8Ν
8ξ
!#$'
Θ 5Γ
5τ
"
Ξ
Z
பைடு நூலகம்
8 Ν!55γτ" Ν ]
5ω ιγ 5Νι
'_ς Ν
!##1J'
Θ 5Η
5τ
"
Ξ
Z
8 ξ !55ητ" ξ ]
5ϖιη 5ξ ι
'_ς
ξ
!##1G'
311 质量守恒方程!连续性方程'
考察一连续体"用 58 Ξ !58 ξ 和 58 Ν 分别表示物质域 8 Ξ !空间域 8 ξ 和参考域 8 Ν的边界"而用 Θ$!Θ和 Θδ 分别表示连续体各构形的密度∀在不同构形中该连续体的质量 Μ 可以写为Y
ϑZ
"
5ξ 5Ξ
ι"
ϕ
Ξ
$"
ϑδ Z " 55ξΝϕι" Ξ $
!/'
根据定义 τ时刻某质点 Ξ 在空间中的运动速度 ϖ等于质点 Ξ 在空间的位置矢量 ξ !Ξ "τ' 对
时间的导数"即Y
ϖZ 5ξ !5Ξτ"τ'" Ξ
!-'
参考构形中某点 Ν在空间中的运动速度!也就是网格点的运动速度'ϖδ则等于参考点!即网
本文论述了 ( )∞ 法的基本理论和研究进展及其在流体动力学!流体2结构相互作用!加工 成型!碰撞!接触等大变形问题中的应用∀
/ ( )∞ 运动学描述
将连续体在初始时刻 $ 的构形!称初始构形' 记为 8 Ξ "将 τ时刻的构形!称现时构形' 记为 8 ξ ∀为了确定各物质点的位置"引入 )JSHJMSF 或物质坐标系 ΟΞ #Ξ /Ξ -!简记为 Λ 坐标系'"它 是和物体固结在一起并随物体一起运动的∀质点在初始时刻 το 的位置由其在 Λ 坐标系中的位 置矢量 Ξ 完全确定"因此 Ξ 也可用来识别物体中的不同质点∀为了描述现时构形"引入空间坐 标系或 Ε υλερ 坐标系 οξ #ξ /ξ -!简记为 Ε 坐标系'"它是和空间固结在一起的∀空间中各点的位置 由其在 Ε 坐标系中的位置矢量 ξ 确定"因此可用 ξ 表示空间中的几何点∀
因此各物理量的物质导数!即某固定质点 ÷ 的物理量对时间的变化率' 应通过该物理量的参
考导数!即某固定参考点 Ν处的物理量对时间的变化率"也称混合导数' 来计算∀为了表达方
便"本文采用指标求和法则∀对 ¬ 求物质导数"得Y
ϖι Z
ϖδι ]
5ξ
ι!Ν"τ' 5Νϕ
#
ω
ϕ
!4'
将式!4'整理得Y
χι Z 5ξ ι!5ΝΝϕ"τ'ω ϕ
#引 言
非线性连续介质力学的有限元分析方法已取得了很多令人瞩目的进展"其中在大变形问 题中绝大多数研究工作都采用拉格朗日!主要用在固体力学中' 或欧拉!主要用在流体力学 中' 描述方法∀
在拉格朗日描述!即物质描述或 Λ 描述' 中"计算网格固定在物体上随物体一起运动"即 网格点与物质点在物体的变形过程中始终保持重合"因此物质点与网格点之间不存在相对运 动!即迁移运动"也称对流运动'∀这大大地简化了控制方程的求解过程"而且能准确描述物体 的移动界面"并可跟踪质点的运动轨迹∀但在涉及到特大变形的问题中"物质的扭曲将导致计 算网格的畸形而使得计算失败∀
第 #& 卷第 # 期 #33: 年 / 月
计算力学学报
≤ < =0∞≥∞ ?, @ A 0( ) , ƒ ≤ , C °@ × ( × ), 0( ) C ∞≤ < ( 0=≤ ≥
∂ 189#& 019# ƒ FGH∏JHK #33:
任意拉格朗日2欧拉描述法研究进展Ξ
张 雄 陆明万 王建军
!清华大学工程力学系"北京"#$$$%&'
利用散度定理以及式!##1J'!式!#-1J' 可将上式写成Y
Θδ 55ϖτι" Ν ]
Θδω
ϕ
5ϖι 5Νϕ
Z
5Τδϕι 5Νϕ
]
Θδφ ι
在 8 Ν 域中
!#.1J'
#期
张 雄等Y任意拉格朗日2欧拉描述法研究进展
3.
式中 Τδϕι 是定义在参考构形下的第一类皮奥拉 2 克希荷夫应力张量!即拉格朗日应力张量'∀它 与真实应力!即柯西应力' 张量 Ριϕ 之间的关系为Y
!%'
3&
计算力学学报
#& 卷
利用式!:' 可将任意物理量 Φ!Ν"τ' 的物质导数表示为Y
5Φ 5τ
"
Ξ
Z
5Φ!5Ντ"τ'" Ν ]
χι
5Φ 5ξ ι
!3'
式!3'是 ( )∞ 描述中很重要的关系式"利用它可将控制方程转化到参考坐标系下求解∀ ( )∞ 描述的重要特征是可以根据需要给定合适的网格运动速度 ϖδι"以维持计算网格的合
( )∞ 描述中"参考构形是已知的"而初始构形和现时构形都是待求解的∀ 因此 ( )∞ 法尤 其适合于在初始构形和现时构形都未知的问题中使用"如接触问题等∀
目前"( )∞ 法已被应用于固体力学领域中求解大变形问题"如碰撞≈-46≈..6! 接触和弹性断 裂力学≈-#∗ -&6!路径相关材料!如弹塑性材料等'≈-.∗ -36!加工成型≈&$∗ &:6等∀ 另外对网格的运动以 及控制方程的求解策略等问题也进行了多方面的研究≈.-∗ .%6∀
拉格朗日描述以初始构形为参考构形研究物质点 Ξ 在空间中的运动规律"即Y
ξ Z ξ !Ξ "τ'
!#1α'
而欧拉描述则是以现时构形为参考构形来研究空间点 ξ 上物质点的运动规律"即研究Y
Ξ Z Ξ !ξ "τ'
!#1β'
式!#1α'描述了同一质点 Ξ 在不同时刻的空间位置"而式!#1β'则描述了同一空间点 ξ 在
摘 要 任意拉格朗日2欧拉!( )∞'描述综合了纯拉格朗日和纯欧拉描述的优点"克服了各 自的缺点"成为非线性连续介质力学中大变形分析的非常有效的方法∀ 本文论述了 ( )∞ 法 的研究进 展及其在流 体动力学! 流体2结构 相互作用! 加工成型! 碰撞! 接触 等大变形 问 题 中 的 应用∀
关键词 非线性连续介质力学+任意拉格朗日2欧拉描述+非线性有限元 分类号 , --+, -.+, /&/1/#
5Θδω 5Νι
ι
Z
$
在 8 Ν 域中
!#-1J'
类似于修正拉格朗日!@ )'方法"有时在空间域 8 ¬ 中求解控制方程是比较方便的∀利用式
!#/'! 式!##1 G'和式!3'可以把质量守恒方程写成Y
55Θτ" Ν ]
χι
5Θ 5ξ ι
]
Θ
5ϖι 5ξ ι
Z
$
在 8 ξ 域中
!#-1G'
312 动量守恒方程!平衡方程'
( )∞ 最早是为了解决流体动力学问题而引入的"并且使用有限差分法≈#2%6∀ 由于核反应堆 结构安全分析的需要"⁄1MFJ≈#/6!NF8KOPQ2R1≈#-"#&6等人分别将 ( )∞ 法引入有限元法中"用以求 解流体与结构相互作用问题≈#/2-$6∀< ∏S2FP等人≈36建立了 ( )∞ 描述的运动学理论"并使用有限 元法解决了粘性不可压缩流体流动和自由表面流动问题∀
Τδιϕ Z
ϑδ
5Νι 5ξ κ
Ρκϕ
!#4'
如果在参考域中直接求解式!#.1J'"Τδϕι 的不对称性将会给求解带来一定的困难∀由于定
动量守恒定律表明"在 τ 时刻占据参考域 8 Ν 的物体的总动量的整体变化率等于施加在物
体上的外力之和"即Y
Θ Θ Θ 5
5τ
"
Ξ
Θδϖι_ςΝ Z
8Ν
τδι_Σ Ν ]
58 Ν
Θδφ ι_ςΝ
8 Ν
!#&'
式中 τδι 是作用在参考域 8 Ν 的边界 58 Ν 的单位表面上的力"φ ι 是作用于物体中单位质量的体力∀
不同时刻被物质点占据的情况∀
不同于拉格朗日和欧拉描述"( )∞ 描述另外引入了一个可以独立于初始构形和现时构形
运动的参考构形"记为 8 Ν∀在物体的变形过程中"观察者始终跟随参考构形运动"因而对观察
者而言参考构形是固定不动的"而初始构形和现时构形则都相对于参考构形运动∀为了确定参
考构形中各参考点的位置"引入参考坐标系 οΝ#Ν/Ν-"参考构形中各点的位置由其在参考坐标系
空间坐标系运动∀这样通过规定合适的网格运动形式可以准确地描述物体的移动界面"并维持 单元的合理形状∀ 类似于 Ε 描述"在 ( )∞ 描述下的控制方程中也将出现对流项"因此也可能 得到振荡解"需要进行相应的数值处理∀ 纯拉格朗日和纯欧拉描述实际上是 ( )∞ 描述的两个 特例"即当网格的运动速度等于物体的运动速度时就退化为拉格朗日描述"而当网格固定于空 间不动时就退化为欧拉描述∀