橡胶参数.docx
《橡胶参数.docx》由会员分享,可在线阅读,更多相关《橡胶参数.docx(74页珍藏版)》请在冰豆网上搜索。
橡胶参数
超弹性(2008-07-1715:
32:
31)
标签:
ansys 教育
分类:
ANSYS学习
铁板材料特性:
Ex=2e5 (泊松比)=0.3 橡胶体材料特性:
(泊松比)=0.499
单轴拉伸时的实验数据:
应变
0.0
0.2
0.4
0.6
0.8
1.0
1.2
1.6
2
3
应力
0.0
1
1.5
2.0
2.9
3.6
5
7.5
9.7
17
单轴压缩时的实验数据:
应变
0.45
0.4
0.35
0.3
0.25
0.2
0.15
0.1
0.05
应力
256
128
64
32
16
8
4
2
1
剪切时的实验数据:
应变
0.0
0.05
0.2
0.3
0.4
0.5
0.6
0.8
1
1.5
应力
0.0
0.7
1.5
2.0
2.9
3.6
5
7.5
9.7
17
fini
/cle
r=180
l1=185
l2=74
h1=6
h2=50
h3=182
r1=10
d=50
/prep7
et,1,56,,,1
et,2,42,,,1
et,3,48,,1
et,4,58
et,5,45
et,6,49,,1
keyopt,6,7,1
keyopt,3,7,1
r,1,10000,,0.5
rect,0,l1,0,h1
rect,0,l2,0,h2
cyl4,0,h3,r,-90,,0
aovlap,all
asel,s,loc,y,0,h1
asel,r,loc,x,0,l2
aadd,all
alls
asel,s,loc,y,h1,h3
aadd,all
alls
lsel,s,loc,x,0
lsel,r,loc,y,0,h1
lcom,all
alls
lsel,s,loc,x,0
lsel,r,loc,y,h1,h3
lcom,all
alls
lsel,s,loc,y,h1
lsel,r,loc,x,0,l2
lcom,all
alls
lsel,s,loc,y,h1,h2
lsel,r,loc,x,l2
*get,line1,line,,num,max
alls
lsel,s,loc,x,l2,l1
lsel,r,loc,y,h2,h3
*get,line2,line,,num,max
alls
lfillt,line1,line2,r1
alls
al,1,4,3
asel,s,loc,y,h1,h3
aadd,all
alls
lsel,s,loc,x,l2,l1
lsel,r,loc,y,h1+0.1,h3-0.1
lcom,all
alls
lsel,s,loc,x,l2,l1
lsel,r,loc,y,h1+0.1,h3-0.1
*get,line3,line,,num,max
alls
kl,line3,0.12
lsel,s,loc,x,0
lsel,r,loc,y,h1,h3
*get,line4,line,,num,max
kl,line4,0.4
kl,line4,0.7
alls
lstr,5,7
lstr,13,8
asel,s,loc,y,h1,h3
lsel,s,,,3,4
asbl,all,all
alls
asel,s,loc,y,h1,h3
aatt,1,1,1
asel,s,loc,y,0,h1
aatt,2,1,2
alls
lesize,3,,,10
lesize,8,,,12,2
lesize,7,,,6
lesize,12,,,1
lesize,19,,,8
mshkey,1
amesh,5
amesh,3
amesh,2
amesh,1
amesh,6
alls
lsel,s,,,19
nsll,s,1
cm,targ,node
alls
lsel,s,,,6
nsll,s,1
cm,cont1,node
lsel,s,,,8
nsll,s,1
nsel,r,loc,y,h1,130
cm,cont2,node
alls
cmsel,s,cont1
cmsel,a,cont2
!
cmsel,a,cont3
cm,cont,node
alls
type,3
mat,1
real,1
gcgen,cont,targ
alls
save,hypelastic,db
resume,hypelastic,db
mp,nuxy,1,0.499
mp,ex,2,2e5
mp,nuxy,2,0.3
m=1.95
m1=2.03
n=1.05
*dim,strn,,10,3
*dim,strss,,10,3
*dim,const,,5
*dim,calc,,30
*dim,sortss,,30
*dim,sortsn,,30
*dim,ffx,table,30,2
*dim,ecalc,table,100
*dim,xval,table,100
strn(1,2)=-0.45,-0.4,-0.35,-0.3,-0.25,-0.2,-0.15,-0.1,-0.05
strss(1,2)=-256,-128,-64,-32,-16,-8,-4,-2,-1
*do,i,1,9
strss(i,2)=strss(i,2)*m
strn(i,2)=strn(i,2)*m1
*enddo
strn(1,1)=0.0,0.2,0.4,0.6,0.8,1.0,1.2,1.6,2,3
strss(1,1)=0.0,1,1.5,2.0,2.9,3.6,5,7.5,9.7,17
strn(1,3)=0.0,0.05,0.2,0.3,0.4,.5,.6,.8,1,1.5
strss(1,3)=0.0,0.7,1.5,2.0,2.9,3.6,5,7.5,9.7,17
*do,i,1,10
strss(i,3)=strss(i,3)*n
*enddo
tb,mooney,,,,1
*mooney,strn(1,1),strss(1,1),,const
(1),calc
(1),sortsn
(1),sortss
(1)
*vfun,ffx(1,1),copy,sortss
(1)
*vfun,ffx(1,2),copy,calc
(1)
*vplot,strn
(1),ffx
(1),2
*eval,1,2,const
(1),-0.2,0,xval
(1),ecalc
(1)
*vplot,xval
(1),ecalc
(1)
fini
/solu
alls
nsel,s,loc,y,0
d,all,all,0
alls
nsel,s,loc,x,0
d,all,ux,0
d,all,uz,0
alls
nsel,s,loc,y,h3
d,all,ux,0
d,all,uz,0
d,all,uy,-d
alls
antype,static
nlgeom,on
nropt,,,on
outpr,all,all
outres,all,all
autots,on
time,1
deltim,0.03,0.01,0.3
cnvtol,f,,0.02,2
lnsrch,on
pred,on
alls
solve
fini
ansys-Beam3二维弹性单元特性翻译
工程应力与真实应力(2008-07-3113:
47:
36)
标签:
ansys 教育
分类:
ANSYS学习
fini
/cle
/PREP7
ET,1,plane182
KEYOPT,1,3,1
R,1,0.001,,,,,,
MP,EX,1,2.1E11 !
STEEL
MP,NUXY,1,0.3
TB,BKIN,1,1 !
DEFINENON-LINEARMATERIALPROPERTYFORSTEEL
TBTEMP,0
TBDATA,1,210e6,8.6e9
BLC4,0,,0.03,0.03
NUMCMP,ALL
AESIZE,1,0.003,
alls
amesh,all
nsel,s,loc,y,0
nsel,a,loc,y,0.05
d,all,uy
nsel,s,loc,y,0.03
nsel,a,loc,y,0.08
D,all,,0.00009,,,,Uy !
if>0.00003materialisyield;
alls
/SOLU
NLGEOM, ON
!
Accordingsmallstraintheory0.005cause0.3%=(0.00009/0.03)strain;butifwetrunNLGEOMON,thestrainis0.2996%=ln(1+0.00009/0.03)
NSUBST, 40, 100, 40
OUTRES, ALL, 1
SOLVE
/POST1
SET, LAST
PLNSOL,EPTO,Y,0,1.0 !
themaxtotalstrainvalueis0.2996%
/repl
/POST26
RFORCE, 2, 22, f, y, FY_2
PLVAR, 2
ANSOL,4,22,EPEL,y,EPELy_2
ANSOL,5,22,EPPL,y,EPPLy_4
ANSOL,6,22,S,y,Sy_4
ADD,7,4,5,,,,,1,1,1,
/AXLAB,X,DEFLECTION
/AXLAB,Y,Stress
/GRID,1
XVAR,7
PLVAR,6
超级大变形(2008-07-1515:
56:
19)
标签:
ansys
分类:
ANSYS学习
fini
/cle
/PREP7
lsize=600
hsize=24
l=135e-3
h=6e-3
p=-10e-3
!
定义单元
!
ET,1,VISCO108
ET,1,SHELL181
ET,2,MASS21
KEYOPT,2,1,0
KEYOPT,2,2,0
KEYOPT,2,3,2
KEYOPT,1,1,0
KEYOPT,1,2,0
KEYOPT,1,3,0
KEYOPT,1,5,0
KEYOPT,1,6,0
KEYOPT,1,7,0
KEYOPT,1,8,0
KEYOPT,1,9,0
KEYOPT,1,11,0
!
实常数
R,1,0.3e-3,,,,,, !
板厚
RMORE,,,,
RMORE
RMORE,,
R,2,(5.3e-3)/(hsize-1), !
集中质量
!
定义材料
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,EX,1,,2.06e11
MPDATA,PRXY,1,,0.3
MPTEMP,,,,,,,,
MPTEMP,1,0
MPDATA,DENS,1,,7900
!
板尺寸
BLC4,p,,l,h
!
单元大小
LESIZE,1,,,lsize,,,,,1
LESIZE,2,,,hsize,,,,,1
MSHAPE,0,2D
MSHKEY,1
AMESH,1
TYPE,2
REAL,2
nsel,s,loc,x,l+p
nsel,r,ext
nsel,u,loc,y,0
nsel,u,loc,y,h
*get,n_min,node,,num,min
n_num=n_min
*do,i,1,hsize-1
E,n_num
n_num=ndnext(n_num)
*enddo
!
约束
DL,4,,ALL
EPLOT
alls
save
!
求解
/SOLU
ANTYPE,4 !
瞬态大变形
TRNOPT,FULL
LUMPM,0
NLGEOM,1
NSUBST,20,100,10 !
子步数
OUTRES,ERASE
OUTRES,esol,LAST
OUTRES,nsol,LAST
AUTOTS,1
lnsrch,1
PRED,ON
SSTIF,1
KBC,0
!
cnvtol,f,0.05,,, !
收敛容差
!
cnvtol,u,0.05,,,
!
冲击波加载
J=1
*do,i,1.1e-5,1.1e-3,1.1e-4 !
载荷步数可在此改
time,i
acel,,,9810*sin(2854.5*i) !
冲击波形
lswrite,j
j=j+1
*enddo
*do,i,1.1e-3+1e-3,30e-3,1e-3 !
载荷步数可在此改
time,i
NSUBST,50,100,10
acel,,,-361.94*sin(108.7*(i-1.1e-3)) !
冲击波形
lswrite,j
j=j+1
*enddo
acel,,,0
*do,i,30e-3+1e-3,0.05,1e-3 !
载荷步数可在此改
time,i
NSUBST,50,100,10
lswrite,j
j=j+1
*enddo
j=j-1
save
!
求解
lssolve,1,J,1,
Ansys疲劳算例(2008-02-2213:
38:
45)
标签:
ansys 教育
!
***************环境设置***************
/units,si
/title,Fatigueanalysisofcylinderwithflathead
!
***************参数设定***************
Di=1000 !
筒体内径
t=20 !
筒体厚度
hc=nint(4*sqrt(Di/2*t)/10)*10 !
模型中筒体长度
tp=60 !
平板封头厚度
r1=10 !
平板封头外测过渡圆弧半径
r2=10 !
平板封头内侧应力释放槽圆弧半径
exx=2e5 !
材料弹性模量
mu=0.3 !
材料泊松比
p1=2 !
最高工作压力
p3=2.88 !
水压试验压力
n1=2e4 !
最高/最低压力循环次数
n2=5 !
水压试验次数
!
***************前处理***************
/prep7
et,1,82 !
设定单元类型
keyopt,1,3,1 !
设定周对称选项
mp,ex,1,exx !
定义材料弹性模量
mp,nuxy,1,mu !
定义材料泊松比
!
******* 建立模型 *******
k,1,0,0 !
定义关键点
k,2,Di/2+t,,
k,3,Di/2+t,-(tp+hc)
k,4,Di/2,-(tp+hc)
k,5,Di/2,-tp
k,6,Di/2-r2,-tp !
定义应力释放槽圆弧中心关键点
k,7,0,-tp
l,1,2 !
生成线
l,2,3
l,3,4
l,4,5
l,5,7
l,7,1
LFILLT,1,2,r1 !
生成外测过渡圆弧
al,all !
生成子午面
CYL4,kx(6),ky(6),r2,180 !
生成应力释放槽面域
ASBA,1,2 !
面相减
wprot,,,90 !
旋转工作平面
wpoff,,,kx(6)-3*r2 !
移动工作平面
asbw,all !
用工作平面切割子午面
wprot,,90 !
旋转工作平面
wpoff,,,tp+r2 !
移动工作平面
asbw,all !
用工作平面切割子午面
esize,5 !
设定单元尺寸
MSHKEY,1 !
设定映射剖分
amesh,1 !
映射剖分面1
amesh,3 !
映射剖分面3
esize,2 !
设定单元尺寸
MSHKEY,0 !
设定自由剖分
amesh,4 !
自由剖分面4
fini !
退出前处理
!
***************求解***************
/solu !
筒体端部施加轴向约束
dl,3,,uy !
筒体端部施加轴向约束
dl,6,,symm !
平板封头对称面施加对称约束
time,1 !
载荷步1
lsel,s,,,8 !
选择内表面各线段
lsel,a,,,11,13
lsel,a,,,15
cm,lcom1,line !
生成内表面线组件
SFL,all,PRES,p1, !
内表面施加内压
alls !
全选
solve !
求解
fini !
退出求解器
!
***************后处理***************
/post1 !
进入后处理
FTSIZE,1,2,2, !
设定疲劳评定的位置数、事件数及载荷数
FP,1,1e1,2e1,5e1,1e2,2e2,5e2 !
根据疲劳曲线输入S-N数据
FP,7,1e3,2e3,5e3,1e4,2e4,5e4
FP,13,1e5,2e5,5e5,1e6,,
FP,19,,
FP,21,4000,2828,1897,1414,1069,724
FP,27,572,441,331,262,214,159
FP,33,138,114,93.1,86.2,,
FP,39,,
!
****** 水压试验循环******
fs,4760,1,1,1,0,0,0,0,0,0 !
储存节点4760对应其第一载荷的应力
set,1,last !
读入第一载荷步数据
FSNODE,4760,1,2 !
储存节点4760对应其第二载荷的应力
fe,1,n2,p3/p1 !
设定事件循环次数及载荷比例系数
!
****** 最高/最低压力循环******
fs,4760,2,1,1,0,0,0,0,0,0 !
储存节点4760对应其第一载荷的应力
set,1,last !
读入第一载荷步数据
FSNODE,4760,2,2 !
储存节点4760对应其第二载荷的应力
FE,2,n1,1, !
设定事件循环次数及载荷比例系数
FTCALC,1 !
进行疲劳计算(并记录使用系数)
fini
!
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
计算结果如下:
PERFORMFATIGUECALCULATIONATLOCATION 1 NODE 0
***POST1FATIGUECALCULATION***
LOCATION 1 NODE 4760
事件1:
****** 水压试验循环******
EVENT/LOADS 1 1 AND 1