第 7 章 · 数值积分：高斯点数的选择

这一节把上一章埋的伏笔展开：积分阶数不是小事，它同时影响精度、稳定性与计算量，是单元工程的核心决策。

精度账：需要几个点？

单元刚度阵 $\int B^T D B\,d\Omega$ 中，线性四边形的 $B$ 含 $\xi,\eta$ 一次项，被积函数最高四次（两个方向各二次）。精确积分需每方向三点。但实践常取 $2\times2$——原因在于：精确积分的减缩版本已能通过分片检验，且更柔、更好。这说明"积分阶"标准不是数学精确性，而是"能否通过分片检验 + 单元性能"。

锁死：全积分的代价

• 剪切锁死：薄梁/薄板弯曲时，线性单元的全积分点处寄生剪切应变远大于真实值，单元刚硬到几乎不弯。八节点二次单元全积分仍轻度锁死，减缩积分后成为最优秀实体单元之一。
• 体积锁死：不可压材料（泊松比 → 0.5，橡胶、饱和土）时，常压力约束在积分点上过约束，单元刚死。解法是 B-bar（应变均值法）或混合格式。

零能模式：减缩的代价

$1\times1$ 单点积分的双线性单元有多个零能量模式（沙漏）。诊断方法：看应变能中沙漏能占比，超过 5-10% 即不可信。控制手段：人工沙漏刚度（小而准）或物理稳定化。

实用决策表

• 三维一般问题：八节点六面体 + $2\times2\times2$ 积分（入门首选）或 20 节点全积分（精度基准）。
• 薄壁弯曲：二次单元 + 减缩积分，或专门的壳单元。
• 不可压：混合格式或 B-bar，不用普通位移格式硬算。

老师的提醒

• 单元测试三板斧：常应变（分片检验）、悬臂梁（锁死）、单轴拉伸不可压（体积锁死）。新单元先过这三关再上工程。
• 高阶积分不等于高精度：锁死路径上多花的每个积分点都在帮倒忙。

自检

• 为什么 $2\times2$ 积分对双线性单元"不精确却合法"？
• 剪切锁死与体积锁死的机理分别是什么？
• 沙漏能占比多大时应怀疑结果？

更深一层：积分阶数是"单元设计"的最后一块拼图

把第 4 节的等参映射与本节的积分阶数放在一起看，一个单元的完整"配方"就齐了：插值阶数（形函数）+ 几何映射（等参）+ 积分阶数（高斯点数）。三者构成一个设计空间，工业单元是这个空间里经过几十年基准赛淘洗出的几个甜点。最著名的甜点是"八节点六面体 + $2\times2\times2$ 减缩积分"（C3D8R 的配方）：全积分会让它剪切锁死，单点积分会沙漏，$8$ 点减缩积分恰好落在"通过分片检验、足够柔、沙漏可控"的三角平衡点上。这个配方撑起了全球碰撞与冲压仿真的主体算力——单元设计的成就最终以"谁在支撑工业算例"来计量。

积分阶数还有一面少有人提的作用：它是误差的采样站。应力在积分点上超收敛（比其他位置更接近真解），所以后验误差估计、应力恢复、自适应加密的判据全部读积分点数据。换句话说，高斯点既是算积分的账房，又是误差估计的哨兵——理解这一点，第 7 节的 ZZ 恢复法就毫无神秘可言。

推导细读：锁死的能量语言

锁死（locking）听起来玄，用能量语言一句话说穿：单元的伪应变能吃掉了真实应变能。薄梁纯弯曲时，真实应变能几乎全是弯曲能；双线性单元全积分点上除了弯曲还被迫产生寄生剪切应变（线性位移场无法表达"弯曲时截面保持平面且垂直"的纯弯变形），寄生剪切能 $\propto$ 剪切模量 $\times$ 剪应变的平方——单元为了存这点伪能量，把本该变成弯曲变形的功全吞了，表现为"怎么压都压不动"。数值诊断：把网格加密，伪剪切不收敛（它是离散产物），挠度被系统性低估且不随加密改善——这是锁死与普通离散误差的决定性区别（后者加密必收敛）。

