箱梁桥ansys命令流

/TITLE, BRIDGE XX
/PREP7
! 定义单元类型
ET,1,SHELL63
ET,2,SOLID45
ET,3,LINK8
!定义1,钢，2，混凝土，3，无限刚度
MP,EX,1,2.1E11
MP,DENS,1,7.8E3
MP,PRXY,1,0.3
mp,alpx,1,1.0e-5
MP,EX,2,3.5E10
MP,DENS,2,2.6E3
MP,PRXY,2,0.1667
MP,EX,3,3.5E15
MP,DENS,3,7.8E3
MP,PRXY,3,0.3
*set,nz,154
!定义箱梁N1
R,1,7,7,7,7

!定义箱梁N2
R,2,8,8,8,8

!定义箱梁N3
R,3,4,4,4,4

!定义混凝土solid65的实常数
R, 4

!定义钢绞线的实常数
R, 5, 140

!降温=预加力/线膨胀系数*钢筋截面积*弹性模量
tem1=-1190.4e6/(1.0e-5*2.11e11)
tem2=-1339.2e6/(1.0e-5*2.11e11)
!定义纵向梯度参数表
*dim, z, array, nz
*SET,z(1),130,260,371,482,593,702,811,920,1029,1138,1247,1356,1465,1574,1683,1787.19,1891.38,1995.57
*SET,z(19),2099.76,2203.95,2308.14,2412.33,2516.52,2620.71,2724.9,2829.09,2933.28,3037.47,3141.66,3245.85,
*SET,z(31),3350.04,3454.23,3558.42,3662.61,3766.8,3870.99
*SET,z(37),3975.18,4079.37,4183.56,4287.75,
*SET,z(41),4391.94,4496.13,4600.32,4704.51,4808.7,4912.89,5017.08,5121.27,5225.46,5329.65,5433.84,5538.03,5642.22,5746.41
*SET,z(55),5850.6,5954.79,6058.98,6163.17,6267.36,6371.55,6475.74,6579.93,6682.47,6785.01,
*SET,z(65),6887.55,6990.09,7092.63,7195.17,7297.71,7400.25,7502.79,7605.33
*SET,z(73),7707.87,7810.41,7912.95,8024.28,8135.61,8246.94,8358.27,
*SET,z(80),8469.6,8580.93,8680.93,8780.93,8880.93,8980.93,9080.93,9180.93,9280.93,9380.93,9480.93
*SET,z(91),9580.93,9691.93,9802.93,9913.93,10024.93,10135.93,10246.93,10346.93,10446.93,
*SET,z(100),10546.93,10646.93,10746.93,10846.93,10946.93,11046.93,11146.93,11246.93,11358.26
*SET,z(109),11469.59,11580.92,11692.25,11803.58,11914.91,12005.91,
*SET,z(115),12096.91,12187.91,12278.91,12369.91,12460.91,12551.91,12642.91,12733.91,12824.91,12915.91,13006.91,13125.91
*SET,z(127),13244.91,13363.91,13482.91,13601.91,13720.91,13839.91,13939.91,14039.91,
*SET,z(135),14139.91,14239.91,14339.91,14439.91,14539.91,14639.91,14739.91,14839.91,14942.41,15044.91
*SET,z(145),15147.41,15249.91,15352.41,15454.91,15557.41,15659.91,15770.91,15881.91,15992.91,16109.91
!
!定义横隔梁的梯度参数表
*set,nd,18
*dim, del, array, nd
*SET,del(1), 2,5,15,31,47,63,76,82,92,95,98,108,114,125,132,142,
*SET,del(17),150,153
!
!定义初始面座标参数表
*set,ni,36
*dim, initial, array,ni, 3
*SET,initial(1,1),-633,-533,-399.6667,-266.333,-133,0,133.333,266.6666667,400,533.333,633,-633,-533,-399.667,-266.3333333,-133,
*SET,initial(17,1),0,133.3333333,266.6666667,400,533.3333333,633,-633,-633,633,633,-766,-699.5,-633,-566.5,-500,500,566.5,633,
*SET,initial(35,1),699.5,766
*SET,initial(1,2),0,0,0,0,0,0,0,0,0,0,0,0
*SET,initial(12,2),60,60,60,60,60,60,60,60,60,60,60,
*SET,Initial(23,2),184,276,184,276,368,368,368,368,368,368,368,368,368,368,
*SET,initial(1,3),115.82 ,97.53 ,73.13 ,48.73 ,24.34 ,0.00 ,-24.40 ,-48.79 ,-73.19 ,-97.59 ,-115.82 ,
*SET,initial(

12,3),115.82 ,97.53 ,73.13 ,48.73 ,24.34 ,0.00 ,-24.40 ,-48.79 ,-73.19 ,-97.59 ,-115.82 ,115.82 ,
*SET,initial(24,3),115.82 ,-115.82 ,-115.82 ,140.16 ,127.99 ,115.82 ,103.66 ,91.49 ,-91.49 ,-103.66 ,-115.82 ,-127.99 ,-140.16

