finish
/clear,start
/filname,truss
/prep7
et,1,link8
et,2,beam44
mp,ex,1,2.1e11 !钢
mp,prxy,1,0.3
mp,dens,1,7850
r,1,0.10488E-01
!SECTYPE, 1, BEAM, I, , 0
!SECOFFSET, CENT
!SECDATA,0.302,0.302,0.294,0.012,0.012,0.012,0,0,0,0 !上弦杆
r,2,0.70080E-02
!SECTYPE, 2, BEAM, I, , 0
!SECOFFSET, CENT
!SECDATA,0.204,0.204,0.200,0.012,0.012,0.012,0,0,0,0 !下弦杆
r,3,0.14256E-01
!SECTYPE, 3, BEAM, I, , 0
!SECOFFSET, CENT
!SECDATA,0.348,0.348,0.344,0.016,0.016,0.010,0,0,0,0 !斜杆
r,4,0.10488E-01
!SECTYPE, 4, BEAM, I, , 0
!SECOFFSET, CENT
!SECDATA,0.302,0.302,0.294,0.012,0.012,0.012,0,0,0,0 !吊杆
!
SECTYPE, 5, BEAM, I, , 0
SECOFFSET, CENT
SECDATA,0.200,0.200,0.294,0.012,0.012,0.008,0,0,0,0 !纵梁
SECTYPE, 6, BEAM, I, , 0
SECOFFSET, CENT
SECDATA,0.302,0.302,0.294,0.012,0.012,0.012,0,0,0,0 !横梁
r,7,0.49975E-02
!SECTYPE, 7, BEAM, I, , 0
!SECOFFSET, CENT
!SECDATA,0.175,0.175,0.175,0.011,0.011,0.0075,0,0,0,0 !横撑
r,8,0.49975E-02
!SECTYPE, 8, BEAM, I, , 0
!SECOFFSET, CENT
!SECDATA,0.175,0.175,0.175,0.011,0.011,0.0075,0,0,0,0 !上平联
r,9,0.29455E-02
!SECTYPE, 9, BEAM, I, , 0
!SECOFFSET, CENT
!SECDATA,0.125,0.125,0.125,0.009,0.009,0.0065,0,0,0,0 !下平联
n,1,0,0,0
ngen,13,1,1,,,5,,,
n,14,5,5
n,15,10,6.3
n,16,15,7.3
n,17,20,8
n,18,25,8.4
n,19,30,8.5
n,20,35,8.4
n,21,40,8
n,22,45,7.3
n,23,50,6.3
n,24,55,5
n,1000,0,100,0
!real,1
type,1
mat,1
!secnum,2
real,2
*do,i,1,12,1
e,i,i+1,1000
*enddo
!secnum,1
real,1
e,1,14,1000
*do,i,14,23,1
e,i,i+1,1000
*enddo
e,24,13,1000
!secnum,3
real,3
*do,i,3,11,2
e,i,i+13,1000
e,i,i+11,1000
*enddo
!secnum,4
real,4
*do,i,2,12,1
e,i,i+12,1000
*enddo
/view,1,1,1,1
allsel
egen,2,45,1,45,1,,,,,,,,5.6
ngen,7,69,1,13,1,,,0.8
!关键点
nsel,s,,,1,13,1
nsel,a,,,70
nsel,a,,,139
nsel,a,,,208
nsel,a,,,277
nsel,a,,,346
nsel,a,,,415
nsel,a,,,46
ngen,2,1000,all,,,,30,
type,2
mat,1
secnum,5
real,5
*do,i,70,81,1
e,i,i+1,1070
e,i+69,(i+69)+1,1139
e,i+69*2,(i+69*2)+1,1208
e,i+69*3,(i+69*3)+1,1277
e,i+69*4,(i+69*4)+1,1346
e,i+69*5,(i+69*5)+1,1415
*enddo
secnum,6
real,6
*do,i,1,13,1
e,i,i+69,1000+i
e,i+69,i+69*2,1000+i
e,i+69*2,i+69*3,1000+i
e,i+69*3,i+69*4,1000+i
e,i+69*4,i+69*5,1000+i
e,i+69*5,i+69*6,1000+i
e,i+69*6,i+45,1000+i
*enddo
type,1
mat,1
real,7
*do,i,14,24,1
e,i,i+45
*enddo
real,8
*do,i,14,23,1
e,i,i+46
*enddo
*do,i,15,24,1
e,i,i+44
*enddo
real,9
*do,i,1,12,1
e,i,i+46
*enddo
*do,i,2,13,1
e,i,i+44
*enddo
D,1,UX
D,1,UY
D,1,UZ
!D,1,ROTX
!D,1,ROTY
D,13,UY
D,13,UZ
!D,13,ROTX
!D,13,ROTY
d,46,UX
d,46,UY
d,58,UY
save
finish
/solu
antype,static
allsel
acel,,9.8,
solve
finish
/post1
etable,fx,ls,1
PLLS,FX,FX,0.1,0
pretab,fx
etable,myi,smisc,5
etable,myj,smisc,11
plls,myi,myj,0.5
pretab,myi,myj
这座桥的模型就是这样了,但是还是不完善:我想用虚拟纵梁的办法加活载(比如人群荷载),但却不知如何实现,希望斑竹以及各位高手能够继续完善,帮我一下。 |