基尔霍夫-洛夫薄壳理论
基尔霍夫-洛夫薄壳理论的理论基础
薄壳的经典理论
老师，基尔霍夫-洛夫壳理论是基尔霍夫板理论的曲面版本吗？
正是如此。基尔霍夫板处理平板的弯曲，而基尔霍夫-洛夫壳处理曲面的弯曲与膜力的耦合。洛夫（1888）将基尔霍夫假设扩展到了曲面。
基本假设
洛夫假设（一阶近似理论）：
1. 壳的板厚 $t$ 相对于其他尺寸（半径 $R$ 等）足够薄
2. 位移相对于板厚很小
3. 忽略板厚方向的法向应力 $\sigma_z$
4. 忽略板厚方向的剪切变形 —— 与基尔霍夫假设相同
5. 忽略 $t/R$ 的高阶项
假设4和基尔霍夫板一样呢。剪切变形为零。
是的。与板（平面）的本质区别在于曲率引起的膜力-弯曲耦合。平板中膜力与弯曲是独立的，但在壳中，由于曲率，两者发生耦合。这正是壳理论复杂性的根源。
膜力与弯曲的耦合
能再详细说明一下"膜力-弯曲的耦合"吗？
考虑一个受压力的球壳。膜理论中会产生 $\sigma = pR/(2t)$ 的均匀拉应力。但如果球壳上有孔洞，或者板厚发生变化，仅靠膜应力无法满足变形的协调条件。这个不足的部分由弯曲力矩来弥补。
和板的"不连续应力"一样呢。
正是相同的机理。压力容器筒体与封头连接处产生不连续应力，这可以用基尔霍夫-洛夫理论中曲率变化引起的膜-弯曲耦合来解释。
FEM中的实现
有基尔霍夫-洛夫壳的FEM单元吗？
由于和基尔霍夫板同样的原因，需要$C^1$ 连续性，因此实现困难。现代FEM中，包含剪切变形的明德林-赖斯纳壳是主流。
不过IGA（等几何分析）可以利用NURBS基函数的 $C^1$ 连续性，直接离散化基尔霍夫-洛夫壳。IGA壳单元研究活跃正是因为这个原因。
总结
我来整理一下基尔霍夫-洛夫壳理论。
要点：
- 忽略剪切变形的薄壳经典理论 —— Love (1888)
- 膜力与弯曲的耦合 —— 壳的曲率引起耦合
- FEM需要 $C^1$ 连续性 —— 实现困难
- 实务中用明德林-赖斯纳壳替代 —— 薄壳时收敛于K-L理论
- IGA正在复兴K-L壳 —— NURBS基函数的 $C^1$ 连续性
和基尔霍夫板是同样的模式呢。理论优美但FEM实现困难，主角让给了明德林系。不过IGA有复兴的迹象。
正是如此。结构力学理论的深度与FEM实现难度之间的鸿沟，正在被IGA这一新范式所填补。
基尔霍夫-洛夫理论的假设
基尔霍夫-洛夫（Kirchhoff-Love）壳理论由古斯塔夫·基尔霍夫于1850年作为板理论提出，1888年由奥古斯塔斯·洛夫扩展到曲面壳。主要假设是"壳中面的法线在变形后仍保持为法线（法线不变）"，此假设使得可以忽略剪切变形。适用于板厚/跨度比1/20以下的薄壁结构。
数值解法与实现
K-L壳的数值解法
请告诉我用FEM处理基尔霍夫-洛夫壳的方法。
有三种方法。
1. 用明德林-赖斯纳壳单元替代
最实用。明德林系壳单元（S4R, CQUAD4等）在薄壳极限下会自动收敛于K-L壳。若有剪切锁定对策（MITC法等）则没问题。
2. DKT壳单元
离散基尔霍夫三角形（DKT）是将板理论的DKT扩展到壳的单元。Nastran的CTRIA3壳和Abaqus的STRI3/STRI65属于此类。专用于薄板，精度高。
3. IGA壳单元
等几何分析（IGA）利用NURBS基函数的 $C^1$ 连续性直接离散化K-L壳。无需旋转自由度，可减少DOF。
IGA壳没有旋转自由度？
在K-L壳中，转角由挠度的微分决定，所以独立变量只有位移 $(u, v, w)$ 三个。比明德林壳的5自由度（$u, v, w, \theta_x, \theta_y$）或6自由度（含drilling）要少。减少DOF数在大规模问题中是一大优势。
各方法比较
| 方法 | 连续性 | DOF/节点 | 薄板精度 | 厚板对应 |
|---|---|---|---|---|
| 明德林系 | $C^0$ | 5〜6 | ○（带锁定对策） | ○ |
| DKT壳 | $C^0$（离散K-L） | 5〜6 | ◎ | × |
| IGA K-L壳 | $C^1$ 以上 | 3 | ◎ | × |
IGA的DOF效率最高呢。
但IGA在商用求解器中的实现还不普及。LS-DYNA中实现了一部分IGA壳，但Abaqus/Nastran/Ansys中还不是标准。目前明德林系是实务标准。
总结
我来整理一下K-L壳的数值方法。
要点：
- 用明德林系壳替代是实务标准 —— 薄壳时自动收敛于K-L
- DKT壳是薄板专用的高精度单元 —— STRI3/STRI65
- IGA使K-L壳的直接离散化成为可能 —— DOF为3，效率高
- IGA的商用实现仍有限 —— 未来的主力候选
KL壳的C1连续性要求
基尔霍夫-洛夫理论中位移的二阶导数出现在弱形式中，因此有限元需要在相邻单元间满足C1连续性（位移及其一阶导数的连续性）。满足此要求的单元设计很复杂，1960~70年代许多研究者曾挑战过。Bogner-Fox-Schmit（BFS）矩形单元是C1单元的经典例子，但由于难以适应任意网格，目前主流已由明德林-赖斯纳壳（C0连续）替代。