基尔霍夫-洛夫薄壳理论
理论与物理
薄壳的经典理论
老师，基尔霍夫-洛夫壳理论是基尔霍夫板理论的曲面版本吗？
正是如此。基尔霍夫板处理平板的弯曲，而基尔霍夫-洛夫壳处理曲面的弯曲与膜力的耦合。洛夫（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以下的薄壁结构。
各项的物理意义
- 惯性项（质量项）：$\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统一 |
数值解法与实现
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连续）替代。
线性单元（1次单元）
节点间线性插值。计算成本低，但应力精度低。注意剪切锁定（可通过减缩积分或B-bar法缓解）。
2次单元（带中间节点）
可表现曲线变形。应力精度大幅提高，但自由度约增加2~3倍。推荐：应力评估重要时使用。
完全积分 vs 减缩积分
完全积分：过约束（锁定）风险。减缩积分：沙漏模式（零能量模式）风险。根据情况选择。
自适应网格
基于误差指标（ZZ估计量等）的自动细化。高效提高应力集中部位的精度。有h法（单元细分）和p法（增加阶次）。
牛顿·拉夫逊法
非线性分析的标准方法。每次迭代更新切线刚度矩阵。在收敛半径内具有二次收敛性，但计算成本高。
修正牛顿·拉夫逊法
切线刚度矩阵使用初始值或每隔几次迭代更新。每次迭代成本低，但收敛速度为线性。
收敛判定标准
力残差范数: $||R|| / ||F_{ext}|| < \epsilon$（通常 $\epsilon = 10^{-3}$〜$10^{-6}$）。位移增量范数: $||\Delta u|| / ||u|| < \epsilon$。能量范数: $\Delta u \cdot R < \epsilon$
载荷增量法
不一次性施加全部载荷，而是分小步增加。弧长法（Riks法）可以越过载荷-位移关系的极值点进行追踪。
直接法 vs 迭代法的比喻
直接法是“用笔算精确解联立方程”的方法——可靠但大规模问题太耗时。迭代法是“反复猜测逼近正确答案”的方法——最初答案粗糙，但每次迭代精度提高。查字典时，与其从第一页开始按顺序找（直接法），不如先估计大概位置（迭代法）更快。FEM的迭代求解器（CG法、GMRES法）就是这种“聪明的猜测”。
なった
詳しく
報告