flac3d基坑开挖支护例子

;part1
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;一区
gen zone radcyl p0 0,0,13.4 p1 6,0,12.41 p2 0,20,13.4 p3 0,0,16.7 &
p4 6,20,12.41 p5 0,20,16.7 p6 6,0,15.13 p7 6,20,15.13 &
p8 2.1,0,13.4 p9 0,0,15.5 p10 2.1,20,13.4 p11 0,20,15.5 &
size 4 20 10 8 ratio 1 1 1 1 group n1
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;二区
gen zone brick p0 2.1,0,12 p1 6,0,11.02 p2 2.1,20,12 p3 2.1,0,13.4 &
p4 6,20,11.02 p5 2.1,20,13.4 p6 6,0,12.41 p7 6,20,12.41 &
size 8 20 5 ratio 1 1 1 group n2
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;三区
gen zone brick p0 -2.1,0,12 p1 6,0,9.78 p2 -2.1,20,12 p3 -1.5,0,12 &
p4 6,20,9.78 p5 -1.5,20,12 p6 6,0,9.935 p7 6,20,9.935 &
size 8 20 1 ratio 1 1 1 group n3
gen zone brick p0 -1.5,0,12 p1 6,0,9.935 p2 -1.5,20,12 p3 1.5,0,12 &
p4 6,20,9.935 p5 1.5,20,12 p6 6,0,10.865 p7 6,20,10.865 &
size 8 20 6 ratio 1 1 1 group n3
gen zone brick p0 1.5,0,12 p1 6,0,10.865 p2 1.5,20,12 p3 2.1,0,12 &
p4 6,20,10.865 p5 2.1,20,12 p6 6,0,11.02 p7 6,20,11.02 &
size 8 20 1 ratio 1 1 1 group n3
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;四区
gen zone brick p0 -6,0,12.99 p1 -2.1,0,12 p2 -6,20,12.99 p3 -6,0,14.39 &
p4 -2.1,20,12 p5 -6,20,14.39 p6 -2.1,0,13.4 p7 -2.1,20,13.4 &
size 8 20 5 ratio 1 1 1 group n14
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;五区
gen zone radcyl p0 0,0,13.4 p1 0,0,16.7 p2 0,20,13.4 p3 -6,0,14.39 &
p4 0,20,16.7 p5 -6,20,14.39 p6 -6,0,18.35 p7 -6,20,18.35 &
p8 0,0,15.5 p9 -2.1,0,13.4 p10 0,20,15.5 p11 -2.1,20,13.4 &
size 4 20 10 8 ratio 1 1 1 1 group n11
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;六区 （？？顶 喷层？？？）
gen zon cshell p0 0,0,13.4 p1 2.1,0,13.4 p2 0,20,13.4 p3 0,0,15.5 &
dim 2.0 2.0 2.0 size 1 20 10 group line
gen zon cshell p0 0,0,13.4 p1 0,0,15.5 p2 0,20,13.4 p3 -2.1,0,13.4 &
dim 2.0 2.0 2.0 size 1 20 10 group line
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;七区（？？侧帮 喷层？？？）
gen zone brick p0 -2.1,0,12 p1 -2.0,0,12 p2 -2.1,20,12 p3 -2.1,0,13.4 &
p4 -2.0,20,12 p5 -2.1,20,13.4 p6 -2.0,0,13.4 p7 -2.0,20,13.4 &
size 1 25 5 ratio 1 1 1 group line
gen zone brick p0 2.0,0,12 p1 2.1,0,12 p2 2.0,20,12 p3 2.0,0,13.4 &
p4 2.1,20,12 p5 2.0,20,13.4 p6 2.1,0,13.4 p7 2.1,20,13.4 &
siz

