工程上普遍采用的薄壁压力容器几何形状主要由轴对称壳体如圆柱形、球形、圆锥形、椭圆形和蝶形等壳体组成，薄壁容器是指壳体外径与内径之比≤1.2（Do/Di≤1.2）的容器。

本次以薄壁球壳为例，采用无力矩理论，进行ANSYS应力分析。

一、问题描述

薄壁球壳，球体直径D =1200mm，壁厚t =10mm。内压p =2MPa。弹性模量E =200GPa，泊松比μ =0.3。计算薄壁球壳上的应力。

问题分析：

（1）因为是回转薄壳，选取任意圆周角度。本算例选取1/4球壳建模，角度为90°，在球坐标系下180°~270°建模。

（2）采用SHELL41薄膜单元，也叫做膜单元。该单元仅有面内膜刚度而无面外弯曲刚度，该单元的每个节点有3个自由度，即沿节点坐标系x、y和z方向的平动位移。

（3）坐标系

（4）将节点坐标系转换到球坐标系[CSYS,2]。

耦合所有节点的径向UX自由度。

模型底端线的节点约束经向UZ自由度。

两侧边线上的节点约束纬向UY自由度。

尖角处的节点同时约束UZ和UY自由度。

（4）在所有单元上施加内压2MPa。根据SHELL41薄膜单元示意图，内压Loadkey为①，在单元施加压力时，Loadkey选择1。

（5）后处理时，将结果坐标系转换到结果球坐标系，结果球坐标系X为径向、Y为纬向，Z为经向。

计算结果：ANSYS数值解结果纬向应力（地球经纬的纬）SY为59.95MPa，经向应力SZ（地球经纬的经）为59.97MPa。解析解结果纬向和经向应力都为60MPa。ANSYS数值解和解析解结果吻合很好。

二、解析解 

主要参考大学本科过程装备与控制工程的《过程设备设计》教材——压力容器应力分析>回转薄壳应力分析>薄壁圆筒的应力。

参考文献：   

[1] 郑津洋，桑芝富. 过程设备设计（第4版）[M]. 北京: 化学工业出版社, 2015.

[2] 王心明，W.Z.麦克. 工程压力容器设计与计算（第2版）[M]. 北京: 国防工业出版社, 2011.

文章来源：申江储气罐  http://www.jhztgz.cn