我现在正在用ＡＮＳＹＳ的ＡＰＤＬ编非稳定渗流的计算程序，遇到了不少问题饱和－非饱和非稳定渗流时渗流微分控制方程式右端项C(h)＝d θ/d hc以及时段的选取对结果的影响很大，我编的程序与赤井浩一的砂槽模型结果相差很大，水位上升的结果还可以，水位下降的结果相差太多请大家给出出主意．欢迎大家在此发表ＡＮＳＹＳ与渗流方面的观点与意见
我现在正在用ＡＮＳＹＳ的ＡＰＤＬ编非稳定渗流的计算程序，遇到了不少问题
饱和－非饱和非稳定渗流时渗流微分控制方程式右端项C(h)＝d θ/d hc以及时段的选取对结果的影响很大，我编的程序与赤井浩一的砂槽模型结果相差很大，水位上升的结果还可以，水位下降的结果相差太多
请大家给出出主意．
欢迎大家在此发表ＡＮＳＹＳ与渗流方面的观点与意见
2楼
堤坝饱和_非饱和渗流的数值分析---高冀
回复
3楼
这个我一直就没搞清楚。
回复
4楼
渗流控制方程式是高度非线性的，容水度与渗透系数在非饱和区均为负压的非线性函数，并且对于土石坝还存在着浸润线与溢出面边界位置事先不确定，必需根据每一次计算结果重新调整不断迭代才能接近真实值，两个不确定因素结合在一起，如何在编程时将两种循环迭代是一个难点．所取时间步过长而计算精度与参数选取精度不够，步长过小又很导致过大的累积误差，模拟计算不能收敛于实验结果的原因也可能是这个原因造成的．
对于土石坝上游水位骤升时由于下游溢出点即使在下游存在，也不会离水面太高，这时的模拟结果已经与实验结果基本吻合；上游水位骤降时，溢出点的位置则对模拟结果影响非常大，我现在模拟出的结果很不理想．
回复
5楼
作研究，关键问题在于溢出点的确定问题，渗透参数与溢出点的交叉迭代可能最后收敛不稳定，还有就是对初始步长的选择难以确定．你能不能对别的算法给出些建议？
回复
6楼
我怎么下不了呀，急！
回复
7楼
谢谢！可以下载！
回复
8楼
渗流与边坡稳定分析，群号39170907
有兴趣的朋友加入讨论问题，用容重替换作为考虑渗流的影响进行边坡稳定分析，不论是极限平衡法还是强度折减有限元，所得最小安全系数，均是在水位降至坝高1/3左右时，但采用有限元法计算，用直接施加节点力的方法作为考虑渗透力影响，却得到最小边坡稳定安全系数在水位降至最低点时，有做过类似研究的吗？
回复