e 1 20 5 ratio 1 1 1 group line
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;八区
gen zone brick p0 -2.0,0,12 p1 2.0,0,12 p2 -2.0,20,12 p3 -2.0,0,13.4 &
p4 2.0,20,12 p5 -2.0,20,13.4 p6 2.0,0,13.4 p7 2.0,20,13.4 &
size 8 20 5 ratio 1 1 1 group t
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;九区
gen zon cshell p0 0,0,13.4 p1 1,0,13.4 p2 0,20,13.4 p3 0,0,14.4 &
dim 0.5 0.5 0.5 size 1 20 10 group t
gen zon cshell p0 0,0,13.4 p1 0,0,14.4 p2 0,20,13.4 p3 -1,0,13.4 &
dim 0.5 0.5 0.5 size 1 20 10 group t
gen zon cshell p0 0,0,13.4 p1 1.5,0,13.4 p2 0,20,13.4 p3 0,0,14.9 &
dim 1 1 1 size 1 20 10 group t
gen zon cshell p0 0,0,13.4 p1 0,0,14.9 p2 0,20,13.4 p3 -1.5,0,13.4 &
dim 1 1 1 size 1 20 10 group t
gen zon cshell p0 0,0,13.4 p1 2,0,13.4 p2 0,20,13.4 p3 0,0,15.4 &
dim 1.5 1.5 1.5 size 1 20 10 group t
gen zon cshell p0 0,0,13.4 p1 0,0,15.4 p2 0,20,13.4 p3 -2,0,13.4 &
dim 1.5 1.5 1.5 size 1 20 10 group t

gen zon cylinder p0 0,0,13.4 p1 0.5,0,13.4 p2 0,20,13.4 p3 0,0,13.9 &
dim 0.5 0.5 0.5 size 1 20 10 group t
gen zon cylinder p0 0,0,13.4 p1 0,0,13.9 p2 0,20,13.4 p3 -0.5,0,13.4 &
dim 0.5 0.5 0.5 size 1 20 10 group t
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;十区
gen zone brick p0 6,0,9.78 p1 15,0,7.36 p2 6,20,9.78 p3 6,0,11.02 &
p4 15,20,7.36 p5 6,20,11.02 p6 15,0,8.6 p7 15,20,8.6 &
size 12 20 8 ratio 1.05 1 1 group n3
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;十一区
gen zone brick p0 6,0,11.02 p1 15,0,8.6 p2 6,20,11.02 p3 6,0,12.41 &
p4 15,20,8.6 p5 6,20,12.41 p6 15,0,10 p7 15,20,10 &
size 12 20 5 ratio 1.05 1 1 group n2
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;十二区
gen zone brick p0 6,0,12.41 p1 15,0,10 p2 6,20,12.41 p3 6,0,15.13 &
p4 15,20,10 p5 6,20,15.13 p6 15,0,12.72 p7 15,20,12.72 &
size 12 20 5 ratio 1.05 1 1 group n1
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;十三区
gen zone brick p0 -15,0,15.4 p1 -6,0,12.99 p2 -15,20,15.4 p3 -15,0,16.8 &
p4 -6,20,12.99 p5 -15,20,16.8 p6 -6,0,14.39 p7 -6,20,14.39 &
size 8 25 5 ratio 0.95 1 1 group n3
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;十四区
gen zone brick p0 -15,0,16.8 p1 -6,0,14.39 p2 -15,20,16.8 p3 -15,0,18.04 &
p4 -6,20,14.39 p5 -15,20,18.04 p6 -6,0,15.182 p7 -6,20,15.182 &
size 8 20 1 ra

