寒区高铁路桥过渡段冻结特征及热影响区域数值分析
1.
2.
3.
4.
Numerical analysis on freezing characteristics and thermal influence areas of roadbed-bridge transition section of high-speed railway in cold region
1.
2.
3.
4.
编委: 武俊杰
收稿日期: 2018-07-19 修回日期: 2018-09-26
基金资助: |
|
Received: 2018-07-19 Revised: 2018-09-26
作者简介 About authors
刘华(1983-),男,陕西延安人,副教授,2013年在中国科学院寒区旱区环境与工程研究所获博士学位,从事寒区岩土工程研究.E-mail:
关键词:
Keywords:
本文引用格式
刘华, 牛泽林, 牛富俊, 张楠, 鲁洁.
LIU Hua, NIU Zelin, NIU Fujun, ZHANG Nan, LU Jie.
0 引言
随着我国高速铁路网络的逐步扩展, 在2016年7月多部委联合发布的《中长期铁路网规划》中, 提出了规划建设“八纵八横”的高速铁路网络格局, 高速铁路的总长度要达到3.8×104 km。其中部分高速铁路线路穿越了我国不同冻结深度与气候条件的广大寒区, 面临着不同工程性质的土质条件。为保障线路的安全运营, 无砟轨道路基的轨面变形控制要求极为严格[1-2], 因此工程上要面对和亟需解决大量的冻害问题。同时, 为保证列车高速运行的安全舒适性, 线路中大量采用了高架桥的形式跨越地形和协调纵向坡率, 从而出现相对应的路桥过渡段, 比如我国第一条寒区高铁哈大高铁的线路中桥梁比约为73%[3]。考虑到线路走向、 桥梁净空和局地地形等工程条件, 不同气候分区的路桥过渡段呈现出不同的冻结特征及演变规律, 进而表现出不同的冻土力学行为, 因此需要对路桥过渡段的冻结特征分布特点和工程影响进行深入分析。
冻土的力学行为具有显著的热敏性, 以及动力荷载对冻土的扰动等诸多因素, 结合高速铁路路桥过渡段差异变形的高要求(工后沉降差异≤5 mm, 且折角≤0.001 rad)[1-2], 使冻土地区路桥过渡带的差异变形问题更为复杂。张玉芝等[4]在哈大高速铁路德惠试验点选取了一个典型的路桥过渡段, 实测了2个冻融循环期间的地温和路基变形数据, 探讨了由于冻融过程诱发的差异性沉降演变规律。田亚护等[5]根据青藏高原多年冻土区公路路桥过渡段路基在竣工后3年内的现场试验数据, 分析了路基不同位置的地温变化、 路基基底沉降变形和路基不同位置沿横向及纵向的沉降变形规律。牛富俊等[6]通过对青藏铁路西大滩至尺曲谷地164座桥梁路桥过渡段的沉降病害调查及相关因素分析, 发现过渡段路基沉降与桥梁走向的南北端、 路基坡向、 路基高度、 多年冻土类型(含冰量)、 地温、 路基结构以及地质条件等因素密切相关。金会军等[7]通过对比讨论国内外几条重要的冻土地区温热油管在极端环境下的热稳定性, 提出了在此类工程中需面对的差异性沉降工程问题和解决思路。由此可见, 路桥过渡段的差异性沉降是铁路、 公路、 输油管线等线性工程中普遍存在的问题[8-9], 尤其是运营时在动力荷载作用下的变形表现得更为明显[10-11]。
基于此, 专家学者和工程师们对寒区高铁路基[12-16]以及路桥过渡段路基进行了针对性的研究。刘建坤等[17]基于青藏铁路清水河地区的工程特点, 提出了一种新型的路桥过渡段结构形式, 即通过在路堤底层堆砌一定厚度的片石填料形成通风型路基, 片石区以上采用倒梯形的过渡段填筑形式, 并于路基两侧加设保温护道, 通过实测发现该种结构的过渡段型式呈现出良好的工程效果。考虑到冻土路基变形对高速列车的动力响应, 盛岱超等[18]、 Sheng等[19]基于列车的动荷载、 冻土路基和孔隙水压力的交互影响, 提出了一种可能的冻胀机理, 分析了差异性变形存在的可能原因。在路桥过渡段工程问题的研究和设计中, 由冻融环境引发的差异性变形及其热效应至关重要。而探寻路桥过渡段路基热稳定性的关键性指标之一是地温与变形的差异, 寻找纵向结构的差异性响应引发的影响区域对设计、 运营具有至关重要的意义。
哈大高铁桥梁数目众多, 全线共有中型以上桥梁162座, 路桥过渡段路基的稳定性严重制约着整条线路的安全运营。一旦出现差异性变形过大, 轨枕系统增大的附加应力极有可能出现桥头跳车(落车)、 扣件脱扣、 安全储备降低和拉裂轨道板的事故。影响冻土路基热稳定性的一个关键因素是冻结特征, 尤其在春季融化期, 路桥过渡段三向融化, 使路基下出现冻结夹层, 更进一步扩大了这种影响。哈大线地处我国东北地区, 一般在第二年的3月底4月初开始解冻, 融化期持续长达一个月以上。路桥过渡段的路基填土在差异性地气换热冻结融化过程中, 热影响域的演变过程和结果对填土的热稳定性和力学特征影响甚大, 有必要对其热影响的几何边界和时间尺度进行分析讨论。本文基于有限元和热传导的基本原理, 建立了不同工况下的路桥过渡段数值计算模型, 目的在于评价路桥过渡段路基热影响域范围和探寻较优的台后填筑方式, 分析不同工况下路桥过渡段的冻结特征, 为工程设计和决策提供数据支持。
1 计算模型
1.1 控制方程
线路走向、 倾向和倾角的不同诱发的热不对称性, 使过渡段上的路基稳定性呈现出极大的差异[20-23], 在冻土路基工程中称为“阴阳坡效应”。考虑到阴阳坡效应在数值模型中的实际叠加方法, 本文中仅将不同朝向坡体的温度差异作为阴阳坡效应的体现, 以给定不同的边界温度值作为实际边界条件。且由于路基中粗颗粒填料的含水率和细粒含量较低, 在分析路桥过渡段的路基温度场数值反演过程中, 对过渡段这一线性工程的传热过程作如下假设:
该温度场演变问题是一个带相变的瞬态温度场问题, 加上冻融界面的位置不固定, 因此界面的能量守恒条件是非线性的, 该问题在数学上是一个强非线性问题。本文采用Galerkin法, 将本计算模型热传导控制方程和有限元计算公式导出为
1.2 模型与计算参数
以哈大高铁和哈齐高铁路基为研究对象, 路桥过渡段计算模型如图1所示。计算区域中的地层从上至下依次为: 混凝土强化表层、 级配碎石砂砾土、 非冻胀性A、 B组填料、 普通A、 B组填料、 褥垫层(3层碎石+2层中粗砂)、 黏质黄土层以及桥头边坡上的浆砌块石护坡。同时, 忽略路基主体中的两层复合土工膜和土工格栅的传热效果; 并假设护坡浆砌块石只发生传导传热, 忽略大孔隙的对流传热; 忽略路面上方的轨道板覆盖影响; 忽略路基边坡上防护对热传导的影响。
图1
在路桥过渡段的端部处, 在填土压实的过程中并无侧向限制或者具有较弱的侧限, 路基土体压实并不能达到理想的压实效果, 同时该处土体又受到三个方向的负温冻结作用。因此, 考虑将锥台护坡取消, 换成倒T型的悬臂式钢筋混凝土挡土墙, 而且也可在挡土墙后铺设保温材料, 以便于机械化施工。
表1 三种模拟工况
Table 1
工况 | 编号 | 保温板厚度/cm | 保温板布设方式 | 台前填土形式 | 台后填土形式 |
---|---|---|---|---|---|
现有工况 | I0C1 | 0 | 无 | 有锥台护坡 | 一段过渡 |
模拟工况1 | I0C0 | 0 | 无 | 无锥台护坡 | 分层 |
模拟工况2 | I1C0 | 10 | 紧贴台壁后方, 全高 | 无锥台护坡 | 分层 |
模拟工况3 | I1C1 | 10 | 紧贴台壁后方, 全高 | 有锥台护坡 | 一段过渡 |
图2
表2 模拟计算时的各土层热物理参数
Table 2
土性/材质 | 位置 | 干密度/(kg·m-3) | 含水率/% | 导热系数/(J·m-1·h-1·℃-1) | 比热/(J·kg-1·℃-1) | ||
---|---|---|---|---|---|---|---|
冻 | 融 | 冻 | 融 | ||||
桥台基座C30 | 基床表层上部 | 2 550 | 0 | 5 580 | 5 580 | 976 | 976 |
10 cm强化表层 | 基床表层上部 | 2 450 | 0 | 5 580 | 5 580 | 976 | 976 |
级配砂砾土 | 基床表层 | 1 800 | 7 | 4 789 | 3 988 | 1 731 | 2 033 |
非冻胀A、 B料 | 基床底层上部 | 2 060 | 6 | 5 040 | 4 140 | 707 | 862 |
普通A、 B料 | 基床底层下部 | 2 080 | 6 | 4 830 | 3 680 | 605 | 751 |
褥垫层 | 基床地基交界 | 1 850 | 0 | 1 246 | 1 246 | 649 | 649 |
黏质黄土 | 地基 | 1 540 | 30 | 5 688 | 4 057 | 1 158 | 1 466 |
保温材料 | 挡土墙后 | 30 | 0 | 126 | 126 | 1 250 | 1 250 |
浆砌块石 | 护坡 | 2 480 | 0 | 8 740 | 8 740 | 1 466 | 1 466 |
1.3 边界条件与初始条件
边界条件采用长春地区的气象统计资料和实测地温资料, 计算方法和思路基于文献[26]。边界条件为:
模型四周,
模型下部,
模型上表面,
式中: Tm为地表温度的年平均值, 实测值为9.32 ℃; ΔT为考虑气候变暖的外界气温升高幅度, 东北地区约为0.04 ℃·a-1; A为地表温度波动变化年振幅的一半, 实测值为19.49 ℃; a为路基填筑时刻的时间初始相位, 随路基填筑时间不同而变化, 本例中为
模拟计算采用通用软件ANSYS对模型进行网格划分。路基主体填筑完成于2010年10月30日, 故将该日的实测温度场数值作为初始温度场。路基温度场初始温度条件为
2 结果与分析
2.1 典型路桥过渡段影响区域
图3
图3
不同工况下路桥过渡段地温场分布(20年后4月15日)
Fig.3
The thermal regime of the four different working conditions in roadbed-bridge transition section on April 15, 20th year after the project completion
2.2 不同台后距离的地温场
图4
图4
I0C0不同的台后距离处竖向地温分布
Fig.4
The vertical temperature distribution in different distance aside abutment at I0C0 on March 30 (a) and April 15 (b), 20th year after the project completion
图5
图5
I1C0不同的台后距离处竖向地温分布
Fig.5
The vertical temperature distribution in different distance aside abutment at I1C0 on March 30 (a) and April 15 (b), 20th year after the project completion
图6
图6
I1C1不同的台后距离处竖向地温分布
Fig.6
The vertical temperature distribution in different distance aside abutment at I1C1 on March 30 (a) and April 15 (b), 20th year after the project completion
2.3 不同台后高度的地温场
图7
图7
I0C0不同的台后距离处地温分布
Fig.7
The transverse temperature profiles of the original surface and 2 m above the original surface at I0C0 on March 30 (a) and April 15 (b), 20th year after the project completion
图8
图8
I1C0不同的台后距离处地温分布
Fig.8
The transverse temperature profiles of the original surface and 2 m above the original surface at I1C0 on March 30 (a) and April 15 (b), 20th year after the project completion
图9
图9
I1C1不同的台后距离处地温分布图
Fig.9
The transverse temperature profiles of the original surface and 2 m above the original surface at I1C1 on March 30 (a) and April 15 (b), 20th year after the project completion
2.4 台后热影响区域定量分析
将数值计算得到的地温场数据, 通过插值处理绘制成等温线分布图, 并提取冻结圈图像的数字化数据, 借用计算机图像统计的方法和程序, 获取了不同工况下冻结圈的位置, 以及面积、 周长、 长宽比和横向热扰动距离等数据。图10为提取数字化信息的流程图。图11为路基施工完成后第20年4月15日冻结圈的面积、 周长和所在区域的长宽比(横向长度与纵向高度之比)。由图可见, 在没有保温板情况下, 无锥台和有锥台时结果差异较大, 无论是面积、 周长还是长宽比均存在下降趋势, 这是由于在春融期锥台增大了热交换的幅值。在有保温板情况下, 无锥台和有锥台时结果同前者恰恰相反, 各项几何指标均呈现增大趋势。对比分析可知, 台后保温板的存在直接影响了热交换过程, 保温板比锥台对融化的影响程度大。保温板的隔热效果非常显著, 但同时叠加锥台的影响后, 冻结圈更加扁平, 使长宽比变大, 进而可知纵向设置的保温板隔热效果要远大于横向路基受到大气热扰动的影响。在有锥台情况下, 不论有无保温板, 冻结圈的面积都差别不大, 但周长和长宽比差异较大, 可见同时选取纵向设置的保温板和锥台护坡的情况下, 这种保温效果并不十分显著。该场地的天然冻结深度为1.37 m, 但桥台属于三向冻结, 横向的热扰动深度要大。尽管台后设置了纵向的保温板, 其横向扰动距离也有3.18 m, 为冻深的2.32倍。而在有锥台情况下, 横向扰动距离可达7 m以上。
图10
图11
图11
不同工况下的冻结圈面积、 周长、 长宽比及横向热扰动距离
Fig.11
Frozen area (a), perimeter (b), ratio of length and width (c), and transverse thermal disturbance distance (d) of the freezing ring under the four different working conditions on April 15, 20th year after the project completion
综合分析结果, 结合现场实际工况并考虑工程经济因素, 有保温板无锥台(I1C0)的结构型式热影响区域最小, 冻结圈的几何指标最优。
3 讨论
(1) 桥 - 堤 - 堑连续过渡与阴阳坡效应的耦合
当考虑线路纵向坡率的平顺需求和出现局地过大高差时, 线路设计中需开挖土体, 以至出现桥 - 堤 - 堑的连续过渡段, 这种情形导致在线路一侧方向出现桥梁 - 阳坡(路堤) - 阴坡(路堑)或是桥梁 - 阴坡(路堤) - 阳坡(路堑)的热环境条件。这种工况下, 过渡段的热影响区域变化更加复杂。结合阴阳坡效应, 极有可能在路桥过渡段中, 线路纵向出现刚度先变大后减小, 或是先减小后增大。
(2) 桥台及台后填土模式与冻融作用的耦合
桥台前锥台护坡和台后的填土模式对路桥过渡段热影响区域的判定, 不仅直接受填筑模式、 桥台型制和动力荷载的综合影响, 同时也受热力环境的间接影响。《高速铁路设计规范》[1]推荐一段过渡的倒梯形填筑型式, 可以较好解决变刚度问题, 但此种情况是否具有最佳的热学结果仍不确定。因此, 在力学和热学的影响域中找到平衡点, 仍需要大量论证和分析。
4 结论
通过构建寒区高铁路桥过渡段的三维几何模型, 应用非线性热传导数值分析的手段, 得到了以下结论:
(1) 现有工况(I0C1)下, 台后的冻融影响区域基本在4 m左右, 考虑到台高为4.31 m, 可将春融期热影响区域的横向长度扩展为1.0H~1.5H, 纵向影响深度为该路基条件下的最大冻结深度。
(2) 提出了冻结圈面积、 周长、 长宽比和横向热扰动距离等几何指标, 对比分析了不同工况下的热影响区域和冻结圈的几何特征。
(3) 对比现有路桥过渡段结构和三种优化后结构型式, 发现应用保温板、 倒T型挡土墙的结构形式, 在融期呈现出较小的冻结圈和较快的融化速率, 这让路基纵向刚度变化更加平顺。
参考文献
General description of key technologies for construction of high speed railway in cold area
[J]. ,
严寒地区高速铁路关键施工技术综述
[J]. ,
Roadbed deformation characteristics of roadbed-bridge transition section for high-speed railway in seasonal frozen soil region
[J]. ,
季节性冻土地区高速铁路路桥过渡段路基变形特征分析
[J]. ,
Experiment on ground temperature field and deformation of bridge-subgrade transition section in permafrost regions
[J]. ,
多年冻土区路桥过渡段变形及地温场试验
[J]. ,
Study of the influencing factors of roadbed settlement in embankment-bridge transition section along Qinghai-Tibet Railway
[J]. ,
青藏铁路路桥过渡段沉降变形影响因素分析
[J]. ,
(Differential) frost heave and thaw settlement in the engineering design and construction of oil pipelines in permafrost regions: a review
[J]. ,
多年冻土区输油管道工程中的差异性融沉和冻胀问题
[J]. ,
Study of thermal stability in embankment-bridge transition section on Qinghai-Tibet Railway under climate warming
[J]. ,
气候变暖条件下青藏铁路路桥过渡段长期热稳定性研究
[J]. ,
Analysis and management methods for the road and bridge transition’s sedimentation on the soft clay ground in Shanghai
[J]. ,
上海软土地基路桥过渡段沉降分析及治理方法研究
[J]. ,
Research on the mechanism of differential settlement in embankment-bridge transition section of permafrost regions and prevention methods of the freezing damage
[D].
冻土地区路桥过渡段不均匀沉降的产生机理及冻害防治技术研究
[D].
Dynamic performance analyses and experiment study on bridge/approach embankment of high-speed railway
[D].
高速铁路路桥过渡段动力学特性分析及工程试验研究
[D].
Characteristics of subgrade’s temperature, moisture and frost heave deformation in Haomen, passenger railway line from Lanzhou to Xinjiang
[J]. ,
兰新客运专线浩门区间路基温度、 水分及冻胀变形特征
[J]. ,
The frost heaving deformation of high-speed railway subgrades in cold regions: monitoring and analyzing
[J]. ,
严寒地区高速铁路路基冻胀变形监测分析
[J]. ,
Numerical model of frost heaving and anti-frost heave measures of high-speed railway subgrade in cold regions
[J]. ,
寒区高速铁路路基冻胀数值模型及防冻胀措施
[J]. ,
Characteristics and engineering significance of frost heaving in subgrade of Harbin-Dalian high-speed railway
[J]. ,
哈大高铁路基面冻胀变形特征及工程意义
[J]. ,
Study of water thermal characteristics and frost heave process monitoring in the roadbed of cutting section in seasonally frozen soil regions
[J]. ,
我国季节冻土区路堑段路基水热特性及冻胀发展过程监测研究
[J]. ,
Observation and analysis of a new type embankment-bridge transition section in permafrost regions
[J]. ,
多年冻土区路桥过渡段一种新结构的试验观测与分析
[J]. ,
Effects of train loads on frost heave of embankments
[J]. ,
高速列车与路基冻胀相互作用机理
[J]. ,
A potential new frost heave mechanism in high-speed railway embankments
[J]. ,
Calculation of difference in temperature between sunny slope and shady slope along railways in permafrost regions in Qinghai-Tibet Plateau
[J]. ,
青藏高原多年冻土区铁路路基阴阳坡表面温差的计算
[J]. ,
Temperature and deformation differences between southern and northern slopes of highway embankment on permafrost
[J]. ,
多年冻土区公路路基阴阳坡温度及变形差异分析
[J]. ,
Sunny-shady slope effect on the thermal and deformation stability of the highway embankment in warm permafrost regions
[J]. ,
Analysis of difference in thermal state between south faced slope and north faced slope of railway embankment in permafrost region
[J]. ,
多年冻土区铁路路基阴阳坡面热状况差异分析
[J]. ,
Numerical solution of phase-change problems
[J]. ,
Three-dimensional nonlinear analysis for temperature characteristic of ventilated embankment in permafrost regions
[J]. ,
Numerical analysis for cooling effect of open boundary ripped-rock embankment on Qinghai-Tibetan Railway
[J]. ,
/
〈 | 〉 |