钢管混凝土拱桥ansys命令流很实用很经典Word文档下载推荐.docx

钢管混凝土拱桥ansys命令流很实用很经典Word文档下载推荐.docx

《钢管混凝土拱桥ansys命令流很实用很经典Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《钢管混凝土拱桥ansys命令流很实用很经典Word文档下载推荐.docx（40页珍藏版）》请在冰豆网上搜索。

k,9,0.75,0.25

k,10,0.5,0

a,1,2,3,4,5,6,7,8,9,10

secwrite,hl1,sect,,1

sectype,4,beam,mesh

hl1'

建立横梁标准截面b-b

k,2,-0.5,1

k,3,-0.75,1.25

k,4,-0.75,1.5

k,5,0.75,1.5

k,6,0.75,1.25

k,7,0.5,1

k,8,0.5

a,1,2,3,4,5,6,7,8

secwrite,hl2,sect,,1

sectype,5,beam,mesh

hl2'

建立横梁标准截面c-c

k,2,-0.5,1.04

k,3,-0.75,1.04+0.25

k,4,-0.75,1.04+0.5

k,5,-0.55,1.04+0.5

k,6,-0.55,1.04+0.5+0.35

k,7,0.55,1.04+0.5+0.35

k,8,0.55,1.04+0.5

k,9,0.75,1.04+0.5

k,10,0.75,1.04+0.25

k,11,0.5,1.04

k,12,0.5

k,13,-0.2,0.58

k,14,-0.2,1.04+0.5-0.3

k,15,0.2,1.04+0.5-0.3

k,16,0.2,0.58

a,1,2,3,4,5,6,7,8,9,10,11,12

a,13,14,15,16

asba,1,2

secwrite,hl3,sect,,1

sectype,6,beam,mesh

hl3'

建立横梁标准截面d-d

ldele,all

kdele,all

k,2,-0.5,1.23

k,3,-0.75,1.23+0.25

k,4,-0.75,1.23+0.5

k,5,-0.55,1.23+0.5

k,6,-0.55,1.23+0.5+0.35

k,7,0.55,1.23+0.5+0.35

k,8,0.55,1.23+0.5

k,9,0.75,1.23+0.5

k,10,0.75,1.23+0.25

k,11,0.5,1.23

k,13,-0.2,0.38

k,14,-0.3,0.38+0.2

k,15,-0.3,1.23+0.5-0.1-0.2

k,16,-0.2,1.23+0.5-0.1

k,17,0.2,1.23+0.5-0.1

k,18,0.3,1.23+0.5-0.1-0.2

k,19,0.3,0.38+0.2

k,20,0.2,0.38

a,13,14,15,16,17,18,19,20

secwrite,hl4,sect,,1

sectype,7,beam,mesh

hl4'

建立横梁标准截面e-e

secwrite,hl5,sect,,1

sectype,8,beam,mesh

hl5'

etdele,1!

!

将辅助单元删除

***********************

定义所有材料特性

et,1,beam44!

钢管特性

mp,ex,1,2.1e11

mp,dens,1,7800

mp,prxy,1,0.3

n,90000,0,0,30!

参考点

et,2,beam44!

钢管内50＃混凝土特性

mp,ex,2,3.5e10

mp,dens,2,2600

mp,prxy,2,0.1667

et,3,beam44!

纵梁30＃混凝土钢管特性

mp,ex,3,3.0e10

mp,dens,3,2600

mp,prxy,3,0.1667

et,4,beam44!

横梁30＃混凝土钢管特性

mp,ex,4,3.0e10

mp,dens,4,2600

mp,prxy,4,0.1667

et,5,beam44!

风撑特性

mp,ex,5,2.1e11

mp,dens,5,7800

mp,prxy,5,0.3

et,6,link10!

吊杆特性（钢绞线）

mp,ex,6,1.9e11

mp,dens,6,7800

mp,prxy,6,0.3

keyopt,6,3,0!

只拉吊杆

et,7,beam44!

盖梁30＃混凝土特性

