近几天在学习摩擦生热遇到一些问题学到了一些知希望与大家共勉~
1.关于各个单位之间的换算,有些题目中给出的单位不见得会是我们使用ansys的默认单位,需要换算一下,比如,比热和密度在mm单位中是需要换算的
2.关于接触单元的属性要设置为UX/UY/TEMP,因为在我看过的许多书中这个并没有提到 下边是我试验成功的摩擦声热的例子:
物块在以固定物块上滑动,求1m/s滑动15mm由于摩擦引起的升温
/PREP7
!*
ET,1,PLANE13
!*
ET,2,TARGE169
!*
ET,3,CONTA171
!*
KEYOPT,1,1,4
KEYOPT,1,3,0
!*
KEYOPT,3,1,1
MP,EX,1,69000
MP,PRXY,1,0.3
MP,DENS,1,2.7e-9
MP,ALPX,1,2.39e-5
MP,MU,1,0.2
MP,KXX,1,150
MP,C,1,9e8
!*
RECTNG,0,20,0,5,
RECTNG,0,5,5,10,
!*
LESIZE,ALL,1, , , ,1, , ,1,
AMESH,all
!*
LPLOT
LSEL,S, , , 3
NSLL,S,1
TYPE, 2
MAT, 1
R,1, , ,
!*
ESURF
ALLSEL,ALL
LSEL,S, , , 5
NSLL,S,1
TYPE, 3
ESURF
ALLSEL
FINISH
ALLSEL,ALL
/SOL
!*
ANTYPE,4
!*
TRNOPT,FULL
LUMPM,0
!*
!*
NLGEOM,1
SSTIF,0
NROPT,FULL, ,
TRNOPT,FULL
EQSLV, , ,0, ,DELE
PRECISION,0
MSAVE,0
TOFFST,460,
!*
TUNIF,0,
TIMINT,0,STRUCT
TIMINT,1,THERM
TIMINT,1,MAG
TINTP,0.005, , ,1,0.5,0.2,
!*
NSUBST,30,10000,30
OUTRES,ERASE
OUTRES,ALL,5
AUTOTS,1
TIME,0.015
APLOT
LSEL,S, , , 7
NSLL,S,1
D,all,UY,-0.1
ALLSEL,ALL
LSEL,S, , , 8
NSLL,S,1
D,all,UX,15
ALLSEL,ALL
LSEL,S, , , 1
NSLL,S,1
D,all, ,0, , , ,UX,UY, , , ,
ALLSEL,ALL
SOLVE
FINISH
上面是稍加修改的结果~
后来我又做了几个例子总是不成功,明明方法是一直的请大家鉴证以及批评指正~
两者不同之处大约也就是后者是两种材料罢了~~
/PREP7
!*
ET,1,PLANE13
!*
ET,2,TARGE169
!*
ET,3,CONTA171
!*
KEYOPT,1,1,4
KEYOPT,1,2,0
KEYOPT,1,3,0
KEYOPT,1,4,0
KEYOPT,1,5,0
!*
KEYOPT,3,1,1
KEYOPT,3,2,0
KEYOPT,3,3,0
KEYOPT,3,4,0
KEYOPT,3,5,0
KEYOPT,3,6,0
KEYOPT,3,7,0
KEYOPT,3,8,0
KEYOPT,3,9,0
KEYOPT,3,10,0
KEYOPT,3,11,0
KEYOPT,3,12,0
KEYOPT,3,14,0
!*
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,7.8e-9
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,1,,66.6
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,1,,4.6e8
MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,1,REFT,,,30
MPDATA,ALPX,1,,1.06e-5
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,206e5
MPDATA,PRXY,1,,0.3
*CSET,1,2, 1,2,
MPCOPY, ,1,2
TBCOPY,ALL,1,2
!*
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDE,DENS,2
MPDATA,DENS,2,,8.9e-9
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDE,C,2
MPDATA,C,2,,3.9e8
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDE,C,2
MPDATA,C,2,,3.9E+008
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDE,EX,2
MPDE,PRXY,2
MPDATA,EX,2,,2.06E+007
MPDATA,PRXY,2,,0.3
MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,2,REFT,,,20
MPDE,ALPX,2
MPDATA,ALPX,2,,1.75e-5
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDE,KXX,2
MPDATA,KXX,2,,383
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDE,EX,2
MPDE,PRXY,2
MPDATA,EX,2,,103e9
MPDATA,PRXY,2,,0.3
RECTNG,,5,,1.25,
RECTNG,,1.25,1.25,2.5,
ESIZE,0.25,0,
MSHAPE,0,2D
MSHKEY,0
!*
CM,_Y,AREA
ASEL, , , , 1
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
AMESH,_Y1
TYPE, 1
MAT, 2
REAL,
ESYS, 0
SECNUM,
CM,_Y,AREA
ASEL, , , , 2
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
AMESH,_Y1
!*
TYPE, 2
MAT, 2
REAL,
ESYS, 0
SECNUM,
TSHAP,LINE
!*
LSEL,S, , , 3
NSLL,S,1
TYPE, 2
MAT, 2
REAL,
ESYS, 0
SECNUM,
TSHAP,LINE
!*
FLST,5,21,1,ORDE,3
FITEM,5,22
FITEM,5,27
FITEM,5,-46
CM,_Y,NODE
NSEL, , , ,P51X
CM,_Y1,NODE
CMSEL,S,_Y
!*
CMSEL,,_Y1
ESURF, ,TOP,
CMSEL,,_Y
CMDELE,_Y
CMDELE,_Y1
!*
TYPE, 3
MAT, 2
REAL, 1
ESYS, 0
SECNUM,
TSHAP,LINE
!*
LSEL,S, , , 5
NSLL,S,1
FLST,5,6,1,ORDE,2
FITEM,5,127
FITEM,5,-132
CM,_Y,NODE
NSEL, , , ,P51X
CM,_Y1,NODE
CMSEL,S,_Y
!*
CMSEL,,_Y1
ESURF, ,TOP,
CMSEL,,_Y
CMDELE,_Y
CMDELE,_Y1
!*
FINISH
/SOL
!*
ANTYPE,4
!*
TRNOPT,FULL
LUMPM,0
!*
!*
NLGEOM,1
SSTIF,0
NROPT,FULL, ,
TRNOPT,FULL
EQSLV, , ,0, ,DELE
PRECISION,0
MSAVE,0
TOFFST,273,
!*
TUNIF,30,
TIMINT,0,STRUCT
TIMINT,1,THERM
TIMINT,1,MAG
TINTP,0.005, , , , , ,
!*
!*
NLGEOM,1
SSTIF,0
NROPT,FULL, ,
TRNOPT,FULL
EQSLV, , ,0, ,DELE
PRECISION,0
MSAVE,0
TOFFST,273,
!*
NSUBST,100,10000,30
OUTRES,ERASE
OUTRES,ALL,10
AUTOTS,1
TIME,3.75e-3
*DIM,press,TABLE,2,1,1,time, ,
!*
*SET,PRESS(1,0,1) , 0
*SET,PRESS(1,1,1) , 10000
*SET,PRESS(2,0,1) , 10
*SET,PRESS(2,1,1) , 10000
/PNUM,KP,0
/PNUM,LINE,1
/PNUM,AREA,0
/PNUM,VOLU,0
/PNUM,NODE,0
/PNUM,TABN,0
/PNUM,SVAL,0
/NUMBER,0
!*
/PNUM,ELEM,0
/REPLOT
!*
LPLOT
/AUTO,1
/REP,FAST
/AUTO,1
/REP,FAST
APLOT
LPLOT
APLOT
ALLSEL,ALL
FLST,2,1,4,ORDE,1
FITEM,2,7
/GO
!*
!*
SFL,P51X,PRES, %PRESS%
ALLSEL,ALL
FLST,2,1,4,ORDE,1
FITEM,2,6
!*
/GO
DL,P51X, ,UX,3.75
FLST,2,1,5,ORDE,1
FITEM,2,1
!*
/GO
DA,P51X,UX,
FLST,2,1,5,ORDE,1
FITEM,2,1
!*
/GO
DA,P51X,UY, |