基于FLAC3D的边坡稳定性动力分析——地震工况

基于FLAC3D的边坡稳定性动力分析——地震工况

参考张迎宾教授的文章——论文复现

Paper : Numerical Simulation of Seismic Slope Stability Analysis
Based on Tension-Shear Failure Mechanism

new
config dyn
set dynamic off
restore zhang.sav
def vel
	E=77.88e6
	mu=0.3
	G=E/(2*(1+mu))
	K=E/(3*(1-2*mu))
	rho=2e3
cp=sqrt((K+3/4*G)/rho)
cs=sqrt(G/rho)
thetas_coff=-2*rho*cs*10
end
@vel
model elas 
prop bulk 1e10 shear 3e9 dens 2000
set grav=0,0,-10
fix x range x 0
fix x range x 120
fix y
fix x y z range z 0
solve	
table 1 read 'kobe_correct.txt'
;table 1 read '18s_acc.txt'
call INT.FIS
@integrate(1,2)
call baseline.fis
set @itab_unc=2 @itab_corr=10 @drift=2.60407009e-02 @ttime=15 @itab_cvel=5
@baseline
;table 5 is correct velocity
config dyn
model mech mohr 
prop  bulk @K shea @G cohes=4e4 tens 0 fric=20 dilation 20
ini dens 2000
ini xdis=0 ydis=0 zdis=0 xvel 0 yvel=0 zvel=0
free x y z ran z -0.1 0.1
;apply stress
range name bottom z=-.1 .1
apply sxz @thetas_coff hist table 5 syz 0.0 szz 0.0 range nrange bottom
apply nvel 0 plane norm 0,0,1 range nrange bottom
apply dquiet squiet nquiet range nrange bottom
apply ff
set dyn damp local 0.157
set dyn time = 0
plot create view shear_incr
plot add zcont simaxshear




;;hist set
;plot create view xdisp
;hist reset
;hist add gp xdisp 60 0 40
;hist add gp xdisp 70 0 30
;hist add gp xdisp 79 0 21
;hist add gp xdisp 40 0 20
;hist dytime
;plot add hist 1 2 3 4 vs 5
;plot add zone colorby state
;
;;max_min_priciple_stress
;hist add zone smax 20 0 35
;hist add zone smin 20 0 35
;plot create view max_min_priciple_stress
;plot add hist 6 7 vs 5
;
;
;hist add gp zdisp 60 0 40
;hist add gp zdisp 70 0 30
;hist add gp zdisp 79 0 21
;hist add gp zdisp 40 0 20
;plot create view zdisp
;plot add hist 8 9 10 11 vs 5
;
;;max_min_priciple_stress
;hist add zone smin 40 0 39
;hist add zone smin 40 0 37
;hist add zone smin 40 0 35
;hist add zone smin 40 0 33
;hist add zone smin 40 0 31
;plot create view depth_stress
;plot add hist 12 13 14 15 16 vs 5
;
;;accelerate field
;hist add gp xacc 60 0 40
;hist add gp xacc 70 0 30
;hist add gp xacc 79 0 21
;hist add gp xacc 40 0 20
;plot create view x_accelerate
;plot add hist 17 18 19 20 vs 5
;
;;velocity field
;hist add gp xvel 60 0 40
;hist add gp xvel 70 0 30
;hist add gp xvel 79 0 21
;hist add gp xvel 40 0 20
;plot create view x_velocity
;plot add hist 21 22 23 24 vs 5


set dyn multi on
solve age 18


;def time_calc
;loop i(1,6)
;m=i*3
;file_name=string(m)
;command
;solve age @m
;save @file_name
;endcommand
;endloop
;end
;@time_calc






;
;
;;;buttom_stress
;;hist add zone sxz 60 0 0
;;plot create view bottom_sxz
;;plot add hist 1
;;;plot add table 5
;
;
;;plot create view zone_show
;;plot add zone
;;plot add history location
;hist unbalance

1. 地震波及其基线校正