根据《FLAC3D在岩土工程中的应用》书中的命令流进行模拟，为什么得到的图形和书上的不符合，求高人指点。命令流如下：
new
def _sol
step1=50000
flag01=0
;loop while flag01<=9
flag001=flag01*10
file00=string(int(flag01))+'0'+'.flac3d'
file01=string(int(flag01))+'.sav'
command
;delete
impgrid file00
group rock range group 1
group weakplane range group 2
endcommand
; 弹性模量(Pa)
E1=0.20e9 ;节理
E3=1.63e9
;
; 泊松比
poi1=0.30
poi3=0.21
;
; 粘结力(Pa)
coh1=100.0e3
coh3=2000.0e3
;
; 容重(N/m3)
weight1=21.0e3
weight3=26.80e3
;
; 膨胀角(度)
dila1=5.0
dila3=10.0
;
; 内摩擦角(度)
fri1=23.0
fri3=41.0
;
;抗拉强度
ten1=0.10e6
ten3=1.41e6
grav0=-9.80 ; 重力加速度(N/kg)
dens1=-weight1/grav0
dens3=-weight3/grav0
K1=E1/(3*(1-2*poi1))
K3=E3/(3*(1-2*poi3))
G1=E1/(2*(1+poi1))
G3=E3/(2*(1+poi3))
command
;model null ;挖空
model ss ;采用Mohr-Coulomb应变软化模型
pro bulk K1 she G1 dens dens1 coh coh1 dilation dila1 friction fri1 tens ten1 ftab 1 ctab 2 dtab 3 range group weakplane
pro bulk K3 she G3 dens dens3 coh coh3 dilation dila3 friction fri3 tens ten3 ftab 4 ctab 5 dtab 6 range group rock
table 1 0 23 0.05 20 0.1 18 1 18
table 4 0 41 0.05 38 0.1 35 1 35
table 2 0 100.0e3 0.05 50.0e3 0.1 25.0e3 1 25.0e3
table 5 0 2000.0e3 0.05 1000.0e3 0.1 500.0e3 1 500.0e3
table 3 0 5.0 0.05 3.0 0.1 1.0 1 1.0
table 6 0 10.0 0.05 6.0 0.1 3.0 1 3.0
;
fix y range y -.1 .1
fix y range y 3.9 4.1
ini yvel 2.e-6 range y -.1 .1
ini yvel -2.e-6 range y 3.9 4.1
endcommand
end
_sol
def ax_str
str = 0
pnt = gp_head
loop while pnt # null
if gp_ypos(pnt) < 0.1 then
str = str + gp_yfunbal(pnt)
endif
pnt = gp_next(pnt)
endloop
ax_str = str / pi ; cylinder radius = 1
end
hist n 1
hist gp ydisp 0,0,0
hist ax_str
hist gp xdisp 1,1,0
hist id 17 unbal
plot hist -2 vs 1 ;axial stress vs axial disp.
step step1
def _sav
command
save file01
endcommand
end
_sav
;plo crea qqq
;plo add surf red range group yield
;plo add axes green
;plo set rot 123 313 3
;plo set mag 1.5
hist write -2 vs 1 skip 200 file 0.txt
得到的图形如下：
为什么位移这么小，数量级和书上的差一个，应力的数量级也和书上的差一个？另外绘图时坐标的间距是如何选取的？就像图中X轴的2.0 4.0 6.0
