资源描述：
《压力容器开孔接管区的应力分析》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、压力容器开孔接管区的应力分析finish/clear/filn,E322/title,FEAofconnectingzoneofnozzletocylinder!*************参数设定*******************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
2、=30!焊缝外侧过渡圆角半径rr2=15!焊缝内侧过渡圆角半径pi=1.2!内压pc=pi*Rci**2/(Rco**2-Rci**2)!筒体端部轴向平衡面载荷pn=pi*Rni**2/(Rno**2-Rni**2)!接管端部纵向平衡面载荷!****************前处理***************************/prep7et,1,45!定义单元类型mp,ex,1,2e5!定义材料的弹性模量mp,nuxy,1,0.3!定义材料的泊松比!****************建立模型*******************
3、********cylind,Rco,Rci,0,-Lc/2,90,270,!生成筒体wpoff,0,0,-Lc/2!将工作面沿-Z向移动Lc/2wprot,0,90,!将工作面沿yz旋转90ocylind,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,
4、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,3
5、6!切割筒体wprot,0,,90-ang1!旋转坐标系wprot,0,90,0!旋转坐标系wpoff,0,0,-2*Rni!移动坐标系至接管区轴向切割位置vsbw,9!切割筒体vsbw,1!切割筒体!*************划分网格******************************lsel,s,loc,y,Ln+rci+tc!选择接管端部各线段lsel,u,length,,0,tc!去掉长度为0至tc的线段lesize,all,,,20,,1!设定所选线段的剖分数alls,alllesize,60,,,4,,1!设定外圆
6、角某一线段的剖分数lesize,45,,,12,,1!设定内接管分段数lesize,57,,,12,,1asel,s,,,45,46,1,1!选择筒体上接管连接影响区外侧两截面accat,all!连接相邻面域alls,all!全选vsweep,all!映射划分allsnummrg,all!合并所有相同项numcmp,all!压缩fini!****************求解***************************/soluasel,s,loc,z,0!选择筒体端面SFA,all,1,PRES,-pc!施加端面平衡面载荷a
7、sel,s,loc,x,0!选择对称面asel,a,loc,z,-Lc/2!选择对称面DA,all,SYMM!施加对称约束asel,s,loc,y,Ln+Rci+Tc!选择接管端面DA,all,UY,!约束轴向位移asel,s,,,6!选择所有内表面asel,a,,,7asel,a,,,15asel,a,,,17asel,a,,,19asel,a,,,23asel,a,,,29asel,a,,,42asel,a,,,47asel,a,,,50SFA,all,1,PRES,pi!施加内压allssolve!求解fini!********
8、********后处理***************************/post1PLNSOL,S,INT,0,1!显示应力云图fini
