第 5–6 章 · 等参单元与数值积分

真实结构是曲线的，直边三角形拼不出圆孔。等参思想用一个映射解决：在"标准正方形单元"上定义形函数，再用同一套形函数把正方形映射成任意曲边四边形。这是有限元几何能力的来源。

等参映射

父单元（自然坐标 $(\xi,\eta) \in [-1,1]^2$）到物理单元的映射用同一形函数：

"等参"即"同一个参数化"：几何与位移用同一组 $N_a$。微分要经过雅可比矩阵换算：

$\det\mathbf{J} \le 0$ 意味着单元翻转，一切完蛋——网格质量检查的数学出处就在这。

高斯积分：用很少的点算准积分

单元刚度阵里的积分在父单元上做，被积函数是多项式（线性单元）或有理函数（映射后）。高斯求积用 $n$ 个点精确积分 $2n-1$ 次多项式：

二维四边形单元用 $2\times2$ 高斯点（四个点）处理双线性单元就够了——这四个高斯点后来还成了应力输出与误差估计的"采样站"。

减缩积分与沙漏

全积分的二次单元在弯曲问题时过于刚硬（剪切锁死、体积锁死）。减缩积分（少取点）能解放柔度，却引入零能量模式——沙漏模式：节点位移上下交错但应变处处为零，结构像手风琴一样假变形。补救是沙漏控制（加微小人工刚度）。锁死与沙漏这对矛盾是单元技术几十年的主线。

老师的提醒

• 看到"hourglass energy"警告别忽略：沙漏模式不控制，结果形似合理、本质全错。
• 曲边单元在尖角处（裂纹尖端）雅可比恶化——裂纹问题要用专门的奇异单元，不是加密普通单元。

自检

• 等参映射里"等参"两个字指什么？
• 二点高斯积分能精确积分几次多项式？
• 减缩积分省了什么、换来什么风险？

更深一层：等参映射是有限元的"几何自由"

直线单元拼不出圆孔，曲边单元的诞生让有限元获得了几何自由——等参映射的妙处在于一套形函数两用：既插值位移，又插值几何。这个"同一参数化"的选择不是唯一可能（超参单元几何阶高于位移阶、次参反之），但等参的对称性让推导统一、性质优雅。映射的代价是雅可比：所有对物理坐标的微分都要经过 $\mathbf{J}^{-1}$，被积函数从多项式变成有理函数（分母含 $\det J$）——这就是为什么"精确积分"不再可行、数值积分成为标配。

雅可比的三个工程解读值得背下来：$\det J > 0$（单元不翻转，合法性的底线）、$\det J$ 接近常数（单元形状方正，精度高）、$\partial$ 网格 = $\det J$ 剧烈变化（应力精度差）。网格质量检查报告里的 skewness、aspect ratio、Jacobian ratio，全都是这三条的图形化翻译。

推导细读：高斯积分为什么"以少胜多"

同样 $n$ 个点，等距的牛顿-科茨公式精确到 $n-1$ 次，高斯积分精确到 $2n-1$ 次——两倍效率的来源是节点位置的优化。直觉版：等距点浪费了端点附近的采样（端点权重分布不经济），高斯点恰好是"勒让德多项式 $P_n$ 的零点"，让误差多项式与前 $n$ 个正交多项式全部正交，从而把可精确积分的次数顶到 $2n-1$。数值上二点高斯积分 $\int_{-1}^1\xi^2d\xi$ 用两个点就精确命中 $2/3$（精确到 3 次）。

工程取舍的核心是减缩积分：双线性单元被积函数名义上 4 次，$3\times3$ 才"精确"，但 $2\times2$ 已能通过分片检验且更柔（不易锁死）；单点积分则引入沙漏。积分阶数的最终裁判是分片检验 + 基准算例，不是数学上的"精确"——这个反直觉的结论是第 7 章全部讨论的种子。

在真实世界里：沙漏事故的著名教训

历史上因沙漏模式失控导致的仿真事故不止一例：某车身碰撞分析中薄壁件全部采用单点积分且未开沙漏控制，结果整车吸能曲线"看起来"正常，但后处理动画里车架呈手风琴状锯齿变形——沙漏能占系统应变能 40%，结论全部作废重算。行业由此形成的检查清单沿用至今：单点/减缩积分单元必须启用沙漏控制；报告中必须给出沙漏能与内能之比（警戒线 5%）；动画检查变形形态有无规律性锯齿。沙漏是"便宜的单元"暗藏的账单，收银台在事后验证环节。

常见误区与辨析

• "曲边单元一定更准"：几何拟合好不等于力学精度高；在应力光滑区，直边单元加密比曲边单元稀疏更划算。曲边的真正价值在高梯度区与曲面贴合（圆孔、翼型）。
• "雅可比为负软件会报错"：多数软件只检查积分点处的 $\det J$；畸形单元在积分点之间可能翻面而漏检——关键区域人工检查网格质量仍是必要的。
• "积分点越多越准"：对锁死问题恰恰相反（全积分加剧锁死）；积分阶数是"防锁死 vs 防沙漏"的天平，不是精度旋钮。

专题深入：减缩积分与沙漏的"军备竞赛"

减缩积分与沙漏的军备竞赛是单元技术几十年的主线：减缩积分解放柔度（剪切锁死、体积锁死的解药）但引入零能量模式（沙漏）；沙漏控制（人工刚度/粘性）压制沙漏但损害精度。现代求解器的标准配置：减缩积分 + 沙漏控制 + 沙漏能监控（占比 < 5-10% 可信）——三件套的配套使用。沙漏模式的图像：节点位移上下交错但应变处处为零——结构像手风琴一样假变形，形似合理、本质全错。这就是为什么"动画检查变形形态有无规律性锯齿"是显式仿真的纪律。

数学细节展开：B-bar 与混合格式——不可压问题的出路

体积锁死（不可压材料 $\nu \to 0.5$）的解法值得展开：B-bar（应变均值法——把体积应变的加权平均替换逐点值——"在约束最紧的维度上放松插值空间"）与混合格式（位移-压力双场插值——独立的静水压力场）。两者的理论根源是Ladyzhenskaya-Babuška-Brezzi（LBB）条件（inf-sup 条件）——位移与压力的插值空间必须满足"不相容的匹配"（否则棋盘压力——与 CFD 的交插网格同源）。Taylor-Hood 单元（二次位移 + 线性压力）是满足 LBB 的经典组合——流体与不可压固体的通用选择。

工程案例细节：橡胶与饱和土的"不可压陷阱"

不可压问题的工程陷阱实录：橡胶密封件（$\nu \approx 0.499$——普通位移单元锁死到几乎不变形，混合格式立刻恢复）、饱和土固结（Biot 耦合 + 不可压孔隙水——双重约束）、金属塑性（J2 塑性在塑性区不可压——大变形塑性必须用 B-bar 或混合）。"泊松比接近 0.5 就换配方"是选型纪律——不是加密能解决的问题，是离散结构性缺陷。