ANSYS地铁车站横断面内力命令流New.docx
- 文档编号:4367363
- 上传时间:2022-12-01
- 格式:DOCX
- 页数:17
- 大小:15.95KB
ANSYS地铁车站横断面内力命令流New.docx
《ANSYS地铁车站横断面内力命令流New.docx》由会员分享,可在线阅读,更多相关《ANSYS地铁车站横断面内力命令流New.docx(17页珍藏版)》请在冰豆网上搜索。
ANSYS地铁车站横断面内力命令流New
ANSYS地铁车站横断面内力命令流
ANSYS地铁车站横断面内力命令流
finish$/clear$/prep7
!
①定义单元类型、实常数、材料属性、截面
ET,1,BEAM3
ET,2,LINK10
et,3,combin14
SECTYPE,1,BEAM,RECT,,0
SECOFFSET,CENT
SECDATA,1.1,1,0,0,0,0,0,0,0,0
SECTYPE,2,BEAM,RECT,,0
SECOFFSET,CENT
SECDATA,0.4,1,0,0,0,0,0,0,0,0
SECTYPE,3,BEAM,RECT,,0
SECOFFSET,CENT
SECDATA,1,1,0,0,0,0,0,0,0,0
SECTYPE,4,BEAM,RECT,,0
SECOFFSET,CENT
SECDATA,0.8,1,0,0,0,0,0,0,0,0
SECTYPE,5,BEAM,RECT,,0
SECOFFSET,CENT
SECDATA,0.23375,1,0,0,0,0,0,0,0,0
R,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,30e9
mp,prxy,1,0.2
mp,dens,1,2.5e3
mp,ex,2,32.5e9
mp,prxy,2,0.2
mp,dens,2,2.5e3
!
②建立几何模型
k,1,0,0
k,2,7.7,0
k,3,14.8,0
k,4,22.5,0
k,5,0,7.275
k,6,7.7,7.275
k,7,14.8,7.275
k,8,22.5,7.275
k,9,0,14.05
k,10,7.7,14.05
k,11,14.8,14.05
k,12,22.5,14.05
l,1,2
l,2,3
l,3,4
l,5,6
l,6,7
l,7,8
l,9,10
l,10,11
l,11,12
l,1,5
l,5,9
l,2,6
l,6,10
l,3,7
l,7,11
l,4,8
l,8,12
!
③单元网格划分
lsel,s,line,,1,2,1
lesize,all,1
LATT,1,1,1,,,,1
lmesh,all
lsel,s,line,,2,3,1
lesize,all,1
LATT,1,1,1,,,,1
lmesh,all
lsel,s,line,,4,5,1
lesize,all,1
LATT,1,2,1,,,,2
lmesh,all
lsel,s,line,,5,6,1
lesize,all,1
LATT,1,2,1,,,,2
lmesh,all
lsel,s,line,,7,8,1
lesize,all,1
LATT,1,3,1,,,,3
lmesh,all
lsel,s,line,,8,9,1
lesize,all,1
LATT,1,3,1,,,,3
lmesh,all
lsel,s,line,,10,11,1
lesize,all,1
LATT,1,4,1,,,,4
lmesh,all
lsel,s,line,,16,17,1
lesize,all,1
LATT,1,4,1,,,,4
lmesh,all
lsel,s,line,,12,13,1
lesize,all,1
LATT,2,5,1,,,,5
lmesh,all
lsel,s,line,,14,15,1
lesize,all,1
LATT,2,5,1,,,,5
lmesh,all
!
④创建弹簧单元
nplot
/pnum,node,1
/replot
type,3
real,6
nsel,s,loc,y,0
ngen,2,500,all,,,,-1
*do,i,1,25,1
e,i,i+500
*enddo
!
⑤施加边界条件
allsel,all
nsel,s,loc,y,-1
d,all,ux,0,,,,uy
allsel,all
d,2,uy,0
d,10,uy,0
d,14,ux,0
!
⑥施加荷载
/solu
q1=134.5e3
q2=16.6e3
q3=268.8e3
e1=97.7e3
e2=403.5e3
allsel,all
esel,s,real,,3
sfbeam,all,1,pres,q1
allsel,all
esel,s,real,,1
sfbeam,all,1,pres,-q3
allsel,all
esel,s,real,,2
sfbeam,all,1,pres,q2
allsel,all
esel,s,real,,4
nsle,s
nsel,r,loc,x,0
esln,r,1
sfbeam,pres,
*get,Nelm,elem,0,count
*get,Emin,elem,0,num,min
*dim,P,array,Nelm,3,1
k=Emin
*do,i,1,Nelm
P(i,1)=k
k=elnext(k)
*enddo
n=1
*do,i,n,Nelm
p(i,2)=e2-(e2-e1)*(i-n)*1/14.05
P(i,3)=e2-(e2-e1)*(i-n+1)*1/14.05
*enddo
*do,i,1,Nelm
sfbeam,P(i,1),1,pres,P(i,2),P(i,3)
*enddo
allsel,all
esel,s,real,,4
nsle,s
nsel,r,loc,x,22.5
esln,r,1
*get,Nelm,elem,0,count
*get,Emin,elem,0,num,min
*dim,P,array,Nelm,3,1
k=Emin
*do,i,1,Nelm
P(i,1)=k
k=elnext(k)
*enddo
n=1
*do,i,n,Nelm
p(i,2)=e2-(e2-e1)*(i-n)*1/14.05
P(i,3)=e2-(e2-e1)*(i-n+1)*1/14.05
*enddo
*do,i,1,Nelm
sfbeam,P(i,1),1,pres,-P(i,2),-P(i,3)
*enddo
allsel,all
acel,,9.8!
自重应该加9.8*1.485
finish
/SOL
/STATUS,SOLU
SOLVE
finish
/POST1
ETABLE,NI,SMISC,1
ETABLE,NJ,SMISC,7
ETABLE,QI,SMISC,2
ETABLE,QJ,SMISC,8
ETABLE,MI,SMISC,6
ETABLE,MJ,SMISC,12
PLLS,NI,NJ,1,0
PLLS,QI,QJ,1,0
PLLS,MI,MJ,-1,0
ESEL,U,TYPE,,2
!
显示变形图
pldisp,1
finish
2.标准组合命令流
FINISH/CLEAR/title,tiyuzhongxinnanzhan
!
①定义单元类型、实常数、材料属性、截面
/prep7
ET,1,BEAM3
ET,2,LINK10
keyopt,2,3,1
et,3,combin14
SECTYPE,1,BEAM,RECT,,0
SECOFFSET,CENT
SECDATA,1.1,1,0,0,0,0,0,0,0,0
SECTYPE,2,BEAM,RECT,,0
SECOFFSET,CENT
SECDATA,0.4,1,0,0,0,0,0,0,0,0
SECTYPE,3,BEAM,RECT,,0
SECOFFSET,CENT
SECDATA,1,1,0,0,0,0,0,0,0,0
SECTYPE,4,BEAM,RECT,,0
SECOFFSET,CENT
SECDATA,0.8,1,0,0,0,0,0,0,0,0
SECTYPE,5,BEAM,RECT,,0
SECOFFSET,CENT
SECDATA,0.23375,1,0,0,0,0,0,0,0,0
R,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,30e9
mp,prxy,1,0.2
mp,dens,1,2.5e3
mp,ex,2,32.5e9
mp,prxy,2,0.2
mp,dens,2,2.5e3
!
②建立几何模型
k,1,0,0
k,2,7.7,0
k,3,14.8,0
k,4,22.5,0
k,5,0,7.275
k,6,7.7,7.275
k,7,14.8,7.275
k,8,22.5,7.275
k,9,0,14.05
k,10,7.7,14.05
k,11,14.8,14.05
k,12,22.5,14.05
l,1,2
l,2,3
l,3,4
l,5,6
l,6,7
l,7,8
l,9,10
l,10,11
l,11,12
l,1,5
l,5,9
l,2,6
l,6,10
l,3,7
l,7,11
l,4,8
l,8,12
!
③单元网格划分
lsel,s,line,,1,2,1
lesize,all,1
LATT,1,1,1,,,,1
lmesh,all
lsel,s,line,,2,3,1
lesize,all,1
LATT,1,1,1,,,,1
lmesh,all
lsel,s,line,,4,5,1
lesize,all,1
LATT,1,2,1,,,,2
lmesh,all
lsel,s,line,,5,6,1
lesize,all,1
LATT,1,2,1,,,,2
lmesh,all
lsel,s,line,,7,8,1
lesize,all,1
LATT,1,3,1,,,,3
lmesh,all
lsel,s,line,,8,9,1
lesize,all,1
LATT,1,3,1,,,,3
lmesh,all
lsel,s,line,,10,11,1
lesize,all,1
LATT,1,4,1,,,,4
lmesh,all
lsel,s,line,,16,17,1
lesize,all,1
LATT,1,4,1,,,,4
lmesh,all
lsel,s,line,,12,13,1
lesize,all,1
LATT,2,5,1,,,,5
lmesh,all
lsel,s,line,,14,15,1
lesize,all,1
LATT,2,5,1,,,,5
lmesh,all
!
④创建弹簧单元
nplot
/pnum,node,1
/replot
type,3
real,6
nsel,s,loc,y,0
ngen,2,500,all,,,,-1
*do,i,1,25,1
e,i,i+500
*enddo
!
⑤施加边界条件
allsel,all
nsel,s,loc,y,-1
d,all,ux,0,,,,uy
allsel,all
d,2,uy,0
d,10,uy,0
d,14,ux,0
!
⑥施加荷载
/solu
q1=94.6e3
q2=12.0e3
q3=181e3
e1=67.7e3
e2=273.6e3
allsel,all
esel,s,real,,3
sfbeam,all,1,pres,q1
allsel,all
esel,s,real,,1
sfbeam,all,1,pres,-q3
allsel,all
esel,s,real,,2
sfbeam,all,1,pres,q2
allsel,all
esel,s,real,,4
nsle,s
nsel,r,loc,x,0
esln,r,1
*get,Nelm,elem,0,count
*get,Emin,elem,0,num,min
*dim,P,array,Nelm,3,1
k=Emin
*do,i,1,Nelm
P(i,1)=k
k=elnext(k)
*enddo
n=1
*do,i,n,Nelm
p(i,2)=e2-(e2-e1)*(i-n)*1/14.05
P(i,3)=e2-(e2-e1)*(i-n+1)*1/14.05
*enddo
*do,i,1,Nelm
sfbeam,P(i,1),1,pres,P(i,2),P(i,3)
*enddo
allsel,all
esel,s,real,,4
nsle,s
nsel,r,loc,x,22.5
esln,r,1
*get,Nelm,elem,0,count
*get,Emin,elem,0,num,min
*dim,P,array,Nelm,3,1
k=Emin
*do,i,1,Nelm
P(i,1)=k
k=elnext(k)
*enddo
n=1
*do,i,n,Nelm
p(i,2)=e2-(e2-e1)*(i-n)*1/14.05
P(i,3)=e2-(e2-e1)*(i-n+1)*1/14.05
*enddo
*do,i,1,Nelm
sfbeam,P(i,1),1,pres,-P(i,2),-P(i,3)
*enddo
allsel,all
acel,,9.8
finish
/SOL
/STATUS,SOLU
SOLVE
finish
/POST1
ETABLE,NI,SMISC,1
ETABLE,NJ,SMISC,7
ETABLE,QI,SMISC,2
ETABLE,QJ,SMISC,8
ETABLE,MI,SMISC,6
ETABLE,MJ,SMISC,12
PLLS,NI,NJ,1,0
PLLS,QI,QJ,1,0
PLLS,MI,MJ,-1,0
ESEL,U,TYPE,,2
!
显示变形图
pldisp,1
finish
3.纵梁命令流
/prep7
ET,1,BEAM3
ET,2,link10
KEYOPT,2,3,1
r,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,3e10
MP,PRXY,1,0.2
MP,DENS,1,2500
MP,EX,2,3.35e10
MP,PRXY,2,0.2
MP,DENS,2,2500
MP,EX,3,1.6e8
*DO,i,1,41,1
N,i,i-1,18.5
N,i+41,i-1,8.5
N,i+82,i-1,0
*ENDDO
TYPE,1
REAL,1
Mat,1
*DO,i,1,3,1
REAL,i
*DO,j,1,40,1
E,j+41*(i-1),j+1+41*(i-1)
*ENDDO
*ENDDO
REAL,4
Mat,2
*DO,i,1,41,8
*DO,j,1,80,41
E,j+(i-1),j+(i-1)+41
*ENDDO
*ENDDO
eplot
/pnum,elem,1
ESEL,S,,,1,40
SFBEAM,ALL,1,PRES,1161000
ESEL,S,,,41,80
SFBEAM,ALL,1,PRES,167000
ESEL,S,,,81,120
SFBEAM,ALL,1,PRES,-1398900
ESEL,ALL
*DO,i,1,6,1
D,83+8*(i-1),UX,0,,,,,UY
*ENDDO
ACEL,0,10,0
EPLOT
/PNUM,ELEM,1
/SOL
/STATUS,SOLU
SOLVE
finish
/POST1
ETABLE,NI,SMISC,1
ETABLE,NJ,SMISC,7
ETABLE,QI,SMISC,2
ETABLE,QJ,SMISC,8
ETABLE,MI,SMISC,6
ETABLE,MJ,SMISC,12
PLLS,NI,NJ,1,0
PLLS,QI,QJ,1,0
PLLS,MI,MJ,-1,0
ESEL,U,TYPE,,2
/PNUM,ELEM,0
/PNUM,NODE,0
/DIST,1,1.08222638492,1
/REP,FAST
/DIST,1,0.924021086472,1
/REP,FAST
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- ANSYS 地铁 车站 横断面 内力 命令 New