基尔霍夫板理论
理论与物理
什么是基尔霍夫板理论？
老师，基尔霍夫板理论是欧拉-伯努利梁理论的二维版本吗？
正是如此。正如欧拉-伯努利梁假设“横截面始终垂直于中性轴”一样，基尔霍夫板理论假设“板厚方向的直线在变形后仍保持为直线，且垂直于中面”。
基本假设
基尔霍夫假设：
1. 直法线假设 — 变形前垂直于中面的直线，变形后仍垂直于中面
2. 法线无伸缩假设 — 板厚方向的应变 $\varepsilon_{zz} = 0$
3. 板厚方向剪切应变为零 — $\gamma_{xz} = \gamma_{yz} = 0$
假设3和欧拉-伯努利梁一样呢。忽略了剪切变形。
是的。这个假设使得转角由面外挠度的微分决定：
自由度只有挠度 $w(x,y)$ 这一个，转角不是独立变量。
控制方程
板挠度 $w(x,y)$ 的双调和方程：
其中 $D = Et^3/(12(1-\nu^2))$ 是板的弯曲刚度，$q$ 是面外分布载荷。
$\nabla^4$ 是四阶微分算子。是梁的 $EI w'''' = q$ 的二维版本呢。
没错。因为 $\nabla^4 = \nabla^2(\nabla^2)$，所以：
弯矩与剪力
内力分量：
弯矩由 $w$ 的二阶微分决定。挠度求两次导就是弯矩。和梁的结构一样呢。
是的。板的弯曲应力为：
在板表面（$z = \pm t/2$）达到最大应力。和梁的 $\sigma = My/I$ 结构相同。
适用范围
基尔霍夫板理论能用到多厚的板？
前提是薄板。经验法则是 $b/t > 20$（$b$: 板短边，$t$: 板厚）。
与铁木辛柯梁结构类似：
- $b/t > 20$: 基尔霍夫板足够
- $10 < b/t < 20$: 考虑使用明德林板
- $b/t < 10$: 明德林板或实体单元
明德林板是铁木辛柯梁的二维版本吗？
正是如此。基尔霍夫板 = EB梁的二维版，明德林板 = 铁木辛柯梁的二维版。区别在于是否考虑剪切变形。
FEM中的基尔霍夫板单元
在FEM中实现基尔霍夫板困难吗？
实际上非常困难。基尔霍夫板理论包含 $w$ 的四阶微分，因此在FEM中实现需要 $C^1$ 连续性（位移和转角在单元间都连续）。通常的FEM（$C^0$ 连续性）无法满足这一点。
$C^1$ 连续性很难吗？
在二维中创建满足 $C^1$ 连续性的多项式单元很困难。历史上开发了Argyris三角形（21自由度）或Bell三角形（18自由度），但自由度多，不实用。因此明德林板理论（只需 $C^0$ 连续性）在FEM中成为了主流。
总结
我来整理一下基尔霍夫板理论。
要点：
- 忽略剪切变形的薄板弯曲理论 — EB梁的二维版
- $D\nabla^4 w = q$ — 双调和方程
- 适用于 $b/t > 20$ — 仅限于薄板
- FEM中需要 $C^1$ 连续性 — 实现困难
- 实际应用中明德林板（只需 $C^0$ 连续性）是主流 — 基尔霍夫板单元很少见
理论虽然优美，但因为FEM实现困难，所以把主角位置让给了明德林板呢。
是的。不过基尔霍夫板理论是理论解的基础，作为验证FEM结果的参考解是不可或缺的。纳维解（矩形板的双重傅里叶级数解）是基尔霍夫理论的经典解法。
基尔霍夫板理论的起源
古斯塔夫·基尔霍夫在1850年的论文“Über das Gleichgewicht und die Bewegung einer elastischen Scheibe”中确立了板的弯曲理论。他的假设（薄板、直法线保持、中面应变忽略）作为“经典板理论（CPT）”至今仍通用。数值实验证实，在轧制铝薄板（t<2mm）的成形分析中，该理论的误差小于1%。
各项的物理意义
- 惯性项（质量项）：$\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统一 |
数值解法与实现
FEM中基尔霍夫板的实现
$C^1$ 连续性问题是如何解决的？
历史上有三种方法。
1. 高阶协调单元
2. DKT/DKQ单元（离散基尔霍夫）
DKT（离散基尔霍夫三角形），在明德林板理论框架内离散化，在高斯积分点“离散地”满足基尔霍夫约束（剪切应变 = 0）。
用明德林板单元，之后再施加基尔霍夫条件…很巧妙呢。
DKT由Batoz, Bathe, Ho于1980年提出。3节点9自由度（每个节点 $w, \theta_x, \theta_y$），数量少，精度也高。Nastran的CTRIA3（弯曲）内部使用了DKT系的公式。
DKQ（离散基尔霍夫四边形）是4节点的四边形版本。同样实用。
3. 将明德林板单元作为薄板使用
最实用的方法。使用明德林板单元（包含剪切变形），如果板很薄，剪切变形会自动变小。对于薄板，能得到与基尔霍夫板相同的结果。
那么，现代FEM中，主流是不使用基尔霍夫板的“专用单元”，而是用明德林板单元替代吗？
没错。Abaqus、Ansys、Nastran的通用壳单元都是基于明德林（Reissner-Mindlin）的。对于薄板，会收敛到基尔霍夫的理论解。
理论解的活用
基尔霍夫板的理论解如何用于FEM验证？
最著名的理论解是纳维解。四边简支矩形板在均布载荷 $q$ 作用下的挠度：
双重傅里叶级数…
なった
詳しく
報告