基于ANSYS模拟非稳定饱和_非饱和渗流场_岳庆河
基于ANSYS模拟非稳定饱和-非饱和渗流场
ｓ 为储水率 ， 其值较小 ， 一般忽略不计 ； ａ在饱 和区取 １ 非饱和 ，
区取 ０ 为已知水头 ； 边界流量 ； 为土体空 隙 比；ｏ （ ； ｑ为 ｒ ｔ ｃｓ ，
通讯作者 ： 刘福胜 （９ ４ ） 教授 。 １６ 一 ，
Ｅ ｍａｌ ｉ ｆｈ ｓ ａ ． ｄ ． ｎ － ｉ ：ｌ ｓ ＠ ｄ ｕ ｅ ｕ ｃ ｕ
虑饱和区渗 流的方法不 能有效地 解决 水位变化 和降雨 等引 起
的非稳 定渗 流问题 。 已有不少学者应 用 Ａ Ｓ Ｓ进行 了渗流场 数值模 拟研 究 。 ＮＹ 许玉景 等心 首次使用 Ａ Ｓ Ｓ热分析模块结合单元 生死技术进 ＮＹ 行了稳 定渗 流模拟 ， 富强 对此 进行 了进 一步 充实 ， 蒙 罗启北 等Ｈ 使用 ＡＮ Ｙ Ｓ Ｓ进行了导流洞有压渗 流的动态模拟 。为解 决 存在浸润面的非 稳定渗 流场模 拟 问题 ， 者对 Ａ Ｓ Ｓ进行 二 笔 ＮＹ 次开发 ， 编制 了非稳定 饱和 一非饱 和渗 流程序 ， 并通过 实例 验
分析奠定 了基础心“］ 。对二维渗流微分方 程应用变分法 ， 并对
去 ） ＋ Ｈ＋ ［ Ｏ ＝ ［ ］ （ ］啬 ） ［ Ｏ （。 Ｈ ｊ
Ｅ ｈ Ｃ（ ）＋ａ 。 ｓ］
Ｈ（ ｙｚｔ ＝ Ｈｔ ｙｚｔ ， ｆ 上 ，，， ） （ ，，） 在 ｌ ， 。（ ，）＋ｋ ｓｆ ｘ ｒ ０（ ） ｓ
（） ３
２ 非稳定饱和 一 非饱和渗流场数值模拟
２ １ 时 间步长 的 自动选择 ．
非稳定 渗流场数值 模拟 中 ， 选取合 适 的时 间步长 ， 既能保
收 稿 日期 ：０ ７ ０－ ２０ —１－ ９ ０
基于水气二相流的稳定饱和-非饱和渗流模拟研究
基于水气二相流的稳定饱和-非饱和渗流模拟研究
孙冬梅;朱岳明;许艳杰
【期刊名称】《大连理工大学学报》
【年(卷),期】2006(046)0z1
【摘要】非饱和带的渗流过程实质上是水、气两种流体在土壤孔隙中相互替代的过程,因此采用水-气二相流模型求解饱和-非饱和渗流问题更加合理.根据水、空气的质量守恒定律和达西定律,结合多相流理论建立水-气二相流模型,采用高效的积分有限差分法求解,给出精确模拟水相、气相边界的处理方法.通过求解Muskat稳定渗流问题得到逸出面长度与解析解基本一致,验证了水气二相流模型的有效性;由孔隙水压力、孔隙气压力和毛细压力的分布可知,稳定渗流中气相的影响几乎可以忽略,而非稳定渗流中气相的影响有待进一步研究.
【总页数】6页(P213-218)
【作者】孙冬梅;朱岳明;许艳杰
【作者单位】河海大学,水利水电工程学院,江苏,南京,210098;河海大学,水利水电工程学院,江苏,南京,210098;沈阳市水利建筑勘测设计院,辽宁,沈阳,110015
【正文语种】中文
【中图分类】TU457
【相关文献】
1.基于双剪理论的饱和/非饱和渗流边坡稳定性分析 [J], 马宗源;任;党发宁;廖红建
2.基于非饱和-饱和渗流的降雨入渗边坡稳定性分析 [J], 董建军;王思萌;杨晓萧;聂
兰磊
3.饱和-非饱和渗流条件下土堤边坡渗流稳定分析 [J], 石敏香;付国栋
4.考虑饱和-非饱和渗流的土石坝渗流及稳定性计算 [J], 张守仁
5.非饱和地表径流-渗流和流固体耦合条件下降雨入渗对路堤边坡稳定性研究 [J], 刘俊新;刘育田;胡启军
因版权原因，仅展示原文概要，查看原文内容请购买。
饱和-非饱和非稳定渗流的数值模拟
饱和-非饱和非稳定渗流的数值模拟
张培文;刘德富;黄达海;宋玉普
【期刊名称】《岩土力学》
【年(卷),期】2003(24)6
【摘要】介绍了饱和-非饱和渗流的计算程序数值模拟方法,该方法模拟了降雨过程中边坡内孔隙水压力的变化,克服了传统降雨边界处理方法的缺陷,考虑了降雨过程中地表边界的条件转化,为以后的径流渗流耦合的模拟提供了良好的基础,为降雨诱发滑坡研究提供了定量分析手段。
【总页数】4页(P927-930)
【关键词】饱和—非饱和渗流;数值模拟;降雨边界;滑坡;孔隙水压力
【作者】张培文;刘德富;黄达海;宋玉普
【作者单位】大连理工大学;三峡大学水电学院
【正文语种】中文
【中图分类】P642.22;TV139.1
【相关文献】
1.三维饱和非饱和稳定非稳定渗流场的有限元模拟 [J], 张家发
2.饱和-非饱和非稳定渗流有限元分析方法的改进 [J], 周桂云
3.饱和-非饱和非稳定渗流数值分析中初始状态的研究 [J], 袁俊平;褚飞飞;季李通
4.大型贮灰场的三维饱和–非饱和瞬变流数值模拟及非饱和参数试验 [J], 丁家平;
徐辉;顾国新;李勤明;王富庆
5.饱和非饱和土体非稳定渗流数值分析 [J], 吴梦喜;高莲士
因版权原因，仅展示原文概要，查看原文内容请购买。
浅谈基于ANSYS的土石坝渗流场模拟
介 质 的分界线应 该与 非均质 岩层 的分界 线也 要保 持
一
情况 ， 由于 Ａ Ｓ Ｓ软件 没有 提供模 拟 图的给水 度 的 ＮＹ 功能 ， 不能直接 分析 非 稳 定 渗 流 的情 况 和渗 流 问题
本身 的复杂性 ， 采用 数值 模拟方 法 ， 如何真 实准确 的 模 拟实 际渗流是 一个 需要 不断完 善 的过 程 。从 而研
参考文 献 ：
［ ］ 毛昶熙 ． 流计算分析与控制 ［ ． １ 渗 Ｍ］ 北京 ： 中国水利水
电 出 版社ห้องสมุดไป่ตู้，ｏ３ ２ｏ．
４ 工 程 实例 分 析
本 文利用 西安市 金盆水 利枢 纽工 程为例 。该 工 程 为大 型二等工 程 ， 于黏 土心墙 土石 坝 ， 址位 于 属 坝
西安 市 的周至县境 内 ， 周至 县城 约 １ ｋ 距 ３ ｍ。该 枢纽
有重要 地位 。对 土石 坝 而 言 ， 透水 流 除浸 湿 土 体 渗 降低其 强度指 标外 ， 当渗 透 力大 到 一定 程 度 时将 导 致坝坡 滑动 、 防渗 体 被击 穿 、 坝基 管 涌 、 土 等重 大 流 渗流事 故 ， 接 影 响大 坝 的运 行安 全 。对 于 混凝 土 直
受力安 全 ； 岸坝 肩渗透 压力 （ 下水 位 ） 两 地 的高 地关 系到坝 肩岸坡 岩体 的抗 滑稳定安 全 。带 有浸蚀 性 的
渗流对 建筑 物和 坝基 的可溶 性 物质 造 成 浸蚀 ， 响 影
结构安 全 问题 。此 外 ， 大 的渗 透 损 失也 将 减少 工 过
程效益 。
（ ） 散化 ： 求 解 域 划 分 为具 有 一 定 几 何形 ２离 将 状 的单元 进行单 元 编号 并 确定 插 值 函数 ， 对结 点 进 行 总体 编号和单元 上 的局部 编号并 给 出结 点局部 标
基于ANSYS的复杂地质条件下大型渗流场计算
基于ANSYS的复杂地质条件下大型渗流场计算水电站设计DHPS第26卷第4期2010年12月基于ANSYS的复杂地质条件下大型渗流场计算段斌,何江达,王瑶(1.国电大渡河流域水电开发有限公司,四川成都610041;2.四川大学水利水电学院,四川成都610065;3.四川电力进出口公司,四川成都610061)摘要:利用ANSYS温度场分析模块计算大岗山水电站复杂地质条件下的大型渗流场,取得了较好的效果.关键词:ANSYS软件;渗流场;计算方法;应用软件;大岗山水电站中图法分类号:TV223.6文献标识码:B文章编号:1003—9805(2010)04—0042—04 1ANSYS进行渗流分析的理论基础1.1ANSYS软件基本情况ANSYS软件是融结构,热,流体,电磁,声学于一体的大型通用有限元分析软件.该软件具有强大的前处理及后处理功能,它的图形界面和交互式操作大大简化了计算模型的创建过程,同时在计算之前,可通过图形显示来验证模型的几何形状,材料及边界条件;在后处理中,其计算结果可以采用多种方式输出,比如计算结果排序和检索,彩色云图,等值线,动画显示等等.与其他通用软件比较,其前,后处理功能优于同类型的软件.ANSYS为我们提供了一个优良的工作环境,使我们从繁琐,单调的常规有限元编程中解脱出来.该软件可以进行热,电, 磁,流体和结构等有限元分析,并可以进行多物理场耦合分析.由于渗流场和温度场可以互相比拟,故采用ANSYS的温度场分析功能进行渗流场的计算,可以取得很好的效果¨.1.2ANSYS温度场分析模块分析渗流场理论基础温度场与渗流场之间具有理论基础相似,微分方程相似,初始条件与边界条件的相似等特点. ANSYS中渗流场与温度场各种相应量比较见表1[一422ANSYS参数化设计语言及其在本文中的运用ANSYS参数化设计语言(ANSYSParameter DesignLanguage,简称APDL)是一种可用来自动完成有限元常规分析或通过参数化变量方式建立分析模型的脚本语言,用建立智能化分析的手段为用户提供自动完成有限元分析过程,即程序的输入可设定为根据制定的函数,变量以及选用的分析类型来做决定,是完成优化设计和自适应网格的最主要的基础_5J.APDL允许复杂的数据输入,使用户实际上对任何设计或分析属性有控制权,例如尺寸,材料,荷载,约束位置和网格密度等...APDL有以下功能:(1)参数;(2)表达式和函数;(3)分支和循环;(4)重复功能和复写;(5)宏;(6)用户子程序.APDL有很多优点,其中之一是可以利用AP-DL从事二次开发.本文借助APDL编制了相应的APDL模块来自动迭代计算出渗流边界和自由面.为了便于与结构计算网格相对应,计算中采取了固定网格法(不变网格法).其中干燥区(自由面以上的区域)的渗透系数取很低的值(本文在具体计算时,干燥区渗透系数取1.0X10cm/s),从而可以实现干燥区和饱和区进行联立统一求解,避免了干燥区渗透系数取为0时导致渗透矩阵奇异而无法进表1ANSYS中渗流场与温度场各种相应量的比较收稿日期:2009一O2—25作者简介:段斌(1980一),男,四川绵阳人,工程师,主要从事水电工程技术和管理工作.行求解.具体计算时,首先将整个计算域视为饱和区,已知的边界条件以实际作用范围输入到计算模型中,通过反复迭代便可最终确定自由面位置.3工程应用3.1工程概况大岗山水电站是大渡河干流规划中的第十四级电站,位于四川省雅安市石棉县境内.工程等级为一等大(1)型工程,开发任务以发电为主,电站总装机容量2600MW,设计多年平均发电量114.3亿kW?h.水库正常蓄水位1130m,总库容7.42亿Ill,调节库容1.17亿IIl,具有日调节能力.电站枢纽工程由最大坝高210m的混凝土双曲拱坝,左岸地下厂房,右岸泄洪洞等组成.3.2地质条件大岗山水电站坝址区河段呈"Q"形嵌入河曲,两岸山体雄厚,谷坡陡峻,河谷狭窄且对称.坝址区地层岩性较为单一,主要为澄江期酸性花岗岩,辉绿岩脉(13)等穿插发育于花岗岩中.厂坝区发育的岩脉,断层众多,其产状以近河流向(NE,SN或NW向)中陡倾角发育为主,垂河流方向(近EW向)发育较少.另外,坝址区还发育有多组裂隙.坝址区发育的岩脉(断层)和多组裂隙共同构成了厂坝区地下水的渗流通道.3.3天然渗流场反演天然渗流场反演分析的目的主要为了确定大岗山水电站三维渗流场的边界水位及岩体的渗透参数, 作为运行期渗流场研究的依据.大岗山水电站天然渗流场反演分析的方法和主要结论参见文献[7].通过反演分析,厂坝区各类岩体的渗透张量见表2.3.4运行期渗流场计算3.4.1计算模型根据大岗山水电站厂坝区主要水工枢纽布置,防渗帷幕布置,排水幕和排水廊道布置情况以及地质地形条件,并考虑计算域的边界效应,三维渗流场有限元计算范围上游取至导流洞进El位置,下游取至导流洞出口位置,顺河向长度约883.0m;横河向左侧边界取至地下厂房洞室群所在的地形分水岭附近,右侧边界取至靠近地形分水岭位置,横河向长度约1423.8m;垂直方向底面取至海拔615.0m.有限元计算坐标定义为:轴:顺河流方向,从上游指向下游为正,轴向方位为SW238.;y轴:垂直河流方向,从右岸指向左岸为正,轴向方位为SE148.;z轴:与x和y垂直,且Z=x×Y,铅直向上.表2厂坝区各向异性岩体的渗透张量雾[芝主.../萎cm~曩s一-1.01×10_41.一等透水组裂隙【对称L822.96.L56m.:3.57x10-4:2.04x10-5.圳.,圳一弱透水组裂隙【对称…;:J………s整个计算域大部分采用空间八节点等参单元进行离散,局部考虑材料介质过渡和地形变化等因素退化为三棱柱和四面体单元,其中运行期计算域共离散为53959个单元和54957个节点.运行期三维网格图见图1,2.3.4.2排水洞(廊道),排水幕的模拟及有关材料的渗透系数在计算模型中,排水洞(廊道)的位置和范围都以节点形式详细模拟.由于排水洞(廊道)在正常运行工况基本不存在失效和出现阻塞的情况,因此其水头按排水洞(廊道)所在高程约束,在迭代计算时,根据排水洞(廊道)是否有排水流量和其水头与位置势是否相等这两个条件来计算渗流自由面. 排水幕的排水效果是通过等效结构面的等效渗透系数来体现的,等效结构面在计算模型中被准确模拟出来,排水幕的等效渗透系数,见表3,主要图1大岗山厂坝区运行期渗流场三维有限元网格43图2大岗山厂坝区运行期坝体+灌浆帷幕+地下厂房三维网格受排水孔尺寸,间距和排数的影响.防渗帷幕在计算模型中以单元形式模拟,有关材料的渗透系数见表4.表3各渗控方案排水幕的等效渗透系数表4设计提供运行期厂坝区部分材料的渗透系数cm/s3.4.3运行期厂区渗流场分析厂区典型剖面的地下水位等势线图,渗透压力等值线图和厂坝区枢纽整体模型的地下水水位(自由面)等值线分别见图3,4.经分析可知:(1)从图3可以看出,由于厂区顺河向防渗帷幕,帷幕后的排水幕(两排)和4层排水廊道(高程1032m,987m,957m,931.4m)的"前堵后排"作用,库水向厂区渗流的自由面在防渗帷幕之后附近显着降低.另外,由于B岩脉的渗透系数远大于周围岩体的渗透系数,防渗帷幕在靠近库水一侧的自由面也有一定的降低.厂区三大洞室中,地下厂房,主变室周边地下水位线较低;对应的外水压力较小,而调压室周边地下水位线相对较高;外水压力相对较大,主要受靠山内一侧的边界水位和运行期调压室水位影响所致.(2)从图4可以看出,靠山内侧的调压室,由于防渗帷幕没有延伸至调压室,而仅有987m高程的排水廊道排水作用,在山体内较高的天然地下水渗流作用,调压室上游侧的地下水位线较高;渗透压力较大,而主变室和地下厂房所在区,虽然山体天然地下水位较高,但是由于上游侧的防渗帷幕,灌浆廊道和4层排水廊道(1032m,987m,957m,931.4m)的"前堵后排"作用,在主变室和地下厂房上游侧地下水位线显着降低,渗透压力也较小.(3)计算成果表明,设计方案的防渗,排水措施能有效地降低厂区三大洞室(尤其是地下厂房和主13t512i51Il5鲁1015\9159157l56I50lO02002004005006007008009000lO0200300400500600700800900 图3运行期厂区三大洞室横剖面水位/渗压等值线∞撕∞啪鲁\13151215llt5吕]015\,158I57】56佰埘5I￡l5【llj蛊i.15\N翱58巧7156】5图4运行期厂区主机间纵剖面水位/渗压等值线变室)周边的地下水位线,厂区三大洞室周边承受的外水压力不大,调压室底部最大压力水头约50m, 地下厂房最大压力水头约10m,而主变室基本位于干燥区,表明设计方案的渗控效果显着.(4)计算成果还表明,厂区1032m高程的排水廊道基本处于干燥区,其对厂区高高程的排水作用不大,可以考虑取消.3.4.4运行期坝区渗流场分析坝区典型剖面的地下水位等势线和渗透压力等值线见图5.经分析可知:(1)由于坝基防渗帷幕对岩体,尤其是坝肩较高高程的弱透水岩体段的阻渗作用,以及帷幕之后排水幕和坝肩排水洞的强排水作用,坝基帷幕前后的水头损失很大,排水幕之后的坝底扬压力很小,坝基建基面925m高程排水幕处扬压力水头减小约90m,位于两岸坝肩中部的坝底扬压力水头减小约120～150m,位于两岸坝肩上部的坝底扬压力水头O10O啪300枷5006o0700800减小约80—100m.(2)坝后水垫塘左岸抗力体排水洞区的自由面位置较低,右岸抗力体排水洞区的自由面位置较高. 由于坝体部位的帷幕,排水幕和排水洞的"前堵后排"作用,使得坝后水垫塘两岸抗力体所在区域的地下水渗流主要受两岸岩体中天然地下水位(右岸高,左岸低)控制.3.4.5运行期厂坝区渗流量通过厂坝区三维渗流场计算,可以获得前述不同计算方案通过厂区排水廊道,地下洞室,防渗帷幕,坝区排水廊道,抗力体排水洞的渗流量,见表5. 4结论(1)本文以大型通用有限元计算软件ANSYS为平台,使用ANSYS温度场分析模块,并用APDL 编制了相应的计算程序,计算了三维渗流场.该方11151015g们5\B157l66150lo0瑚300400500600700800X/mX/m(a)渗压等值线(b)水位等势线图5运行期坝体最大横剖面水位和渗压等值线表5运行期厂坝区各部位渗流量m./d(下转第54页)45图5材料赋值,施加约束和荷载后的有限元模型(a)主压应力等值线p0¨E~zaphic￡图6典型等值线esel,S,elem,,1,1081单元选择sfgrad,pres,,Z,2520.00,一9810.001水荷载施加sfe,all,1,pres(7)计算成果检查.根据已经建好的三维拱坝模型进行线弹性有限元计算,典型的位移,应力等值线见图6.,经与拱梁分载法计算成果比较,其分布规律和拱梁分载法吻合良好,极值差别不大,有限元计算结果可靠.5结语根据ANSYS有限元软件的内在规定,用编程语言编制大坝ANSYS---维模型自动建立的方法和(b)顺河向位移等值线p.'地0r●phjc●思路,能大幅节约建模时间,提高工作效率.(1)直接进行节点,单元建立的三维模型建立方法,节点,单元可根据需要随意调整,且不易发生单元奇异,有利于后面的计算收敛,减少奇异单元检查工作.(2)合理的单元规划很重要,根据建筑物的重要性和受力特点进行单元划分分组,并以每一组单元的交界作为单元控制面进行不同组单元的衔接, 单元的疏密能满足计算的要求,单元编号有序,更易控制.(3)前处理采用生成ANSYS命令流文件的方式进行模型材料的赋值,约束及荷载的施加,方便,准确,且便于修改.(上接第45页)法将ANSYS软件与复杂渗流场计算结合起来,方便人们进行大型复杂模型的建模和后处理工作,所建立的有限元计算模型不仅适用于渗流场的计算,还可用于应力场的计算.(2)大岗山水电站地质情况复杂,运行期厂坝区渗流场计算模型复杂,单元众多,渗控措施模拟难度大,采用ANSYS温度场分析模块后,由于其强大的前处理,后处理,二次开发的功能,可以较方便地对大岗山水电站厂坝区渗流场进行模拟计算,并取得较好的效果.参考文献:54[1]ANSYS中国公司.ANSYS高级技术分析指南[M].北京:美国ANSYS北京办事处,2005.[2]任辉启.ANSYS7.0工程分析实例详解[M].北京:人民邮电出版社,2004.[3]李景涌.有限元法[M].北京:北京邮电大学出版社,1999.[4]李军华.大坝渗流监测系统设计及渗流计算机模拟[D].郑州:郑州大学,2004.[5]龚曙光,谢桂兰.ANSYS操作命令与参数化编程[M].北京: 机械工业出版社,2004.[6]刘涛,杨凤鹏.精通ANSYS[M].北京:清华大学出版社,20o5.[7]唐正州,段斌,何江达,等.大岗山水电站厂坝区三维天然渗流场反演分析[J].四川水利,2007,(1).。
ANSYS非稳态热分析及实例详解
第7 章非稳态热分析及实例详解本章向读者介绍非稳态热分析的基本知识，主要包括非稳态热分析的应用、非稳态热分析单元、非稳态热分析的基本步骤。
本章要点非稳态导热的基本概念非稳态热分析的应用非稳态热分析单元分析的基本步骤本章案例钢球非稳态传热过程分析不同材料金属块水中冷却的非稳态传热过程分析高温铜导线冷却过程分析7.1 非稳态热分析概述物体的温度随时间而变化的导热过程称为非稳态导热。
根据物体温度随着时间的推移而变化的特性可以区分为两类非稳态导热：物体的温度随时间的推移逐渐趋于恒定的值以及物体的温度随时间而作周期性的变化。
无论在自然界还是工程实际问题中，绝大多数传热过程都是非稳态的。
许多工程实际问题需要确定物体内部的温度场随时间的变化，或确定其内部温度达到某一限定值所需要的时间。
例如：在机器启动、停机及变动工况时，急剧的温度变化会使部件因热应力而破坏，因此需要确定物体内部的瞬时温度场；钢制工件的热处理是一个典型的非稳态导热过程，掌握工件中温度变化的速率是控制工件热处理质量的重要因素。
再例如，金属在加热炉内加热时，需要确定它在加热炉内停留的时间，以保证达到规定的中心温度。
可见，非稳态热分析是有相当大的应用价值的。
ANSYS 11.0及其相关的下属产品均支持非稳态的热分析。
非稳态热分析确定了温度以及其它随时间变化的热参数。
7.1.1 非稳态热分析特性瞬态热分析用于计算一个系统的随时间变化的温度场及其它热参数。
在工程上一般用瞬态热分析计算温度场，并将之作为热载荷进行应力分析。
瞬态热分析的基本步骤与稳态热分析类似。
主要的区别是瞬态热分析中的载荷是随时间变化的。
为了表达随时间变化的载荷，首先必须将载荷-时间曲线分为载荷步。
对于每一个载荷步，必须定义载荷值及时间值，同时必须选择载荷步为渐变或阶越。
7.1.2 非稳态热分析的控制方程热储存项的计入将稳态系统变为非稳态系统，计入热储存项的控制方程的矩阵形式如下：[]{}[]{}{}C T K T Q +=其中，[]{}C T 为热储存项。
饱和—非饱和渗流作用下边坡稳定分析的混合法
ｆ ）作 用 在 单 兀 土 体 之 间 的 土 压 力 大 小 相 ３ 等 向 相 反 ． 滑 动 体 整 体 束 看 ，】 为 相 互 抵 方 从 ｌ Ｊ
消， 对单 元体 分析 时也 不 考虑
和 渗 流 作 用 下 的 边 坡 稳 定 问 题 进 行 Ｉ研 究 ， 广 利
２ 滑 坡 计 算 的 有 限 元 法
２ １ 方 法 简 述 ．
用 渗流 计算 时 的剖分 网格 （ 四边 形等参 元 ） 和计算
结 果 ， 接 连 续 进 行 渗 流 作 用 下 的 边 坡 稳 定 分 析 直 在 寻找 最 小 安 全 系 数 及 其 相 应 滑 弧 位 置 时 ， 用 采 数 学 规 划 中 的 单 纯 形 法 ， 现 了 自动 查 找 的 日的 实 计 算 结 槊 表 明 ， 法 为 饱 和 一非 饱 和 不 稳 定 渗 流 该 作 用 下 的 边 坡 稳 定 分 析 问 题 提 供 Ｉ一 个 实
文 章 编 号 ：６ ｌ ６ ３ （ ０ ２ ０ ｏ ２ ０ １７ ８ ３ ２０ ） １ ０ ５ ３
饱 和 一非 饱 和 渗 流 作 用 下 边 坡 稳 定 分 析 的混 合 法
汪 自力 ，朱 明霞 ，高 青 伟 ，岳瑜 素
ｔ ｌ黄 州 水 利 委 员 会 黄 河 水 利 科 学 研 究 院 ． 南 部 州 ４ ０ ０ ：２ 郭 州 大 学 ｌ程 力 学 系 ， 南 部 州 河 ５０３ 河 ４ ０ ０ ； ． 『 建筑 ｌ 公 司 。 南 郑 州 ４ ０ ０ ） ５０２ ３ ｉ南 ｂ 程 河 ５０２
维普资讯
２） ４ ２年 ￣ ３ Ｈ 第 ２ ３卷 第 ｌ期
郑 卅 ｌ＾ 学 学 报 【 ［ 学 版 ｊ
ｈＬ １ｌ ｆ ｈ ｎ ｚ ｎ ｖ ｒｉ （￣ｇ ｎ ｅ ｉｇ Ｓ ｉ ｃ ） ｌ ａ ｏ ｅ ｇ ｝ Ｕ ｉ ｅｓ ｙ Ｆ ｉｅ ｒ ｃｅ ｅ ｎ Ｚ ｍｕ ｔ ｎ ｎ ｎ
超低渗透储层饱和_非饱和渗流启动压力梯度试验研究_谢全
静态毛管压力 /kPa
120.58 125.40 －327.28 －119.20 －51.68
第 30 卷 增 1
谢 全等：超低渗透储层饱和/非饱和渗流启动压力梯度试验研究
上游端压力传感器
下游端压力传感器
恒压注入
上游端死体积
岩心夹持器
下游端死体积
• 3131 •
高精度压差传感器
图 1 非稳态法启动压力梯度测试原理 Fig.1 Test principle of threshold pressure gradient with unsteady state method
本文通过非稳态法，当岩心渗流达到稳态后(岩 心两端压力差在 1.4 kPa 范围内波动，流量波动范围 小于 1%)，关闭岩心出口端阀门，同时保持岩心入 口端恒压。在压力差条件下，岩心出口端将有流体 流入下游端死体积(管线死体积为 6.4～6.5 mL)，由 于出口端死体积部分阀门关闭，因此出口端管线死 体积部分压力将逐渐增加。因为液体压缩系数小 (2×10－4～6×10－4 MPa－1)，即便有微小流量(10－6～ 10－7 mL/min)从岩心出口端端面流出，也可以导致 10 kPa 左右的压力变化。因此，采用非稳态法完全 可以通过岩心两端的压力差变化来判断是否存在启 动压力梯度。如果超低渗饱和渗流存在启动压力梯 度，那么岩心下游端压力将永远低于岩心上游端压 力。相反，如果一定时间后岩心两端压力差为 0， 就可以说明超低渗饱和渗流不存在真实的启动压力 梯度。非稳态法启动压力梯度测试原理见图 1。 2.2 试验设备
PERMEABILITY RESERVOIRS
XIE Quan1，2，HE Shunli2，JIAO Chunyan2，ZHU Huayin3，FU Chunyan4，LEI Gang2
基于ANSYS的土石坝渗流与稳定分析研究的开题报告
基于ANSYS的土石坝渗流与稳定分析研究的开题报告一、研究背景和意义土石坝是一种重要的水利工程结构，其安全稳定性直接关系到人们的生命财产安全和社会经济发展。
而渗流问题是土石坝安全稳定性研究的重要内容之一。
在土石坝工程设计、施工和运行过程中，渗流问题一直是困扰工程师的难题，如何在渗流对土石坝安全稳定性产生影响的情况下，保证土石坝的安全运行是当前亟需解决的问题。
本研究利用ANSYS软件，通过有限元数值模拟方法，研究土石坝内部的渗流分析及土石坝的稳定性分析，旨在探讨土石坝渗流及其对稳定性的影响规律，为土石坝设计、施工和运行提供技术支撑和依据。
二、研究内容和方法本研究的主要内容有两个方面：一是土石坝内部渗流模拟及分析；二是土石坝稳定性分析。
具体通过以下步骤实现：1. 确定研究对象：本研究以某一具体土石坝为研究对象，对其渗流分析及稳定性进行模拟和分析。
2. 建立土石坝模型：根据实际情况建立土石坝三维有限元模型，包括坝体、坝基、边坡等，考虑土、石材料的物理力学特性。
3. 渗流模拟：对建立的土石坝模型进行渗流模拟，通过ANSYS中的多孔介质渗流模型，对土石坝内部流场进行数值计算和分析。
4. 渗流分析：根据渗流模拟结果，分析产生渗流的原因，判断坝体、坝基是否产生渗漏现象，并分析渗漏现象的破坏机理。
5. 稳定性分析：根据建立的土石坝模型，通过ANSYS有限元分析软件对土石坝的稳定性进行数值计算，分析坝体的变形、破坏状况，确定安全系数，预测土石坝的破坏条件。
三、预期成果和意义本研究通过ANSYS软件，对土石坝内部的渗流分析及稳定性分析进行研究，预期取得以下成果：1. 对土石坝内部渗流模拟及分析技术的研究与应用，提高土石坝设计、施工和运行的水平，为工程师在实际工程中提供技术支撑和依据；2. 对土石坝安全稳定性分析方法的探究和应用，为土石坝的安全设计和管理提供科学依据，提高工程的安全性和经济效益；3. 深入了解土石坝渗流及其对稳定性的影响规律，为水工、环境等领域的科研人员提供参考，促进相关学科的发展。
基于饱和-非饱和理论的尾矿库渗流场模拟分析
90工业安全与环保2013年第39卷第10期I ndus t r i al Saf et y a nd Envi w nnm [1t al Pr o t ect i onO ct ober 2013基于饱和一非饱和理论的尾矿库渗流场模拟分析朱毅川1徐克2程镇2杨楠2(1．中国地质大学工程学院武汉430074；2．湖北省安全科学技术研究中心武汉430070)摘要在稳定性分析过程中，尾矿库力学参数选取困难，尾矿库边界条件复杂，需要能够较为精确预测尾矿坝浸润线位置和稳定性计算的方法。
在饱和一非饱和渗流理论研究的基础上，对尾矿库渗流场的分布进行了数值模拟，通过对比研究，饱和一非饱和模型比饱和渗流模型更加符合尾矿库的渗流规律。
研究了初期坝透水性与干滩长度对尾矿库浸润线的影响，对比分析了影响尾矿库渗流的一些关键因素，以期能更好地解决尾矿库渗流稳定性评估问题。
关键词尾矿库稳定性分析数值模拟饱和一非饱和理论T he S i m ul a t i on of T ai l i ngs Seepage Fi el d B as edont heSaturat ed —-嘲turat edT heor yZ H U Y i c hua nlX U K e2C H EN G Zben2Y A N G N an2(1．Eng ／neer ／ngFacul t y ，C hi na U ni ver si t yofC ．eos ci encesW uba n430074)A bst r actI n t he s t abil i t y anal y si s of t he r ai l i ng s ，t he sel e ct i on of t he m ec ha nl e al para m et e r s of t he t ai li ngs is usua l l y i na c-cu _rl l t ea ndSOi t i s ba dl y i n need ofa m or e ac cu r at e m et h od t o f or ecastt aat ,峥dam s ee pa ge l l ne l ocat i on a nd s t abi l it y ca l cu-l at i on becaus e of t he bou ndar y condi t i ons of t he t a i l i n 铲．Base d o nt he sa t ura t e d —unsa t ura t e d s ee pa ge t heor y ，t he di st r i bu-t i m 0f t he t ai li ng pond s ee pa ge f i el d is si m ul a t e d i n t his paper ．．I l II oL 曲com parat i ve gt udy ，i tis f ound t hat sat ur at e d —t i nsa t ura t ed m odelc anbet t er agr e et o t het aaa 啦seepage r ul e t han sat ur at ed s ee pa ge m odel ；t he i nf l uence s of t he i ni ti al da ml m m eabi l i t y a nd dr y l e ngt ho nt het aj l i I 骖i nf i l t r at i ng l i near er es ea r ched ．The paper ha sal so 舡l a 】咖脚80nl e keyf act or s af-f ect i ngthe 咖seepagei n or de r t o bet t er addr es s t he r ai li ngs s ee pa ge st abi l i t y asse ssm en t ．K ey Wor ds t ai l ings ponds t abi l it ymla ：Iys i 8num er i c alsi m ul a t i onsa t ura t e d —unsa t ura t e d t heor y0引盲尾矿库浸润线的高低是影响尾矿堆积坝稳定性的主要因素之一。
ANSYS软件在平原水库土工膜气胀计算中应用论文
ANSYS软件在平原水库土工膜气胀计算中的应用摘要：本文分析了热传导问题温度场与土工膜下非饱和土气场的基本理论、微分方程以及边界条件，将ansys软件的热分析模块应用于膜下气场的计算，提出了应用温度场理论解决非饱和土气场问题的求解方法。
该方法可以计算复杂边界和多种介质的稳态或非稳态气场问题，为工程设计提供了极大便利。
关键词：ansys；热分析；土工膜；气胀abstract: this paper analyzes the heat conduction problem temperature field and earthwork film the unsaturated rustic field, the basic theory of differential equation and boundary condition, the ansys software of thermal analysis module used in the calculation of the gas membrane, and puts forward the application temperature field theory to solve a problem of unsaturated rustic method. this method can be calculated boundary and a variety of medium complex steady or unsteady gas problem, for engineering design provides great convenience.keywords: ansys; thermal analysis; geomembrane; flatulence中图分类号：tv62 文献标识码：a文章编号：气胀问题是近年来造成平原水库土工膜防渗体系破坏和失效的主要原因之一[1]，一旦发生将使水库产生严重的集中渗漏，因此合理地布置排气设施是土工膜防渗体系设计的核心问题。
用ansys做非稳定渗流与边坡稳定分析第10页均为动画表示的渗流场随时步变化的情况_后就能看到了
基本理论
非稳定饱和—非饱和渗流控制方程
x
kx
h
H x
y
ky
h
H y
z
kz
h
H z
C
h
Ss
H t
非饱和土抗剪强度
c tan ua uw tanb
基本理论
二、软件之间是相互联通的，如CAD，3D在建模 上与ANSYS的相似性。可以在不同软件间取长补 短，结合运用。
三、对ANSYS错误记录文件是操作失误时的重要 提示，为如何改正指出了方向。
二、浅谈ANSYS学习方法
四、试验性。从最简单的例子入手，有什么想法， 可以做一最简单的例子，使用软件进行验证。
非稳定渗流程序结构框图
非稳定渗流程序结构框图
！调整单元渗透系数的语句
/PREP7 MATNUM=2 *DO,I,1,E_MAX ePRESS=BBB(i)-ccc(i) *IF,ePRESS,GT,0,THEN MPCHG,1,I *ELSEIF,ePRESS,LT,-13.812,THEN MPCHG,2,I *ELSE MP,KXX,MATNUM+1,kxxbh*hsl(ePRESS) MPCHG,MATNUM+1,I MATNUM=MATNUM+1 *ENDIF *ENDDO
ANSYS使用经验交流
作者邮箱：erep@
汇报安排
一、简介使用ANSYS对课题的研究情况 二、浅谈ANSYS学习方法 三、相互交流学习
一、简介使用ANSYS对课题的研究情况
基于ANSYS,开发出了土石坝渗流与边坡稳 定分析程序.
基于ANSYS的土石坝渗流数值模拟研究
防 渗方 案 的科学性 与 合理性 ， 从 而 为土石 坝 工程 设计
０ 引言
十 石坝 是一种 常 见 的水 工建 筑物 ， ＿ 它 对地质 要 求
和 建设提 供 了一 定 的参考和 借 鉴 。
１ 基于 Ａ ＮＳ ＹＳ的渗流计算的基本原理
相 对较 低 、施 工简单 ，能就地 取材 ，造 价 相对 经济 ，
（ １ ．湖南省 水运 建设投 资集团有 限公 司，湖 南 长沙 ４ １ ０ ０ ０ ０ ；２ ．大连 海洋大学，辽 宁 大连 １ １ ６ ０ ２ ３）
摘 要 ：基 于有限元原理 ，采用 Ａ ＮＳ ＹＳ软件 中的热处理模 块 ，建立 了坝体数值 分析模 型，针 对湖南省 长沙县五 龙 山 土石坝的渗流过程进行 了数值模拟 ，得到 了在三种 工况下坝体加 固前后 的渗 流速度 、渗流量 、渗透坡 降、选 出点 高 程 、渗透比降等特征值 ，通过比较 可知 ：在 正常蓄水位 、设计 水位及校核 洪水位条件 下，采用 粘土防渗墙进行坝体 的加 固，其渗流速度 、单宽渗 流量、逸出点 高程 以及 逸出点渗透 比降均 小于坝体加 固前计算相应值 ，平 均渗透坡 降 高于加 固前计算相应值 。计算结果表 明采用 粘土防渗处理 方法进行坝体加 固，对保证 土石坝 的渗流稳定和坝体安全
土石 坝加 固前 后渗 流 的数值模 拟 结果 比较 ， 肯定 了该
洮河海甸峡堆石面板坝在稳定状态和非稳定状态下渗流分析研究报告大学方案
1绪论1.1 本文研究得意义大坝作为一种挡水构筑物一般都会涉及到渗流稳定性问题,并且它们的破坏主要都是从渗流开始。
如若不及时引起重视进行处理,必将危及大坝的整体稳定性。
坝体发生渗流破坏的原因,通常都是由于水力梯度引起的。
因此，要研究确定水工构筑物的水力梯度的分布 ,同时必须通过其进行渗流分析 ,确定浸润自由面，再对坝体进行渗流稳定分析，以此对大坝的安全运行进行作出合理评估。
水利枢纽工程在蓄水期，由于库区内水位的上升，在水压力的作用下会在大坝的坝体内部和大坝基底产生渗流的现象，这对水工建筑物的正常运行会产生比较不利的影响，可是又是不可避免的。
根据统计，我国已建成上万座各种类型水库，经过长年累月的运行，很多坝体在不同程度上都存在严重的病害或者隐患，这种病险危害的存在会从根本上影响水工建筑物的渗流稳定性及其变形稳定。
近年来，在水利枢纽工程的正常运行过程中发生了很多具有灾难性的的破坏，经过广泛分析研究后，分析得出在可能对坝体造成破坏的所有模式中，渗流都在其中扮演着重要的作用。
因此，对大坝渗流的研究是大坝安全检测的一项重要内容，深入的分析和研究坝体渗流状况，这对于我们水工建筑物的安全运行以及应对具有破坏性事故的威胁和保护人民生命财产具有重大而深远的意义。
目前较为理想的一种计算渗流场的方法是有限单元法。
在渗流场的基础上对其进行整体稳定性分析 ,才能得到较为可靠的结果。
因此 ,用有限单元法计算其渗流场，对挡水构筑物安全、维护、治理等都有巨大的社会效益和经济效益。
有限单元法的基本思路和原理早就有人提出，它是一套求解微分方程的数值计算法，与传统的解题法相比，有物理意义直观明确，解题效能高方法本身的适应性强，解题步骤易规范化，电子计算机配合使用后，有限单元法得到了迅速的应用和发展。
1.2 国内外渗流计算研究的发展现状法国工程师达西 (H.Dracy>，在 1986 年通过实验研究提出了线性渗透定律，为渗流理论的发展莫定了基础。
应用ANSYS方法揭示海水循环泵蜗壳流道内部非稳态流场的研究
及 叶 轮流 道进 行 了 数 值 分析 。 分析 研 究其 内 部流 动
特性 为 判定 其水 力 设计 的 合 理性 、 模 型 的改 进 和优
化提 供 了 依 据 。
１ 计算 输 入
１ ． １ 物 理 模 型 的 建 立
应用 Ｕ Ｇ 软 件 对海 水 循环 泵 进 水 流 道 蜗 壳 流 、
．１
ｍ，
ｎ ＝ １ ６９
ｍ ｒ／
ｎ ｉ
，
在 该 况 下 进 行 边 条 件 设 置 ｐ＝
１
０３
０ｋ
ｇ
３
／ｍ
。
工
界
的
，
进 边 界 条 件设 置 为 流 量进 ｎ
口 ，
４ （＾ ＝２ ９ ７
１ ． ２ ５ ｋ ｇ ／ｓ 。
３ １ ．
基 本 方 程
假 定 海 水 循 环 泵 内 流 体 流 动 为 不 可 压 缩 、 非 定
４ ３ １
．
．
壁面 边 界 条 件
壁 面 条 件用 来 限 制 流体 的 边 界 计算 中 取粘 性 ，
流 体无 滑 移 条件 即 时均 速 度 和 脉动 速 度 的 各 分
为 ６０／（
１
６ ９ ｘ １
００
）
ｓ ，
是 时 步 算 间 间 隔 即
间
长 为 计 ，
０
００３
．
５ ５ ０２ ９
５ ８ ６ｓ。
叶
轮旋转 ８ 周 的情 况
总时 组 则
间 为 两 （ ８ ｘ ６ ０ ）／ １ ６ ９ Ｓ ｄ
，
接触 面 和 设 置 为 Ｄ ｒ ｆ
基于van Genuchten-Mualem模型的饱和-非饱和介质流动随机数值分析
了更为全面的掌握土壤参数空间变异性对地下水运动的影响 还必须进一步探索多随机输入情况下的
随机模型分析方法
参考文献
杨金忠 蔡树英 黄冠华 叶自桐 多孔介质中水分及溶质运移的随机理论 北京 科学出版社
责任编辑 吕斌秀
上面第二个例子中 动量方法求解摄动方程大约 次 计算耗时 而谱方法仅求解了约 次 计
算耗时 很明显 基于 展开的谱方法的计算效率高于动量方法
饱和 非饱和土壤中的水流 模拟区域为宽
长
的剖面 剖分为
个单元 上边
界给定入渗率 ? 下边界为隔水边界 左右边界的下部分为定水头边界 水头分别为
和
上部分为无流量边界 其它参数为 饱和水力传导度的几何平均
法分解 把压力水头表示为多项式 通过摄动方法得到一系列关于水头展开式的偏微分方程 用有限差分法
进行求解 获得了压力水头的随机描述 并计算其均值和方差 应用本文的随机模型研究了二维非饱和以及饱和
非饱和介质流动的实例 结果与动量方法的计算结果一致 而且计算效率高于传统的动量方法
关键词 随机模型 饱和 非饱和流
均水头相同 图
两者计算的水头标准差存在一点偏差 图
计算水头标准差 谱方法只需
要展开前 项就可获得与动量方法近似的结果 项以后的展开项对结果没有太大影响
将上边界的入渗率提高到 ? 同时增大 和
的变异性 取
计
算结果如图 所示 可以看到 计算水头标准差需要 展开前 项 在该例中 谱方法和动量方法的
计算结果也有一点偏差 通过模拟 可以发现水头标准差在离开下边界后快速增加 到达重力占优区后
项都需要求解 只需
求解关于
项的摄动方程 这些性质使得数值计算具有较高的效率
实例分析
基于ANSYS的复杂地质条件下大型渗流场计算
ＡＤ Ｐ Ｌ 以下功 能 ：
（） １ 参数 ；２ 表达 式 和 函数 ； ３ 分 支 和循 环 ； （） （）
（） ４ 重复 功能 和复 写 ； ５ 宏 ； ６ 用户 子程序 。 （） （） ＡＤ Ｐ Ｌ有很 多 优 点 ， 中 之 一 是 可 以利 用 Ａ － 其 Ｐ
实现 干燥 区和饱 和 区进 行 联 立 统一 求 解 ， 免 了 干 避
方 程相 似 、 始条 件 与边 界条 件 的相 似等 特点。 初 Ａ Ｙ ＮＳ Ｓ中渗 流 场 与 温 度 场 各 种 相 应 量 比较 见 表
１ 一 ［
燥 区渗透 系数取 为 ０时导致 渗透 矩 阵奇异 而无法进
ＡＰ ＤＬ模 块 来 自动 迭代 计 算 出渗流 边 界 和 自由面 。
温 度 场 与渗 流场 之 间具有 理 论基 础 相 似 、 微分
为 了便 于 与结构 计算 网格 相 对 应 ， 算 中采 取 了 固 计 定 网格 法 （ 不变 网格法 ） 。其 中干燥 区（ 自由面 以上 的 区域 ） 的渗 透 系 数 取 很低 的值 （ 文 在具 体 计 算 本 时 ， 燥 区渗透 系数 取 １０Ｘ１ ｃ ／ ） 从 而 可 以 干 ． ０ ｍ ｓ，
３ ．四 川 电力 进 出 口公 司 ， 四川 成 都 ６ ０ ６ ） １０ １
摘
要： 利用 ＡＮ ＹＳ温度场分析模块计算大岗山水 电站复杂地质条件下的大型渗流场 ， Ｓ 取得了较好的效果 。
关键词 ： ＮＳ Ｓ软件； Ａ Ｙ 渗流场 ； 计算方法 ； 应用软件 ； 岗山水 电站 大
中图 法 分 类 号 ： Ｖ ２ ． Ｔ ２３６ 文 献 标 识码 ： Ｂ 文章 编 号 ：０ ３—９０ （ ００）４—０ ２— １０ ８５ ２ １ ０ ４ ０ ４ ０
基于饱和-非饱和渗流的煤层注水数值模拟及应用
基于饱和-非饱和渗流的煤层注水数值模拟及应用吴兵;于振江;王紫薇;江泽泳【摘要】针对当前同忻煤矿井下粉尘浓度偏大、工作条件恶劣以至严重影响该矿安全生产的问题,根据非饱和土渗流理论,提出采用 comsol 有限元数值模拟软件进行煤层注水方案优化,最后进行现场应用及效果测定。
通过数值模拟得出同忻矿8103工作面的煤层湿润半径约为15 m,在保持4°倾角条件下,其湿润范围最大。
采用了新的注水方案后,工作面的粉尘浓度明显降低,其中采煤机中部和前滚筒下风侧10 m 处的降尘率分别达到了19.55％、30.62％。
%For larger concentration of dust in the current Tongxin Coal Mine and the poor working conditions,which seriously affects mine production safety issues,based on unsaturated flow theory,this article uses comsol finite element simulation software to optimize seam water injection program,and make the final determination of field application and effect. The results shows:the moist radius of 8103 working face in Tongxin Coal Mine is about 15 m,and in keeping + 4°inclination conditions,th e moist range max. After the introduction of a new water injection program,the dust concentration was signifi-cantly lower. The rate at which the dust before the shearer drum central and leeward side of the 10 m position reaches 19.55%、30.62% respectively.【期刊名称】《中国煤炭》【年(卷),期】2015(000)006【总页数】4页(P102-105)【关键词】不饱和土;非饱和流动;煤层注水;湿润半径【作者】吴兵;于振江;王紫薇;江泽泳【作者单位】中国矿业大学北京资源与安全工程学院，北京市海淀区，100083;中国矿业大学北京资源与安全工程学院，北京市海淀区，100083;中国矿业大学北京资源与安全工程学院，北京市海淀区，100083;中国矿业大学北京资源与安全工程学院，北京市海淀区，100083【正文语种】中文【中图分类】TD714.41煤层注水就是通过钻孔将压力水注入到尚未开采的煤层中,增加煤的水分,使煤体预先得到湿润,减少煤尘的生成和飞扬,是采煤工作面最有效的防尘措施,能够大幅度减少采煤工作面产尘量。
堤坝饱和-非饱和渗流场数值模拟及稳定性分析的开题报告
堤坝饱和-非饱和渗流场数值模拟及稳定性分析的开题报告一、选题背景与意义随着我国经济和社会的快速发展，对水利水电建设的需求也越来越大。
而水利工程中的堤坝工程是一个关键的领域，其稳定性问题一直是一个热点问题。
在堤坝设计和施工过程中，要考虑到多方面因素，其中之一是堤坝内的渗流问题。
如果渗流问题得不到解决，将会对整个工程的安全性产生极大的影响。
因此，对堤坝饱和-非饱和渗流场数值模拟及稳定性分析具有重要的理论和实际意义。
二、研究内容与方法本课题将以一座某水库工程的堤坝为例，开展以下研究内容：（1）对该堤坝进行渗透试验，获取堤坝内的水文地质参数；（2）建立该堤坝的饱和-非饱和渗流场数值模型；（3）在ANSYS等软件中进行数值模拟，并进行模型验证；（4）对模型的稳定性进行分析。
具体方法包括：（1）对堤坝结构、土质组成、水文地质条件进行调研和实地勘察；（2）结合试验数据建立堤坝的饱和-非饱和数值模型；（3）利用有限元方法在ANSYS等软件中进行模拟计算，并比对实验数据以验证模型的准确性；（4）通过分析模型结果，预测堤坝在满足设计要求的前提下的稳定性。
三、研究意义本研究对于提高水利工程建设质量、优化设计方案、保障工程安全等方面有着积极的促进作用。
具体如下：（1）促进水利工程领域渗透试验方法的改进与完善；（2）提高堤坝饱和-非饱和渗流场数值模拟技术水平，并在实际工程中得到应用；（3）为堤坝设计和施工提供科学依据，减轻水利工程的风险；（4）为今后更深入地研究水利工程的渗流问题提供理论基础和经验参考。
四、研究进度安排本研究将于XX年XX月开始，具体进度安排如下：第一阶段（XX年XX月-XX年XX月）：调研和实地勘察，获取实验试验数据和工程参数。
第二阶段（XX年XX月-XX年XX月）：分析和处理数据，建立数值模型。
第三阶段（XX年XX月-XX年XX月）：进行数值模拟计算，并与实验数据进行比对。
第四阶段（XX年XX月-XX年XX月）：对模型结果进行分析和总结，撰写研究报告。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
Ss为储水率 , 其值较小 , 一般忽略不计 ;α在饱和区取 1, 非饱和 区取 0;H0 为已知水头 ;q为 边界流 量 ;n为土 体空隙 比 ;cos(n,
x)、cos(n, y)、cos(n, z)为 边 界 面 外 法 线 方 向 余 弦 ;t为 时 间
变量 。
其中式 (1)为连 续性方 程 , 式 (2)～ 式 (4)为边界 条件 ,
图 3 非稳定渗流模拟流程
(下转第 80页 )
· 80·
人 民 黄 河 2008年
tem[ J] .OperRes, 1995(3):187 -197. [ 2] LoucksD P.Somecommentsonlineardecisionrulesand
2.2 渗透参数调整与逸出点位置的确定
渗透 参数是 未知的 , 需 要根据 计算结果 不断调 整 , 以往进 行调整时只使用上一 时步末的水头场分布 , 会产生过 大的累积 误差 , 因此考虑时步初 (即为上一时步末 )与时步末水头场分布 来确定时步内的渗透 参数更为合理 [ 8] 。
逸出 按如 下步骤确定逸出点位 置 [ 9 -10] :①第 1次计算 , 对于初始时步 , 将 可能的逸出面节点 全部作 为未 知点 , 对于 其余 时步 , 将时 步初 的逸出面节点作为 已知 点 , 经 过有 限元 计算求 出渗 流场分 布 。 ②将水面以上的正压 边界点全部按逸出点处 理 , 求出 渗流场分 布后 , 再对水面以上零 压边界 点进 行判断 , 若 流通 量大于 等于 0, 即
第 30卷第 4期 人 民 黄 河 Vol.30, No.4 2008 年 4月 YELLOW RIVER Apr., 2008
【水利水电工程 】
渗流计算是堤 坝设 计及 除险加 固中 的一项 重要 内容 。 以 往的渗流计算 , 一般忽 略非饱 和区 渗流的 影响 , 但 是土的 基质 势和重力势的作用 , 使土 体中孔 隙水 产生非 饱和 流动 [ 1] , 只考
虑饱和区渗流的方法 不能 有效地 解决 水位变 化和 降雨等 引起 的非稳定渗流问题 。
该方法在判断可能逸出面 节点流通量后 , 从逸出 面节点中 成批排除流通量小于 0的节点 , 没有考虑 施加在假定 逸出面上 的固定水 头对节点 流通量的 影响 , 现节合图 2 进行说明 。 图 2 为取水位降落 后部 分坝 坡面 , 曲 线 A为上 一 时步 求得 的 浸润 线 , 曲线 B为此时步待求的浸润线 , 节点 2、6则分别为上一时步 求 得的逸出点及此时步待求的逸出点 。求解时 , 对节点 2 ～ 9施 加 固定水头 (H =z), 节点 2 ～ 6所加固定水头将影响与其相邻 单元的水头值 , 水头与渗透参数的高度非 线性关系可 能导致节 点 2 ～ 7在此时步流通量均流向域内 , 而排除在逸出面之外 , 从 而使逸出点错误地收 敛于节点 8, 浸 润线相 应地 由实际 曲线 B 变为曲线 C。
2.3 非线性迭代的程序结构
非饱和土的渗流 问题是非线性问题 , 饱和度 、渗透系 数 、逸 出面范围等只能在计 算过程中迭代确定 。 迭代分 3层进 行 (见 图 3):最内层迭代确定渗流参数 , 迭代中 边界条 件及时 步长不 变 , 渗流参数采用单元孔隙水压力值所对 应的参数不 断进行调 整 , 直到前后两次计算 各单元 孔压 变化均 小于 某一 定值 ;次外 层迭代确定渗流逸 出面边 界 , 迭代 时步长 不变 , 根 据内层 的计 算结果调整逸出面 边界范 围 , 直至 迭代收 敛 ;最外 层迭代 调整 时间步长 , 只有相邻时 步间的 水头 差在允 许范 围之 内时 , 才进 行下一时步计算 , 否则 此步计 算结 果作废 , 更 改时 间步长 重新 计算 。
1 饱和 -非饱和渗流数学模型及有限元解法
Neuman将非饱和区 渗透 系数 看做 孔隙 水负 压的 函数 , 使
饱和区与非饱和区 渗流控 制方 程得到 统一 [ 5] 。 不 考虑土 骨架
变形及水压缩性时 , 非稳定饱和 -非饱和渗流数学模型如下 :
x[ kx(h) Hx] + y[ ky(h) Hy] + z[ kz(h) Hz] =
逐步积分 , 这些积分只要设定时间步长及 相应的边界 条件与初
始条件后 , 即可由 ANSYS自动实现 。
2 非稳定饱和 -非饱和渗流场数值模拟
2.1 时间步长的自动选择
非稳 定渗流 场数值 模拟中 , 选取 合适的 时间步 长 , 既能保
收稿日期 :2007-10-09 作者简介 :岳庆河 (1981— ), 男 , 山东济南人 , 硕士研究 生 , 主 要研究方向为水工结 构 。 通讯作者 :刘福胜 (1964— ), 教授 。 E-mail:liufsh@
Г 1
为已知水头边界 , Г2 为已知流量边界 , Г3 为逸出段边界 (渗流场
边界见图 1), 式 (5)为初始状态 , H0 即为渗流场的初始水头 。
图 1 渗流场边界条件示意
热传导基本微分 方程为
kx
2xT2 +ky
2yT2 +kz 2zT2 =Cρ
T t
(6)
比较渗流微分方 程与热传导微分方程 , 可见二者 控制方程
已有不少学者 应用 ANSYS进 行了 渗流 场数 值模 拟研 究 。 许玉景等 [ 2] 首次使用 ANSYS热分析模块 结合单元生 死技术进 行了稳定渗流模拟 , 蒙富强 [ 3] 对此 进行了 进一 步充 实 , 罗 启北 等 [ 4] 使用 ANSYS进行了导流洞有压渗 流的动态 模拟 。 为解决 存在浸润面的 非稳定 渗流 场模 拟问题 , 笔 者对 ANSYS进 行二 次开发 , 编制了非稳定饱 和 -非饱 和渗流 程序 , 并 通过实 例验 证了程序的正确性 。
(4)
H(x, y, z, t0) =H0(x, y, z, t0 )
(5)
式中 :H =h+z, H为总水头 ;h为压力水头 , 在饱和区为 正值 ,
非 饱和区为负值 ;kx、ky、kz分别为沿坐标 x、y、z的孔隙水主渗透 系 数 , 是压力水头 h的函数 ;C= θh为容水度 , θ为介质含水率 ;
余时步时 , 将上一时步 确定的 逸出 面节点 作为 已知 点 , 计 算出 渗流场分布 。 ②不 改变 边界条 件 , 对渗 透参数 进行 循环迭 代 , 直到达到足够精度 后 , 再判断 逸出 点流通 量 , 若流 通量大 于等 于零 , 则该点为 逸出 点 ;反 之 , 在 可能 逸出 面 节点 中删 除 此节 点 。 重复这 2个步骤 , 直到找出流通量大 于等于 0的 可能逸出 点为止 。
基于 ANSYS模拟非稳定饱和 -非饱和渗流场
岳庆河 1 , 刘福胜1 , 刘佩玺2
(1.山东农业大学 水利土木工程学院 , 山东 泰安 271018;2.济宁市水利局 , 山东 济宁 272019)
摘 要 :堤坝非稳定渗流数值模拟中 , 浸润面 、逸出点事先未知 , 需要迭代计算 , 渗透参数需根据 计算结果不断调整 , 都增 加了数值模拟的难度 。 选用 Neuman饱和 -非饱和渗 流数学 模型进行 非稳定 渗流场 数值模 拟 , 分析了数 值模拟 中逸出 面一般求解方法存在的 问题并加以改进 , 使用 ANSYS开发出渗流数值模拟程序 , 经实例验证表明该法是实用有效的 。 关 键 词 :饱和 -非饱和 ;渗流 ;逸出面 ;ANSYS 中图分类号 :TV139.14 文献标识码 :A 文章编号 :1000-1379(2008)04-0076-02
-kx Hxcos(n, x)-ky Hycos(n, y)-kz Hzcos(n, z)≥ 0 说明该节点流向域外 , 则该点为逸出面节 点 , 否 则 , 该 点不是逸 出面节点 。 ③将前 一次 确 定的 逸 出 面节 点 作为 已 知点 处 理 。 重复这 3个步骤 , 直到计算收敛 , 再进行下一时步 计算 。
的形式完全相 同 , 这就为 应用 ANSYS热分 析模 块进 行渗 流场 分析奠定了基础 [ 2 -4] 。 对二维渗流微分 方程应用 变分法 , 并对
时间取隐式差分 , 可得相应的有限元格式 。 对于存在 时 间变 量 的瞬 态 热 分析 问 题的 求 解 [ 6] , ANSYS
在空间域上按有限元 方法求逐步积分 , 在 时间域上用 差分方法
[ C(h)+αSs]
H t
(1)
H(x, y, z, t) =H1 (x, y, z, t), 在 Г1 上
(2)
kx Hxcos(n, x)+ky Hycos(n, y)+kz Hzcos(n, z)=q, 在 Г2 上
(3)
H(x, y, z, t) =z(x, y, z, t), t> 0, 在 Г3 上
图 2 逸出面调整示意 笔者对此加以改 近 , 只判 断可能 逸出 点 (也就 是逸出 面上 最高位置节点 )的流通量 , 如果流向域外则 计算收敛 , 否则从可 能逸出面节点中删除此节点 , 进行下一 次循环 。 这样 在每次循 环中只排除一 个可 能 逸出 面节 点 , 避 免了 上 述错 误情 况 的发 生 。 具体 做法为 :①第 1次 计算 , 对于初 始时步 将可能 的逸出 面节点全部作为已 知水头 节点 (H =z), 求 出渗流 场分布 , 其
chanceconstraints[ J] .WaterResour, 1970, 6(2):668 671. [ 3] 赵鸣雁 , 程春田 , 李刚 .水 库群系 统优化 调度新进 展 [ J] . 水文 , 2005, 25(6):18 -23. [ 4] AhmedI.Onthedeterminationofmultireservoiroperation policyunderuncertainty[ D] .TucsonArizona:TheuniversityofArizona, 2001. [ 5] 谭维炎 , 黄守信 .应用随机动态规 划进行水电站水库的优 化调度 [ J] .水利学报 , 1982(7):1 -7. [ 6] 李爱玲 .梯级水电 站水库 群兴利随 机优化 调度数 学模型 与方法研究 [ J] .水利学报 , 1998(5):71 -74. [ 7] DorfmanR.Themulti-structureapproachindesignofwaterresourcesystems[ M] .HarvardUniversityPress, 1962. [ 8] WindsorJS.Optimizationmodelforreservoirfloodcontrol [ J] .WaterResour, 1973, 9(5):1103 -1114. [ 9] BarrosM, TsaiF, YangS, etal.OptimizationofLargescalehydropowersystem operations[ J] .JournalofWater Resour, 2003, 129(3):178 -188. [ 10] PengCS, BurasN.Practicalestimationofinflowintomultireservoirsystem[ J] .JournalofWaterResour, 2000, 126 (5):35 -40. [ 11] 黄强 .黄河干流水库联合实施调 度及智能决策支持系统 研究 [ D] .西安 :西安理 工大学 , 1995.