壳体屈曲
理论与物理
壳屈曲的特殊性
与板的屈曲相比，壳（曲面结构）的屈曲有什么不同？
壳屈曲是结构力学中最危险且最难预测的屈曲现象。与板屈曲的根本区别有三点。
1. 膜力-弯曲耦合。 板的面内力和面外弯曲是独立的，但壳由于曲率导致面内和面外耦合。这使得曲面本身成为刚度的来源，理论屈曲应力非常高。
2. 极端的缺陷敏感性。 圆柱壳在轴压作用下，实验只能得到理论屈曲载荷的20〜40%。这种差异在1930年代被发现，并已研究了70多年。
3. 后屈曲路径不稳定。 板在屈曲后仍能承载载荷，但许多壳在屈曲瞬间载荷会急剧下降。这就是Koiter理论所说的“不稳定对称分岔”（$b < 0$）。
实验值只有理论值的2成…这意味着如果直接使用特征值屈曲分析结果，会高估5倍，处于危险侧？
正是如此。因此，在壳屈曲中，折减系数（理论值的降低系数）成为设计的核心概念。从NASA的SP-8007（1968年）开始，到现在的ESA ECSS标准，如何合理确定这个折减系数是最大的课题。
圆柱壳的经典屈曲理论
请告诉我圆柱壳的理论屈曲应力。
承受轴压的薄壁圆柱壳的经典屈曲应力为：
（$\nu = 0.3$ 的情况）
与 $t/R$ 成比例…也就是说只由板厚和半径的比值决定。长度没有影响！
这是从短圆柱到长圆柱都广泛成立的“经典解”。屈曲模式是轴向和周向具有多个半波的菱形图案。由于许多模式都密集在相同的特征值（$0.605 Et/R$）附近，因此会发生模式密集。
模式密集和缺陷敏感性有关吗？
有直接关系。由于许多模式处于相同的能量水平，微小的缺陷就会引起模式间的相互作用，导致屈曲载荷急剧下降。这就是壳屈曲缺陷敏感性的物理机制。
承受外压的圆柱壳
外压的情况呢？
承受外压的长圆柱壳的屈曲类似于柱屈曲。临界外压为：
其中 $n$ 是周向波数。给出最小 $p_{cr}$ 的 $n$ 即为屈曲模式。
$t/R$ 的3次方！与轴压的1次方相比，对板厚的敏感性高得惊人。
是的。外压屈曲是弯曲主导的，因此依赖于 $t^3$。轻微的减薄（如腐蚀）就会导致屈曲载荷急剧下降，因此压力容器的检查中板厚管理极为重要。
折减系数
请告诉我关于NASA SP-8007的折减系数。
SP-8007是NASA于1968年发布的火箭结构设计指南，其中轴压圆柱的推荐折减系数 $\gamma$ 为：
$R/t = 500$ 的薄壳时 $\gamma \approx 0.25$…是理论值的四分之一。
相当保守。实际上，SP-8007是基于1930〜60年代实验数据下限包络线的。当时制造精度低，缺陷大。现代制造技术（精密轧制、机械加工）可以实现更高的屈曲载荷，因此SP-8007被批评为“过度保守”。
有新的标准吗？
有ESA的ECSS-E-HB-32-24A（2010年修订）以及NASA/ESA联合的DESICOS项目成果。这些标准基于实测缺陷数据进行个体评估、概率论屈曲评估，并利用VCT（振动相关技术）进行无损屈曲预测，从而推导出比SP-8007更合理的折减系数。
总结
我深刻理解了壳屈曲理论的可怕之处。
要点：
- 壳屈曲是最具缺陷敏感性的屈曲现象 — 可能降至理论值的20〜40%
- 经典解 $\sigma_{cr} \approx 0.605 Et/R$ — 优美但无法直接用于实际结构
- 模式密集是缺陷敏感性的物理原因 — 众多模式处于相同的能量水平
- 折减系数是设计的核心 — SP-8007保守，新标准正在开发中
- 仅靠特征值屈曲分析完全不够 — 必须进行考虑缺陷的非线性分析
也就是说，在壳屈曲中，绝对不能直接将特征值屈曲分析结果用作设计值，对吧？
绝对不能。如果仅凭特征值分析就对壳屈曲设计放行，那意味着结构力学素养的缺失。折减系数或非线性分析，两者必选其一。
薄壁壳屈曲的折减系数问题
圆柱薄壁壳的实验屈曲载荷仅为线性理论值的25〜80%，这种差异自1930年代起就以“折减系数”问题而闻名。Lorenz、Timoshenko、Southwell的理论公式预测值是实验值的3〜4倍，因此NASA SP-8007（1968年）采用γ=0.6的经验公式修正，大幅降低了设计载荷。
各项的物理意义
- 惯性项（质量项）：$\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分析策略
用FEM分析壳屈曲时，与板屈曲有什么不同？
在三个方面比板屈曲困难得多。
1. 网格要求严格。 屈曲波长可能达到板厚量级，需要非常精细的网格。
2. 初始缺陷的建模主导结果。 板受初始缺陷影响小，但壳的屈曲载荷会因缺陷模式和振幅而改变数倍。
3. 非线性分析收敛困难。 伴随回弹行为，即使使用Riks法也容易追踪失败。
特征值屈曲分析的用法
您说壳不能用特征值屈曲分析，是完全不用吗？
只是“不能直接用作设计值”，特征值分析本身扮演重要角色：
1. 确认屈曲载荷的上限值 — 如果非线性分析结果超过此值，说明有问题
2. 利用模态形状作为初始缺陷 — 非线性分析的预处理
3. 识别临界部位 — 筛选何处容易屈曲
4. 折减系数的起点 — 用 $P_{cr,FEM} \times \gamma$ 估算设计载荷
缺陷建模方法
壳的初始缺陷怎样引入最好？
有四种方法。按可信度从高到低排列：
1. 实测缺陷
对实际结构进行激光扫描，将理想形状与实际形状的差值直接作为初始缺陷输入。可信度最高，但需要试制品。
2. 特征值模态叠加（EIMP法）
最常用的方法。叠加线性屈曲的低阶模态来创建初始形状：
振幅通常为板厚 $t$ 的0.5〜2.0倍。问题在于最不利的振幅组合未知。
3. SPLA（单扰动载荷法）
在壳的一点施加集中载荷，人为制造凹陷，然后对该凹陷结构进行非线性分析。由DLR（德国航空航天中心）的Hühne（2008）提出。物理意义明确，再现性高。
4. 随机缺陷
使用傅里叶级数的随机系数生成缺陷。与蒙特卡洛模拟结合用于概率论评估。
相关主题
なった
詳しく
報告