预应力混凝土模型

chinshlia

本模型为预应力混凝土连三跨续梁桥，采用188梁单元，8杆单元，采用约束方程法将杆单元与梁单元连接，采用初应变实现预应力施加。
/prep7
!武汉理工大学桥梁工程09级
*dim,hw,array,9,1
*dim,Tw,array,9,1
*dim,H,array,9,1
*dim,nn,array,24,1
HW(1)=0.35,0.41,0.47,0.52,0.57,0.61,0.66,0.70,0.70   !腹板厚
TW(1)=0.39,0.43,0.47,0.505,0.54,0.57,0.60,0.80,1.60  !底板厚
H(1)=2.90,2.96,3.14,3.41,3.76,4.14,4.58,5.10,5.10    !梁高
nn(1)=72.6,72,69,67,63,59,55,51,47.5,44              !节点坐标
nn(11)=41,38,34.5,32,29.5,26,23,20,16.5,13
nn(21)=9,5,1,0
et,1,beam188!主梁
et,2,plane82
et,3,link8!预应力钢筋
mp,ex,1,3.55e10
mp,prxy,1,0.167
mp,dens,1,2650
*do,i,1,9
k,1,3.2,0
k,2,3.2,h(i)-0.65
k,3,4.68,h(i)-0.35
k,4,4.68,h(i)
ksymm,x,1,4,1,4
k,9,3.2-0.3-hw(i),tw(i)
k,10,3.2-hw(i),tw(i)+0.3
k,11,3.2-hw(i),h(i)-0.65
k,12,3.2-hw(i)-0.6,h(i)-0.35
ksymm,x,9,12,1,4
a,1,2,3,4,8,7,6,5
a,9,10,11,12,16,15,14,13
asba,1,2
wprota,0,90
wpoffs,0,0,-tw(i)
asbw,all
wpoffs,0,0,-0.3
asbw,all
wpoffs,0,0,0.95+tw(i)-h(i)
asbw,all
wpcsys
kwpave,6
wprota,0,0,90
asbw,all
wpoffs,0,0,hw(i)
asbw,8
wpoffs,0,0,0.6
asbw,9
kwpave,12
asbw,10
wpoffs,0,0,0.6
asbw,11
wpoffs,0,0,hw(i)
asbw,12
kwpave,13
asbw,1
kwpave,9
asbw,14
wpcsys
wpstyl
mshkey,1
mshape,0,2d
gplot
/pnum,line,1
/rep
lsel,all
lplot
lsel,s,line,,5,27,22
lesize,all,,,2
lsel,s,line,,4,29,25
lsel,a,line,,31,33,2
lsel,a,line,,3,35,32
lesize,all,,,2
lsel,s,line,,22,30,8
lsel,a,line,,21,39,18
lesize,all,,,2
lsel,s,line,,13,32,19
lsel,a,line,,11,36,25
lesize,all,,,2
lsel,s,line,,6,28,22
lsel,a,line,,2,38,36
lesize,all,,,3
lsel,s,line,,12,34,22
lesize,all,,,4
lsel,s,line,,1,23,22
lesize,all,,,2
lsel,s,line,,19,20
lesize,all,,,2
lsel,s,line,,37,43,6
lesize,all,,,2
lsel,s,line,,17,18
lesize,all,,,2
lsel,s,line,,8,40,32
lesize,all,,,2
lsel,s,line,,15,24,9
lsel,a,line,,7,9,2
lesize,all,,,1
lsel,s,line,,16,42,26
lesize,all,,,6
lsel,s,line,,14,26,12
lsel,a,line,,10,25,15
lesize,all,,,6
amesh,all
gplot
allsel,all
secwrite,%i%,sect,,2
aclear,all
adele,all
ldele,all
kdele,all
*enddo
*do,i,1,24
n,i,0,0,nn(i)
*enddo
/pnum,node,1
nsym,z,24,1,24
nummrg,all
numcmp,all
n,1000,0,50
!******直线段******************
*do,i,1,9
sectype,i,beam,mesh
secread,%i%,sect,,mesh
*enddo
type,1
mat,1
secnum,1
*do,i,1,5
e,i,i+1,1000
e,i+24,i+25,1000
*enddo
e,22,23,1000
e,23,24,1000
e,46,47,1000
e,47,24,1000
secnum,9
e,13,14,1000
e,14,15,1000!根部块
e,37,38,1000
e,38,39,1000
!********************************
!******曲线段********************
*do,i,1,7
sectype,i,beam,mesh
secread,%i%,sect,,mesh
secoffset,user,0,h(i)-h(1)
sectype,i+1,beam,mesh
secread,%i+1%,sect,,mesh
secoffset,user,0,h(i+1)-h(1)
sectype,10+i,taper
secdata,i,0,0,nn(5+i)
secdata,i+1,0,0,nn(i+1+5)
type,1
mat,1
secnum,10+i
e,i+5,i+1+5,1000
*enddo
*do,i,1,7
sectype,20+i,taper
secdata,i,0,0,nz(23-i)
secdata,i+1,0,0,nz(22-i)
type,1
mat,1
secnum,20+i
e,23-i,22-i,1000
*enddo
!*******************
!****对称跨*********
*do,i,1,7
sectype,30+i,taper
secdata,i,0,0,nz(47-i)
secdata,i+1,0,0,nz(46-i)
type,1
mat,1
secnum,30+i
e,46-i,47-i,1000
*enddo
*do,i,1,7
sectype,40+i,taper
secdata,i,0,0,nz(29+i)
secdata,i+1,0,0,nz(30+i)
type,1
mat,1
secnum,40+i
e,29+i,30+i,1000
*enddo
/eshape,1
eplot
gplot
!****建立预应力筋***********
mp,ex,2,195e9
mp,prxy,2,0.333
r,1,0.000165*9,0.00667923
!****1号到6号钢筋***********
!1号钢筋
n,51,3.04,-0.09,nz(12)
n,52,3.04,1.51,nz(13)
n,53,3.04,1.51,nz(14)
n,54,3.04,1.51,nz(15)
n,55,3.04,-0.09,nz(16)
type,3
mat,2
real,1
e,51,52
*repeat,4,1,1
!2号钢筋
n,56,3.04,0.19,nz(11)
n,57,3.04,(0.19+1.69)/2,nz(12)
n,58,3.04,1.69,nz(13)
n,59,3.04,1.69,nz(14)
n,60,3.04,1.69,nz(15)
n,61,3.04,(0.19+1.69)/2,nz(16)
n,62,3.04,0.19,nz(17)
real,1
e,56,57
*repeat,6,1,1
!3号钢筋
n,63,3.04,0.47,nz(10)
n,64,3.04,1.87,nz(11)
n,65,3.04,1.87,nz(12)
n,66,3.04,1.87,nz(13)
n,67,3.04,1.87,nz(14)
n,68,3.04,1.87,nz(15)
n,69,3.04,1.87,nz(16)
n,70,3.04,1.87,nz(17)
n,71,3.04,0.47,nz(18)
real,1
e,63,64
*repeat,8,1,1
!4号钢筋
n,72,3.04,0.75,nz(9)
n,73,3.04,2.05,nz(10)
n,74,3.04,2.05,nz(11)
n,75,3.04,2.05,nz(12)
n,76,3.04,2.05,nz(13)
n,77,3.04,2.05,nz(14)
n,78,3.04,2.05,nz(15)
n,79,3.04,2.05,nz(16)
n,80,3.04,2.05,nz(17)
n,81,3.04,2.05,nz(18)
n,82,3.04,0.75,nz(19)
real,1
e,72,73
*repeat,10,1,1
!5号钢筋
n,83,3.04,1.03,nz(8)
n,84,3.04,2.23,nz(9)
n,85,3.04,2.23,nz(10)
n,86,3.04,2.23,nz(11)
n,87,3.04,2.23,nz(12)
n,88,3.04,2.23,nz(13)
n,89,3.04,2.23,nz(14)
n,90,3.04,2.23,nz(15)
n,91,3.04,2.23,nz(16)
n,92,3.04,2.23,nz(17)
n,93,3.04,2.23,nz(18)
n,94,3.04,2.23,nz(19)
n,95,3.04,1.03,nz(20)
real,1
e,83,84
*repeat,12,1,1
!6号钢筋
n,96,3.04,1.21,nz(7)
n,97,3.04,2.41,nz(8)
n,98,3.04,2.41,nz(9)
n,99,3.04,2.41,nz(10)
n,100,3.04,2.41,nz(11)
n,101,3.04,2.41,nz(12)
n,102,3.04,2.41,nz(13)
n,103,3.04,2.41,nz(14)
n,104,3.04,2.41,nz(15)
n,105,3.04,2.41,nz(16)
n,106,3.04,2.41,nz(17)
n,107,3.04,2.41,nz(18)
n,108,3.04,2.41,nz(19)
n,109,3.04,2.41,nz(20)
n,110,3.04,1.21,nz(21)
real,1
e,96,97
*repeat,14,1,1
!7号钢筋
*do,i,111,115
n,i,2.24,2.59,nz(i-99)
*enddo
real,1
e,111,112
*repeat,4,1,1
*do,i,111,115
n,i+5,2.64,2.59,nz(i-99)
*enddo
e,116,117
*repeat,4,1,1
*do,i,111,115
n,i+10,3.04,2.59,nz(i-99)
*enddo
e,121,122
*repeat,4,1,1
*do,i,111,115
n,i+15,3.44,2.59,nz(i-99)
*enddo
e,126,127
*repeat,4,1,1
!8号钢筋
*do,i,132,138
n,i,2.44,2.59,nz(i-121)
*enddo
real,1
e,132,133
*repeat,6,1,1
*do,i,132,138
n,i+7,2.84,2.59,nz(i-121)
*enddo
real,1
e,139,140
*repeat,6,1,1
*do,i,132,138
n,i+14,3.24,2.59,nz(i-121)
*enddo
real,1
e,146,147
*repeat,6,1,1
*do,i,132,138
n,i+21,3.64,2.59,nz(i-121)
*enddo
real,1
e,153,154
*repeat,6,1,1
!9号钢筋
*do,i,160,168
n,i,1.04,2.77,nz(i-150)
*enddo
real,1
e,160,161
*repeat,8,1,1
*do,i,160,168
n,i+9,1.84,2.77,nz(i-150)
*enddo
e,169,170
*repeat,8,1,1
!10号钢筋
*do,i,178,188
n,i,0.64,2.77,nz(i-169)
*enddo
real,1
e,178,179
*repeat,10,1,1
*do,i,178,188
n,i+11,2.84,2.77,nz(i-169)
*enddo
e,189,190
*repeat,10,1,1
!11号钢筋
*do,i,200,212
n,i,0.44,2.77,nz(i-192)
*enddo
real,1
e,200,201
*repeat,12,1,1
*do,i,200,212
n,i+13,1.44,2.77,nz(i-192)
*enddo
e,213,214
*repeat,12,1,1
*do,i,200,212
n,i+26,2.44,2.77,nz(i-192)
*enddo
e,226,227
*repeat,12,1,1
!12号钢筋
*do,i,239,253
n,i,1.64,2.77,nz(i-232)
*enddo
real,1
e,239,240
*repeat,14,1,1
*do,i,239,253
n,i+15,3.04,2.77,nz(i-232)
*enddo
e,254,255
*repeat,14,1,1
*do,i,239,253
n,i+30,3.44,2.77,nz(i-232)
*enddo
e,269,270
*repeat,14,1,1
!13号钢筋
*do,i,284,300
n,i,0.84,2.77,nz(i-278)
*enddo
real,1
e,284,285
*repeat,16,1,1
*do,i,284,300
n,i+17,2.04,2.77,nz(i-278)
*enddo
e,301,302
*repeat,16,1,1
*do,i,284,300
n,i+34,3.24,2.77,nz(i-278)
*enddo
e,318,319
*repeat,16,1,1
!14号钢筋
*do,i,335,353
n,i,1.24,2.77,nz(i-330)
*enddo
real,1
e,335,336
*repeat,18,1,1
*do,i,335,353
n,i+19,2.24,2.77,nz(i-330)
*enddo
e,354,355
*repeat,18,1,1
*do,i,335,353
n,i+38,2.64,2.77,nz(i-330)
*enddo
real,1
e,373,374
*repeat,18,1,1
!15号钢筋
*do,i,392,396
n,i,1.04,2.77,nz(i-388)
*enddo
real,1
e,392,393
*repeat,4,1,1
*do,i,392,396
n,i+5,1.84,2.77,nz(i-388)
*enddo
e,397,398
*repeat,4,1,1
!16号钢筋
*do,i,402,407
n,i,0.64,2.77,nz(i-401)
*enddo
real,1
e,402,403
*repeat,5,1,1
*do,i,402,407
n,i+6,1.44,2.77,nz(i-401)
*enddo
e,408,409
*repeat,5,1,1
!17号钢筋
*do,i,414,420
n,i,1.34,0.285,nz(i-413)
*enddo
real,1
e,414,415
*repeat,6,1,1
*do,i,414,420
n,i+7,2.15,0.285,nz(i-413)
*enddo
e,421,422
*repeat,6,1,1
!18号钢筋
*do,i,428,434
n,i,1.07,0.285,nz(i-427)
*enddo
n,435,1.07,0.285-0.24,nz(8)
*do,i,436,442
n,i,1.88,0.285,nz(i-435)
*enddo
n,443,1.88,0.285-0.24,nz(8)
real,1
e,428,429
*repeat,7,1,1
e,436,437
*repeat,7,1,1
!19号钢筋
*do,i,444,450
n,i,0.8,0.285,nz(i-443)
*enddo
n,451,0.8,0.285-0.24,nz(8)
n,452,0.8,0.285-0.51,nz(9)
real,1
e,444,445
*repeat,8,1,1
*do,i,453,459
n,i,1.61,0.285,nz(i-452)
*enddo
n,460,1.61,0.285-0.24,nz(8)
n,461,1.61,0.285-0.51,nz(9)
e,453,454
*repeat,8,1,1
!20号钢筋
*do,i,462,468
n,i,1.07,0.11,nz(i-461)
*enddo
n,469,1.07,0.11-0.24,nz(8)
n,470,1.07,0.11-0.51,nz(9)
n,471,1.07,0.11-0.86,nz(10)
real,1
e,462,463
*repeat,9,1,1
*do,i,472,478
n,i,1.61,0.11,nz(i-471)
*enddo
n,479,1.61,0.11-0.24,nz(8)
n,480,1.61,0.11-0.51,nz(9)
n,481,1.61,0.11-0.86,nz(10)
real,1
e,472,473
*repeat,9,1,1
!21号钢筋
*do,i,482,488
n,i,1.34,0.11,nz(i-481)
*enddo
n,489,1.34,0.11-0.24,nz(8)
n,490,1.34,0.11-0.51,nz(9)
n,491,1.34,0.11-0.86,nz(10)
n,492,1.34,0.11-1.25,nz(11)
real,1
e,482,483
*repeat,10,1,1
*do,i,493,499
n,i,1.88,0.11,nz(i-492)
*enddo
n,500,1.88,0.11-0.24,nz(8)
n,501,1.88,0.11-0.51,nz(9)
n,502,1.88,0.11-0.86,nz(10)
n,503,1.88,0.11-1.25,nz(11)
real,1
e,493,494
*repeat,10,1,1
!22号钢筋
n,504,1.84,2.77,nz(22)
n,505,1.84,2.77,nz(23)
n,506,1.84,2.77,nz(24)
real,1
e,504,505
e,505,506
!23号钢筋
*do,i,507,510
n,i,1.01,0.285,nz(i-486)
*enddo
real,1
e,507,508
*repeat,3,1,1
*do,i,511,514
n,i,1.41,0.285,nz(i-490)
*enddo
e,511,512
*repeat,3,1,1
*do,i,515,518
n,i,1.81,0.285,nz(i-494)
*enddo
e,515,516
*repeat,3,1,1
*do,i,519,522
n,i,2.21,0.285,nz(i-498)
*enddo
e,519,520
*repeat,3,1,1
!24号钢筋
n,523,0.81,0.285-0.24,nz(20)
*do,i,524,527
n,i,0.81,0.285,nz(i-503)
*enddo
real,1
e,523,524
*repeat,4,1,1
n,528,1.21,0.285-0.24,nz(20)
*do,i,529,532
n,i,1.21,0.285,nz(i-508)
*enddo
real,1
e,528,529
*repeat,4,1,1
n,533,1.61,0.285-0.24,nz(20)
*do,i,534,540
n,i,1.61,0.285,nz(i-513)
*enddo
real,1
e,533,534
*repeat,4,1,1
n,541,2.01,0.285-0.24,nz(20)
*do,i,542,545
n,i,2.01,0.285,nz(i-521)
*enddo
e,541,542
*repeat,4,1,1
!25号钢筋
n,546,1.01,0.11-0.51,nz(19)
n,547,1.01,0.11-0.24,nz(20)
*do,i,548,551
n,i,1.01,0.11,nz(i-527)
*enddo
e,546,547
*repeat,5,1,1
n,552,1.61,0.11-0.51,nz(19)
n,553,1.61,0.11-0.24,nz(20)
*do,i,554,557
n,i,1.61,0.11,nz(i-533)
*enddo
e,552,553
*repeat,5,1,1
n,558,2.21,0.11-0.51,nz(19)
n,559,2.21,0.11-0.24,nz(20)
*do,i,560,563
n,i,2.21,0.11,nz(i-539)
*enddo
e,558,559
*repeat,5,1,1
!26号钢筋
n,564,0.81,0.11-0.86,nz(18)
n,565,0.81,0.11-0.51,nz(19)
n,566,0.81,0.11-0.24,nz(20)
*do,i,567,570
n,i,0.81,0.11,nz(i-546)
*enddo
e,564,565
*repeat,6,1,1
n,571,1.41,0.11-0.86,nz(18)
n,572,1.41,0.11-0.51,nz(19)
n,573,1.41,0.11-0.24,nz(20)
*do,i,574,577
n,i,1.41,0.11,nz(i-553)
*enddo
e,571,572
*repeat,6,1,1
n,578,2.01,0.11-0.86,nz(18)
n,579,2.01,0.11-0.51,nz(19)
n,580,2.01,0.11-0.24,nz(20)
*do,i,581,584
n,i,2.01,0.11,nz(i-560)
*enddo
e,578,579
*repeat,6,1,1
!27号钢筋
n,585,1.21,0.11-1.25,nz(17)
n,586,1.21,0.11-0.86,nz(18)
n,587,1.21,0.11-0.51,nz(19)
n,588,1.21,0.11-0.24,nz(20)
*do,i,589,592
n,i,1.21,0.11,nz(i-568)
*enddo
e,585,586
*repeat,7,1,1
n,593,1.81,0.11-1.25,nz(17)
n,594,1.81,0.11-0.86,nz(18)
n,595,1.81,0.11-0.51,nz(19)
n,596,1.81,0.11-0.24,nz(20)
*do,i,597,600
n,i,1.81,0.11,nz(i-576)
*enddo
e,593,594
*repeat,7,1,1
!**********钢筋对称*******
nsel,u,loc,x,0
nsym,x,1000,all
esel,s,mat,,2
esym,,1000,all
allsel,all
numcmp,node
nsel,u,loc,x,0
nsym,z,3000,all
esel,s,mat,,2
esym,,3000,all
nummrg,node
numcmp,node
gplot
!约束耦合
allsel,all
*do,i,1,23
nsel,s,loc,z,nz(i)
cerig,i,all,all
*enddo
*do,i,25,47
nsel,s,loc,z,nz(i)
cerig,i,all,all
*enddo
nsel,s,loc,z,0
nsel,r,loc,y,0,ny(597)-10
cerig,24,all,all
allsel,all
/solu
!*****约束****************
d,2,uy,0,,,,ux
d,26,uy,0,,,,ux
d,14,ux,0,,,,uy,uz
d,38,ux,0,,,,uy,uz
!****自重*****************
allsel,all
acel,0,10,0
!****二期恒载*************
esel,s,mat,,1
sfbeam,all,,pres,187000
!****中活载***************
sfcum,pres,add
*do,i,26,32
sfbeam,i,,pres,91100*2
*enddo
sfbeam,16,,pres,91100*2
sfbeam,11,,pres,91100*2
sfbeam,12,,pres,91100*2
sfbeam,13,,pres,91100*2
sfbeam,14,,pres,91100*2
*do,i,33,39
sfbeam,i,,pres,91900*2
*enddo
sfbeam,18,,pres,91900*2
!**************************
allsel,all
solve
/post1
esel,s,mat,,1
pldisp
plnsol,s,x

dahair

学习一下！！！谢谢

leading0716

好的，很好的学习模型。

xiaod5285

多谢版主帮助！

星期六

很实用，谢谢

stanomo

发错地儿了吧 这是ABA呢

liudapengsnd

楼主支持你，多发点这方面的资料