含节理黄土渗流数值模型研究-工程地质学报.pdfVIP

－ ＪｏｕｒｎａｌｏｆＥｎｇｉｎｅｅｒｉｎｇＧｅｏｌｏｇｙ　　工程地质学报　２０１４／２２（６）１１１５０８ ＤＯＩ：１０．１３５４４／ｊ．ｃｎｋｉ．ｊｅｇ．２０１４．０６．０１５ 含节理黄土渗流数值模型研究 罗　扬　王铁行　王娟娟 （西安建筑科技大学土木工程学院　西安　７１００５５） 摘　要　本文建立了在考虑黄土中有节理存在时，水分场数值计算的有限元方法。采用质量守恒的观点推导了质点元中饱 和度的变化与流速的关系。进而利用达西定律得到以水头为变量的渗流基本方程。针对黄土垂直节理的渗流特点，确定了 节理渗流基本方程中的参数。采用四边形等参元，利用Ｇａｌｅｒｋｉｎ加权余量法建立考虑节理影响的黄土非饱和渗流的有限元形 式。对局部水头边界条件下的黄土节理二维水分入渗问题进行了数值分析。结果表明，节理对黄土场地湿润峰的迁移有很 大影响。 关键词　黄土节理　非饱和　Ｇａｌｅｒｋｉｎ法　有限元　数值分析 ＋ 中图分类号：Ｐ６４２１３１　　文献标识码：Ａ ＦＩＮＩＴＥＥＬＥＭＥＮＴＳＥＥＰＡＧＥＦＬＯＷ ＭＯＤＥＬＦＯＲＵＮＳＡＴＵＲＡＴＥＤＬＯＥＳＳ ＷＩＴＨＪＯＩＮＴＳ ＬＵＯＹａｎｇ　ＷＡＮＧＴｉｅｈａｎｇ　ＷＡＮＧＪｕａｎｊｕａｎ （ＣｏｌｌｅｇｅｏｆＣｉｖｉｌＥｎｇｉｎｅｅｒｉｎｇ，Ｘｉ′ａｎＵｎｉｖｅｒｓｉｔｙｏｆＡｒｃｈｉｔｅｃｔｕｒｅａｎｄＴｅｃｈｎｏｌｏｇｙ，Ｘｉ′ａｎ　７１００５５） Ａｂｓｔｒａｃｔ　Ａｆｉｎｉｔｅｅｌｅｍｅｎｔｍｅｔｈｏｄｉｓｅｓｔａｂｌｉｓｈｅｄｆｏｒｃａｌｃｕｌａｔｉｎｇｔｈｅｓｅｅｐａｇｅｆｌｏｗｉｎｌｏｅｓｓｗｉｔｈｃｏｎｓｉｄｅｒａｔｉｏｎｏｆｉｔｓ ｊｏｉｎｔｓｉｎｆｌｕｅｎｃｅ．Ｔｈｅｂａｓｉｃｓｅｅｐａｇｅｆｌｏｗｅｑｕａｔｉｏｎｉｎｕｎｓａｔｕｒａｔｅｄｍｅｄｉａｉｓｗｏｒｋｅｄｏｕｔｗｉｔｈｔｈｅｖｉｅｗｏｆｃｏｎｓｅｒｖａｔｉｏｎ ｏｆｍａｓｓ．Ｉｔｉｓａｂｏｕｔｔｈｅｖｅｌｏｃｉｔｙａｎｄｔｈｅｓａｔｕｒａｔｉｏｎｓｅｘｐｒｅｓｓｉｏｎ．ＡｐｐｌｉｎｇｔｈｅＤａｒｃｙｌａｗ，ｔｈｅｖａｒｉａｂｌｅｃａｎｂｅｃｈａｎｇｅｄ ｆｒｏｍｖｅｌｏｃｉｔｙｔｏｈｙｄｒａｕｌｉｃｈｅａｄｉｎｔｈｅｅｑｕａｔｉｏｎ．Ｔｈｅｐａｒａｍｅｔｅｒｉｎｔｈｅｂａｓｉｃｆｌｏｗｅｑｕａｔｉｏｎｏｆｌｏｅｓｓｊｏｉｎｔｓｉｓ ｄｅｔｅｒｍｉｎｅｄａｃｃｏｒｄｉｎｇｔｏｔｈｅｃｈａｒａｃｔｅｒｉｓｔｉｃｏｆｌｏｅｓｓｖｅｒｔｉｃａｌｊｏｉｎｔｓ．Ｕｓｉｎｇｔｈｅｑｕａｄｒｉｌａｔｅｒａｌｉｓｏｐａｒａｍｅｔｒｉｃｅｌｅｍｅｎｔ，ｔｈｅ ｆｉｎｉｔｅｅｌｅｍｅｎｔｆｏｒｍｏｆｔｈｅｅｑｕａｔｉｏｎｉｓｇｉｖｅｎｗｉｔｈｔｈｅＧａｌｅｒｋｉｎｗｅｉｇｈｔｅｄｒｅｓｉｄｕａｌｍｅｔｈｏｄ．ＵｎｄｅｒｌｏｃａｌＤｉｒｉｃｈｌｅｔ ｂｏｕｎｄａｒｙｃｏｎｄｉｔｉｏｎ，ｔｈｅｌｏｅｓｓ２Ｄｓｅｅｐａｇｅｐｒｏｂｌｅｍｉｓｎｕｍｅｒｉｃａｌｌｙｓｔｕｄｉｅｄｃｏｎｓｉｄｅｒｉｎｇｔｈｅｉｎｆｌｕｅｎｃｅｏｆｌｏｅｓｓｊｏｉｎｔ． Ｔｈｅｒｅｓｕｌｔｓｓｈｏｗｔｈａｔｊｏｉｎｔｈａｓａｈｕｇｅｉｎｆｌｕｅｎｃｅｏｎｔｈｅｍｉｇｒａｔｉｏｎｏｆｗｅｔｔｉｎｇｆｒｏｎｔ． Ｋｅｙｗｏｒｄｓ　Ｌｏｅｓｓｊｏｉｎｔ，Ｕｎｓａｔｕｒａｔｅｄ，Ｇａｌｅｒｋｉｎｍｅｔｈｏｄ，ＦＥＭ，Ｎｕｍｅｒｉｃａｌａｎａｌｙｓｉｓ 的促进作用。对于岩石节理来说，由于岩石的低透 ０　引　言 水性，节理岩体的水分运动基本上是通过节理通道 完成的。节理对岩石的渗流起到了重要作用。岩石 ［１～５］ 节理作为优势渗流通道，对水分的入渗有巨大 节理的渗流问题历来得到了充分的重视 。 － － － － 　 收稿日期收到修改稿日期  基金项目：国家自然科学基金，西安建筑科技