压力容器中ansys分析支座反力的命令流.doc

FINISH/CLEAR/FILN, E43 /TITLE, FEA OF CONNECTING ZONE OF SUPPORT OF EVAPERATOR! * 参数设定*DCI=5000 !筒体内径TC=40 !筒体壁厚TCON=40 !锥形封头壁厚ALFA=30 !锥形封头半锥角LC=2500 !筒体长度*AFUN,DEG !角度单位设定LCON=DCI/2/TAN(ALFA) !锥形封头全锥长度LRAT=0.5 !模型中锥形封头经线方向长度系数LCON=LCON*LRAT !模型中锥形封头经线方向长度NS=4 !支座个数H0=1500 !支座垂直方向安装位置DIS=DCI+TC !支座水平方向安装位置TS=20 !支座筋板、腹板厚度WS=500 !筋板宽度HS=1300 !筋板高度BB=800 !底板长度W1=250 !腹板位置PI=0.1 !内压值GRA=9.81 !重力加速度，M/S2DENSTELL=7850 !钢材密度，KG/M3DENWATER=1000 !水的密度，KG/M3M1=58.1E3 !容器自重，KGHY1MAX=13 !锥形封头部分平均液面高度HY0MAX=15.0 !最大液面高度PY1=HY1MAX*DENWATER*GRA/1E6 !锥形封头部分平均液柱静压力PY0=HY0MAX*DENWATER*GRA/1E6 !容器最大液柱静压力SCON=(DCI/2*LRAT*2+TCON*COS(ALFA)*TCON!模型中锥形封头端面面积/?PCON=(PI+PY0)*(DCI/2*LRAT)*2/SCON/COS(ALFA) !模型中锥形封头端面平衡面载荷PCM1=M1*GRA/(3.1415927*(DCI+2*TCON)*2-DCI*2) !自重在筒体端面引起的面载荷PC=PI*DCI*2/(DCI+2*TCON)*2-DCI*2)-PCM1 !筒体端面平衡面载荷!*前处理*/PREP7ET,1,45 !定义单元类型MP,EX,1,2E5 !定义材料的弹性模量MP,NUXY,1,0.3 !定义材料的泊松比!*建立模型*K,1,0,0,0 !定义回转轴中心关键点K,2,0,LC,0 !定义回转轴中心关键点BLC4,DCI/2,0,TCON,-LCON !生成锥形封头子午面母体LOCAL,11,1,DCI/2 !定义局部坐标系AGEN,2,1,-ALFA,1 !旋转母体得到旋转锥形封头子午面BLC4,DCI/2,0,TCON,LC !生成筒体子午面NUMMRG,ALL, , , ,LOW !合并相同项KSEL,S,LOC,Y,-TCON,0 !选择筒体与封头焊接区关键点*GET,K1,KP,NUM,MIN !提取焊接区关键点KSEL,U,K1*GET,K2,KP,NUM,MIN*GET,K3,KP,NUM,MAXKSEL,A,K1A,K1,K2,K3 !生成焊接区子午面ALLS !全选VROTAT,ALL,1,2,360/2/NS !生成筒体、封头及焊缝区1/8回转体WPOFF,0,0,-BB/2 !移动工作面BLC4,-H0,DIS/2,HS,TS !生成筋板母体BLC4,DIS/2-W1,-H0,-TS,HS,BB/2 !生成腹板母体BLC4,-H0,DIS/2,TS,BB/2 !生成底板母体VSEL,S,4,6 !选择筋板、腹板及底板母体VSBA,ALL,5 !用锥形封头外表面切割筋板、腹板及底板VDELE,7,1 !删除筋板母体多余部分VDELE,10,1 !删除腹板母体多余部分VDELE,11,1 !删除底板母体多余部分WPOFF,DIS/2-WS !移动工作面WPROT,90 !旋转工作面VSBW,ALL !再次切割筋板与底板VDELE,4,1 !再次删除筋板多余部分VDELE,6,1 !再次删除底板多余部分VPTN,ALL !筋板、腹板及底板互分ALLS!*划分网格*KWPAVE,23 !将坐标平面移至关键点23WPROT,90 !旋转工作平面VSBW,ALL !采用工作平面剖分体KWPAVE,19 !将坐标平面移至关键点19VSBW,ALL !采用工作平面剖分体 WPROT,90, !旋转工作平面 VSBW,ALL !采用工作平面剖分体 VPTN,ALL !体互分ALLS !全选择LGLUE,ALL !粘合各线LESIZE,36,3, !设定线分段数LESIZE,55,3, !设定线分段数LESIZE,45,3, !设定线分段数ESIZE,40 !设定单元大小VSWEEP,ALL !映射划分NUMMRG,ALL !合并所有项NUMCMP,ALL !压缩编号FINISH!*求解*/SOLUWPCSYS,-1,0 !转换到原始坐标系ASEL,S,LOC,Z,0 !选择在Z=0的面DA,ALL,SYMM !施加对称约束LOCAL,12,1,-90 !定义局部柱坐标系ASEL,S,LOC,Y,360/NS/2 !选择1/8位置截面DA,ALL,SYMM !施加对称约束CSYS,0 !激活整体直角坐标系ASEL,S,LOC,Y,-H0 !选择支座底面ASEL,R,LOC,X,DIS/2-WS,DIS/2DA,ALL,ALL !约束底面各项位移ASEL,S,LOC,Y,LC !选择筒体端面SFA,ALL,1,PRES,-PC !施加平衡面载荷ASEL,S,LOC,Y,-LCON*COS(ALFA)-TCON,-LCON*COS(ALFA) !选择锥形封头端面SFA,ALL,1,PRES,-PCON !施加平衡面载荷ASEL,S,AREA,83,87,4 !选择筒体内表面ASEL,A,AREA,42,SFA,ALL,1,PRES,PI !在筒体内表面施加内压ASEL,S,AREA,44,46,2,1 !选择封头内表面