new
generate zone radcyl p0 = (0,0,0) p1 = (20,0,0.996) &
p2 = (0,1,0) p3 = (0,0,20) p4 = (20,1,0.996) p5 = (0,1,20) &
p6 = (20,0,20) p7 = (20,1,20) p8 = (2.4,0,0.996) p9 = (0,0,2.596) p10 = (2.4,1,0.996) &
p11 = (0,1,2.596) size 10 5 10 10 dimension 2.596 2.596 2.596 2.596 ratio 1 1 1 1 fill group tunnel
gen zone b p0 2.4 0 0 p1 20 0 0 p2 2.4 1 0 p3 2.4 0 0.996 size 10 10 1
gen zone c p0 0 0 0 p1 2.4 0 0 p2 0 1 0 p3 2.4 0 0.996 p4 2.4 1 0 p5 2.4 1 0.996 fill group tunnel2
gen zone radtunnel p0 = (0,0,0) p1 = (0,0,-20) &
p2 = (0,1,0) p3 = (20,0,0) &
size 10 5 10 10 dimension 1.204,2.4,1.204 fill group tunnel3
gen zone reflect normal (-1,0,0) origin =(0,0,0)
group section1 range y 0 1 group tunnel
group section range y 0 1 group tunnel2
group section range y 0 1 group tunnel3
save model.sav
;德鲁克-普拉格模型的材料参数
model dr
prop bulk 1.41e9 shear 0.97e9 kshear 1.603e7 qvol 0.917
;给定边界条件
fix x range x -20.1 -19.9
fix x range x 19.9 20.1
fix z range z -20.1 -19.9
fix y range y -0.1 0.1
fix y range y 0.9 1.1
;初始应力赋值
set grav 0 0 -10
ini density 2750
ini szz -17.6e6 grad 0 0 27500
ini sxx -17.6e6 grad 0 0 27500
ini syy -17.6e6 grad 0 0 27500
mod null range group section1
solve
save kaiwa1.sav
pau
;第一排
sel cable id=1 begin 0 0 2.596 end 0 0 4.596 nseg 20 ;建立锚杆自由段
sel cable id=1 begin 0 0 4.596 end 0 0 5.096 nseg 5 ;建立锚杆锚固段
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1 gr_k 1 gr_p 0.157 range cid 2 20 ;设置自由端参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 21 25 ;设置锚固段参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 1 1 ;设置端头参数模拟托盘
sel cable id=1 pretension 100e3 range cid 1 20 ;在锚杆自由段施加预应力
;第二排
sel cable id=2 begin 0.75 0 2.48 end 1.33 0 4.40 nseg 20
sel cable id=2 begin 1.33 0 4.40 end 1.48 0 4.88 nseg 5
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1 gr_k 1 gr_p 0.157 range cid 27 45 ;设置自由端参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 46 50 ;设置锚固段参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 26 26 ;设置端头参数模拟托盘
sel cable id=2 pretension 100e3 range cid 26 45 ;在锚杆自由段施加预应力
;第三排
sel cable id=3 begin 1.44 0 2.16 end 2.55 0 3.82 nseg 20
sel cable id=3 begin 2.55 0 3.82 end 2.83 0 4.24 nseg 5
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1 gr_k 1 gr_p 0.157 range cid 52 70 ;设置自由端参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 71 75 ;设置锚固段参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 51 51 ;设置端头参数模拟托盘
sel cable id=3 pretension 100e3 range cid 51 70 ;在锚杆自由段施加预应力
;第四排
sel cable id=4 begin 2 0 1.65 end 3.55 0 2.92 nseg 20
sel cable id=4 begin 3.55 0 2.92 end 3.93 0 3.24 nseg 5
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1 gr_k 1 gr_p 0.157 range cid 77 95 ;设置自由端参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 96 100 ;设置锚固段参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 76 76 ;设置端头参数模拟托盘
sel cable id=4 pretension 100e3 range cid 76 95 ;在锚杆自由段施加预应力
;第五排
sel cable id=5 begin 2.4 0 0.996 end 4.4 0 0.996 nseg 20
sel cable id=5 begin 4.4 0 0.996 end 4.9 0 0.996 nseg 5
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1 gr_k 1 gr_p 0.157 range cid 102 120 ;设置自由端参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 121 125 ;设置锚固段参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 101 101 ;设置端头参数模拟托盘
sel cable id=5 pretension 100e3 range cid 101 120
;第六排
sel cable id=6 begin 2.4 0 0.196 end 4.4 0 0.196 nseg 20
sel cable id=6 begin 4.4 0 0.196 end 4.9 0 0.196 nseg 5
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1 gr_k 1 gr_p 0.157 range cid 127 145 ;设置自由端参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 146 150 ;设置锚固段参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 126 126 ;设置端头参数模拟托盘
sel cable id=6 pretension 100e3 range cid 126 145
;第七排
sel cable id=7 begin 2.4 0 -0.604 end 4.4 0 -0.604 nseg 20
sel cable id=7 begin 4.4 0 -0.604 end 4.9 0 -0.604 nseg 5
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1 gr_k 1 gr_p 0.157 range cid 152 170 ;设置自由端参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 171 175 ;设置锚固段参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 151 151 ;设置端头参数模拟托盘
sel cable id=7 pretension 100e3 range cid 151 170
;第八排
sel cable id=8 begin -0.75 0 2.48 end -1.33 0 4.40 nseg 20
sel cable id=8 begin -1.33 0 4.40 end -1.48 0 4.88 nseg 5
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1 gr_k 1 gr_p 0.157 range cid 177 195 ;设置自由端参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 196 200 ;设置锚固段参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 176 176 ;设置端头参数模拟托盘
sel cable id=8 pretension 100e3 range cid 176 195
;第九排
sel cable id=9 begin -1.44 0 2.16 end -2.55 0 3.82 nseg 20
sel cable id=9 begin -2.55 0 3.82 end -2.83 0 4.24 nseg 5
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1 gr_k 1 gr_p 0.157 range cid 202 220 ;设置自由端参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 221 225 ;设置锚固段参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 201 201 ;设置端头参数模拟托盘
sel cable id=9 pretension 100e3 range cid 201 220
;第十排
sel cable id=10 begin -2 0 1.65 end -3.55 0 2.92 nseg 20
sel cable id=10 begin -3.55 0 2.92 end -3.93 0 3.24 nseg 5
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1 gr_k 1 gr_p 0.157 range cid 227 245 ;设置自由端参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 246 250 ;设置锚固段参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 226 226 ;设置端头参数模拟托盘
sel cable id=10 pretension 100e3 range cid 226 245
;第十一排
sel cable id=11 begin -2.4 0 0.996 end -4.4 0 0.996 nseg 20
sel cable id=11 begin -4.4 0 0.996 end -4.9 0 0.996 nseg 5
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1 gr_k 1 gr_p 0.157 range cid 252 270 ;设置自由端参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 271 275 ;设置锚固段参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 251 251 ;设置端头参数模拟托盘
sel cable id=11 pretension 100e3 range cid 251 270
;第十二排
sel cable id=12 begin -2.4 0 0.196 end -4.4 0 0.196 nseg 20
sel cable id=12 begin -4.4 0 0.196 end -4.9 0 0.196 nseg 5
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1 gr_k 1 gr_p 0.157 range cid 277 295 ;设置自由端参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 296 300 ;设置锚固段参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 276 276 ;设置端头参数模拟托盘
sel cable id=12 pretension 100e3 range cid 276 295
;第十三排
sel cable id=13 begin -2.4 0 -0.604 end -4.4 0 -0.604 nseg 20
sel cable id=13 begin -4.4 0 -0.604 end -4.9 0 -0.604 nseg 5
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1 gr_k 1 gr_p 0.157 range cid 302 320 ;设置自由端参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 321 326 ;设置锚固段参数
sel cable prop emod 20e10 yten 600e10 xcarea 1.9625e-3 &
gr_coh 1e20 gr_k 2e10 range cid 301 301 ;设置端头参数模拟托盘
sel cable id=13 pretension 100e3 range cid 301 320
solve
save maogan13.save
plot set rotation 20 0 30
plot set center auto
plot add surface blue
plot add axes
plot show
这个是个巷道的开挖支护模型，出来没有塑性区，然后锚杆也不起作用，谁能告诉我怎么回事？