mp,ex,7,3.0e10

mp,dens,7,2600

mp,prxy,7,0.1667

et,8,beam44!

墩柱40＃混凝土特性

mp,ex,8,3.3e10

mp,dens,8,2600

mp,prxy,8,0.1667

et,9,beam44!

横梁30＃混凝土特性

mp,ex,9,3.0e10

mp,dens,9,2600

mp,prxy,9,0.1667

et,10,link10!

系杆特性（钢绞线）

mp,ex,10,2.0e11

mp,dens,10,7800

mp,prxy,10,0.3

keyopt,10,3,0!

只拉系杆

et,11,beam44!

承台桩基基础25＃混凝土特性

mp,ex,11,2.85e10

mp,dens,11,2500

mp,prxy,11,0.1667

##############################

开始建立模型

0号横梁（端横梁）

n,1,,2.025

n,2,6.25,2.025

n,3,6.25+5.1,2.025

n,4,6.25+5.1+0.2+1.9/2,2.025!

12.5m位置

n,5,6.25+5.1+2.3,2.025

n,6,6.25+5.1+2.3+2.23,2.025

nsym,x,10,2,6,1

type,4

r,1

real,1

mat,4

secnum,7

单元1到10号为0号横梁

e,1,2

secnum,6

e,2,3

secnum,5

e,3,4

e,4,5

secnum,4

e,5,6

e,1,12

e,12,13

e,13,14

e,14,15

e,15,16

1号端横梁

egen,2,100,1,10,1,,,,,,0,4.7+1.25/2+1.5/2

其余端横梁2～17号

egen,17,100,11,20,1,,,,,,0,5.3+1.5

18号端横梁

egen,2,100,171,180,1,,,,,,0,4.7+1.25/2+1.5/2

cm,hl,elem

esel,none

********************************

纵梁单元

type,3

r,2

real,2

mat,3

sectype,10,beam,rect,zl

secdata,0.4,1.2

secnum,10

0号纵梁

e,1,101

e,2,102

e,3,103

e,5,105

e,12,112

e,13,113

e,15,115

中间纵梁1～16号

*do,i,0,15,1

e,101+i*100,201+i*100

e,102+i*100,202+i*100

e,104+i*100,204+i*100

e,112+i*100,212+i*100

e,114+i*100,214+i*100

*enddo

17号纵梁（端横梁）

e,1701,1801

e,1702,1802

e,1703,1803

e,1705,1805

e,1712,1812

e,1713,1813

e,1715,1815

cm,zl,elem

****************************

定义拱轴线

*dim,y,array,16!

将纵桥向距离定义为数组

*dim,cs,array,16!

定义分值cs＝（1－2x/l）

*dim,z,array,16!

将拱高度定义为数组

吊杆上的点

首先定义吊杆的y向（顺桥向）坐标

y

（1）=8.1

y

（2）=14.9

y（3）=21.7

y（4）=28.5

y（5）=35.3

y（6）=42.1

y（7）=48.9

y（8）=55.7

y（9）=62.5

风撑上的点

首先定义风撑的y向（顺桥向）坐标

y（11）=21.7

y（12）=28.5

y（13）=42.1

y（14）=48.9

y（15）=59.1

0号横梁的y向坐标

y（16）=2.025

fl=25!

矢高

ml=1.1

k=log（ml+sqrt（ml*ml-1））!

k=ch-1（m）

*do,i,1,16,1

cs（i）=1-y（i）*2/125

z（i）=25-fl/（ml-1）*（cosh（k*cs（i））-1）!

拱轴高度

上述定义的拱轴线可参照姚玲森《桥梁工程》p1299

*dim,ygd,array,23!

定义预工拱度（是在后续工况中反复调算的结果）

!

初始运行中可将预拱度全部设为0

ygd

（1）=0.4000e-02+1.29e-02

ygd

（2）=0.10685e-01+2.22e-02

ygd（3）=0.19483e-01+3.03e-02

ygd（4）=0.29219e-01+3.69e-02

