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