關于橡膠等超彈性體的有限元分析ANSYS命令流

2017-01-16 by:CAE仿真在線來源:互聯網

版本不近相同

(1)

fini
/clear
/prep7
rid=1.35
rsd=0.12
et,1,56,,,1
et,2,169
et,3,171
tb,hyper,1,1,2,mooney
tbdata,1,80.0,20,0.002
mp,nuxy,1,.49
mp,mu,2,0.0
r,2,,,1,0.1
r,3,,,1,0.1
!*
local,11,1,rid+rsd/2
wpcsys,11
pcirc,rsd/2,,0,90
pcirc,rsd/2,,90,180
pcirc,rsd/2,,180,270
pcirc,rsd/2,,270,360
nummrg,kp
eshape,2
esize,.015
amesh,all
!
csys,0
k,20,rid-rsd,-rsd/2
k,21,rid+(2*rsd),-rsd/2
k,22,rid+(2*rsd),rsd/2
k,23,rid+(2*rsd),rsd/2
k,24,rid-rsd,rsd/2
l,21,20
l,22,21
l,24,23
!*
type,2
mat,2
real,2
lmesh,6
lmesh,9
real,3
lmesh,11
kmesh,24
csys,11
nsel,s,loc,x,rsd/2
type,3
real,2
esurf
real,3
esurf
esel,all
nsel,all
d,all,uz,0
n_load=node(kx(24),ky(24),0)
nplot

/user
save,oring,db
fini

/solu
nlgeom,on
solc,on
time,0.06
autots,off
nsubst,10,100,5
monitor,var3,n_load,fy
outres,all,all
dk,24,uy,-0.06
solve
time,0.085
dk,24,uy,-0.085
solve

(2)

*SET,r,1.55
/PREP7
ET,1,HYPER56 ! 2-D 4-NODE MIXED U-P HYPERELASTIC
ET,2,CONTA175
R,2, , ,-6000
R,3, , ,-6000,,-1
R,4, , ,-6000
ET,3,TARGE169
et,4,plane182
MP,EX,1,2.82 ! YOUNG'S MODULUS [MPA]
MP,NUXY,1,0.49967 ! POISSON'S RATIO
TB,MOONEY,1 ! MOONEY-RIVLIN CONSTANTS [MPA]
TBDATA,1,2.571
TBDATA,2,-0.818
mptemp,,,,,,,,
mptemp,1,0
MPdata,MU,1,,0.2
mptemp,,,,,,,,
mptemp,1,0
MPdata,EX,2,,206000
MPdata,PRXY,2,,0.3
CSYS,1 ! SWITCH TO CYLINDRICAL C.S.
K,1 ! DEFINE KEYPOINTS
K,2,R,-90
K,3,R
K,4,(0.5*R),-90
K,5,(0.6*R),-45
K,6,(0.5*R)
K,7,R,-45
L,2,7
L,7,3
CSYS,0 ! SWITCH TO CARTESIAN C.S.
A,2,7,5,4
A,7,3,6,5
A,4,5,6,1
ARSYM,X,all
ARSYM,Y,all
NUMMRG,KP
ESIZE,,6 ! SET ELEMENT DIVISION SIZE
AMESH,ALL
csys,0
rectng,-2*r,2*r,r,2*r
rectng,-2*r,2*r,-r,-2*r

rectng,1.8,2*r,-r,0.5
AADD,14,15 !
lsel,s,line,,13
lesize,all,,,30,-0.5
lsel,s,line,,42
lesize,all,,,28,-0.5
lsel,s,line,,40
lesize,all,,,20,2
LSEL,S, , ,33,36,3
lesize,all,,,8
LSEL,S, , ,18,26,8
lesize,all,,,8
LSEL,S, , , 39
lesize,all,,,5
LSEL,S, , ,41
lesize,all,,,8

