需要模拟好几十排的预应力锚杆, 按照陈育民提供的三种方法的第二种,用到cid,编一个变量放在cid后边,识别不了。其他几种也要用cid,怎么办??? 所有的几十排上百根锚杆id都是1有没有问题?
def dingzhihu
;左侧第一根锚杆起点、锚固点、终点的x、z值
x_d11=-1.1804
x_d12=-1.7252
x_d13=-2.27
z_d11=2.3166
z_d12=3.3858
z_d13=4.455
;左侧第二根锚杆起点、锚固点、终点的x、z值
x_d21=-0.4067
x_d22=-0.5945
x_d23=-0.7822
z_d21=2.568
z_d22=3.7532
z_d23=4.9384
;右侧第三根锚杆起点、锚固点、终点的x、z值
x_d31=-x_d21
x_d32=-x_d22
x_d33=-x_d23
z_d31=z_d21
z_d32=z_d22
z_d33=z_d23
;右侧第四根锚杆起点、锚固点、终点的x、z值
x_d41=-x_d11
x_d42=-x_d12
x_d43=-x_d13
z_d41=z_d11
z_d42=z_d12
z_d43=z_d13
loop k(1,37)
y_d=k*0.8-0.4
command
sel cable id=1 beg x_d11 y_d z_d11 end x_d12 y_d z_d12 nseg 3 ;自由段,下同
sel cable id=1 beg x_d12 y_d z_d12 end x_d13 y_d z_d13 nseg 3 ;锚固段,下同
sel cable id=1 beg x_d21 y_d z_d21 end x_d22 y_d z_d22 nseg 3
sel cable id=1 beg x_d22 y_d z_d22 end x_d23 y_d z_d23 nseg 3
sel cable id=1 beg x_d31 y_d z_d31 end x_d32 y_d z_d32 nseg 3
sel cable id=1 beg x_d32 y_d z_d32 end x_d33 y_d z_d33 nseg 3
sel cable id=1 beg x_d41 y_d z_d41 end x_d42 y_d z_d42 nseg 3
sel cable id=1 beg x_d42 y_d z_d42 end x_d43 y_d z_d43 nseg 3
endcommand
end_loop
end
def dingfuzhi
loop i(1,148)
cid_1=6*i-5
cid_2=6*i-4
cid_3=6*i-3
cid_4=6*i-2
cid_5=6*i-1
cid_6=6*i
;设置锚杆自由段参数
sel cable id=1 prop emod 2e10 yten 1.03e10 xcarea 3.8e-4 gr_k 1 gr_coh 1 range cid 6*i-4,6*i-3
;设置锚杆锚固段参数
sel cable id=1 prop emod 2e10 yten 1.03e10 xcarea 3.8e-4 gr_k 3.3e6 gr_coh 2.7e5 range cid 6*i-2,6*i
;设置端头参数模拟托盘
sel cable id=1 prop emod 2e10 yten 1.03e10 xcarea 3.8e-4 gr_k 2e10 gr_coh 1e10 range cid 6*i-5,6*i-5
;在锚杆自由段施加预紧力
sel cable id=1 pretension 3e5 range cid 6*i-5,6*i-3
end_loop
end
求大神解怎么办???
例题是这样的:
sel cable id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10
sel cable id=1 beg 0,29,0 end 0,35,0 nseg 6
sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 gr_coh 1 gr_k 1 gr_per 0.0785 range cid 2,10
sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 gr_coh 10e5 gr_k 2e7 range cid 11,17
sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 gr_coh 10e8 gr_k 2e10 range cid 1,1
sel cable id=1 pretension 60e3 range cid 1,10