python自动化算法_基于Python语言和Abaqus平台的边坡可靠度计算自动化算法开发

2.2 求解过程

求解过程包括7步，如

图 2

图 2

自动化程序的求解过程

Fig. 2

The solving part of the automation program

1) Python形成初始Inp文件。将边坡的几何参数、材料信息及随机场数据写入Inp文件，该文件称为初始Inp文件。

2) Inp文件进行初次运算。该步骤的主要目的是平衡地应力。在初始Inp文件中，有施加土体重力的分析步骤。边坡在初始状态下，由于自重作用，存在与重力相平衡的应力状态，因此，在进行数值模拟时，需要在边坡开始运算之前建立相应的应力场。

3) 提取初始应力生成Rpt文件。经过步骤2)的初始运算，得到一系列包含各个单元应力的Job文件，然后使用Python编写的脚本文件，提取各个单元的内力，并生成包含各个应力提取代码的Rpt文件。

4) 地应力平衡的Csv文件。在Abaqus平台上运行步骤3)得到Rpt文件，得到与每种情况相对应的Csv文件，以便平衡地应力。

5) 地应力平衡。模型的地应力平衡结果满足要求后，程序自动调用提前编写的命令读取Csv文件。

6) 得到最终的Inp文件。在初始Inp文件中加入强度折减法的分析步，得到最终的Inp文件。

7) Abaqus强度折减法运算。调用Abaqus求解器得到最终包含边坡变形、应力和场变量等信息的Job文件。

8) 计算边坡的失效概率。根据Pf=Nfs < 1/N(Pf表示边坡的失效概率; Nfs < 1表示安全系数小于1的数量; N表示总的计算次数)输出边坡的失效概率。

3 算例验证

为验证编写的自动化算法程序的精度，采用经典边坡算例。边坡尺寸如ρ为2 000 kg/m3，变形模量为10 MPa，泊松比v=0.3[。为简化计算，只考虑黏聚力生成的随机场，内摩擦角为0°。

图 3

图 3

典型边坡几何尺寸

Fig. 3

Size of slope

边坡采用平面应变单元CPE4，共划分910个单元，971个单元节点，土体失效模式采用Mohr-Coulomb屈服准则。边界条件为约束边界的侧向位移及底部的水平及竖向位移。Der Kiureghian等[和Huang等[指出单元尺寸与相关距离之比应小于0.25。单元水平长度为2 m，高度为0.5 m，其中，单元水平长度/水平相关距离=2/38=0.05 < 0.25，单元高度/竖向相关距离=0.5/3.8=0.13 < 0.25，单元尺寸符合要求。

地应力平衡是岩土工程数值模拟过程中