弱问：我要在洞库里打锚杆，怎么实现

caikewei1984

我刚学，对锚杆这一块感觉很不好理解，我是不是要把锚杆的间距找出来，然后再打入啊，书中fish语言建的数组理解不了。
急求

caikewei1984

zenm mei ren hui da

lookcity

计算好锚杆群起点和终点的坐标存到数组中，在大锚杆的时候从数组中读取坐标数据给beg 和 end就可以了。需要用到循环语句。

caikewei1984

3# lookcity 谢谢了，非常感谢，我建了数组，但是锚杆大部分受压，怎么回事呢

li_2001

计算好锚杆群起点和终点的坐标存到数组中，在大锚杆的时候从数组中读取坐标数据给beg 和 end就可以了。需要用到循环语句。

hondu5

那你得注意调整参数，以及边界条件了。 4# caikewei1984

chaylo

3# lookcity  谢谢了，非常感谢，我建了数组，但是锚杆大部分受压，怎么回事呢
caikewei1984 发表于 2009-3-28 17:48

应该是隧道上部的锚杆受拉下部的受压吧

ca1b000

def ZCFJBCS
;圆弧半径
  R_O=17.5
;锚杆间距
  B_O=3
  B_YM=1
  B_YWM=1.5
;钢筋长度
  L_O=6
  L_L=L_O-1
  L_BM=8
  L_MB=L_BM-1
  L_JM=9
  L_XJM=L_JM-1
  L_WO=5
  L_WL=L_WO-1
  L_QM=7
  L_QBM=L_QM-1
  L_SWM=3.5
  
;定义根数
  AAAA=110
  BBBB=4
;端点1 坐标(拱顶)
  X1_1=2.4
  Y1_1=72.2
  Z1_1=76.9
;锚杆截面面积和外周长
  
  D_25=0.00031416
  KLQD_25=180000.0
  D_28=0.0004146912
  KLQD_28=230000.0
  D_32=0.0004901
  KLQD_32=250000.0
  D_P22=0.0003801
  KLQD_P22=335.0*1000000*D_P22
  D_P25=0.0004909
  KLQD_P25=335.0*1000000*D_P25
  D_P28=0.00061544
  KLQD_P28=335.0*1000000*D_P28
  D_P32=0.00080384
  KLQD_P32=335.0*1000000*D_P32
  
  S_25=0.0785
  S_28=0.0879648
  S_32=0.1005312
  S_P22=0.0691152
  S_P25=0.07854
  S_P28=0.0879648
  S_P32=0.1005312

;水泥浆的刚度
  SNJQD_O=20000000
;水泥浆的粘结力
  SNJNJL_O=1000000
;预应力值
  yuyinli_O=120000
;初始锚杆编号和节点编号
  SN_ID=0
  SD_ID=0
end
ZCFJBCS


def get_cid1
now_cid_1=s_head
now_cid_2=0
loop while now_cid_1 # null
if s_cid(now_cid_1) > now_cid_2
now_cid_2=s_cid(now_cid_1)
endif
now_cid_1=s_next(now_cid_1)
endloop
SD_ID=now_cid_2
end
get_cid1




;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;28(6M)中间没有（偶数个），弧形发散。;预应力中空注浆锚杆 28  @x@'=1.5m x 3.0m  T=120kN,L=6m
def caibo1
LOOP AB(1,AAAA)
A1_Z=Z1_1-B_O*(AB-1)
A2_Z=Z1_1-B_O*(AB-1)
A3_Z=Z1_1-B_O*(AB-1)
LOOP CD(1,BBBB)
;坐标
A1_X=X1_1-R_O*SIN((CD-0.5)*B_O/R_O)
A2_X=A1_X-L_L*SIN((CD-0.5)*B_O/R_O)
A3_X=A1_X-L_O*SIN((CD-0.5)*B_O/R_O)

A1_Y=Y1_1-(R_O-R_O*COS(((CD-0.5)*B_O)/R_O))
A2_Y=A1_Y+L_L*COS(((CD-0.5)*B_O)/R_O)
A3_Y=A1_Y+L_O*COS(((CD-0.5)*B_O)/R_O)
;编号
N_ID=(AB-1)*BBBB+CD+SN_ID
D_ID=(AB-1)*L_O*BBBB+(CD-1)*L_O+SD_ID
B_1=D_ID+2
B_2=D_ID+L_O-1
B_3=D_ID+L_O
B_4=D_ID+L_O
B_5=D_ID+1