tio 0.95 1 1 group n2
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;十五区
gen zone brick p0 -15,0,18.04 p1 -6,0,15.182 p2 -15,20,18.04 p3 -15,0,20.76 &
p4 -6,20,15.182 p5 -15,20,20.76 p6 -6,0,18.35 p7 -6,20,18.35 &
size 8 20 4 ratio 0.95 1 1 group n1
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;十六区
gen zone brick p0 -6,0,9.42 p1 0,0,7.82 p2 -6,20,9.42 p3 -6,0,12.99 &
p4 0,20,7.82 p5 -6,20,12.99 p6 -2.1,0,12 p7 -2.1,20,12 &
size 8 20 4 ratio 1 1 0.9 group n12
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;十七区
gen zone brick p0 0,0,7.82 p1 6,0,6.21 p2 0,20,7.82 p3 -2.1,0,12 &
p4 6,20,6.21 p5 -2.1,20,12 p6 6,0,9.78 p7 6,20,9.78 &
size 8 20 4 ratio 1 1 0.9 group n13
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;十八区
gen zone brick p0 6,0,6.21 p1 15,0,3.8 p2 6,20,6.21 p3 6,0,9.78 &
p4 15,20,3.8 p5 6,20,9.78 p6 15,0,7.36 p7 15,20,7.36 &
size 12 20 4 ratio 1.05 1 0.9 group n4
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;十九区
gen zone brick p0 -15,0,11.8 p1 -6,0,9.42 p2 -15,20,11.8 p3 -15,0,15.4 &
p4 -6,20,9.42 p5 -15,20,15.4 p6 -6,0,12.99 p7 -6,20,12.99 &
size 8 20 4 ratio 0.95 1 0.9 group n4
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;二十区
gen zone brick p0 -6,0,0 p1 0,0,0 p2 -6,20,0 p3 -6,0,9.42 &
p4 0,20,0 p5 -6,20,9.42 p6 0,0,7.82 p7 0,20,7.82 &
size 8 20 6 ratio 1 1 0.9 group n5
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;二十一区
gen zone brick p0 0,0,0 p1 6,0,0 p2 0,20,0 p3 0,0,7.82 &
p4 6,20,0 p5 0,20,7.82 p6 6,0,6.21 p7 6,20,6.21 &
size 8 20 6 ratio 1 1 0.9 group n5
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;二十二区
gen zone brick p0 6,0,0 p1 15,0,0 p2 6,20,0 p3 6,0,6.21 &
p4 15,20,0 p5 6,20,6.21 p6 15,0,3.8 p7 15,20,3.8 &
size 12 20 6 ratio 1.05 1 0.9 group n5
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;二十三区
gen zone brick p0 -15,0,0 p1 -6,0,0 p2 -15,20,0 p3 -15,0,11.8 &
p4 -6,20,0 p5 -15,20,11.8 p6 -6,0,9.42 p7 -6,20,9.42 &
size 8 20 6 ratio 0.95 1 0.9 group n5
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;二十四区
gen zone brick p0 -6,0,19.6 p1 0,0,17.95 p2 -6,20,19.6 p3 -6,0,23.93 &

p4 0,20,17.95 p5 -6,20,23.93 p6 0,0,22.33 p7 0,20,22.33 &
size 5 20 4 ratio 1 1 1.1 group n6
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;二十五区
gen zone brick p0 0,0,17.95 p1 6,0,16.38 p2 0,20,17.95 p3 0,0,22.33 &
p4 6,20,16.38 p5 0,20,22.33 p6 6,0,20.72 p7 6,20,20.72 &
size 5 20 4 ratio 1 1 1.1 group n6
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;二十六区
gen zone brick p0 6,0,16.38 p1 15,0,13.97 p2 6,20,16.38 p3 6,0,20.72 &
p4 15,20,13.97 p5 6,20,20.72 p6 15,0,18.31 p7 15,20,18.31 &
size 12 20 4 ratio 1.05 1 1.1 group n6
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;二十七区
gen zone brick p0 -15,0,22.01 p1 -6,0,19.6 p2 -15,20,22.01 p3 -15,0,26.35 &
p4 -6,20,19.6 p5 -15,20,26.35 p6 -6,0,23.93 p7 -6,20,23.93 &
size 8 20 4 ratio 0.95 1 1.1 group n6
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;c1煤层
gen zone brick p0 -6,0,18.35 p1 0,0,16.7 p2 -6,20,18.35 p3 -6,0,19.6 &
p4 0,20,16.7 p5 -6,20,19.6 p6 0,0,17.95 p7 0,20,17.95 &
size 5 20 2 ratio 1 1 1.1 group c1
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;c1
gen zone brick p0 0,0,16.7 p1 6,0,15.13 p2 0,20,16.7 p3 0,0,17.95 &
p4 6,20,15.13 p5 0,20,17.95 p6 6,0,16.38 p7 6,20,16.38 &
size 5 20 2 ratio 1 1 1.1 group c1
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;c1
gen zone brick p0 6,0,15.13 p1 15,0,12.72 p2 6,20,15.13 p3 6,0,16.38 &
p4 15,20,12.72 p5 6,20,16.38 p6 15,0,13.97 p7 15,20,13.97 &
size 12 20 2 ratio 1.05 1 1.1 group c1
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;c1
gen zone brick p0 -15,0,20.76 p1 -6,0,18.35 p2 -15,20,20.76 p3 -15,0,22.01 &
p4 -6,20,18.35 p5 -15,20,22.01 p6 -6,0,19.6 p7 -6,20,19.6 &
size 8 20 2 ratio 0.95 1 1.1 group c1
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;二十八区
gen zone brick p0 -6,0,23.93 p1 0,0,22.33 p2 -6,20,23.93 p3 -6,0,30 &
p4 0,20,22.33 p5 -6,20,30 p6 0,0,30 p7 0,20,30 &
size 5 20 6 ratio 1 1 1.1 group n9
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;二十九区
gen zone brick p0 0,0,22.33 p1 6,0,20.72 p2 0,20,22.33 p3 0,0,30 &
p4 6,20,20.72 p5 0,20,30 p6 6,0,30 p7 6,20,30 &
size 5 20 6 ratio 1 1

