/TITLE, BRIDGE XX/PREP7! 定义单元类型ET,1,SHELL63ET,2,SOLID45ET,3,LINK8!定义1,钢，2，混凝土，3，无限刚度MP,EX,1,MP,DENS,1,MP,PRXY,1,mp,alpx,1,MP,EX,2,MP,DENS,2,MP,PRXY,2,MP,EX,3,MP,DENS,3,MP,PRXY,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=*tem2=*!定义纵向梯度参数表*dim, z, array, nz*SET,z(1),130,260,371,482,593,702,811,920,1029,1138,1247,1356,146 5,1574,1683,,,*SET,z(19),,,,,,,,,,,,,*SET,z(31),,,,,,*SET,z(37),,,,,*SET,z(41),,,,,,,,,,,,,,*SET,z(55),,,,,,,,,,,*SET,z(65),,,,,,,,*SET,z(73),,,,,,,,*SET,z(80),,,,,,,,,,,*SET,z(91),,,,,,,,,,*SET,z(100),,,,,,,,,*SET,z(109),,,,,,,*SET,z(115),,,,,,,,,,,,*SET,z(127),,,,,,,,,*SET,z(135),,,,,,,,,,*SET,z(145),,,,,,,,,,!!定义横隔梁的梯度参数表*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,,,-133,0,,,400,,633,-633,-533,,,-133, *SET,initial(17,1),0,,,400,,633,-633,-633,633,633,-766,,-633,,-50 0,500,,633,*SET,initial(35,1),,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,36 8,368,368,*SET,initial(1,3), , , , , , , , , , , ,*SET,initial(12,3), , , , , , , , , , , , ,*SET,initial(24,3), , , , , , , , , , , , ,!定义横隔梁参数表*set,nv,25*dim,vertical,array,nv,3*SET,vertical(1 ,3), , , , , , , , , , , , ,*SET,vertical(14,3), , , , , , , , , , , ,*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, , ,-133 ,0 , , , , , , , ,-133*SET,vertical(14,1),0, , , , , , , , , , ,!!定义混凝土的节点参数表*set,nc,135*dim,concrete,array,nc,3*SET,concrete(1,1), , , , , , , , , ,*SET,concrete(11,1), , , , , , , , , ,*SET,concrete(21,1), , , , , , , , , ,*SET,concrete(31,1), , , , , , , , , , *SET,concrete(41,1), , , , , , , , , , *SET,concrete(51,1), , , , , , , , , , *SET,concrete(61,1), , , , , , , , , ,*SET,concrete(71,1), , , , , , , , , , *SET,concrete(81,1), , , , , , , , , ,*SET,concrete(91,1), , , , , , , , , ,*SET,concrete(101,1), , , , , , , , , , *SET,concrete(111,1), , , , , , , , , ,*SET,concrete(121,1), , , , , , , , , , *SET,concrete(131,1), , , , ,*SET,concrete(1, 2), , , , , , , , , ,*SET,concrete(11, 2), , , , , , , , , , *SET,concrete(21, 2), , , , , , , , , , *SET,concrete(31, 2), , , , , , , , , , *SET,concrete(41, 2), , , , , , , , , , *SET,concrete(51, 2), , , , , , , , , , *SET,concrete(61, 2), , , , , , , , , , *SET,concrete(71, 2), , , , , , , , , , *SET,concrete(81, 2), , , , , , , , , , *SET,concrete(91, 2), , , ,,, ,,,368,*SET,concrete(101, 2), ,368,,,368, ,,,,*SET,concrete(111, 2), ,,,,, ,,,,*SET,concrete(121, 2), ,,,,,,,,,*SET,concrete(131, 2), ,,,,,*SET,concrete(1, 3), , , , , , , , , ,*SET,concrete(11, 3), , , , , , , , , ,*SET,concrete(21, 3), , , , , , , , , ,*SET,concrete(31, 3), , , , , , , , , ,*SET,concrete(41, 3), , , , , , , , , ,*SET,concrete(51, 3), , , , , , , , , ,*SET,concrete(61, 3), , , , , , , , , ,*SET,concrete(71, 3), , , , , , , , , ,*SET,concrete(81, 3), , , , , , , , , ,*SET,concrete(101, 3), ,,,,, ,,,,*SET,concrete(111, 3), ,,,,, ,,,,*SET,concrete(121, 3), ,,,,,,,,,*SET,concrete(131, 3), ,,,,,!节点编号，生成箱梁的节点*do, i, 1,1,1*do, j, 1, nin, j+ni*i, initial(j,1),initial(j,2),initial(j,3)+z(i) *enddo*enddo*do, j, 1, nin, 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, nin, 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, 32n, infl+i*nfl+j, vertical(j-7,1), vertical(j-7,2), vertical(j-7,3)+z(k-1)*enddon, 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, ncn, incr+j+nc*i, concrete(j,1),concrete(j,2),concrete(j,3)+z(i)*enddo*enddo*do, i, 2, nz*do, j, 1, ncn, 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, 1n, 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+in, instb1+i, , +12*i, +z(k)n, instb2+i, , +12*i, +z(k)*enddo*do, i, 1, 4*SET,k,27+in, instb1+i+4, 0 , +12*i, z(k) n, instb2+i+4, 0 , +12*i, z(k) *enddo*do, i, 1, 4*SET,k,27+in, instb1+i+8, , +12*i, +z(k) n, instb2+i+8, , +12*i, +z(k)*enddo*do, i, 1, 4*SET,k,128+i+6n, instb1+i+12, , *(i+6), +z(k) n, instb2+i+12, , *(i+6), +z(k) *enddo*do, i, 1, 4*SET,k,128+i+6n, instb1+i+16, 0 , *(i+6), z(k) n, instb2+i+16, 0 , *(i+6), z(k) *enddo*do, i, 1, 4*SET,k,128+i+6n, instb1+i+20, , *(i+6), +z(k) n, instb2+i+20, , *(i+6), +z(k) *enddo!Y2*do, i, 1, 4*SET,k,56+in, instb1+i+24, , +12*i, +z(k)n, instb2+i+24, , +12*i, +z(k)*enddo*do, i, 1, 4*SET,k,56+in, instb1+i+28, ,+12*i, +z(k)n, instb2+i+28, ,+12*i, +z(k)*enddo*do, i, 1, 4*SET,k,115+i+6n, instb1+i+32, ,*(i+6), +z(k) n, instb2+i+32, ,*(i+6), +z(k) *enddo*do, i, 1, 4*SET,k,115+i+6n, instb1+i+ni, ,*(i+6), +z(k)n, instb2+i+ni, ,*(i+6), +z(k)*enddo!生成箱梁的单元*do, i, 1,nz-1,1!底板，shell63，钢，厚度8type,1mat,1real,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,1mat,1real,1*do, j, 27, 30, 1e, i*ni+j, i*ni+j+1, (i+1)*ni+j+1, (i+1)*ni+j, *enddo*do, j, 32, 35, 1e, i*ni+j, i*ni+j+1, (i+1)*ni+j+1, (i+1)*ni+j, *enddo!腹板、底板加劲肋，shell63，钢，厚度4type,1mat,1real,3*do, j, 1, 11, 1e, i*ni+j, (i+1)*ni+j, (i+1)*ni+j+11, i*ni+j+11*enddoe, i*ni+12, (i+1)*ni+12, (i+1)*ni+23, i*ni+23e, i*ni+23, (i+1)*ni+23, (i+1)*ni+24, i*ni+24e, i*ni+24, (i+1)*ni+24, (i+1)*ni+29, i*ni+29e, i*ni+22, (i+1)*ni+22, (i+1)*ni+25, i*ni+25e, i*ni+25, (i+1)*ni+25, (i+1)*ni+26, i*ni+26e, i*ni+26, (i+1)*ni+26, (i+1)*ni+34, i*ni+34*enddo!生成横隔板单元type,1mat,1real,3*do, i, 1, 18, 1!生成横隔板单元*do, l, 0, 9, 1e, 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,1e, 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, 1e, 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, 1e, infl+i*nfl+l+19, infl+i*nfl+l+20, infl+i*nfl+l+27, infl+i*nfl+l+26e, infl+i*nfl+l+33, infl+i*nfl+l+34, infl+i*nfl+l+20, infl+i*nfl+l+19*enddoe, 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,2mat,2real, 4*do, i, 1, 153, 1*do, j, 0, 130, 3*do, k, 2, 3, 1e,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,3mat,1real, 5*do, i, 37, 126, 1e, 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 *enddoe, 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*enddofinish/SOL!进入求解器施加荷载ESEL,S,TYPE,,3 !选择现有单元类型为LINK8的单元EPLOTBFE,all,TEMP,1,tem1 !对目前显示的所有单元施加温度荷载，大小为-tem1 ESEL,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)+141-66+1, concrete(80,1),concrete(80,2),concrete(80,3)+z(i)*enddo*do, i, 66, 113, 1e, insty2+i-65,insty2+i-65+1e, insty2+(i-65)+141-66+1, insty2+(i-65)+141-66+2*enddoe, instb1+28, insty2+1e, instb1+32, insty2+1+141-66+1e, instb1+33, insty2+(113-65)e, instb1+37,insty2+(113-65)+141-66+1*do, i, 1, 3*do, j ,0, 3e, instb1+i+j*4+24, instb1+i+j*4+1+24*enddo*enddoFINISH/SOL!进入求解器施加荷载BFE,all,TEMP,1,tem2 !对目前选择的单元即Y2施加温度荷载，大小为-tem2 ESEL,all !选择所有的单元FINISH/PREP7!生成钢筋锚固区type,2mat,2real, 4*do, i, 1, 3e, 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)+69e, 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, 3e, 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,8d, ni*i+j, all,*enddo*enddo*do, i, 94, 96 *do, j, 4,8d, ni*i+j, all, *enddo*enddo*do, i, 153, 154 *do, j, 4,8d, ni*i+j, ux, *enddo*enddo*do, i, 153, 154 *do, j, 4,8d, ni*i+j, uz, *enddo*enddoFINISH/PREP7cpintf, ux, cpintf, uy, cpintf, uz, nsel, alleplot/solu solve。