初始状态地下水位以上空隙水压力定义问题

beiximan

在SEEP中模拟膨胀土降雨入渗问题时
计算区域高度取10米，宽度11米，土性函数取Silt(Tailings)-Function#22,地下水位距离地表10米，设地表总水头值10.1米，两侧边不透水，模拟地表径流情况，请问如何实现初始状态地表吸力值为30kpa，且随深度呈线性减小？

basara

参考算例Kisch liners and caps设置初始水压力场

basara

自己做了一个小例子，得到了比较漂亮的结果。
总结模拟的两个关键点：1.初始水力场的设定；2.模拟地表降雨边界的瞬态分析需要在地表表层设定一个surface layers，这个是必须的，否则无法得到收敛的结果的。
下图是设定的满足楼主要求的初始水力场
这个水力场是通过keyin 菜单的 Spatial Functions 设定，定义高程0的两个点压力水头为0，定义高程为10的两个点压力水头为-3，其它点水压力软件自动线性差值得到。
我定义了降水边界为Unit Flux (q) 值3e-5m/sec。分析时间设定为1天
以下分别是1小时、5.5、10、19.5小时、1天的水压力云图

对此问题感兴趣的朋友可以与我交流。

beiximan

谢谢basara,我好好琢磨下，对于水位以上合理的吸力分布问题，不知basara及各位同仁有何看法？

beiximan

我用的是2004版，是不是没有“keyin 菜单的 Spatial Functions ”？？:)

basara

关于水位以上的吸力分布跟非饱和区土的体积含水量密切相关。基质吸力与体积含水量是个对应的关系，这个关系由试验得出，所以精确的定义应该由试验得出。
如果04版本菜单中没有Spatial Functions 估计是07版本的新增的功能。