5.1 拉格朗日、欧拉与 ALE：网格的三种活法

5.1 拉格朗日、欧拉与 ALE：网格的三种活法

拉格朗日网格跟着材料走（界面精确、大变形时网格扭曲失真），欧拉网格固定在空间（任意流动都不怕、物质界面弥散），ALE 介于两者之间（网格可动、按需重划，兼得精度与稳健，代价是算法复杂与计算成本）。爆炸模拟的第一道选型题，就是决定网格怎么动。

本章开篇：坐标描述方式是所有离散方法的母题。第 1 章的控制方程组是连续的，离散之前先要回答"方程在谁身上写"。

两种极端：贴身跟随与定点观察

拉格朗日描述把网格"焊"在材料上，每个单元一辈子装同一团物质。好处直接：材料界面永远与网格边界重合，自由面、材料边界一清二楚；历史量（应变、损伤）天然跟随质点累计，第 4 章的层裂判据、损伤演化写在单元里就行。坏处同样直接：材料大变形时网格跟着扭曲，单元被拉成扁条后精度崩塌、时间步长被压到不可用——典型的失败场景是爆炸焊接的射流区、爆破抛掷的土石飞散。

欧拉描述把网格固定在空间，材料从中"流过"。单元永远规则，任意大变形都不怕，时间步长稳定——处理爆炸产物气体膨胀、空中冲击波传播是它 的主场。代价是界面：两种材料流过同一单元时只能记录体积分数，界面被抹成几个单元宽的过渡带。对气体动力学这无所谓，对"钢板在爆轰产物推动下变形"这种界面主导的问题，界面模糊直接毁掉精度。

图：三种描述方式对同一压缩过程的处理。中间绿色区域是欧拉法里"材料流过单元"的体积分数表示。

ALE：折中出来的工业主力

ALE（任意拉格朗日-欧拉）让网格以任意速度运动：材料变形不大时网格跟随（保界面），临近扭曲时网格相对材料滑动（保形状）。每个时间步分三拍执行：拉格朗日拍（网格随材料动）、重划拍（按质量准则调整网格）、输运拍（物质在调整后的网格间重新分配）。爆炸驱动的结构响应、水中爆炸气泡这类"固有大变形但界面又重要"的问题，ALE 是工业默认选择。

主流软件里，LS-DYNA 与 AUTODYN 都提供完整三套方案，选型经验成表：

💡 关键直觉：选型的本质是"把网格畸变和界面模糊这两个误差，挪到你更不在乎的那个上"。结构响应问题界面是命根子，选拉格朗日；流场问题界面只是过渡，选欧拉。

⚠️ 常见坑：侵蚀单元（单元变形超标即删除）被当成万能保命符。删掉的单元带着质量与能量一起消失，若侵蚀阈值设置过松，破片数量和飞散速度系统性偏低——勘测复算时必须检查侵蚀单元的质量丢失率，超过百分之几就要存疑。

单元尺寸与时间步长的硬约束

无论选哪种活法，两个数值参数卡住所有模拟的咽喉。单元尺寸：波头要抹在三个单元以上（人工粘性的工作宽度），因此关键波阵面前沿的单元尺寸应小于脉冲空间宽度的十分之一——按第 4 章波速乘脉宽估计空间宽度，再反推网格尺寸，这一步在建模前完成，而不是跑完发现波形不对再返工。时间步长：显式格式的稳定性要求步长小于声波穿越最小单元的时间（Courant 条件），最小单元决定全局步长——一处畸变的小单元能把千万单元的模型拖慢数倍，这解释了拉格朗日网格畸变"杀步长"的机理，也是 ALE 重划拍存在的经济理由。

工程折中的常见做法是分区网格：关心波形的接触区细网格、远场粗网格，中间用过渡带连接。过渡带的网格尺寸渐变率要控制在一到一点二之间，突变过渡等于人为制造反射界面——又一条"数值界面即物理界面"的纪律。

一个选型实例

以水下接触爆炸攻击舰船壳板为例走一遍选型：炸药与水用欧拉（大范围流场、水下气泡脉动），壳板与骨架用拉格朗日（结构响应与塑性历史），两者之间流固耦合传力；壳板破口附近的飞溅碎片若需要追踪，局部转 SPH。单一方法吃不下这个问题，混合建模才是常态——而混合的耦合界面正是误差高发区，对账时耦合参数（罚刚度、耦合点密度）要单独做敏感性检查。

网格质量的勘测判据

复算他人模型时，网格质量可以快速体检三项：一是长宽比，畸形单元的长短轴比超过五比一是预警线；二是最小单元尺寸，它决定步长与总耗时，也暴露建模者是否在关键区加密；三是界面两侧网格尺寸过渡率，超过一点二的跳变等于埋下反射界面。三项指标多数软件都能一键输出，五分钟体检能省五天复算——勘测者的时间应花在物理上，而不是替人收拾网格。

计算成本的量级感

给一组参考量级建立预算直觉：十万单元级的显式模型在单工作站跑数小时，百万单元级要数天，千万单元级进集群。预算与精度之争永远存在，混合方案与分区加密是主要省法，而"为百分之五的精度多花十倍机时是否值得"应回到证据需求回答——勘测复算要的是区间内的结论，不是排行榜上的分辨率。

再补一句关于文件与工况管理：每个模型版本连同输入卡与后处理脚本一起打包归档，命名含日期与工况号——勘测复算追溯到"哪一版模型出的哪一条曲线"靠的就是这份纪律，多数返工的根源是版本混乱而非物理错误。

选型一节的小结话术也值得统一：技术报告里写"采用某方案"要附一句理由（本问题大变形且界面重要），让后来的复算者不必反推你的思路——可复算性从选型说明就开始了。

网格三活法的比喻收束：拉格朗日是贴身保镖、欧拉是岗亭哨兵、ALE是巡巡逻队——各守其位，安防才完整。

选型的最后一句话送给初学者：不确定时先用最简单方案跑通全流程，再换精细方案——粗糙但完整的流程，永远优于精细但断在半路的尝试。

本节要点回顾

• 拉格朗日：界面精确、历史量天然跟随，大变形畸变是死穴；
• 欧拉：任意流动不怕、界面弥散是代价，流场问题主力；
• ALE 三拍：拉格朗日动、重划、输运，工业界的大变形默认解；
• 侵蚀单元：救命与失真一线之隔，质量丢失率是必查项。

下一节看网格彻底用不了的场景：无网格方法如何让每个粒子自带"身份证"地参与计算。