预应力混凝土连续梁模型模型

/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

全桥模型

全桥模型

耦合

耦合

预应力钢筋

预应力钢筋

楼主，想请问下，这个模型跟midas对照过吗？   这种模拟方法验证过吗？。。谢谢！！

太牛了 值得学习啊 做这个模型可是要用很多时间的 加油

挺不错的，值得学习！！！支持！！版主应该给加分。

楼主啊，你也太牛了把，顶你啊，这么好的资料，谢谢分享。

楼主的功底确实不错
但是，我有个疑问
你用ANSYS的梁单元建这么个空间模型，跟midas的结果对比过吗？
我倒觉得这种梁单元分析，桥博和midas足矣了。
我没仔细看，只看到你的预应力是耦合上去的，损失啊，收缩徐变啊，恐怕还midas和桥博的计算好吧？

回复 5# seulqs


    你说的很对
对于这种模型用ANSYS很麻烦，而且由于模型处理不当的话计算结果往往会错误，用MIDAS和桥博则是非常理想的，对于此模型用ANSYS建模时是很麻烦的，而且要计算预应力损失的话，需要另外编写相关程序计算，所以这个模型给大家作为学习ANSYS之用，对于以后工作中计算桥梁模型，肯定是用MIDAS和桥博了，ANSYS仅适用于科学研究。

lz，真的很牛哦，做这个。呵呵~~~~~~~~~~~~

lz,你好，其中有这块我看的不是很明白
!约束耦合
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
这个约束耦合是耦合哪些单元，能不能给我解释下，谢谢了，请加我QQ，932977873，非常希望能和你交流更多。

回复 8# jiangbo1026 的帖子


    是用于耦合主梁节点与钢筋节点，即钢筋节点从属与主梁节点，实现刚臂，将预应力筋的力施加于主梁上，实现预应力效应

回复 9# CHINSHLIA 的帖子


    恩，这想法很好哦，哈哈，先前没想到过，能不能加我QQ932977873，想对于ansys 的有些问题向你请教！！

看楼主所贴的图，似乎将一个范围内的预应力钢束的节点都耦合于一个固定节点，这与实际受力不符，请楼主指教。

耦合节点模拟预应力和混凝土连接，强

楼主，您好！“武汉理工大学桥梁工程09级”，却神牛！！！不得不佩服，希望可以多交流：QQ：1039002445。

神牛！！！不得不佩服

太牛了 值得学习啊 做这个模型可是要用很多时间的 加油