《压力容器分析设计-直接法》L.4.1命令流

!加氢裂化反应器的总体塑性变形设计校核GPD-DC,附录L.4.1!!!!!!!!环境设置!!!!!!!!!
FINISH
/CLEAR
/filnam,model,1
/PREP7
pi=acos(-1)
*afun,deg

!geomety cylindrical shell
riz=2515/2
raz=2515/2+190
hz=2000
k,1,riz,0,0
k,2,raz,0,0
k,3,raz,hz,0
k,4,riz,hz,0

!geometry hemispherical end
rik=1297.5
rak=1297.5+110
local,11,1,,hz
k,5,rak,acos(riz/rik),0
k,6,rik,acos(riz/rik),0

!geometry nozzle
ris=520/2
ras=520/2+130
rar=520/2+90
hsv=295
hsr=150
rf=50
k,7,rak,acos(ras/rak),0
k,8,rik,acos(ris/rik),0
L,5,7
L,8,6
CSYS,0
K,9,RAS,hz+(sqrt(rak**2-ras**2))+hsv
K,10,RiS,hz+(sqrt(rak**2-ras**2))+hsv
K,11,RAr,hz+(sqrt(rak**2-ras**2))+hsv+(ras-rar)/tan(30)
K,12,Ris,hz+(sqrt(rak**2-ras**2))+hsv+(ras-rar)/tan(30)
K,13,Rar,hz+(sqrt(rak**2-ras**2))+hsv+(ras-rar)/tan(30)+hsr
K,14,Ris,hz+(sqrt(rak**2-ras**2))+hsv+(ras-rar)/tan(30)+hsr
L,7,9
L,10,8
LFILLT,1,3,RF,0
LFILLT,2,4,RF,0

!geometry flange
hf=200
raf=1045/2
rda=700/2
K,19,Raf,hz+(sqrt(rak**2-ras**2))+hsv+(ras-rar)/tan(30)+hsr
K,20,Raf,hz+(sqrt(rak**2-ras**2))+hsv+(ras-rar)/tan(30)+hsr+hf

K,21,Rda,hz+(sqrt(rak**2-ras**2))+hsv+(ras-rar)/tan(30)+hsr+hf
K,22,Ris,hz+(sqrt(rak**2-ras**2))+hsv+(ras-rar)/tan(30)+hsr+hf

a,1,2,3,4
a,4,3,5,6
a,10,9,11,12
a,12,11,13,14
lang,2,15,90,,0.1
lang,4,16,90,,0.7
al,1,22,21,12
al,5,24,23,6,2,22
al,3,14,4,24
a,14,13,21,22
a,13,19,20,21

csys,0
k,100,0,0,0
k,101,0,1,0
vrotat,all,,,,,,100,101,180,2

!defining elements
ET,1,plane42
ET,2,SOLID45

!material properties
!material 1 cylindrical shell & end
MP,EX,1,191915    !e-modul 305
MP,NUXY,1,0.3
TBDE,BKIN,1,,,
!TBTEMP,400
tb,bkin,1,2,2,1
TBDATA,,135,0

!material 2 nozzle
MP,EX,2,191915.001    !e-modul 305
MP,NUXY,2,0.3
TBDE,BKIN,2,,,
!TBTEMP,400
tb,bkin,2,2,2,1
TBDATA,,135.001,0

!material assignment
VSEL,S,VOLU,,1,2,1
VSEL,A,VOLU,,5,14,9
VSEL,A,VOLU,,10,11,1
VATT,1,,2,0
ALLS
VSEL,U,VOLU,,1,2,1
VSEL,U,VOLU,,5,14,9
VSEL,U,VOLU,,10,11,1
VATT,2,,2,0
ALLS

!meshing
lesize,5,,,11
lesize,2,,,4
lesize,23,,,4
lesize,6,,,3
lesize,7,,,5
lesize,9,,,5
lesize,12,,,5
lesize,14,,,5
lesize,16,,,5
lesize,19,,,5
lesize,25,,,5
lesize,27,,,5
lesize,28,,,5
lesize,22,,,5
lesize,24,,,5
lesize,26,,,5
lesize,29,,,5
lesize,30,,,5
lesize,3,,,10
lesize,4,,,10
lesize,15,,,5
lesize,17,,,5
lesize,18,,,5
lesize,20,,,5
lesize,1,,,25
lesize,21,,,25
lesize,11,,,8
lesize,13,,,8
lesize,8,,,30
lesize,10,,,30

