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

2017-01-16  by:CAE仿真在線  來源:互聯(lián)網(wǎng)

版本不近相同

(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
!!!參數(shù)設定
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
!!!網(wǎng)格劃分
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
!!!參數(shù)設定
fe=0.9 !!!摩擦系數(shù)
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 !!!無槽側(cè)總高
oex=2.82
c10=0.293
c01=0.177
nu1=0.49967
DD=2*(1-2*nu1)/(c10+c01)
!!!環(huán)境設置
/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 !!!實常數(shù)1
R,2,,,fk1,ftoln !!!實常數(shù)1
R,3,,,fk1,ftoln !!!實常數(shù)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
!!!!翻蓋本體側(cè)
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
!!!!!筒體側(cè)
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
!!!!創(chuàng)建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
!!!創(chuàng)建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
!!!!創(chuàng)建O圈與槽右側(cè)接觸
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

開放分享:優(yōu)質(zhì)有限元技術文章,助你自學成才

相關標簽搜索:關于橡膠等超彈性體的有限元分析ANSYS命令流 Ansys有限元培訓 Ansys workbench培訓 ansys視頻教程 ansys workbench教程 ansys APDL經(jīng)典教程 ansys資料下載 ansys技術咨詢 ansys基礎知識 ansys代做 Fluent、CFX流體分析 HFSS電磁分析 Abaqus培訓 

編輯
在線報名:
  • 客服在線請直接聯(lián)系我們的客服,您也可以通過下面的方式進行在線報名,我們會及時給您回復電話,謝謝!
驗證碼

全國服務熱線

1358-032-9919

廣州公司:
廣州市環(huán)市中路306號金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓QQ咨詢:點擊咨詢 點擊咨詢
項目QQ咨詢:點擊咨詢
email:kf@1cae.com