基于二维拉格朗日有限差分法的土石坝渗流和稳定性分析.docxVIP

基于二维拉格朗日有限差分法的土石坝渗流和稳定性分析 0 土石坝渗流流动规律 土石坝广泛应用于中国的节水工程中。根据统计数据，2010年中国共有8399座水库，其中83万多座水库超过83万人。根据研究，30%的洪水是由土石坝和水库底板的不稳定气候造成的。洪水是指水在岩石和裂缝之间流动的现象。其流动性质与岩石和水流的性质有关。现在，假设土石坝的渗透规律符合西方定律。由于渗流骨架中的不同大小和形状，分布也非常复杂，因此很难获得水流的实际流量。在渗透理论上，比较参数和综合参数是难以替代的。当稳定水位时，土石渗透量随着时间的推移而减少。当变化水的作用下，土石渗透规律不符合西方定律。 综合国内外土石坝数值模拟研究现状, 基于二维拉格朗日有限差分法, 建立了二维土石坝有限元分析模型.拉格朗日分析法求解中一般采用3种方法:离散模型法、有限差分法、动态松弛方法.采用ANSYS有限元软件, 通过理论研究和数值模拟相结合的方法, 对土石坝进行渗漏分析、渗透分析、坝坡稳定性分析, 并给出了土石坝加固方案. 1 数值结果分析 达西定律是法国学者达西在1856年提出的, 达西定律只适用于线性阻力关系的测流, 对于惯性力和粘性力支配的渗流, 并不能很好的进行求解, 其认为渗流量: 其中:k—土粒渗透常数;A—断面面积; h—测压管的水头;L—渗流长度. 土粒渗透常数k表征了土粒的渗透强度, 其与土石料材质、土颗粒级配、密实度、水的动力粘度、温度等因素有关. 渗流运动的液体受到质量力和表面力作用, 因此可利用N-S方程、连续性方程、动量方程和能量方程对其进行分析.对于二维空间内的渗流, 如果认为介质和水体是可压的, 则推导出非稳定渗流微分方程为: 其中:S—储水系数;kx、ky、kz—x、y、z方向的质量力. 通过ANSYS有限元软件建立数学模型, 采用具有渗流特性的实体模拟流体流动.本文建立了土石坝流固耦合模型, 流固耦合中孔隙水压力消失导致大坝土体发生沉降, 其中主要包含两个力学变化: (1) 孔隙水压力变化导致大坝有效应力变化, 可能导致塑性区的出现; (2) 孔隙水压力变化导致某一区域内的流体运动状态变化. 渗流数值计算分为各向同性、各向异性两种, 分别采用相应的计算模型分析.非渗流材料采用null进行模拟, 不同区域内可以采用不同的流动模型.流体源采用点源插入到材料中, 对于饱和流动可以采用隐式算法.计算中主要根据达西定律的线性理论, 假设材料参数恒定, 并且不考虑对流.力学过程的特征时间为: 其中:Ku—不排水体积模量;G—剪切模量; ρ—密度;Lc—模型的平均尺寸.流固耦合计算中的储水系数为: 非饱和渗流的流体扩散率为: 其中:M—节点质量;α—扩散系数. 在渗流模型下, 完全流固耦合模型分析需要从以下两方面进行: (1) 孔隙压力改变影响体积应变和应力变化. (2) 体积应变变化会引起孔隙压力变化. 渗流力学扰动都是瞬间出现的, 但是扰动导致的渗流作用时间往往相当长, 例如土体中的孔隙水压力耗散和估计周期往往是几十小时.因此将扩散过程特征时间与流体力学特征时间的比值作为完全流固耦合模型分析的重点. 其中, tcf为扩散过程特征时间. 2 土石坝坝体渗流路径模型 依据渗流理论, 土石坝病害主要是由于渗流和渗漏引起的, 因此有必要建立大坝有限元模.土石坝土体一般采用低渗透粘砂和碎石填筑而成, 中间有一道砖结构的心墙, 心墙的渗透系数远低于土石料, 因此可近似认为心墙的渗透系数为0.将模型坝顶高度设置为1.2 m, 坝基和坝体的参数 (见表1) . 建立大坝的初始状态, 假设初始时大坝土体中不存在超孔隙水压力.利用GAMBIT进行网格划分, 模型为对称结构, 砖砌心墙深入坝基, 有限元模型网格划分 (见图1) . 首先计算初始压力下的孔隙水压和竖向应力情况, 孔隙水压云图中可以清楚的看到心墙上的孔压情况, 有些心墙单元处的孔压不为0, 是因为ANSYS软件对节点信息进行了内插拟合.土石坝初始孔隙压力云图 (见图2) . 逐步上升水位, 计算此情况在水压力的变化, 导出的实验数据显示, 初始孔压为1 000 Pa, 水位上升后的孔压为10 k Pa, 模型内部的竖向应力为50 k Pa, 最大主应力为25 k Pa.水位上升后的孔压云图 (见图3) . 由图3可以看出, 孔压分布自上而下依次增加, 同一深度下孔压云图分布较为均匀, 从坝坡面开始到心墙处终止.利用软件模拟大坝内部的渗流路径, 推断出坝体内浸润线的大体位置.水位上升后的渗流路径分布 (见图4) . 研究表明, 设置的渗流通道不同, 渗流量不同模拟出的渗流路径也不同.水位上升后, 坝体内部孔隙压力、应力都明显增加, 因此水位高度是坝体稳定性的重要影响因素. 3 土石坝防渗设计 通过对土石