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

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

《ANSYS地铁车站横断面内力命令流.docx》由会员分享，可在线阅读，更多相关《ANSYS地铁车站横断面内力命令流.docx（19页珍藏版）》请在冰豆网上搜索。

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

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