cc0714yantu |
2016-11-17 13:52 |
结果是这个样子,是不是少东西啊,求指导???[attachment=305671]
命令流如下:
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
|
|