堤坝渗漏及结构破坏条件下渗流传热特性的数值分析.pdf
《堤坝渗漏及结构破坏条件下渗流传热特性的数值分析.pdf》由会员分享,可在线阅读,更多相关《堤坝渗漏及结构破坏条件下渗流传热特性的数值分析.pdf(8页珍藏版)》请在文库网上搜索。
1、水利水电科技进展,2023,43(6)摇 Tel:02583786335摇 E鄄mail:jz 摇 http:/第 43 卷第 6 期Vol.43 No.6水 利 水 电 科 技 进 展Advances in Science and Technology of Water Resources2023 年 11 月Nov.2023基金项目:湖南省水利科技项目(XSKJ2021000鄄9,XSKJ2021000鄄1)作者简介:蒋中明(1969),男,教授,博士,主要从事水工岩土工程安全评价研究。E鄄mail:DOI:10.3880/j.issn.10067647.2023.06.010堤坝渗漏及结
2、构破坏条件下渗流传热特性的数值分析蒋中明1,2,杨江寅1,黄湘宜1,廖峻慧1,肖喆臻1(1.长沙理工大学水利与环境工程学院,湖南 长沙摇 410114;2.长沙理工大学洞庭湖水环境治理与生态修复湖南省重点实验室,湖南 长沙摇 410114)摘要:以洞庭湖区典型土质堤防为对象,采用 FLAC3D软件开展了渗流 温度 应力耦合计算,分析了堤身及堤基在渗漏及结构破坏条件下的渗流场、温度场和损伤场的时空演变特性,探明了温度与渗流、温度与结构破坏之间的内在联系。结果表明,堤坝发生局部渗透破坏乃至整体结构破坏过程中,下游堤脚及地表温度的演变呈阶段性变化特征;地表温度的变化特征可作为预测或判断堤坝内部渗透破
3、坏及结构损伤程度的间接指标。关键词:渗流传热特性;堤坝渗漏;渗透破坏;结构破坏;FLAC3D中图分类号:TV222.2摇 摇 摇 文献标志码:A摇 摇 摇 文章编号:10067647(2023)06006608Numerical analysis of seepage and heat transfer characteristics under dam leakage and structure failureconditions/JIANG Zhongming1,2,YANG Jiangyin1,HUANG Xiangyi1,LIAO Junhui1,XIAO Zhezhen1(1.Sch
4、ool ofHydraulic and Environmenal Engineering,Changsha University of Science&Technology,Changsha 410114,China;2.KeyLaboratory of Dongting Lake Aquatic Eco鄄Environmental Control and Restoration of Hunan Province,Changsha University ofScience&Technology,Changsha 410114,China)Abstract:Taking typical soi
5、l dikes in the Dongting Lake area as the object,the coupling calculation of seepage,temperature and stress was carried out by using FLAC3Dsoftware.The temporal鄄spatial evolution characteristics of theseepage field,temperature field and damage field of dikes and foundations were analyzed under the co
6、nditions of seepageand structural failure.The internal relationships between temperature and seepage and between temperature and structuralfailure were investigated.The results show that the evolution of dike foot and surface temperature in the lower reachesreveals the characteristics of stage chang
7、e during the process of local penetration failure and even the whole structurefailure.The change characteristics of surface temperature can be used as an indirect index to predict or judge the degree ofseepage damage and structural damage in dams.Key words:seepage heat transfer characteristics;dike
8、leakage;seepage failure;structural failure;FLAC3D摇 摇渗透破坏既是威胁堤坝安全的主要问题,也是堤坝破坏的主要形式1。准确探测渗透破坏的产生及扩展,可为防汛抢险的决策提供强有力的科学依据。近年来,以温度为对象的堤坝渗漏隐患探测技术成为研究热点2鄄4。美国、德国等国家最早将温度作为天然示踪剂,通过捕捉温度变化来研究堤坝渗漏情况5。该方法于 20 世纪 80 年代引入国内,后被称为温度示踪渗漏监测技术6。早期温度示踪渗漏监测技术需在水工建筑物或其基础内埋设大量热敏温度计,由于埋设点有限,其对温度场的监测不够准确7。随后兴起的分布式光纤温度测量技术克服
9、了前者的不足,但需在堤坝内埋设大量光纤,工作量依然巨大,且在已建工程中难以实现。目前,红外热成像法凭借其非接触、无污染、速度快、可连续扫描等优点,成为堤坝渗漏探测技术研究的重点发展方向8。深入分析渗漏条件下堤坝内部温度场的变化规律,探明渗漏出口处温度的演变特性是红外热成像堤坝渗漏探测技术的基础。国内外研究人员针对渗漏条件下堤坝温度场变化规律等方面开展了大量的研究,并取得了丰硕的研究成果。理论分析方面,罗日洪等9建立了稳定渗流影响下堤坝渗漏热流耦合模型,分析了坝体中渗流场对温度场的影响;苏怀智等8认为当渗流存在时,水体迁移将迫使土体温度与水温相适应,从而引起温度场的局部不规则变化。室内试验方面,
10、倪枫10通过物理模型试验研究了不66水利水电科技进展,2023,43(6)摇 Tel:02583786335摇 E鄄mail:jz 摇 http:/同入渗水头、渗漏通道及入渗水温等条件下堤坝温度场与渗流场的变化规律;马佳佳等11开展了多种工况下土坝渗漏的红外热像探测试验;Song 等12采用水热耦合分析对堤坝渗流状况进行监测,提出温度可成为判断堤坝前期渗流状况的有效工具。数值模拟方面,陈建生等13提出了利用虚拟热源法研究坝基裂隙岩体中存在的集中渗漏通道,通过温度场来判定钻孔周围是否存在渗漏,并定量计算出渗漏量等参数;Yousefi 等14以 Shamil 大坝为对象,研究了坝体内渗流对压力及温
11、度场的影响,表明温度比压力更能反映渗流状态;张启义等15通过对双层堤基异常渗漏险情下的温度场进行计算分析,在渗漏通道出口处观测到明显的低温区;Cheng 等16基于渗流传热理论,采用数值模拟法对土石坝渗漏过程中的传热特性进行分析,提出孔隙内渗流稳定后,饱和区与非饱和区之间会形成温差,进而可通过温度判断坝体内浸润线位置。堤坝渗漏是一个十分复杂的非稳态过程,尽管堤坝工程渗漏问题已从单一场的研究向多场耦合方向发展,但在数值模拟方面依旧存在考虑的环境条件和边界荷载不够全面等问题。为此,本文在现有研究基础上,采用数值模拟方法,考虑堤坝工程实际边界情况,研究堤坝渗漏乃至结构破坏过程中温度场的演变特性,以期
12、探明堤坝渗漏及结构破坏与温度场演变的内在联系,为基于红外热成像的堤坝渗漏隐患探测技术提供一条有益的实现途径。1摇 数值模拟理论1.1摇 多孔介质渗流传热理论多孔介质渗流传热过程中涉及复杂的耦合作用。渗流方面,堤内及地下水的流动符合达西定律。假设流体和固体不可压缩,流体运动方程可参考文献17。传热方面,堤坝渗流是一个十分缓慢的过程,可忽略土体与渗水之间的能量交换,基于局部热平衡理论的对流 扩散热传导能量平衡方程可参考文献18。1.2摇 多孔介质热流固耦合理论多孔介质中土体、水及温度之间相互影响。基于多孔介质线弹性假设,当产生非等温渗流时,土体本构方程表述为19滓ij=2G着ij+姿着kk+琢P啄
13、ij-3茁TbKb(T-T0)啄ij(1)式中:滓 为应力;着 为应变;i、j、k 代表坐标轴方向;G为剪切模量;姿 为拉梅常数;琢 为比奥系数;P 为孔隙水压力;啄 为 Kronecher 符号;茁Tb为多孔介质热膨胀系数;Kb为多孔介质体积模量;T 为温度;T0为参考温度。1.3摇 土体结构损伤分析理论沈珠江20基于连续介质力学基本概念,提出了考虑黏土结构破坏损伤过程的损伤模型,后被广泛应用于工程界。该土体损伤模型通过定义损伤系数,对各项力学及流体计算参数进行折减,以反映土体的损伤程度。损伤系数和损伤过程中任意时刻土体参数的损伤值为21棕=1-e-(a着v+b着s)(2)E=棕E0+(1-
14、棕)Ed(3)K=棕K0+(1-棕)Kd(4)式中:棕 为损伤系数;a、b 为拟合参数;着v为单元体积应变;着s为单元剪切应变;E、K 分别为损伤后的弹性模量和渗透系数;E0为单元初始时刻弹性模量;Ed为单元定义损伤弹性模量;K0为单元初始时刻渗透系数;Kd为单元定义损伤渗透系数。利用内嵌 FISH 语言实现基于 FLAC3D的土体损伤计算,具体步骤如下:a.读取单元初始参数,并保存在单元额外变量中,作为后续损伤计算的基准。b.读取各单元的体积应变及剪切应变,采用式(2)计算各个单元的损伤系数值,并保存于第二个单元额外变量中。c.根据各个单元的损伤系数值,采用式(3)(4)分别计算损伤过程中各
15、个单元的参数值,并保存于第三个单元额外变量中。d.运用 FISH 的赋值功能,将第三个单元额外变量赋值给模型的对应参数,由此实现对单元参数的损伤。e.通过 FISH 的循环调用功能,将上述步骤反复执行,以模拟计算过程中土体单元实时损伤的效果。2摇 模型有效性验证2.1摇 饱和多孔介质流体热对流计算饱和多孔介质中,流体流动时会携带热量而发生强制对流传热,流体静止时会由于内部温度差异导致密度不均匀,进而发生自由对流传热。Zhao等22提出了一种求解稳态饱和多孔介质流体自由对流问题的渐进逼近方法,解决了采用传统有限元方法难以得到准确对流解的弊端。本文通过引入文献22中饱和多孔介质中流体的热对流计算,
16、验证FLAC3D软件求解多孔介质中渗流传热问题的有效性及准确性。文献22计算模型及监测点布置如图 1所示,模型高度为10km,宽度为10km,厚度为1m,共计网格节点169 个,单元总数144 个,监测点5 个。76水利水电科技进展,2023,43(6)摇 Tel:02583786335摇 E鄄mail:jz 摇 http:/图 1摇 热对流计算模型(单位:km)渗流计算边界:左右两侧及上下两侧均为不透水边界。传热计算边界:顶面为固定温度边界,大小为 20益;底面也是固定温度边界,其值为 220益;左右两侧均为绝热边界。具体计算参数如下:渗透率 为 1伊10-8m/s,孔 隙 率 为 0.1,
17、导 热 率 为0.335 W/(mK),比热容为803.2J/(kgK),热膨胀系数为 7伊10-5K-1。图 2 为初始温度分布和初始孔隙压力分布图。图 2摇 初始温度和初始孔隙压力分布给模型施加顺时针旋转的流速,驱动流体发生强制对流。随着强制对流现象的产生,多孔介质内的温度发生改变。当瑞利数大于临界值时,自由对流现象产生。最后,将初始设置的倾斜加速度去除,由流体内部温度差驱动流体发生自由对流。饱和多孔介质中流体自由对流的温度分布稳态解如图 3 所示(图 3(b)为采用文献22的渐进逼近方法计算得到的解析解)。图 3摇 稳态温度分布(单位:益)通过对比图 2(a)及图 3(a)可知,在自由对
18、流过程中,流体携带热量运动,极大地改变了整体温度分布;对比图3(a)及图3(b)可知,基于 FLAC3D软件计算得到的温度数值解与采用文献22方法计算得到的解析解几乎吻合。为进一步判断两者拟合度,提取监测点温度分布并进行比对,结果见表 1。由表 1 可知,数值解与解析解最大温度差值出现在1 号监测点,为3.4益,误差仅为3.89%。因此,基于FLAC3D软件计算得到的数值解与解析解高度吻合,渗流传热计算的有效性及准确性得以验证。表 1摇 监测点稳态温度值监测点数值解/益解析解/益187.484258.8563119.91204181.11845152.71562.2摇 饱和土热固结计算饱和土在
19、热固结过程中受外部荷载、孔隙压力、温度荷载的同时作用,伴随着孔隙水的排出、温度的扩散,土体内部的有效应力变化较为复杂。白冰23基于饱和多孔介质热 水 力完全耦合的控制方程,通过有限傅立叶变换及其逆变换,推导了土柱内部温度、孔隙压力和位移演化过程的解析表达式,并对一维热弹性固结模型进行解析求解;孙致学等24利用该解析解进行了热流固耦合模型的数值验证。为验证基于 FLAC3D软件求解热流固耦合问题的有效性及准确性,本文引入文献24中经典一维饱和土的热弹性固结算例进行分析。计算模型及监测点布置情况如图 4 所示,模型高度为 7 m,宽度为 3 m,厚度为 1 m,共计网格节点 1 240 个,单元总
20、数570 个,监测点 4 个,分 别 设 置 在 z 坐 标 为 0、3.5、5.0、7.0 m 处。图 4摇 热弹性固结计算模型模型初始孔隙压力为 10 kPa,初始温度为 10益。力学计算边界:模型顶面受向下的均布荷载作用,大小为 10 kPa;模型左右两侧为法向位移约束;底面为固定约束。渗流计算边界:左右两侧及底面均为不透水边界;顶面为自由渗漏边界。传热计算边界:顶面受固定温度荷载作用,大小为 60益;左右两侧及底面均为绝热边界。具体计算参数如下:弹性模量86水利水电科技进展,2023,43(6)摇 Tel:02583786335摇 E鄄mail:jz 摇 http:/为60 MPa,泊
21、松比为 0.4,渗透率为 4伊10-6m/s,孔隙率为 0.2,导热率为 0.836 W/(mK),比热容为167.2 J/(kgK),热膨胀系数为 3伊10-7K-1,密度为2000 kg/m3,流体密度为 1000 kg/m3。图 5 为土体不同位置处竖向位移、孔隙压力及温度的解析解和数值解对比图,可见基于 FLAC3D软件计算得到的数值解与文献24中的解析解高度吻合,其有效性及准确性得以验证。图 5摇 各状态变量解析解和数值解对比表 2摇 计算参数土样类型弹性模量/MPa泊松比密度/(kg/m3)黏聚力/kPa内摩擦角/(毅)渗透率/(10-4cm/s)孔隙率热膨胀系数/(10-6K-1
22、)导热系数/(mK)比热容/(kgK)素填土4.5090.20190030.517.09.810.44812.281.39962粉质黏土3.9450.30186030.215.03.110.48513.411.851650淤泥粉质黏土4.5460.28186017.518.03.010.53811.751.911010粉细砂3.4250.25190009.53.690.41111.381.451450弱透水层38.4600.22200060.530.00.100.21113.631.33656摇 摇 注:表中各参数来自室内试验及护城垸工程地质勘探报告。3摇 模型应用与结果分析3.1摇 工程背景
23、根据湖南省水利水电勘测设计研究总院 2021年 4 月的护城垸工程地质勘察报告,护城垸位于洞庭湖区华容县,西南临藕池河东支,东靠华容河,垸内总面积 364.6 km2,总人口 37.8 万人,是洞庭湖区重点堤垸。该垸一线大堤分华容河段和藕池河段,总长 82.816 km,其堤防级别为域级,建筑物级别为域级。其中,华容河堤段桩号为 70+030 76+790的断面堤身土质量较差,密实程度不均。1998 年、2016 年汛期该堤段堤身、堤脚出现散浸及渗水险情。3.2摇 计算模型以华容河堤段桩号为 71+590 的断面为研究对象,建立堤防计算网格模型。堤顶高程为 35.97 m,宽6.12m,堤身高
24、6.47m,内坡比为1 4,外坡比为1 颐 2.58。堤身为素填土,堤基由粉质黏土、淤泥粉质黏土及粉细砂层组成。计算模型堤身及堤基采用八节点六面体单元划分,共计网格节点数 27882 个,单元总数 13640 个。计算网格及分组情况如图 6 所示,计算参数见表 2。图 6摇 计算网格及模型分组为深入分析下游堤基及地表处渗透破坏及其传热特征,设置 7 个监测点用于监测各状态变量随时间变化的规律。测点1 位于下游堤脚处,测点2、3、4位于下游堤脚附近地表处,测点 5、6、7 位于下游地表以下 0.75 m 处的堤基表层,具体编号及位置如图 7所示。图 7摇 模型监测点位置3.3摇 计算方案及边界条
25、件洞庭湖水位受季节及气候的影响,水位高程变化复杂。为模拟堤防渗流实际情况,拟定计算工况为天然工况、初始工况及高水位工况。其中天然工况及初始工况为稳态计算工况,高水位工况为瞬态计算工况。96水利水电科技进展,2023,43(6)摇 Tel:02583786335摇 E鄄mail:jz 摇 http:/a.天然工况。为模拟堤防枯水期渗流及温度分布,拟定堤防天然工况为上游无水的情况。渗流方面,根据实测资料拟定地下水位。传热方面,模型底部温度为 20益,顶部由于受到气温的作用,拟定温度为 25益。b.初始工况。渗流方面,拟定上游水位为2020 年汛期华容河水文站平均水位,高程为30.8 m。传热方面,
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 堤坝 渗漏 结构 破坏 条件下 渗流 传热 特性 数值 分析