!定义横隔梁参数表
*set,nv,25
*dim,vertical,array,nv,3
*SET,vertical(1 ,3),97.53 ,73.13 ,48.73 ,24.34 ,0.00 ,-24.40 ,-48.79 ,-73.19 ,-97.59 ,97.53 ,73.13 ,48.73 ,24.34
*SET,vertical(14,3),0.00 ,-24.40 ,-48.79 ,-73.19 ,-97.59 ,73.19 ,48.73 ,24.34 ,0.00 ,-24.40 ,-48.79 ,-73.19
*SET,vertical(1,2),184,184,184,184,184,184,184,184,184,
*SET,vertical(10,2),276,276,276,276,276,276,276,276,276,
*SET,vertical(19,2),368,368,368,368,368,368,368
*SET,vertical(1,1),-533,-399.667 ,-266.333 ,-133 ,0 ,133.333 ,266.667 ,400.000 ,533.333 ,-533.000 ,-399.667 ,-266.333 ,-133
*SET,vertical(14,1),0,133.333 ,266.667 ,400.000 ,533.333 ,-400.000 ,-266.333 ,-133.000 ,0.000 ,133.333 ,266.667 ,400.000
!
!定义混凝土的节点参数表
*set,nc,135
*dim,concrete,array,nc,3
*SET,concrete(1,1),-1301.625 ,-1301.625 ,-1301.625 ,-1248.225 ,-1248.225 ,-1248.225 ,-1194.825 ,-1194.825 ,-1194.825 ,-1141.425
*SET,concrete(11,1),-1141.425 ,-1141.425 ,-1088.025 ,-1088.025 ,-1088.025 ,-1034.625 ,-1034.625 ,-1034.625 ,-981.225 ,-981.225
*SET,concrete(21,1),-981.225 ,-927.825 ,-927.825 ,-927.825 ,-874.425 ,-874.425 ,-874.425 ,-821.025 ,-821.025 ,-821.025
*SET,concrete(31,1),-766.000 ,-766.000 ,-766.000 ,-699.500 ,-699.500 ,-699.500 ,-633.000 ,-633.000 ,-633.000 ,-566.500
*SET,concrete(41,1),-566.500 ,-566.500 ,-500.000 ,-500.000 ,-500.000 ,-450.000 ,-450.000 ,-450.000 ,-400.000 ,-400.000
*SET,concrete(51,1),-400.000 ,-333.167 ,-350.000 ,-333.167 ,-266.333 ,-250.000 ,-266.333 ,-199.667 ,-199.667 ,-199.667
*SET,concrete(61,1),-133.000 ,-133.000 ,-133.000 ,-66.500 ,-66.500 ,-66.500 ,0.000 ,0.000 ,0.000 ,66.500
*SET,concrete(71,1),66.500 ,66.500 ,133.000 ,133.000 ,133.000 ,199.667 ,199.667 ,199.667 ,266.333 ,250.000
*SET,concrete(81,1),266.333 ,333.167 ,350.000 ,333.167 ,400.000 ,400.000 ,400.000 ,450.000 ,450.000 ,450.000
*SET,concrete(91,1),500.000 ,500.000 ,500.000 ,566.500 ,566.500 ,566.500 ,633.000 ,633.000 ,633.000 ,699.500
*SET,concrete(101,1),699.500 ,699.500 ,766.000 ,766.000 ,766.000 ,821.025 ,821.025 ,821.025 ,874.425 ,874.425
*SET,concrete(111,1),874.425 ,927.825 ,927.825 ,927.825 ,981.225 ,981.225 ,981.225 ,1034.625 ,1034.625 ,1034.625
*SET,concrete(121,1),1088.025 ,1088.025 ,1088.025 ,1141.425 ,1141.425 ,1141.425 ,1194.825 ,1194.825 ,1194.825 ,1248.225
*SET,concrete(131,1),1248.225 ,1248.225 ,1301.625 ,1301.625 ,1301.625
*SET,concrete(1, 2), 468.25,443.219 ,418.188 ,468.250 ,440.716 ,413.181 ,468.250 ,443.219 ,408.175 ,468.250
*SET,concrete(11, 2), 443.219 ,403.169 ,468.250 ,430.703 ,398.163 ,468.250 ,430.703 ,393.156 ,468.250 ,418.188
*SET,concrete(21, 2), 388.150 ,468.250 ,418.188 ,383.144 ,468.250 ,418.188 ,378.138 ,468.250 ,418.188 ,373.131
*SET,concrete(31, 2), 468.250 ,418

