找回密码
 注册
Simdroid-非首页
查看: 61|回复: 1

[接触分析] 为何我施加的是压力,经过接触分析反而会使得物件飞起来

[复制链接]
发表于 2011-11-16 22:51:24 | 显示全部楼层 |阅读模式 来自 清华大学
我想做一个一直有压力存在,多个板摩擦的模型。以后在此基础上改进。
由于三层接触问题就更难了,所以我从简单的做起,把上面两个直接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




 楼主| 发表于 2011-11-16 22:55:10 | 显示全部楼层 来自 清华大学
Simdroid开发平台
求大家帮帮忙,谢谢了!~~~~
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Simapps系列直播

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-9-29 17:25 , Processed in 0.026941 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表