压力容器线性分析
理论与物理
压力容器力学
老师，压力容器的应力分析需要特殊的理论吗？
压力容器以旋转壳体力学为基础。对于薄壁圆筒或球壳受内压的问题，通常分两个阶段考虑：膜理论（membrane theory）和弯曲理论（bending theory）。
薄壁理论（膜理论）
请告诉我薄壁圆筒的应力公式。
对于内压 $p$、内径 $D$、壁厚 $t$ 的薄壁圆筒：
周向应力是轴向应力的2倍！所以压力容器才会沿周向开裂啊。
正是如此。这和煮香肠时纵向裂开的道理一样。压力容器设计中，周向应力起主导作用，壁厚主要由周向应力决定。
对于薄壁球壳：
球壳的应力更低（是圆筒周向应力的一半）。所以球形储罐比圆筒形更高效，但制造更困难。
厚壁理论（拉梅问题）
薄壁公式能用到多大壁厚程度？
大致以 $D/t > 20$ 为基准。当 $D/t < 20$ 时，需要采用厚壁理论（拉梅公式）：
厚壁情况下 $\sigma_r$ 也不能忽略呢。内表面处 $\sigma_r = -p$（压缩）。
是的。薄壁理论假设 $\sigma_r \approx 0$，但厚壁情况下壁厚方向的应力梯度变得重要。内表面的周向应力比薄壁公式计算值更高。
不连续应力
压力容器的“不连续应力”是什么？
这是压力容器分析的核心。如果只是均匀的圆筒或球壳，膜理论就足够了，但实际的压力容器存在：
- 筒体与封头的连接处
- 接管连接处
- 壁厚变化处
- 支撑部位
在这些几何不连续处，为了满足仅靠膜理论无法满足的协调条件，会产生局部弯曲应力。这就是不连续应力。
膜应力之上，叠加了局部弯曲应力呢。
是的。不连续应力在远离不连续处时会迅速衰减。衰减距离大致是 $\sqrt{Rt}$ 的量级（$R$: 半径，$t$: 壁厚）。在这个“影响范围”内，仅用膜理论是不准确的，需要借助FEM进行弯曲分析。
ASME规范的应力分类
ASME规范是如何对应力进行分类的？
ASME BPVC Section VIII Div. 2 Part 5 的应力分类：
| 应力类别 | 符号 | 含义 | 许用值 |
|---|---|---|---|
| 总体薄膜应力 | $P_m$ | 作用于容器整体的薄膜应力 | $S$（许用应力） |
| 局部薄膜应力 | $P_L$ | 不连续处的局部薄膜应力 | $1.5S$ |
| 一次弯曲应力 | $P_b$ | 为平衡载荷所需的弯曲应力 | $1.5S$（对于 $P_L + P_b$） |
| 二次应力 | $Q$ | 具有自限性的位移引起的应力 | $3S$（对于 $P_L + P_b + Q$） |
| 峰值应力 | $F$ | 局部应力集中 | 用于疲劳评估 |
薄膜应力要求最严格，二次应力允许到3倍啊。
二次应力具有“自限性”（self-limiting）。即使发生局部屈服，由于变形受到约束，也不会导致塑性坍塌。因此设定了比薄膜应力更大的许用值。大部分不连续应力被归类为二次应力。
总结
我来整理一下压力容器的理论。
要点：
- 薄壁公式 — $\sigma_\theta = pD/(2t)$。周向应力是轴向的2倍
- 厚壁用拉梅公式 — $D/t < 20$ 时需要。壁厚方向的应力梯度
- 不连续应力 — 几何不连续处的局部弯曲。FEM是必须的
- ASME应力分类 — $P_m, P_L, P_b, Q, F$ 这5个类别
- 二次应力允许到3S — 具有自限性的应力
压力容器设计的关键是正确分类“哪里是薄膜应力，哪里是二次应力”啊。
正是如此。应力分类错误会导致结果偏于安全或偏于危险。对于压力容器工程师来说，不仅要求能用FEM计算出应力，更要求具备正确分类这些应力的能力。
压力容器设计规范的历史
压力容器的安全设计标准体系化始于1914年AMSE（美国机械工程师学会）发布《锅炉及压力容器规范（BPVC）》初版。其契机是1905年马萨诸塞州布罗克顿的锅炉爆炸事故（死亡58人）。现行的ASME BPVC Section VIII Division 2（2017版）正式认可了基于有限元法的“设计-分析（DBA）”。
各项的物理意义
- 惯性项（质量项）：$\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分析压力容器时，使用什么单元？
建模选择有三种。
| 模型 | 单元类型 | 适用场景 |
|---|---|---|
| 轴对称 | CAX8R等 | 圆筒筒体、封头、轴向接管 |
| 壳 | S4R, S8R等 | 接管连接处、整体模型 |
| 实体 | C3D20R等 | 接管连接处细节、焊接部位 |
能用轴对称的话效率最高呢。
是的。因为压力容器大部分是轴对称的，所以通常先进行轴对称分析，仅对非轴对称部分（接管、支腿）进行三维分析，这是标准做法。
应力分类线（Stress Classification Line, SCL）
如何根据FEM结果进行ASME应力分类？
设定应力分类线（SCL），将线上的应力分解为薄膜成分和弯曲成分。SCL是沿壁厚方向画一条直线，对该线上的应力分布进行积分。
对于壁厚方向的应力分布 $\sigma(y)$（$y$: 壁厚位置，$-t/2$ 〜 $t/2$）：
总应力减去薄膜和弯曲应力后剩下的就是峰值应力。应力集中和缺口效应都包含在这里。
SCL的设置位置由分析者决定。在不连续处附近画一条垂直于壁厚方向的线。ASME指南（Div. 2, Part 5, Annex 5-A）规定了SCL的设置方法。
网格要求
压力容器分析中，网格需要注意什么？
壁厚方向的网格最为重要。为了进行应力分类，需要准确捕捉壁厚方向的应力梯度。
- 轴对称模型 — 壁厚方向至少4个单元（二次单元则至少2个）
- 实体模型 — 同上
- 壳模型 — 壁厚方向自动考虑（在积分点输出应力）
不连续处的网格密度：
- 转角过渡区（封头与筒体连接处）: 单元尺寸小于壁厚的1/2
- 接管连接处: 小于接管壁厚的1/2
- 焊接部位: 小于焊脚长度的1/3
相关主题
なった
詳しく
報告