压力容器设计 — CAE术语解说

很好的问题。膜应力是沿容器壁厚方向均匀分布，主要承受内压的拉伸应力。而弯曲应力发生在喷嘴、法兰等不连续部位，沿壁厚方向符号相反。从危险性角度讲，弯曲应力更容易成为脆性破裂的起点。ASME Boiler and Pressure Vessel Code, Section VIII, Division 2 规定，膜应力许容值为屈服强度的 2/3，弯曲应力为屈服强度的 1.0 倍。例如，SM490 钢（屈服强度 325 MPa）的膜应力许可值约为 217 MPa，弯曲应力许可值为 325 MPa。

实务中通常以内半径 r 与厚度 t 的比值 r/t > 10 作为薄壁的标准。例如内径 1000 mm、厚度 20 mm，则 r/t=25，属于薄壁。对于厚壁情况，内壁和外壁应力不同，需用 Lamé 公式。内压 p 作用下厚壁圆筒内壁切向应力为

典型产生部位有：①壳体与端盖的接合部，②喷嘴安装部，③法兰周围，④鞍座支持部。CAE 中这些局部区域需细化网格分划。例如喷嘴周围，在壳体厚度的 1.5 倍范围内，网格尺寸应小于厚度的 1/4。Ansys Workbench 的"Sphere of Influence"功能可提高效率。不连续应力衰减快速（"edge effect"），仅需局部评估，但应力集中系数（SCF）常超过 3.0。

判断标准如下：
壳单元：薄壁结构（r/t > 10）、整体应力评估、初期设计阶段。计算成本低，膜应力/弯曲应力分量易于分解。典型产品有 Abaqus 的 S4R 和 Ansys 的 Shell181。
实体单元：厚壁结构、不连续部详细评估、应力集中、疲劳分析。能捕捉喷嘴圆角部等三维应力状态。建议使用二阶单元（20 节点六面体）。但网格数常增加 10 倍以上。

受弯曲的部位，厚度方向至少需 3 层，理想情况 4～6 层二阶单元。一层线性实体单元会产生"剪切锁定"问题，刚性被过估计。例如 30 mm 厚平板弯曲，一层模型会低估挠度 50% 以上。ASME Code Case 2695 等也暗示，为执行沿应力分类线（Stress Classification Line）的应力线性化，需要厚度方向有足够的节点数。

接触必要部位：①法兰与垫圈界面，②多层结构层间，③鞍座与壳体接触。用简单约束（CP 或 Coupling）表示"螺栓联接"会使法兰旋转刚性过估计，螺栓轴力和垫圈面压易被低估，存在泄漏评估风险。真实建模螺栓，设置初始预紧力（例如按 JIS B 2251 的 70% 为参考）和法兰面摩擦接触（摩擦系数 μ=0.15 左右），可得到更现实的评估。Ansys 的"Bolt Pretension"功能或 Abaqus 的"Bolt Load"可用。

至少要确认以下内容：
1. 设计压力、温度：不只是常规运转条件，还要明确设计条件（通常为 1.1 倍）和水压试验条件（设计压力的 1.3～1.5 倍）。
2. 材料特性：设计温度下的屈服强度、拉伸强度、弹性模量、泊松比。例如 JIS G 3106 SM490，室温屈服强度 325 MPa，但 150℃ 时约降至 310 MPa。
3. 荷载工况：内压、自重、内装物重量、风荷载、地震荷载的组合。ASME VIII Div.2 中定义为"Load Case"。
4. 模型范围：能否利用对称性，需要全模型还是 1/2、1/4 模型。

判断结果"可信性"的三步法：
首先是1. 网格敏感性分析：特别在应力集中部，将网格细化 1.5 倍，最大应力变化在 5% 以内即可收敛，超过则未收敛。
其次是2. 全体行为确认：变形图物理上是否合理（内压下膨胀等）。反力合计是否与内压×投影面积平衡（例如内压 1 MPa、内径 1000 mm 圆筒，端盖应受约 785 kN 轴力）。
最后是3. 单纯部位验算：圆筒部膜应力手算（

ASME VIII Div.2, Annex 5-A 有粗略指针，但具体位置由工程师判断。不能完全依赖软件。关键是"不连续部充分远离、仅一次应力主导的区域"和"评估对象结构的全厚度"相交的直线。例如喷嘴接合部，在壳体侧应至少离不连续点

专业度从高到低排序：Ansys Mechanical + ASME 规范评估模块功能最完善。特别是"ACP（ASME Code Postprocessor）"可自动化执行 ASME VIII Div.2 符合的应力线性化、分类、与许可值对比。次之是Abaqus/CAE。接触算法强大，Python 脚本自动化能力强，特别适合复杂法兰紧固分析。而COMSOL Multiphysics在热流体与结构耦合（热应力）及腐蚀等化学现象与结构结合的"多物理场"分析中强大，但原生 ASME 规范评估功能较弱。

有的。PV Elite 和 COMPRESS 等压力容器专用设计软件擅长 ASME VIII Div.1 面积补强法等简易计算。实务做法是：用这类软件确定全体布局和主要尺寸，不连续部和复杂形状再用 Ansys 或 Abaqus 做详细 FEA。意义在于：FEA 易遗漏"所有喷嘴"的网格检查，简易法可全面覆盖。例如容器有 20 个喷嘴，FEA 仅对最严酷的 2～3 个，其余用简易法验证。

CalculiX（PrePoMax 前端）或 Code_Aster（Salome-Meca 环境）的线性静力分析核心精度与商用软件相当。弹性范围应力分布精度足够。但致命差异在"规范评估自动化和验证"。商用软件的 ASME 模块经数十年实务应用和规范委员会对话验证。免费软件要达同等信度，需用户自己用脚本实现应力线性化算法到许可应力适用的全过程，验证工数和风险很大。学习和研究可用，但设计认证使用现阶段困难。

非线性分析（接触或大变形）中常见原因三点：
1. 接触设置问题：初始接触状态"分离过远"或"嵌入过深"。特别是垫圈模型要准确设置初始间隙。用"Adjust to Touch"功能或手动调整。
2. 荷载陡峭加载：内压或螺栓预紧力一次满载会不收敛。应分 10～20 步加载，减小增量步（Ansys 中开启"Auto Time Stepping"）。
3. 刚性矩阵奇异性：存在未约束刚体运动。如鞍座支持只约束 Z 向，但 X 向热膨胀未约束等，需确保边界条件物理正确。

尖锐凹角（缺口）的弹性分析中，应力在理论上趋于无穷（"奇异点"）。这是正确的数学现象，但实际材料会用塑性变形缓解应力。CAE 对策有二：
1. 设计改进：反映实际加工，必须加圆角（例如 R≥5 mm）。分析模型中也加相同圆角。
2. 评估方法改变：不评估奇异点本身应力，而评估略远处（例如圆角根部外厚度 1/4）的应力。或用"结构应力法"、"Boltzmann 应力"等避开奇异点影响的方法。ASME VIII Div.2 规定，局部应力集中除疲劳评估外可忽略。

仅壳单元难以准确评估喷嘴接合部。标准工作流是：
1. 全体模型用壳单元，把握大致应力分布和荷载路径。
2. 采用子模型（切边界）技术。从全体模型中切出要评估的喷嘴周围部分，将全体模型得到的位移作为强制位移加到切边界。
3. 子模型用实体单元详细网格分划，含圆角形状。
这样以计算成本低且局部三维应力精度高。Ansys 的"Submodeling"、Abaqus 的"Submodel"功能对应。