.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250
*SET,concrete(41, 2), 418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188
*SET,concrete(51, 2), 368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000
*SET,concrete(61, 2), 468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250
*SET,concrete(71, 2), 418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188
*SET,concrete(81, 2), 368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000 ,468.250 ,418.188 ,368.000
*SET,concrete(91, 2), 468.25,418.188 ,368.000 ,468.25,418.188,368.000 ,468.25,418.188,368,468.25
*SET,concrete(101, 2), 418.188,368,468.25,418.188,368,468.250 ,418.188,373.131,468.25,418.188
*SET,concrete(111, 2), 378.138,468.25,418.188,383.144,468.25,418.188 ,388.15,468.25,430.703,393.156
*SET,concrete(121, 2), 468.25,430.703,398.163,468.25,443.219,403.169,468.25,443.219,408.175,468.25
*SET,concrete(131, 2), 440.716,413.181,468.25,443.219,418.188,
*SET,concrete(1, 3), 238.1649941,238.165 ,238.165 ,228.394 ,228.394 ,228.394 ,218.623 ,218.623 ,218.623 ,208.852
*SET,concrete(11, 3), 208.852 ,208.852 ,199.082 ,199.082 ,199.082 ,189.311 ,189.311 ,189.311 ,179.540 ,179.540
*SET,concrete(21, 3), 179.540 ,169.769 ,169.769 ,169.769 ,159.998 ,159.998 ,159.998 ,150.227 ,150.227 ,150.227
*SET,concrete(31, 3), 140.159 ,140.159 ,140.159 ,127.991 ,127.991 ,127.991 ,115.823 ,115.823 ,115.823 ,103.655
*SET,concrete(41, 3), 103.655 ,103.655 ,91.488 ,91.488 ,91.488 ,82.339 ,82.339 ,82.339 ,73.190 ,73.190
*SET,concrete(51, 3), 73.190 ,60.961 ,64.041 ,60.961 ,48.732 ,45.744 ,48.732 ,36.534 ,36.534 ,36.534
*SET,concrete(61, 3), 24.336 ,24.336 ,24.336 ,12.168 ,12.168 ,12.168 ,0.000 ,0.000 ,0.000 ,-12.168
*SET,concrete(71, 3), -12.168 ,-12.168 ,-24.336 ,-24.336 ,-24.336 ,-36.534 ,-36.534 ,-36.534 ,-48.732 ,-45.744
*SET,concrete(81, 3), -45.744 ,-60.961 ,-64.041 ,-60.961 ,-73.190 ,-73.190 ,-73.190 ,-82.339 ,-82.339 ,-82.339
*SET,concrete(91, 3), -91.48756137,-91.488 ,-91.488 ,-103.65,-103.65,-103.655 ,-115.8232527,-115.8232527,-115.82,-127.99
*SET,concrete(101, 3), -127.991,-127.9910,-140.159,-140.159,-140.159,-150.227 ,-150.23,-150.23,-159.9980217,-159.9980217
*SET,concrete(111, 3), -159.998,-169.768,-169.768,-169.768,-179.539,-179.540 ,-179.5397,-189.3106,-189.3106364,-189.3106364
*SET,concrete(121, 3), -199.0815,-199.0815,-199.0815,-208.852,-208.8523,-208.8523,-218.62,-218.623,-218.623251,-228.3941226
*SET,concrete(131, 3), -228.3941226,-228.3941226,-238.1649941,-238.1649941,-238.1649941,