LSEL,S, , ,24,31,7
lesize,all,,,20
TYPE,4
MAT,2
ASEL, , , ,13,16,3
AMESH,all
ALLSEL,ALL
lsel,s,line,,19,28,9
NSLL,S,1
NPLOT
type,2
real,2
esurf
lsel,s,line,,13
NSLL,S,1
NPLOT
type,3
real,2
esurf
lsel,s,line,,2,23,21
NSLL,S,1
csys,1
NSEL,U,LOC,Y,-45,-22.5
type,2
real,3
esurf
lsel,s,line,,40
NSLL,S,1
type,3
real,3
esurf
lsel,s,line,,1,10,9
lsel,a,line,,2
NSLL,S,1
NSEL,U,LOC,Y,-22.5,0
type,2
real,4
esurf
lsel,s,line,,42
NSLL,S,1
type,3
real,4
esurf
FINISH
/SOL
lsel,s,line,,31
NSLL,S,1
d,all,all
ALLSEL,ALL
csys,1
nsel,s,loc,y,270
nsel,r,loc,x,r
d,all,uy
lsel,s,line,,24
NSLL,S,1
d,all,ux
d,all,uy,-0.62
ALLSEL,ALL
NSEL,s,LOC,y,90
d,all,ux
ALLSEL,ALL
ANTYPE,0
NLGEOM,1
NSUBST,100,0,0
OUTRES,ALL,ALL
LSWRITE,1,
lsel,s,line,,10,14,4
lsel,a,line,,28,32,4
NSLL,R,1
NPLOT
NSEL,U,LOC,Y,90
NSEL,U,LOC,Y,270
SF,all,PRES,0.8
NSEL,s,LOC,y,90
ddele,all,all
nsel,s,loc,y,270
nsel,r,loc,x,r
ddele,all,all
ALLSEL,ALL
NSUBST,40,0,0
LSWRITE,2,
LSSOLVE,1,2,1,

(2)

fini
/cle
/Title,O-ring analysis
!!!參數設定
r1=0.2 !!!圓角
r2=0.6 !!!圓角
dia=7 !!!O圈直徑
w1=9.5 !!!槽寬
h1=5 !!!槽高
wlr=w1/2 !!!左右寬
wt=2*w1 !!!總寬
gap=dia-h1 !!!間隙
/prep7
et,1,plane182
mp,ex,1,2.82e+6
mp,nuxy,1,0.49
mp,mu,1,0.1
c10=0.293
c01=0.177
nu1=0.49967
DD=2*(1-2*nu1)/(c10+c01)
tb,hyper,2,1,2,mooney
tbdata,1,c10,c01,DD
!!!!!o圈模型
csys,1
k,1,dia/2
kgen,36,1,,,,10
/pnum,kp,1
kplot
*do,j,1,35,1
l,j,j+1
*enddo
l,36,1
lsel,all
al,all
!!!網格劃分
type,1
mat,1
lesize,all,,,4
mshkey,0
amesh,1
save,oring_1,db
!!!上下接觸邊
csys,0
k,40,-wt/2,dia/2
k,41,wt/2,dia/2
k,42,-wt/2,h1-dia/2
k,43,-wt/4,h1-dia/2
k,44,-wt/4,-dia/2
k,45,wt/4,-dia/2
k,46,wt/4,h1-dia/2
k,47,wt/2,h1-dia/2
l,40,41
l,42,43
*repeat,5,1,1
/pnum,line,1
/rep
lplot
lfil,38,39,r1 !43
lfil,39,40,r2 !44
lfil,40,41,r2 !45
lfil,41,42,r1 !46
!!!!create the contact pair,define element and contact parameter
et,2,169
et,3,171
keyopt,3,5,4
keyopt,3,9,1
keyopt,3,10,1
keyopt,3,12,2
mp,mu,1,0
r,1,,,0.1,0.1
r,2,,,0.1,0.1
r,3,,,0.1,0.1
!!!pair1
type,2
real,1
mat,1
lesize,37,,,1
lmesh,37
kmesh,40
!!!contact face
csys,1
lsel,s,loc,x,dia/2
lsel,r,loc,y,-225,-315
type,3
nsll,s,1
esln,s,0
esurf,all
alls
lplot
!!!!contact pair2
type,2
real,2
csys,0
lsel,s,line,,38,46
lesize,all,,,1
lmesh,all
!!!contact face 2
csys,1
lsel,s,loc,x,dia/2
lsel,r,loc,y,-45,-135
type,3
nsll,s,1
esln,s,0
esurf,all
alls
!!!!apply load
n_load=node(kx(40),ky(40),0)
fini
/solu
nlgeom,on
solc,on
time,0.85
d,n_load,uy,-gap
csys,1
lsel,s,loc,x,dia/2
nsll,s,1
nsel,r,loc,y,90
nsel,r,loc,x,dia/2
d,all,ux
lsel,s,loc,x,dia/2
nsll,s,1
nsel,r,loc,y,270
nsel,r,loc,x,dia/2
d,all,uy
nsubst,25,500,10
outres,all,all
monitor,var3,n_load,fy
alls
solve

