| 
除去蠕变,这个模型的结果可靠性是不错的。作了一系列接触问题,通过试验验证符合的很好。 
模型解释:(1)一个弹性结构受压(接触)变形,到发生塑性变形。(2)拿开压缩板,结构回弹,但不会回到原始位置。(3)这时计算蠕变,释放掉应力。(4)再压弹性结构到开始压缩位置。比较这四步的接触力。结果:第二,三步当然没有接触力,(若没有应力释放,第一、第四步接触力应一样,)有了应力释放,第四步接触力比第一步减小。
 这个模型中的蠕变没用太好。用的是隐式6号蠕变方程,蠕变是时间和应力的函数,参数是乱定的(应力释放太快)。
 想请教有关蠕变方面的资料,尤其是材料蠕变方程选用及参数方面的资料。
 /prep7
 !------------CuSn8----------
 ET,1,182,,,3
 mp,ex,1,115e9
 mp,prxy,1,0.3
 r,1,0.3
 TB,BKIN,1
 TBDATA,1,470E6,0
 tm=100
 *SET,C1,1.5625E-14!ASSIGN VALUE
 *SET,C2,1.5!ASSIGN VALUE
 *SET,C3,!ASSIGN VALUE
 *SET,C4,0!ASSIGN VALUE
 TB,CREEP,1,,,6!ACTIVATE DATA TABLE
 TBDATA,1,C1,C2,C3,C4!DEFINE DATA FOR TABLE
 !-----------contact-----------------
 ET,9,169
 ET,10,171
 R,9,,,0.1,0.1,,
 !RMORE,,,1.0E20,0.0,1.0,
 !RMORE,0.0,0,1.0,0,0,0.5
 !RMORE,,,1.0,0.0
 MP,MU,9,0.0
 !----------------geometry
 k,,2
 k,,2,0.2
 k,,,0.2
 k,,-0.2
 k,,-0.2,1.2
 k,,,1
 k,,2,1.2
 k,,1,1
 k,,1.25,1
 k,,2,1
 L,8,9,
 k,,1.5,1.2
 k,,1.75,1.45
 L, 1, 2
 L, 1, 4
 L, 4, 5
 L, 5, 11
 larc,7,12,11,0.25
 larc,11,12,7,0.25
 L, 7, 10
 L, 10, 9
 L, 8, 6
 L, 6, 3
 L, 3, 2
 LFILLT,11,10,0.3, ,
 !*
 LFILLT,4,5,0.5, ,
 !*
 LFILLT,11,12,0.3, ,
 !*
 LFILLT,4,3,0.5, ,
 FLST,2,16,4
 FITEM,2,12
 FITEM,2,15
 FITEM,2,11
 FITEM,2,13
 FITEM,2,10
 FITEM,2,1
 FITEM,2,9
 FITEM,2,8
 FITEM,2,7
 FITEM,2,6
 FITEM,2,5
 FITEM,2,14
 FITEM,2,4
 FITEM,2,16
 FITEM,2,3
 FITEM,2,2
 AL,P51X
 rect,1,3,1.45+0.001,1.5
 type,1
 mat,1
 esize,0.05
 amesh,all
 !---------contact------------
 alls
 type,10
 mat,9
 real,9
 lsel,s,,,6,7
 nsll,s,1
 esln,s,0
 esurf,all
 type,9
 mat,9
 real,9
 lsel,s,,,17
 nsll,s,1
 esln,s,0
 esurf,all
 !------boundary
 lsel,s,,,3
 nsll,,1
 d,all,ux
 d,all,uy
 lsel,s,,,19
 nsll,,1
 cp,11,uy,all
 cplgen,11,ux
 *get,nmin,node,,num,min
 d,nmin,ux
 ksel,s,,,10
 nslk
 *get,ndis,node,,num,min
 fini
 /solu
 antype,static
 nlgeom,on
 autots,on
 alls
 save
 rate,off
 time,1e-8
 d,nmin,uy,-0.3
 nsub,20
 outres,all,all
 solve
 *get,rf1,node,nmin,rf,fy
 *get,dis1,node,ndis,u,y
 time,2e-8
 d,nmin,uy,0.0
 nsub,20
 outres,all,all
 solve
 *get,rf2,node,nmin,rf,fy
 *get,dis2,node,ndis,u,y
 !BFUNIF,TEMP,90
 rate,on
 TIME,tm
 !NSUBST,10
 OUTPR,BASIC,10 ! PRINT BASIC SOLUTION FOR EVERY 10TH SUBSTEP
 OUTRES,ESOL,1 ! STORE ELEMENT SOLUTION FOR EVERY SUBSTEP
 SOLVE
 *get,rf3,node,nmin,rf,fy
 *get,dis3,node,ndis,u,y
 rate,off
 time,tm+1e-8
 d,nmin,uy,-0.3
 nsub,20
 outres,all,all
 solve
 *get,rf4,node,nmin,rf,fy
 *get,dis4,node,ndis,u,y
 /EOF
 time,11
 d,nmin,uy,-0.0
 nsub,20
 outres,all,all
 solve
 *get,rf11,node,nmin,rf,fy
 *get,dis11,node,ndis,u,y
 /eof
 fini
 /post1
 *get,rf2,node,nmin,rf,fy
 fini
 /eof
 
   |