马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
- finish
- /clear
- /prep7
- /title,the analyse of steel arch bridge
- 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
- asel,all
- 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
- /replot
- 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
- /replot
- k,1,-0.5,0
- 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
- /replot
- k,1,-0.5,0
- 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
- /replot
- k,1,-0.5,0
- 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,0
- 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
- /replot
- ldele,all
- kdele,all
- k,1,-0.5,0
- 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,0
- k,13,-0.2,0.38
- k,14,-0.3,0.38+0.2
- k,15,-0.3,1.23+0.5-0.3
- 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
-
- /replot
- k,1,-0.5,0
- 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,0
- 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
-
- /replot
- etdele,1
- /replot
- 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
- mp,ex,2,3.5e10
- mp,dens,2,2600
- mp,prxy,2,0.1667
- et,3,beam44
- mp,ex,3,3.0e10
- mp,dens,3,2600
- mp,prxy,3,0.1667
- et,4,beam44
- 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
- mp,ex,7,3.0e10
- mp,dens,7,2600
- mp,prxy,7,0.1667
- et,8,beam44
- mp,ex,8,3.3e10
- mp,dens,8,2600
- mp,prxy,8,0.1667
- et,9,beam44
- 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
- mp,ex,11,2.85e10
- mp,dens,11,2500
- mp,prxy,11,0.1667
- n,90000,0,0,30
- n,1,0,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
- 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
- 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
- egen,2,100,1,10,1,,,,,,0,4.7+1.25/2+1.5/2
- egen,17,100,11,20,1,,,,,,0,5.3+1.5
- egen,2,100,171,180,1,,,,,,0,4.7+1.5/2+1.25/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
- e,1,101
- e,2,102
- e,3,103
- e,5,105
- e,12,112
- e,13,113
- e,15,115
- *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
- 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
- *dim,z,array,16
- 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(11)=21.7
- y(12)=28.5
- y(13)=42.1
- y(14)=48.9
- y(15)=59.1
- y(16)=2.025
- f1=25
- m1=1.1
- k=log(m1+sqrt(m1*m1-1))
- *do,i,1,16,1
- cs(i)=1-y(i)*2/125
- z(i)=25-f1/(m1-1)*(cosh(k*cs(i))-1)
- *enddo
- *dim,ygd,array,23
- ygd(1) =0.0040004 +0.0129
- ygd(2) =0.010685 +0.0222
- ygd(3) =0.019483 +0.0303
- ygd(4) =0.029219 +0.0369
- ygd(5) =0.039022 +0.0423
- ygd(6) =0.048243 +0.0465
- ygd(7) =0.05639 +0.0495
- ygd(8) =0.06307 +0.0515
- ygd(9) =0.068014 +0.0521
- ygd(10) =0.071003 +0.0515
- ygd(11) =0.071739 +0.0495
- ygd(12) =0.069962 +0.0465
- ygd(13) =0.065514 +0.423
- ygd(14) =0.058289 +0.0369
- ygd(15) =0.048248 +0.0304
- ygd(16) =0.035507 +0.0224
- ygd(17) =0.20415 +0.0123
- ygd(18) =0.06577 +0.040234
- ygd(19) =0.069769 +0.041215
- 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)
- dgarea=109*3.1415926*0.25*0.007*0.007
- dgyl=0.3*1670*1e6
- dgl=dgyl*dgarea
- edg=1.95e11
- *dim,L0,array,17
- *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)
- L0(18-i)=z(i)
- *enddo
- fdg(1 )=dgl
- fdg(2 )=dgl
- fdg(3 )=dgl
- fdg(4 )=dgl
- fdg(5 )=dgl
- fdg(6 )=dgl
- fdg(7 )=dgl
- fdg(8 )=dgl
- fdg(9 )=dgl
- fdg(10)=dgl
- fdg(11)=dgl
- fdg(12)=dgl
- fdg(13)=dgl
- fdg(14)=dgl
- fdg(15)=dgl
- fdg(16)=dgl
- fdg(17)=dgl
- 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,1+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, -12.5 ,42.1 ,z(13)
- n,20024, -12.5 ,48.9 ,z(14)
- n,20025, -12.5 ,59.1 ,z(15)+ygd(18)
- n,20026, -12.5 ,65.9 ,z(15)+ygd(19)
- n,20027, -12.5 ,76.1 ,z(14)
- n,20028, -12.5 ,82.9 ,z(13)
- n,20029, -12.5 ,96.5 ,z(12)
- n,20030, -12.5 ,103.3 ,z(11)
- n,20000, 0 ,62.5 ,z(15)
- type,5
- mat,5
- sectype,500,beam,ctube,windc
- secdata,0.486,0.50
- secnum,500
- e,20002,10404,90000
- e,20002,10414,90000
- e,20002,10304,90000
- e,20002,10314,90000
- e,20004,10704,90000
- e,20004,10714,90000
- e,20004,10604,90000
- e,20004,10614,90000
- e,20000,20015,90000
- e,20000,20016,90000
- e,20000,20025,90000
- e,20000,20026,90000
- e,20007,11104,90000
- e,20007,11114,90000
- e,20007,11204,90000
- e,20007,11214,90000
- e,20009,11404,90000
- e,20009,11414,90000
- e,20009,11504,90000
- e,20002,11514,90000
- cm,fc,elem
- esel,none
- n,30000,12.5,0,0
- n,30001,12.5,125,0
- n,30002,12.5,2.025,z(16)
- n,30003,12.5,125-2.025,z(16)
- n,40000,12.5,0,0
- n,40001,12.5,125,0
- n,40002,12.5,2.025,z(16)
- n,40003,12.5,125-2.025,z(16)
- *dim,zgzx,array,23
- *dim,ygzx,array,23
- zgzx(1 )=30000
- zgzx(2 )=30002
- zgzx(3 )=10104
- zgzx(4 )=10204
- zgzx(5 )=10304
- zgzx(6 )=10404
- zgzx(7 )=10504
- zgzx(8 )=10604
- zgzx(9 )=10704
- zgzx(10 )=10804
- zgzx(11 )=20015
- zgzx(12 )=10904
- zgzx(13 )=20016
- zgzx(14 )=11004
- zgzx(15 )=11104
- zgzx(16 )=11204
- zgzx(17 )=11304
- zgzx(18 )=11404
- zgzx(19 )=11504
- zgzx(20 )=11604
- zgzx(21 )=11704
- zgzx(22 )=30003
- zgzx(23 )=30001
- zgzx(1 )=40000
- zgzx(2 )=40002
- zgzx(3 )=10114
- zgzx(4 )=10214
- zgzx(5 )=10314
- zgzx(6 )=10414
- zgzx(7 )=10514
- zgzx(8 )=10614
- zgzx(9 )=10714
- zgzx(10 )=10814
- zgzx(11 )=20025
- zgzx(12 )=10914
- zgzx(13 )=20026
- zgzx(14 )=11014
- zgzx(15 )=11114
- zgzx(16 )=11214
- zgzx(17 )=11314
- zgzx(18 )=11414
- zgzx(19 )=11514
- zgzx(20 )=11614
- zgzx(21 )=11714
- zgzx(22 )=40003
- zgzx(23 )=40001
- n,90001,12.5,0,30
- n,90002,-12.5,0,30
- type,1
- mat,1
- r,2000
- secnum,1
- *do,i,1,22
- e,zgzx(i),zgzx(i+1),90001
- e,ygzx(i),ygzx(i+1),90002
- *enddo
- cm,gzgg,elem
- esel,none
- type,2
- mat,2
- r,2001
- secnum,2
- *do,i,1,22
- e,zgzx(i),zgzx(i+1),90001
- e,ygzx(i),ygzx(i+1),90002
- *enddo
- cm,gzhnt,elem
- esel,none
- type,9
- r,6
- real,6
- mat,9
- sectype,70,beam,rect,dizuo
- secdata,1.4,4
- secnum,70
- e,30001,30003,90001
- e,40001,40003,90002
- e,30000,30002,90001
- e,30000,40002,90002
- cm,gzdz,elem
- esel,none
- ngen,2,2000,1,16,1,0,-2.025,0
- ngen,2,100,1801,1816,1,0,2.025,0
- type,7
- mat,7
- r,7
- real,7
- sectype,60,beam,rect,gail
- secdata,4,1.6
- !!!!!!!!!!*****************************以上已经校核过!
复制代码 |