ygd（5）=0.39022e-01+4.23e-02

ygd（6）=0.48243e-01+4.65e-02

ygd（7）=0.56390e-01+4.95e-02

ygd（8）=0.63070e-01+5.15e-02

ygd（9）=0.68014e-01+5.21e-02

ygd（10）=0.71003e-01+5.15e-02

ygd（11）=0.71739e-01+4.95e-02

ygd（12）=0.69962e-01+4.65e-02

ygd（13）=0.65514e-01+4.23e-02

ygd（14）=0.58289e-01+3.69e-02

ygd（15）=0.48248e-01+3.04e-02

ygd（16）=0.35507e-01+2.24e-02

ygd（17）=0.20415e-01+1.23e-02

ygd（18）=0.65770e-01+0.40234e-01

ygd（19）=0.69769e-01+0.41215e-01

生成各吊杆的竖向节点（z向）

ngen,2,10000,104,114,10,0,0,z

（1）+ygd

（1）

ngen,2,10000,204,214,10,0,0,z

（2）+ygd

（2）

ngen,2,10000,304,314,10,0,0,z（3）+ygd（3）

ngen,2,10000,404,414,10,0,0,z（4）+ygd（4）

ngen,2,10000,504,514,10,0,0,z（5）+ygd（5）

ngen,2,10000,604,614,10,0,0,z（6）+ygd（6）

ngen,2,10000,704,714,10,0,0,z（7）+ygd（7）

ngen,2,10000,804,814,10,0,0,z（8）+ygd（8）

ngen,2,10000,904,914,10,0,0,z（9）+ygd（9）

ngen,2,10000,1004,1014,10,0,0,z（8）+ygd（10）

ngen,2,10000,1104,1114,10,0,0,z（7）+ygd（11）

ngen,2,10000,1204,1214,10,0,0,z（6）+ygd（12）

ngen,2,10000,1304,1314,10,0,0,z（5）+ygd（13）

ngen,2,10000,1404,1414,10,0,0,z（4）+ygd（14）

ngen,2,10000,1504,1514,10,0,0,z（3）+ygd（15）

ngen,2,10000,1604,1614,10,0,0,z

（2）+ygd（16）

ngen,2,10000,1704,1714,10,0,0,z

（1）+ygd（17）

ngen,2,3000,104,114,10,0,0,0

ngen,2,3000,204,214,10,0,0,0

ngen,2,3000,304,314,10,0,0,0

ngen,2,3000,404,414,10,0,0,0

ngen,2,3000,504,514,10,0,0,0

ngen,2,3000,604,614,10,0,0,0

ngen,2,3000,704,714,10,0,0,0

ngen,2,3000,804,814,10,0,0,0

ngen,2,3000,904,914,10,0,0,0

ngen,2,3000,1004,1014,10,0,0,0

ngen,2,3000,1104,1114,10,0,0,0

ngen,2,3000,1204,1214,10,0,0,0

ngen,2,3000,1304,1314,10,0,0,0

ngen,2,3000,1404,1414,10,0,0,0

ngen,2,3000,1504,1514,10,0,0,0

ngen,2,3000,1604,1614,10,0,0,0

ngen,2,3000,1704,1714,10,0,0,0

ngen,2,6000,104,114,10,0,0,z

（1）+ygd

（1）

ngen,2,6000,204,214,10,0,0,z

（2）+ygd

（2）

ngen,2,6000,304,314,10,0,0,z（3）+ygd（3）

ngen,2,6000,404,414,10,0,0,z（4）+ygd（4）

ngen,2,6000,504,514,10,0,0,z（5）+ygd（5）

ngen,2,6000,604,614,10,0,0,z（6）+ygd（6）

ngen,2,6000,704,714,10,0,0,z（7）+ygd（7）

ngen,2,6000,804,814,10,0,0,z（8）+ygd（8）

ngen,2,6000,904,914,10,0,0,z（9）+ygd（9）

ngen,2,6000,1004,1014,10,0,0,z（8）+ygd（10）

