算弯矩的命令流

finish/CLEAR/Filname,Ⅲ,1!指定工作名称/title,LANXI!指定图形显示标题/units,si!国际单位制/PREP7!**********************尺寸参数的设定***********************h=1.5!钢轨间距!轨道板+自密实混凝土层参数kjjj=0.625!扣件间距kjn=9!扣件对数bw=2.5!板宽bl=5.6!板长bn=3!板数bf=0.07!板缝bh=0.21!板厚bbkjj=0.3!板边到第一组扣件的距离!底座板参数dzl=bl*bn+(bn-1)*bfdzw=2.74dzh=0.3!找平层参数zpl=bl*bn+(bn-1)*bfzpw=3.3zph=0.12!***************定义单元类型并设置材料属性***********************et,1,shell63et,2,beam3et,3,combin14keyopt,3,2,3mp,ex,1,3.65e10!轨道板mp,ey,1,3.65e10mp,ez,1,3.65e10mp,prxy,1,0.2mp,pryz,1,0.2mp,prxz,1,0.2mp,alpx,1,1e-5mp,alpy,1,1e-5mp,alpz,1,1e-5mp,dens,1,2.5e3mp,ex,2,3.25e10!底座板mp,ey,2,3.25e10mp,ez,2,3.25e10mp,prxy,2,0.2mp,pryz,2,0.2mp,prxz,2,0.2mp,alpx,2,1e-5mp,alpy,2,1e-5mp,alpz,2,1e-5mp,dens,2,2.5e3mp,ex,3,3e10!找平层mp,ey,3,3e10mp,ez,3,3e10mp,prxy,3,0.2mp,pryz,3,0.2mp,prxz,3,0.2mp,alpx,3,1e-5mp,alpy,3,1e-5mp,alpz,3,1e-5mp,dens,3,2.5e3mp,ex,4,2.06e11!钢轨材料属性mp,alpx,4,1.2e-5mp,gxy,4,7.9e10mp,dens,4,7.8e3mp,prxy,4,0.3r,1,bhr,2,dzhr,3,zphr,4,77.45e-4,3217e-8,0.176r,5,46.4156e7/1479!橡胶垫层刚度r,6,2.5e7!扣件刚度r,7,14e9/1305!自密实层与底座板的连接弹簧刚度r,8,1e20!路基弹簧刚度!*****************建立模型*****************!*****************建立轨道板***************blc4,-bbkjj,-(bw-h)/2,bl,bwagen,3,all,,,bl+bf!*****************建立底座板***************wpoffs,,,-(bh+dzh)/2blc4,-bbkjj,-(dzw-h)/2,dzl,dzw!*****************建立找平层***************wpoffs,,,-(dzh+zph)/2blc4,-bbkjj,-(zpw-h)/2,zpl,zpw!*****************切分轨道板**********WPCSYS,-1asel,allwprota,0,0,90asbw,all*do,i,1,kjn-1,1wpoffs,,,kjjjasbw,all*enddowpoffs,,,bbkjjasbw,allwpoffs,,,bfasbw,allwpoffs,,,bbkjjasbw,all*do,j,1,kjn-1,1wpoffs,,,kjjjasbw,all*enddowpoffs,,,bbkjjasbw,allwpoffs,,,bfasbw,allwpoffs,,,bbkjjasbw,all*do,k,1,kjn-1,1wpoffs,,,kjjjasbw,all*enddowprota,0,90asbw,allwpoffs,,,(bw-h)/2asbw,allwpoffs,,,(dzw-bw)/2asbw,allwpoffs,,,-dzwasbw,allwpoffs,,,(dzw-bw)/2asbw,allwpoffs,,,(bw-h)/2asbw,all!划分轨道板单元wpcsys,-1lsel,s,loc,z,0lsel,r,length,,bbkjjlesize,all,,,2lsel,s,loc,z,0lsel,r,length,,(bw-h)/2lesize,all,,,2lsel,s,loc,z,0lsel,r,length,,kjjjlesize,all,,,3lsel,s,loc,z,0lsel,r,length,,hlesize,all,,,10asel,s,loc,z,0type,1mat,1real,1mshape,0,3dmshkey,1amesh,all!*****************生成底座板单元**********lsel,s,loc,z,-(bh+dzh)/2lsel,r,length,,bbkjjlesize,all,,,2lsel,s,loc,z,-(bh+dzh)/2lsel,r,length,,(bw-h)/2lesize,all,,,2lsel,s,loc,z,-(bh+dzh)/2lsel,r,length,,kjjjlesize,all,,,3lsel,s,loc,z,-(bh+dzh)/2lsel,r,length,,hlesize,all,,,10lsel,s,loc,z,-(bh+dzh)/2lsel,r,length,,(dzw-bw)/2lesize,all,,,1lsel,s,loc,z,-(bh+dzh)/2lsel,r,length,,bflesize,all,,,1asel,s,loc,z,-(bh+dzh)/2type,1mat,2real,2mshape,0,3dmshkey,1amesh,all!*****************生成找平层单元**********lsel,s,loc,z,-(bh+zph)/2-dzhlsel,r,length,,bbkjjlesize,all,,,2lsel,s,loc,z,-(bh+zph)/2-dzhlsel,r,length,,(bw-h)/2lesize,all,,,2lsel,s,loc,z,-(bh+zph)/2-dzhlsel,r,length,,kjjjlesize,all,,,3lsel,s,loc,z,-(bh+zph)/2-dzhlsel,r,length,,hlesize,all,,,10lsel,s,loc,z,-(bh+zph)/2-dzhlsel,r,length,,(dzw-bw)/2lesize,all,,,1lsel,s,loc,z,-(bh+zph)/2-dzhlsel,r,length,,(zpw-dzw)/2lesize,all,,,2lsel,s,loc,z,-(bh+zph)/2-dzhlsel,r,length,,bflesize,all,,,1asel,s,loc,z,-(bh+zph)/2-dzhtype,1mat,3real,3mshape,0,3dmshkey,1amesh,all!*****************轨道板与底座板连接************nsel,s,loc,z,0nsel,a,loc,z,-(bh+dzh)/2type,3real,7tshap,lineeintf,0.0001,,low,,,,-(bh+dzh)/2!*****************橡胶垫层************nsel,s,loc,z,-(bh+dzh)/2nsel,a,loc,z,-(bh+zph)/2-dzhtype,3real,5tshap,line!确定目标面的几何形状eintf,0.0001,,low,,,,-(dzh+zph)/2!*****************钢轨************wpcsys,-1!右股钢轨节点alls*do,m,1,bn,1*do,l,1,kjn,1allsnsel,s,loc,z,0nsel,r,loc,y,0nsel,r,loc,x,(m-1)*(bl+bf)+kjjj*(l-1)ngen,2,4999+l+10000*(m-1),all,,,,,0.04+bh/2*enddo*enddoallswpcsys,-1alls!左股钢轨节点*do,n,1,bn,1*do,o,1,kjn,1allsnsel,s,loc,z,0nsel,r,loc,y,hnsel,r,loc,x,(n-1)*(bl+bf)+kjjj*(o-1)ngen,2,5999+o+10000*(n-1),all,,,,,0.04+bh/2*enddo*enddoallsnsel,s,loc,z,0.04+bh/2!连成钢轨type,2mat,4real,4tshap,lineeintf,0.0001,,low,,kjjjallsnsel,s,loc,z,0.04+bh/2type,2mat,4real,4tshap,lineeintf,0.0001,,low,,2*bbkjj+bfallsnsel,s,loc,y,0!右股钢轨扣件nsel,r,loc,z,0,0.04+bh/2type,3real,6tshap,lineeintf,0.0001,,high,,,,0.04+bh/2allsnsel,s,loc,y,h!左股钢轨扣件nsel,r,loc,z,0,0.04+bh/2type,3real,6tshap,lineeintf,0.0001,,high,,,,0.04+bh/2alls!*******************路基弹簧*******************nsel,s,loc,z,-(bh+zph)/2-dzhngen,2,10000,all,,,,,-0.2nsel,s,loc,z,-(bh+zph)/2-dzh!路基弹簧nsel,a,loc,z,-(bh+zph)/2-dzh-0.2type,3real,8tshap,lineeintf,0.0001,,high,,,,0.2alls!*******************约束*******************allsnsel,s,loc,z,-(bh+zph)/2-dzh-0.2d,all,allallseplotNSEL,S,,,NODE((3*bl/2+bf),0,0.04+bh/2)F,ALL,fz,-300000/2ALLSEL,ALLNSEL,S,,,NODE((3*bl/2+bf),h,0.04+bh/2)F,ALL,fz,-300000/2ALLSEL,ALLLOCAL,11,0,,h/2,GD_D/2ALLSEL,ALLEMODIF,ALL,ESYS,11ALLSEL,ALLNUMCMP,NODEFINISH/solusolvefinish/post1etable,wjx,smisc,4etable,wjy,smisc,5