creo液体有限元分析

我做过一个储液罐的静水压力分析项目，客户用 Creo Simulate 想验证罐壁在装满液体时的强度和变形。我这个项目虽然简单，但把静水压力分布、壁面变形和 CAD-FEA 一体化流程讲得很清楚。

左侧图是储液罐的静水压力云图。罐内液体自由液面在顶部，压力从液面向下线性增加，颜色从顶部的浅黄到底部的深蓝，对应压力从 0 增加到约 45 kPa。静水压力的公式是 P = ρgh，其中 h 是深度。我假设液体密度 1000 kg/m³、罐高 4 m，底部压力就是 ρgH ≈ 39.2 kPa，加上顶部可能还有小气相压力，图中最大值 45 kPa 合理。

右侧图有两条曲线：蓝色实线是压力 P 随深度 h 线性增加，从 0 到 40 kPa；红色虚线是罐壁水平变形 u 随深度变化，顶部变形小，底部变形大，最大约 0.8 mm。变形的分布和弯矩分布有关：底部压力大，但罐壁底部通常有加劲或固定约束，所以变形最大值不一定在最底部，而可能在离底部一定高度处，具体取决于边界条件。

Creo Simulate 的操作流程：先在 Creo Parametric 里建好储液罐三维模型，定义材料（比如不锈钢，E=210 GPa，ν=0.3）；然后进入 Simulate，定义静水压力载荷，选择罐体内表面，设置液面高度和液体密度；底部约束固定或简支；运行静态分析，输出 von Mises 应力和变形。网格用自动划分，对薄壁结构可以抽中面做壳单元分析，效率更高。

踩坑经验：第一个是单位制：Creo 默认可能是英寸-磅，做静水压力时一定要注意统一成 m-Pa-kg 或 mm-MPa-tonne。第二个是压力方向：静水压力垂直于壁面指向外侧，如果方向反了，罐壁会被压瘪而不是膨胀。第三个是薄壁结构：如果罐壁很薄，大变形效应明显，需要打开几何非线性。

工程经验：储液罐设计通常以底部环向应力控制。我帮客户算了最大环向应力 σ_θ = ρgHR/t，其中 R 是半径、H 是液高、t 是壁厚。如果算出来应力接近材料屈服强度，就需要加筋或增加壁厚。这个案例里最大应力远小于许用应力，但底部变形 0.8 mm 对客户密封结构有影响，建议他在底部加环形加劲肋。

更多 Creo/Simulate 结构与流体静压力分析的实战案例，我整理在 [creo液体有限元分析](https://www.keyanxueshu.com/category/fea/) 栏目里。[creo液体有限元分析](https://www.keyanxueshu.com/)

薄壁储液罐的应力分析还要区分薄膜应力和弯曲应力。对于圆柱形罐壁，环向薄膜应力 σ_θ = ρgHR/t 是主要分量，轴向薄膜应力 σ_z = ρgH/(2t) 是次要分量；罐壁与底板连接处由于约束转折，会产生弯曲应力集中。我们这张图主要展示环向应力引起的膨胀变形，最大变形出现在中下部。如果客户关心焊缝或底板连接，需要用实体单元局部细化，看弯曲应力峰值。Creo Simulate 里可以用壳单元做整体分析，再用子模型技术把关键连接区切出来做实体分析，兼顾效率和精度。另一个设计要点是加劲肋间距，通常取 2-3 倍罐壁厚度对应的名义环向应力波长。