COMMAND
sel cable id=N_ID beg A1_X A1_Y A1_Z end A2_X A2_Y A2_Z nseg L_L
sel cable id=N_ID beg A2_X A2_Y A2_Z end A3_X A3_Y A3_Z nseg 1
sel cable prop emod 0.2e12 yten KLQD_28 xcarea D_28 gr_per S_28  gr_k SNJQD_O gr_coh SNJNJL_O range cid B_1 B_2
sel cable prop emod 0.2e12 yten KLQD_28 xcarea D_28 gr_per S_28  gr_k SNJQD_O gr_coh SNJNJL_O range cid B_3 B_4
sel cable prop emod 0.2e12 yten KLQD_28 xcarea D_28 gr_per S_28  gr_k SNJQD_O gr_coh SNJNJL_O range cid B_5 B_5
sel cable id N_ID pretension yuyinli_O range cid B_5 B_3
END_COMMAND
   end_loop
   end_loop
SN_ID=SN_ID+AAAA*BBBB
SD_ID=SD_ID+AAAA*BBBB*L_O
end
CAIBO1

def caibo2
LOOP AB(1,AAAA)
A1_Z=Z1_1-B_O*(AB-1)
A2_Z=Z1_1-B_O*(AB-1)
A3_Z=Z1_1-B_O*(AB-1)
LOOP CD(1,BBBB)
;坐标
A1_X=X1_1+R_O*SIN((CD-0.5)*B_O/R_O)
A2_X=A1_X+L_L*SIN((CD-0.5)*B_O/R_O)
A3_X=A1_X+L_O*SIN((CD-0.5)*B_O/R_O)
A1_Y=Y1_1-(R_O-R_O*COS(((CD-0.5)*B_O)/R_O))
A2_Y=A1_Y+L_L*COS(((CD-0.5)*B_O)/R_O)
A3_Y=A1_Y+L_O*COS(((CD-0.5)*B_O)/R_O)
;编号
N_ID=(AB-1)*BBBB+CD+SN_ID
D_ID=(AB-1)*L_O*BBBB+(CD-1)*L_O+SD_ID
B_1=D_ID+2
B_2=D_ID+L_O-1
B_3=D_ID+L_O
B_4=D_ID+L_O
B_5=D_ID+1

COMMAND
sel cable id=N_ID beg A1_X A1_Y A1_Z end A2_X A2_Y A2_Z nseg L_L
sel cable id=N_ID beg A2_X A2_Y A2_Z end A3_X A3_Y A3_Z nseg 1
sel cable prop emod 0.2e12 yten KLQD_28 xcarea D_28 gr_per S_28  gr_k SNJQD_O gr_coh SNJNJL_O range cid B_1 B_2
sel cable prop emod 0.2e12 yten KLQD_28 xcarea D_28 gr_per S_28  gr_k SNJQD_O gr_coh SNJNJL_O range cid B_3 B_4
sel cable prop emod 0.2e12 yten KLQD_28 xcarea D_28 gr_per S_28  gr_k SNJQD_O gr_coh SNJNJL_O range cid B_5 B_5
sel cable id N_ID pretension yuyinli_O range cid B_5 B_3
END_COMMAND
   end_loop
   end_loop
SN_ID=SN_ID+AAAA*BBBB
SD_ID=SD_ID+AAAA*BBBB*L_O
   end
CAIBO2
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;



抛砖引玉吧。。。。大概就这个意思。。。学会上面的，锚杆你基本就会了

ca1b000

上面的命令流，我本身是想做预应力锚索的，后来改锚杆的，所以稍微写复杂了。。。写的不好的地方，大家就将就下吧。。。不要喷我哦

lookcity

8#楼的命令流不错呀
就是代号太多 看的有点晕 不过修改参数方便些

ca1b000

上面的参数修改非常方便，直接改第一个函数就可以了。而且里面的id 和cid 是自动编号的。。可以随便插入。算是半自动吧。。。

abot866

学习了，代号太多 看的有点晕 不过修改参数方便

chengsulei

9# ca1b000 请问高手，我想做锚杆支护的命令流应当怎么编写呢》？

漫步听海

8# ca1b000
打个包分享下吧，
zhanggp1987@126.com 谢谢！

要不是人类

高手呀！！现在还看不懂...努力学习中！

chengsulei

可以再解释的详细一些么？