锚杆巷道ansys命令流

fini/clear/title,jiamaogan4/filename,jiamaogan4,0!指定工作名称/units,si!国际单位制!参数设置/prep7*set,hezai,9e5!荷载*set,mglength,4!锚杆长度4m*set,mgshu,54!锚杆数*set,mgdanyuan,16!锚杆单元数16*set,mgbshu,18!定义材料et,1,solid45et,2,link8!定义材料属性!岩石mp,ex,1,6.4e9mp,prxy,1,0.25mp,dens,1,2410tb,dp,1TBDATA,,1.55e6,44,0!锚杆mp,ex,2,230e9mp,prxy,2,0.26mp,dens,2,5800r,1r,2,3.14*(0.02/2)**2,2.6e-3!矩形5m建体k,1,-20,0,0k,2,-3,0,0k,3,3,0,0k,4,20,0,0k,5,-20,10,0k,6,-3,10,0k,7,3,10,0k,8,20,10,0k,9,-20,18,0k,10,-3,18,0k,11,3,18,0k,12,20,18,0k,13,-20,20,0k,14,-3,20,0k,15,3,20,0k,16,20,20,0k,17,-20,24.8,0k,18,-3,24.8,0k,19,3,24.8,0k,20,20,24.8,0k,21,-20,27,0k,22,-3,27,0k,23,3,27,0k,24,20,27,0k,25,-20,32,0k,26,-3,32,0k,27,3,32,0k,28,20,32,0k,29,-20,45,0k,30,-3,45,0k,31,3,45,0k,32,20,45,0a,1,2,6,5a,2,3,7,6a,3,4,8,7a,5,6,10,9a,6,7,11,10a,7,8,12,11a,9,10,14,13a,10,11,15,14a,11,12,16,15a,13,14,18,17a,14,15,19,18a,15,16,20,19a,17,18,22,21a,18,19,23,22a,19,20,24,23a,21,22,26,25a,22,23,27,26a,23,24,28,27a,25,26,30,29a,26,27,31,30a,27,28,32,31!dragareastovolumsallsel,allk,100,-20,0,-5l,1,100vdrag,all,,,,,,53ldele,53nummrg,all,,,,lownumcmp,alllsel,s,,,27lsel,a,,,25lsel,a,,,28lsel,a,,,30lesize,all,,,10lsel,s,,,5lsel,a,,,7lsel,a,,,15lsel,a,,,22lsel,a,,,29lsel,a,,,36lsel,a,,,43lsel,a,,,50lesize,all,,,12lsel,s,,,4lsel,a,,,2lsel,a,,,6lsel,a,,,9lesize,all,,,10lsel,s,,,13lsel,a,,,11lsel,a,,,14lsel,a,,,16lesize,all,,,8lsel,s,,,20lsel,a,,,18lsel,a,,,21lsel,a,,,23lesize,all,,,4lsel,s,,,34lsel,a,,,32lsel,a,,,35lsel,a,,,37lesize,all,,,4lsel,s,,,41lsel,a,,,39lsel,a,,,42lsel,a,,,44lesize,all,,,5lsel,s,,,48lsel,a,,,46lsel,a,,,49lsel,a,,,51lesize,all,,,13allsel,alllsel,s,,,54lesize,all,,,5!体厚的单元划分allsel,allvatt,1,1,1vsweep,allallsel,all!定义锚杆k,1001,-3,20.45,-1k,1002,-3,21.35,-1k,1003,-3,22.25,-1k,1004,-3,23.15,-1k,1005,-3,24.05,-1k,1006,-3,24.8,-1k,1007,-2.3,24.8,-1k,1008,-1.4,24.8,-1k,1009,-0.5,24.8,-1k,1010,0.5,24.8,-1k,1011,1.4,24.8,-1k,1012,2.3,24.8,-1k,1013,3,24.8,-1k,1014,3,24.05,-1k,1015,3,23.15,-1k,1016,3,22.25,-1k,1017,3,21.35,-1k,1018,3,20.45,-1k,1019,-7,20.45,-1k,1020,-7,21.35,-1k,1021,-7,22.25,-1k,1022,-7,23.15,-1k,1023,-7,24.05,-1k,1024,-7,24.8,-1k,1025,-2.3,28.8,-1k,1026,-1.4,28.8,-1k,1027,-0.5,28.8,-1k,1028,0.5,28.8,-1k,1029,1.4,28.8,-1k,1030,2.3,28.8,-1k,1031,7,24.8,-1k,1032,7,24.05,-1k,1033,7,23.15,-1k,1034,7,22.25,-1k,1035,7,21.35,-1k,1036,7,20.45,-1*do,ii,1001,1000+mgbshu,1l,ii,ii+mgbshu*enddoallsellgen,4,137,154,1,,,-1lsel,s,line,,137,154,1!需要改动cm,line0,linelsel,s,,,line0latt,2,2,2lesize,all,,,mgdanyuanlmesh,allallsel,all!锚杆节点耦合lsel,s,,,line0nsll,s,1cm,cmljnod,node!定义锚杆节点名称为cmljnod*get,max1,node,0,count!节点的总数为MAX1*dim,ojd,,max1!ojd最大下标*dim,jd,,max1!jd最大下标*get,nod1,node,0,num,min!nod1为最小节点号ojd(1)=nod1!ojd(1)为最小节点号*do,i,2,max1ojd(i)=ndnext(ojd(i-1))!当前节点号存入ojd(i)*enddoallsel,allnsel,allcmsel,u,cmljnod*do,i,1,max1nod1=ojd(i)j=nnear(nod1)nsel,u,,,jjd(i)=j!将与锚杆节点最近的节点存入jd(i)*enddonsel,allji=1*do,i,1,max1cp,ji,ux,ojd(i),jd(i)cp,ji+1,uy,ojd(i),jd(i)cp,ji+2,uz,ojd(i),jd(i)ji=ji+3*enddoallsel,allfinish!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!1/config,nproc,2!设置处理器数/soluallselantype,staticasel,s,loc,x,-20asel,a,loc,x,20da,all,ux,0allsel!以下约束地面y方向的约束asel,s,loc,y,0da,all,uy,0allsel!施加重力加速度acel,,9.8asel,s,area,,86asel,a,area,,90asel,a,area,,93SFA,all,,pres,hezaiesel,s,mat,,2ekill,allallsel,all!设置非线性选项!nropt,Modi!用改变的牛顿-拉普拉斯方程求解!nlgeom,on!打开大变形选项!autots,off!打开自动时间步!lnsrch,on!自动线性搜索!pred,on!打开预测求解器!outres,all,10!TIME,1!kbc,0!NSUBST,10!allsel,allsolve/soluallsel,allvsel,s,,,11eslv,s,1ekill,allesel,s,mat,,2ealive,allallsel,allsolve!savefini