火车过桥全过程仿真(64m钢桁梁)
这是本人的毕业设计,大家凑合着看看/view,1,1,1,1
/eshape,1
cj=1.27
/prep7
mp,ex,1,2e11
mp,prxy,1,0.3
mp,dens,1,780
et,1,beam188
!上下弦杆
sectype,1,beam,I
secdata,0.46,0.46,0.46,0.02,0.02,0.012
!端斜杆
sectype,2,beam,I
secdata,0.6,0.6,0.46,0.02,0.02,0.012
!斜杆
sectype,3,beam,I
secdata,0.44,0.44,0.46,0.012,0.012,0.01
!竖杆
sectype,4,beam,I
secdata,0.26,0.26,0.46,0.012,0.012,0.01
!上下纵联
sectype,5,beam,I
secdata,0.12,0.12,0.24,0.012,0.012,0.01
!纵梁
sectype,6,beam,I
secdata,0.24,0.24,1.29,0.016,0.016,0.01
!横梁
sectype,7,beam,I
secdata,0.24,0.24,1.29,0.024,0.024,0.012
n,
n,,6
n,,0,0,8
n,,6,0,8
n,,6,8,8
n,,0,8,8
n,,0,0,16
n,,6,0,16
n,,6,8,16
n,,0,8,16
n,,0,0,24
n,,6,0,24
n,,6,8,24
n,,0,8,24
n,,0,0,32
n,,6,0,32
n,,6,8,32
n,,0,8,32
n,,0,0,40
n,,6,0,40
n,,6,8,40
n,,0,8,40
n,,0,0,48
n,,6,0,48
n,,6,8,48
n,,0,8,48
n,,0,0,56
n,,6,0,56
n,,6,8,56
n,,0,8,56
n,,0,0,64
n,,6,0,64
n,500,0,10,0
n,508,0,10,64
fill
n,600,2,10,0
n,601,4,10,0
n,602,6,10,0
!上下弦杆
secnum, 1
e,1,3,500
e,3,7,500
e,7,11,500
e,11,15,500
e,15,19,500
e,19,23,500
e,23,27,500
e,27,31,500
e,2,4,602
e,4,8,602
e,8,12,602
e,12,16,602
e,16,20,602
e,20,24,602
e,24,28,602
e,28,32,602
e,5,9,602
e,9,13,602
e,13,17,602
e,17,21,602
e,21,25,602
e,25,29,602
e,6,10,500
e,10,14,500
e,14,18,500
e,18,22,500
e,22,26,500
e,26,30,500
!端斜杆
secnum,2
e,1,6
e,2,5
e,29,32
e,30,31
!斜杆
secnum,3
e,6,7
e,7,14
e,14,15
e,15,22
e,22,23
e,23,30
e,5,8
e,8,13
e,13,16
e,16,21
e,21,24
e,24,29
!竖杆
secnum,4
e,3,6
e,4,5
e,7,10
e,8,9
e,11,14
e,12,13
e,15,18
e,16,17
e,19,22
e,20,21
e,23,26
e,24,25
e,27,30
e,28,29
!上纵联
secnum,5
e,5,6,501
e,5,10,502
e,6,9,501
e,9,10,502
e,9,14,503
e,10,13,502
e,13,14,503
e,13,18,504
e,14,17,503
e,17,18,504
e,17,22,505
e,18,21,504
e,21,22,505
e,21,26,506
e,22,25,505
e,25,26,506
e,25,30,506
e,26,29,506
e,29,30,506
!下纵联
e,1,4,500
e,2,3,501
e,3,8,501
e,4,7,502
e,7,12,502
e,8,11,503
e,11,16,503
e,12,15,504
e,15,20,504
e,16,19,505
e,19,24,505
e,20,23,506
e,23,28,506
e,24,27,507
e,27,32,507
e,28,31,508
!纵梁
secnum,6
n,33,2,0,0
n,161,2,0,64
fill
e,33,34,600
*repeat,128,1,1
n,162,4,0,0
n,290,4,0,64
fill
e,162,163,601
*repeat,128,1,1
!横梁
secnum,7
e,1,33,500
e,33,162,500
e,162,2,500
e,3,49,501
e,49,178,501
e,178,4,501
e,7,65,502
e,65,194,502
e,194,8,502
e,11,81,503
e,81,210,503
e,210,12,503
e,15,97,504
e,97,226,504
e,226,16,504
e,19,113,505
e,113,242,505
e,242,20,505
e,23,129,506
e,129,258,506
e,258,24,506
e,27,145,507
e,145,274,507
e,274,28,507
e,31,161,508
e,161,290,508
e,290,32,508
d,1,all
d,2,all
ddele,1,rotx
ddele,2,rotx
d,31,all
d,32,all
ddele,31,rotx
ddele,32,rotx
ddele,31,uz
ddele,32,uz
/solu
antype,static
i=0
*do,i,0,128,1
fdele,all,all
sfedele,all,all,pres
nsel,all
nsel,s,loc,x,2,4,2
nsel,r,loc,y,0
nsel,r,loc,z,i*0.5
f,all,fy,-cj*110000
nsel,all
nsel,s,loc,x,2,4,2
nsel,r,loc,y,0
nsel,r,loc,z,(i-3)*0.5
f,all,fy,-cj*110000
nsel,all
nsel,s,loc,x,2,4,2
nsel,r,loc,y,0
nsel,r,loc,z,(i-6)*0.5
f,all,fy,-cj*110000
nsel,all
nsel,s,loc,x,2,4,2
nsel,r,loc,y,0
nsel,r,loc,z,(i-9)*0.5
f,all,fy,-cj*110000
nsel,all
nsel,s,loc,x,2,4,2
nsel,r,loc,y,0
nsel,r,loc,z,(i-12)*0.5
f,all,fy,-cj*110000
nsel,all
!
*if,i,le,15,then
*elseif,i,le,75
nsel,s,loc,x,2,4,2
nsel,r,loc,y,0
nsel,r,loc,z,0,(i-15)*0.5,1
esln
esel,u,sec,,7
sfbeam,all,1,pres,cj*46000
*elseif,i,gt,75
nsel,s,loc,x,2,4,2
nsel,r,loc,y,0
nsel,r,loc,z,(i-75)*0.5,(i-15)*0.5,1
esln
esel,u,sec,,7
sfbeam,all,1,pres,cj*46000
nsel,all
nsel,s,loc,x,2,4,2
nsel,r,loc,y,0
nsel,r,loc,z,0,(i-75)*0.5,1
esln
esel,u,sec,,7
sfbeam,all,1,pres,cj*40000
*endif
allsel
outres,all,all
solve
*enddo
finish
/post26
numvar,50
timerange,1,130
!下弦
esol,2,1,1,smisc,1,fx1
esol,3,2,3,smisc,1,fx2
esol,4,3,7,smisc,1,fx3
esol,5,4,11,smisc,1,fx4
esol,6,5,15,smisc,1,fx6
esol,7,6,19,smisc,1,fx7
esol,8,7,23,smisc,1,fx8
esol,9,8,27,smisc,1,fx9
!上弦
esol,10,23,6,smisc,1,fx23
esol,11,24,10,smisc,1,fx24
esol,12,25,14,smisc,1,fx25
esol,13,26,18,smisc,1,fx26
esol,14,27,22,smisc,1,fx27
esol,15,28,26,smisc,1,fx28
!端斜杆
esol,16,29,6,smisc,1,fx29
esol,17,32,30,smisc,1,fx32
!斜杆
esol,18,33,6,smisc,1,fx33
esol,19,34,7,smisc,1,fx34
esol,20,35,14,smisc,1,fx35
esol,21,36,15,smisc,1,fx36
esol,22,37,22,smisc,1,fx37
esol,23,38,30,smisc,1,fx38
!竖杆
esol,24,45,6,smisc,1,fx45
esol,25,47,7,smisc,1,fx47
esol,26,49,14,smisc,1,fx49
esol,27,51,15,smisc,1,fx51
esol,28,53,22,smisc,1,fx53
esol,29,55,26,smisc,1,fx55
esol,30,57,30,smisc,1,fx57
prval,2,3,4,5,6,7,8,9
prval,10,11,12,13,14,15
prval,16,17
prval,18,19,20,21,22,23
prval,24,25,26,27,28,29,30
/view,1,1,1,1
/eshape,1
cj=1.27
/prep7
mp,ex,1,2e11
mp,prxy,1,0.3
mp,dens,1,780
et,1,beam188
!上下弦杆
sectype,1,beam,I
secdata,0.46,0.46,0.46,0.02,0.02,0.012
!端斜杆
sectype,2,beam,I
secdata,0.6,0.6,0.46,0.02,0.02,0.012
!斜杆
sectype,3,beam,I
secdata,0.44,0.44,0.46,0.012,0.012,0.01
!竖杆
sectype,4,beam,I
secdata,0.26,0.26,0.46,0.012,0.012,0.01
!上下纵联
sectype,5,beam,I
secdata,0.12,0.12,0.24,0.012,0.012,0.01
!纵梁
sectype,6,beam,I
secdata,0.24,0.24,1.29,0.016,0.016,0.01
!横梁
sectype,7,beam,I
secdata,0.24,0.24,1.29,0.024,0.024,0.012
n,
n,,6
n,,0,0,8
n,,6,0,8
n,,6,8,8
n,,0,8,8
n,,0,0,16
n,,6,0,16
n,,6,8,16
n,,0,8,16
n,,0,0,24
n,,6,0,24
n,,6,8,24
n,,0,8,24
n,,0,0,32
n,,6,0,32
n,,6,8,32
n,,0,8,32
n,,0,0,40
n,,6,0,40
n,,6,8,40
n,,0,8,40
n,,0,0,48
n,,6,0,48
n,,6,8,48
n,,0,8,48
n,,0,0,56
n,,6,0,56
n,,6,8,56
n,,0,8,56
n,,0,0,64
n,,6,0,64
n,500,0,10,0
n,508,0,10,64
fill
n,600,2,10,0
n,601,4,10,0
n,602,6,10,0
!上下弦杆
secnum, 1
e,1,3,500
e,3,7,500
e,7,11,500
e,11,15,500
e,15,19,500
e,19,23,500
e,23,27,500
e,27,31,500
e,2,4,602
e,4,8,602
e,8,12,602
e,12,16,602
e,16,20,602
e,20,24,602
e,24,28,602
e,28,32,602
e,5,9,602
e,9,13,602
e,13,17,602
e,17,21,602
e,21,25,602
e,25,29,602
e,6,10,500
e,10,14,500
e,14,18,500
e,18,22,500
e,22,26,500
e,26,30,500
!端斜杆
secnum,2
e,1,6
e,2,5
e,29,32
e,30,31
!斜杆
secnum,3
e,6,7
e,7,14
e,14,15
e,15,22
e,22,23
e,23,30
e,5,8
e,8,13
e,13,16
e,16,21
e,21,24
e,24,29
!竖杆
secnum,4
e,3,6
e,4,5
e,7,10
e,8,9
e,11,14
e,12,13
e,15,18
e,16,17
e,19,22
e,20,21
e,23,26
e,24,25
e,27,30
e,28,29
!上纵联
secnum,5
e,5,6,501
e,5,10,502
e,6,9,501
e,9,10,502
e,9,14,503
e,10,13,502
e,13,14,503
e,13,18,504
e,14,17,503
e,17,18,504
e,17,22,505
e,18,21,504
e,21,22,505
e,21,26,506
e,22,25,505
e,25,26,506
e,25,30,506
e,26,29,506
e,29,30,506
!下纵联
e,1,4,500
e,2,3,501
e,3,8,501
e,4,7,502
e,7,12,502
e,8,11,503
e,11,16,503
e,12,15,504
e,15,20,504
e,16,19,505
e,19,24,505
e,20,23,506
e,23,28,506
e,24,27,507
e,27,32,507
e,28,31,508
!纵梁
secnum,6
n,33,2,0,0
n,161,2,0,64
fill
e,33,34,600
*repeat,128,1,1
n,162,4,0,0
n,290,4,0,64
fill
e,162,163,601
*repeat,128,1,1
!横梁
secnum,7
e,1,33,500
e,33,162,500
e,162,2,500
e,3,49,501
e,49,178,501
e,178,4,501
e,7,65,502
e,65,194,502
e,194,8,502
e,11,81,503
e,81,210,503
e,210,12,503
e,15,97,504
e,97,226,504
e,226,16,504
e,19,113,505
e,113,242,505
e,242,20,505
e,23,129,506
e,129,258,506
e,258,24,506
e,27,145,507
e,145,274,507
e,274,28,507
e,31,161,508
e,161,290,508
e,290,32,508
d,1,all
d,2,all
ddele,1,rotx
ddele,2,rotx
d,31,all
d,32,all
ddele,31,rotx
ddele,32,rotx
ddele,31,uz
ddele,32,uz
/solu
antype,static
i=0
*do,i,0,128,1
fdele,all,all
sfedele,all,all,pres
nsel,all
nsel,s,loc,x,2,4,2
nsel,r,loc,y,0
nsel,r,loc,z,i*0.5
f,all,fy,-cj*110000
nsel,all
nsel,s,loc,x,2,4,2
nsel,r,loc,y,0
nsel,r,loc,z,(i-3)*0.5
f,all,fy,-cj*110000
nsel,all
nsel,s,loc,x,2,4,2
nsel,r,loc,y,0
nsel,r,loc,z,(i-6)*0.5
f,all,fy,-cj*110000
nsel,all
nsel,s,loc,x,2,4,2
nsel,r,loc,y,0
nsel,r,loc,z,(i-9)*0.5
f,all,fy,-cj*110000
nsel,all
nsel,s,loc,x,2,4,2
nsel,r,loc,y,0
nsel,r,loc,z,(i-12)*0.5
f,all,fy,-cj*110000
nsel,all
!
*if,i,le,15,then
*elseif,i,le,75
nsel,s,loc,x,2,4,2
nsel,r,loc,y,0
nsel,r,loc,z,0,(i-15)*0.5,1
esln
esel,u,sec,,7
sfbeam,all,1,pres,cj*46000
*elseif,i,gt,75
nsel,s,loc,x,2,4,2
nsel,r,loc,y,0
nsel,r,loc,z,(i-75)*0.5,(i-15)*0.5,1
esln
esel,u,sec,,7
sfbeam,all,1,pres,cj*46000
nsel,all
nsel,s,loc,x,2,4,2
nsel,r,loc,y,0
nsel,r,loc,z,0,(i-75)*0.5,1
esln
esel,u,sec,,7
sfbeam,all,1,pres,cj*40000
*endif
allsel
outres,all,all
solve
*enddo
finish
/post26
numvar,50
timerange,1,130
!下弦
esol,2,1,1,smisc,1,fx1
esol,3,2,3,smisc,1,fx2
esol,4,3,7,smisc,1,fx3
esol,5,4,11,smisc,1,fx4
esol,6,5,15,smisc,1,fx6
esol,7,6,19,smisc,1,fx7
esol,8,7,23,smisc,1,fx8
esol,9,8,27,smisc,1,fx9
!上弦
esol,10,23,6,smisc,1,fx23
esol,11,24,10,smisc,1,fx24
esol,12,25,14,smisc,1,fx25
esol,13,26,18,smisc,1,fx26
esol,14,27,22,smisc,1,fx27
esol,15,28,26,smisc,1,fx28
!端斜杆
esol,16,29,6,smisc,1,fx29
esol,17,32,30,smisc,1,fx32
!斜杆
esol,18,33,6,smisc,1,fx33
esol,19,34,7,smisc,1,fx34
esol,20,35,14,smisc,1,fx35
esol,21,36,15,smisc,1,fx36
esol,22,37,22,smisc,1,fx37
esol,23,38,30,smisc,1,fx38
!竖杆
esol,24,45,6,smisc,1,fx45
esol,25,47,7,smisc,1,fx47
esol,26,49,14,smisc,1,fx49
esol,27,51,15,smisc,1,fx51
esol,28,53,22,smisc,1,fx53
esol,29,55,26,smisc,1,fx55
esol,30,57,30,smisc,1,fx57
prval,2,3,4,5,6,7,8,9
prval,10,11,12,13,14,15
prval,16,17
prval,18,19,20,21,22,23
prval,24,25,26,27,28,29,30
火车头粘在桥面下了:)
很好的算例,运行很流畅。叫板板多给点分哦! 这是下承式钢桁梁,用的是“中活载”,计算结果比影响线法手算的稍微小点,感觉还行,让大家见笑了 mp,dens,1,7800
不知是不是大意了
牛呀
很好呀 这么厉害啊 牛人真多啊:victory: 谢谢分享,好人啊,呵呵 真是太强了 牛人真多啊 太强了 粘下来,看看运行结果... 好是好就是看不懂 粘下来,看看运行结果... 好东西,看看结果先! 真是一代比一代强啊:P 1# haojingwei1651 学习学习,。 学习了,正好需要 学习了,xiexie 很贴近生活,楼主强悍 向你致敬,要好好学习下