考虑太阳热辐射的梯形渠道温...动态数值模拟和衬砌冻胀分析_王建有.pdf
《考虑太阳热辐射的梯形渠道温...动态数值模拟和衬砌冻胀分析_王建有.pdf》由会员分享,可在线阅读,更多相关《考虑太阳热辐射的梯形渠道温...动态数值模拟和衬砌冻胀分析_王建有.pdf(5页珍藏版)》请在咨信网上搜索。
1、第4 1卷第5期2023年5月水 电 能 源 科 学W a t e r R e s o u r c e s a n d P o w e rV o l.4 1 N o.5M a y 2 0 2 3D O I:1 0.2 0 0 4 0/j.c n k i.1 0 0 0-7 7 0 9.2 0 2 3.2 0 2 2 1 3 9 5考虑太阳热辐射的梯形渠道温度场动态数值模拟和衬砌冻胀分析王建有1,石鹏超1,禄林方2,王 娟1,梁太利1(1.郑州大学黄河实验室,河南 郑州 4 5 0 0 0 1;2.赵口引黄灌区二期工程建设管理局,河南 开封 4 7 5 0 0 8)摘要:针对因缺乏渠道表面实测温
2、度而难以准确模拟渠道温度场的问题,考虑太阳热辐射对渠道阴坡、阳坡和底板温度的影响,建立渠道表面增温计算模型,采用有限元软件A B AQU S,对经历三次降温情况下渠道的温度场变化过程进行数值模拟并分析渠道衬砌冻胀情况。结果表明,在太阳热辐射影响下渠道阴坡冻深显著大于阳坡和底板冻深,计算冻深和衬砌破坏规律也与实际情况吻合较好,在降温过程中阴坡的混凝土衬砌相较于阳坡和底板更容易发生冻胀破坏,按照多年最大冻深利用规范计算时也有相同规律,研究成果可为渠道混凝土衬砌温控防裂设计提供新思路。关键词:渠道;混凝土衬砌;太阳热辐射;数值模拟中图分类号:T V 9 1;TU 5 2 8 文献标志码:A 文章编号
3、:1 0 0 0-7 7 0 9(2 0 2 3)0 5-0 1 0 9-0 4收稿日期:2 0 2 2-0 7-0 9,修回日期:2 0 2 2-0 8-0 6基金项目:国家自然科学基金项目(5 2 1 7 9 1 4 5);河南省2 0 2 1年水利科技攻关项目(Z K E Q-Q D C QMN F K)作者简介:王建有(1 9 7 7-),男,博士、教授级高级工程师,研究方向为水工建筑物的安全性检测与评价,E-m a i l:w a n g j i a n y o u z z u.e d u.c n1 引言混凝土衬砌冻胀问题是渠道工程中亟待解决的问题。目前,关于渠道冻胀问题的研究成果较
4、丰富,但现有研究大多将冻胀的动态过程看作稳态传热过程,无法准确反映温度波动下的冻胀问题;另外,在温度边界条件的取值上,由于缺乏长时段渠表温度数据,通常直接利用气温数据预测渠道冻胀破坏,误差较大。因此,合理设定温度场的温度边界条件成为渠道温度场数值模拟难题。对此,根据混凝土水平表面太阳辐射的增温计算方法1,从太阳热辐射角度出发,结合气象学知识,对渠道阴坡、底板和阳坡赋予不同太阳热辐射增温,建立渠道表面增温计算模型。在此基础上,基于有限元软件A B AQU S对赵口引黄灌区的梯形渠道进行8 5 d温度场仿真计算,模拟渠道冻深演化过程,以期掌握外界温度变化过程中渠道冻胀的动态变化规律,以期为渠道抗冻
5、胀设计提供理论依据。2 太阳热辐射增温模型太阳辐射能量中,一部分到达地表前被大气吸收;另一部分到达地表后,又有一部分被吸收,余下部分则被反射,其中被地表吸收的能量引起了地表温度的升高2。设在单位时间内单位面积上,水平混凝土板表面太阳辐射热均值为S,其中被混凝土吸收的部分为R。水平混凝土板表面的太阳辐射增温公式1为:T=R/(1)R=sS(2)S=S0(1-b n)(3)=2 1.8+1 3.5va(4)式中,T为混凝土表面温度与大气温 度的差值,;R为 混 凝 土 表 面 吸 收 的 太 阳 辐 射 能,k J/(m2h);为 混 凝 土 表 面 放 热 系 数,k J/(m2h);s为吸收系
6、数,混凝土表面为0.6 5;S为太阳辐射热均值,k J/(m2h);S0为晴天太阳辐射热,k J/(m2h);b为云度系数;n为云量;va为风速,m/s。式(1)(4)为理想情况下的太阳辐射增温计算,即无遮挡物的理想水平面。对于渠道工程,尤其是东西走向渠道,阴坡和渠底会遮蔽一部分阳光,使得阴坡和渠底的太阳辐射增温比阳坡小。不同渠道断面形状不同,渠道边坡板坡度越大,渠道所在地纬度越高,这一差异越显著。为考虑这部分增温差异,提出太阳辐射增温系数k,定义为任意坡面任意坡比的太阳辐射增温与理想水平平面太阳辐射增温的比值。不同坡面的太阳辐射增温计算公式为:k=k(5)=s i ni(6)k=(m a x
7、-i)/m a x(7)T=kT(8)式中,为太阳光照效率,与太阳光线和受光面夹角的正弦值成正比;k 为角度系数,理论意义为实际可曝晒角度与理论可曝晒角度的比值;i为太阳光线和受光面的夹角;m a x为正午太阳高度角;i为不同坡面开始接收太阳辐射的角度(图1),由于阴坡坡顶的遮挡,当太阳光线与地面夹角大于1时阳面中部开始接收太阳光辐射,当太阳光线与地面夹角大于2时底板中部开始接收太阳光辐射,当太阳光线与地面夹角大于3时阴面整体开始接收太阳光辐射;T 为渠道混凝土衬砌不同坡面的太阳辐射增温值。ABC平行光线平行光线平行光线阴坡阳坡底板北东南西321321321图1 太阳辐射规律F i g.1 S
8、 o l a r r a d i a t i o n p a t t e r n s3 计算模型3.1 基本假定(1)冻土和未冻土均为各向同性材料,在不同方向上热传导性能相同。(2)由于渠道断面在渠长方向基本无变化,因此采用二维模型进行模拟研究。3.2 热力学温度场方程由于计算区域内存在一个随时间变化的冻融两相分界面,在考虑水冰相变释放潜热时,需在冻土材料属性中设置潜热参数,张玉芝等3认为水冰相变区间为0,-1。因此本研究中考虑土骨架、介质传导及水冰相变3个主要因素的影响,渠道基土中的温度分布满足方程:cT/t=(-T)+LVi/t(9)式中,为材料密度,k g/m3;c为材料质量比热容,J/
9、(k g);T为温度,;t为计算时间,s;为哈密顿算子;为材料导热系数,W/(m);L为水冰相变潜热,k J/k g;Vi为体积含冰率。4 有限元数值模拟4.1 原型渠道基本情况以赵口灌区石岗分干桩号1+0 5 0断面渠道为模拟对象,渠道为东西走向,渠道断面尺寸见图2,边坡系数12,衬砌板厚度为1 0 c m。工程位于河南省开封市,在2 0 2 1年1月出现了大范围渠道衬砌 开 裂 现 象,据 气 象 资 料 显 示,该 工 程 在2 0 2 0-1 2-1 42 0 2 0-1 2-1 7、2 0 2 0-1 2-2 92 0 2 0-1 2-3 1、2 0 2 1-0 1-0 62 0 2
10、 1-0 1-0 8先后经历3次强冷空气,开封通许县具体气温变化过程见图3。3004 4004 0004 40030012 800100310300100500300300100设计渠底2 2001 99010.451:21:2图2 石岗分干1+0 5 0断面图(单位:mm)F i g.2 S h i g a n g s u b-s t e m 1+0 5 0 c r o s s-s e c t i o n a l v i e w温度/日最高温度日最低温度日平均温度2020-11 01-2020 11 11-2020-11-212020-12-012020-12-112020-12-21202
11、0-12-312021-01-102021-01-20年-月-日图3 通许县2 0 2 0-1 1-0 12 0 2 1-0 1-2 4气温变化过程 F i g.3 T e m p e r a t u r e c h a n g e p r o c e s s o f T o n g x u C o u n t y f r o m 2 0 2 0-1 1-0 1 t o 2 0 2 1-0 1-2 44.2 有限元模型4.2.1 有限元几何建模及单元网格划分在有限元建模时,取渠底衬砌中点为整体坐标原点,渠道断面向两侧沿伸3 m作为模型左右边界,渠底向下取5 m作为模型下边界。有限元模型采用四边
12、形网格划分单元,各部分单元保证连续。渠道左右边界和下边界在空间上无限延伸,认为其边界热通量为零。渠道上边界采用温度边界条件。阳坡、底板和阴坡考虑太阳辐射增温折算后的日最低温度作为其温度边界条件。梯形渠道有限元网格见图4。4.2.2 材料参数选择与分析步设定分析步选择瞬态热传导进行计算,根据原型渠道气象资料确定瞬态求解总时间为8 5 d(2 0 2 0/1 1/12 0 2 1/1/2 4)。在参数选取方面,参考王正中等4,5关于混凝土和土体材料的取值,见表1。011水 电 能 源 科 学 2 0 2 3年 第4 1卷第5期王建有等:考虑太阳热辐射的梯形渠道温度场动态数值模拟和衬砌冻胀分析图4
13、石岗分干1+0 5 0梯形渠道网格划分F i g.4 T r a p e z o i d a l c h a n n e l g r i d d i v i s i o n o f S h i g a n g s u b-s t e m 1+0 5 0 s e c t i o n表1 材料参数表T a b.1 T a b l e o f m a t e r i a l p a r a m e t e r s材料类别密度/(k gm-3)导热系数/(Wm-1-1)比热容/(Jk g-1-1)热线膨胀系数/-1混凝土2 4 0 01.2 89 7 01.1 1 0-5未冻土1 3 0 01.5 4
14、1 7 8 00冻土1 3 0 02.4 01 7 8 0-1.4 6 1 0-34.2.3 阳坡、底坡与阴坡的温度边界条件设定理想水平混凝土平面的太阳辐射增温值可根据式(1)(4)得出。根据气象资料,计算周期内平均晴天太阳辐射热S0为6 1 1.3 7 k J/(m2h),云度系数b为0.6 8,平均云量n为0.2,平均风速Va为2.3 9 m/s。由此可计算出赵口灌区所在地水平混凝土表面吸收的太阳辐射能R为5 2 8.2 2 k J/(m2h),T为6.3 4,该值即为理想水平混凝土平面下太阳辐射的增温值。考虑到该工程为东西走向渠道,存在与太阳光线夹角不同的阳坡、底板和阴坡,因此根据式(6
15、),阳坡的太阳光照效率阳为0.8 4 3,底板的太阳光照效率底为0.5 1 5,阴坡的太阳光照效率阴为0.1 7 8。赵口灌区冬季正午太阳高度角即m a x为3 1,根 据 式(7),k阴为0.1 4 5,k 底为0.3 8 7,k 阳为0.8 0 8,因此根据式(5)得到太阳辐射增温系数k阴为0.0 1 1,k底为0.2,k阳为0.6 8 1。结合式(9)可得阴坡增温为0.0 7,底板增温为1.2 7,阳坡增温为4.3 2。4.3 结果与分析4.3.1 温度场分析从气温变化图(图3)可知,在2 0 2 0-1 1-0 12 0 2 1-0 1-2 4日渠道所在地一共经历了三次强冷空气。温度场
16、随时间的变化过程见图5,为方便观察测量冻深,温度场云图中将冻土区域处理为白色显示。在第一次强冷空气影响下,最大冻深出现在2 0 2 0-1 2-2 2日,见图5(a),其中左岸为阴坡,右岸为阳坡。此时阴坡冻深为1 4 c m,渠底冻深为7 c m,阳坡未出现冻深。在第二次强冷空气影响下,最大冻深出现在+1.60010+1.533+1.466+1.400+1.333+1.266+1.200+1.133+1.066+9.997+9.331+8.664+7.998+7.331+6.665+5.998+5.332+4.665+3.999+3.332+2.666+1.999+1.333+6.665+0.
17、000-1.7301101010101010101010101010101010101010101010101010101111111100000000000000-100+1.59710+1.530+1.46411010+1.39710+1.33410+1.26410+1.19810+1.13110+1.06510+9.98010+9.31510+8.64910+7.98410+7.31910+6.65310+5.98810+5.32310+4.65710+3.99210+3.32710+2.66110+1.99610+1.33110+6.65310+0.00010-6.7301011111
18、11100000000000000-100(a)第一次强冷空气下温度场云图(b)第二次强冷空气下温度场云图(c)第三次强冷空气下温度场云图(d)三次强冷空气结束时温度场云图+1.59710+1.527+1.46111010+1.39410+1.32810+1.26210+1.19510+1.12910+1.06210+9.96110+9.29710+8.63210+7.96810+7.30410+6.64010+5.97610+5.31210+4.64810+3.98410+3.32010+2.65610+1.99210+1.32810+6.64010+0.00010-9.7901011111
19、11100000000000000-101+1.59110+1.524+1.45811010+1.39210+1.32610+1.25910+1.19310+1.12710+1.06010+9.94210+9.27910+8.61610+7.95310+7.29110+6.62810+5.96510+5.30210+4.63910+3.97710+3.31410+2.65110+1.98810+1.32610+6.62810+0.00010-9.516101111111100000000000000-10-1温度/温度/温度/温度/图5 三次冷空气下温度场云图F i g.5 T e m p e
20、 r a t u r e f i e l d c l o u d s u n d e r t h r e e c o l d s p e l l s2 0 2 1-0 1-1 0日,见图5(b),此时阴坡冻深 为2 3 c m,渠底冻深为1 4 c m,阳坡冻深为7.5 c m。在第三次强冷空气影响下,最大冻深出现在2 0 2 1-0 1-1 9日,见 图5(c),此 时 阴 坡 冻 深 为3 0 c m,渠底冻深为1 6 c m,阳坡未出现冻深。在2 0 2 1-0 1-2 02 0 2 1-0 1-2 4,此时渠道冻深趋于稳定状态,在2 0 2 1-0 1-2 4日,见图5(d),此时阴坡
21、冻深为3 0 c m,渠底冻深为1 7 c m,阳坡未出现冻深。在渠道阴坡出现冻胀裂缝后,地质人员于2 0 2 1-0 1-2 4日在石岗分干测得实地冻深结果,在渠道阴坡有约2 7 c m冻深,阳坡未见冻结层,与温度场模拟结果吻合较好。温度场计算结果表明,在渠道走向影响下,温度场呈不对称分布,阳坡温度明显大于阴坡,温度场的分布规律决定了冻深情况,这与现场检测中阴坡冻深大于阳坡冻深的规律相吻合。主要原因在于该渠道为东西走向渠道,不同坡面太阳热辐射增温的不同导致阳坡接受的太阳热辐射大于底111板与阴坡,结合现场检测结果可知,数值结果与实际结果冻深最大差异为3 c m且分布规律基本吻合。这表明此模型
22、可较好模拟太阳辐射作用下的渠道温度场。4.3.2 数值模拟结果分析此工程所在地地下水位埋深约2 0 0 c m。依据 渠系工程抗冻胀设计规范63.2.2.1中冻深与冻胀量的关系曲线中查得冻深为1 0 c m时冻胀量约为0.2 c m,冻深2 5 c m时冻胀量约为0.5 c m,冻深4 0 c m时冻胀量约为1 c m。规范6规定对于梯形混凝土衬砌渠道,当衬砌板单块长边尺寸5 m时允许冻胀量为0.5 c m。说明当冻深大于2 5 c m时产生的冻胀量大于规范6允许冻胀量。温度场模拟结果显示,渠道阴坡在第一、二、三次强冷空气后最大冻深分别为1 4、2 3、3 0 c m,若不采取防冻胀措施,阴坡
23、渠道衬砌板易发生冻胀破坏。渠道阳坡经历三次冷空气降温后最大冻深不超过7.5 c m,远小于规范6允许冻深值,因此阳坡衬砌板不会冻胀破坏。渠道底板最大冻深出现在第三次强冷空气后,达到了1 6 c m,也未达到规范6允许冻深值,因此渠道阳坡和底板均不会发生冻胀破坏。综上,此次强冷空气作用下,阴坡混凝土衬砌板相较于阳坡和底板更容易发生冻胀破坏。在赵口引黄灌区石岗分干桩号1+0 5 0断面渠道现场勘察中,发现阳坡和底板部位的混凝土衬砌未出现冻胀开裂现象,在阴坡的混凝土衬砌出现了长距离沿渠长方向的冻胀开裂现象。这与数值模拟结果的冻胀破坏规律吻合较好,验证了数值模拟分析方法的合理性。4.3.3 依据3 0
24、年最大冻深参考规范计算依据 民用建筑供暖通风与空气调节设计规范7对中国各大城市最大冻深的汇总可知,河南开封 历 史 记 录1 9 7 12 0 0 0年 最 大 冻 深 为2 6 mm。将其作为历史最大冻深,按照 渠系工程抗冻胀设计规范6计算渠道各部位的冻深和冻胀量,与本次极端强冷空气下的数值模拟结果具体对比见表2。由表2可知,当按照历史最大冻深参照 渠系工程抗冻胀设计规范6推算冻胀破坏结果时,渠道阴坡混凝土衬砌会发生冻胀破坏,渠道阳坡和底板混凝土衬砌均不会发生冻胀破坏,这与本次数值模拟分析结果一致。即在面临极端强冷空气或历史记录多年最大冻深时,阴坡混凝土衬砌相较于底板和阳坡的混凝土衬砌,更容
25、易产生冻胀破坏。表2 石岗分干桩号1+0 5 0冻胀量计算表T a b.2 C a l c u l a t i o n t a b l e f o r f r o s t h e a v e f o r S t o n e G a n g S u b-d r y 1+0 5 0c m 数值模拟结果与规范计算结果对比计算值阳坡阴坡底板规范计算结果渠道各部位计算冻深8.5 1 3 0.7 7 2 3.7 5冻胀量0.1 70.6 90.4 8数值模拟结果冷空气作用下最大冻深7.5 0 3 0.0 0 1 6.0 0冻胀量0.1 50.6 70.3 2是否破坏依据规范允许冻胀量0.5 05 结论a.
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 考虑 太阳 热辐射 梯形 渠道 动态 数值 模拟 衬砌 冻胀 分析 王建有
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。