同感，我也不知道。不过有个程序希望对你有用
关于橡胶等超弹性体的有限元分析ANSYS命令流
版本不近相同
（１）
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
（２）
*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,
（２）
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
（３）
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 !!!实常数１
R,2,,,fk1,ftoln !!!实常数１
R,3,,,fk1,ftoln !!!实常数１
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 |
IP卡
狗仔卡
发表于 2014-10-17 15:16:23
提升卡
沉默卡
喧嚣卡
变色卡
抢沙发
千斤顶
显身卡
