用Ansys分析有浸润线的土石坝平面渗流问题
anays渗流
第2章水工流体力学问题的ANSYS模拟ANSYS软件具有专门针对流体的计算模块，可以模拟理想流体的流动、稳定和非稳定渗流。
本章主要介绍利用ANSYS进行简单的渗流分析和流体分析。
2.1 水工渗流场模拟需要指出的是，在ANSYS中并不存在专门的渗流分析模块，但由于渗流场与温度场的有限元计算公式相同，因此可以利用ANSYS中的热分析模块进行分析，只需相应参数对应采取即可。
鉴于问题的复杂性，这里仅用一个重力坝的例子进行分析。
2.1.1 数值模拟对象考虑混凝土大坝下水的渗流。
假定坝下土壤的渗流系数K约为每天15米，即K＝15m/天，试确定该土壤的渗流速度分布。
本问题为一个稳态渗流问题，可以利用ANSYS的稳态热传导进行分析。
2.1.2 有限元模型本问题中，考虑混凝土坝体不透水，渗流主要发生在坝基部分，因此，模型仅涉及坝基，具体建模过程如下：1．模型的定义启动ANSYS，设置好文件夹及文件名；2．设置分析类型点击Preferences，在弹出的对话框中选择Thermal，选择热分析，点击OK；随后进入前处理模块：Main Menu>PreProcessor；3．设置单元类型对于二维渗流（热）分析，采用Plane55单元。
Element Type>Add/Edit/Delete…，弹出对话框：4．定义材料参数稳态渗流分析中，主要设置材料的渗透系数（本问题中，坝基材料的渗透系数为K ＝15m/天，其对应热分析中的热传导系数）。
设置方法如下：Material Props>Material Models …，弹出对话框中：5．建立数学模型本问题中，通过创建关键点，再由关键点直接生成二维（坝基）面。
●在命令行中，分别输入关键点命令，生成关键点：K,1,0,0；K,2,5,0；K,3,9,0；K,4,16,0；K,5,16,4；K,6,16,5；K,7,9,5；K,8,9,4；K,9,5,4；K,10,5,5；K,11,0,5；K,12,0,4●连接关键点生成面Modelling>Creat>Areas>Arbitary>Throuth KPs…，弹出点选择对话框：8，9）、（3，4，5，8）、（8，5，6，7）、（12、9，10，11），生成其它四个面。
ANSYS软件在土坝渗流稳定计算中的应用_许玉景
文章编号:0559-9342(2003)04-0069-03ANSYS软件在土坝渗流稳定计算中的应用许玉景1,孙克俐2,黄福才1(1.天津市水利勘测设计院,天津 300204;2.天津大学建工学院,天津 300072)关键词:ANSYS软件;温度场;渗流场;死活单元;浸润线;渗流稳定计算;土坝摘 要:根据基本方程及定解条件的比较分析,将ANSYS软件的温度场分析功能应用于渗流场的分析,并采用死活单元技术,通过迭代算法计算自由水面位置(浸润线),解决了土坝渗流稳定问题的求解。
该方法可以解决复杂边界、多种介质的渗流问题,为工程应用提供极大的便利。
Applicatio n of ANSYS in the Earth dam seepage flow stability computatio nXU Yu-jing1,SUN Ke-li2,HUANG Fu-cai1(1.Tianjin Water Conservancy Survey&Design Institute,Tianjin300204;2.Civil Engineering School of the Tianjin University,Tianjin300072)Key Words:ANSYS software,temperature field,seepage flow field,element birth or death,saturation line,seepage flow stability calculation,earth damA bstract:According to comparison of the fundamental equations and boundary conditions,the ANSYS function of temperaturefield analysis is applied to that of seepage flow field,the element birth or death technique with overlap method is adopted to calculate the free water surface site(saturation line),in order to solve the problem of Earth dam seepage flow stability.This method can solve the seepage flow problem of complex boundaries and mediums.中图分类号:TV139.14;TV641 文献标识码:B渗流会对土坝稳定产生严重的危害。
浅谈基于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二次开发的心墙堆石坝渗流分析
确定 复合单 元 的等效渗 透矩 阵 。该方 法需要 判 断 自由 面与 复合单 元 的交截情 况 ， 对于 交截情 况较 多 的问题 ，
程序 处理上 有 些不便 。 本文基 于 Ａ Ｎ Ｓ Ｙ Ｓ热分 析 模 块 进 行 二 次 开发 应 用
前人 基于 Ａ Ｎ Ｓ Ｙ Ｓ热 分析 模 块 进 行 二 次 开 发 用 于 渗 流稳定 分析 开展 了很多 研究 。许玉 景运 用死 活单元
第４ ４卷 第 １ ９期
２ ０ １ ３年 １ ０月
人 民 长 江
Ｙａ ｎｇ ｔ ｚ ｅ Ｒｉ ｖ ｅ ｒ
Ｖｏ １ ． ４ ４． ＮＯ ．１ ９
０ｅ ｔ ．． ２ ０１ ３
文章编号 ： １ ０ ０ １ — ４ １ ７ ９ （ ２ ０ １ ３ ） １ ９— ０ ０ １ ２— ０ ４
关 键 词： 浸润线 ； Ａ Ｎ Ｓ Ｙ Ｓ二 次 开 发 ； 渗 流 稳 定 ；心墙 堆石 坝 文献标志码 ： Ａ
中 图 法分 类 号 ： Ｔ Ｖ ６ ４ １ ． ４ １
据不 完全 统计 ， 我 国已建 大坝有 ８ ． ７万 余 座 ， ９ ５ ％ 以上为 土石坝 ， 且在 我 国高坝 中所 占的 比重 呈 逐 渐增
的权 系数 ， 然 后按 加权 平均确 定单 元 的等效渗 透 系数 ， 从 而计 算 出相 关 渗 流 物 理 量 。这 种 方 法 的 核 心 是
坝 渗流稳 定进 行分 析和评 价 。 Ａ Ｎ Ｓ Ｙ Ｓ作为一 种 大型通 用 数值 分 析 软件 ， 拥 有 强 大 的分 析及 前后处 理 功 能 ， 但 在水 工 结 构 方 面 的实 际 运用 中缺乏 处理具 有 自由面 的无 压 渗 流 问题 的能 力 ，
基于ANSYS的大坝渗流分析研究
总体而言，ANSYS在大坝渗流分析中具有重要的应用价值和潜力。未来可以 进一步探索ANSYS在解决实际工程问题中的其他应用，为水利工程领域的科学研 究和技术发展提供更全面的支持。
谢谢观看
然而，ANSYS在应用于大坝渗流分析时也存在一定的局限性。例如，对于复 杂的三维模型，网格划分和计算量可能会变得非常大，需要较高的计算资源和时 间成本。此外，ANSYS的数值计算结果也会受到许多因素的影响，如模型的简化 程度、参数设置是否合理等。因此，在使用ANSYS进行大坝渗流分析时，需要充 分考虑这些因素，并进行必要的实验验证和对比分析，以确保分析结果的准确性 和可靠性。
结论
本次演示介绍了基于ANSYS的大坝渗流分析方法及其应用。通过文献综述， 总结了ANSYS在水利工程领域的应用及研究现状，指出了现有研究的不足之处。 在此基础上，阐述了ANSYS进行大坝渗流分析的基本原理、方法和流程，并通过 实例分析展示了ANSYS在解决实际工程问题中的应用。最后，对实例分析的结果 进行了讨论和评价，总结了ANSYS在大坝渗流分析中的优势和适用性，并指出了 其局限性。
结果表明，ANSYS作为一种强大的数值计算和可视化工具，在应用于大坝渗 流分析时能够有效地模拟渗流场的分布情况和影响因素。同时，ANSYS还具有广 泛的适用性和灵活性，可以针对不同的工程实际问题进行模型的灵活调整和优化 计算。然而，对于复杂的三维模型和特定的工程问题，仍需充分考虑ANSYS的局 限性，并进行必要的实验验证和对比分析。
实例分析
以某实际大坝渗流问题为例，运用ANSYS进行渗流分析。首先，根据大坝的 结构特征和实际运行环境，建立大坝及周围区域的几何模型；然后，根据工程实 际情况，将模型划分为适当的计算网格；接着，设定合理的边界条件和材料属性， 进行数值计算；最后，根据计算结果进行后处理，得到渗流场的分布情况。
基于ANSYS软件的堤防渗流分析
的渗透 系数 ， 坐标轴 方向与渗透主 方向一致 ； 为 渗流 区域 ； Ｓ 为已知水头值 的边 界 曲面 ； Ｓ ２ 为给定
流量边界曲面 ； Ｓ ， 为浸润面 ； Ｓ ． 为逸 出段 ； ｇ为边界
上 的单 位 面积 流 量 ，这 里 ｑ ＝ Ｏ表 示 为 无流 量交 换
由于有 限元法能够建立各 种复杂 的几 何模 型 、 具 有强大的前 、 后处理 功能等特点 ， 其越来越 多的被
堤防前后均取 ２ ０ ｍ， 沿堤防轴线方 向取 ３ ０ ｍ。模 拟 中依次建立管道 、 回填土体 （ 形状按开挖尺寸定
头及等势线云 图， 根据计算结果 ， 堤防 的最大渗透
比降为 ０ ． ４ ７ ６ ， 管道 穿越处筑堤材料 的允许 渗透 比 降为 ０ ． ５ Ｏ ． ８ ｏ ， 即发 生 ５ Ｏ 年 一遇洪水 时 ， 堤 防的
１ 计算原 理
对于各 向异性 、 非均质的连续介质 ， 服从 达西
定律 的稳定渗流问题 ， 可归结为下列定 解问题 ：
軎 鲁 ａ （ ｋ ， 等） ＋ 軎 誓） ＝ ０ 在 以上的单元 的渗透系数 ，来调 整这 部分单元 对总
，
单元渗透矩阵调整法的原理 上 ，通过 改变 自由面
【 摘 要】 某管道 穿越 西里尼 西 河的计 算 实例表 明， Ａ Ｎ Ｓ Ｙ Ｓ 软件 中的热分析模 块 ，可以 用来进
行堤防的渗流分析。对计算结果的分析发现 ， 回填土体的渗透系数， 可以作为衡量管道穿越 工
程对堤 防影响的 工程 质量控 制指标 。
【 关键词 】 ／ ￣ Ｎ Ｓ Ｙ Ｓ 堤 防； 渗 流分析 【 中图分类号 】 Ｔ Ｖ ２ ２ ３ ． ４
后 ｃ ０ ｎ （ ） ＋ 后 ｃ ｏ ｎ （ ｎ， ） ＋
ANSYS有限元生死单元技术在砼面板堆石坝渗流计算中的研究与应用
甘 肃 农 业 大 学 学 报 J OU RNAL O F GANSU A GRICUL TU RAL UN 单元技术在砼面板堆石坝 渗流计算中的研究与应用
吴建东 ,许健
(甘肃农业大学工学院 ,甘肃 兰州 730070)
WU J ian2do ng , XU J ian
(College of Engineering , Gansu Agricult ural U niversity , Lanzho u 730070 , China)
Abstract :Based o n t he basic t heory of seepage and co nt rol equatio ns , a new finite2element seepage cal2 culatio n mat hematical model of co ncrete face rockfill dam was p ropo sed. By utilizing t he similarit y of tem2 perat ure and seepage field and large2scale finite2element co mp utatio n soft ware AN S YS , ekill and ealive technology co uld be used to kill t he unit above t he infilt ratio n line and activate t he unit below t he infilt ra2 tio n line while impo sing correspo nding bo rder co nditio n , and finally o btained t he infilt ratio n line and t he seepage speed of t he dam t hro ugh simulatio n. The result s suggested t hat t he model accuracy meet t he p ro2 ject requirement s , at t he same time t he model was simple and t he result s were int uitive.
基于ANSYS的土石坝稳定渗流场的数值模拟
基于ANSYS的土石坝稳定渗流场的数值模拟一、本文概述随着水利工程的日益发展，土石坝作为一种重要的水利结构，其稳定性与安全性受到了广泛关注。
渗流是土石坝中普遍存在的物理现象，对坝体的稳定性产生深远影响。
因此，对土石坝稳定渗流场的深入研究和分析具有重要的工程实践意义。
本文旨在利用ANSYS这一强大的工程模拟软件，对土石坝的稳定渗流场进行数值模拟，以期更准确地理解渗流对土石坝稳定性的影响，并为土石坝的设计、施工和维护提供理论支持和实践指导。
本文将简要介绍土石坝及其渗流现象的基本概念，阐述稳定渗流场研究的重要性和必要性。
然后，详细介绍ANSYS软件在水利工程中的应用，以及其在土石坝稳定渗流场数值模拟中的优势。
接下来，本文将详细描述数值模拟的过程，包括模型的建立、边界条件的设定、计算参数的选择等。
通过对模拟结果的分析和讨论，揭示土石坝稳定渗流场的特征和规律，为土石坝的安全稳定运行提供理论支撑。
本文的研究不仅有助于深化对土石坝渗流规律的理解，也有助于提升水利工程的设计水平和施工质量，为保障水利工程的安全运行提供有力支持。
二、土石坝渗流基本理论土石坝是一种利用当地石料、土料或混合料，经过抛填、碾压等方法堆筑成的挡水建筑物。
在土石坝的运行过程中，渗流是一个不可忽视的物理过程，它关系到坝体的稳定性和安全性。
因此，对土石坝渗流的基本理论进行深入研究，对于保障坝体安全、优化坝体设计具有重要意义。
渗流是指液体在固体骨架中通过孔隙或裂隙流动的现象。
在土石坝中，渗流主要受到重力、孔隙水压力、坝体材料性质以及边界条件等因素的影响。
当库水通过坝体向下游渗流时，会形成一定的渗流场。
这个渗流场是一个三维的空间分布，其中包含了渗流速度、渗流压力、渗流量等多个物理量。
土石坝的渗流场分析通常采用达西定律来描述渗流速度与渗流压力梯度之间的关系。
达西定律表达式为：v = -k * (dP/dx)，其中v为渗流速度，k为渗透系数，dP/dx为渗流压力梯度。
基于ANSYS的土坝渗流分析及其应用
４７ ６
基 于 Ａ Ｓ Ｓ的土 坝 渗 流 分 析及 其 应 用 ＮＹ
许 尚杰 ，党发 宁 ，田威
（． １ 西安理工 大学 岩 土工程研究所 ， 陕西 西安 ７ ０ ４ ； ． １０ ８ ２ 山东省水利科学研究 院， 山东 济南 ２ ０ １ ） ５ ０ ３
摘要 ： 利用热传导和水渗流原理的相似性 ， 采用 Ａ Ｓ Ｓ软件提供 的热分析模块， ＮＹ 开发编制 了以总 水 头为 因变量 的 多孔 隙介 质 的饱 和． 非饱 和渗 流程 序 。该程序 可用来模 拟复 杂 渗流 场 ， 而且 能够检
西安理工 大学 学报 Ｊｕ ａ ｏ ｉａ ｎｖｒ ｔ ｏ ｅｈ ｏ ｇ （０ ８ ｏ．４Ｎ ． ｏｒ ｌ ｆ ’ｎＵ ｉ ｓ ｙ ｆ ｃ ｎｌ ｙ ２０ ）Ｖ １ ｏ４ ｎ Ｘ ｅ ｉ Ｔ ｏ ２
文章编 号：１０ －７ ０ ２ ０ ４０ ６ －５ ０ ６４ １ （０ ８ ０ －４ ７０
ｐｏ ｅ ｓ ｒ ｃ ｓ ．Ａｌｏ ｈ ｓｓ ｆ ｒ ａ ｅ ｕ ｅ ｏ ｅ ｔ ｂ ｉ ｈ ｄ ｌｏ ｌ ｙ ｓ ｐ － ｌ ｄ ｍ ｆＬ ｏ ｈ ｎ Ｒ ｓ ｔ ｓ ，ｔ ｉ ｏｔ ｅ ｃ ｎ ｂ ｓ ｄ ｔ ｓａ ｌ ｈ ｔ ｅ ｍｏ ｅ ｆｃａ ｌ ｅ ｗａｌ ａ ｏ ａ ｓ ａ ｅ ｅ — ｗａ ｓ ｏ
中图分 类号 ： Ｕ ３ Ｔ ４ 文献标 识码 ： Ａ
Ａｎ ｌ ｓ ｎ ｐ ｉａｉ ｎ ｏ ｍ ｅ ｐ ｇ ｓ ｄ ｏ ａｙ ｉ ａ ｄ Ａｐ ｌ ｔ ｆＤａ Ｓ ｅ ａ ｅＢａ ｅ ｎ ＡＮＳ ｓ ｃ ｏ ＹＳ
Ｘ ｈ ｎ －ｅ ， Ａ Ｇ Ｆ —ｉｇ ， Ｉ Ｎ Ｗｅ ＵＳａｇ ｉ ｊ Ｄ Ｎ ａｎｎ ＴＡ ｉ
查每 一迭代 过程 的计 算结 果并 查看 收敛过 程 。应 用该软件 建 立 了崂 山水 库粘 土斜墙 坝模 型并进行
基于ANSYS的土坝三维渗流场模拟
h z
）＝0
h│t=t0＝h（x，y，z）
温度场
温度 T
热流量 Qr 热传导系数 kr 热传导速度 qr
温度梯度 E
热传导定律：Qr＝AKr
dT dn
x（krx
T ）+ x
y（kry
T ）+ y
z（krz
T ）＝0 z
T│t=t0＝T（x，y，z）
h（M，t）│Γ＝φ（s M，t）
T（M，t）│Γ＝φ（s M，t）
2 渗流场计算理论及方法
由于渗流场与温度场有极大的相似性，基于此可利用
ANSYS 热分析模块来分析渗流问题。
用有限元方法求解渗流场时最核心、最复杂的问题便是
渗流自由面的确定。目前常用的方法有固定网格法和变动网
格法，这两种方法均可以用 ANSYS 的 APDL 参数化编程语言
来实现，变动网格法需要借助 ANSYS 的生死单元功能来实
现，而固定网格法计算相对简单实用，本文采用 Bathe(1979)年
提出的单元传导矩阵调整法来进行计算，将整个渗流域 R 划
分为两个子域 R1 和 R2。在自由面以上的子域 R1 中节点流速 为 0，令其渗透系数 k 降低至很小的值；在自由面以下子域 R2 中，流速大于 0，其渗透系数 k 仍保持不变，即：
1 引言
在水库除险加固工程中，渗流问题一直是关系到土石坝 安全的关键。目前渗流分析计算领域中应用比较广泛的有数 学解析法、电模拟法和数值分析法。数值分析方法主要有：有 限差分法、有限单元发和边界元法。其中，有限元法随着各种 大型有限元软件快速发展而得到了广泛的应用，常用的大型 通用有限元分析软件有 ANSYS、ADINA、ABAQUS 等，也有一 些专门针对渗流场的计算软件，如 Seep3D、Autobank、STSE、 Modflow 等。各种计算方法和软件各有优缺点，对于复杂边界 条件和复杂渗流介质的空间渗流问题很多平面渗流计算软件 都无法准确计算，本文提出利用 ANSYS 的热分析模块来分析 土石坝的三维渗流问题，并将 ANSYS 三维分析成果中的典型 断面和 Autobank 计算结果进行对比，检验计算结果的准确 性。本文的方法为解决复杂边界、多种介质、多物理场耦合的 三维渗流问题提供有效的方法。
基于ANSYS的土石坝渗流与稳定分析研究的开题报告
基于ANSYS的土石坝渗流与稳定分析研究的开题报告一、研究背景和意义土石坝是一种重要的水利工程结构，其安全稳定性直接关系到人们的生命财产安全和社会经济发展。
而渗流问题是土石坝安全稳定性研究的重要内容之一。
在土石坝工程设计、施工和运行过程中，渗流问题一直是困扰工程师的难题，如何在渗流对土石坝安全稳定性产生影响的情况下，保证土石坝的安全运行是当前亟需解决的问题。
本研究利用ANSYS软件，通过有限元数值模拟方法，研究土石坝内部的渗流分析及土石坝的稳定性分析，旨在探讨土石坝渗流及其对稳定性的影响规律，为土石坝设计、施工和运行提供技术支撑和依据。
二、研究内容和方法本研究的主要内容有两个方面：一是土石坝内部渗流模拟及分析；二是土石坝稳定性分析。
具体通过以下步骤实现：1. 确定研究对象：本研究以某一具体土石坝为研究对象，对其渗流分析及稳定性进行模拟和分析。
2. 建立土石坝模型：根据实际情况建立土石坝三维有限元模型，包括坝体、坝基、边坡等，考虑土、石材料的物理力学特性。
3. 渗流模拟：对建立的土石坝模型进行渗流模拟，通过ANSYS中的多孔介质渗流模型，对土石坝内部流场进行数值计算和分析。
4. 渗流分析：根据渗流模拟结果，分析产生渗流的原因，判断坝体、坝基是否产生渗漏现象，并分析渗漏现象的破坏机理。
5. 稳定性分析：根据建立的土石坝模型，通过ANSYS有限元分析软件对土石坝的稳定性进行数值计算，分析坝体的变形、破坏状况，确定安全系数，预测土石坝的破坏条件。
三、预期成果和意义本研究通过ANSYS软件，对土石坝内部的渗流分析及稳定性分析进行研究，预期取得以下成果：1. 对土石坝内部渗流模拟及分析技术的研究与应用，提高土石坝设计、施工和运行的水平，为工程师在实际工程中提供技术支撑和依据；2. 对土石坝安全稳定性分析方法的探究和应用，为土石坝的安全设计和管理提供科学依据，提高工程的安全性和经济效益；3. 深入了解土石坝渗流及其对稳定性的影响规律，为水工、环境等领域的科研人员提供参考，促进相关学科的发展。
ANSYS在土石坝防渗加固分析中的应用
ANSYS 在土石坝防渗加固分析中的应用摘要：为了评价多头小直径深搅桩防渗墙的防渗效果，通过有限元软件 ANSYS 热分析模块 APDL 参数化语言结合反演分析理论，分析了大坝加固前后的渗流特性，对比其浸润面、逸出点、渗流流量的变化。
计算结果表明 ANSYS 热分析能准确得模拟大坝渗流状态，多头小直径深搅防渗墙能有效降低大坝渗流量，加强坝体稳定性，该防渗加固方法可广泛推广应用。
关键词：多头小直径；深搅桩防渗墙；防渗加固；APDL ；反演分析1引言多头小直径深搅桩防渗墙技术在近几年得到了较大范围的应用，但对其防渗效果工程界还无准确的定论。
本文运用 ANSYS 热分析模块，分析某土石坝防渗前后渗流量的变化，为多头小直径深搅桩防渗墙技术的推广提供一定的理论依据。
2 ANSYS 分析渗流原理渗流本是复杂的三维应力场问题，但考虑到坝身长度远大于坝宽，可将其简化为二维问题加以研究。
假设坝体每层土体为各向同性并不可压缩，在稳定流作用下，渗流控制方程为：式中：、分别为、方向的渗流系数；为水头。
而热力学二维微分控制方程为：式中：、分别为、方向的热传导系数；T 为温度。
由（ 1）式和（ 2）式可以看出，渗流场和热力场的控制方程相同。
同样，在边界条件、坡降、流量计算等方面两者也具有高度的一致性。
由此，不难看出渗流场其实是热力场的一种特殊情况，通过ANSYS 热力学模块仿真模拟大坝渗流是合理的。
3有限元模型建立及渗流系数反演3.1 有限元模型建立为了分析多头小直径深搅桩防渗墙技术防渗效果，选取某水库经多头小直径深搅桩防渗墙技术防渗后的典型断面CS01 进行分析。
坝基土层分布如下：①层为粉质黏土，主要是人工填筑而成。
由于多次加高，碾压不实渗透严重；②层为粉质黏土，塑性较好，层内发现多处小孔，粉性局部略大。
上述两层为渗流发生主要区域；③层为粉质黏土夹粉土；④层为粉质黏土，其下部为完整基岩，可视作不透水层。
通过 ANSYS的APDL命令输入各关键点坐标，然后分别创建线和面建立大坝平面模型。
用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的土石坝渗流数值模拟研究
防 渗方 案 的科学性 与 合理性 ， 从 而 为土石 坝 工程 设计
０ 引言
十 石坝 是一种 常 见 的水 工建 筑物 ， ＿ 它 对地质 要 求
和 建设提 供 了一 定 的参考和 借 鉴 。
１ 基于 Ａ ＮＳ ＹＳ的渗流计算的基本原理
相 对较 低 、施 工简单 ，能就地 取材 ，造 价 相对 经济 ，
（ １ ．湖南省 水运 建设投 资集团有 限公 司，湖 南 长沙 ４ １ ０ ０ ０ ０ ；２ ．大连 海洋大学，辽 宁 大连 １ １ ６ ０ ２ ３）
摘 要 ：基 于有限元原理 ，采用 Ａ ＮＳ ＹＳ软件 中的热处理模 块 ，建立 了坝体数值 分析模 型，针 对湖南省 长沙县五 龙 山 土石坝的渗流过程进行 了数值模拟 ，得到 了在三种 工况下坝体加 固前后 的渗 流速度 、渗流量 、渗透坡 降、选 出点 高 程 、渗透比降等特征值 ，通过比较 可知 ：在 正常蓄水位 、设计 水位及校核 洪水位条件 下，采用 粘土防渗墙进行坝体 的加 固，其渗流速度 、单宽渗 流量、逸出点 高程 以及 逸出点渗透 比降均 小于坝体加 固前计算相应值 ，平 均渗透坡 降 高于加 固前计算相应值 。计算结果表 明采用 粘土防渗处理 方法进行坝体加 固，对保证 土石坝 的渗流稳定和坝体安全
土石 坝加 固前 后渗 流 的数值模 拟 结果 比较 ， 肯定 了该
基于ANSYS的大伙房水库土坝渗透系数反演
基于ANSYS的大伙房水库土坝渗透系数反演佟晓娜【摘要】根据温度场和渗流场的相似性,利用有限元分析软件ANSYS的热分析功能计算土石坝渗流问题.通过ANSYS的优化设计功能,建立了基于ANSYS的渗透系数反演模型,反演得到更符合大伙房水库大坝实际工况的渗透系数,利用该渗透系数进行渗流数值计算,通过对计算值与实测值的对比表明,有限元计算值与实测值相对误差较小,拟合结果能够满足精度要求.说明所求得的渗透系数比较合理,可以用于渗流数值计算.【期刊名称】《黑龙江水利科技》【年(卷),期】2012(040)012【总页数】3页(P47-49)【关键词】ANSYS;土坝;渗透系数;反演;有限元模型;误差较小【作者】佟晓娜【作者单位】辽宁省大伙房水库管理局,辽宁抚顺113007【正文语种】中文【中图分类】TV2141 概述由于大坝边界条件的复杂性和坝体、坝基渗流的非均匀性，用渗流基本方程计算成果难以与大坝实际运行情况相一致，因此有必要通过渗流监测资料反演其渗透系数。
即在渗流资料分析基础上，选取一段库水位相对稳定时间的大坝渗流监测资料，开展大坝稳定渗流反演分析。
本文根据温度场和渗流场的相似性，利用有限元分析软件ANSYS的热分析功能计算土石坝渗流问题，通过ANSYS的优化设计功能，建立了基于ANSYS的渗透系数反演模型，反演得到大坝稳定渗流后的心墙料的渗透系数。
2 工程概况大伙房水库是我国“一五”期间修建的重点水利工程，具有防洪、供水、灌溉、养鱼、发电等多种功能，是新中国成立后我国第一座自行设计、自行施工的大型水利枢纽工程。
水库主坝为碾压式黏土心墙坝，心墙顶部宽度为2.25 m，心墙底部宽度为15.0 m，心墙顶部高程为136.8 m。
黏土心墙底部为混凝土截水墙，底座宽1.0 m，墙底宽0.7 m，顶宽0.5 m，墙身高2.0 m。
坝顶宽8 m;坝顶高程139.8 m，防浪墙顶部高程141.0 m。
3 稳态温度场与稳定渗流场的相似性土石坝渗流基本微分控制方程如下:式中:H为总渗透水头;kx，ky，kz为以 x，y，z轴为主方向的渗透系数;μs=ρg(α+nβ)为单位贮水量或贮存率。
基于ANSYSAPDL语言的黄壁庄水库土坝渗流稳定分析(PDF X页)
最 臀」 景 臀」 母 哥」
二
,
相应 的边 界条件进行分析 元的死活
计 算精度
。
相应修 正边界条件后重新计算直到达到
由于 这种
。
“
定解 条件 为
住 一
。 。
,
死活
单 元技术的基本 对象
,
是单元
几
任
,
若 希望足够精度 和足 够光滑的浸润线
网
。‘ 。
格划分必 须足够密
卜
兀
单元 用
黄壁庄水库土 坝计算模型采用
,
随后于
,
世纪
年代逐 渐被
,
向渗流速度为 水量的函 数 所以在 流 问题
,
而 在浸润面 以上
,
引用到流体力学的领域 作为 二 向渗流问题考虑
土 体的压 缩性时
,
最早就是引用到稳定渗流 的分布为研 究对象 垂直剖面
二
,
一 般 可 假设为
。
即浸润线 以上 部分
,
一 般 以测压 管水头
不应参于计算
由于浸润线的位置 一 般是未知 的
,
、
城 市供水
,
、
养 鱼等综合 利
月开 始 兴
。 ,
集结构
。
热
、
流体
,
、
电磁
、
声学
用 的大型水利枢纽
,
总库容
亿 立 方米 黄壁庄 年
,
在
多年 的发展过程 中 是
,
不 断改
水库主 坝为水中填 土 均质坝 建
最 大坝高 长
, 。 ,
进 提高 功能不断增强
,
年 又 将主 坝 加高 至现坝 顶高程 坝顶宽
Ansys分析有浸润线的土石坝平面渗流问题
Ansys分析有浸润线的土石坝平面渗流问题土石坝渗流分析，采用非饱和土渗流参数，迭代计算浸润线，根据前次计算结果，不断修改单元的渗透系数和浸润逸出点位置，直到满足精度要求。
本算例的土石坝体型比较简单.采用非饱和渗流计算.即渗透系数为空隙压力的函数.首先建立一个数据文件PPPP.TXT，存储渗透系数函数关系，如下。
第一列为空隙压力值（水头M)，第二列为渗透系数指数，渗透系数等于10^A(M/D)。
-10.00 -4.0E+00-9.00 -3.6E+00-8.00 -3.2E+00-7.00 -2.8E+00-6.00 -2.4E+00-5.00 -2.0E+00-4.00 -1.6E+00-3.00 -1.2E+00-2.00 -8.0E-01-1.00 -4.0E-010.00 0.0E+00土坝顶宽4M，上下游坡比均为1:2，总高12M，底宽52M。
上游水深8M，下游无水。
FINI/TITLE, EARTHDAM SEEPAGE/FILNAME,SEEPAGE5/PLOPTS,DATE,0*DIM,TPRE,TABLE,11,1,1,PRESS,KKPE ! 定义水压与渗透系数的关系数组*TREAD,TPRE,PPPP,TXT ! 读入数组*DIM,NCON,ARRAY,4 ! 定义数组，用于存贮单元四个节点号/PREP7SMRT,OFFANTYPE,STATIC ! THERMAL ANALYSISET,1,PLANE55MP,KXX,1,1 ! 饱和状态下的渗透系数MP,KXX,2,1E-4 ! 完全干燥下的渗透系数，假设空隙水压力小于-10M时K,1,24,12K,2,24,0K,3,0,0K,4,28,12K,5,28,0K,6,52,0L,1,3L,3,2L,1,2L,4,5L,5,6L,4,6LESIZE,ALL,,,24A,1,3,2A,1,2,5,4A,4,5,6MSHK,2 ! MAPPED AREA MESH IF POSSIBLEMSHA,0,2D ! USING QUADSAMESH,ALL ! MESH AREASNUMMRG,NODE ! MERGE NODES AT BOTTOM OF CAISSON *GET,N_MAX,NODE,,NUM,MAX ! 获得最大节点号*GET,E_MAX,ELEM,,NUM,MAX ! 获得最大单元号*DIM,N_TEMP,ARRAY,N_MAX ! 定义节点温度变量-总水头*DIM,N_PRE,ARRAY,N_MAX ! 定义节点压力水头变量!定义上游面总水头值LSEL,S,LINE,,1NSLL,S,1NSEL,R,LOC,Y,0,8D,ALL,TEMP,8 !定义上游面总水头值!定义下游面总水头值LSEL,S,LINE,,6NSLL,S,1*GET,N_NUM2,NODE,,COUNT*DIM,N_NO2,ARRAY,N_NUM2II=0*DO,I,1,N_MAX*IF,NSEL(I),EQ,1,THEN ! 判断节点是否选中II=II+1N_NO2(II)=I ! 存储渗流可能逸出点节点编号*ENDIF*ENDDONSEL,R,LOC,Y,0,8 ! 第一次计算，假设浸润线逸出点在8M高位置,与上游同高*GET,N_NUM,NODE,,COUNT ! 获得渗流出口节点总数*DIM,N_NO,ARRAY,N_NUM ! 定义变量,存储渗流出口节点编号II=0*DO,I,1,N_MAX*IF,NSEL(I),EQ,1,THEN ! 判断节点是否选中II=II+1N_NO(II)=I ! 存储渗流出口节点编号*ENDIF*ENDDO*DO,I,1,N_NUMD,N_NO(I),TEMP,NY(N_NO(I)) ! 定义下游面总水头值*ENDDOALLSEL,ALLFINISH/SOLUSOLVEFINISH!第一次计算完毕!-------------------------------------------------------------------------!迭代计算CONUTT=20 ! 最大循环次数DD_HEAT=0.001 ! 前后两次计算，总水头最大允许计算差CHUK_ST=3 ! 出口边界条件重新设定的起始点CHUK_MAXY2=10E5 ! 临时变量，用于存储浸润线出口坐标*DO,COM_NUM,1,CONUTTDD_H=0/POST1SET,1*DO,I,1,N_MAX*IF,COM_NUM,GT,CHUK_ST+1,THENDD1=N_TEMP(I)*IF,ABS(DD1-TEMP(I)),GT,DD_H,THENDD_H=ABS(DD1-TEMP(I))*ENDIF*ENDIFN_TEMP(I)=TEMP(I) ! 计算节点温度（总水头）N_PRE(I)=N_TEMP(I)-NY(I) ! 计算节点压力，总水头-Y坐标*ENDDO*IF,COM_NUM,GT,CHUK_ST+1,THEN*IF,DD_H,LE,DD_HEAT,THEN*EXIT*ENDIF*ENDIF/PREP7! 重新给每个单元设定材料MATNUM=2*DO,I,1,E_MAX*DO,KK,1,4*GET,NCON(KK),ELEM,I,NODE,KK ! 获取单元四个节点编号*ENDDOTEMP_Y=(N_TEMP(NCON(1))+N_TEMP(NCON(2))+N_TEMP( NCON(3))+N_TEMP(NCON(4)))/4 !计算单元中心点平均温度RESS_T=TEMP_Y-CENTRY(I)*IF,PRESS_T,GT,0,THENRESS_T=0MPCHG,1,I*ELSEIF,PRESS_T,LT,-10,THENRESS_T=-10MPCHG,2,IMP,KXX,MATNUM+1,10**TPRE(PRESS_T)MPCHG,MATNUM+1,IMATNUM=MATNUM+1*ENDIF*ENDDO! 重新设定出口边界条件*IF,CONUTT,GT,CHUK_ST,THEN !前CHUK_ST次采用原边界条件LSEL,S,LINE,,6NSLL,S,1DDELE,ALL,TEMP ! 删除原边界条件II=0CHUK_MAXY=0*DO,JJ,1,N_NUM2*IF,N_TEMP(N_NO2(JJ)),GE,NY(N_NO2(JJ)),THEND,N_NO2(JJ),TEMP,NY(N_NO2(JJ)) ! 总水头=Y坐标*IF,NY(N_NO2(JJ)),GT,CHUK_MAXY,THENCHUK_MAXY=NY(N_NO2(JJ))*ENDIF*ENDIF*ENDDO*IF,CHUK_MAXY2,NE,CHUK_MAXY,THEN ! 判断前后两次计算的浸润线出口位置是否相同NSEL,R,LOC,Y,CHUK_MAXY ! 选择最高节点*IF,CHUK_MAXY,GT,0,THENDDELE,ALL,TEMP ! 删除出口最高节点边界条件*ENDIFCHUK_MAXY2=CHUK_MAXY*ENDIFALLSEL,ALLFINI/SOLUSOLVEFINISH*ENDDOSAVE!迭代计算完毕，进入后处理FINISH/POST1/CLABEL,,1/EDGE,,0/CONTOUR,,8,0,1,8PLNSOL,TEMP ! 显示总水头云图PLVECT,TF, , , ,VECT,ELEM,ON,0PLVECT,TF, , , ,VECT,NODE,ON,0LSEL,S,LINE,,6NSLL,S,1PRRSOL,HEAT ! PRINT FLOWRATE THROUGH SOIL FSUM,HEAT ! 计算渗流量*GET,Q_DAY,FSUM,0,ITEM,HEATALLSEL,ALLSAVE*DO,I,1,N_MAXN_TEMP(I)=TEMP(I) ! 计算节点总水头（温度）N_PRE(I)=N_TEMP(I)-NY(I) ! 计算节点压力，总水头-Y坐标DNSOL,I,TEMP,,N_PRE(I) ! 将压力水头值复制到节点*ENDDOPLNSOL,TEMP ! 显示压力水头云图FINI渗流量就是热流量，通过选择边界上的节点，各节点热流量的总和就是流量。
ANSYS在土石坝防渗加固分析中的应用
ANSYS在土石坝防渗加固分析中的应用土石坝是一种常见的水利工程结构，用于蓄水或防洪目的。
土石坝防渗加固是土石坝设计中非常重要的一环，其主要目的是防止坝体内的水从坝体内部泄漏出去，影响坝体的稳定性和使用寿命。
在土石坝防渗加固设计中，ANSYS软件可以发挥重要的作用。
1.坝体内水压分析：土石坝在蓄水或雨水侵蚀的情况下，坝体内部会产生水压力。
通过ANSYS软件可以对坝体内部的水压力进行分析，包括水压的大小、分布情况等。
这些信息对于设计合适的防渗加固措施非常重要。
2.渗流路径分析：土石坝防渗加固的设计需要准确地预测坝体内水的渗流路径。
通过ANSYS软件可以建立数值模型，模拟坝体内水通过不同部位的渗流路径，找出可能的渗漏点和路径，为加固设计提供参考。
3.加固材料力学性能分析：土石坝防渗加固通常涉及到各种加固材料，如混凝土、聚合物和土工合成材料等。
通过ANSYS软件可以对这些材料的力学性能进行分析，包括强度、刚度、变形等，为加固设计提供依据。
4.加固结构设计优化：在土石坝防渗加固设计过程中，需要确定合适的加固结构形式和参数。
通过ANSYS软件可以建立不同加固结构的有限元模型，进行受力分析和设计优化，找出最优的加固方案。
5.加固效果评估：在加固完成后，需要对加固效果进行评估。
通过ANSYS软件可以对加固后的土石坝进行受力分析和模拟，评估加固效果，判断是否满足设计要求。
综上所述，ANSYS在土石坝防渗加固设计中的应用具有重要意义。
通过ANSYS软件的分析和模拟，可以帮助工程师更准确地了解土石坝的受力和渗流情况，设计出更科学、安全、经济的防渗加固方案。
在实际工程应用中，工程师们可以借助ANSYS软件的强大功能，提高土石坝防渗加固设计的效率和可靠性，确保土石坝的安全运行和长期稳定性。
ANSYS在隧道渗流计算中的应用_郝晨琦
开 发 应 用ANSYS在隧道渗流计算中的应用郝晨琦１刘福胜２（１．深圳高速工程顾问有限公司，广东 深圳 ５１８０３４；２．华南理工大学，广东 广州 ５１０６４１） 摘 要：针对无压稳定渗流场，根据温度场与渗流场在理论基础、微分方程和边界条件方面具有极大的相似性，可以利 用ＡＮＳＹＳ热分析功能对其进行模拟计算。
以方斗山隧道为实例，运用ＡＮＳＹＳ参数化语言编程，通过程序迭代计算，得出了 渗流自由面，整个模型和衬砌的渗流量随注浆圈渗透系数的变化规律，并证实围岩注浆圈可使衬砌渗流量明显减少 ，对 隧道的设计施工具有一定的参考意义。
关键词：ＡＮＳＹＳ；渗流场；渗流量；浸润面；注浆圈 ＤＯＩ：１０．３９６９／ｊ．ｉｓｓｎ．１６７１－６３９６．２０１０．１３．０１８ Application of ANSYS in Tunnel Seepage Calculation １ ２ ＨＡＯ Ｃｈｅｎ－ｑｉ ，ＬＩＵ Ｆｕ－ｓｈｅｎｇ （１．ＳｈｅｎＺｈｅｎ Ｅｘｐｒｅｓｓｗａｙ Ｅｎｇｉｎｅｅｒｉｎｇ Ｃｏｎｓｕｌｔａｎｔｓ Ｃｏ．，Ｌｔｄ，ＳｈｅｎＺｈｅｎ，Ｇｕａｎｇｄｏｎｇ ５１８０３４；２．Ｓｏｕｔｈ Ｃｈｉｎａ Ｕｎｉｖｅｒｓｉｔｙ ｏｆ Ｔｅｃｈｎｏｌｏｇｙ，ＧｕａｎｇＺｈｏｕ，Ｇｕａｎｇｄｏｎｇ ５１０６４１） Abstract:Ａｃｃｏｒｄｉｎｇ ｔｏ ｔｈｅ ｇｒｅａｔ ｓｉｍｉｌａｒｉｔｙ ｂｅｔｗｅｅｎ ｔｅｍｐｅｒａｔｕｒｅ ｆｉｅｌｄ ａｎｄ ｓｅｅｐａｇｅ ｆｉｅｌｄ ｉｎ ｔｈｅ ｔｈｅｏｒｅｔｉｃａｌ ｆｏｕｎｄａｔｉｏｎ，ｄｉｆｆｅｒｅｎｔｉａｌ ｅｑｕａｔｉｏｎｓ ａｎｄ ｂｏｕｎｄａｒｙ ｃｏｎｄｉｔｉｏｎｓ，ｔｈｅ ｔｈｅｒｍａｌ ａｎａｌｙｓｉｓ ｆｕｎｃｔｉｏｎ ｏｆ ＡＮＳＹＳ ｃｏｕｌｄ ｂｅ ｕｓｅｄ ｔｏ ｓｉｍｕｌａｔｉｏｎ ｔｈｅ ｎｏｎ－ｐｒｅｓｓｕｒｅ ｓｔａｂｌｅ ｓｅｅｐａｇｅ．Ｔｈｅ ｓａｔｕｒａｔｅｄ ｓｕｒｆａｃｅ ａｎｄ ｔｈｅ ｖａｒｉａｔｉｏｎ ｌａｗ ｏｆ ｗｈｏｌｅ ｍｏｄｅｌ ａｎｄ ｔｈｅ ｌｉｎｉｎｇ ｓｅｅｐａｇｅ ｆｌｏｗ ｗｉｔｈ ｄｉｆｆｅｒｅｎｔ ｇｒｏｕｔｉｎｇ ｃｉｒｃｌｅ ｐｅｒｍｅａｂｉｌｉｔｙ ｃｏｅｆｆｉｃｉｅｎｔ ｉｓ ｏｂｔａｉｎｅｄ ｂｙ ｔａｋｅ Ｆａｎｇｄｏｕｓｈａｎ ｔｕｎｎｅｌ ａｓ ａｎ ｅｘａｍｐｌｅ，ｕｎｄｅｒ ｕｓｉｎｇ ＡＮＳＹＳ ｐａｒａｍｅｔｅｒ ｌａｎｇｕａｇｅ．Ｉｔ ｐｒｏｖｅｓ ｔｈａｔ ｔｈｅ ｇｒｏｕｔｉｎｇ ｃｉｒｃｌｅ ｃａｎ ｏｂｖｉｏｕｓｌｙ ｄｅｃｒｅａｓｅ ｔｈｅ ｓｅｅｐａｇｅ ｆｌｏｗ ｔｈｒｏｕｇｈ ｔｈｅ ｌｉｎｉｎｇ，ａｎｄ ｈａｓ ａ ｃｅｒｔａｉｎ ｒｅｆｅｒｅｎｃｅ ｖａｌｕｅ ｔｏ ｔｈｅ ｄｅｓｉｇｎ ａｎｄ ｃｏｎｓｔｒｕｃｔｉｏｎ ｏｆ ｔｕｎｎｅｌ． Key words: ＡＮＳＹＳ；Ｓｅｅｐａｇｅ ｆｉｅｌｄ；Ｓｅｅｐａｇｅ ｆｌｏｗ；Ｓａｔｕｒａｔｅｄ ｓｕｒｆａｃｅ；Ｇｒｏｕｔｉｎｇ ｃｉｒｃｌｅ１引言边界条件的相似性，为利用ＡＮＳＹＳ热分析模块来求解渗流问 题提供了依据。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
用Ansys分析有浸润线的土石坝平面渗流问题
土石坝渗流分析，采用非饱和土渗流参数，迭代计算浸润线，根据前次计算结果，不断修改单元的渗透系数和浸润线出口位置，直到满足精度要求。
本算例的土石坝体型比较简单.采用非饱和渗流计算.即渗透系数为空隙压力的函数.首先建立一个数据文件PPPP.TXT，存储渗透系数函数关系，如下。
第一列为空隙压力值（水头M)，第二列为渗透系数指数，渗透系数等于10^A(M/D)。
-10.00 -4.0E+00
-9.00 -3.6E+00
-8.00 -3.2E+00
-7.00 -2.8E+00
-6.00 -2.4E+00
-5.00 -2.0E+00
-4.00 -1.6E+00
-3.00 -1.2E+00
-2.00 -8.0E-01
-1.00 -4.0E-01
0.00 0.0E+00
土坝顶宽4M，上下游坡比均为1:2，总高12M，底宽52M。
上游水深8M，下游无水。
FINI
/TITLE, EARTHDAM SEEPAGE
/FILNAME,SEEPAGE5
/PLOPTS,DATE,0
*DIM,TPRE,TABLE,11,1,1,PRESS,KKPE ! 定义水压与渗透系数的关系数组
*TREAD,TPRE,PPPP,TXT ! 读入数组
*DIM,NCON,ARRAY,4 ! 定义数组，用于存贮单元四个节点号
/PREP7
SMRT,OFF
ANTYPE,STATIC ! THERMAL ANALYSIS
ET,1,PLANE55
MP,KXX,1,1 ! 饱和状态下的渗透系数
MP,KXX,2,1E-4 ! 完全干燥下的渗透系数，假设空隙水压力小于-10M时
K,1,24,12
K,2,24,0
K,3,0,0
K,4,28,12
K,5,28,0
K,6,52,0
L,1,3
L,3,2
L,1,2
L,4,5
L,5,6
L,4,6
LESIZE,ALL,,,24
A,1,3,2
A,1,2,5,4
A,4,5,6
MSHK,2 ! MAPPED AREA MESH IF POSSIBLE
MSHA,0,2D ! USING QUADS
AMESH,ALL ! MESH AREAS
NUMMRG,NODE ! MERGE NODES AT BOTTOM OF CAISSON
*GET,N_MAX,NODE,,NUM,MAX ! 获得最大节点号
*GET,E_MAX,ELEM,,NUM,MAX ! 获得最大单元号
*DIM,N_TEMP,ARRAY,N_MAX ! 定义节点温度变量-总水头
*DIM,N_PRE,ARRAY,N_MAX ! 定义节点压力水头变量
!定义上游面总水头值
LSEL,S,LINE,,1
NSLL,S,1
NSEL,R,LOC,Y,0,8
D,ALL,TEMP,8 !定义上游面总水头值
!定义下游面总水头值
LSEL,S,LINE,,6
NSLL,S,1
*GET,N_NUM2,NODE,,COUNT
*DIM,N_NO2,ARRAY,N_NUM2
II=0
*DO,I,1,N_MAX
*IF,NSEL(I),EQ,1,THEN ! 判断节点是否选中
II=II+1
N_NO2(II)=I ! 存储渗流可能出口节点编号
*ENDIF
*ENDDO
NSEL,R,LOC,Y,0,8 ! 第一次计算，仅假设浸润线出口在8M高位置,与上游同高*GET,N_NUM,NODE,,COUNT ! 获得渗流出口节点总数
*DIM,N_NO,ARRAY,N_NUM ! 定义变量,存储渗流出口节点编号
II=0
*DO,I,1,N_MAX
*IF,NSEL(I),EQ,1,THEN ! 判断节点是否选中
II=II+1
N_NO(II)=I ! 存储渗流出口节点编号
*ENDIF
*ENDDO
*DO,I,1,N_NUM
D,N_NO(I),TEMP,NY(N_NO(I)) ! 定义下游面总水头值
*ENDDO
ALLSEL,ALL
FINISH
/SOLU
SOLVE
FINISH
!!!第一次计算完毕
!-------------------------------------------------------------------------
!迭代计算
CONUTT=20 ! 最大循环次数
DD_HEAT=0.001 ! 前后两次计算，总水头最大允许计算差
CHUK_ST=3 ! 出口边界条件重新设定的起始点
CHUK_MAXY2=10E5 ! 临时变量，用于存储浸润线出口坐标
*DO,COM_NUM,1,CONUTT
DD_H=0
/POST1
SET,1
*DO,I,1,N_MAX
*IF,COM_NUM,GT,CHUK_ST+1,THEN
DD1=N_TEMP(I)
*IF,ABS(DD1-TEMP(I)),GT,DD_H,THEN
DD_H=ABS(DD1-TEMP(I))
*ENDIF
*ENDIF
N_TEMP(I)=TEMP(I) ! 计算节点温度（总水头）
N_PRE(I)=N_TEMP(I)-NY(I) ! 计算节点压力，总水头-Y坐标
*ENDDO
*IF,COM_NUM,GT,CHUK_ST+1,THEN
*IF,DD_H,LE,DD_HEAT,THEN
*EXIT
*ENDIF
*ENDIF
/PREP7
! 重新给每个单元设定材料
MATNUM=2
*DO,I,1,E_MAX
*DO,KK,1,4
*GET,NCON(KK),ELEM,I,NODE,KK ! 获取单元四个节点编号
*ENDDO
TEMP_Y=(N_TEMP(NCON(1))+N_TEMP(NCON(2))+N_TEMP(NCON(3))+N_TEMP(NCON(4)))/4 !计算单元中心点平均温度
RESS_T=TEMP_Y-CENTRY(I)
*IF,PRESS_T,GT,0,THEN
RESS_T=0
MPCHG,1,I
*ELSEIF,PRESS_T,LT,-10,THEN
RESS_T=-10
MPCHG,2,I
*ELSE
MP,KXX,MATNUM+1,10**TPRE(PRESS_T)
MPCHG,MATNUM+1,I
MATNUM=MATNUM+1
*ENDIF
*ENDDO
! 重新设定出口边界条件
*IF,CONUTT,GT,CHUK_ST,THEN !前CHUK_ST次采用原边界条件
LSEL,S,LINE,,6
NSLL,S,1
DDELE,ALL,TEMP ! 删除原边界条件
II=0
CHUK_MAXY=0
*DO,JJ,1,N_NUM2
*IF,N_TEMP(N_NO2(JJ)),GE,NY(N_NO2(JJ)),THEN
D,N_NO2(JJ),TEMP,NY(N_NO2(JJ)) ! 总水头=Y坐标
*IF,NY(N_NO2(JJ)),GT,CHUK_MAXY,THEN
CHUK_MAXY=NY(N_NO2(JJ))
*ENDIF
*ENDIF
*ENDDO
*IF,CHUK_MAXY2,NE,CHUK_MAXY,THEN ! 判断前后两次计算的浸润线出口位置是否相同 NSEL,R,LOC,Y,CHUK_MAXY ! 选择最高节点
*IF,CHUK_MAXY,GT,0,THEN
DDELE,ALL,TEMP ! 删除出口最高节点边界条件
*ENDIF
CHUK_MAXY2=CHUK_MAXY
*ENDIF
*ENDIF
ALLSEL,ALL
FINI
/SOLU
SOLVE
FINISH
*ENDDO
SAVE
!迭代计算完毕，进入后处理
FINISH
/POST1
/CLABEL,,1
/EDGE,,0
/CONTOUR,,8,0,1,8
PLNSOL,TEMP ! 显示总水头云图
PLVECT,TF, , , ,VECT,ELEM,ON,0
PLVECT,TF, , , ,VECT,NODE,ON,0
LSEL,S,LINE,,6
NSLL,S,1
PRRSOL,HEAT ! PRINT FLOWRATE THROUGH SOIL FSUM,HEAT ! 计算渗流量
*GET,Q_DAY,FSUM,0,ITEM,HEAT
ALLSEL,ALL
SAVE
*DO,I,1,N_MAX
N_TEMP(I)=TEMP(I) ! 计算节点总水头（温度）
N_PRE(I)=N_TEMP(I)-NY(I) ! 计算节点压力，总水头-Y坐标
DNSOL,I,TEMP,,N_PRE(I) ! 将压力水头值复制到节点
*ENDDO
PLNSOL,TEMP ! 显示压力水头云图
FINI。