!节点编号，生成箱梁的节点
*do, i, 1,1,1
*do, j, 1, ni
n, j+ni*i, initial(j,1),initial(j,2),initial(j,3)+z(i)
*enddo
*enddo
*do, i, nz,nz,1
*do, j, 1, ni
n, j+ni*i, initial(j,1),initial(j,2),initial(j,3)+z(i)
*enddo
*enddo
*do, i, 2,nz-1,1
*do, j, 1, ni

n, j+ni*i, initial(j,1),initial(j,2),initial(j,3)+z(i)
*enddo
*enddo
!生成横隔板节点

*set, infl, ni*(nz+1)
*set, nfl, 39
*do, i, 1, 18, 1
*SET,k,del(i)
*do, j, 1, 25, 1
n, infl+i*nfl+j, vertical(j,1), vertical(j,2), vertical(j,3)+z(k)
*enddo
*do, j, 26, 32
n, infl+i*nfl+j, vertical(j-7,1), vertical(j-7,2), vertical(j-7,3)+z(k-1)
*enddo
*do, j, 33, nfl
n, infl+i*nfl+j, vertical(j-14,1), vertical(j-14,2), vertical(j-14,3)+z(k+1)
*enddo
*enddo

!节点编号，生成混凝土的节点
*set,incr,infl+18*nfl+1
*do, i, 1, 1,
*do, j, 1, nc
n, incr+j+nc*i, concrete(j,1),concrete(j,2),concrete(j,3)+z(i)
*enddo
*enddo
*do, i, 2, nz
*do, j, 1, nc
n, incr+j+nc*i, concrete(j,1),concrete(j,2),concrete(j,3)+z(i)
*enddo
*enddo

!生成钢绞线的节点
*set,insty1,incr+nc*(nz+1)
!Y1
*do, i, 37, 127, 1
n, insty1+(i-36), concrete(53,1),concrete(53,2),concrete(53,3)+z(i)
n, insty1+(i-36)+127-37+1, concrete(68,1),concrete(68,2),concrete(68,3)+z(i)
n, insty1+(i-36)+2*(127-37)+2, concrete(83,1),concrete(83,2),concrete(83,3)+z(i)
*enddo
!Y2
*set,insty2,insty1+3*(127-37+1)
*do, i, 66,114,1
n, insty2+(i-65), concrete(56,1),concrete(56,2),concrete(56,3)+z(i)
n, insty2+(i-65)+141-66+1, concrete(80,1),concrete(80,2),concrete(80,3)+z(i)
*enddo


