火车过桥全过程仿真
- /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
复制代码 |