練練手的，和大家交流一下：
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