体积锁死同理：不可压材料（$\nu \to 0.5$）的静水约束在每个积分点上过约束位移场，单元"胀不动也压不动"。解药 B-bar 方法的思路漂亮：把体积应变的加权平均替换逐点值——在约束最紧的维度上放松插值空间。这个思想后来长成混合格式（位移-压力插值）的大家族。

在真实世界里：单元选型的一条快速决策链

工程现场遇到"该用什么单元"，可以走这条十秒决策链： 1. 几何厚度远小于其他尺寸（< 1/10）？→ 壳或梁，别用实体硬算薄壁（锁死 + 网格浪费）； 2. 材料接近不可压（橡胶、饱和土）？→ 混合格式 / B-bar / C3D10H，普通位移元直接出局； 3. 弯曲主导 + 精度要求高？→ 二次单元 + 减缩积分，或厚板理论壳； 4. 接触 + 大变形 + 显式动力学？→ 减缩积分实体 + 沙漏控制（显式显式界的工业标配）； 5. 只算刚度不算应力细节？→ 线性单元粗网格足够。

把这条链贴在工位上，能省掉 80% 的"单元该怎么选"讨论。

常见误区与辨析

• "减缩积分会不准"：通过分片检验的减缩积分在光滑问题上精度不降反升（更柔、更接近真解）；不准的是沙漏失控时的解。
• "锁死加密网格就好"：锁死是离散结构性缺陷，加密只是花钱买同样的锁死；换配方（减缩/混合/高阶）才是正解。
• "沙漏能小就没事"：沙漏能占比小仍可能有局部沙漏变形污染应力细节；看占比之外还要看变形动画。

专题深入：高斯积分的"勒让德零点"来历

高斯积分的"两倍效率"来历值得从正交多项式角度说透：高斯点是勒让德多项式 $P_n$ 的零点，权重由插值公式唯一确定。误差的分析：用 $n$ 个点精确积分 $2n-1$ 次多项式——这个"双倍"来自节点位置的优化（等距点的牛顿-科茨只能到 $n-1$ 次）。工程上的"甜点"清单：双线性单元用 $2\times2$、双二次用 $3\times3$、四面体与三角形有专门的对称点集。"积分阶数 = 精度与锁死的权衡"——减缩积分的合法性由分片检验担保，锁死风险由基准算例暴露。

数学细节展开：剪切锁死的"寄生应变"解剖

剪切锁死的解剖值得做一次：薄梁纯弯曲时，双线性单元的位移场是线性的——但纯弯曲要求位移含二次项（挠度抛物线），线性单元被迫用"阶梯状线性逼近抛物线"，寄生剪切应变在积分点上产生（不真实的剪切），单元"假刚"到几乎不弯。加密网格不收敛（锁死是离散结构性缺陷）、减缩积分立刻恢复——"锁死 vs 离散误差"的决定性区别：加密是否收敛。这个判据在工程诊断里是分水岭：加密收敛 = 普通离散误差、加密不收敛 = 锁死。

工程案例细节：悬臂梁基准的"单元体检"

悬臂梁基准（端部集中力、理论挠度已知）是单元体检的标准：双线性全积分 → 挠度只有理论的 38%（剪切锁死现场）；双线性减缩积分 → 93%（解锁但轻微沙漏）；二次单元全积分 → 99.5%（二次位移场表达弯曲）；二次减缩积分 → 99.9%（最优秀实体单元之一）。这个四行表格就是单元选型的"成绩单"——每个新单元先在悬臂梁上跑一遍，成绩立刻可读。

常见误区与辨析（补充）

• "沙漏能小就没事"：占比小仍可能有局部沙漏污染应力细节——看占比之外还要看变形动画。
• "积分点越多越准"：对锁死问题恰恰相反——全积分加剧锁死；积分阶数是"防锁死 vs 防沙漏"的天平。