渗流稳定性分析（MATLAB实现）

function SSA;
%%渗流稳定性分析（Seepage stability analysis）
%%  C.0.2 悬挂式截水帷幕底端位于碎石土、砂土或粉土含水层时，对均质含水层，地下水渗流的流土稳定性应符合下式规定（图C.0.2）：
%% （2D+0.8D1）γ'/(Δhγw) ≥Kse      （C.0.2）
%%  止水帷幕底面至坑底的土层厚度：D≥ (KseΔhγw/γ'-0.8D1)

%%  式中： Kse──流土稳定性安全系数；安全等级为一、二、三级的支护结构，Kse分别不应小于1.6、1.5、1.4；
%%  D──截水帷幕底面至坑底的土层厚度(m)；
%%  D1──潜水水面或承压水含水层顶面至基坑底面的土层厚度(m)；
%%  γ＇──土的浮重度(kN/m3) γ＇=gamma_s；
%%  Δh──基坑内外的水头差(m) Δh =delta_h；
%%  γw──水的重度(kN/m3)  γw=gamma_w。

syms Kse delta_h gamma_s gamma_w D1;
syms d D;
Kse=3;
delta_h=4.8;
gamma_w=10;
gamma_s=9.3;
D1=delta_h;
d=1.0/2*(Kse*delta_h*gamma_w/gamma_s-0.8*D1);
D=ceil(d)
%% https://wenku.baidu.com/view/b869b69c29ea81c758f5f61fb7360b4c2f3f2a67