haojingwei1651 发表于 2006-9-17 14:27:46

火车过桥全过程仿真(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


inkey 发表于 2006-9-19 20:27:55

火车头粘在桥面下了:)
很好的算例,运行很流畅。叫板板多给点分哦!

haojingwei1651 发表于 2006-9-20 09:43:14

这是下承式钢桁梁,用的是“中活载”,计算结果比影响线法手算的稍微小点,感觉还行,让大家见笑了

inkey 发表于 2006-9-23 16:10:13

mp,dens,1,7800
不知是不是大意了

rende1982 发表于 2006-9-24 17:15:24

牛呀

很好呀

850776307 发表于 2010-11-11 22:24:50

这么厉害啊

sihansmart 发表于 2010-11-16 13:50:17

牛人真多啊:victory:

tangwu169 发表于 2010-11-17 17:36:29

谢谢分享,好人啊,呵呵

haycgz 发表于 2010-11-29 19:05:51

真是太强了

zhangxin1210 发表于 2010-11-30 16:10:12

牛人真多啊 太强了

chajian0704 发表于 2011-4-21 21:38:53

粘下来,看看运行结果...

muweijie 发表于 2011-5-18 14:41:36

好是好就是看不懂

dechong 发表于 2011-5-19 07:31:56

粘下来,看看运行结果...

wurong 发表于 2011-5-25 12:24:58

好东西,看看结果先!

徐明海 发表于 2011-6-15 19:56:21

真是一代比一代强啊
:P

fqh888 发表于 2011-6-16 09:45:33

1# haojingwei1651 学习学习,。

webccbkky 发表于 2011-6-18 08:28:13

学习了,正好需要

wqqzgd 发表于 2011-7-25 13:14:29

学习了,xiexie

chenjinliang 发表于 2011-8-16 09:33:14

很贴近生活,楼主强悍

carville1981 发表于 2011-9-5 08:14:19

向你致敬,要好好学习下
页: [1] 2 3
查看完整版本: 火车过桥全过程仿真(64m钢桁梁)