HSS模型在FLAC3D中的实现及改进 [复制链接]

HSS模型在FLAC3D中的实现及改进
简介：将吸力考虑入经典的饱和土HSS模型，提出了非饱和土的HSS模型，本人命名为HSS_Unsat模型，同时基于该模型提出了非饱和土强度折减理论。
论文在构思之中，先贴出部分计算结果。

图片:no-suction-state.jpg

单元的状态采用FLAC3D的指标习惯，引入了hs屈服，cap屈服以及mc屈服（破坏）

非饱和强度折减法计算的边坡稳定系数。
图为破坏工况的网格变形和位移场。

命令流思路：




new;
title
Slop Stability Analyze by HSS_Unsat

;-------------Config for HSS_Unsat Model----------
config cppudm
model load HSS_UNSAT.dll

;-------------Config Zone Extra memory of saturation----------
Config zextra 1

;-------------The Core of SSR should be loaded earlier----------
call SSRUnsat.dat
call LoadSat.dat

set gravity 0,-9.8,0
set mech ratio 1.0e-5

;------------Import the mesh----------------
IMPGRID grid.Flac3D

;------------Assign the soil material----------------
model HSS_UNSAT
;layer 2
prop  dens 1.85e3 cohesion 16.0e3 fric_q 32.7
prop fric_b 11.7 dialancy 0.0
prop E50_ref 4.8e6. Eur_ref 30.0e6 Ei_ref 150.0e6 Eoed_ref 4.8e6
pro poisson_ur 0.20 door .0001 m_index 1.0 dm 1.0
prop rf 0.75 suction 0.0
prop pc0 50e3 mtype 3


;------------Load Saturation information----------------
load_sat
;----Config the suction by saturation----
set_suction

fix z
fix x range x -0.1 0.1
fix x rang  x 29.9 30.1
fix x y z range y -0.1 0.1

sol

;------------Initialize State----------------
def zero_de1
new_disp = 0.0
p_z = zone_head
loop while p_z # null
if z_model(p_z) = 'HSS_UNSAT'
     z_prop(p_z,'strain1') = 0.0
     z_prop(p_z,'gamma_s') = 0.0
     z_prop(p_z,'e1_p')    = 0.0
endif
p_z = z_next(p_z)
endloop
end

zero_de1

ini xdisp 0.0 ydisp 0.0 zdisp 0.0
ini state 0

def ini_SSRF
_delta = 0.02
_upperK = 5.
_lowerK = 1.
SSR_ratio = 1.0e-4
end
ini_SSRF

;------------Model Null--------------
mod null rang group 3

;------------Apply the surcharge------------
apply syy -20.e3 range x 16.0 21.5 y 14.9 15.1


;----------save the final step information of the model-------------
put_state

;----------get the characteristic step of the system-------------
get_step

_SSRUnsat

save final2.sav

贴一个ssi表示的滑裂面可能更加形象一些。

其他部分计算结果

特别感谢ssw1216 博士提供的强度折减法方面的例子，对本人有很大的启发！

强度折减法还没有在机理上根本解决非饱和土的问题，不过针对工程上的应用来说是可以解决问题的，支持！

引用第6楼geofem于2009-01-17 22:54发表的  :
特别感谢ssw1216 博士提供的强度折减法方面的例子，对本人有很大的启发！

是不是这个

是孙博士的原创，我引用了一下思路。

漂亮
大家一起学习flac

赞。学习。