明德林-赖斯纳板理论
理论与物理
铁木辛柯梁的二维版
老师，明德林板是给基尔霍夫板加上剪切变形吗？
正是如此。正如铁木辛柯梁给欧拉-伯努利梁添加了剪切变形一样，明德林板（Reissner-Mindlin板）是给基尔霍夫板添加了剪切变形的理论。
基本假设
明德林板的假设：
1. 板厚方向的直线在变形后仍保持直线（但不一定与中面垂直）
2. 板厚方向的应变 $\varepsilon_{zz} = 0$
3. 考虑剪切变形 — $\gamma_{xz} \neq 0, \gamma_{yz} \neq 0$
与基尔霍夫板的区别在于假设3呢。剪切应变不为零。
这个区别使得转角 $\theta_x, \theta_y$ 与挠度 $w$ 的微分相互独立：
自由度不仅是 $w$，$\theta_x, \theta_y$ 也独立。每个节点有3个自由度呢。
是的。基尔霍夫板中有 $\theta = -\partial w / \partial x$ 的约束，但明德林板没有这个约束。由于这种独立性，使得FEM的实现只需 $C^0$ 连续性即可。
控制方程
明德林板的平衡方程是三个联立偏微分方程：
其中 $Q_x, Q_y$ 是横向剪切力，$M_x, M_y, M_{xy}$ 是弯曲/扭转力矩。
基尔霍夫板是 $w$ 的一个四阶偏微分方程，而明德林板是联立方程…。
但相应地，每个方程只包含二阶及以下的微分。这就是为什么可以用 $C^0$ 连续的FEM单元进行离散化。
剪切闭锁的问题
明德林板单元也有剪切闭锁问题吧。
会发生与铁木辛柯梁完全相同的问题。对于薄板（$b/t > 20$），理论上剪切变形几乎为零，但在通常的FEM单元中剪切变形不会消失，导致单元变得过于刚硬。
已开发出多种对策：
- MITC单元（Mixed Interpolation of Tensorial Components） — Bathe-Dvorkin的方法。独立假设剪切应变的插值
- DSG法（Discrete Shear Gap） — 剪切间隙法
- 减缩积分 — 适用于一阶单元
- Assumed Natural Strain (ANS) — 在自然坐标系中假设剪切应变
我听说过MITC单元。
MITC4（4节点）和MITC9（9节点）是Bathe教授开发的单元，在消除剪切闭锁的同时也能通过分片试验。许多商用求解器的壳单元都基于MITC。Abaqus的S4R、Nastran的CQUAD4（壳单元）、Ansys的SHELL181都包含了MITC系的技术。
总结
我来整理一下明德林板理论。
要点：
- 考虑剪切变形的板理论 — 铁木辛柯梁的二维版
- $w, \theta_x, \theta_y$ 是独立变量 — 可用 $C^0$ FEM单元离散化
- 薄板时收敛于基尔霍夫板 — $b/t > 20$ 时一致
- 剪切闭锁是最大课题 — 用MITC法、DSG法、减缩积分等对策
- 现代的FEM壳单元全部基于明德林理论 — 事实上的标准理论
既然使用FEM的壳单元，明德林板理论就是无法避开的基础知识呢。
正是如此。壳单元设置中“为什么推荐减缩积分”、“为什么薄板也能使用壳单元”等问题，都可以用明德林板理论和剪切闭锁的知识来解释。
明德林板理论的成立
Raymond D. Mindlin在1951年的论文《Influence of Rotatory Inertia and Shear on Flexural Motions of Isotropic, Elastic Plates》中发表了考虑剪切变形的板理论。E. Reissner也在1944年独立发表了类似的公式化，现在常并称为“Mindlin-Reissner理论”。与基尔霍夫理论相比，该理论在t/L>0.05的板中精度提高了10~20%。
各项的物理意义
- 惯性项（质量项）：$\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统一 |
数值解法与实现
MITC单元的原理
请讲解一下MITC单元消除剪切闭锁的机制。
通常的明德林板单元用相同的形函数对 $w$ 和 $\theta$ 进行插值。对于薄板，应有 $\gamma = \partial w / \partial x + \theta_y \approx 0$，但同阶插值在数值上不会变为零。
MITC法是独立地插值剪切应变。不是从 $w$ 和 $\theta$ 计算剪切应变，而是在边上的“连接点”处假设剪切应变，并在单元内进行内插。
“连接点”是什么？
在单元边上的特定点评估剪切应变，并将该值内插到整个单元。这样就能消除寄生剪切（闭锁的原因）。因为是“连接”（tie）边上的点，所以称为连接点。
各求解器的壳/板单元
请介绍一下各求解器的明德林板（壳）单元。
| 单元 | 求解器 | 节点数 | 闭锁对策 | 推荐度 |
|---|---|---|---|---|
| S4R | Abaqus | 4 | 减缩积分+沙漏控制 | ○ 通用推荐 |
| S4 | Abaqus | 4 | 完全积分+非协调模式 | ○ 精密分析 |
| S8R | Abaqus | 8 | 减缩积分 | ◎ 最高精度 |
| CQUAD4 | Nastran | 4 | MITC系+非协调 | ◎ 行业标准 |
| CQUAD8 | Nastran | 8 | 二阶单元 | ○ |
| SHELL181 | Ansys | 4 | MITC系 | ○ 通用推荐 |
| SHELL281 | Ansys | 8 | MITC系 | ◎ 最高精度 |
Abaqus的S4R和S4有什么区别？
なった
詳しく
報告