0

GEO5 software

Online nápověda

Stromeček
Nastavení
Program:
Jazyk:

渗流分析

非稳定流

非饱和介质中的非稳定流问题可以由Richard方程(连续性方程)得到:

其中:

n

-

材料孔隙率

-

饱和度变化率

Kr

-

相对渗透系数

-

饱和介质渗透系数矩阵

-

总水头梯度

Richard方程的时间离散采用的是显式Picard迭代法[1]。这样,就能构造一个可以保证质量守恒的杂交元列式。由于非稳定渗流问题是非线性的,所以分析过程往往是递增的。同时,采用了标准的Newton-Raphson迭代法,以满足平衡条件。

应当注意的是,迭代过程的速度和稳定性在很大程度上受到所选材料模型的影响,换句话说,也就是受到计算相对渗透系数Kr、饱和度S的方法和相关参数取值的影响。例如,对于典型的砂性土,当设置的初始条件不恰当时,可能导致数值分析问题。更多细节和信息请参考文献[2, 3]。

稳定流

稳定流分析假设岩土体的饱和度随时间是不变的。相应的计算公式如下:

和非稳定渗流分析不同,稳定渗流分析和时间无关,仅仅需要输入渗流边界条件。但是,这仍是一个一般意义上的非线性问题(例如,非承压流分析),需要采用Newton-Raphson法进行迭代计算。更多细节和信息请参考文献[2, 3]。

参考文献:

[1] M. A. Celia and E. T. Bouloutas, A general mass-conservative numerical solutionfor the unsaturated flow equation, Water Resources Research 26 (1990), no. 7, 1483-1496.

[2] M. Šejnoha, Finite element analysis in geotechnical design, to appear (2015).

[3] M. Šejnoha, T. Janda, H. Pruška, M. Brouček, Metoda konečných prvků v geomechanice: Teoretické základy a inženýrské aplikace, předpokládaný rok vydání (2015).

Vyzkoušejte si programy GEO5. Zdarma.