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+28)+84,

incr+nc*(i+27)+69,incr+nc*(i+27)+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+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+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

instb2+29+i,

instb2+37+i,

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