ex9.3.inp

ansys有限元分析程序

! 9.3  利用8节点2D实体单元PLANE82计算压力容器
! 本程序来源于邢静忠编著的《ANSYS应用实例与分析》，科学出版社，2006年
!
Finish
/Clear
/FileName, EX9.3
/Title, EX9.3, Long Cyclinder Pressure Vessel for CNG with PLANE82.
! （1）定义基本参数
h1=298.5 		! 封头的高度坐标298.5 mm
b1=44.5  		! 封头的内半径44.5 mm
b2=b1+23 		! 封头的外半径44.5+23mm厚度
r1=185   		! 罐体内半径185mm
r2=203   		! 罐体外半径203mm
! （2）定义剖面区域边界上的关键点
/Prep7   		! 进入前处理模块
k,1,b1,h1		! 封头起始端内侧的点
k,2,b2,h1		! 封头起始端外侧的点
k,3,b1,sqrt(r1*r1-b1*b1) 	! 封头内侧剖面线和罐体半球内剖面线的交点
k,4,b2,sqrt(r2*r2-b2*b2) 	! 封头外侧剖面线和罐体半球外剖面线的交点
k,5,r1,0    		! 罐体半球和罐体主体的内侧交点
k,6,r2,0    		! 罐体半球和罐体主体的外侧交点
k,7,0,0     		! 坐标原点
k,8,r1,-1000+h1 	! 罐体主体截断内侧点
k,9,r2,-1000+h1 	! 罐体主体截断外侧点
K,10,0,h1   		! 封头起始端对称线上的点
! （3）定义围成结构剖面区域的线，并定义剖面
L,1,2      		! 封头起始端水平线
L,1,3      		! 封头内侧竖直线
L,2,4      		! 封头外侧竖直线
L,5,8      		! 罐体主体内侧线
L,8,9      		! 罐体主体截断位置水平线
L,6,9      		! 罐体主体外侧线
Larc,5,3,7,r1	! 过5号，3号关键点，以7号关键点为圆心，r1为半径做半球内圆弧线
Larc,6,4,7,r2	! 过6号，4号关键点，以7号关键点为圆心，r2为半径做半球外圆弧线
LFILLT,2,7,63, , 	! 过2号，7号线，以半径63mm做半球内侧位置的倒圆
LFILLT,3,8,40, , 	! 过3号，8号线，以半径40mm做半球外侧位置的倒圆
AL,all       		! 将所有线围成的区域定义为一个面
! （4）定义单元类型、材料参数和实常数
ANTYPE,STATIC   	! 定义分析类型为静力分析
ET,1,PLANE82		! 定义第1类单元是8节点二维实体单元(用轴对称选项)
KEYOPT,1,3,1		! 定义单元的第3个可选项，打开轴对称选项
MP,EX,1,207E3   	! 定义第1类材料的弹性模量为207 GPa
MP,NUXY,1,0.3   	! 定义第1类材料的泊松比NUXY=0.3
! （5）定义位移约束和荷载 
DL,5,,SYMM  		! 给5号线定义对称位移条件
SFL,2,PRES,23, , , , 	! 给剖面线的内侧线2施加内压荷载23MPa
SFL,4,PRES,23, , , , 	! 给剖面线的内侧线4施加内压荷载23MPa
SFL,7,PRES,23, , , , 	! 给剖面线的内侧线7施加内压荷载23MPa
SFL,9,PRES,23, , , , 	! 给剖面线的内侧线9施加内压荷载23MPa
SFL,1,PRES,-17.681, 	! 给端部封头水平线施加封头传递过来的-17.681MPa
Esize,5         		! 指定单元尺寸大小为5mm
Amesh,1         		! 将罐体剖面1号面剖分为面单元
FINISH          		! 退出前处理模块
! （6）进入求解模块求解
/SOLU      		! 进入求解模块
Solve    		! 开始求解
FINISH   		! 退出求解模块
! （7）进入后处理模块处理计算结果
/POST1  		! 进入后处理模块
Pldisp,2		! 显示变形图
Plnsol,s,1  	! 显示最大拉应力分布
Plnsol,s,2  	! 显示第2主应力分布
Plnsol,s,3  	! 显示最大压应力分布
Plnsol,s,EQV	! 显示最大Mises应力分布
Plnsol,s,Z   	! 显示环向应力分布
Finish