!生成钢绞线锚固区节点
*set,instb1,insty2+2*(141-66)
*set,instb2,instb1+40
!Y1
*do, i, 1, 4
*SET,k,27+i
n, instb1+i, -350.000 , 298.188+12*i, 64.0413+z(k)
n, instb2+i, -350.000 , 298.188+12*i, 64.0413+z(k)
*enddo
*do, i, 1, 4
*SET,k,27+i
n, instb1+i+4, 0 , 298.188+12*i, z(k)
n, instb2+i+4, 0 , 298.188+12*i, z(k)
*enddo
*do, i, 1, 4
*SET,k,27+i
n, instb1+i+8, 350.000 , 298.188+12*i, -64.0413+z(k)
n, instb2+i+8, 350.000 , 298.188+12*i, -64.0413+z(k)
*enddo
*do, i, 1, 4
*SET,k,128+i+6
n, instb1+i+12, -350.000 , 418.188-12*(i+6), 64.0413+z(k)
n, instb2+i+12, -350.000 , 418.188-12*(i+6), 64.0413+z(k)
*enddo
*do, i, 1, 4
*SET,k,128+i+6
n, instb1+i+16, 0 , 418.188-12*(i+6), z(k)
n, instb2+i+16, 0 , 418.188-12*(i+6), z(k)
*enddo
*do, i, 1, 4
*SET,k,128+i+6
n, instb1+i+20, 350.000 , 418.188-12*(i+6), -64.0413+z(k)
n, instb2+i+20, 350.000 , 418.188-12*(i+6), -64.0413+z(k)
*enddo
!Y2
*do, i, 1, 4
*SET,k,56+i
n, instb1+i+24, -250.000 , 298.188+12*i, 45.743+z(k)
n, instb2+i+24, -250.000 , 298.188+12*i, 45.743+z(k)
*enddo
*do, i, 1, 4
*SET,k,56+i
n, instb1+i+28, 250.000 ,298.188+12*i, -45.743+z(k)
n, instb2+i+28, 250.000 ,298.188+12*i, -45.743+z(k)
*enddo
*do, i, 1, 4
*SET,k,115+i+6
n, instb1+i+32, -250.000 ,418.188-12*(i+6), 45.743+z(k)
n, instb2+i+32, -250.000 ,418.188-12*(i+6), 45.743+z(k)
*enddo
*do, i, 1, 4
*SET,k,115+i+6
n, instb1+i+ni, 250.000 ,418.188-12*(i+6), -45.743+z(k)
n, instb2+i+ni, 250.000 ,418.188-12*(i+6), -45.743+z(k)
*enddo

!生

成箱梁的单元
*do, i, 1,nz-1,1
!底板，shell63， 钢， 厚度8
type,1
mat,1
real,2
*do, j, 1, 10, 1
e, i*ni+j, i*ni+j+1, (i+1)*ni+j+1, (i+1)*ni+j,
*enddo
!顶板，shell63， 钢， 厚度7
type,1
mat,1
real,1
*do, j, 27, 30, 1
e, i*ni+j, i*ni+j+1, (i+1)*ni+j+1, (i+1)*ni+j,
*enddo
*do, j, 32, 35, 1
e, i*ni+j, i*ni+j+1, (i+1)*ni+j+1, (i+1)*ni+j,
*enddo
!腹板、底板加劲肋，shell63， 钢， 厚度4
type,1
mat,1
real,3
*do, j, 1, 11, 1
e, i*ni+j, (i+1)*ni+j, (i+1)*ni+j+11, i*ni+j+11
*enddo
e, i*ni+12, (i+1)*ni+12, (i+1)*ni+23, i*ni+23
e, i*ni+23, (i+1)*ni+23, (i+1)*ni+24, i*ni+24
e, i*ni+24, (i+1)*ni+24, (i+1)*ni+29, i*ni+29
e, i*ni+22, (i+1)*ni+22, (i+1)*ni+25, i*ni+25
e, i*ni+25, (i+1)*ni+25, (i+1)*ni+26, i*ni+26
e, i*ni+26, (i+1)*ni+26, (i+1)*ni+34, i*ni+34
*enddo
!生成横隔板单元
type,1
mat,1
real,3

