大侠能帮忙看看这个程序对么？打锚杆

276774950

;定义顶板锚杆二维数组函数，便于确定顶锚杆的施加位置，顶6根锚杆，端部锚杆斜打2根
def maoshuzu1
    array bw3(2,6) ew3(2,6)
    bw3(1,1)=135.105
    bw3(1,2)=136.055
    bw3(1,3)=137.005
    bw3(1,4)=137.995
    bw3(1,5)=138.945
    bw3(1,6)=139.895
    bw3(2,1)=13.4
    bw3(2,2)=13.4
    bw3(2,3)=13.4
    bw3(2,4)=13.4
    bw3(2,5)=13.4
    bw3(2,6)=13.4
    ew3(1,1)=134.7577
    ew3(1,2)=136.055
    ew3(1,3)=137.005
    ew3(1,4)=137.995
    ew3(1,5)=138.945
    ew3(1,6)=140.2421
    ew3(2,1)=15.3696
    ew3(2,2)=15.4
    ew3(2,3)=15.4
    ew3(2,4)=15.4
    ew3(2,5)=15.4
    ew3(2,6)=15.3696
end

;左帮锚杆函数数组,3根
def maoshuzu2
   array bw(2,3)  ew(2,3)
   bw(1,1)=135
   bw(1,2)=135
   bw(1,3)=135
   bw(2,1)=11.6
   bw(2,2)=12.4
   bw(2,3)=13.2
   ew(1,1)=133.6
   ew(1,2)=133.6
   ew(1,3)=133.6
   ew(2,1)=11.6
   ew(2,2)=12.4
   ew(2,3)=13.2
end
;右帮锚杆函数数组,3根
def maoshuzu3
   array bw2(2,3)  ew2(2,3)
   bw2(1,1)=140
   bw2(1,2)=140
   bw2(1,3)=140
   bw2(2,1)=11.6
   bw2(2,2)=12.4
   bw2(2,3)=13.2
   ew2(1,1)=141.4
   ew2(1,2)=141.4
   ew2(1,3)=141.4
   ew2(2,1)=11.6
   ew2(2,2)=12.4
   ew2(2,3)=13.2
end
;调用数组函数，并作索引确定循环坐标值
def dingbanzuobian
  maoshuzu1
  dy_bw=bw3(1,index)
  dz_bw=bw3(2,index)
  dy_ew=ew3(1,index)
  dz_ew=ew3(2,index)
end
def zuobangzuobiao
  maoshuzu2
  ly_bw=bw(1,index)
  lz_bw=bw(2,index)
  ly_ew=ew(1,index)
  lz_ew=ew(2,index)
end
def youbangzuobiao
  maoshuzu3
  ry_bw=bw2(1,index)
  rz_bw=bw2(2,index)
  ry_ew=ew2(1,index)
  rz_ew=ew2(2,index)
end
set large
;开挖巷道80m加锚杆
model null range group f6 x 20 100 y 135 140

;打顶板锚杆
def maodingban1
   loop n1(1,30)
          x_d=2*n1+20
      loop index(1,6)
        dingbanzuobiao
        command
          sel cable id x_d beg dy_bw x_d dz_bw end dy_ew x_d dz_ew
          sel cable prop emod 2e10 ytens 3e8 xcarea 2.5434e-4 gr_k 1.28e10 gr_coh 3e5 gr_per 0.088
        end_command
      end_loop
   end_loop
end
;打左帮锚杆
def maozuobang1
   loop n1(1,30)
          x_d=2*n1+20
      loop index(1,3)
        zuobangzuobiao
        command
          sel cable id x_d beg ly_bw x_d lz_bw end ly_ew x_d lz_ew
          sel cable prop emod 2e10 ytens 3e8 xcarea 2.0096e-4 gr_k 1.28e10 gr_coh 3e5 gr_per 0.088
        end_command
      end_loop
   end_loop
end

;打右帮锚杆
def maoyoubang1
   loop n1(1,30)
          x_d=2*n1+20
      loop index(1,3)
        youbangzuobiao
        command
          sel cable id x_d beg ry_bw x_d rz_bw end ry_ew x_d rz_ew
          sel cable prop emod 2e10 ytens 3e8 xcarea 2.0096e-4 gr_k 1.28e10 gr_coh 3e5 gr_per 0.088
        end_command
      end_loop
   end_loop
end

maodingban1
maozuobang1
maoyoubang1
step 2000
save 全支护80m.sav