finish
/clear
/filename,S-1
/title,S-1
/prep7
!——————————————建立模型N-mm————————————————
!混凝土
fc=60$ft=6
et,1,solid65
mp,ex,1,5.4e4
mp,prxy,1,0.2
mp,dens,1,2.3e-5 !混凝土密度取23kN/m3
r,1
tb,concr,1 !定义混凝土(C50)失效准则
tbdata,,0.35,0.95,ft,-1
tb,miso,1,,10 !多线性等向强化模型
tbpt,,0.0004,fc*0.36
tbpt,,0.0008,fc*0.64
tbpt,,0.0012,fc*0.84
tbpt,,0.0016,fc*0.96
tbpt,,0.002,fc
tbpt,,0.0023,fc*0.9775
tbpt,,0.0026,fc*0.91
tbpt,,0.0029,fc*0.7975
tbpt,,0.0031,fc*0.6975
tbpt,,0.0033,fc*0.5775
!钢绞线并考虑预应力损失
area=274.9
f=1023
E=1.95e5
istrn=f/E
et,2,link10
mp,ex,2,1.95e5
mp,prxy,2,0.3
r,2,area,istrn
tb,bkin,2
tbdata,,1860
!箍筋
et,3,link8
mp,ex,3,1.9e5
mp,prxy,3,0.3
r,3,28.27
tb,bkin,3
tbdata,,335
!纵向钢筋
et,4,link8
mp,ex,4,1.9e5
mp,prxy,4,0.3
r,4,254.47$r,5,113.1
tb,bkin,4
tbdata,,400
!..................................工字梁部分
type,1
real,1
mat,1
*dim,x,,94
x(1)=-300,-260,-220,-180,-140,-100,-60,-20,20,60
x(11)=100,140,180,220,260,300,-300,-260,-220,-180
x(21)=-140,-100,-60,-20,20,60,100,140,180,220
x(31)=260,300,-300,-260,-220,-180,-140,-100,-60,-20
x(41)=20,60,100,140,180,220,260,300,-60,-20
x(51)=20,60,-60,-20,20,60,-60,-20,20,60
x(61)=-60,-20,20,60,-60,-20,20,60,-60,-20
x(71)=20,60,-60,-20,20,60,-100,-60,-20,20
x(81)=60,100,-100,-60,-20,20,60,100,-100,-60
x(91)=-20,20,60,100
*dim,y,,94
y(1)=500,500,500,500,500,500,500,500,500,500
y(11)=500,500,500,500,500,500,450,450,450,450
y(21)=450,450,450,450,450,450,450,450,450,450
y(31)=450,450,400,400,400,400,400,400,400,400
y(41)=400,400,400,400,400,400,400,400,360,360
y(51)=360,360,320,320,320,320,280,280,280,280
y(61)=240,240,240,240,200,200,200,200,160,160
y(71)=160,160,120,120,120,120,80,80,80,80
y(81)=80,80,40,40,40,40,40,40,0,0
y(91)=0,0,0,0
*do,i,1,94,1
n,i,x(i),y(i),285
*enddo
ngen,2,100,1,94,1,,,35
ngen,28,100,101,194,1,,,40
ngen,11,100,2801,2894,1,,,50
*do,i,0,37,1 !控制z向
*do,j,0,14,1 !控制x向
*do,m,0,1,1 !控制y向
e,1+100*i+j+16*m,2+100*i+j+16*m,18+100*i+j+16*m,17+100*i+j+16*m,101+100*i+j+16*m,102+100*i+j+16*m,118+100*i+j+16*m,117+100*i+j+16*m
*enddo
*enddo
*enddo
*do,i,0,37,1
*do,j,0,2,1
e,39+100*i+j,40+100*i+j,50+100*i+j,49+100*i+j,139+100*i+j,140+100*i+j,150+100*i+j,149+100*i+j
*enddo
*enddo
*do,i,0,37,1
*do,j,0,2,1
*do,m,0,5,1
e,49+100*i+j+4*m,50+100*i+j+4*m,54+100*i+j+4*m,53+100*i+j+4*m,149+100*i+j+4*m,150+100*i+j+4*m,154+100*i+j+4*m,153+100*i+j+4*m
*enddo
*enddo
*enddo
*do,i,0,37,1
*do,j,0,2,1
e,73+100*i+j,74+100*i+j,79+100*i+j,78+100*i+j,173+100*i+j,174+100*i+j,179+100*i+j,178+100*i+j
*enddo
*enddo
*do,i,0,37,1
*do,j,0,4,1
*do,m,0,1,1
e,77+100*i+j+6*m,78+100*i+j+6*m,84+100*i+j+6*m,83+100*i+j+6*m,177+100*i+j+6*m,178+100*i+j+6*m,184+100*i+j+6*m,183+100*i+j+6*m
*enddo
*enddo
*enddo
!..............................边横隔板部分
*dim,x1,,208
x1(1)=-300,-260,-220,-180,-140,-100,-60,-20,20,60
x1(11)=100,140,180,220,260,300,-300,-260,-220,-180
x1(21)=-140,-100,-60,-20,20,60,100,140,180,220
x1(31)=260,300,-300,-260,-220,-180,-140,-100,-60,-20
x1(41)=20,60,100,140,180,220,260,300,-300,-260
x1(51)=-220,-180,-140,-100,-60,-20,20,60,100,140
x1(61)=180,220,260,300,-300,-260,-220,-180,-140,-100
x1(71)=-60,-20,20,60,100,140,180,220,260,300
x1(81)=-300,-260,-220,-180,-140,-100,-60,-20,20,60
x1(91)=100,140,180,220,260,300,-300,-260,-220,-180
x1(101)=-140,-100,-60,-20,20,60,100,140,180,220
x1(111)=260,300,-300,-260,-220,-180,-140,-100,-60,-20
x1(121)=20,60,100,140,180,220,260,300,-300,-260
x1(131)=-220,-180,-140,-100,-60,-20,20,60,100,140
x1(141)=180,220,260,300,-300,-260,-220,-180,-140,-100
x1(151)=-60,-20,20,60,100,140,180,220,260,300
x1(161)=-300,-260,-220,-180,-140,-100,-60,-20,20,60
x1(171)=100,140,180,220,260,300,-300,-260,-220,-180
x1(181)=-140,-100,-60,-20,20,60,100,140,180,220
x1(191)=260,300,-300,-260,-220,-180,-140,-100,-60,-20
x1(201)=20,60,100,140,180,220,260,300
*dim,y1,,208
y1(1)=500,500,500,500,500,500,500,500,500,500
y1(11)=500,500,500,500,500,500,450,450,450,450
y1(21)=450,450,450,450,450,450,450,450,450,450
y1(31)=450,450,400,400,400,400,400,400,400,400
y1(41)=400,400,400,400,400,400,400,400,360,360
y1(51)=360,360,360,360,360,360,360,360,360,360
y1(61)=360,360,360,360,320,320,320,320,320,320
y1(71)=320,320,320,320,320,320,320,320,320,320
y1(81)=280,280,280,280,280,280,280,280,280,280
y1(91)=280,280,280,280,280,280,240,240,240,240
y1(101)=240,240,240,240,240,240,240,240,240,240
y1(111)=240,240,200,200,200,200,200,200,200,200
y1(121)=200,200,200,200,200,200,200,200,160,160
y1(131)=160,160,160,160,160,160,160,160,160,160
y1(141)=160,160,160,160,120,120,120,120,120,120
y1(151)=120,120,120,120,120,120,120,120,120,120
y1(161)=80,80,80,80,80,80,80,80,80,80
y1(171)=80,80,80,80,80,80,40,40,40,40
y1(181)=40,40,40,40,40,40,40,40,40,40
y1(191)=40,40,0,0,0,0,0,0,0,0
y1(201)=0,0,0,0,0,0,0,0
*do,i,1,208,1
n,4000+i,x1(i),y1(i)
*enddo
ngen,6,300,4001,4208,1,,,50
*do,i,0,4,1
*do,j,0,14,1
*do,m,0,11,1
e,4001+300*i+j+16*m,4002+300*i+j+16*m,4018+300*i+j+16*m,4017+300*i+j+16*m,4301+300*i+j+16*m,4302+300*i+j+16*m,4318+300*i+j+16*m,4317+300*i+j+16*m
*enddo
*enddo
*enddo
*do,j,0,14,1
*do,m,0,1,1
e,5501+j+16*m,5502+j+16*m,5518+j+16*m,5517+j+16*m,1+j+16*m,2+j+16*m,18+j+16*m,17+j+16*m
*enddo
*enddo
*do,j,0,2,1
e,5539+j,5540+j,5556+j,5555+j,39+j,40+j,50+j,49+j
*enddo
*do,j,0,2,1
*do,m,0,5,1
e,5555+j+16*m,5556+j+16*m,5572+j+16*m,5571+j+16*m,49+j+4*m,50+j+4*m,54+j+4*m,53+j+4*m
*enddo
*enddo
*do,j,0,2,1
e,5651+j,5652+j,5668+j,5667+j,73+j,74+j,79+j,78+j
*enddo
*do,j,0,4,1
*do,m,0,1,1
e,5666+j+16*m,5667+j+16*m,5683+j+16*m,5682+j+16*m,77+j+6*m,78+j+6*m,84+j+6*m,83+j+6*m
*enddo
*enddo
!..............................中横隔板部分
*do,i,1,208,1
n,8000+i,x1(i),y1(i),2000
*enddo
ngen,4,300,8001,8208,1,,,-25
*do,i,0,2,1
*do,j,0,14,1
*do,m,0,11,1
e,8001+300*i+j+16*m,8002+300*i+j+16*m,8018+300*i+j+16*m,8017+300*i+j+16*m,8301+300*i+j+16*m,8302+300*i+j+16*m,8318+300*i+j+16*m,8317+300*i+j+16*m
*enddo
*enddo
*enddo
*do,j,0,14,1
*do,m,0,1,1
e,8901+j+16*m,8902+j+16*m,8918+j+16*m,8917+j+16*m,3801+j+16*m,3802+j+16*m,3818+j+16*m,3817+j+16*m
*enddo
*enddo
*do,j,0,2,1
e,8939+j,8940+j,8956+j,8955+j,3839+j,3840+j,3850+j,3849+j
*enddo
*do,j,0,2,1
*do,m,0,5,1
e,8955+j+16*m,8956+j+16*m,8972+j+16*m,8971+j+16*m,3849+j+4*m,3850+j+4*m,3854+j+4*m,3853+j+4*m
*enddo
*enddo
*do,j,0,2,1
e,9051+j,9052+j,9068+j,9067+j,3873+j,3874+j,3879+j,3878+j
*enddo
*do,j,0,4,1
*do,m,0,1,1
e,9066+j+16*m,9067+j+16*m,9083+j+16*m,9082+j+16*m,3877+j+6*m,3878+j+6*m,3884+j+6*m,3883+j+6*m
*enddo
*enddo
!..............................建立半跨纵筋
type,4
real,5
mat,4
allsel
*do,i,0,3,1 !顶板纵筋
*do,j,0,2,1
e,4318+300*i+3*j,4618+300*i+3*j
*enddo
*enddo
*do,i,0,3,1
*do,j,0,2,1
e,4325+300*i+3*j,4625+300*i+3*j
*enddo
*enddo
*do,j,0,2,1
e,5518+3*j,18+3*j
*enddo
*do,j,0,2,1
e,5525+3*j,25+3*j
*enddo
*do,i,0,37,1
*do,j,0,2,1
e,18+100*i+3*j,118+100*i+3*j
*enddo
*enddo
*do,i,0,37,1
*do,j,0,2,1
e,25+100*i+3*j,125+100*i+3*j
*enddo
*enddo
*do,j,0,2,1
e,3818+3*j,8918+3*j
*enddo
*do,j,0,2,1
e,3825+3*j,8925+3*j
*enddo
*do,i,0,2,1
*do,j,0,2,1
e,8918-300*i+3*j,8618-300*i+3*j
*enddo
*enddo
*do,i,0,2,1
*do,j,0,2,1
e,8925-300*i+3*j,8625-300*i+3*j
*enddo
*enddo
type,4
real,4
mat,4
*do,i,0,4,1 !底板纵筋
*do,j,0,3,1
e,4183+300*i+j,4483+300*i+j
*enddo
*enddo
*do,j,0,3,1
e,5683+j,84+j
*enddo
*do,i,0,37,1
*do,j,0,3,1
e,84+100*i+j,184+100*i+j
*enddo
*enddo
*do,j,0,3,1
e,3884+j,9083+j
*enddo
*do,i,0,2,1
*do,j,0,3,1
e,9083-300*i+j,8783-300*i+j
*enddo
*e