本文探討了輥壓接觸分析的相關內容。

　　練練手的，和大家交流一下：

　　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