钢管混凝土拱桥ansys命令流很实用很经典.docx
《钢管混凝土拱桥ansys命令流很实用很经典.docx》由会员分享,可在线阅读,更多相关《钢管混凝土拱桥ansys命令流很实用很经典.docx(40页珍藏版)》请在冰豆网上搜索。
钢管混凝土拱桥ansys命令流很实用很经典
fini
/clear
/title,concrete-filledtubearchbridge
/prep7
截面1
et,1,82
cyl4,0.6,0.6,0.6,,0.586
cyl4,0.6,2.4,0.6,,0.586
rectng,0.35,0.364,0.8,2.2
rectng,0.85,0.836,0.8,2.2
allsel
aptn,all
adele,5,8,1
allsel
aadd,all
smrtsize,5
amesh,all
secwrite,gg,sect,,1
sectype,1,beam,mesh
secoffset,cent,,,
secread,'gg','sect',,mesh
asel,all
aclear,all
adele,all,,,1
截面2
cyl4,0.6,0.6,0.586
cyl4,0.6,2.4,0.586
rectng,0.364,0.836,0.8,2.2
allsel
aadd,all
smrtsize,5
amesh,all
secwrite,hnt,sect,,1
sectype,2,beam,mesh
secoffset,cent,,,
secread,'hnt','sect',,mesh
asel,all
aclear,all
adele,all,,,1
建立横梁标准截面a-a
k,1,-0.5
k,2,-0.75,0.25
k,3,-0.75,0.5
k,4,-0.55,0.5
k,5,-0.55,0.85
k,6,0.55,0.85
k,7,0.55,0.5
k,8,0.75,0.5
k,9,0.75,0.25
k,10,0.5,0
a,1,2,3,4,5,6,7,8,9,10
smrtsize,5
amesh,all
secwrite,hl1,sect,,1
sectype,4,beam,mesh
secoffset,cent,,,
secread,'hl1','sect',,mesh
asel,all
aclear,all
adele,all,,,1
建立横梁标准截面b-b
k,1,-0.5
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
smrtsize,5
amesh,all
secwrite,hl2,sect,,1
sectype,5,beam,mesh
secoffset,cent,,,
secread,'hl2','sect',,mesh
asel,all
aclear,all
adele,all,,,1
建立横梁标准截面c-c
k,1,-0.5
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
smrtsize,5
amesh,all
secwrite,hl3,sect,,1
sectype,6,beam,mesh
secoffset,cent,,,
secread,'hl3','sect',,mesh
asel,all
aclear,all
adele,all,,,1
建立横梁标准截面d-d
ldele,all
kdele,all
k,1,-0.5
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,12,0.5
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,1,2,3,4,5,6,7,8,9,10,11,12
a,13,14,15,16,17,18,19,20
asba,1,2
smrtsize,5
amesh,all
secwrite,hl4,sect,,1
sectype,7,beam,mesh
secoffset,cent,,,
secread,'hl4','sect',,mesh
asel,all
aclear,all
adele,all,,,1
建立横梁标准截面e-e
k,1,-0.5
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,12,0.5
a,1,2,3,4,5,6,7,8,9,10,11,12
smrtsize,5
amesh,all
secwrite,hl5,sect,,1
sectype,8,beam,mesh
secoffset,cent,,,
secread,'hl5','sect',,mesh
asel,all
aclear,all
adele,all,,,1
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
secnum,7
e,1,12
secnum,6
e,12,13
secnum,5
e,13,14
e,14,15
secnum,4
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
esel,none
****************************
!
定义拱轴线
*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)!
拱轴高度
*enddo
!
上述定义的拱轴线可参照姚玲森《桥梁工程》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)!
吊杆节点
!
*enddo
!
!
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
!
*enddo
!
!
为了对以后吊杆转动进行放松,重复吊杆节点
!
*do,i,1,9,1
!
ngen,2,3000,i*100+4,i*100+14,10,0,0,z(i)+ygd(i)!
吊杆下节点
!
*enddo
!
!
jj=1
!
*do,i,10,17,1
!
ngen,2,3000,i*100+4,i*100+14,10,0,0,z(i-2*jj)+ygd(i)!
吊杆节点
!
jj=jj+1
!
*enddo
!
!
*do,i,1,10,1
!
ngen,2,6000,i*100+4,i*100+14,10,0,0,z(i)+ygd(i)!
吊杆节点
!
*enddo
!
!
jj=1
!
*do,i,10,17,1
!
ngen,2,6000,i*100+4,i*100+14,10,0,0,z(i-2*jj)+ygd(i)!
吊杆节点
!
jj=jj+1
!
*enddo
!
定义吊杆的参数(每根吊杆均采用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!
定义在吊杆上施加的力数组
*do,i,1,9,1
L0(i)=z(i)!
0应力时单元长度
L0(18-i)=z(i)!
吊杆沿中轴对称
*enddo
*do,i,1,17,1
fdg(i)=dgl!
各吊杆施加的预应力
*enddo
type,6
mat,6
*do,i,1,17,1
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)!
定义吊杆的实常数特性
*enddo
*do,i,1,17,1
real,i+1000!
赋予吊杆的实常数特性
e,i*100+4+3000,i*100+4+6000!
生成吊杆单元
e,i*100+10+4+3000,i*100+10+4+6000
*enddo
cm,dg,elem
esel,none
allsel
!
!
*********************************
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,-