马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
<P>我用ansys二维接触模拟抗浮锚杆<br>怎么锚干一拔就出来了<br>我不知道是接触的问题还是没有加初始应力<br>下面是我的命令流<br>敬请高手指点,我是哪里出问题了?<br><br>/prep7 <br>et,1,solid45<br>et,2,plane42<br>mp,ex,1,2.5e10 !锚杆的弹性模量<br>mp,nuxy,1,0.2 !锚杆的泊松比<br>mp,dens,1,2500 !锚杆的密度<br>mp,ex,2,2.5e8 !土的弹性模量<br>mp,nuxy,2,0.3 !土的泊松比<br>mp,dens,2,2000 !土的密度 <br>tb,dp,2 <br>tbdata,1,1900,32,30 !粘聚力c为19,摩擦角为32度,膨胀角为30</P>
<P>RECTNG,0,0.5,0,14, !面1<br>RECTNG,-8,8.5,-5,12 !面2<br>ASBA,2,1 <br>RECTNG,0,0.5,0,12.5, !面1 <br>/pnum,area,1<br>/pnum,line,1<br>numcmp,all<br>allsel<br>aplot</P>
<P>lsel,s,,,2,4,2 !划分面1 <br>lesize,all,,,1 <br>lsel,s,,,3,7,2 <br>lesize,all,0.5 <br>mat,1<br>amesh,1 </P>
<P>lsel,s,,,6,8,2 !划分面2<br>lesize,all,0.5<br>lsel,s,,,9 <br>lesize,all,0.5 <br>lsel,s,,,10 <br>lesize,all,0.5 <br>lsel,s,,,11 <br>lesize,all,0.5 <br>lsel,s,,,12 <br>lesize,all,0.5 <br>lsel,s,,,1<br>lesize,all,,,2 <br>lsel,s,,,5<br>lesize,all,0.5 <br>mat,2<br>amesh,2</P>
<P>!施加初始应力</P>
<P>!接触单元设置<br>allsel<br>et,3,169 !定义接触单元<br>et,4,171<br>asel,s,,,1<br>nsla,s,1<br>nsel,r,loc,x,0.5<br>r,10<br>real,10<br>mp,mu,2,0.3 !摩擦系数<br>type,3<br>esurf<br>allsel<br>asel,s,,,2</P>
<P>nsla,s,1<br>nsel,r,loc,x,0.5<br>type,4<br>esurf<br>allsel,all</P>
<P>asel,s,,,1<br>nsla,s,1<br>nsel,r,loc,x,0<br>r,10<br>real,10<br>mp,mu,2,0.8<br>type,3<br>esurf<br>allsel<br>asel,s,,,2</P>
<P>nsla,s,1<br>nsel,r,loc,x,0<br>type,4<br>esurf<br>allsel,all</P>
<P>!底端接触<br>asel,s,,,1<br>nsla,s,1<br>nsel,r,loc,y,0<br>r,10<br>real,10<br>mp,mu,2,0.3 !摩擦系数<br>type,3<br>esurf<br>allsel<br>asel,s,,,2</P>
<P>nsla,s,1<br>nsel,r,loc,y,0<br>type,4<br>esurf<br>allsel,all</P>
<P><br>allsel !施加约束及荷载<br>/solu<br>asel,s,,,1<br>nsla,s,1<br>sfl,2,pres,20000<br>!dl,2,1,uy,0.8<br>asel,s,,,2<br>nsla,s,1<br>dl,8,2,ux<br>dl,6,2,ux<br>dl,5,2,all,all<br>asel,s,,,1</P>
<P>!sfl,9,pres,-200000<br>!sfl,10,pres,-20000</P>
<P><br>alls !选择所有实体<br>antype,static<br>!nlgeom,on<br>nsub,100,10000,100<br>outres,all,all<br>allsel<br>solve<br>finish<br></P>
[此贴子已经被作者于2006-4-27 21:15:22编辑过]
|