本人最近做了个预应力锚杆支护的例子，发现支护前后巷道变形量基本没变，求各位大神指导一下，代码如下：
def zhihu
loop a (1,60)
command
;--------------------------------------左帮锚杆--------------------------------------
sel cable id 1 begin 12,a,78.85 end 10.5511,a,78.4618 nseg 10
sel cable id 2 begin 10.5511,a,78.4618 end 10.2613,a,78.3841 nseg 6
sel cable id 1 begin 12,a,79.85 end 10.5,a,79.85 nseg 10
sel cable id 2 begin 10.5,a,79.85 end 10.2,a,79.75 nseg 6
sel cable id 1 begin 12,a,80.85 end 10.5511,a,81.2832 nseg 10
sel cable id 2 begin 10.5511,a,81.2832 end 10.2613,a,81.3159 nseg 6
;--------------------------------------右帮锚杆----------------------------------
sel cable id 1 begin 16.5,a,82.05 end 1,a,82.05 nseg 10
sel cable id 2 begin 18,a,82.05 end 18.3,a,82.05 nseg 6
sel cable id 1 begin 16.5,a,81.05 end 18,a,81.05 nseg 10
sel cable id 2 begin 18,a,81.05 end 18.3,a,81.05 nseg 6
sel cable id 1 begin 16.5,a,80.05 end 18,a,80.05 nseg 10
sel cable id 2 begin 18,a,80.05 end 18.3,a,80.05 nseg 6
sel cable id 1 begin 16.5,a,79.05 end 18,a,79.05 nseg 10
sel cable id 2 begin 18,a,79.05 end 18.3,a,79.05 nseg 6
;--------------------------------------顶板锚杆--------------------------------------
sel cable id 1 begin 12.5424,a,81.1808 end 11.0574,a,82.6657 nseg 10
sel cable id 2 begin 11.0574,a,82.6657 end 10.845,a,82.8778 nseg 6
sel cable id 1 begin 13.3962,a,81.4654 end 12.7321,a,83.4576 nseg 10
sel cable id 2 begin 12.7321,a,83.4576 end 12.6372,a,83.7422 nseg 6
sel cable id 1 begin 14.25,a,81.75 end 13.5859,a,83.7422 nseg 10
sel cable id 2 begin 13.5859,a,83.7422 end 13.4911,a,84.0268 nseg 6
sel cable id 1 begin 15.1038,a,82.0346 end 14.4397,a,82.0268 nseg 10
sel cable id 2 begin 14.4397,a,82.0268 end 14.3449,a,84.3114 nseg 6
sel cable id 1 begin 15.9576,a,82.3192 end 17.0076,a,84.1397 nseg 10
sel cable id 2 begin 17.0076,a,84.1397 end 17.1576,a,84.3977 nseg 6
end_command
end_loop
loop a (1,720)
b=16*a-15
g=16*a-14
h=16*a-6
d=16*a-5
f=16*a
command
sel cable id 1 prop emod 2e10 ytens 3.1e7 xcarea 4.9e-4 gr_coh 10e5 gr_k 2e7 ran cid d f ;设置锚杆锚固段参数
sel cable id 1 prop emod 2e10 ytens 3.1e7 xcarea 4.9e-4 gr_coh 1 gr_k 1 gr_per 0.1 ran cid g h ;设置锚杆自由端参数
sel cable id 1 prop emod 2e10 ytens 3.1e7 xcarea 4.9e-4 gr_coh 10e8 gr_k 2e10 ran cid b b ;模拟托盘
sel cable id 1 pretension 60e3 ran cid b h
end_command
end_loop
end
zhihu
