压力容器开孔接管区的应力分析.doc

压力容器开孔接管区的应力分析 finish /clear /filn, E322 /title, FEA of connecting zone of nozzle to cylinder ! *************参数设定******************* Rci=1000 !筒体内半径 tc=30 !筒体厚度 Rco=Rci+tc !筒体外半径 Lc=4000 !筒体长度 Rno=530 !接管外半径 tn=15 !接管厚度 Rni=Rno-tn !接管内半径 Hmin=sqrt(Rci**2-Rni**2) !接管最,,,,内伸长度 Lni=50 !接管内伸长度 Ln=500 !接管外伸长度 rr1=30 !焊缝外侧过渡圆角半径 rr2=15 !焊缝内侧过渡圆角半径 pi=1.2 !内压 pc=pi*Rci**2/(Rco**2-Rci**2) !筒体端部轴向平衡面载荷 pn=pi*Rni**2/(Rno**2-Rni**2) !接管端部纵向平衡面载荷 !****************前处理*************************** /prep7 et,1,45 !定义单元类型 mp,ex,1,2e5 !定义材料的弹性模量 mp,nuxy,1,0.3 !定义材料的泊松比 !****************建立模型*************************** cylind,Rco,Rci,0,-Lc/2,90,270, !生成筒体 wpoff,0,0,-Lc/2 !将工作面沿-Z向移动Lc/2 wprot,0,90, !将工作面沿yz旋转90o cylind,Rno,Rni,-Ln-Rci-tc,-Hmin+Lni,90,180, !生成接管 vovlap,all !体overlap布尔运算 vdel,3,,,1 !删除筒体上开孔失去的部分 vsel,s,,,7 !选择筒体 *afun,deg !设定角度函数中单位为角度 ang1=2*nint(asin(Rno/Rci)) !计算接管区切割角度 wprot,0,0,-90+ang1 !旋转坐标系 vsbw,all !切割筒体 afillt,21,12,rr1 !筒体与接管外表面圆角 al,33,60,5 !生成圆角区域所需的面 al,37,61,7 !生成圆角区域所需的面 va,12,14,16,34,35 !生成圆角区域体 afillt,23,13,rr2 !筒体内表面与接管外表面圆角 al,24,40,68 !生成圆角区域所需的面 al,25,41,69 !生成圆角区域所需的面 va,13,23,31,39,40 !生成圆角区域体 alls,all !全选 askin,62,70 !根据内外圆角边界蒙皮生成切割面 vsba,2,36 !切割筒体 wprot,0,,90-ang1 !旋转坐标系 wprot,0,90,0 !旋转坐标系 wpoff,0,0,-2*Rni !移动坐标系至接管区轴向切割位置 vsbw,9 !切割筒体 vsbw,1 !切割筒体 !*************划分网格****************