橡胶接触例子

meclbs

/PREP7  
!*  
ET,1,SOLID45
!*  
ET,2,HYPER58
MP,EX,1,2e11
MP,nuXY,1,0.3
mp,ex,2
MP,NUXY,2,0.4997   ! NUXY should be almost equal to, but less than 0.5
TB,MOONEY,2,1
TBDATA,1,525000
TBDATA,2,161000

!建立橡胶模型
csys,0
CYL4, , ,0.4,45,0.55,135,0.1
CYL4, , ,0.4,135,0.55,405,0.1
vglue,1,2
!划分网格
lesize,2, , ,20,-0.25, , , ,1
lesize,4, , ,20,-0.25, , , ,1
lesize,5, , ,20,-0.25, , , ,1
lesize,7, , ,20,-0.25, , , ,1

lesize,25, , ,25,-2, , , ,1
lesize,26, , ,25,-2, , , ,1
lesize,27, , ,25,-2, , , ,1
lesize,28, , ,25,-2, , , ,1

lesize,8, , ,5, , , , ,1
lesize,10, , ,3, , , , ,1

csys,1
type,2
mat,2
mshk,1
msha,0,3d
vmesh,1,3
!建立圆筒模型
local,11,0,,2.05-1e-6,-0.05
WPAVE,0,0,0
CYL4, , ,1.4,-110,1.5,-70,0.2

lesize,14, , ,20,-0.25, , , ,1
lesize,16, , ,20,-0.25, , , ,1
lesize,17, , ,20,-0.25, , , ,1
lesize,19, , ,20,-0.25, , , ,1

lesize,18, , ,3, , , , ,1
lesize,21, , ,5, , , , ,1

type,1
mat,1
mshk,1
msha,0,3d
vmesh,2
!建立接触
/COM, CONTACT PAIR CREATION - START
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,
MAT,1   
MP,EMIS,1,7.88860905221e-031
R,3
REAL,3  
ET,3,170
ET,4,174
R,3,,,1.0,0.1,0,
RMORE,,,1.0E20,0.0,1.0,
RMORE,0.0,0,1.0,,1.0,0.5
RMORE,0,1.0,1.0,0.0,,1.0
KEYOPT,4,4,0
KEYOPT,4,5,0
KEYOPT,4,7,0
KEYOPT,4,8,0
KEYOPT,4,9,0
KEYOPT,4,10,1   
KEYOPT,4,11,0   
KEYOPT,4,12,0   
KEYOPT,4,2,0
KEYOPT,3,5,0
! Generate the target surface   
ASEL,S,,,9  
CM,_TARGET,AREA
TYPE,3  
NSLA,S,1
ESLN,S,0
ESURF   
CMSEL,S,_ELEMCM
! Generate the contact surface  
ASEL,S,,,3  
CM,_CONTACT,AREA
TYPE,4  
NSLA,S,1
ESLN,S,0
ESURF   
ALLSEL  
ESEL,ALL
ESEL,S,TYPE,,3  
ESEL,A,TYPE,,4  
ESEL,R,REAL,,3  
/PSYMB,ESYS,1   
/PNUM,TYPE,1
/NUM,1  
EPLOT   
ESEL,ALL
ESEL,S,TYPE,,3  
ESEL,A,TYPE,,4  
ESEL,R,REAL,,3  
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  

/SOL
CSYS,1  
ALLSEL,ALL  
NSEL,S,LOC,X,0.4
d,all,all

local,12,1,,2.05-1e-6
nsel,s,loc,x,1.4
nsel,r,loc,y,-90
csys,0
cm,nc,node

allsel,all
csys,12
nsel,s,loc,x,1.4
csys,0
d,all,ux
d,all,uz
allsel,all

NLGEOM,1
NROPT,AUTO, ,   
LUMPM,0
EQSLV, , ,0,
PRECISION,0
MSAVE,0
PIVCHECK,1  
DSOPT,LOCAL,AUTO,   
DSPROC,1,   
SSTIF   
PSTRES  
TOFFST,0,   
!*  
!*  
TIME,1  
AUTOTS,-1   
DELTIM,0.02,0.00001,0.5,1
!cnvtol,f,500,0.01
neqit,100

csys,0
d,nc,uy,-0.001
solve
ddele,nc,uy
SFA,10,1,PRES,3e4
CSYS,1  

solve
finish

cao001

什么软件？
ANSYS？

byhyj

是一个好例子，但我提两个建议，一：命令流的使用很好，但是你对命令流的掌握还需要进一步加强。二：对橡胶分析来说，你对橡胶单元的选择需要进一步更新，现大有更好的单元给你选择如185单元，另外你的材料特性的定义也需要更新，现在有更好的定义材料特性的功能。尽量不要用原来的MOONEY 材料选项。

弄影清风

太好了,非常感谢,我正需要这些东西
两位高人

luxiaobin001

上面钢体的应力值大体相等，且在最大值附近，吗