lsel,s,line,,2,6,4
lsel,a,line,,23
lccat,all
alls

asel,s,area,,1,9,1
amesh,all
alls

extopt,aclear,1
extopt,vswe,auto,0
extopt,esize,22,0

vsweep,1,1,14
vsweep,10,14,53
vsweep,2,2,18
vsweep,11,18,57
vsweep,5,5,31
vsweep,14,31,70
vsweep,6,6,37
vsweep,15,37,76
vsweep,7,7,40
vsweep,16,40,79
vsweep,3,3,23
vsweep,12,23,62
vsweep,4,4,27
vsweep,13,27,66
vsweep,8,8,44
vsweep,17,44,83
vsweep,9,9,48
vsweep,18,48,87

!defining beam element
et,3,beam189
save
fini



/filnam,beamsection,1
/CLEAR
/PREP7
local,11,1,,,,180
*set,ri,520/2
*set,ra,1045/2
k,1,ri,0
k,2,ra,0
k,3,ra,180
k,4,ri,180
a,1,2,3,4

et,1,mesh200,7
esize,ra-ri
mshkey,1
amesh,all
secwrite,beamsection,sect,,1
save
fini

/filnam,gpd_mises,1
resume,MODEL,db,D:\YINGLIFENXI\     !见注
/prep7
sectype,1,beam,mesh,hrohr0
secRead,beamsection,,,mesh
secoffset,origin
K,200,0,hz+(sqrt(rak**2-ras**2))+hsv+(ras-rar)/tan(30)+hsr+hf,0
K,201,0,hz+(sqrt(rak**2-ras**2))+hsv+(ras-rar)/tan(30)+hsr+hf+200,0
L,200,201
LSEL,S,,,136
LATT,1,1,3
LESIZE,136,,,1
LMESH,136
ALLS

!boundary condition
csys,0
nrotat,all
nsel,s,loc,y,0
d,all,uy,0
nsel,s,loc,z,0
d,all,uz,0
nsel,s,loc,y,0
nsel,r,loc,x,0
d,all,ux,0

*set,h,hz+(sqrt(rak**2-ras**2))+hsv+(ras-rar)/tan(30)+hsr+hf
nsel,s,loc,x,0
nsel,r,loc,z,0
nsel,r,loc,y,h
d,all,rotx,0
d,all,roty,0

!rigid region
csys,0
*set,h,hz+(sqrt(rak**2-ras**2))+hsv+(ras-rar)/tan(30)+hsr+hf
nsel,s,loc,y,h,h+1
nsel,r,loc,x,0
nsel,r,loc,z,0
*get,ncerig,node,,num,max
nsel,s,loc,y,h,h+1
csys,5
nsel,r,loc,x,0,rda
csys,0
cerig,ncerig,all,uy
nsel,s,loc,y,h,h+1
nsel,r,loc,x,0
csys,5
nsel,r,loc,x,0,rda
csys,0
cerig,ncerig,all,ux
alls
save
fini

/SOLU
ANTYPE,0,NEW
ERESX,NO
PR=21.6
MS=313600000
FS=115200

ASEL,S,,,13,17,4
ASEL,A,,,22,26,4
ASEL,A,,,30,34,4
ASEL,A,,,35,36,1
ASEL,A,,,39,43,4
ASEL,A,,,52,56,4
ASEL,A,,,61,65,4
ASEL,A,,,69,73,4
ASEL,A,,,74,75,1
ASEL,A,,,78,82,4
NSLA,S,1
SF,ALL,PRES,PR
ALLS

F,NCERIG,FY,PR*RIS**2*PI/2-FS/2
F,NCERIG,MZ,-MS/2

TIME,1
DELTIM,0.05,0.0001,0.05
NEQIT,15
OUTRES,ALL,ALL
SOLV
fini

/post26

ANSOL,2,19362,EPEL,1,EPEL1_2
STORE,MERGE

!注：resume,MODEL,db,D:\YINGLIFENXI\ ，要改成你自己的ansys工作目录

牛逼的人生不需要解释！！吊炸天了

大师稍微介绍一下啊。

就是把《压力容器分析设计-直接法》中的一个分析实例敲了下来。是一个筒体加球封头加接管的一个结构，根据en13445对其进行总体塑性变形设计校核。

学习一下