ansys建模实例-混凝土箱梁(命令流).doc
-/TITLE, BRIDGE XX /PREP7 ! 定义单元类型ET,1,SHELL63ET,2,SOLID45ET,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-5MP,EX,2,3.5E10 MP,DENS,2,2.6E3 MP,PRXY,2,0.1667MP,EX,3,3.5E15 MP,DENS,3,7.8E3 MP,PRXY,3,0.3 *set,nz,154!定义箱梁N1R,1,7,7,7,7!定义箱梁N2R,2,8,8,8,8!定义箱梁N3R,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, 1n, 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,1n, 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+6n, 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+6n, 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+6n, 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+6n, 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+6n, 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， 钢， 厚度8type,1 mat,1 real,2 *do, j, 1, 10, 1e, i*ni+j, i*ni+j+1, (i+1)*ni+j+1, (i+1)*ni+j, *enddo !顶板，shell63， 钢， 厚度7type,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, 1e, 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+23e, 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+1e, insty1+i-36+127-37+1, insty1+i-ni+127-37+2e, 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+1e, instb1+12, insty1+1+2*(127-37)+2e, 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, 5e, instb1+i+j*4, instb1+i+j*4+1*enddo *enddo finish/SOL!进入求解器施加荷载ESEL,S,TYPE,3 !选择现有单元类型为LINK8的单元EPLOT BFE,all,TEMP,1,tem1 !对目前显示的所有单元施加温度荷载，大小为-tem1ESEL,none !全不选择 /PREP7 !进入前处理器,生成Y2钢绞线的单元*set,insty2,insty1+3*(127-37+1)*do, i, 66,114,1n, insty2+(i-65), concrete(56,1),concrete(56,2),concrete(56,3)+z(i) n, insty2+(i-65)
