Mindlin-Reissner板理论
Mindlin-Reissner板理论的理论基础
Timoshenko梁的二维版本
老师，Mindlin板是Kirchhoff板加上剪切变形吗？
没错。正如Timoshenko梁在Euler-Bernoulli梁基础上添加了剪切变形，Mindlin板（Reissner-Mindlin板）在Kirchhoff板基础上添加了剪切变形理论。
基本假设
Mindlin板的假设：
1. 板厚方向的直线在变形后仍为直线（但不一定垂直于中面）
2. 板厚方向应变 $\varepsilon_{zz} = 0$
3. 考虑剪切变形 — $\gamma_{xz} \neq 0, \gamma_{yz} \neq 0$
与Kirchhoff板的区别是假设3对吧。剪切应变不为零。
这个区别使得旋转角 $\theta_x, \theta_y$ 与挠度 $w$ 的微分独立：
自由度不仅是 $w$，而是 $\theta_x, \theta_y$ 也是独立的。每个节点有3个自由度。
对的。在Kirchhoff板中，$\theta = -\partial w / \partial x$ 有约束，但在Mindlin板中没有这个约束。这种独立性意味着FEM实现可以只用 $C^0$ 连续单元。
控制方程
Mindlin板的平衡方程是三个联立偏微分方程：
其中 $Q_x, Q_y$ 是横向剪切力，$M_x, M_y, M_{xy}$ 是弯曲/扭转矩。
Kirchhoff板是关于 $w$ 的4阶偏微分方程，而Mindlin板是联立方程…
相反地，每个方程只包含2阶及以下的微分。这就是为什么可以用 $C^0$ 连续的FEM单元进行离散化。
剪切自锁问题
Mindlin板单元也有剪切自锁问题吗？
Timoshenko梁的问题完全相同。对于薄板（$b/t > 20$），理论上剪切变形应接近零，但普通FEM单元剪切变形不消失，导致单元过于刚硬。
已开发多种对策：
- MITC单元（混合张量分量插值） — Bathe-Dvorkin方法。独立假设剪切应变补间
- DSG法（离散剪切间隙） — 剪切间隙法
- 低减积分 — 应用于一次单元
- 假设自然应变 (ANS) — 在自然坐标系中假设剪切应变
我听过MITC单元。
MITC4（4节点）或MITC9（9节点）是Bathe教授开发的单元，消除了剪切自锁同时通过了补丁测试。许多商用求解器的壳单元基于MITC技术。Abaqus的S4R、Nastran的CQUAD4（壳）、Ansys的SHELL181都包含MITC系技术。
总结
总结一下Mindlin板理论。
要点：
- 考虑剪切变形的板理论 — Timoshenko梁的二维版本
- $w, \theta_x, \theta_y$ 是独立变量 — 可用 $C^0$ FEM单元离散化
- 薄板时收敛到Kirchhoff板 — 当 $b/t > 20$ 时一致
- 剪切自锁是最大课题 — 用MITC法、DSG法、低减积分应对
- 现代FEM壳单元都基于Mindlin — 事实上的标准理论
只要使用FEM壳单元，就避免不开Mindlin板理论这个基础知识。
没错。壳单元设置中"为什么推荐低减积分""为什么薄板也能用壳单元"等问题，都可以用Mindlin板理论和剪切自锁知识解释。
Mindlin板理论的确立
Raymond D. Mindlin在1951年的论文"Influence of Rotatory Inertia and Shear on Flexural Motions of Isotropic, Elastic Plates"中发表了考虑剪切变形的板理论。E. Reissner也在1944年独立发表了类似的公式，所以现在称为"Mindlin-Reissner理论"。与Kirchhoff理论相比，当t/L>0.05时，该理论的精度提高10～20%。
Mindlin-Reissner板理论的数值计算方法
MITC单元的原理
请教我MITC单元如何消除剪切自锁。
在普通Mindlin板单元中，$w$ 和 $\theta$ 用相同的形状函数插值。对于薄板，$\gamma = \partial w / \partial x + \theta_y \approx 0$ 应成立，但同次插值不能数值上为零。
MITC法独立插值剪切应变。不是从 $w$ 和 $\theta$ 计算剪切应变，而是在"约束点"处假设边上的剪切应变，然后内插到单元内。
"约束点"是什么意思？
单元边上特定点处评估剪切应变，然后在整个单元内插该值。这样消除了寄生剪切（自锁的原因）。从"结合"（tie）这些点的概念来的"约束点"名称。
各求解器的壳/板单元
请教我各求解器的Mindlin板（壳）单元。
| 单元 | 求解器 | 节点数 | 自锁对策 | 推荐度 |
|---|---|---|---|---|
| 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有什么区别？
实务推荐：
- 一般分析 → S4R（速度和精度的平衡）
- 精密应力评估 → S8R（二次单元。最高精度）
- 座屈分析 → S4R（也支持显式法结合）
板厚方向的积分点
壳单元在板厚方向也有积分点对吧。
Mindlin板单元在板厚方向有多个Simpson积分点。默认通常为5点（Abaqus）或5～7点（Nastran/Ansys）。
板厚方向的积分点数产生影响的场景：
- 弹性分析 → 3点足够
- 弹塑性分析 → 5点以上。追踪板厚方向屈服的进展
- 复合材料分层 → 每层最少3点。全层15～20点以上
总结
整理一下Mindlin板的数值方法。
要点：
- MITC法是剪切自锁的主要对策 — 已内置于主要求解器
- S4R/CQUAD4/SHELL181是通用推荐 — 4节点低减积分＋自锁对策
- S8R/CQUAD8/SHELL281是精密分析 — 8节点二次单元
- 板厚方向的积分点 — 弹性3点，弹塑性5点以上
- 壳单元设置因求解器而异 — 查阅手册确认
剪切自锁与DSG法
Mindlin板单元4节点版在薄板极限会出现剪切自锁。1985年Bleylock和Hughes提出DSG I（离散剪切间隙）法，通过在单元边上独立评估剪切应变作为独立变量来消除自锁。该方法从1990年代起在NX Nastran的CQUAD4（BCSCALE=1.0）中实现，大幅改善了剪切变形和弯曲的双重精度。
Mindlin-Reissner板理论的实务应用
壳单元的实务应用
Mindlin板基础的壳单元在实务中如何应用？
这是FEM结构分析中最常用的单元类型。薄板、钢结构、车体、飞机、船体…几乎所有薄壁结构都用壳单元分析。
壳单元的应用范围
| $b/t$ 范围 | 推荐单元 | 理由 |
|---|---|---|
| > 100 | 壳（S4R等） | 薄板。剪切变形可忽略 |
| 20 〜 100 | 壳（S4R等） | 标准应用范围 |
| 10 〜 20 | 壳或实体 | 两者都可以 |
| < 10 | 实体 | 厚板。壳的假设失效 |
当 $b/t < 10$ 时应该用实体单元？
当"板厚方向的应力梯度很重要"的问题中需要用实体。螺栓紧固部、厚翼缘接合、接触压力评估部等。
壳单元的偏移
壳单元的"偏移"是什么意思？
壳单元网格定义在中面上，但实结构有时希望网格在板的表面（上面或下面）。偏移指定从中面的距离。
典型用法：
- T字接合部 — 翼缘和腹板的连接。翼缘的中面和腹板的中面有偏移
- 复合材分层 — 每层的中面不同
- 焊接部 — 焊线位置与壳单元中面不同
偏移设置错误会怎样？
对结果有很大影响。特别是T字接合部，不设置偏移会低估翼缘的弯曲刚性。
网格大小的确定
壳单元的网格大小：
- 板厚的5倍以上 — 通常最小尺寸（$h_{elem} > 5t$）
- 特征长度 — 孔的半径/8，圆角半径/3
- 座屈半波长的6分之1以上 — 表示座屈模式所必需
单元大小至少是板厚的5倍…那1 mm板厚的板需要至少5 mm的单元吗？
没错。壳单元的网格大小对板厚有约束。小于板厚的单元在物理上无意义（违反壳理论假设）。需要比板厚更精细的精度时用实体单元。
实务检查清单
请提供壳单元的检查清单。
壳单元是结构FEM的主力，所以这个检查清单是最常用的。
没错。能否正确设置壳单元是结构FEM工程师的基本技能。
船体甲板的结构分析
大型LNG运输船（载货量26万立方米）的甲板结构分析中，对板厚20～35 mm的甲板钢板应用Mindlin-Reissner板单元。2016年三菱重工技术报告显示，COMSOL的Mindlin板单元模型（约50万单元）计算波浪荷载下的主应力，与应变计测量值的误差在±8%以内。
Mindlin-Reissner板理论的软件比较
壳单元的求解器比较
请比较各求解器的壳单元。
为什么Nastran的CQUAD4在航空航天领域压倒性占优？
Nastran的CQUAD4以1980年代MacNeal-Harder单元为基础，经过多年改进。在航空航天认证（DO-160等）中有大量验证实绩。特别是与PCOMP的组合在复合材面板座屈、振动中积累了无法复制的经验。
选择指南
壳单元在各求解器上都是成熟单元。
没错。壳单元是各求解器投入最多开发力量的单元之一。求解器间的差异在逐年缩小。
Mindlin板单元的实现差异比较
Abaqus S4R、ANSYS SHELL181、NX Nastran CQUAD4（默认）、LS-DYNA ELFORM=2各自采用不同的剪切自锁对策。S4R采用4点积分＋沙漏控制，SHELL181采用B-bar法，NX Nastran采用DSG法。2019年SolverBench均一弯曲问题中S4R和CQUAD4的最大弯矩差为0.7%，实务上基本无差异。
Mindlin-Reissner板理论的先端研究
Mindlin板的先端研究
壳单元的最新研究动向是什么？
壳单元是研究最活跃的领域之一。
实体壳单元
实体壳单元具有普通实体单元（HEX8等）的形状，但达到壳单元的精度。板厚方向可用1个单元表示弯曲。
有什么优点？
Abaqus的SC8R（8节点实体壳）或Ansys的SOLSH190属于这类。
IGA壳
等几何分析（IGA）基础的壳单元直接使用CAD的NURBS曲面。无需网格生成，形状近似误差为零。薄板壳座屈分析是IGA优势最明显的领域。
连续统壳
用超薄三维实体单元的方法。板厚方向1～2个实体单元来解板弯曲。用增强假设应变（EAS）或ANS法回避自锁。未来壳单元的概念可能本身就会消失。
总结
总结壳单元的先端研究。
只要有"薄壁结构的高效建模"这个本质需求，壳单元就会继续是FEM的核心。
Mindlin理论对复合材分层板的扩展
纤维增强复合材（CFRP）分层板分析使用扩展的Mindlin理论、一阶剪切变形理论（FSDT）。2010年代波音787机身的碳纤维分层面板（t=12 mm、拟各向同性层压）座屈分析中，Abaqus S4R计算的座屈荷载与实验值的差为98.7%，与层状理论（更精密）的差仅1.3%。
Mindlin-Reissner板理论的故障处理
壳单元的故障
（Mindlin板基础）壳单元常见的故障有哪些？
壳单元是最广泛使用的单元，所以故障也很多样。
剪切自锁
薄板的挠度比理论值小。
完全积分的壳单元（S4等）偶尔会发生剪切自锁。对策：
沙漏模式
S4R的变形呈锯齿形。
低减积分的壳单元激发了沙漏模式。
原因：
- 集中荷载
- 网格过粗
- 荷载偏向单元的一条边
对策：
- 分散荷载（用RBE3等）
- 细化网格
- 增加Abaqus沙漏刚性
- 切换到S4（完全积分）
法向方向不一致
相邻单元的应力结果符号相反。
单元法向方向反向。壳单元有"上面""下面"的区别，法向方向（正z方向）每个单元不同时，弯曲应力符号反向。
对策：
- 用前处理器统一全部单元的法向方向
- Abaqus中用 *NORMAL 明确指定法向
- Nastran中用 PARAM,SNORM 自动调整法向
这是容易被忽视的故障。
法向不统一表现为力矩图的不连续。力矩符号在单元间反向就该检查法向。
板厚变化部位的应力
板厚变化处的应力有问题。
壳单元的板厚在各单元定义。板厚急变部位会出现：
- 膜力不连续（相同膜应力但板厚不同时膜力不连续）
- 中面偏移（板厚变化使中面偏移）
对策：
- 在板厚变化部正确设置偏移
- 使板厚变化平缓过渡（渐变）
- 必要时切换为实体单元
总结
整理壳单元的故障处理。
法向不一致是盲点。应力符号反向这样…
壳单元的"表裏"特性与2D平面单元或实体单元本质不同。常记"表裏"是掌握壳单元的关键。
旋转DOF问题的处理
Mindlin板单元缺面内旋转（旋转DOF）的单元用在非平面网格时，节点在面外方向"摇晃"的数值奇点导致的面外奇点问题会发生。ANSYS SHELL181的KEYOPT(3)=0自动补强旋转DOF，但过度补强导致高纵横比单元的面内应力偏差5～15%的案例记录在2008年ANSYS检证事例集。
相关主题
价值
更详
要报告