压力容器的线性解析
压力容器线性理论基础
压力容器的力学
老师，压力容器的应力解析需要特殊的理论吗？
压力容器是回转壳的力学基础。薄肉的圆筒或球壳在内压作用下的问题，分两阶段考虑：膜理论（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$ 五个分类
- 二次应力许可3S — 具有自制性的应力
压力容器设计的关键是"正确区分哪些是膜应力，哪些是二次应力"。
恰恰如此。应力分类错误时，可能既过度安全也可能不安全。FEM输出应力只是第一步，正确分类应力的能力是压力容器工程师的必备技能。
压力容器设计规范的历史
压力容器安全设计基准的体系化开始于1914年，当时美国机械学会（ASME）发表了《锅炉和压力容器规范》（Boiler and Pressure Vessel Code, BPVC）的初版。其背景是1905年美国马萨诸塞州Brockton锅炉爆炸事故（58人死亡）。现在的ASME BPVC Section VIII Division 2（2017版）正式承认基于有限要素法的"设计与分析"（Design by Analysis, DBA）方法。
压力容器线性数值计算方法
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
Abaqus
```
*STEP
*STATIC
*DLOAD
cylinder, P, 10.0 $ 内压 10 MPa
*END STEP
```
内压通过 *DLOAD 的 P（压力）选项直接施加到面。方向为单元法向。
Nastran
```
PLOAD4, 100, elem_id, 10.0, , , , , 1
```
用PLOAD4卡定义面压。最后的标志位指定法向方向。
Ansys
```
SF, area_num, PRES, 10.0
```
或在Workbench中用Pressure荷载选择面并输入压力值。
应力线性化工具
应力分类是手动做吗？
有应力线性化功能的后处理器很方便：
- Abaqus/CAE — 用Path功能定义SCL，自动分离膜/弯曲
- Ansys Workbench — Linearized Stress功能直接输出SCL上的应力
- Nastran — 没有直接功能，从f06应力数据手动或脚本计算
Ansys Workbench的功能最方便似的。
Ansys Workbench的Linearized Stress在压力容器工程师中很受欢迎。可视化设置SCL，自动输出膜、弯曲、峰值各成分和等效应力（von Mises/Tresca）。
总结
整理压力容器的数值方法。
要点：
- 轴对称+3D子模型最高效 — 不需要全部3D
- 板厚方向最少4个单元 — 为应力分类捕捉梯度
- SCL分离膜、弯曲、峰值 — ASME Div. 2应力分类必需
- 利用后处理器的应力线性化功能 — Ansys Workbench方便
- 不连续部网格要充分细 — 皱纹、喷嘴、焊接部
拉梅应力式和FEM验证
厚肉圆筒的内压应力存在拉梅（Gabriel Lamé, 1852年）的严格解，历来作为FEM验证问题经典使用。内径50mm、外径100mm、内压10MPa的钢管，最大径向应力由拉梅解为13.3MPa。此问题用Abaqus CAX4单元求解，单元数20时误差0.3%以内收敛，已列入著名FEM教科书（Zienkiewicz 2000版）作为基准。
压力容器线性实务应用
压力容器设计的实务流程
压力容器设计在实际中怎样进行？
ASME BPVC Section VIII 中有两种设计方法。
Division 1（按规则设计）
基于公式和图表的设计。无需FEM。各主要部件（胴、镜板、喷嘴、法兰）都有设计式。简单形状足够了。
Division 2（按分析设计）
基于FEM应力解析的设计。两阶段：Part 4（按规则设计）和Part 5（按分析设计）。Part 5中：
- 弹性应力解析法 — FEM线性弹性结果应力分类后与许可值比较
- 弹塑性解析法 — FEM非线性弹塑性解析直接评估塑性崩溃荷载
Division 1和Division 2的选择？
Division 1简单安全，但板厚倾向于厚。Division 2用FEM精密评估，可减薄板厚，但解析成本高。高压、大型容器或需要轻量化时选Division 2。
喷嘴接续部的解析
喷嘴接续部是压力容器最重要的部位吧。
喷嘴接续部是应力集中最大的部位，容易成为破坏和疲劳的起源。
解析的要点：
- WRC 297/537 — 喷嘴接续部应力的简化计算方法。提供局部应力系数
- FEM详细解析 — 实体单元含焊接形状模型
- 应力分类 — 喷嘴接续部的SCL位置特别重要。遵循ASME Annex 5-A
WRC 297/537与FEM结果一致吗？
大体一致，但常有10-20%差异。WRC基于无限圆筒壳上接喷嘴的理论解，当有端部效应或喷嘴间干涉时不精确。重要容器应用FEM验证。
镜板的设计
镜板（头部）有哪些种类？
| 镜板类型 | 形状 | 应力特性 | 成本 |
|---|---|---|---|
| 半球形 | 完全球壳 | 应力最小 | 最高 |
| 半椭圆形（2:1） | 椭圆旋转体 | 皱纹部有不连续应力 | 中等 |
| 碟形（三球面） | 球壳+圆锥切线 | 皱纹部高应力 | 相对便宜 |
| 平板 | 平面 | 应力最大（弯曲主导） | 最便宜 |
半椭圆形的"皱纹部不连续应力"是什么？
皱纹（knuckle）是镜板与胴的接续处，曲率急剧变化的地点。膜理论无法处理这种曲率变化，因此产生局部弯曲应力（不连续应力）。用轴对称FEM可精确评估。
2:1半椭圆镜板的皱纹部应力可达胴环向应力的1.5-2.0倍。这是Division 2中FEM被要求的典型情况。
实务检查清单
请给压力容器解析的检查清单。
应力分类的正确性最难判断。
是的。应力分类不能自动执行。工程师的判断必需。特别是"这个应力是一次还是二次"的判断需要经验，是压力容器工程师最重要的技能。
氢燃料汽车高压储罐解析
丰田MIRAI（2014年上市）搭载的70MPa高压氢储罐采用CFRP缠绕+铝内衬结构，其内压解析使用轴对称FEM与平面应力壳的复合模型。按ASME VIII Div.3基准设计破裂压力（175MPa）的FEM预测为176.2MPa，达到了实验破裂压力的99.3%。
压力容器线性软件比较
压力容器设计工具
压力容器设计用什么工具？
规格计算软件和FEM二阶段组成很一般。
规格计算软件（按规则设计）
| 软件 | 支持规格 | 特点 |
|---|---|---|
| PVElite | ASME VIII Div.1/2, EN 13445 | 使用最广。代码计算+简易FEM |
| Compress | ASME VIII Div.1/2 | ASME专长。代码检查自动化 |
| PASSA/Vessel | EN 13445 | 欧洲规格强 |
| AutoPIPE Vessel | ASME, EN | 与管道解析统合 |
PVElite是业界标准？
在全球是这样。胴、镜板、喷嘴、法兰、支撑脚的规格计算一贯进行。内置简易FEM（Nozzle Pro等）也可评估喷嘴不连续应力。
FEM（按分析设计）
用FEM解析压力容器选什么求解器？
Ansys的应力线性化功能对压力容器有利。
Ansys Workbench的Linearized Stress结果在压力容器工程师中很受欢迎。SCL视觉化设置，膜、弯曲、峰值各成分自动按ASME格式输出。Abaqus也可以，但需要更多手动操作。
特殊工具
选型指南
总结下？
规格计算软件覆盖大部分，FEM作为验证和复杂部位使用。
是的。压力容器设计中FEM是"辅助工具"，规格理解和应力分类技能最核心。工具再好，如果工程师应力分类错误，设计就成不了。
压力容器解析求解器的认证
原子炉压力容器等核压力容器的解析要求使用美国核管理委员会（NRC）认证的验证过的FEM代码。Abaqus（Dassault Systèmes）、ANSYS Mechanical、ADINA（Adina R&D）都取得了NRC NUREG/CR-7103（2012年）认证，在ASME VIII验证问题中精度99%以上，都有文件。Nastran虽包括于结构规范但无原子力专用认证。
压力容器线性先端研究
压力容器解析的先端课题
压力容器解析技术怎样发展？
有三个活跃方向。
弹塑性崩溃解析（DBA-Plastic）
听说不是弹性应力分类，而是直接塑性解析做设计的手法。
ASME Div. 2 Part 5.2.3 的弹塑性解析法。用弹完全塑性模型（屈服应力 = 1.5S）施加2倍设计荷载，解收敛则合格。不再需要应力分类的划时代手法。
不需要应力分类！这是大优点。
应力分类的模糊性（SCL位置、一次/二次判断）排除掉了。缺点是需要非线性FEM分析，计算成本比弹性解析大。但近年计算机性能提高，实用化进展。
适用性评估（Fitness-for-Service, FFS）
已有压力容器的健全性评估怎么做？
FEM应用：
- Level 3评估 — 详细FEM解析的健全性评估
- 有裂纹容器 — FEM计算J积分或SIF，破坏力学评估
- 腐蚀减厚 — 实测板厚分布映射到FEM模型进行应力评估
腐蚀使板厚不均的容器直接用FEM解析。
用实测数据（超声波测厚或激光扫描）映射FEM网格。板厚分布不一致的模型进行应力解析，最薄部应力在许可值以内确认。
氢气环境下的设计
氢储罐（Type I-IV）是压力容器最热的领域。
| 类型 | 结构 | 最高压力 | 用途 |
|---|---|---|---|
| Type I | 全金属 | ～200 bar | 工业用 |
| Type II | 金属内衬+部分FRP缠绕 | ～300 bar | 公交车 |
| Type III | 金属内衬+全面FRP缠绕 | ～700 bar | FCV |
| Type IV | 树脂内衬+全面FRP缠绕 | ～700 bar | FCV（最轻） |
Type IV几乎全部荷载由FRP（碳纤）负担。
总结
整理压力容器的先端课题。
压力容器看似"已成熟技术"，但氢社会和高龄装置的课题中产生了新解析需求。
应力线性化和疲劳评估
压力容器疲劳寿命按ASME VIII应力线性化（Stress Linearization）程序评估。断面上分解为膜应力、弯曲应力、峰值应力，各有许可值。1990年代ABB原子能的PWR压力容器喷嘴解析用Abaqus自动处理200条应力线性化路径的脚本，手动评估比工程工数削减90%的事例已发表在核工学期刊。
压力容器线性故障排除
压力容器解析的故障
压力容器FEM解析的常见故障？
压力容器特有的故障多与应力分类有关。
环向应力与薄肉公式不符
FEM的环向应力与 $pD/(2t)$ 不符。
检查项目：
1. 端部条件 — 开放端或闭合端，轴向应力不同，泊松效应影响环向应力
2. 不连续部影响范围内 — 镜板接续部距离 $\sqrt{Rt}$ 以内有不连续应力。离足够远处比较
3. 厚肉效应 — $D/t < 20$ 时需与拉梅式比较。内表面应力高于薄肉公式
4. 压力施加面 — 确认内表面压力向外（法向）
端部离足够远处比较是关键。
胴中央部（端离足够远的一致区域）与薄肉公式比较。不符则建模有问题。
应力分类困惑
"这个应力是一次还是二次"不知道。
最难的判断。基本原则：
一次应力（$P$） — 荷载平衡必需的应力。移除则不平衡。
- 例：内压膜应力、自重应力
二次应力（$Q$） — 变位协调条件产生的应力。具有自制性。
- 例：不连续应力（胴与镜板接续）、热应力
判断测试："移除此应力结构不崩溃？"→ 是则二次。否则一次。
不连续应力移除，结构虽变形但应力自重新配分，不崩溃，所以是二次。
SCL设置位置影响结果
SCL位置稍微变化，应力分类结果也变。
这是SCL法的根本性弱点。对策：
- 遵循ASME Annex 5-A指南 — SCL位置规定
- SCL板厚方向垂直 — 曲面上是局部法线方向
- 多个SCL位置确认结果 — 采最严格的
- 弹塑性解析法检验 — 不需SCL的手法验证
弹塑性解析法彻底解决SCL问题。
是的。如果为应力分类困扰，应该考虑弹塑性解析法（ASME Div. 2 Part 5.2.3）。计算成本增加，但判断模糊排除。
压力方向错误
施加了内压，应力反而是压缩。
压力方向反向了。FEM中压力沿面法向作用。法向内向时，正压作为外压。
检查：
- 显示单元法向确认
- 内压应变位向外（半径增大）
- 反力与压力×面积一致
实务检查清单补充
压力容器FEM解析的最终检查清单。
"薄肉公式整合确认"是出发点。不符不能推进。
完全同意。简单理论解对齐是FEM解析的铁则。压力容器中薄肉公式是最基本的验证手段。
喷嘴孔口应力集中误差
压力容器喷嘴取合部（贯通孔）应力集中系数Kt=2.5～4，但粗壳网格严重低估集中应力。ASME PTC 60指南推荐孔口半径r对网格尺寸 ≤ r/5，且用二阶单元。如不守此基准用3倍粗网格，峰值应力预测可低约40%的实例存在。
相关课题
更详细
报告