我想做一个一直有压力存在,多个板摩擦的模型。以后在此基础上改进。
由于三层接触问题就更难了,所以我从简单的做起,把上面两个直接GLUE在了一起,可是现在还是无法运行,这是为什么?
为什么我施加压力,反而会使得物件飞起来?
finish
/clear
/filn,2D-1116-04F
/title,2D-shiyan
/PREP7
!*选取单元,2D耦合场单元
ET,1,PLANE13
ET,2,TARGE169
ET,3,CONTA171
KEYOPT,1,1,4
KEYOPT,1,3,0
KEYOPT,3,1,1
!*输入材料参数
!*材料1,铝 AL3003-H18-Foil
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,6.89e10
MPDATA,PRXY,1,,0.33
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,2.73e3
MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,1,REFT,,,293
MPDATA,ALPX,1,,2.5e-5
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,1,,154
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,1,,0.893e3
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDE,C,1
MPDATA,C,1,,8.93E+002
!*材料2,钛 Ti-6Al-4V -- sonotrode
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,2,,1.138e11
MPDATA,PRXY,2,,0.342
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,2,,4.43e3
MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,2,REFT,,,
MPDATA,ALPX,2,,9.2e-6
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,2,,6.7
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,2,,0.5263e3
!*材料3,AL3003-H16-baseplate
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,3,,6.89e10
MPDATA,PRXY,3,,0.33
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,3,,2.73e3
MPTEMP,,,,,,,,
MPTEMP,1,0
UIMP,3,REFT,,,293
MPDE,ALPX,3
MPDATA,ALPX,3,,2.5E-005
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,KXX,3,,157
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,C,3,,0.893e3
!*建立模型
RECTING,-30,30,,5
RECTING,0,5,5,6
RECTING,0,RECTING,-30,30,,5
!*粘贴电极和F板
RECTING,0,5,5,6
RECTING,0,5,6,10
FLST,2,2,5,ORDE,2
FITEM,2,2
FITEM,2,-3
AGLUE,P51X
)/GOP !5,6,10
!*划分网格
FLST,5,2,4,ORDE,2
FITEM,5,13
FITEM,5,-14
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,8, , , , ,1
!*
FLST,5,2,4,ORDE,2
FITEM,5,6
FITEM,5,8
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,4, , , , ,1
!*
FLST,5,2,4,ORDE,2
FITEM,5,2
FITEM,5,4
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,8, , , , ,1
!*
FLST,5,3,4,ORDE,3
FITEM,5,5
FITEM,5,7
FITEM,5,11
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,5, , , , ,1
!*
FLST,5,2,4,ORDE,2
FITEM,5,1
FITEM,5,3
CM,_Y,LINE
LSEL, , , ,P51X
CM,_Y1,LINE
CMSEL,,_Y
!*
LESIZE,_Y1, , ,60, , , , ,1
!*
TYPE, 1
MAT, 2
REAL,
ESYS, 0
SECNUM,
!*
MSHAPE,0,2D
MSHKEY,0
!*
CM,_Y,AREA
ASEL, , , , 4
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
TYPE, 1
MAT, 1
REAL,
ESYS, 0
SECNUM,
!*
MSHKEY,1
!*
CM,_Y,AREA
ASEL, , , , 2
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
TYPE, 1
MAT, 3
REAL,
ESYS, 0
SECNUM,
!*
/AUTO,1
/REP,FAST
APLOT
CM,_Y,AREA
ASEL, , , , 1
CM,_Y1,AREA
CHKMSH,'AREA'
CMSEL,S,_Y
!*
AMESH,_Y1
!*
CMDELE,_Y
CMDELE,_Y1
CMDELE,_Y2
!*
!*建立接触
!*
CM,_NODECM,NODE
CM,_ELEMCM,ELEM
CM,_KPCM,KP
CM,_LINECM,LINE
CM,_AREACM,AREA
CM,_VOLUCM,VOLU
/GSAV,cwz,gsav,,temp
MP,MU,1,0.2
MAT,1
R,3
REAL,3
ET,4,169
ET,5,172
KEYOPT,5,9,0
KEYOPT,5,10,2
R,3,
RMORE,
RMORE,,0
RMORE,0
KEYOPT,5,1,1
! Generate the target surface
LSEL,S,,,5
CM,_TARGET,LINE
TYPE,4
NSLL,S,1
ESLN,S,0
ESURF
CMSEL,S,_ELEMCM
! Generate the contact surface
LSEL,S,,,3
CM,_CONTACT,LINE
TYPE,5
NSLL,S,1
ESLN,S,0
ESURF
!* Create Companion Pair - Start
R,4
REAL,4
ET,6,169
ET,7,172
KEYOPT,7,9,0
KEYOPT,7,10,2
R,4,
RMORE,
RMORE,,0
RMORE,0
KEYOPT,7,1,1
TYPE,6
ESEL,S,TYPE,,5
NSLE,S
ESLN,S,0
ESURF
TYPE,7
ESEL,S,TYPE,,4
NSLE,S
ESLN,S,0
ESURF
!* Create Companion Pair - End
ALLSEL
ESEL,ALL
ESEL,S,TYPE,,4
ESEL,A,TYPE,,5
ESEL,R,REAL,,3
/PSYMB,ESYS,1
/PNUM,TYPE,1
/NUM,1
EPLOT
ESEL,ALL
ESEL,S,TYPE,,4
ESEL,A,TYPE,,5
ESEL,R,REAL,,3
ESEL,A,TYPE,,6
ESEL,A,TYPE,,7
CMSEL,A,_NODECM
CMDEL,_NODECM
CMSEL,A,_ELEMCM
CMDEL,_ELEMCM
CMSEL,S,_KPCM
CMDEL,_KPCM
CMSEL,S,_LINECM
CMDEL,_LINECM
CMSEL,S,_AREACM
CMDEL,_AREACM
CMSEL,S,_VOLUCM
CMDEL,_VOLUCM
/GRES,cwz,gsav
CMDEL,_TARGET
CMDEL,_CONTACT
/MREP,EPLOT
allsel
!*施加结构载荷
*DIM,PRESS,TABLE,2,1,1,TIME, ,
!*
*SET,PRESS(1,0,1) , 0
*SET,PRESS(1,1,1) , 400
*SET,PRESS(2,0,1) , 10
*SET,PRESS(2,1,1) , 400
FINISH
/SOL
FLST,2,1,4,ORDE,1
FITEM,2,11
/GO
!*
!*
SFL,P51X,PRES, %PRESS%
!*其它方向的力约束
/SOL
FLST,2,6,4,ORDE,6
FITEM,2,2
FITEM,2,4
FITEM,2,6
FITEM,2,8
FITEM,2,13
FITEM,2,-14
!*
/GO
DL,P51X, ,UX,0
FLST,2,1,4,ORDE,1
FITEM,2,1
!*
/GO
DL,P51X, ,UY,0
)/GOP !
!*初始温度场
SFL,P51X,CONV,40, ,20,
TUNIF,20,
!*设置动力学求解选项,rampped加载方式
!*
ANTYPE,4
!*
TRNOPT,FULL
LUMPM,0
!*
ANTYPE,4
NLGEOM,1
NSUBST,200,2000,200
OUTRES,ERASE
OUTRES,ALL,1
KBC,0 !*rampped加载
TIME,4
!*全选求解
allsel
solve
|