非饱和土渗流变形耦合的数值分析-重庆大学学报.PDF

第 卷第 期 土 木 建 筑 与 环 境 工 程 ３３ ３ Ｖｏｌ．３３Ｎｏ．３ 年 月 ， ２０１１ ６ ＪｏｕｒｎａｌｏｆＣｉｖｉｌＡｒｃｈｉｔｅｃｔｕｒａｌ＆ＥｎｖｉｒｏｎｍｅｎｔａｌＥｎｉｎｅｅｒｉｎ Ｊｕｎ．２０１１ ｇ ｇ 非饱和土渗流 变形耦合的数值分析  吴礼舟，黄润秋 （成都理工大学 地质灾害防治与地质环境保护国家重点实验室，成都 ６１００５９） 摘 要：基于一维非饱和土的渗流 变形控制方程，采用 （ ）软  ＦｌｅｘＰＤＥＰａｒｔｉａｌｄｉｆｆｅｒｅｎｔｉａｌｅｕａｔｉｏｎ 　 ｑ 件对该耦合方程组进行求解分析。该方法突破了解析法对非饱和土导水系数函数的特殊限定，适 用于任意的土 水特征曲线表达式；还可考虑到饱和时的渗透系数以及孔隙率是变量。与解析解相  比，该数值解表现较高的精度，具有解决非饱和土耦合问题的可行性。计算分析表明，非饱和土渗 流 变形耦合作用对暂态孔隙水压力分布产生重要的影响，在降雨入渗过程中需考虑土体渗流 变形   耦合的影响。降雨初期，位移随着时间明显增大，地表出现下沉，考虑耦合效应的孔隙水压力慢于 非耦合情况，原因是 值为正的。随着降雨持续时间的增大，地表下沉的速度减缓，到最后变形开 犎 始稳定。位移的变化快慢与孔隙水压力变化规律相同。地表沉降量还与初始孔隙水压力分布以及 犎值密切相关。饱和时的渗透系数以及孔隙率对非饱和土降雨入渗以及稳态流的分布产生影响， 但对其地表变形产生的影响微弱。 关键词：非饱和土；渗流 变形耦合；降雨入渗；数值分析  中图分类号： 文献标志码： 文章编号： （ ） Ａ １６７４４７６４２０１１０３００６３０５ 　　 　　 犖狌犿犲狉犻犮犪犾犃狀犪犾狊犻狊狅犳犛犲犲犪犲犪狀犱犇犲犳狅狉犿犪狋犻狅狀犻狀犝狀狊犪狋狌狉犪狋犲犱犛狅犻犾狊 狔 狆犵 ， 犠犝犔犻狕犺狅狌犎犝犃犖犌犚狌狀犻狌 狇 （ ， ＳｔａｔｅＫｅ ＬａｂｏｒａｔｏｒｏｆＧｅｏｌｏｉｃａｌＨａｚａｒｄＰｒｅｖｅｎｔｉｏｎａｎｄＧｅｏｌｏｉｃａｌＥｎｖｉｒｏｎｍｅｎｔＰｒｏｔｅｃｔｉｏｎ ｙ ｙ ｇ ｇ ， ，