ANSYS软件在平原水库土工膜气胀计算中的应用测试

发布时间：2014-07-27 11:04

　　气胀问题是近年来造成平原水库土工膜防渗体系破坏和失效的主要原因之一[1]，一旦发生将使水库产生严重的集中渗漏，因此合理地布置排气设施是土工膜防渗体系设计的核心问题。以往排气设施的布置主要依赖于经验，缺乏理论指导，且对于复杂边界和膜下多种介质等问题往往无能为力。
　　国内外对于土工膜膜下气体传递和聚集的研究多集中在垃圾填埋场的防渗体系[2]及其排气状态的稳态分析[3]，对于平原水库土工膜的气胀探索尚不充分。ANSYS软件发展至今已相当完善，可以进行结构、热学、电磁场以及各种耦合场的稳态和瞬态有限元分析，其强大的前、后处理功能可以方便地表述模型各点状态沿时间变化的过程，并且利用物理场相似的原理已经在土石坝渗流中取得了应用[4]。本文将从非饱和土气场和温度场的可比拟性出发，用ANSYS热分析模块解决膜下气场的瞬态数值模拟。
　　1. 非饱和土气场与温度场的相似
　　非饱和土中气场与热传导问题的温度场的相似性存在于基本理论、微分方程以及边界条件等三方面。
　　1.1 基本理论相似性
　　热传导问题的傅里叶定律[5]为：
　　（1）
　　式中为单位时间内单位面积的热流量；为导热系数；为温度；为方向上的温度梯度。非饱和土中的气体扩散遵从的Fick定律[6]的表达形式为：
　　（2）
　　式中为通过单位面积土的空气流量；为土中空气流动的传导系数；为空气浓度；为方向上的浓度梯度。由式（1）、（2）可知，非饱和土中气体传导和温度场中热传导的驱动势能都是浓度梯度，在数值计算中具有相同的表达形式。
　　1.2 微分方程相似性
　　三维热传导方程的表达形式为[7]：
　　（3）
　　式中为介质比热容；为密度；为导热系数；为热源强度。
　　在介质各向同性和不考虑热源的情况下（3）式可以简化为：
　　（4）
　　非饱和土中非稳态气流的气象偏微分方程如下：
　　（5）
　　式中为气体孔隙压力；为孔隙水压力；为时间；为气象固结系数，与传导系数有关；为常数；为传导系数。在各向同性问题中，笔耕论文，传导系数不沿程变化，则、和为零。只取地下水位之上的土体作为研究对象，则可假定此部分的孔隙水压力为零，亦为零，则（5）式化简为
　　（6）
　　与各向同性的热传导方程（4）一致。
　　由式（4）、（6）可知，非饱和土中气体传导和温度场中热传导可以用相同的微分方程描述，因此在数值计算中属于同一类问题。
　　1.3 边界条件相似性
　　对于温度场
　　①初始条件： ，即初始时刻温度场各点的温度分布；
　　②第一类边界条件 ： 边界上点在时刻的温度；
　　③第二类边界条件： 边界上点在时刻热流强度函数，特别地，对于绝热边界。
　　对于非饱和土气场
　　①初始条件： ，即初始时刻非饱和土气场各点的孔隙气压分布；
　　②第一类边界条件：
　　边界上点在时刻的孔隙气压；
　　③第二类边界条件： ，边界上点在时刻气流函数，特别地，对于不透气边界。
　　由温度场和非饱和土气场的边界条件可知，两种物理场的边界条件在本质上也是相等的。因此在利用温度场分析非饱和土气体问题时，可以将各种参数对应输入进行求解。
　　2. 应用ANSYS计算非饱和气场
　　2.1 前处理
　　ANSYS中热分析模块专为解决热传导问题而设置，在计算非饱和土气场问题时还需做一定的转化。在前处理中，热分析模块要求输入材料的导热系数（Conductivity）、密度（Density）以及比热容（Specific Heat）。比热容的物理意义为单位质量物质升高单位温度所需热量，其表达式为
　　（7）
　　因此在输入材料的比热容时，需要考虑的是单位质量的土体在气体输入时孔隙气压升高的能力，由理想气体状态方程
　　（8）
　　可知，在体积不变的情况下土体的孔隙气压与体积内空气的物质的量成正比，而该方程也可以写成类似于比热容公式的表达形式
　　（9）
　　孔隙气压在计算初始状态下一般为标准大气压101.325kPa。在不考虑耦合场的情况下，为简便起见，初始状态下单位体积内空气的物质的量取亦取101.325，则土体“比气容”可取1，其余各项赋值1。在ANSYS定义材料性质、输入边界和荷载条件时，各项参数的对比和输入如下
　　表1. ANSYS分析非饱和土气场输入的对应参数
　　在建模和划分网格阶段，需利用对称原理将原本庞大的计算模型分割。在处理对称图形时可以将模型分割成1/2或1/4，在前处理中可使模型细节所占比例扩大和方便网格细化，在后处理中也大大节约了计算时间，同时对于一些关键节点的状态追踪也更为直观。
　　2.2 后处理
　　ANSYS软件包含了强大的后处理功能，为研究非饱和土气体运动的规律提供了强有力的技术支持。对于平原水库气胀现象的研究，在排气设施布置、地下水位上升速率和土体性质各不相同的情况下利用ANSYS后处理可以从不同的模型计算实例中找出各方面影响因素的重点，为之后的工程设计提供理论依据。对于非稳态分析，可以通过对某一特定节点或单元的追踪了解一个特定模型可能出现的最不利状况及其发生位置，并对提出的防治措施提出方案比较。
　　3. 计算实例
　　某水库地下水位为土工膜下5.0m处，经过2个月上升至距膜下1.3m处，土的渗气系数为0.2m/h，排气盲沟为矩形网格排列，相距20m，膜上覆土压重30kPa，计算简图如图1所示。现利用ANSYS软件分析此问题。
 

 

本文编号：7067