ANSYS地铁车站横断面内力命令流.doc

finish$/clear$/prep7!定义单元类型、实常数、材料属性、截面et,1,beam3et,2,link10et,3,combin14sectype,1,beam,rect,0secoffset,centsecdata, 1.1, 1, 0,0,0,0,0,0,0,0sectype,2,beam,rect,0secoffset,centsecdata, 0.4, 1, 0,0,0,0,0,0,0,0sectype,3,beam,rect, ,0secoffset,centsecdata, 1, 1, 0,0,0,0,0,0,0,0sectype,4, beam,rect,0secoffset,centsecdata, 0.8, 1, 0,0,0,0,0,0,0,0sectype,5,beam,rect,0secoffset,centsecdata, 0.23375, 1, 0,0,0,0,0,0,0,0r,1, 1*1.1, 1*1.1*1.1*1.1/12, 1.1, , , ,r,2, 1*0.4, 1*0.6*0.6*0.6/12, 0.4, , , ,r,3, 1*1.0, 1*0.8*0.8*0.8/12, 1.0, , , ,r,4, 1*0.8, 1*0.8*0.8*0.8/12, 0.8, , , ,r,5, 0.23375*1, 0.23375*1*1*1/12, 1, , , ,r,6, 160*1*1*1e6,r,7, 61.63*1*1*1e6,mp,ex,1,30e9mp,prxy,1,0.2mp,dens,1,2.5e3mp,ex,2,32.5e9mp,prxy,2,0.2mp,dens,2,2.5e3!建立几何模型k,1, 0,0k,2, 7.7,0k,3, 14.8,0k,4, 22.5,0k,5, 0,7.275k,6, 7.7,7.275k,7, 14.8,7.275k,8, 22.5,7.275k,9, 0,14.05k,10, 7.7,14.05k,11, 14.8,14.05k,12, 22.5,14.05l,1,2l,2,3l,3,4l,5,6l,6,7l,7,8l,9,10l,10,11l,11,12l,1,5l,5,9l,2,6l,6,10l,3,7l,7,11l,4,8l,8,12!单元网格划分lsel,s,line,1,2,1lesize,all,1latt,1,1,1,1lmesh,alllsel,s,line,2,3,1lesize,all,1latt,1,1,1,1lmesh,alllsel,s,line,4,5,1lesize,all,1latt,1,2,1,2lmesh,alllsel,s,line,5,6,1lesize,all,1latt,1,2,1,2lmesh,alllsel,s,line,7,8,1lesize,all,1latt,1,3,1,3lmesh,alllsel,s,line,8,9,1lesize,all,1latt,1,3,1,3lmesh,alllsel,s,line,10,11,1lesize,all,1latt,1,4,1,4lmesh,alllsel,s,line,16,17,1lesize,all,1latt,1,4,1,4lmesh,alllsel,s,line,12,13,1lesize,all,1latt,2,5,1,5lmesh,alllsel,s,line,14,15,1lesize,all,1latt,2,5,1,5lmesh,all!创建弹簧单元nplot/pnum,node,1/replottype,3real,6nsel,s,loc,y,0ngen,2,500,all,-1*do,i,1,25,1e,i,i+500*enddo!施加边界条件allsel,allnsel,s,loc,y,-1d,all,ux,0,uyallsel,alld,2,uy,0d,10,uy,0d,14,ux,0!施加荷载/soluq1=134.5e3q2=16.6e3q3=268.8e3e1=97.7e3e2=403.5e3allsel,allesel,s,real,3sfbeam,all,1,pres,q1allsel,allesel,s,real,1sfbeam,all,1,pres,-q3allsel,allesel,s,real,2sfbeam,all,1,pres,q2allsel,allesel,s,real,4nsle,snsel,r,loc,x,0esln,r,1sfbeam,pres,*get,nelm,elem,0,count*get,emin,elem,0,num,min*dim,p,array,nelm,3,1k=emin*do,i,1,nelmp(i,1)=kk=elnext(k)*enddon=1*do,i,n,nelmp(i,2)=e2-(e2-e1)*(i-n)*1/14.05p(i,3)=e2-(e2-e1)*(i-n+1)*1/14.05*enddo*do,i,1,nelmsfbeam,p(i,1),1,pres,p(i,2),p(i,3)*enddoallsel,allesel,s,real,4nsle,snsel,r,loc,x,22.5esln,r,1*get,nelm,elem,0,count*get,emin,elem,0,num,min*dim,p,array,nelm,3,1k=emin*do,i,1,nelmp(i,1)=kk=elnext(k)*enddon=1*do,i,n,nelmp(i,2)=e2-(e2-e1)*(i-n)*1/14.05p(i,3)=e2-(e2-e1)*(i-n+1)*1/14.05*enddo*do,i,1,nelmsfbeam,p(i,1),1,pres,-p(i,2),-p(i,3)*enddoallsel,allacel,9.8 !自重应该加9.8*1.485finish/sol/status,solusolvefinish/post1etable,ni,smisc,1etable,nj,smisc,7etable,qi,smisc,2etable,qj,smisc,8etable,mi,smisc,6etable,mj,smisc,12plls,ni,nj,1,0plls,qi,qj,1,0plls,mi,mj,-1,0esel,u,type,2！显示变形图pldisp,1finish2.标准组合命令流finish/clear/title,tiyuzhongxinnanzhan!定义单元类型、实常数、材料属性、截面/prep7et,1,beam3et,2,link10keyopt,2,3,1et,3,combin14sectype,1,beam,rect,0secoffset,centsecdata, 1.1, 1, 0,0,0,0,0,0,0,0sectype,2,beam,rect,0secoffset,centsecdata, 0.4, 1, 0,0,0,0,0,0,0,0sectype,3,beam,rect, ,0secoffset,centsecdata, 1, 1, 0,0,0,0,0,0,0,0sectype,4, beam,rect,0secoffset,centsecdata, 0.8, 1, 0,0,0,0,0,0,0,0sectype,5,beam,rect,0secoffset,centsecdata, 0.23375, 1, 0,0,0,0,0,0,0,0r,1, 1*1.1, 1*1.1*1.1*1.1/12, 1.1, , , ,！底板r,2, 1*0.4, 1*0.6*0.6*0.6/12, 0.4, , , ,！中板r,3, 1*1.0, 1*0.8*0.8*0.8/12, 1.0, , , ,！顶板r,4, 1*0.8, 1*0.8*0.8*0.8/12, 0.8, , , ,！侧墙r,5, 0.23375*1, 0.23375*1*1*1/12, 1, , , ,！中柱r,6, 160*1*1*1e6,r,7, 61.63*1*1*1e6,mp,ex,1,30e9mp,prxy,1,0.2mp,dens,1,2.5e3mp,ex,2,32.5e9mp,prxy,2,0.2mp,dens,2,2.5e3!建立几何模型k,1, 0,0k,2, 7.7,0k,3, 14.8,0k,4, 22.5,0k,5, 0,7.275k,6, 7.7,7.275k,7, 14.8,7.275k,8, 22.5,7.275k,9, 0,14.05k,10, 7.7,14.05k,11, 14.8,14.05k,12, 22.5,14.05l,1,2l,2,3l,3,4l,5,6l,6,7l,7,8l,9,10l,10,11l,11,12l,1,5l,5,9l,2,6l,6,10l,3,7l,7,11l,4,8l,8,12!单元网格划分lsel,s,line,1,2,1lesize,all,1latt,1,1,1,1lmesh,alllsel,s,line,2,3,1lesize,all,1latt,1,1,1,1lmesh,alllsel,s,line,4,5,1lesize,all,1latt,1,2,1,2lmesh,alllsel,s,line,5,6,1lesize,all,1latt,1,2,1,2lmesh,alllsel,s,line,7,8,1lesize,all,1latt,1,3,1,3lmesh,alllsel,s,line,8,9,1lesize,all,1latt,1,3,1,3lmesh,alllsel,s,line,10,11,1lesize,all,1latt,1,4,1,4lmesh,alllsel,s,line,16,17,1lesize,all,1latt,1,4,1,4lmesh,alllsel,s,line,12,13,1lesize,all,1latt,2,5,1,5lmesh,alllsel,s,line,14,15,1lesize,all,1latt,2,5,1,5lmesh,all!创建弹簧单元nplot/pnum,node,1/replottype,3real,6nsel,s,loc,y,0ngen,2,500,all,-1*do,i,1,25,1e,i,i+500*enddo!施加边界条件allsel,allnsel,s,loc,y,-1d,all,ux,0,uyallsel,alld,2,uy,0d,10,uy,0d,14,ux,0!施加荷载/soluq1=94.6e3q2=12.0e3q3=181e3e1=67.7e3e2=273.6e3allsel,allesel,s,real,3sfbeam,all,1,pres,q1allsel,allesel,s,real,1sfbeam,all,1,pres,-q3allsel,allesel,s,real,2sfbeam,all,1,pres,q2allsel,allesel,s,real,4nsle,snsel,r,loc,x,0esln,r,1*get,nelm,elem,0,count*get,emin,elem,0,num,min*dim,p,array,nelm,3,1k=emin*do,i,1,nelmp(i,1)=kk=elnext(k)*enddon=1*do,i,n,nelmp(i,2)=e2-(e2-e1)*(i-n)*1/14.05p(i,3)=e2-(e2-e1)*(i-n+1)*1/14.05*enddo*do,i,1,nelmsfbeam,p(i,1),1,pres,p(i,2),p(i,3)*enddoallsel,allesel,s,real,4nsle,snsel,r,loc,x,22.5esln,r,1*get,nelm,elem,0,count*get,emin,elem,0,num,min*dim,p,array,nelm,3,1k=emin*do,i,1,nelmp(i,1)=kk=elnext(k)*enddon=1*do,i,n,nelmp(i,2)=e2-(e2-e1)*(i-n)*1/14.05p(i,3)=e2-(e2-e1)*(i-n+1)*1/14.05*enddo*do,i,1,nelmsfbeam,p(i,1),1,pres,-p(i,2),-p(i,3)*enddoallsel,allacel,9.8finish/sol/status,solusolvefinish/post1etable,ni,smisc,1etable,nj,smisc,7etable,qi,smisc,2etable,qj,smisc,8etable,mi,smisc,6etable,mj,smisc,12plls,ni,nj,1,0plls,qi,qj,1,0plls,mi,mj,-1,0esel,u,type,2！显示变形图pldisp,1finish3.纵梁命令流/prep7et,1,beam3et,2,link10keyopt,2,3,1r,1,1.2*2,1.2*2*2*2/12,2 ！顶纵梁r,2,1.2*1.2,1.2*1.2*1.2*1.2/12,1.2 ！中纵梁r,3,1.2*2.5,1.2*2.5*2.5*2.5/12,2.5 ！底纵梁r,4,0.8*1.2,0.8*1.2*1.2*1.2/12,1.2 ！中柱r,5,1*1 ！侧墙mp,ex,1,3e10mp,prxy,1,0.2mp,dens,1,2500mp,ex,2,3.35e10mp,prxy,2,0.2mp,dens,2,2500mp,ex,3,1.6e8*do,i,1,41,1n,i,i-1,18.5n,i+41,i-1,8.5n,i+