1.1 group n9
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;三十区
gen zone brick p0 6,0,20.72 p1 15,0,18.31 p2 6,20,20.72 p3 6,0,30 &
p4 15,20,18.31 p5 6,20,30 p6 15,0,30 p7 15,20,30 &
size 12 20 6 ratio 1.05 1 1.1 group n9
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;三十一区
gen zone brick p0 -15,0,26.35 p1 -6,0,23.93 p2 -15,20,26.35 p3 -15,0,30 &
p4 -6,20,23.93 p5 -15,20,30 p6 -6,0,30 p7 -6,20,30 &
size 8 20 6 ratio 0.95 1 1.1 group n9
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
mod e
pro bulk 1e10 she 1e10
ini density 2000 range group n1
ini density 2550 range group n2
ini density 1650 range group n3
ini density 2100 range group t
ini density 2620 range group n4
ini density 2680 range group n5
ini density 2500 range group n6
ini density 2650 range group n9
ini density 2680 range group line
ini density 1650 range group c1
ini density 2000 range group n11
ini density 2200 range group n13
ini density 1650 range group n14
set gravity 0 0 -9.81
set large
;;;;;;;;;;;;;;
ini szz=-27.1e6
ini sxx=-30.61e6
ini syy=-29.89e6
; boundary conditions
apply szz -27.1e6 range z 29.9,30.1
fix x range x -15.1,-14.9
fix x range x 14.9,15.1
fix y range y -0.1,0.1
fix y range y 19.9,20.1
fix x range z -0.1,0.1
fix y range z -0.1,0.1
fix z range z -0.1,0.1
; equilibrate
hist unbal
solve
sav 弹性模型.sav


rest 弹性模型.sav
model mohr
pro bulk 4e9 she 2.5e9 fric 18 coh 8e5 ten 5e4 range group n1
pro bulk 3e9 she 2e9 fric 18 coh 2e5 ten 1e3 range group c1
pro bulk 5e9 she 2.3e9 fric 45 coh 20e17 ten 3e5 range group n2
pro bulk 5e9 she 2.3e9 fric 28 coh 1.2e6 ten 3e5 range group t
pro bulk 3e9 she 2e9 fric 20 coh 5e5 ten 1e3 range group n3
pro bulk 5e9 she 4e9 fric 20 coh 4.5e5 ten 6e5 range group n4
pro bulk 6e9 she 4e9 fric 22 coh 3e5 ten 2e6 range group n5
pro bulk 5e9 she 3e9 fric 16 coh 6.0e5 ten 1e4 range group line
pro bulk 5e9 she 3.0e9 fric 15 coh 3e5 ten 8e5 range group n6
pro bulk 6e9 she 3.6e9 fric 20 coh 4e5 ten 1e6 range group n9
pro bulk 4e9 she 2.5e9 fric 17 coh 5e5 ten 5e4 range group n11
pro bulk 3e9 she 2e9 fric 14 coh 1e5 ten 1e3 range group n12
pro bulk 3e9 she 2e9 fric 18 coh 3.5e5 ten 1e3 range group n13
pro bulk 3e9 she 2e9 fric 22 coh 1.5e6 ten 1e3 range group n14
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
ini xd 0 yd 0 zd 0
ini xv 0 yv 0 zv 0
hist unbal
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;开挖;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
model null range group t
;;;;;;;;;;;;;;;;;;;;;锚杆;;;;;;;;;;;;;;;;;;;
def ins_cable1
loop n (0,28)
y1=0.35
y2=y1+0.7*n
command
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;NO 2
sel cable id=1 begin -2.1 y2 13.2118 end -4.1 y2 13.2118 nseg 4
;;;;;;;;;;;;;;
sel ca

