非稳定渗流场数值计算和闸基渗流“异常现象”研究
|
摘要
本文在总结分析前人研究成果的基础上，论述了用有限元法计算非稳定渗流场的基本理论，并将其应用到实际工程中。文中主要做了以下工作：
(1)通过变分原理，运用插值函数和等参变换推导了渗流微分方程的有限元法格式，建立了非稳定渗流数学模型，讨论了有限元方程求解的稳定性，介绍了本文所采用的求解大型稀疏方程组的基于对称逐步超松弛迭代的预处理共轭梯度法(SSOR-PCG)。
(2)求解无压渗流场问题的关键在于确定渗流逸出面和自由面的位置，在分析目前常用的几种求解自由面方法的基础上，重点介绍了本文采用的求解有自由面渗流问题的固定网格的结点虚流量法，并通过算例说明该方法能实现无压渗流场的快速迭代求解。
(3)非稳定渗流场数值求解中的自由面积分项很难精确处理，本文将曲面积分转换为在坐标平面上的积分来近似处理非稳定渗流的自由面边界，编制了非稳定渗流场有限元计算程序，并通过算例验证了程序的正确性。
(4)遗传算法作为一种高效的全局寻优的方法，经过改进其局部寻优能力，能够作为渗流反分析的重要工具。本文结合非稳定渗流计算程序，运用改进的加速遗传算法编制了渗流问题的反分析计算程序，并对本程序进行了算例验证。
(5)针对嶂山闸闸基的渗流“异常现象”，将可能引起这种“异常现象”的主要因素概化为特征参数变量，经过对闸基主要参数的反分析、渗流场的正分析以及计算结果的对比分析，搞清了该闸闸基渗流“异常现象”的具体成因，最后提出了经济可靠的工程加固处理方法。
Based on the previous research findings in hydraulic engineering, the basic theories on caculating the unsteady seepage field with finited element method (FEM) are studied in this paper. The key contents involved are as follows:
(1)Based on the theory of variation, the FEM form of seepage calculation is concluded by using interpolation function and isoparametric transformation, and the mathematical model of unsteady seepage is established. The solution stability of the FEM matrix is studied, and the symmteric successive over relaxation-precondition conjugate gradient method (SSOR-PCG) to solve the large-scale FEM equations is recommended.
(2)After the several methods to solve this problem have been analyzed, the key issue in solving unsteady seepage problem is to identify the location of free table and outflow line. The nodal virtual flux method is used to solve the seepage problem with free table, and the validity of the method is proved by an example.
(3)The approximate method is proposed by changing the surface integral to plane integral, because it is difficult to correctly deal with the integral about the free table. A program for solving the unsteady seepage field is written, and it is proved to be correct by an example.
(4)As a effective method on seeking the reasonable result in complete region, the genetic algorithm (GA) could be a important tool in seepage conter-analysis, after being improved its partial search ability. A program for the seepage conter-analysis is written, based on the improved accelerating genetic algorithm (IAGA), combined with the program of solving the unsteady seepage field, and it is validated by a example.
(5)The primary factors about the abnormal phenomenon of seepage in Zhangshan Lock are summed up as a series of characteristic parameter variable. The causes which bring the abnormal phenomenon are made clear by conter-analysis, seepage analysis and contrast analysis, and then some economical reliable processing measures are proposed.
(1)通过变分原理，运用插值函数和等参变换推导了渗流微分方程的有限元法格式，建立了非稳定渗流数学模型，讨论了有限元方程求解的稳定性，介绍了本文所采用的求解大型稀疏方程组的基于对称逐步超松弛迭代的预处理共轭梯度法(SSOR-PCG)。
(2)求解无压渗流场问题的关键在于确定渗流逸出面和自由面的位置，在分析目前常用的几种求解自由面方法的基础上，重点介绍了本文采用的求解有自由面渗流问题的固定网格的结点虚流量法，并通过算例说明该方法能实现无压渗流场的快速迭代求解。
(3)非稳定渗流场数值求解中的自由面积分项很难精确处理，本文将曲面积分转换为在坐标平面上的积分来近似处理非稳定渗流的自由面边界，编制了非稳定渗流场有限元计算程序，并通过算例验证了程序的正确性。
(4)遗传算法作为一种高效的全局寻优的方法，经过改进其局部寻优能力，能够作为渗流反分析的重要工具。本文结合非稳定渗流计算程序，运用改进的加速遗传算法编制了渗流问题的反分析计算程序，并对本程序进行了算例验证。
(5)针对嶂山闸闸基的渗流“异常现象”，将可能引起这种“异常现象”的主要因素概化为特征参数变量，经过对闸基主要参数的反分析、渗流场的正分析以及计算结果的对比分析，搞清了该闸闸基渗流“异常现象”的具体成因，最后提出了经济可靠的工程加固处理方法。
Based on the previous research findings in hydraulic engineering, the basic theories on caculating the unsteady seepage field with finited element method (FEM) are studied in this paper. The key contents involved are as follows:
(1)Based on the theory of variation, the FEM form of seepage calculation is concluded by using interpolation function and isoparametric transformation, and the mathematical model of unsteady seepage is established. The solution stability of the FEM matrix is studied, and the symmteric successive over relaxation-precondition conjugate gradient method (SSOR-PCG) to solve the large-scale FEM equations is recommended.
(2)After the several methods to solve this problem have been analyzed, the key issue in solving unsteady seepage problem is to identify the location of free table and outflow line. The nodal virtual flux method is used to solve the seepage problem with free table, and the validity of the method is proved by an example.
(3)The approximate method is proposed by changing the surface integral to plane integral, because it is difficult to correctly deal with the integral about the free table. A program for solving the unsteady seepage field is written, and it is proved to be correct by an example.
(4)As a effective method on seeking the reasonable result in complete region, the genetic algorithm (GA) could be a important tool in seepage conter-analysis, after being improved its partial search ability. A program for the seepage conter-analysis is written, based on the improved accelerating genetic algorithm (IAGA), combined with the program of solving the unsteady seepage field, and it is validated by a example.
(5)The primary factors about the abnormal phenomenon of seepage in Zhangshan Lock are summed up as a series of characteristic parameter variable. The causes which bring the abnormal phenomenon are made clear by conter-analysis, seepage analysis and contrast analysis, and then some economical reliable processing measures are proposed.
引文
[1] 魏茂杰，孙雪琦，陶振常，等．基坑防渗帷幕插入相对不透水层中深度的浅析[J]．防渗技术，2001，7(3)：26-27
[2] 赖远明，吴紫汪，朱元林，等．寒区隧道温度场和渗流场耦合问题的非线性分析[J]．中国科学，1999，(29)：21-26
[3] 赵光恒．力学—现代工程技术的支柱[M]．南京：河海大学出版社，2001，12：34-44
[4] 薛禹群．地下水动力学原理[M]．北京：地质出版社，1986
[5] 毛昶熙．电模拟试验与渗流研究[M]．北京：水利出版社，1981
[6] Neuman S. P., T. N. Narasimhan and P. A. Witherspoon. Application of Mixed Explicit-implicit Finite Element Method to Nonlinear Diffusion-type Problems[M], Finite Elements in Water Resources, Pentech Press, London, Plymouth, 1997, (1): 153～186
[7] Bathe K. J. and Khoshgoftan M. R., Finite Element for Surface Seepage Analysis without Mesh Iteration[J], Int. J. Num. Anal. Methods in Geomechanics, 1979, (3): 3-22
[8] 王勖成，劭敏．有限单元法原理[M]．北京：清华大学出版社，1995
[9] 徐芝纶．弹性力学(上册)[M]．北京：高等教育出版社，1994
[10] 张有天，张武功，半无限域渗流问题的边界元法[J]，水利学报，1981(4)，p8-17
[11] Cundall P. A., Formulation of Three-dimensional Distinct Element Model, Part 1, A Scheme to Detect and Represent Contact in System Composed of Many Polyhedral Blocks[J]. Int. J. Rock Mech. Min. Sci. Geomech. Abstr., 1988, 25(3): 107-116
[12] Shi G. H. & Goodman R. E., Generalization of Two-dimensional Discontinuous Analysis for Forward Modeling[J], Int. J. Num. Anal. Methods Geomech., 1989, (13): 359-380
[13] ITASCA Consulting Group, Inc. In: FLAC 3.3 User's Manual[R], 1996
[14] Pokharel G., Honjo Y., Mapped Infinite Elements in Multi-layered Seepage Analysis[C], Computer Methods and Advances in Geo-mechanics, Balkma, Rotterdam, 1994: 1243-1248
[15] 寇晓东，周维垣．无单元法追踪结构开裂[C]．岩土力学新计算方法讲义，中国科学院武汉岩土力学研究所，1999，10：45-72
[16] 王水林．数值流形方法及其应用[C]．岩土力学新计算方法讲义，中国科学院武汉岩土力学研究所，1999，10：178-207
[17] Zienkiewicz O. C. and Cheung, Y. K. Finite Element in Solution of Field Problems[J]. The Engineering, 1965: 220(5710)
[18] 朱伯芳．有限元原理及其应用[M]．北京．中国水利水电出版社，1998
[19] 林绍忠．对称逐步超松弛预处理共轭梯度法的改进迭代格式[J]．数值计算与计算机应用，1997．12：
[20] 程昭，陈胜宏，田肚利，等．大规模p型有限元方程组的修正SSOR-PCG解法[J]．上海交通大学学报，2002．11，36(11)：1608-1611
[21] 金忠青，陈夕庆．用脉冲谱法确定承压水含水层非均质导水系数[J]．河海大学学报．1991，19(3)：53-63
[22] Sagar B., Yabowitz S., Duuckstein L., A Direct Method for the Identification of the Parameter of Dynamic Non-homogeneous Aquifers[J]. Water Resour. Res., 1975, 11 (4): 563-570
[23] Philip J. R., General Method of Exact Solution of the Concentration Dependent Diffusion Equation[J]. Aust.J.Phys., 1960, 13:1-12
[24] 张乾飞，宋一明，吴耿同，等．南水大坝渗透系数反演分析[J]．水利技术监督，2000，8(4)：41-44
[25] 李守巨，刘迎曦，周承芳，等．岩土渗流反分析的非线性最小二乘法[J]．中国有色金属学报，1998，8(增2)：726-728
[26] Neuman S.P., A Statistical Approach to the Inverse Problem of Aquifer Hydrology, 3. Improved Solution Method and Added Perspective[J]. Water Resource Research, 1980,16(2): 331-346
[27] Cheng J.M.,Yeh W. W-G., A Proposed Quasi Newton Method For Parameter Identification in a Flow and Transport System[J]. Advances in Water Resource Research, 1992, 15(4): 239-250
[28] 朱岳明．测压管水位反分析确定性混合模型[C]．大坝安全监测技术国际学术讨论会论文集，1992：150．154
[29] 吴中如，顾冲时．大坝原型反分析及其应用[M]．南京：江苏科学技术出版社，2000
[30] 周志芳．裂隙地下水反分析中若干问题的探析[J]．水文地质工程地质，1993(4)：33-35
[31] 周志芳，朱学愚．确定岩体渗透数的结构面控制反演法[J]．南京大学学报，1999，35(3)
[32] 周志芳，朱学愚．岩体渗透系数张量的半解析计算[J]．水利学报，1997(9)：6-11
[33] 金忠青，周志芳．工程水力学反问题[M]．南京：河海大学出版社，1997
[34] 张奇，夏颂佑，俞国青，等．裂隙基岩渗透张量反分析及等效连续介质模型[J]．河海大学学报，1994，22(3)：74-80
[35] 李晓峰，朱本仁．Randon变换在渗流力学反问题中的应用[J]．计算物理，2001，18(2)：115-118
[36] 郑东健，顾冲时，张涛，等．测压管水位的混合模型研究[J]．河海大学学报，2001，29(6)：61-64
[37] 刘迎曦，李守巨，李正国，等．岩体渗透系数反演的数值方法及其适定性[J]．辽宁工程技术大学学报(自然科学版)，19(4)：375-378
[38] 刘杰，王媛，刘宁，等．岩土工程渗流参数反问题[J]．岩土力学，2002，23(2)：152-161
[39] 张乾飞，王健，吴中如，等．基于人工神经网络的大坝渗透系数分区反演分析[J]．水电能源科学，2001，19(4)：4-7
[40] 张士儒，张越．三河闸闸基渗流观测的初步分析[J]．水利工程管理技术，1990(2)：47-52
[41] 吴明祥．山河闸闸基渗流规律的探索[J]．中国水利，2003．2，(B)：119-120
[42] 吕世德，陈浩然．嶂山闸闸基渗流异常现象初步分析[J]．水电自动化与大坝监测，2004．4，28(2)：57-59
[43] 唐劲松，刘光标．嶂山闸扬压力异常原因分析[J]．水电自动化与大坝监测，2004．12，28(6)：49-52
[43] 吕世德，陈浩然．嶂山闸闸基渗流异常原因[J]．治淮，2004．2：25-27
[44] 毛昶熙．土石坝非稳定渗流和坝坡稳定分析的有限单元法计算[M]．南京水利科学研究所，1979
[45] 李定方，李祖贻，陈平，等．土石坝三向渗流计算[J]．水利水运科学研究，1980(3)：1～15．
[46] 叶兴才，丁家平．澄西船坞井点降水电模拟试验和有限单元法计算[M]．南京水利科学研 院．1979
[47] 谢春红，陈平．有限单元法在不稳定渗流分析中的初步应用．水利水运科技情报，1975(2)
[48] 毛昶熙．渗流计算分析与控制[M]．北京：中国水利电力出版社，1990
[49] 张蔚榛．地下水非稳定计算与地下水资源评价[M]．北京：科学出版社，1983
[50] 吕涛，石济民，林振宝，等．区域分解法—偏微分方程数值解法新技术[M]．北京：科学出版社．1992
[51] 杜延龄，许国安．渗流分析的有限元法和电网络法[M]．北京：水利电力出版社，1992：1-292
[52] 朱岳明，张燎军．渗流场求解的改进排水子结构法．岩土工程学报[J]．1997．
[53] 朱岳明，龚道勇．三维饱和-非饱和渗流场求解及其逸出面边界条件处理[J]．水科学进展，2003(1)
[54] 李春华．固定网络法用于三向稳定渗流有限元的实践．北京：水利电力出版社，第四界全国水利水电工程渗流学术讨论会论文集，1993
[55] 河野伊一郎．用有限元法解坝体的渗流问题[J]．水利水运科技情报(渗流译文汇编(第9辑))，1978：12-20
[56] 陈洪凯．渗流自由面求解的基本方法、修正及应用[J]．重庆交通学院学报，1997．9：
[57] 张有天，陈平，王雷，等．有自由面渗流分析的初流量法[J]．水利学报，1988(8)：18-26
[58] 速宝玉，朱岳明．不变网格法确定渗流自由面的节点虚流量法[J]．河海大学学报，1991，19(5)：113-117
[59] 吴梦喜，张学勤．有自由面渗流分析的虚单元法[J]．水利学报，1994(8)：67-70
[60] 朱军，刘光廷．改进的单元渗透矩阵调整法求解无压渗流场[J]．水利学报，2001．8(8)：49-52
[61] 速宝玉，沈振中，赵坚．用变分不等式理论求解渗流问题的截止负压法[J]．水利学报，1996(3)：22-29
[62] 张乾飞．复杂渗流场演变规律及转异特征研究[D]．南京：河海大学博士学位论文，2002，9：20-30
[63] 王媛．求解有自由面渗流问题的初流量法的改进[J]．水利学报，1998(3)：68-73．
[64] Huang S., Zhou C., Numerical solution of non-steady state porous flows free boundary problems[J]. J.of computational math., 1985,3(1):72-89.
[65] 佘颍禾，孙鹰，郭小明．具有自由边界的二维渗流问题[J]．应用数学和力学，1996，17(6)：523-527
[66] 孙鹰，郭小明，佘颍禾．非稳态渗流的自由边界问题[J]．应用数学和力学，1999，20(7)：756-760
[67] 熊文林，周创兵．有自由面渗流分析的子单元法[J]．水利学报，1997(8)：34-38
[68] 陈洪凯，唐红梅，肖盛燮．求解渗流自由面的复合单元全域迭代法[J]．应用数学和力学，1999，20(10)：1045-1058
[69] 黄蔚，刘迎曦，周承芳．三维无压渗流场的有限元算法研究[J]．水利学报，2001(6)：33-36
[70] 凌道盛．有自由面渗流分析的虚节点法[J]．浙江大学学报(工学版)，2002，36(3)：243-246
[71] 郭洪兴．三维非稳定渗流计算中有限元截止负压法的研究及应朋[D]．南京：河海大学硕士学位论文1997．3
[72] 刘杰．裂隙岩体渗流场及其与应力场耦合的参数反问题研究[D]．南京：河海大学博士学位论文，2002．7
[73] 徐士良．FORTRAN常用算法程序集(第二版)．北京：清华大学出版社，1996
[74] 赤井浩一．用有限单元法分析饱和非饱和渗流问题[A]．日本土木学会论文报告集[C]．1977-1978
[75] 周明，孙树栋．遗传算法原理及应用[M]．北京：国防工业出版社，1999．6
[76] 金菊良，杨晓华，储开凤，等．加速基因算法在海洋环境预报中的应用[J]．海洋环境科学，1997(4)：7-11
[77] 陈建余，朱岳明，陈晓明，等．改进加速遗传算法及其在非稳定渗流场反分析中的应用[J]．水电能源科学，2003．9，21(3)：59-61
[78] 刘杰，王瑗．一种高效混合遗传算法[J]．2002(2)：49-53
[79] 陈建余．非稳定饱和-非饱和渗流场数值计算关键技术及其应用研究[D]．河海大学博士论文，2003．10
[80] 付成华．非均质土石坝稳定_非稳定渗流有限元分析[D]．武汉大学硕士论文，2004．5
[81] 裴利华．有自由面三维非稳定渗流数值模拟若干关键技术研究[D]．河海大学硕士论文．2005．4
[2] 赖远明，吴紫汪，朱元林，等．寒区隧道温度场和渗流场耦合问题的非线性分析[J]．中国科学，1999，(29)：21-26
[3] 赵光恒．力学—现代工程技术的支柱[M]．南京：河海大学出版社，2001，12：34-44
[4] 薛禹群．地下水动力学原理[M]．北京：地质出版社，1986
[5] 毛昶熙．电模拟试验与渗流研究[M]．北京：水利出版社，1981
[6] Neuman S. P., T. N. Narasimhan and P. A. Witherspoon. Application of Mixed Explicit-implicit Finite Element Method to Nonlinear Diffusion-type Problems[M], Finite Elements in Water Resources, Pentech Press, London, Plymouth, 1997, (1): 153～186
[7] Bathe K. J. and Khoshgoftan M. R., Finite Element for Surface Seepage Analysis without Mesh Iteration[J], Int. J. Num. Anal. Methods in Geomechanics, 1979, (3): 3-22
[8] 王勖成，劭敏．有限单元法原理[M]．北京：清华大学出版社，1995
[9] 徐芝纶．弹性力学(上册)[M]．北京：高等教育出版社，1994
[10] 张有天，张武功，半无限域渗流问题的边界元法[J]，水利学报，1981(4)，p8-17
[11] Cundall P. A., Formulation of Three-dimensional Distinct Element Model, Part 1, A Scheme to Detect and Represent Contact in System Composed of Many Polyhedral Blocks[J]. Int. J. Rock Mech. Min. Sci. Geomech. Abstr., 1988, 25(3): 107-116
[12] Shi G. H. & Goodman R. E., Generalization of Two-dimensional Discontinuous Analysis for Forward Modeling[J], Int. J. Num. Anal. Methods Geomech., 1989, (13): 359-380
[13] ITASCA Consulting Group, Inc. In: FLAC 3.3 User's Manual[R], 1996
[14] Pokharel G., Honjo Y., Mapped Infinite Elements in Multi-layered Seepage Analysis[C], Computer Methods and Advances in Geo-mechanics, Balkma, Rotterdam, 1994: 1243-1248
[15] 寇晓东，周维垣．无单元法追踪结构开裂[C]．岩土力学新计算方法讲义，中国科学院武汉岩土力学研究所，1999，10：45-72
[16] 王水林．数值流形方法及其应用[C]．岩土力学新计算方法讲义，中国科学院武汉岩土力学研究所，1999，10：178-207
[17] Zienkiewicz O. C. and Cheung, Y. K. Finite Element in Solution of Field Problems[J]. The Engineering, 1965: 220(5710)
[18] 朱伯芳．有限元原理及其应用[M]．北京．中国水利水电出版社，1998
[19] 林绍忠．对称逐步超松弛预处理共轭梯度法的改进迭代格式[J]．数值计算与计算机应用，1997．12：
[20] 程昭，陈胜宏，田肚利，等．大规模p型有限元方程组的修正SSOR-PCG解法[J]．上海交通大学学报，2002．11，36(11)：1608-1611
[21] 金忠青，陈夕庆．用脉冲谱法确定承压水含水层非均质导水系数[J]．河海大学学报．1991，19(3)：53-63
[22] Sagar B., Yabowitz S., Duuckstein L., A Direct Method for the Identification of the Parameter of Dynamic Non-homogeneous Aquifers[J]. Water Resour. Res., 1975, 11 (4): 563-570
[23] Philip J. R., General Method of Exact Solution of the Concentration Dependent Diffusion Equation[J]. Aust.J.Phys., 1960, 13:1-12
[24] 张乾飞，宋一明，吴耿同，等．南水大坝渗透系数反演分析[J]．水利技术监督，2000，8(4)：41-44
[25] 李守巨，刘迎曦，周承芳，等．岩土渗流反分析的非线性最小二乘法[J]．中国有色金属学报，1998，8(增2)：726-728
[26] Neuman S.P., A Statistical Approach to the Inverse Problem of Aquifer Hydrology, 3. Improved Solution Method and Added Perspective[J]. Water Resource Research, 1980,16(2): 331-346
[27] Cheng J.M.,Yeh W. W-G., A Proposed Quasi Newton Method For Parameter Identification in a Flow and Transport System[J]. Advances in Water Resource Research, 1992, 15(4): 239-250
[28] 朱岳明．测压管水位反分析确定性混合模型[C]．大坝安全监测技术国际学术讨论会论文集，1992：150．154
[29] 吴中如，顾冲时．大坝原型反分析及其应用[M]．南京：江苏科学技术出版社，2000
[30] 周志芳．裂隙地下水反分析中若干问题的探析[J]．水文地质工程地质，1993(4)：33-35
[31] 周志芳，朱学愚．确定岩体渗透数的结构面控制反演法[J]．南京大学学报，1999，35(3)
[32] 周志芳，朱学愚．岩体渗透系数张量的半解析计算[J]．水利学报，1997(9)：6-11
[33] 金忠青，周志芳．工程水力学反问题[M]．南京：河海大学出版社，1997
[34] 张奇，夏颂佑，俞国青，等．裂隙基岩渗透张量反分析及等效连续介质模型[J]．河海大学学报，1994，22(3)：74-80
[35] 李晓峰，朱本仁．Randon变换在渗流力学反问题中的应用[J]．计算物理，2001，18(2)：115-118
[36] 郑东健，顾冲时，张涛，等．测压管水位的混合模型研究[J]．河海大学学报，2001，29(6)：61-64
[37] 刘迎曦，李守巨，李正国，等．岩体渗透系数反演的数值方法及其适定性[J]．辽宁工程技术大学学报(自然科学版)，19(4)：375-378
[38] 刘杰，王媛，刘宁，等．岩土工程渗流参数反问题[J]．岩土力学，2002，23(2)：152-161
[39] 张乾飞，王健，吴中如，等．基于人工神经网络的大坝渗透系数分区反演分析[J]．水电能源科学，2001，19(4)：4-7
[40] 张士儒，张越．三河闸闸基渗流观测的初步分析[J]．水利工程管理技术，1990(2)：47-52
[41] 吴明祥．山河闸闸基渗流规律的探索[J]．中国水利，2003．2，(B)：119-120
[42] 吕世德，陈浩然．嶂山闸闸基渗流异常现象初步分析[J]．水电自动化与大坝监测，2004．4，28(2)：57-59
[43] 唐劲松，刘光标．嶂山闸扬压力异常原因分析[J]．水电自动化与大坝监测，2004．12，28(6)：49-52
[43] 吕世德，陈浩然．嶂山闸闸基渗流异常原因[J]．治淮，2004．2：25-27
[44] 毛昶熙．土石坝非稳定渗流和坝坡稳定分析的有限单元法计算[M]．南京水利科学研究所，1979
[45] 李定方，李祖贻，陈平，等．土石坝三向渗流计算[J]．水利水运科学研究，1980(3)：1～15．
[46] 叶兴才，丁家平．澄西船坞井点降水电模拟试验和有限单元法计算[M]．南京水利科学研 院．1979
[47] 谢春红，陈平．有限单元法在不稳定渗流分析中的初步应用．水利水运科技情报，1975(2)
[48] 毛昶熙．渗流计算分析与控制[M]．北京：中国水利电力出版社，1990
[49] 张蔚榛．地下水非稳定计算与地下水资源评价[M]．北京：科学出版社，1983
[50] 吕涛，石济民，林振宝，等．区域分解法—偏微分方程数值解法新技术[M]．北京：科学出版社．1992
[51] 杜延龄，许国安．渗流分析的有限元法和电网络法[M]．北京：水利电力出版社，1992：1-292
[52] 朱岳明，张燎军．渗流场求解的改进排水子结构法．岩土工程学报[J]．1997．
[53] 朱岳明，龚道勇．三维饱和-非饱和渗流场求解及其逸出面边界条件处理[J]．水科学进展，2003(1)
[54] 李春华．固定网络法用于三向稳定渗流有限元的实践．北京：水利电力出版社，第四界全国水利水电工程渗流学术讨论会论文集，1993
[55] 河野伊一郎．用有限元法解坝体的渗流问题[J]．水利水运科技情报(渗流译文汇编(第9辑))，1978：12-20
[56] 陈洪凯．渗流自由面求解的基本方法、修正及应用[J]．重庆交通学院学报，1997．9：
[57] 张有天，陈平，王雷，等．有自由面渗流分析的初流量法[J]．水利学报，1988(8)：18-26
[58] 速宝玉，朱岳明．不变网格法确定渗流自由面的节点虚流量法[J]．河海大学学报，1991，19(5)：113-117
[59] 吴梦喜，张学勤．有自由面渗流分析的虚单元法[J]．水利学报，1994(8)：67-70
[60] 朱军，刘光廷．改进的单元渗透矩阵调整法求解无压渗流场[J]．水利学报，2001．8(8)：49-52
[61] 速宝玉，沈振中，赵坚．用变分不等式理论求解渗流问题的截止负压法[J]．水利学报，1996(3)：22-29
[62] 张乾飞．复杂渗流场演变规律及转异特征研究[D]．南京：河海大学博士学位论文，2002，9：20-30
[63] 王媛．求解有自由面渗流问题的初流量法的改进[J]．水利学报，1998(3)：68-73．
[64] Huang S., Zhou C., Numerical solution of non-steady state porous flows free boundary problems[J]. J.of computational math., 1985,3(1):72-89.
[65] 佘颍禾，孙鹰，郭小明．具有自由边界的二维渗流问题[J]．应用数学和力学，1996，17(6)：523-527
[66] 孙鹰，郭小明，佘颍禾．非稳态渗流的自由边界问题[J]．应用数学和力学，1999，20(7)：756-760
[67] 熊文林，周创兵．有自由面渗流分析的子单元法[J]．水利学报，1997(8)：34-38
[68] 陈洪凯，唐红梅，肖盛燮．求解渗流自由面的复合单元全域迭代法[J]．应用数学和力学，1999，20(10)：1045-1058
[69] 黄蔚，刘迎曦，周承芳．三维无压渗流场的有限元算法研究[J]．水利学报，2001(6)：33-36
[70] 凌道盛．有自由面渗流分析的虚节点法[J]．浙江大学学报(工学版)，2002，36(3)：243-246
[71] 郭洪兴．三维非稳定渗流计算中有限元截止负压法的研究及应朋[D]．南京：河海大学硕士学位论文1997．3
[72] 刘杰．裂隙岩体渗流场及其与应力场耦合的参数反问题研究[D]．南京：河海大学博士学位论文，2002．7
[73] 徐士良．FORTRAN常用算法程序集(第二版)．北京：清华大学出版社，1996
[74] 赤井浩一．用有限单元法分析饱和非饱和渗流问题[A]．日本土木学会论文报告集[C]．1977-1978
[75] 周明，孙树栋．遗传算法原理及应用[M]．北京：国防工业出版社，1999．6
[76] 金菊良，杨晓华，储开凤，等．加速基因算法在海洋环境预报中的应用[J]．海洋环境科学，1997(4)：7-11
[77] 陈建余，朱岳明，陈晓明，等．改进加速遗传算法及其在非稳定渗流场反分析中的应用[J]．水电能源科学，2003．9，21(3)：59-61
[78] 刘杰，王瑗．一种高效混合遗传算法[J]．2002(2)：49-53
[79] 陈建余．非稳定饱和-非饱和渗流场数值计算关键技术及其应用研究[D]．河海大学博士论文，2003．10
[80] 付成华．非均质土石坝稳定_非稳定渗流有限元分析[D]．武汉大学硕士论文，2004．5
[81] 裴利华．有自由面三维非稳定渗流数值模拟若干关键技术研究[D]．河海大学硕士论文．2005．4