练练手的，和大家交流一下：
THICK=0.1
RADIUS=0.5
LENGTH=1.5
/PREP7
ET,1,182,1,,2
ET,2,169
ET,3,171
R,2,0,0,1,0.1
RMORE,0.1,0.01
R,3,0,0,1,0.1
RMORE,0.1,0.01
!* !* Material Properties !*
MP,EX,1,16E6
MP,NUXY,1,0.33
MP,MU,1,0.3
TB,MISO,1,0,9
TBPT,,0.000625,10000
TBPT,,0.0025,15000
TBPT,,0.005,21000
TBPT,,0.01,29000
TBPT,,0.015,32600
TBPT,,0.02,34700
TBPT,,0.04,36250
TBPT,,0.10,39000
TBPT,,0.20,40250
K,1
K,2,LENGTH
KGEN,2,1,2,1,,THICK,,2
A,1,2,4,3
TYPE,1
MAT,1
LESIZE,1,,,70
LESIZE,2,,,5
LESIZE,3,,,70
LESIZE,4,,,5
AMESH,ALL
NUMMRG,ALL
ASEL,NONE
LSEL,NONE
WPAVE,,THICK+RADIUS
CSWPLA,11,1
K,20,RADIUS,180
K,21,RADIUS,-90
K,22,RADIUS,0
K,23,0,0
k,24,radius,90
l,20,21
L,21,22
l,22,24
l,24,20
TYPE,2
REAL,2
LMESH,ALL
ALLS
CSYS,0
NSEL,S,LOC,Y,THICK
TYPE,3
ESURF
nsel,s,loc,y,0
real,3
esurf
ASEL,NONE
LSEL,NONE
WPAVE,,-THICK
CSWPLA,12,1
k,30
k,31,thick,0
k,32,thick,90
k,33,thick,180
l,31,32
l,32,33
type,2
lmesh,all
KMESH,23
KSEL,S,KP,,23
NSLK,S
*GET,N_LOAD,NODE,,NUM,MAX
csys,0
!* !* Boundary Conditions !*
asel,all
esla,s
nsle,s
nsel,r,loc,x,0
d,all,ux,0
nsel,r,loc,y,thick
d,all,uy,0
alls
FINISH
/SOLU
ANTYPE,STATIC
NLGEOM,ON
SOLC,ON
NSUBST,20,500,10
OUTRES,ALL,ALL
D,N_LOAD, ,2, , , ,ROTZ, , , , ,
SOLVE
FINI
大辊可变形，小辊刚性的辊压练习：
fini
/clear,start
THICK=0.1
RADIUS=0.5
LENGTH=1.5
/PREP7
ET,1,182,1,,2
ET,2,169
ET,3,171
R,2,0,0,1,0.1
RMORE,0.1,0.01
R,3,0,0,1,0.1
RMORE,0.1,0.01
!* !* Material Properties !*
MP,EX,1,16E6
MP,NUXY,1,0.33
MP,MU,1,0.3
TB,MISO,1,0,9
TBPT,,0.000625,10000
TBPT,,0.0025,15000
TBPT,,0.005,21000
TBPT,,0.01,29000
TBPT,,0.015,32600
TBPT,,0.02,34700
TBPT,,0.04,36250
TBPT,,0.10,39000
TBPT,,0.20,40250
MP,EX,2,210E8
MP,NUXY,2,0.27
MP,MU,2,0.3
K,1
K,2,LENGTH
KGEN,2,1,2,1,,THICK,,2
A,1,2,4,3
TYPE,1
MAT,1
LESIZE,1,,,70
LESIZE,2,,,5
LESIZE,3,,,70
LESIZE,4,,,5
AMESH,ALL
asel,none
WPAVE,,THICK+RADIUS
cyl4,,,radius,0,radius/3,90
cyl4,,,radius,0,radius/3,-90
mat,2
esize,0.04
mshkey,1
amesh,all
ARSYM,X,ALL
esla,s
nsle,s
nummerg,node
CSWPLA,11,1
k,30
ASEL,NONE
LSEL,NONE
CSYS,0
WPAVE,,-THICK
CSWPLA,12,1
k,40,thick,0
k,41,thick,90
k,42,thick,180
k,43,thick,-90
l,40,41
l,41,42
l,42,43
l,43,40
type,2
real,2
lmesh,all
kmesh,30
KSEL,S,KP,,30
NSLK,S
*GET,N_LOAD,NODE,,NUM,MAX
asel,s,loc,x,thick/2,thick
esla,s
nsle,s
nsel,r,loc,x,thick/2
type,3
esurf
nsle,s
nsel,r,loc,x,thick
type,2
real,3
esurf
csys,11
asel,s,loc,x,0,radius
esla,s
nsle,s
nsel,r,loc,x,radius
REAL,3
esurf
csys,0
asel,s,loc,y,0,thick
esla,s
nsle,s
nsel,r,loc,y,0
type,3
real,2
esurf
real,3
esla,s
nsle,s
nsel,r,loc,y,thick
esurf
nsle,s
nsel,r,loc,x,0
d,all,ux,0
nsel,r,loc,y,thick
d,all,uy,0
csys,11
nsel,s,loc,x,radius/3
d,all,all
csys,0
alls
FINISH
/SOLU
ANTYPE,STATIC
NLGEOM,ON
SOLC,ON
NSUBST,20,500,10
OUTRES,ALL,ALL
!,N_LOAD, uy,0.01
D,N_LOAD, ,2, , , ,ROTZ, , , , ,
SOLVE
FINI
上一篇： 在ANSYS中实现预应力
下一篇：ansys命令中英文解释