*do, i, 1, 18, 1
!生成横隔板单元
*do, l, 0, 9, 1
e, del(i)*ni+1+l, del(i)*ni+2+l, del(i)*ni+13+l, del(i)*ni+12+l,
*enddo
*do, l, 0, 7,1
e, infl+i*nfl+l+1, infl+i*nfl+l+2, infl+i*nfl+l+11, infl+i*nfl+l+10,
e, del(i)*ni+l+13,del(i)*ni+l+14, infl+i*nfl+l+2, infl+i*nfl+l+1,
*enddo
*do, l, 0, 5, 1
e, infl+i*nfl+l+11, infl+i*nfl+l+12, infl+i*nfl+l+20, infl+i*nfl+l+19,
*enddo
*do, l, 0, 5, 1
e, infl+i*nfl+l+19, infl+i*nfl+l+20, infl+i*nfl+l+27, infl+i*nfl+l+26
e, infl+i*nfl+l+33, infl+i*nfl+l+34, infl+i*nfl+l+20, infl+i*nfl+l+19
*enddo
e, del(i)*ni+31, infl+i*nfl+19, infl+i*nfl+26, (del(i)-1)*ni+31,
e, (del(i)+1)*ni+31, infl+i*nfl+33, infl+i*nfl+19, del(i)*ni+31,
e, infl+i*nfl+25, del(i)*ni+32, (del(i)-1)*ni+32, infl+i*nfl+32,
e, infl+i*nfl+nfl, (del(i)+1)*ni+32, del(i)*ni+32, infl+i*nfl+25,
e, del(i)*ni+12, del(i)*ni+13, infl+i*nfl+1,del(i)*ni+23
e, del(i)*ni+23, infl+i*nfl+1, infl+i*nfl+10, del(i)*ni+24,
e, del(i)*ni+24, infl+i*nfl+10, del(i)*ni+31, del(i)*ni+29,
e, infl+i*nfl+10, infl+i*nfl+11, infl+i*nfl+19, del(i)*ni+31,
e, del(i)*ni+21, del(i)*ni+22, del(i)*ni+25, infl+i*nfl+9,
e, infl+i*nfl+9, del(i)*ni+25, del(i)*ni+26, infl+i*nfl+18,
e, infl+i*nfl+18, del(i)*ni+26, del(i)*ni+34, del(i)*ni+32,
e, infl+i*nfl+17, infl+i*nfl+18, del(i)*ni+32, infl+i*nfl+25,
*enddo


!生成混凝土单元
type,2
mat,2
real, 4
*do, i, 1, 153, 1
*do, j, 0, 130, 3
*do, k, 2, 3, 1
e,incr+nc*i+j+k,incr+nc*(i+1)+j+k,incr+nc*(i+1)+j+k+3,incr+nc*i+j+k+3,incr+nc*i+j+k-1,incr+nc*(i+1)+j+k-1,incr+nc*(i+1)+j+k+2, incr+nc*i+j+k+2,
*enddo
*enddo
*enddo




!生成Y1钢绞线单元
type,3
mat,1
real, 5
*do, i, 37, 126, 1
e, insty1+i-36,insty1+i-ni+1
e, insty1+i-36+127-37+1, insty1+i-ni+127-37+2
e, insty1+(i-36)+2*(127-37)+2,insty1+(i-ni)+2*(127-37)+3
*enddo
e, instb1+4, insty1+1,
e, instb1+8, insty1+1+127-37+1
e, instb1+12, insty1+1+2*(127-37)+2
e, instb1+13, insty1+(127-36)
e, instb1+17, insty1+

