/nerr,0
!将参数读入数据库
finish
/clear
/PREP7
!单元类型编号
et,1,beam44
et,2,shell63
et,3,link8
!材料类型编号 Q420/Q370/拉索-1，混凝土C40/C50-2，刚臂-3
mp,ex,1,2.1e11
mp,prxy,1,0.3
mp,dens,1,7850
mp,alpx,1,1e-5
mp,ex,2,3.5e7
mp,prxy,2,0.16
mp,dens,2,2500
mp,alpx,2,1e-5
mp,ex,3,1e15
mp,prxy,3,1
mp,dens,3,0
mp,alpx,3,0
!读取自定义截面并编号
*dim,name,char,24
name(1)='Sec_1a'
name(2)='Sec_2a'
name(3)='Sec_3a'
name(4)='Sec_4a'
name(5)='Sec_5a'
name(6)='Sec_6a'
name(7)='Sec_7a'
name(8)='Sec_1'
name(9)='Sec_2'
name(10)='Sec_3'
name(11)='Sec_4'
name(12)='Sec_5'
name(13)='Sec_6'
name(14)='Sec_7'
name(15)='Sec_8'
name(16)='Sec_9'
name(17)='Sec_10'
name(18)='Sec_11'
name(19)='Sec_12'
name(20)='Sec_13'
name(21)='Sec_14'
name(22)='Sec_15'
name(23)='Sec_16'
name(24)='Sec_x2'
*do,i,1,24,1
sectype,i,beam,mesh,
secread,name(i),,,mesh
*enddo
!将ta_sec参数读入数据库
!读取自定义截面并编号
*dim,name2,char,15
name2(1)='ta_1'
name2(2)='ta_2'
name2(3)='ta_3'
name2(4)='ta_4'
name2(5)='ta_5'
name2(6)='ta_6'
name2(7)='ta_7'
name2(8)='ta_8'
name2(9)='ta_9'
name2(10)='ta_10'
name2(11)='ta_11'
name2(12)='ta_12'
name2(13)='ta_13'
name2(14)='ta_14'
*do,i,1,14,1
sectype,40+i,beam,mesh,
secread,name2(i),,,mesh
*enddo
SECTYPE, 26, BEAM, HREC, , 0
SECOFFSET, USER, 0.1, 0.266
SECDATA,0.2,0.266,0.006,0.006,0.006,0.006,0,0,0,0
SECTYPE, 27, BEAM, I, , 0
SECOFFSET, USER, 0, 0.584
SECDATA,0.240,0.240,0.6,0.016,0.016,0.012,0,0,0,0
SECTYPE, 28, BEAM, T, T_nin, 0
SECOFFSET, USER, 0, 1.4
SECDATA,0.6,1.4,0.032,0.018,0,0,0,0,0,0
SECTYPE, 29, BEAM, T, T_max, 0
SECOFFSET, USER, 0, 2.1
SECDATA,0.6,2.1,0.032,0.016,0,0,0,0,0,0
SECTYPE, 31, BEAM, HREC, , 0
SECOFFSET, CENT
SECDATA,0.548,0.500,0.016,0.016,0.016,0.016,0,0,0,0
SECTYPE, 32, BEAM, I, , 0
SECOFFSET, CENT
SECDATA,0.560,0.560,0.500,0.028,0.028,0.016,0,0,0,0
!刚臂的截面
SECTYPE, 100, BEAM, CSOLID, , 0
SECOFFSET, CENT
SECDATA,0.06,0,0,0,0,0,0,0,0,0
r,1,
r,2,0.016,
!拉索截面控制
!7-349
r,3,134.31e-4,
!7-301
r,4,115.84e-4,
!7-283
r,5,108.91e-4,
!7-241
r,6,100.000e-4,
finish
!开始建立模型
/prep7
!*********************************************
!
!
!xiaXianGan
!
!
!*********************************************
*SET,j,-14
*do,i,1,(94*25+1)*2,94*25+1
n,i,0,j,0
n,94*25+i,1363.000,j,0
FILL,i,94*25+i,94*25-1,,1,,,
!中桁
*if,i,eq,94*25+2,then
!7截面
type,1
mat,1
secnum,14
*do,a,1,3*25,1
E,i-1+a,i+a
*enddo
*do,b,1,11*25,1
E,35*25-1+i+b,35*25+i+b
*enddo
*do,c,1,8*25,1
E,77*25-1+i+c,77*25+i+c
*enddo
!6截面
type,1
mat,1
secnum,13
*do,a,1,3*25,1
E,i+3*25-1+a,i+3*25+a
*enddo
*do,b,1,3*25,1
E,32*25-1+i+b,32*25+i+b
*enddo
*do,c,1,3*25,1
E,46*25-1+i+c,46*25+i+c
*enddo
*do,d,1,1*25,1
E,i+76*25-1+d,i+76*25+d
*enddo
*do,
e,1,7*25,1
E,87*25-1+i+e,87*25+i+e
*enddo
!5截面
type,1
mat,1
secnum,12
*do,a,1,25,1
E,8*25-1+i+a,8*25+i+a
*enddo
*do,b,1,2*25,1
E,30*25-1+i+b,30*25+i+b
*enddo
*do,c,1,2*25,1
E,49*25-1+i+c,49*25+i+c
*enddo
*do,d,1,2*25,1
E,72*25-1+i+d,72*25+i+d
*enddo
!4截面
type,1
mat,1
secnum,11
*do,a,1,3*25,1
E,9*25-1+i+a,9*25+i+a
*enddo
*do,b,1,4*25,1
E,26*25-1+i+b,26*25+i+b
*enddo
*do,c,1,3*25,1
E,51*25-1+i+c,51*25+i+c
*enddo
*do,d,1,4*25,1
E,68*25-1+i+d,68*25+i+d
*enddo
!3截面
type,1
mat,1
secnum,10
*do,a,1,1*25,1
E,7*25-1+i+a,7*25+i+a
*enddo
*do,b,1,5*25,1
E,12*25-1+i+b,12*25+i+b
*enddo
*do,c,1,4*25,1
E,22*25-1+i+c,22*25+i+c
*enddo
*do,d,1,4*25,1
E,54*25-1+i+d,54*25+i+d
*enddo
*do,e,1,5*25,1
E,63*25-1+i+e,63*25+i+e
*enddo
!2截面
type,1
mat,1
secnum,9
*do,a,1,1*25,1
E,6*25-1+i+a,6*25+i+a
*enddo
*do,b,1,2*25,1
E,17*25-1+i+b,17*25+i+b
*enddo
*do,c,1,1*25,1
E,21*25-1+i+c,21*25+i+c
*enddo
*do,d,1,1*25,1
E,58*25-1+i+d,58*25+i+d
*enddo
*do,e,1,2*25,1
E,61*25-1+i+e,61*25+i+e
*enddo
*do,f,1,2*25,1
E,74*25-1+i+f,74*25+i+f
*enddo
*do,g,1,1*25,1
E,85*25-1+i+g,85*25+i+g
*enddo
!1截面
type,1
mat,1
secnum,8
*do,a,1,2*25,1
E,19*25-1+i+a,19*25+i+a
*enddo
*do,b,1,2*25,1
E,59*25-1+i+b,59*25+i+b
*enddo
*do,c,1,1*25,1
E,86*25-1+i+c,86*25+i+c
*enddo
!左边桁
*elseif,i,eq,(94*25+1)*2+1,then
!7截面
type,1
mat,1
secnum,7
*do,a,1,3*25,1
E,i+a,i-1+a
*enddo
*do,b,1,11*25,1
E,35*25+i+b,35*25-1+i+b
*enddo
*do,c,1,8*25,1
E,77*25+i+c,77*25-1+i+c
*enddo
!6截面
type,1
mat,1
secnum,6
*do,a,1,3*25,1
E,i+3*25+a,i+3*25-1+a
*enddo
*do,b,1,3*25,1
E,32*25+i+b,32*25-1+i+b
*enddo
*do,c,1,3*25,1
E,46*25+i+c,46*25-1+i+c
*enddo
*do,d,1,1*25,1
E,i+76*25+d,i+76*25-1+d
*enddo
*do,e,1,7*25,1
E,87*25+i+e,87*25-1+i+e
*enddo
!5截面
type,1
mat,1
secnum,5
*do,a,1,25,1
E,8*25+i+a,8*25-1+i+a
*enddo
*do,b,1,2*25,1
E,30*25+i+b,30*25-1+i+b
*enddo
*do,c,1,2*25,1
E,49*25+i+c,49*25-1+i+c
*enddo
*do,d,1,2*25,1
E,72*25+i+d,72*25-1+i+d
*enddo
!4截面
type,1
mat,1
secnum,4
*do,a,1,3*25,1
E,9*25+i+a,9*25-1+i+a
*enddo
*do,b,1,4*25,1
E,26*25+i+b,26*25-1+i+b
*enddo
*do,c,1,3*25,1
E,51*25+i+c,51*25-1+i+c
*enddo
*do,d,1,4*25,1
E,68*25+i+d,68*25-1+i+d
*enddo
!3截面
type,1
mat,1
secnum,3
*do,a,1,1*25,1
E,7*25+i+a,7*25-1+i+a
*enddo
*do,b,1,5*25,1
E,12*25+i+b,12*25-1+i+b
*enddo
*do,c,1,4*25,1
E,22*25+i+c,22*25-1+i+c
*enddo
*do,d,1,4*25,1
E,54*25+i+d,54*25-1+i+d
*enddo
*do,e,1,5*25,1
E,63*25+i+e,63*25-1+i+e
*enddo
!2截面
type,1
mat,1
secnum,2
*do,a,1,1*25,1
E,6*25+i+a,6*25-1+i+a
*enddo
*do,b,1,2*25,1
E,17*25+i+b,17*25-1+i+b
*enddo
*do,c,1,1*25,1
E,21*25+i+c,21*25-1+i+c
*enddo
*do,d,1,1*25,1
E,58*25+i+d,58*25-1+i+d
*enddo
*do,e,1,2*25,1
E,61*25+i+e,61*25-1+i+e
*enddo
*do,f,1,2*25,1
E,74*25+i+f,74*25-1+i+f
*enddo
*do,g,1,1*25,1
E,85*25+i+g,85*25-1+i+g
*enddo
!1截面
type,1
mat,1
secnum,1
*do,a,1,2*25,1
E,19*25+i+a,19*25-1+i+a
*enddo
*do,b,1,2*25,1
E,59*25+i+b,59*25-1+i+b
*enddo
*do,c,1,1*25,1
E,86*25+i+c,86*25-1+i+c
*enddo
!右边桁
*else
!7截面
type,1
mat,1
secnum,7
*do,a,1,3*25,1
E,i-1+a,i+a
*enddo
*do,b,1,11*25,1
E,35*25-1+i+b,35*25+i+b
*enddo
*do,c,1,8*25,1
E,77*25-1+i+c,77*25+i+c
*enddo
!6截面
type,1
mat,1
secnum,6
*do,a,1,3*25,1
E,i+3*25-1+a,i+3*25+a
*enddo
*do,b,1,3*25,1
E,32*25-1+i+b,32*25+i+b
*enddo
*do,c,1,3*25,1
E,46*25-1+i+c,46*25+i+c
*enddo
*do,d,1,1*25,1
E,i+76*25-1+d,i+76*25+d
*enddo
*do,e,1,7*25,1
E,87*25-1+i+e,87*25+i+e
*enddo
!5截面
type,1
mat,1
secnum,5
*do,a,1,25,1
E,8*25-1+i+a,8*25+i+a
*enddo
*do,b,1,2*25,1
E,30*25-1+i+b,30*25+i+b
*enddo
*do,c,1,2*25,1
E,49*25-1+i+c,49*25+i+c
*enddo
*do,d,1,2*25,1
E,72*25-1+i+d,72*25+i+d
*enddo
!4截面
type,1
mat,1
secnum,4
*do,a,1,3*25,1
E,9*25-1+i+a,9*25+i+a
*enddo
*do,b,1,4*25,1
E,26*25-1+i+b,26*25+i+b
*enddo
*do,c,1,3*25,1
E,51*25-1+i+c,51*25+i+c
*enddo
*do,d,1,4*25,1
E,68*25-1+i+d,68*25+i+d
*enddo
!3截面
type,1
mat,1
secnum,3
*do,a,1,1*25,1
E,7*25-1+i+a,7*25+i+a
*enddo
*do,b,1,5*25,1
E,12*25-1+i+b,12*25+i+b
*enddo
*do,c,1,4*25,1
E,22*25-1+i+c,22*25+i+c
*enddo
*do,d,1,4*25,1
E,54*25-1+i+d,54*25+i+d
*enddo
*do,e,1,5*25,1
E,63*25-1+i+e,63*25+i+e
*enddo
!2截面
type,1
mat,1
secnum,2
*do,a,1,1*25,1
E,6*25-1+i+a,6*25+i+a
*enddo
*do,b,1,2*25,1
E,17*25-1+i+b,17*25+i+b
*enddo
*do,c,1,1*25,1
E,21*25-1+i+c,21*25+i+c
*enddo
*do,d,1,1*25,1
E,58*25-1+i+d,58*25+i+d
*enddo
*do,e,1,2*25,1
E,61*25-1+i+e,61*25+i+e
*enddo
*do,f,1,2*25,1
E,74*25-1+i+f,74*25+i+f
*enddo
*do,g,1,1*25,1
E,85*25-1+i+g,85*25+i+g
*enddo
!1截面
type,1
mat,1
secnum,1
*do,a,1,2*25,1
E,19*25-1+i+a,19*25+i+a
*enddo
*do,b,1,2*25,1
E,59*25-1+i+b,59*25+i+b
*enddo
*do,c,1,1*25,1
E,86*25-1+i+c,86*25+i+c
*enddo
*endif
*SET,j,j+14
*enddo
!*********************************************
!
!
!T_beam AND shell63's nodes
!
!
!*********************************************
/prep7
ngen,2,7053,1,2351,1,0,0,0.7,
ngen,4,2351,7054,7054+2350,1,0,0.57,0,
ngen,4,2351,7054+3*2351,7054+3*2351+2350,1,0,0.60,0,
ngen,3,2351,7054+6*2351,7054+6*2351+2350,1,0,0.45,0,
ngen,2,2351,7054+8*2351,7054+8*2351+2350,1,0,0.60,0,
ngen,3,2351,7054+9*2351,7054+9*2351+2350,1,0,0.45,0,
ngen,2,2351,7054+11*2351,7054+11*2351+2350,1,0,0.70,0,
local,11,0,0,-7,0
nsym,y,2351*12,7054+2351,32914+2*2351,1
csys,0
nsym,y,58775,7054,65828,1
ngen,2,122252,2352,2352+2350,1,0,0,0.7,
!生成 边上 4条 U肋
type,1
mat,1
secnum,26
*do,h,6,9,1
*do,k,1,2350,1
e,2351*h+k,2351*h+1+k
*end
do
*enddo
!生成T梁1
type,1
mat,1
secnum,27
*do,k,1,2350,1
e,2351*10+k,2351*10+1+k
*enddo
!生成 U肋 2条
type,1
mat,1
secnum,26
*do,h,11,12,1
*do,k,1,2350,1
e,2351*h+k,2351*h+1+k
*enddo
*enddo
!生成T梁2
type,1
mat,1
secnum,27
*do,k,1,2350,1
e,2351*13+k,2351*13+1+k
*enddo
!生成 U肋 另外2条
type,1
mat,1
secnum,26
*do,h,14,15,1
*do,k,1,2350,1
e,2351*h+k,2351*h+1+k
*enddo
*enddo
!镜像四分之一纵向的单元
csys,11
!（最好用*get 最大节点号单元号）
esym,,2351*12,7051,30550,1
!生成横梁
type,1
mat,1
secnum,29
*do,h,7054,9404,25
*do,k,h,h+11*2351,2351
e,k,k+2351
*enddo
*enddo
*do,h,37617,39967,25
*do,k,h,h+10*2351,2351
e,k,k+2351
*enddo
*enddo
*do,h,35266,37616,25
e,h,h+2351*12
*enddo
!生成横肋
type,1
mat,1
secnum,28
*do,ah,7059,9384,25
*do,ac,0,15,5
*do,ak,ah+ac,ah+ac+11*2351,2351
e,ak,ak+2351
*enddo
*enddo
*enddo
*do,ah,37622,39947,25
*do,ac,0,15,5
*do,ak,ah+ac,ah+ac+10*2351,2351
e,ak,ak+2351
*enddo
*enddo
*enddo
*do,ah,35271,37596,25
*do,ac,0,15,5
e,ah+ac,ah+ac+2351*12
*enddo
*enddo
!镜像一半单元
csys,0
esym,,58775,7051,65354,1
!补充横梁
type,1
mat,1
secnum,29
*do,h,37617,39967,25
e,h,h+86987
*enddo
*do,h,96392,98742,25
e,h,h+28212
*enddo
!补充横肋
type,1
mat,1
secnum,28
*do,ah,37622,39947,25
*do,ac,0,15,5
e,ah+ac,ah+ac+86987
*enddo
*enddo
*do,ah,96397,98722,25
*do,ac,0,15,5
e,ah+ac,ah+ac+28212
*enddo
*enddo
!刚臂
type,1
mat,3
secnum,100
*do,h,1,2351,1
e,h,h+7053
*enddo
*do,h,4703,7053,1
e,h,h+61126
*enddo
*do,h,2352,4702,1
e,h,h+124604-2352
*enddo
!*********************************************
!
!
!shell
!
!
!*********************************************
type,2
mat,1
real,2
*do,i,0,11,1
*do,h,7054,7054+2349,1
e,h+i*2351,h+1+i*2351,h+2352+i*2351,h+2351+i*2351
*enddo
*enddo
*do,h,35266,35266+2349,1
e,h,h+1,h+63479-35266,h+63478-35266
*enddo
*do,i,0,-10,-1
*do,h,63478,63478+2349,1
e,h+i*2351,h+1+i*2351,h+i*2351-2350,h+i*2351-2351
*enddo
*enddo
*do,h,37617,37617+2349,1
e,h,h+1,h+1+124604-37617,h+124604-37617
*enddo
csys,0
esym,,58775,131654,190403,1
*do,h,124604,124604+2349,1
e,h,h+1,h+1-(124604-96392),h-(124604-96392)
*enddo
!*********************************************
!
!
!ShangXianGan
!
!
!*********************************************
/prep7
real,1
*SET,Nmax1,300000
*SET,j,-14
*do,i,Nmax1+1,Nmax1+(92*25+1)*2,92*25+1
n,i,14.50,j,15.00
n,92*25+i,1348.50,j,15.00
FILL,i,92*25+i,,,,,,
!中桁
*if,i,eq,Nmax1+92*25+2,then
!13截面
type,1
mat,1
secnum,20
*do,a,1,2*25,1
E,i-1+a,i+a
*enddo
*do,b,1,6*25,1
E,3*25-1+i+b,3*25+i+b
*enddo
*do,c,1,3*25,1
E,70*25-1+i+c,70*25+i+c
*enddo
*do,d,1,10*25,1
E,74*25-1+i+d,74*25+i+d
*enddo
*do,e,1,3*25,1
E,85*25-1+i+e,85*25+i+e
*enddo
*do,f,1,1*25,1
E,91*25-
1+i+f,91*25+i+f
*enddo
!12截面
type,1
mat,1
secnum,19
*do,a,1,1*25,1
E,i+2*25-1+a,i+2*25+a
*enddo
*do,b,1,1*25,1
E,9*25-1+i+b,9*25+i+b
*enddo
*do,c,1,8*25,1
E,35*25-1+i+c,35*25+i+c
*enddo
*do,d,1,1*25,1
E,i+69*25-1+d,i+69*25+d
*enddo
*do,e,1,1*25,1
E,73*25-1+i+e,73*25+i+e
*enddo
*do,f,1,3*25,1
E,88*25-1+i+f,88*25+i+f
*enddo
!11截面
type,1
mat,1
secnum,18
*do,a,1,25,1
E,10*25-1+i+a,10*25+i+a
*enddo
*do,b,1,1*25,1
E,18*25-1+i+b,18*25+i+b
*enddo
*do,c,1,1*25,1
E,26*25-1+i+c,26*25+i+c
*enddo
*do,d,1,4*25,1
E,31*25-1+i+d,31*25+i+d
*enddo
*do,e,1,1*25,1
E,43*25-1+i+e,43*25+i+e
*enddo
*do,f,1,5*25,1
E,45*25-1+i+f,45*25+i+f
*enddo
*do,g,1,1*25,1
E,59*25-1+i+g,59*25+i+g
*enddo
*do,h,1,1*25,1
E,68*25-1+i+h,68*25+i+h
*enddo
!10截面
type,1
mat,1
secnum,17
*do,a,1,1*25,1
E,11*25-1+i+a,11*25+i+a
*enddo
*do,b,1,1*25,1
E,17*25-1+i+b,17*25+i+b
*enddo
*do,c,1,1*25,1
E,20*25-1+i+c,20*25+i+c
*enddo
*do,d,1,1*25,1
E,25*25-1+i+d,25*25+i+d
*enddo
*do,e,1,4*25,1
E,27*25-1+i+e,27*25+i+e
*enddo
*do,f,1,3*25,1
E,50*25-1+i+f,50*25+i+f
*enddo
*do,g,1,1*25,1
E,57*25-1+i+g,57*25+i+g
*enddo
*do,h,1,1*25,1
E,60*25-1+i+h,60*25+i+h
*enddo
*do,hh,1,1*25,1
E,67*25-1+i+hh,67*25+i+hh
*enddo
!9截面
type,1
mat,1
secnum,16
*do,a,1,5*25,1
E,12*25-1+i+a,12*25+i+a
*enddo
*do,b,1,1*25,1
E,19*25-1+i+b,19*25+i+b
*enddo
*do,c,1,4*25,1
E,21*25-1+i+c,21*25+i+c
*enddo
*do,d,1,4*25,1
E,53*25-1+i+d,53*25+i+d
*enddo
*do,e,1,1*25,1
E,58*25-1+i+e,58*25+i+e
*enddo
*do,g,1,6*25,1
E,61*25-1+i+g,61*25+i+g
*enddo
*do,h,1,1*25,1
E,84*25-1+i+h,84*25+i+h
*enddo
!8截面
type,1
mat,1
secnum,15
*do,a,1,1*25,1
E,44*25-1+i+a,44*25+i+a
*enddo
!左边桁
*elseif,i,eq,Nmax1+(92*25+1)*2+1,then
!13截面
type,1
mat,1
secnum,20
*do,a,1,2*25,1
E,i-1+a,i+a
*enddo
*do,b,1,6*25,1
E,3*25-1+i+b,3*25+i+b
*enddo
*do,c,1,3*25,1
E,70*25-1+i+c,70*25+i+c
*enddo
*do,d,1,10*25,1
E,74*25-1+i+d,74*25+i+d
*enddo
*do,e,1,3*25,1
E,85*25-1+i+e,85*25+i+e
*enddo
*do,f,1,1*25,1
E,91*25-1+i+f,91*25+i+f
*enddo
!12截面
type,1
mat,1
secnum,19
*do,a,1,1*25,1
E,i+2*25-1+a,i+2*25+a
*enddo
*do,b,1,1*25,1
E,9*25-1+i+b,9*25+i+b
*enddo
*do,c,1,8*25,1
E,35*25-1+i+c,35*25+i+c
*enddo
*do,d,1,1*25,1
E,i+69*25-1+d,i+69*25+d
*enddo
*do,e,1,1*25,1
E,73*25-1+i+e,73*25+i+e
*enddo
*do,f,1,3*25,1
E,88*25-1+i+f,88*25+i+f
*enddo
!11截面
type,1
mat,1
secnum,18
*do,a,1,25,1
E,10*25-1+i+a,10*25+i+a
*enddo
*do,b,1,1*25,1
E,18*25-1+i+b,18*25+i+b
*enddo
*do,c,1,1*25,1
E,26*25-1+i+c,26*25+i+c
*enddo
*do,d,1,4*25,1
E,31*25-1+i+d,31*25+i+d
*enddo
*do,e,1,1*25,1
E,43*25-1+i+e,43*25+i+e
*enddo
*do,f,1,5*25,1
E,45*25-1+i+f,45*25+i+f
*enddo
*do,g,1,1*25,1
E,59*25-1+i+g,59*25+i+g
*enddo
*do,h,1,1*25,1
E,68*25-1+i+h,68*25+i+h
*enddo
!10截面
type,1
mat,1
secnum,17
*do,a,1,1*25,1
E,11*25-1+i+a,11*25+i+a
*enddo
*do,b,1,1*25,1
E,17*25-1+i+b,17*25+i+b
*enddo
*do,c,1,1*25,1
E,20*25-1+i+c,20*25+i+c
*enddo
*do,d,1,1*25,1
E,25*25-1+i+d,25*25+i+d
*enddo
*do,e,1,4*25,1
E,27*25-1+i+e,27*25+i+e
*enddo
*do,f,1,3*25,1
E,50*25-1+i+f,50*25+i+f
*enddo
*do,g,1,1*25,1
E,57*25-1+i+g,57*25+i+g
*enddo
*do,h,1,1*25,1
E,60*25-1+i+h,60*25+i+h
*enddo
*do,hh,1,1*25,1
E,67*25-1+i+hh,67*25+i+hh
*enddo
!9截面
type,1
mat,1
secnum,16
*do,a,1,5*25,1
E,12*25-1+i+a,12*25+i+a
*enddo
*do,b,1,1*25,1
E,19*25-1+i+b,19*25+i+b
*enddo
*do,c,1,4*25,1
E,21*25-1+i+c,21*25+i+c
*enddo
*do,d,1,4*25,1
E,53*25-1+i+d,53*25+i+d
*enddo
*do,e,1,1*25,1
E,58*25-1+i+e,58*25+i+e
*enddo
*do,g,1,6*25,1
E,61*25-1+i+g,61*25+i+g
*enddo
*do,h,1,1*25,1
E,84*25-1+i+h,84*25+i+h
*enddo
!8截面
type,1
mat,1
secnum,15
*do,a,1,1*25,1
E,44*25-1+i+a,44*25+i+a
*enddo
!右边桁
*else
!13截面
type,1
mat,1
secnum,20
*do,a,1,2*25,1
E,i-1+a,i+a
*enddo
*do,b,1,6*25,1
E,3*25-1+i+b,3*25+i+b
*enddo
*do,c,1,3*25,1
E,70*25-1+i+c,70*25+i+c
*enddo
*do,d,1,10*25,1
E,74*25-1+i+d,74*25+i+d
*enddo
*do,e,1,3*25,1
E,85*25-1+i+e,85*25+i+e
*enddo
*do,f,1,1*25,1
E,91*25-1+i+f,91*25+i+f
*enddo
!12截面
type,1
mat,1
secnum,19
*do,a,1,1*25,1
E,i+2*25-1+a,i+2*25+a
*enddo
*do,b,1,1*25,1
E,9*25-1+i+b,9*25+i+b
*enddo
*do,c,1,8*25,1
E,35*25-1+i+c,35*25+i+c
*enddo
*do,d,1,1*25,1
E,i+69*25-1+d,i+69*25+d
*enddo
*do,e,1,1*25,1
E,73*25-1+i+e,73*25+i+e
*enddo
*do,f,1,3*25,1
E,88*25-1+i+f,88*25+i+f
*enddo
!11截面
type,1
mat,1
secnum,18
*do,a,1,25,1
E,10*25-1+i+a,10*25+i+a
*enddo
*do,b,1,1*25,1
E,18*25-1+i+b,18*25+i+b
*enddo
*do,c,1,1*25,1
E,26*25-1+i+c,26*25+i+c
*enddo
*do,d,1,4*25,1
E,31*25-1+i+d,31*25+i+d
*enddo
*do,e,1,1*25,1
E,43*25-1+i+e,43*25+i+e
*enddo
*do,f,1,5*25,1
E,45*25-1+i+f,45*25+i+f
*enddo
*do,g,1,1*25,1
E,59*25-1+i+g,59*25+i+g
*enddo
*do,h,1,1*25,1
E,68*25-1+i+h,68*25+i+h
*enddo
!10截面
type,1
mat,1
secnum,17
*do,a,1,1*25,1
E,11*25-1+i+a,11*25+i+a
*enddo
*do,b,1,1*25,1
E,17*25-1+i+b,17*25+i+b
*enddo
*do,c,1,1*25,1
E,20*25-1+i+c,20*25+i+c
*enddo
*do,d,1,1*25,1
E,25*25-1+i+d,25*25+i+d
*enddo
*do,e,1,4*25,1
E,27*25-1+i+e,27*25+i+e
*enddo
*do,f,1,3*25,1
E,50*25-1+i+f,50*25+i+f
*enddo
*do,g,1,1*25,1
E,57*25-1+i+g,57*25+i+g
*enddo
*do,h,1,1*25,1
E,60*25-1+i+h,60*25+i+h
*enddo
*do,hh,1,1*25,1
E,67*25-1+i+hh,67*25+i+hh
*enddo
!9截面
type,1
mat,1
secnum,16
*do,a,1,5*25,1
E,12*25-1+i+a,12*25+i+a
*enddo
*do,b,1,1*25,1
E,19*25-1+i+b,19*25+i+b
*enddo
*do,c,1,4*25,1
E,21*25-1+i+c,21*25+i+c
*enddo
*do,d,1,4*25,1
E,53*25-1+i+d,53*25+i+d
*enddo
*do,e,1,1*25,1
E,58*25-1+i+e,58*25+i+e
*enddo
*do,g,1,6*25,1
E,61*25-1+i+g,61*25+i+g
*end
do
*do,h,1,1*25,1
E,84*25-1+i+h,84*25+i+h
*enddo
!8截面
type,1
mat,1
secnum,15
*do,a,1,1*25,1
E,44*25-1+i+a,44*25+i+a
*enddo
*endif
*SET,j,j+14
*enddo
!*********************************************
!
!
!FuGan
!
!
!*********************************************
/prep7
*SET,j,-14
*do,i,1,(94*25+1)*2,94*25+1
k,i,0,j,0
k,94*25+i,1363.000,j,0
kFILL,i,94*25+i,93,,,,,
j=j+14
*enddo
*SET,Nmax1,10000
*SET,j,-14
*do,i,Nmax1+1,Nmax1+(92*25+1)*2,92*25+1
k,i,14.50,j,15.00
k,92*25+i,1348.50,j,15.00
kFILL,i,92*25+i,91,,,,,
j=j+14
*enddo
!连接右斜线1
j=Nmax1+1
*do,n,1,(94*25+1)*2,94*25+1
kk=0
*do,i,n,n-1+20*25,25
lstr,i,j
kk=kk+1
*if,kk,eq,1,then
latt,1,1,1,,,,21
esize,,25
lmesh,all
*elseif,kk,ge,4,and,kk,le,6,then
latt,1,1,1,,,,21
esize,,25
lmesh,all
*elseif,kk,ge,8,and,kk,le,9,then
latt,1,1,1,,,,21
esize,,25
lmesh,all
*elseif,kk,eq,7,then
latt,1,1,1,,,,23
esize,,25
lmesh,all
*else
latt,1,1,1,,,,23
esize,,25
lmesh,all
*endif
j=j+25
*enddo
j=j+2301-(20-0)*25
*enddo
!连接右斜线2
j=Nmax1+1+40*25
*do,n,1,(94*25+1)*2,94*25+1
kk=0
*do,i,n+40*25,n-1+60*25,25
lstr,i,j
kk=kk+1
*if,kk,eq,20,then
latt,1,1,1,,,,21
esize,,25
lmesh,all
*elseif,kk,ge,3,and,kk,le,9,then
latt,1,1,1,,,,21
esize,,25
lmesh,all
*else
latt,1,1,1,,,,22
esize,,25
lmesh,all
*endif
j=j+25
*enddo
j=j+2301-(60-40)*25
*enddo
!连接左斜线1
j=Nmax1+19*25+1
*do,n,1,(94*25+1)*2,94*25+1
kk=0
*do,i,n+21*25,n-1+41*25,25
lstr,i,j
kk=kk+1
*if,kk,eq,1,then
latt,1,1,1,,,,21
esize,,25
lmesh,all
*elseif,kk,ge,13,and,kk,le,18,then
latt,1,1,1,,,,21
esize,,25
lmesh,all
*else
latt,1,1,1,,,,22
esize,,25
lmesh,all
*endif
j=j+25
*enddo
j=j+2301-(41-21)*25
*enddo
!连接左斜线2
j=Nmax1+59*25+1
*do,n,1,(94*25+1)*2,94*25+1
kk=0
*do,i,n+61*25,n-1+95*25,25
lstr,i,j
kk=kk+1
*if,kk,ge,1,and,kk,le,12,then
latt,1,1,1,,,,22
esize,,25
lmesh,all
*elseif,kk,ge,21,and,kk,le,22,then
latt,1,1,1,,,,22
esize,,25
lmesh,all
*elseif,kk,ge,15,and,kk,le,16,then
latt,1,1,1,,,,23
esize,,25
lmesh,all
*elseif,kk,ge,24,and,kk,le,28,then
latt,1,1,1,,,,23
esize,,25
lmesh,all
*elseif,kk,ge,33,and,kk,le,34,then
latt,1,1,1,,,,23
esize,,25
lmesh,all
*else
latt,1,1,1,,,,21
esize,,25
lmesh,all
*endif
j=j+25
*enddo
j=j+2301-(95-61)*25
*enddo
!连接shu线
j=Nmax1+1
*do,n,1,(94*25+1)*2,94*25+1
kk=0
*do,i,n+1*25,n-1+94*25,25
lstr,i,j
kk=kk+1
*if,kk,ge,1,and,kk,le,3,then
latt,1,1,1,,,,22
esize,,25
lmesh,all
*elseif,kk,ge,22,and,kk,le,58,then
latt,1,1,1,,,,22
esize,,25
lmesh,all
*elseif,kk,ge,62,and,kk,le,74,then
latt,1,1,1,,,,22
esize,,25
lmesh,all
*elseif,kk,ge,78,and,kk,le,83,then
latt,1,1,1,,,,22
esize,,25
lmesh,all
*elseif,kk,ge,89,and,kk,le,93,then
latt,1,1,1,,,,22
esize,,25
lmesh,all
*elseif,kk,eq,7,then
latt,1,,1,,,,23
esize,,25
lmesh,all
*elseif,kk,eq,60,then
latt,1,1,1,,,,23
esize,,25
lmesh,all
*elseif,kk,eq,75,then
latt,1,1,1,,
,,23
esize,,25
lmesh,all
*elseif,kk,eq,86,then
latt,1,1,1,,,,23
esize,,25
lmesh,all
*elseif,kk,eq,87,then
latt,1,1,1,,,,23
esize,,25
lmesh,all
*else
latt,1,1,1,,,,21
esize,,25
lmesh,all
*endif
j=j+25
*enddo
j=j+2301-93*25
*enddo
modmsh,detack
ldele,all
kdele,all
!*********************************************
!
!
!2X+上平联
!
!
!*********************************************
/prep7
k,1,14.5,-14,15.00
k,2,14.5,-7,15.00
k,3,14.5,-14,10.2
k,4,14.5,-7,10.2
lstr,1,4
lstr,2,3
lstr,3,4
lgen,2,all,,,,7,
nummrg,kp
numcmp,kp
lstr,2,4
lgen,93,all,,,14.5,,
lgen,2,all,,,,14,,
latt,1,,1,,,,24
esize,,10
lmesh,all
modmsh,detack
ldele,all
kdele,all
*SET,j,-14
*do,i,1,(92*25+1)*2,92*25+1
k,i,14.50,j,15.00
k,92*25+i,1348.50,j,15.00
kFILL,i,92*25+i,91,,,,,
j=j+14
*enddo
*do,i,1,92*25+1,25
lstr,i,i+92*25+1
lstr,i+92*25+1,i+(92*25+1)*2
*enddo
latt,1,,1,,,,31
esize,,20
lmesh,all
modmsh,detack
ldele,all
*do,i,1,92*25,25
lstr,i,i+92*25+1+1*25
lstr,i+92*25+1,i+(92*25+1)*2+1*25
*enddo
*do,i,1+1*25,92*25+1,25
lstr,i,i+92*25+1-25
lstr,i+92*25+1,i+(92*25+1)*2-25
*enddo
nummrg,kp
numcmp,kp
latt,1,,1,,,,32
esize,,20
lmesh,all
modmsh,detack
ldele,all
kdele,all
!************************************************!
!
!1#ta
!
!
!************************************************
/prep7
nn=400000
n,1+nn,290,0,158.9
n,13+nn,290,0,132.5
n,20+nn,290,0,115
n,22+nn,290,0,110.5
n,63+nn,290,0,110.5
n,52+nn,290,-19,-6
n,62+nn,290,-11.5,-34.5
n,93+nn,290,19,-6
n,103+nn,290,11.5,-34.5
n,104+nn,290,-19,-6
n,107+nn,290,19,-6
fill,1+nn,13+nn,11,
fill,13+nn,20+nn,6,
fill,20+nn,22+nn,1,
fill,22+nn,52+nn,29,
fill,52+nn,62+nn,9,
fill,63+nn,93+nn,29,
fill,93+nn,103+nn,9,
fill,104+nn,107+nn,2,
real,1
type,1
mat,2
secnum,53
*do,a,1+nn,12+nn,1
E,a,a+1
*enddo
type,1
mat,2
secnum,52
*do,b,13+nn,19+nn,1
E,b,b+1
*enddo
type,1
mat,2
secnum,51
*do,b,20+nn,21+nn,1
E,b,b+1
*enddo
type,1
mat,2
secnum,46
*do,b,22+nn,51+nn,1
E,b,b+1
*enddo
type,1
mat,2
secnum,42
*do,b,52+nn,61+nn,1
E,b,b+1
*enddo
type,1
mat,2
secnum,46
*do,b,63+nn,92+nn,1
E,b,b+1
*enddo
type,1
mat,2
secnum,42
*do,b,93+nn,102+nn,1
E,b,b+1
*enddo
type,1
mat,2
secnum,54
*do,b,104+nn,106+nn,1
E,b,b+1
*enddo
!************************************************
!
!2#ta
!
!*************************************************
nn=450000
n,1+nn,870,0,158.9
n,13+nn,870,0,132.5
n,20+nn,870,0,115
n,22+nn,870,0,110.5
n,63+nn,870,0,110.5
n,52+nn,870,-19,-6
n,62+nn,870,-11.5,-34.5
n,93+nn,870,19,-6
n,103+nn,870,11.5,-34.5
n,104+nn,870,-19,-6
n,107+nn,870,19,-6
fill,1+nn,13+nn,11,
fill,13+nn,20+nn,6,
fill,20+nn,22+nn,1,
fill,22+nn,52+nn,29,
fill,52+nn,62+nn,9,
fill,63+nn,93+nn,29,
fill,93+nn,103+nn,9,
fill,104+nn,107+nn,2,
real,1
type,1
mat,2
secnum,53
*do,a,1+nn,12+nn,1
E,a,a+1
*enddo
type,1
mat,2
secnum,52
*do,b,13+nn,19+nn,1
E,b,b+1
*enddo
type,1
mat,2
secnum,51
*do,b,20+nn,21+nn,1
E,b,b+1
*enddo
type,1
mat,2
secnum,46
*do,b,22+nn,51+nn,1
E,b,b+1
*enddo
type,1
mat,2
secnum,42
*do,b,52+nn,61+nn,1
E,b,b+1
*enddo
type,1
mat,2
secnum,46
*do,b,63+nn,92+nn,1
E,b,b+1
*enddo
type,1
mat,2
secnum,42
*do,b,93+nn,102+nn,1
E,b,b+1
*enddo
type,1
mat,2
secnum,54
*do,b,104+nn,106+nn,1
E,b,b+1
*enddo
!************************************************
!
!link
!
!************************************************
/prep7
nn=300000
mm=400000
!*****************************ES
*do,j,nn+1,nn+2301*2+1,2301
*do,a,j,10*25+j,25
type,3
mat,1
real,3
E,a,mm+3
mm=mm+1
*enddo
mm=400000
*enddo
*do,j,nn+1+11*25,nn+2301*3,2301
*do,a,j,2*25+j,25
type,3
mat,1
real,4
E,a,mm+14
mm=mm+1
*enddo
mm=400000
*enddo
*do,j,nn+1+14*25,nn+2301*3,2301
*do,a,j,2*25+j,25
type,3
mat,1
real,5
E,a,mm+17
mm=mm+1
*enddo
mm=400000
*enddo
*do,j,nn+1+17*25,nn+2301*3,2301
type,3
mat,1
real,6
E,j,mm+20
mm=400000
*enddo
!****************************EM
mm=400000
*do,j,nn+1+21*25,nn+2301*3,2301
type,3
mat,1
real,6
E,j,mm+20
mm=400000
*enddo
*do,j,nn+1+22*25,nn+2301*3,2301
*do,a,j,2*25+j,25
type,3
mat,1
real,5
E,a,mm+19
mm=mm-1
*enddo
mm=400000
*enddo
*do,j,nn+1+25*25,nn+2301*3,2301
*do,a,j,2*25+j,25
type,3
mat,1
real,4
E,a,mm+16
mm=mm-1
*enddo
mm=400000
*enddo
*do,j,nn+1+28*25,nn+2301*3,2301
*do,a,j,10*25+j,25
type,3
mat,1
real,3
E,a,mm+13
mm=mm-1
*enddo
mm=400000
*enddo
!***************************************WM
mm=450000
*do,j,nn+1+40*25,nn+2301*3,2301
*do,a,j,10*25+j,25
type,3
mat,1
real,3
E,a,mm+3
mm=mm+1
*enddo
mm=450000
*enddo
*do,j,nn+1+51*25,nn+2301*3,2301
*do,a,j,2*25+j,25
type,3
mat,1
real,4
E,a,mm+14
mm=mm+1
*enddo
mm=450000
*enddo
*do,j,nn+1+54*25,nn+2301*3,2301
*do,a,j,2*25+j,25
type,3
mat,1
real,5
E,a,mm+17
mm=mm+1
*enddo
mm=450000
*enddo
*do,j,nn+1+57*25,nn+2301*3,2301
type,3
mat,1
real,6
E,j,mm+20
mm=450000
*enddo
!***************************************WS
mm=450000
*do,j,nn+1+61*25,nn+2301*3,2301
type,3
mat,1
real,6
E,j,mm+20
mm=450000
*enddo
*do,j,nn+1+62*25,nn+2301*3,2301
*do,a,j,2*25+j,25
type,3
mat,1
real,5
E,a,mm+19
mm=mm-1
*enddo
mm=450000
*enddo
*do,j,nn+1+65*25,nn+2301*3,2301
*do,a,j,2*25+j,25
type,3
mat,1
real,4
E,a,mm+16
mm=mm-1
*enddo
mm=450000
*enddo
*do,j,nn+1+68*25,nn+2301*3,2301
*do,a,j,10*25+j,25
type,3
mat,1
real,3
E,a,mm+13
mm=mm-1
*enddo
mm=450000
*enddo
nummrg,node
numcmp,node
FINISH