ble id=1 begin -2.0376 y2 13.9081 end -3.9782 y2 14.3919 nseg 4
sel cable id=1 begin -1.7571 y2 14.5500 end -3.4306 y2 15.6453 nseg 4
;;
sel cable id=1 begin -1.3037 y2 15.0463 end -2.5454 y2 16.6142 nseg 4
;;;;
sel cable id=1 begin -0.6902 y2 15.3833 end -1.3476 y2 17.2722 nseg 4
;;;;;;;;;;;;;;
sel cable id=1 begin 0 y2 15.5 end 0 y2 17.5 nseg 4
;;;;;;;;;;;;;;;;;;
sel cable id=2 begin 0.6902 y2 15.3833 end 1.3476 y2 17.2722 nseg 4
;;;;;;;;;;;;;;;
sel cable id=2 begin 1.3037 y2 15.0463 end 2.5454 y2 16.6142 nseg 4
;;;;;;;;;;;;;;;;;;
sel cable id=2 begin 1.7571 y2 14.5500 end 3.4306 y2 15.6453 nseg 4
;;;;;;;;;;;;;;;;;;;
sel cable id=2 begin 2.0376 y2 13.9081 end 3.9782 y2 14.3919 nseg 4
;;;;;;;;;;;;;;;
sel cable id=2 begin 2.1 y2 13.2118 end 4.1 y2 13.2118 nseg 4
;;;;;;;;;;;;;;;;;
sel cable id=2 begin 2.1 y2 12.28 end 3.1 y2 10.5479 nseg 4
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
sel cable prop emod 6e4 xcarea 1.57e-3 gr_per 1. yten 15e4 &
gr_k 1e5 gr_c 10e6 ran id=1
sel cable prop emod 6e9 xcarea 1.57e-3 gr_per 1. yten 15e4 &
gr_k 6e12 gr_c 10e4 ran id=2
end_command
end_loop
end
ins_cable1
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;锚索;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
def ins_cable2
loop n (0,5)
y1=0.6
y2=y1+3.5*n
command
sel cable id=5 begin -0.9712 y2 15.2619 end -3.4225 y2 19.9610 nseg 1
sel cable id=5 begin 0.9712 y2 15.2619 end 3.4225 y2 19.9610 nseg 1
sel cable id=4 begin -3.4225 y2 19.9610 end -4.3475 y2 21.7342 nseg 5
sel cable id=3 begin 3.4225 y2 19.9610 end 4.3475 y2 21.7342 nseg 5
sel cable prop emod 3e8 xcarea 1.57e-3 gr_per 1. yten 1e-3 &
gr_k 1e-3 gr_c 1e-3 ran id=5
sel cable prop emod 1e8 xcarea 1.57e-3 gr_per 1. yten 25e5 &
gr_k 1e5 gr_c 8e5 ran id=3
sel cable prop emod 8e8 xcarea 1.57e-3 gr_per 1. yten 25e5 &
gr_k 8e5 gr_c 8e5 ran id=4
endcommand
endloop
end
ins_cable2
;;;;;;;;;;;;;;;;;;;;;;;;;;;;喷层;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
sel shell range group line
sel shell prop isotropic 3e9,0.25 thickness 0.1 density 2500
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;;;;;;;;;;;;;;;;;;;;;底板中央控制点位移;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
hist gp zdis 0 0 11.8
hist gp xdis 0 0 11.8
;;;;;;;;;;;;;;;;;;;;;;左下帮控制点位移;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
hist gp zdis -2.35 0 12.2
hist gp xdis -2.35 0 12.2
;;;;;;;;;;;;;;;;;;;;;;顶板右侧控制点侧位移;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
hist gp zdis 1.3 0 15.4
hist gp xdis 1.3 0 15.4
;;;;;;;;;;;;;;;;;;;;;;顶板左侧控制点侧位移;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
hist gp zdis -1.6 0 15.04
hist gp xdis -1.6 0 15.04
step 500
save 500.sav
rest 500.sav

step 500
save 1000.sav
rest 1000.sav
step 500
save 1500.sav
rest 1500.sav
step 500
save 2000.sav
rest 2000.sav
step 500
save 2500.sav
rest 2500.sav
step 500
save 3000.sav
rest 3000.sav
step 500
save 3500.sav
rest 3500.sav
step 500
save 4000.sav
rest 4000.sav
step 500
save 4500.sav
rest 4500.sav
step 500
save 5000.sav
rest 5000.sav
solve
sav 原支护.sav

运行到红色字体的部分时出现了错误“zone 6682 does not conform to node numbering convention”
请高手指点一下，这个程序还很不成熟，有什么建议也请指点，不胜感激！