论坛 › C07：ANSYS--实例赏评 › 预应力混凝土模型

chinshlia 发表于 2010-5-10 23:35:14

预应力混凝土模型

本帖最后由 chinshlia 于 2010-5-10 23:53 编辑

本模型为预应力混凝土连三跨续梁桥，采用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 发表于 2010-5-18 14:15:27

学习一下！！！谢谢

leading0716 发表于 2010-6-8 10:08:29

好的，很好的学习模型。

xiaod5285 发表于 2010-7-11 19:29:10

多谢版主帮助！

星期六 发表于 2010-10-11 16:51:33

很实用，谢谢

stanomo 发表于 2010-11-12 10:35:07

发错地儿了吧 这是ABA呢

liudapengsnd 发表于 2011-3-10 11:21:06

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

查看完整版本: 预应力混凝土模型