1.一种边坡强度折减系数多线程并行计算方法,其特征在于,包括以下步骤:S1、根据边坡岩土体性质确定边坡强度的折减系数上限值Kmax和折减系数下限值Kmin;
S2、通过线性插值形成折减系数等差序列(Kmin,K1,K2,...,KT,Kmax);T为并行计算线程数,T为正整数;
S3、将FLAC3D命令流内嵌于Python脚本中,通过Python脚本同时调用T个FLAC3D应用程序,分别对折减系数为K1,K2,...,KT的数值模拟计算;
S4、根据数值模拟计算结果确定折减系数收敛区域;
‑5
所述折减系数收敛的准则为不平衡力比率小于或等于10 ,所述不平衡力比率为节点平均内力与最大不平衡力的比值;
S5、根据折减系数收敛区域分别对折减系数上限值Kmax和折减系数下限值Kmin进行更新,得到更新后的折减系数上限值K′max和折减系数下限值K′min;
S6、判断是否满足K′max‑K′min<η,若是则进入步骤S7,否则返回步骤S2,将Kmax和Kmin更新为K′max和K′min后进行迭代;η为给定计算误差;
S7、根据更新后的折减系数上限值K′max和折减系数下限值K′min计算得到边坡强度折减系数K。
2.根据权利要求1所述的边坡强度折减系数多线程并行计算方法,其特征在于,所述步骤S1中折减系数上限值Kmax初始设置为20,折减系数下限值Kmin初始设置为0。
3.根据权利要求1所述的边坡强度折减系数多线程并行计算方法,其特征在于,所述步骤S2中并行计算线程数T=N‑1,N为计算机处理器的线程数。
4.根据权利要求1所述的边坡强度折减系数多线程并行计算方法,其特征在于,所述步骤S5具体为:
若K1不收敛,则K′min=Kmin,K′max=K1;
若Kt‑1收敛,Kt不收敛,则K′min=Kt‑1,K′max=Kt,t=2,3,...,T;
若KT收敛,则K′min=KT,K′max=Kmax。
5.根据权利要求1所述的边坡强度折减系数多线程并行计算方法,其特征在于,所述步骤S6中计算误差η根据边坡安全系数的设计值设置为0.001。
6.根据权利要求1所述的边坡强度折减系数多线程并行计算方法,其特征在于,所述步骤S7中边坡强度折减系数K的计算公式为: