结果是这个样子，是不是少东西啊，求指导？？？
命令流如下：
new
config cppudm
gen zon bri p1 .6 0 0 p2 0 .6 0 p3 0 0 .6 size 1 1 1
;model elastic
;prop bulk 3e8 shear 1e8
;model load E:\STUDY\VC\FLAC\DuncanChang-2d\Release\DuncanChang.dll;; Your Dll file location
;model duncan
model load DuncanChang.dll
model duncan
prop cohesion 110e3 friction 48.5 fricDel 0.0 ratiofail 0.79 ke 704 ne 0.38 kb 303 mb 0.18 kur 844.8
;三轴试验尺寸：0.6*0.6*0.6 2个单元
;采用I过度料饱和样的试验结果kur暂时未考虑
fix z ran z -.01 .01
;加300kPa的围压
def sigma3
sigma3=-20e3
end
sigma3
app nstress sigma3 ran x -.01 .01
app nstress sigma3 ran x .59 .61
app nstress sigma3 ran y -.01 .01
app nstress sigma3 ran y .59 .61
app nstress sigma3 ran z .59 .61
ini den 2190
;设置初始应力
ini szz sigma3
ini syy sigma3
ini sxx sigma3
set grav 0 0 -10
;solve
ini xd 0 yd 0 zd 0
ini xv 0 yv 0 zv 0
hist id=1 gp zd 0 0 0.6
hist id=2 gp xd 0 0 0.6
hist id=3 gp yd 0 0 0.6
hist id=11 zone szz .3 .3 .6
hist id=12 zone syy .3 .3 .6
hist id=13 zone sxx .3 .3 .6
;plo hist -11 vs -1
tab 1 name loads
tab 2 name ev_ea
;第1次加卸载
def load1
p_gp = gp_near(0,0,0.6)
loop n(1,50)
zss_load=sigma3 - float(n)*5e3
z_dis = -1*gp_zdisp(p_gp) / 0.6
z_load = (sigma3 - zss_load) / 5e3
command
app nstress zss_load ran z .59 .61
solve
tab 1 z_dis z_load
end_command
end_loop
end
load1
plo tab 1
def unload1
p_gp = gp_near(0,0,0.6)
loop m(1,50)
zss_load=-1600e3 + float(m)*10e3
z_dis = gp_zdisp(p_gp)
command
app nstress zss_load ran z .59 .61
solve
tab 1 zss_load z_dis
endcommand
endloop
end
unload1
;第2次加卸载
def load2
p_gp = gp_near(0,0,0.6)
loop n(1,150)
zss_load=-1100e3 - float(n)*10e3
z_dis = gp_zdisp(p_gp)
command
app nstress zss_load ran z .59 .61
solve
tab 1 zss_load z_dis
end_command
end_loop
end
load2
def unload2
p_gp = gp_near(0,0,0.6)
loop m(1,150)
zss_load=-2600e3 + float(m)*10e3
z_dis = gp_zdisp(p_gp)
command
app nstress zss_load ran z .59 .61
solve
tab 1 zss_load z_dis
endcommand
endloop
end
unload2
;第3次加卸载
def load3
p_gp = gp_near(0,0,0.6)
loop n(1,250)
zss_load=-1100e3 - float(n)*10e3
z_dis = gp_zdisp(p_gp)
command
app nstress zss_load ran z .59 .61
solve
tab 1 zss_load z_dis
end_command
end_loop
end
load3
def unload3
p_gp = gp_near(0,0,0.6)
loop m(1,250)
zss_load=-3600e3 + float(m)*10e3
z_dis = gp_zdisp(p_gp)
command
app nstress zss_load ran z .59 .61
solve
tab 1 zss_load z_dis
endcommand
endloop
end
unload3
;第4次加卸载
def load4
p_gp = gp_near(0,0,0.6)
loop n(1,280)
zss_load=-1100e3 - float(n)*10e3
z_dis = gp_zdisp(p_gp)
command
app nstress zss_load ran z .59 .61
solve
tab 1 zss_load z_dis
end_command
end_loop
end
load4
save sig600_un.sav
