计算固体力学7_ALE公式.ppt

作为一个使用修正的弹性方程的ALE网格更新的例子，考虑一个位于一个矩形流体域中的圆柱绕流的有限元网格。在圆柱沿着y方向移动了一个位移0.25w后的网格见图b。矩形域的边界保持固定。由图可见，圆柱附近的网格精度在更新后的网格中得以保持，并且无明显的单元扭曲发生。 网格更新的例子 6 网格更新算法 7 Petrov-Galerkin方法 伽辽金(Galerkin)方法，是利用满足位移和应力边界条件 的函数寻求积分方程的解答。 关键问题是增加粘性，消除不稳定项，保证数值稳定性。 建立Petrov-Galerkin方法的迎风流线(Streamline Upwind Petrov-Galerkin－SUPG)公式。对流－扩散方程是一个有用的方法，它对应于动量方程的线性化。对于离散的稳态对流—扩散方程，将得到闭合解答。将证明当网格参数(已知的Peclet数)超过临界值时，这个解答在空间是振荡的。通过建立P-G方法以消除这些振荡，即纠正不稳定性。 在一维中，离散方程类似于标准迎风方程。然而在多维中，它们提供了沿着流线引导迎风项的一致理论框架。 迎风格式的基本思想：当用差分方程求解偏微分方程时，利用特征线方向一侧的单边差商来代替空间偏导数。目的是保证数值稳定性。如一阶线性常系数双曲型方程： 7 Petrov-Galerkin方法 在特征线方向一侧的单边差商来代替偏导数，迎风格式为： j-1 n+1 a0 n j j+1 j-1 n+1 a0 n j j+1 中心差分格式需要加入非线性人工粘性项以去掉激波前振荡，同时加入线性数值粘性项以确保计算稳定。 迎风格式依据Euler方程中波传播的信息构造格式，无需添加人工粘性项，迎风格式为空间差分格式(偏心差分格式)，理论上比中心差分格式扎实。 7 Petrov-Galerkin方法 如果差分格式(所用的网格点)与微分方程的特征线方向一致，那么网格比 在满足一定条件下是稳定的，否则，差分格式是不稳定的。 双曲线型微分方程的特征线－存在交叉，导数不连续。 7 Petrov-Galerkin方法 稳态线性对流－扩散方程为 运动粘度 给定速度 对于一维问题，偏微分方程成为常微分方程 应用边界条件 这是在0xL 域上的两点边值问题，容易证明公式的精确解答是 空间非独立变量 7 Petrov-Galerkin方法 应用线性形状函数建立Galerkin离散，并在全域上积分 变分函数 分部积分并且应用散度原理，对流－扩散方程的弱形式是 将域（0，L）划分成相同尺寸的单元 在每个单元上的离散方程给出为 有限元的形状函数 对于第j个节点的内部方程为 7 Petrov-Galerkin方法 上式恰好是中心差分方程，可以方便地重写成为 1 如果Peclet数小于1， 则离散解答是类似于精确解答； 2 如果Peclet数大于1， 则离散解答是正或者负而振荡。 这种不稳定是数值离散的空间不稳定。 Peclet数 精确解 离散解 * * 非线性有限元第7章 任意的Lagrangian和Eulerian公式 计算固体力学 第7章 任意的Lagrangian和Eulerian公式 引言 ALE连续介质力学 ALE守恒规则 ALE控制方程 弱形式 网格更新算法 Petrov-Galerkin方法 1 引言 解决：在发生严重大变形的模拟中，重新划分网格是不可避免的，工作量大，而且由于网格投影引入了误差。 提出：许多问题应用Lagrangian网格不能有效地解决。 问题：当材料严重变形时，Lagrangian单元同样发生严重的扭曲，因为它们随材料一起变形，从而恶化了这些单元的近似精度，特别是对于高阶单元。因此，在积分点的Jacobian行列式可能成为负值，从而使计算中止或者引起严重的局部误差。此外，也恶化了线性化牛顿方程的条件，并且显式稳定时间步长明显地下降。 一个Lagrangian网格像在材料上的蚀刻：当材料变形时，蚀刻(和单元)随着变形。 一个Eulerian网格像放在材料前面一薄片玻璃上的蚀刻：当材料变形时，蚀刻不变形，而材料横穿过网格。 1 引言 Lagrangian网格，材料点与网格点保持重合，单元随材料变形，适合描述固体与结构的变形，但容易严重扭曲。 解决方法：ALE网格(Arbitrary Lagrangian Eulerian) 节点能够有序地任意运动，在边界上的节点保持在边界上运动，内部的节点运动使网格扭曲最小化。 1 引言 1 引言 Mesh adaptivity is based on solution variables as well as minimum elem