(3)

fini
/clear
/units,mpa
/Title,O-ring analysis
!!!參數設定
fe=0.9 !!!摩擦系數
fk1=10 !!!法向接觸剛度
fk2=1 !!!法向接觸剛度
fd=3 !!!切向接觸剛度
ftoln=1 !!!初始滲透因子
r1=0.2 !!!圓角
r2=0.6 !!!圓角
dia=7 !!!O圈直徑
w1=9.5 !!!槽寬
h1=5 !!!槽高
wlr=w1/2 !!!左右寬
wt=2*w1 !!!總寬
gap=dia-h1 !!!間隙
pres=1 !!!壓縮量(mm)
hd=2*h1 !!!O圈槽段總高
ht=hd !!!無槽側總高
oex=2.82
c10=0.293
c01=0.177
nu1=0.49967
DD=2*(1-2*nu1)/(c10+c01)
!!!環境設置
/prep7
et,1,plane182
keyopt,1,6,1
et,2,conta172
et,3,targe169
keyopt,2,5,4
keyopt,2,9,1
keyopt,2,10,1
keyopt,2,12,2
et,4,conta172
et,5,targe169
keyopt,4,5,4
keyopt,4,9,1
keyopt,4,10,1
keyopt,4,12,2
et,6,conta172
et,7,targe169
keyopt,6,5,4
keyopt,6,9,1
keyopt,6,10,1
keyopt,6,12,2
R,1,,,fk1,ftoln !!!實常數1
R,2,,,fk1,ftoln !!!實常數1
R,3,,,fk1,ftoln !!!實常數1
mp,ex,1,210e+6 !!!本體與筒體材料彈模
mp,nuxy,1,0.3 !!!本體與筒體材料彈模
mp,mu,1,fe
mp,ex,2,2.82e+3 !!!O圈彈模
mp,nuxy,2,nu1
mp,mu,2,fe
tb,hyper,2,,2,mooney
tbdata,1,c10,c01,DD

!!!!翻蓋本體側
k,1
k,2,0,hd
k,3,wlr-r1,hd
k,4,wlr,hd-r1
k,5,wlr,hd-h1+r2
k,6,wlr+r2,hd-h1
k,7,wt-wlr-r2,hd-h1
k,8,wt-wlr,hd-h1+r2
k,9,wt-wlr,hd-r1
k,10,wt-wlr+r1,hd
k,11,wt,hd
k,12,wt
k,21,wlr-r1,hd-r1
k,22,wlr+r2,hd-h1+r2
k,23,wt-wlr-r2,hd-h1+r2
k,24,wt-wlr+r1,hd-r1
l,1,2
l,2,3
larc,3,4,21,r1
l,4,5
larc,5,6,22,r2
l,6,7
larc,7,8,23,r2
l,8,9
larc,9,10,24,r1
l,10,11
l,11,12
l,12,1
al,all
/pnum,line,1
lplot
!llist
!!!!!筒體側
k,31,,hd+gap
k,32,,hd+gap+ht
k,33,wt,hd+gap+ht
k,34,wt,hd+gap
l,31,32
l,32,33
l,33,34
l,34,31
lsel,s,loc,y,hd+gap/2,2*hd+gap
lplot
al,13,14,15,16
allsel
aplot
lsel,s,,,12,16,2
lesize,all,,,20
lsel,s,,,11,15,2
lsel,a,,,1
lesize,all,,,10
lsel,s,,,6
lesize,all,,,10
lsel,s,,,2,10,8
lesize,all,,,5
lsel,s,,,4,8,4
lesize,all,,,5
lsel,s,,,3,9,6
lesize,all,,,2
lsel,s,,,5,7,2
lesize,all,,,2
mat,1
asel,all
mshkey,0
amesh,all
!lsel,s,,,2,10,1
!lsel,a,,,16
!lrefine,all,,,1 !!!!細化

