厚壳理论（退化实体）
理论与物理
退化实体壳是什么
老师，“退化实体”是什么？是壳单元和实体单元的混合体吗？
正是如此。退化实体壳单元是通过退化（缩减）三维实体单元在板厚方向上的自由度而得到的壳单元。由Ahmad-Irons-Zienkiewicz（1970）提出。
思路很简单：
1. 拥有三维实体单元的顶面和底面节点
2. 将中面的位移和法线方向的旋转角作为自由度
3. 假设板厚方向的位移分布为线性（Mindlin假设）
从三维出发，加入假设变成二维。这与通常的壳理论是相反的思路呢。
是的。通常的壳理论从二维方程出发，而退化实体壳从三维出发，“退化掉不用的自由度”。结果虽然都归结为相同的Mindlin壳，但实现基于三维，因此更简单。
对厚壳的应对
“厚壳”是指什么样的情况？
指 $R/t$ 大约在10~30左右的中等厚度的壳。既不是薄壳（$R/t > 30$），也不是厚壳（$R/t < 10$，实质上是实体）的中间领域。
在厚壳中：
- 剪切变形不可忽略
- 板厚方向的应力 $\sigma_z$ 不完全为零
- 膜-弯曲耦合较强
能用Mindlin壳单元应对吗？
可以应对剪切变形，但 $\sigma_z \neq 0$ 是通常的壳单元无法处理的。要处理这个，需要使用实体壳单元（在板厚方向也具有位移自由度的壳单元）或实体单元。
实体壳单元
“实体壳”是什么样的单元？
外观是实体单元（HEX8或HEX20），但内部公式针对壳进行了优化。
| 单元 | 求解器 | 特征 |
|---|---|---|
| SC8R | Abaqus | 8节点实体壳。减缩积分+锁定对策 |
| SOLSH190 | Ansys | 实体壳。板厚方向1个单元即可表现弯曲 |
| CHEXA(solid-shell) | LS-DYNA | LSDYNA的实体壳实现 |
实体壳的优点是什么？
总结
我来整理一下厚壳理论。
要点：
- 退化实体 — 通过退化三维实体在板厚方向来创建壳
- $R/t = 10 \sim 30$ 的中间领域 — 既非薄壳也非实体
- 实体壳单元 — 外观是实体，内核是壳。应对接触面和$\sigma_z$
- SC8R（Abaqus）, SOLSH190（Ansys） — 代表性单元
- 板厚方向1个单元即可表现弯曲 — 高效
薄壳→Mindlin壳，中间→实体壳，厚壳→实体单元，是这样区分使用的吧。
基本判断依据是 $R/t$。如果拿不准，可以两种方法都计算并比较结果。
Mindlin-Reissner厚壳理论
作为厚壳基础的Mindlin-Reissner理论由Raymond Mindlin和Eric Reissner在1945年至1951年间独立完成公式化。与Kirchhoff假设不同，它允许“法线因剪切变形而倾斜”，并将横向剪切应变εxz和εyz作为自由度显式处理。这使得其适用范围扩展到板厚/跨度比约1/5的厚板，成为复合材料层合板层间剪切分析不可或缺的理论。
各项的物理意义
- 惯性项（质量项）：$\rho \ddot{u}$，即“质量×加速度”。您有过急刹车时身体被向前甩出的经历吗？那种“被带走的感觉”正是惯性力。物体越重越难启动，一旦启动也越难停止。地震时建筑物摇晃，也是因为地面突然移动而建筑物的质量“被落下”。静力分析中此项设为零，那是基于“缓慢施力，加速度可忽略”的假设。对于冲击载荷或振动问题，此项绝对不能省略。
- 刚度项（弹性恢复力）：$Ku$ 或 $\nabla \cdot \sigma$。拉弹簧时会感觉到“想要恢复的力”吧？那就是胡克定律 $F=kx$，也是刚度项的本质。那么提问——铁棒和橡皮筋，用相同的力拉，哪个伸得更长？当然是橡皮筋。这种“不易伸长性”就是杨氏模量 $E$，它决定了刚度。常见的误解是：“刚度高=强度高”。不对。刚度是“不易变形性”，强度是“不易破坏性”，是不同的概念。
- 外力项（载荷项）：体积力 $f_b$（如重力）和表面力 $f_s$（压力、接触力等）。可以这样想——桥上卡车的重量是“作用在整个内部上的力”（体积力），轮胎压路面的力是“只作用在表面上的力”（表面力）。风压、水压、螺栓紧固力…这些都是外力。这里容易犯的错误是：弄错载荷方向。本想“拉伸”却成了“压缩”——听起来像笑话，但在3D空间中坐标系发生旋转时，确实会发生。
- 阻尼项：瑞利阻尼 $C\dot{u} = (\alpha M + \beta K)\dot{u}$。弹一下吉他的弦试试。声音会一直持续吗？不，会逐渐变小。这是因为振动能量通过空气阻力或弦的内部摩擦变成了热。汽车的减震器也是同样原理——故意吸收振动能量来改善乘坐舒适性。如果阻尼为零会怎样？建筑物在地震后会一直摇晃不停。实际上不会这样，所以设置适当的阻尼很重要。
假设条件与适用范围
量纲分析与单位制
| 变量 | SI单位 | 注意事项·换算备忘 |
|---|---|---|
| 位移 $u$ | m（米） | 输入为mm时，载荷·弹性模量也需统一为MPa/N系 |
| 应力 $\sigma$ | Pa（帕斯卡）= N/m² | MPa = 10⁶ Pa。与屈服应力比较时注意单位制不一致 |
| 应变 $\varepsilon$ | 无量纲（m/m） | 注意工程应变与对数应变的区别（大变形时） |
| 弹性模量 $E$ | Pa | 钢：约210 GPa，铝：约70 GPa。注意温度依赖性 |
| 密度 $\rho$ | kg/m³ | mm系中为tonne/mm³（钢为 = 10⁻⁹ tonne/mm³） |
| 力 $F$ | N（牛顿） | mm系用N，m系也用N统一 |
数值解法与实现
实体壳的实现
请告诉我实体壳在实现上的注意事项。
实体壳是“薄”的实体单元，因此板厚方向的长宽比会变得非常大。通常的实体单元在长宽比 > 5 时精度会下降，但实体壳对此进行了内部修正。
锁定对策
实体壳中发生的锁定：
1. 剪切锁定 — 薄板弯曲时发生。用ANS法对策
2. 体积锁定 — 不可压缩材料时发生。用EAS法或B-bar法
3. 梯形锁定（trapezoidal locking） — 板厚方向单元呈锥形时发生。实体壳特有
4. 曲率厚度锁定 — 曲面中板厚方向单元呈梯形时发生
梯形锁定是实体壳特有的吗？
是的。当板厚方向存在锥度（梯形形状）时，通常的实体单元无法正确表现弯曲。实体壳单元通过EAS（Enhanced Assumed Strain）排除了梯形锁定。
使用要点
使用实体壳单元时的注意事项：
- 板厚方向1个单元 — 不需要2个以上单元（实体壳的设计思想）
- 正确指定单元的“厚度方向” — 在Abaqus的 *SOLID SECTION 中指定 stack direction
- 曲面的网格 — 分别对CAD的顶面和底面进行网格划分，再在板厚方向连接
板厚方向1个单元就够了，这很高效啊。通常的实体HEX8需要在板厚方向有4个以上单元呢。
实体壳的最大优点正在于此。用相当于板厚方向1个单元的HEX8，实现了与通常壳单元同等的弯曲精度。自由度数与壳单元相当，但在接触和板厚变化方面有优势。
总结
我来整理一下实体壳的数值方法。
要点：
- 板厚方向1个单元即可表现弯曲 — 高效
- 锁定对策必不可少 — ANS+EAS的组合
- 梯形锁定是实体壳特有 — 用EAS法对策
- 正确指定stack direction很重要 — 正确定义板厚方向
- 最适合接触面位于上下两面的问题 — 通常壳单元不具备的优点
MITC单元的剪切锁定对策
MITC（Mixed Interpolation of Tensorial Components）法是Bathe和Dvorkin于1986年在MIT开发的厚壳锁定对策手法。通过对剪切应变进行独立插值，确保了从薄板到厚板的均匀精度。MITC4对应4节点壳，MITC9对应9节点壳，即使板厚/跨度比为1/1000，也能将位移误差控制在5%以内，证明了其高性能。
线性单元（1次单元）
节点间进行线性插值。计算成本低，但应力精度低。注意剪切锁定（可通过减缩积分或B-bar法缓解）。
2次单元（带中间节点）
可以表现曲线变形。应力精度大幅提高，但自由度约增加2~3倍。推荐：应力评价很重要的情况。
完全积分 vs 减缩积分
完全积分：过约束（锁定）风险。减缩积分：沙漏模式（零能量模式）风险。根据情况选择。
自适应网格
基于误差指标（ZZ估计量等）的自动细化。高效提高应力集中部位的精度。有h法（单元分割）和p法（增加阶次）。
牛顿·拉夫逊法
非线性分析的标准方法。每次迭代更新切线刚度矩阵。在收敛半径内具有二次收敛性，但计算成本高。
修正牛顿·拉夫逊法
切线刚度矩阵使用初始值或每隔几次迭代更新。每次迭代成本低，但收敛速度是
相关主题
なった
詳しく
報告