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

1、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,0SECO

2、FFSET,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

3、*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,2

4、l,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,

5、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,

6、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

7、,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)=e

8、2-(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)=

9、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

10、,-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,

11、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*

12、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

13、.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,

14、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,alllse

15、l,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.6

16、e3q2=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,

17、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

18、,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

19、,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+