!!!O圈模型
!!!圓心位置
!!!偏移workplane
!local,11,1,wt-wlr-dia/2,hd-h1+dia/2
local,11,1,wt/2,hd-h1+dia/2
csys,11
k,101,dia/2
NUMSTR,kp,102
kgen,36,101,,,,10
/pnum,kp,1
kplot
*do,j,101,135,1
l,j,j+1
*enddo
l,136,101
lsel,s,loc,x,dia/2
lsel,u,line,,6,16,10
lplot
al,all
mat,2
lesize,all,,,2
asel,s,area,,3
amesh,all
!!!!創建O圈與筒體接觸
csys,0
lsel,s,line,,16
type,3
real,1
nsll,s,1
esln,s,0
esurf,all
allsel
csys,11
lsel,s,loc,y,-225,-315
lsel,r,loc,x,dia/2
lsel,u,line,,16
type,2
nsll,s,1
esln,s,0
esurf,all
!!!創建O圈與槽底接觸
csys,0
alls
lplot
lsel,s,line,,6
type,5
real,2
nsll,s,1
esln,s,0
esurf,all
alls
csys,11
lsel,s,loc,y,225,315
lsel,r,loc,x,dia/2
lsel,u,line,,6
type,4
nsll,s,1
esln,s,0
esurf,all
!!!!創建O圈與槽右側接觸
csys,0
alls
lsel,s,line,,8
type,7
real,3
nsll,s,1
esln,s,0
esurf,all
alls
csys,11
lsel,s,loc,y,-45,45
lsel,r,loc,x,dia/2
lsel,u,line,,8
type,6
nsll,s,1
esln,s,0
esurf,all
!!!!loads and solve加載與求解
/solu
csys,11
lsel,s,loc,x,dia/2
nsll,s,1
nsel,r,loc,y,270
nsel,r,loc,x,dia/2
!d,all,ux
lsel,s,loc,x,dia/2
nsll,s,1
nsel,r,loc,y,90
nsel,r,loc,x,dia/2
d,all,ux
csys,0
nsel,s,loc,y
d,all,all
nsel,s,loc,y,2*hd+gap
d,all,ux
outres,all,all
time,1
nsubst,20
nlgeom,on
lnsrch,on
nsel,s,loc,y,2*hd+gap
d,all,uy,-gap
allsel
solve
/solu
outres,all,all
time,2
nsubst,20
nlgeom,on
lnsrch,on
csys,11
lsel,s,loc,x,dia/2
nsll,s,1
nsel,r,loc,y,270
nsel,r,loc,x,dia/2
ddele,all,ux
d,all,uy
lsel,s,loc,x,dia/2
nsll,s,1
nsel,r,loc,y,90
nsel,r,loc,x,dia/2
ddele,all,ux
d,all,uy
lsel,s,loc,x,dia/2
nsll,s,1
nsel,r,loc,y,0
nsel,r,loc,x,dia/2
d,all,uy
lsel,s,loc,x,dia/2
nsll,s,1
nsel,r,loc,y,180
nsel,r,loc,x,dia/2
d,all,uy
lsel,s,loc,y,135,225
lsel,r,loc,x,dia/2
sfl,all,pres,60000
alls
solve

開放分享：優質有限元技術文章,助你自學成才