图9159m定水头设置 175m171m167m163m159m 图10三舟溪滑坡工况1边界条件设置情况
后续在相应位置assign边界条件,注意一定为相应高程处。如果相应高程不能对应,可通过draw points按钮在175m、159m等处添加。然后继续添加边界条件。
单位设置可通过set units and scale 完成,见图11,时间单位一般选择days。
图10 单位的设置
初始地下水位的添加,首先keyin analyses中,选择water table选项,见图11。
图11初始地下水位的添加设置
按照工况1,时间步为123days,设置见图12。一定要设置时间步,否则无法计算完全。
图12时间步的设置
初始地下水位首先是175m出水点,大致平行于地形线,形状一般呈凸形,同时参考统计的水力坡度角大小。
上述参数均设置好后,可通过verify进行模型检查,若无错误,点击solve按钮start求解,求解过程需收敛,convergence tolerance 为0.1,查看方式为solve求解过程,不要勾选close solve after each analysis选项,计算完成后点击graph,查看是否收敛,见图13、14,否则调整水力学参数或边界条件。同时,计算过程中若水位超出地面线,在draw boundary condition选项,点击keyin按钮设置边界条件参数时,勾选Potential seepage face review 选项,超出地面的水沿地表径流;若水位未超出地面线,不要勾选此选项,见图15。
图13计算求解
图14查看是否收敛
图15potential seepage face review
上述计算完成后,结果查看时,在contour状态下draw isolines利用ctrl键可显示不同天地下水浸润线,见图16,draw vectors可显示地下水流向,graph可设置显示一竖剖面的水头或孔隙水压力随时间的变化曲线及数据。。
图16 地下水浸润线显示方法
以三舟溪滑坡为例,选择第0、31、61、92、123天,可显示175、171、167、163、159m时地下水浸润线分布,如下图17、18所示。此过程中,浸润线分布应合理,随着库水位降,浸润线下降,否则调整计算。上述渗流分析到此完成。
250240230220210250公路240230220210Elevation2001901801701601501401300100STK-2STK1WZ03200190175m180171m167m163m159m170160150140200300130400Distance
图16 三舟溪滑坡浸润线分布特征