求教各位,看看是什么问题
代码如下:
1.txt
n
plot add surface green
plot set back white
plot set fore black
plot set center 150 0 50
plot set mag 1
plot add axes black
plot show
gen zone brick p0 0 0 0 p1 220 0 0 p2 0 200 0 p3 0 0 10 &
size 44 50 2 ratio 1 1 1 group ld
attach face range z 9.9 10.1
gen zone brick p0 0 0 10 p1 220 0 10 p2 0 200 10 p3 0 0 11.5 &
size 88 100 2 ratio 1 1 1 group zjdi
attach face range z 11.4 11.6
gen zone brick p0 0 0 11.5 p1 220 0 11.5 p2 0 200 11.5 p3 0 0 16.3 &
size 88 100 6 ratio 1 1 1 group sqcoal
attach face range z 16.2 16.4
gen zone brick p0 0 0 16.3 p1 220 0 16.3 p2 0 200 16.3 p3 0 0 18.3 &
size 88 100 2 ratio 1 1 1 group zjding
attach face range z 18.2 18.4
gen zone brick p0 0 0 18.3 p1 220 0 18.3 p2 0 200 18.3 p3 0 0 32.3 &
size 88 100 7 ratio 1 1 1 group slld
attach face range z 32.2 32.4
gen zone brick p0 0 0 32.3 p1 220 0 32.3 p2 0 200 32.3 p3 0 0 36.3 &
size 44 50 1 ratio 1 1 1 group 6zjd
attach face range z 36.2 36.4
gen zone brick p0 0 0 36.3 p1 220 0 36.3 p2 0 200 36.3 p3 0 0 37.8 &
size 44 50 1 ratio 1 1 1 group 6coal
attach face range z 37.7 37.9
gen zone brick p0 0 0 37.8 p1 220 0 37.8 p2 0 200 37.8 p3 0 0 40.8 &
size 44 50 1 ratio 1 1 1 group 6zjd
attach face range z 40.7 40.9
gen zone brick p0 0 0 40.8 p1 220 0 40.8 p2 0 200 40.8 p3 0 0 49.8 &
size 44 50 1 ratio 1 1 1 group 6ld
attach face range z 49.7 49.9
gen zone brick p0 0 0 49.8 p1 220 0 49.8 p2 0 200 49.8 p3 0 0 58.8 &
size 44 50 1 ratio 1 1 1 group 6ld2
attach face range z 58.7 58.9
gen zone brick p0 0 0 58.8 p1 220 0 58.8 p2 0 200 58.8 p3 0 0 75 &
size 44 50 2 ratio 1 1 1 group fy
sav 1.save
2.txt
res 1.save
m m
set largeprop bulk 1e9 shear 1e9 den 1000 f 10 coh 1e6 range z 11.5 16.3
prop bulk 1e9 shear 1e9 den 1000 f 10 coh 1e6 range z 36.3 37.8
prop bulk 1e9 shear 1e9 den 1000 f 10 coh 1e6 range z 10 11.5
prop bulk 1e9 shear 1e9 den 1000 f 10 coh 1e6 range z 16.3 18.3
prop bulk 1e9 shear 1e9 den 1000 f 10 coh 1e6 range z 32.3 36.3
prop bulk 1e9 shear 1e9 den 1000 f 10 coh 1e6 range z 37.8 40.8prop bulk 1e9 shear 1e9 den 1000 f 10 coh 1e6 range z 0 10
prop bulk 1e9 shear 1e9 den 1000 f 10 coh 1e6 range z 18.3 32.3prop bulk 1e9 shear 1e9 den 1000 f 10 coh 1e6 range z 40.8 49.8
prop bulk 1e9 shear 1e9 den 1000 f 10 coh 1e6 range z 49.8 58.8prop bulk 1e9 shear 1e9 den 1000 f 10 coh 1e6 range z 58.8 75
set g 0 0 -10
fix z range z -.1 .1
fix x range x -.1 .1
fix x range x 219.9 220.1
fix y range y -.1 .1
fix y range y 199.9 200.1apply szz -10e6 range z 74.9 75.1
ini sxx=-10e06 syy=-10e06 szz=-10e06 grad 0 0 2e4
hist unbal
hist gp xdis 155,0,13.25
hist gp xdis 155,20,13.25
hist gp xdis 155,40,13.25
hist gp xdis 155,60,13.25
hist gp xdis 155,80,13.25
hist gp xdis 155,100,13.25
hist gp xdis 155,120,13.25
hist gp xdis 155,140,13.25
hist gp xdis 155,160,13.25
hist gp xdis 155,180,13.25
hist gp xdis 155,200,13.25hist gp zdis 152.5,0,15
hist gp zdis 152.5,20,15
hist gp zdis 152.5,40,15
hist gp zdis 152.5,60,15
hist gp zdis 152.5,80,15
hist gp zdis 152.5,100,15
hist gp zdis 152.5,120,15
hist gp zdis 152.5,140,15
hist gp zdis 152.5,160,15
hist gp zdis 152.5,180,15
hist gp zdis 152.5,200,15
hist gp xdis 150,0,13.25
hist gp xdis 150,20,13.25
hist gp xdis 150,40,13.25
hist gp xdis 150,60,13.25
hist gp xdis 150,80,13.25
hist gp xdis 150,100,13.25
hist gp xdis 150,120,13.25
hist gp xdis 150,140,13.25
hist gp xdis 150,160,13.25
hist gp xdis 150,180,13.25
hist gp xdis 150,200,13.25hist gp zdis 152.5,0,11.5
hist gp zdis 152.5,20,11.5
hist gp zdis 152.5,40,11.5
hist gp zdis 152.5,60,11.5
hist gp zdis 152.5,80,11.5
hist gp zdis 152.5,100,11.5
hist gp zdis 152.5,120,11.5
hist gp zdis 152.5,140,11.5
hist gp zdis 152.5,160,11.5
hist gp zdis 152.5,180,11.5
hist gp zdis 152.5,200,11.5
set mechanical force 100save 2.save
3.txtres 2.save
;形成巷道
ini zdis=0
ini xdis=0
ini ydis=0
m null range x 150 155 y 0 200 z 11.5 15
;巷道支护;锚杆def mgloop a(1,133)
x0=150
x1=147.9
y0=1.5*a-1.5
y1=y0
z0=13.1
z1=z0
command
sel cable begin x0,y0,z0 end x1,y1,z1 nseg 5
sel cable prop emod 4e10 ytens 4e8 &
gr_k=2e7 gr_coh=2.2e7 xcarea=3.14e-4 &
gr_per 6.28e-2
endcommand
endloop
end
mgsolve 运行2小时后,出现上述问题,请教各位如何解决