ngen,2,6000,1104,1114,10,0,0,z（7）+ygd（11）

ngen,2,6000,1204,1214,10,0,0,z（6）+ygd（12）

ngen,2,6000,1304,1314,10,0,0,z（5）+ygd（13）

ngen,2,6000,1404,1414,10,0,0,z（4）+ygd（14）

ngen,2,6000,1504,1514,10,0,0,z（3）+ygd（15）

ngen,2,6000,1604,1614,10,0,0,z

（2）+ygd（16）

ngen,2,6000,1704,1714,10,0,0,z

（1）+ygd（17）

*do,i,1,9,1

ngen,2,10000,i*100+4,i*100+14,10,0,0,z（i）+ygd（i）!

吊杆节点

jj=1

*do,i,10,17,1

ngen,2,10000,i*100+4,i*100+14,10,0,0,z（i-2*jj）+ygd（i）!

jj=jj+1

为了对以后吊杆转动进行放松，重复吊杆节点

ngen,2,3000,i*100+4,i*100+14,10,0,0,z（i）+ygd（i）!

吊杆下节点

ngen,2,3000,i*100+4,i*100+14,10,0,0,z（i-2*jj）+ygd（i）!

*do,i,1,10,1

ngen,2,6000,i*100+4,i*100+14,10,0,0,z（i）+ygd（i）!

ngen,2,6000,i*100+4,i*100+14,10,0,0,z（i-2*jj）+ygd（i）!

定义吊杆的参数（每根吊杆均采用109根直径7mm高强钢丝）

dgarea=109*3.141593*0.25*0.007*0.007

dgyl=0.3*1670*1e6!

吊杆除应力采用0.3倍抗拉强度（n/mm*mm）

dgl=dgyl*dgarea

edg=1.95e11!

吊杆弹模

*dim,L0,array,17!

定义吊杆的初始0应力长度数组

*dim,L,array,17!

定义吊杆的受力后长度数组

*dim,deltL,array,17!

定义吊杆的受力后长度增量数组

*dim,istrn,array,17!

定义吊杆的应力数组

*dim,fdg,array,17!

定义在吊杆上施加的力数组

L0（i）=z（i）!

0应力时单元长度

L0（18-i）=z（i）!

吊杆沿中轴对称

*do,i,1,17,1

fdg（i）=dgl!

各吊杆施加的预应力

type,6

mat,6

deltL（i）=fdg（i）*L0（i）/（edg*dgarea）!

吊杆受力后的增长量

L（i）=L0（i）+deltL（i）

istrn（i）=deltL（i）/L（i）!

吊杆受力后的应变

r,i+1000,dgarea,istrn（i）!

定义吊杆的实常数特性

real,i+1000!

赋予吊杆的实常数特性

e,i*100+4+3000,i*100+4+6000!

生成吊杆单元

e,i*100+10+4+3000,i*100+10+4+6000

cm,dg,elem

*********************************

real,1018!

定义一个全为零的实常数便于利用截面特性时，不产生冲突

生成风撑上的节点

n,20001,0,21.7,z（11）

n,20002,0,28.5,z（12）

n,20003,0,42.1,z（13）

n,20004,0,48.9,z（14）

n,20005,0,59.1,z（15）

n,20006,0,65.9,z（15）

n,20007,0,76.1,z（14）

n,20008,0,82.9,z（13）

n,20009,0,96.5,z（12）

n,20010,0,103.3,z（11）

n,20011,12.5,21.7,z（11）

n,20012,12.5,28.5,z（12）

n,20013,12.5,42.1,z（13）

n,20014,12.5,48.9,z（14）

n,20015,12.5,59.1,z（15）+ygd（18）

n,20016,12.5,65.9,z（15）+ygd（19）

n,20017,12.5,76.1,z（14）

n,20018,12.5,82.9,z（13）

n,20019,12.5,96.5,z（12）

n,20020,12.5,103.3,z（11）

n,20021,-12.5,21.7,z（11）

n,20022,-12.5,28.5,z（12）

n,20023,-