(127-36)+127-37+1,
e, instb1+21, insty1+(127-36)+2*(127-37)+2
*do, i, 1, 3
*do, j ,0, 5
e, instb1+i+j*4, instb1+i+j*4+1
*enddo
*enddo
finish

/SOL　
!进入求解器施加荷载
ESEL,S,TYPE,,3 !选择现有单元类型为LINK8的单元
EPLOT
BFE,all,TEMP,1,tem1 !对目前显示的所有单元施加温度荷载，大小为-tem1
ESEL,none !全不选择

/PREP7
!进入前处理器,生成Y2钢绞线的单元

*set,insty2,insty1+3*(127-37+1)
*do, i, 66,114,1
n, insty2+(i-65), concrete(56,1),concrete(56,2),concrete(56,3)+z(i)
n, insty2+(i-65)+141-66+1, concrete(80,1),concrete(80,2),concrete(80,3)+z(i)
*enddo

*do, i, 66, 113, 1
e, insty2+i-65,insty2+i-65+1
e, insty2+(i-65)+141-66+1, insty2+(i-65)+141-66+2
*enddo
e, instb1+28, insty2+1
e, instb1+32, insty2+1+141-66+1
e, instb1+33, insty2+(113-65)
e, instb1+37,insty2+(113-65)+141-66+1
*do, i, 1, 3
*do, j ,0, 3
e, instb1+i+j*4+24, instb1+i+j*4+1+24
*enddo
*enddo
FINISH

/SOL
!进入求解器施加荷载
BFE,all,TEMP,1,tem2 !对目前选择的单元即Y2施加温度荷载，大小为-tem2
ESEL,all !选择所有的单元
FINISH
/PREP7
!生成钢筋锚固区
type,2
mat,2
real, 4
*do, i, 1, 3
e, instb2+i, instb2+4+i, instb2+5+i, instb2+1+i, incr+54+nc*(i+27),incr+69+nc*(i+27),incr+69+nc*(i+28), incr+54+nc*(i+28),
e, instb2+4+i,instb2+8+i, instb2+9+i,instb2+5+i, incr+nc*(i+27)+69,incr+nc*(i+27)+84, incr+nc*(i+28)+84, incr+nc*(i+28)+69

e, instb2+12+i, instb2+16+i, instb2+17+i, instb2+13+i, incr+54+nc*(i+134),incr+69+nc*(i+134),incr+69+nc*(i+135),incr+54+nc*(i+135)
e, instb2+16+i, instb2+20+i, instb2+21+i, instb2+17+i, incr+69+nc*(i+134),incr+84+nc*(i+134),incr+84+nc*(i+135),incr+69+nc*(i+135)

*enddo


*do, i,1, 3
e, instb2+24+i, instb2+28+i, instb2+29+i, instb2+25+i,incr+57+nc*(i+56),incr+81+nc*(i+56), incr+81+nc*(i+57),incr+57+nc*(i+57)
e, instb2+32+i, instb2+36+i, instb2+37+i, instb2+33+i,incr+57+nc*(i+121),incr+81+nc*(i+121), incr+81+nc*(i+122),incr+57+nc*(i+122)
*enddo


!加载
FINISH
/SOL
*do, i, 1, 2
*do, j, 4,8
d, ni*i+j, all,
*enddo
*enddo
*do, i, 94, 96
*do, j, 4,8
d, ni*i+j, all,
*enddo
*enddo
*do, i, 153, 154
*do, j, 4,8
d, ni*i+j, ux,
*enddo
*enddo
*do, i, 153, 154
*do, j, 4,8
d, ni*i+j, uz,
*enddo
*enddo
FINISH
/PREP7
cpintf, ux, 0.01
cpintf, uy, 0.01
cpintf, uz, 0.01

nsel, all
eplot
/solu
solve


3421,3426,3431,3443,3445,3447,3449,3454,3456,
19139,19141,19145,19146,19147,19160,19162,19175,19177
19190,19192,19205,19206,19207,19211,19213