辊压接触分析

　　本文探讨了辊压接触分析的相关内容。

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

　　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