ANSYS热分析模块分析渗流场问题的探究

1、ANS Y訣分析模块分析渗流场问题的探究摘要：本文分析了温度场与渗流场的基本理论、 微分方程、初始条件 和边界条件的相似性，为后续利用ANSYS勺热分析模块来求解渗流问题提供了依 据。通过计算黑河金盆大坝0+225主断面的在正常蓄水位下的流场分布情况，验 证了热分析模块对渗流分析结果的可行性与正确性。关键词:ANSYS热分析 渗流场 温度场GEO-SEEPWANSYS是一种应用十分广泛的通用的成熟的完备的有限元工程分析软 件。它具有多种多样的分析能力，从简单的线性静态分析到复杂的非线性动态分 析,除此之外还有方便用户设计的参数化设计语言和产品的优化设计等附加的功 能。ANSYS软件能够提供的分

2、析类型有：结构静力分析；结构动力分 析；结构非线性分析；结构屈曲分析；电磁场分析；声场分析； 压电分析；流体动态分析；热力学分析。其中热分析用于计算一个系统或 部件的温度分布及其它热物理参数，如热量的获取或损失、热梯度、热流密度（热 能量）等。ANSY热分析包括热传导、热对流及热辐射三种热传递方式。此外 ,还 可以分析相变、有内热源、接触热阻等问题。但是没有直接可以利用的渗流分析 模块,所以本文提出利用热分析热传导模块来分析渗流问题。虽然已有比较成熟的渗流分析软件 GEO-SEE P但是仅在正分析方面可 以应用,缺少用户自定义参数化语言，不能二次编程，这为渗流参数设计和反分析 问题带来很大限制

3、和不便。ANSYS勺参数化设计语言（APDL）为设计渗流参数和反 分析渗流参数提供了有力的工具。1渗流场与温度场的理论相似下面从基本理论、微分方程、初始边界条件三个方面来证明温度场与渗 流场的相似性。理论基础的相似根据渗流基本理论可知，对于多孔介质满足达西定律：或者其中：一渗流量；一断面面积；一测压管水头；一渗透系数； 渗径长度；一断面平均流速；一渗透坡降。而对于热传导定律（傅里叶假设）为：或者其中：一热（流）量；一断面面积；一温度场梯度值；一热传导热流强度； 传热系数。微分方程相似渗流场微分方程：对于不可压缩各项异性非均质无源稳定渗流微分方程为 ：对于可压缩各项异性非均质非稳定瞬态渗流微分方

4、程为其中：、一、方向的渗透系数；一单位贮存量。温度场微分方程：对于无热源的各项异性非均质稳定热传导微分方程为 ：对于无热源的各项异性非均质瞬态热传导微分方程为其中：、一、方向的热传导传热系数；一比热。初始条件与边界条件的相似(8)渗流场的初始条件：热传导温度场的初始条件：第一类边界条件：渗流场（9）温度场（10）其中：、一时刻点的测压管水头值和温度值；、一边界上给定的已 知测压管水头和温度函数；一边界上的点。第二类边界条件：渗流场（11）温度场（12）其中：、一沿边界法线方向的渗透系数和导热系数；、一渗流场和温 度场沿边界法线方向的梯度值；、一边界上给定的已知流速和热流强度函数 边界上的点。在

5、不透水边界和绝热边界上，则有：和。位函数特征参数源强渗流场 测压管水头 渗透系数 单位贮存量SS 外渗强度温度场 温度 传热系数 比热 热生产率 表1 ANSYS中渗流场与温度场各种相应量比较 2实例验证为了验证热分析模块对渗流问题分析的可行性和正确性，本文利用西安524m,水库总库容亿m3,有效库容亿m3,大坝市金盆水利枢纽工程(大(II)型二等工程，粘土心墙土石坝)为例。该工程位于西 安市的周至县境内，坝址距周至县城约13km工程以城市供水为主，兼有农田灌 溉结合防洪、发电等综合利用。枢纽工程由拦河坝、泄洪洞、引水洞、溢洪洞和 坝后电站及古河道防渗组成。大坝为粘土心墙砂砾石坝，最大坝高，坝

6、顶高程。坝 顶长度440m,宽11m,最大坝底宽度 属1级水工建筑物。10-3cm/s,反滤层的渗透系数为X 10-5cm/s, 正常蓄水位594m,下游为490m砂卵石坝壳的渗透系数为X 粘土心墙的渗透系数为X 10-7cm/s,大坝0+225横断面(如图1)。根据坝体材料的不同，将坝体分为五个部分， 分别为：上游坝壳、心墙上游反滤层、心墙、心墙下游反滤层和下游坝壳，可以建 立ANSYS几何模型(如图2)。由于边界条件比较复杂，本文选用三节点三角形单PLANE35节点总数为3456,单元总数为1643。图10+225横断面图ANSYS几何模型图对于浸润线的求法，本文利用ANSY防的生死单元技

7、术来确定。(1)首先假设浸润线上下的土体渗透系数一样，利用/solu求出各点的总水头。(2) 利用get命令求得各单元节点上的测压管水头，并与各单元节点的 坐标值(高程)比较。如果测压管水头小于 坐标值,说明该点于浸润线以上，杀死 该单元。找出所有浸润线以上单元，将其杀死。(3) 重新加载边界条件并计算。(4) 重复2、3步,直至各点的水头变化很小。迭代计算中，有可能需要激 活己杀死的单元。由于这种“死活”单元技术的基本对象是单元，因此,若希望足够精度和 足够光滑的浸润线，格划分必须足够密，但是相应的计算量就也增大，或者采用 ANSYS勺格自适应功能在浸润线附近以及水头变化较大区域进行局部格加

8、密。图3等势图  4流速矢量图利用ANSYS的后处理结果查看模块，可以实现对计算结果的图形化和具 体数值化查看。图3和图4分别为三角形土坝计算结果的等势图和流速矢量图。 对于数值化查看就要使用*GET等读取ANSY敎据库的命令进行操作。3结论本文通过对比温度场与渗流场的基本理论、微分方程、初始条件和边界 条件的相似性，为后续利用ANSYS的热分析模块来求解渗流场问题提供了依据 为了验证热分析模块对渗流分析结果的可行性与正确性，本文通过金盆水利枢纽 工程大坝来实例计算，并通过ANSY生死单元技术确定浸润线，通过ANSYS!处理 模块/post1来查看渗流场的等势线和流速分布情况。本文的重要意义在于，利用ANSYS勺APDL参数化及其优化设计，突破了渗 流分析软件GEO-SEEP只能分析渗流正问题，以及对用户不能使用用户参数化定 义的限制，为分析研究渗流问题提供了一种有效的渗流参数设计以及反分析的思 路。参考文献强锋科技,李黎明.ANSYS有限元分析实用教程M.北京：清华大学出版 社,20XX.博弈创作室.APDL