CN110390176B - 一种无砟轨道冻结与损伤行为计算方法 - Google Patents
一种无砟轨道冻结与损伤行为计算方法 Download PDFInfo
- Publication number
- CN110390176B CN110390176B CN201910698540.4A CN201910698540A CN110390176B CN 110390176 B CN110390176 B CN 110390176B CN 201910698540 A CN201910698540 A CN 201910698540A CN 110390176 B CN110390176 B CN 110390176B
- Authority
- CN
- China
- Prior art keywords
- water
- freezing
- ice
- porous
- temperature
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000007710 freezing Methods 0.000 title claims abstract description 97
- 230000008014 freezing Effects 0.000 title claims abstract description 90
- 238000004364 calculation method Methods 0.000 title claims description 14
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 120
- 239000011148 porous material Substances 0.000 claims abstract description 97
- 238000000034 method Methods 0.000 claims abstract description 34
- 238000009826 distribution Methods 0.000 claims abstract description 32
- 230000006399 behavior Effects 0.000 claims abstract description 28
- 230000008569 process Effects 0.000 claims abstract description 23
- 238000012360 testing method Methods 0.000 claims abstract description 20
- 238000013508 migration Methods 0.000 claims abstract description 16
- 230000005012 migration Effects 0.000 claims abstract description 16
- 230000008878 coupling Effects 0.000 claims abstract description 13
- 238000010168 coupling process Methods 0.000 claims abstract description 13
- 238000005859 coupling reaction Methods 0.000 claims abstract description 13
- 239000013354 porous framework Substances 0.000 claims abstract description 13
- 238000010257 thawing Methods 0.000 claims abstract description 12
- 238000007906 compression Methods 0.000 claims abstract description 11
- 230000006835 compression Effects 0.000 claims abstract description 10
- 230000005540 biological transmission Effects 0.000 claims abstract description 8
- 230000003993 interaction Effects 0.000 claims abstract description 6
- 238000005516 engineering process Methods 0.000 claims abstract description 5
- 238000002591 computed tomography Methods 0.000 claims abstract description 4
- 239000013078 crystal Substances 0.000 claims description 19
- 230000008859 change Effects 0.000 claims description 18
- 238000002425 crystallisation Methods 0.000 claims description 14
- 230000008025 crystallization Effects 0.000 claims description 14
- 230000035699 permeability Effects 0.000 claims description 10
- 230000000694 effects Effects 0.000 claims description 7
- 230000009471 action Effects 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 5
- 239000007787 solid Substances 0.000 claims description 5
- 238000001179 sorption measurement Methods 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 2
- 230000002706 hydrostatic effect Effects 0.000 claims description 2
- 238000011946 reduction process Methods 0.000 claims description 2
- 238000003860 storage Methods 0.000 claims description 2
- 239000010410 layer Substances 0.000 description 40
- 230000007423 decrease Effects 0.000 description 7
- 241001669679 Eleotris Species 0.000 description 4
- 239000004568 cement Substances 0.000 description 4
- 238000001816 cooling Methods 0.000 description 4
- 238000009413 insulation Methods 0.000 description 4
- 238000004321 preservation Methods 0.000 description 4
- 101100532451 Rattus norvegicus Slc22a17 gene Proteins 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000009467 reduction Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000015556 catabolic process Effects 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000011065 in-situ storage Methods 0.000 description 2
- 238000012423 maintenance Methods 0.000 description 2
- 239000012466 permeate Substances 0.000 description 2
- 239000002002 slurry Substances 0.000 description 2
- 238000005481 NMR spectroscopy Methods 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000005266 casting Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000010485 coping Effects 0.000 description 1
- 230000001808 coupling effect Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000001351 cycling effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002194 freeze distillation Methods 0.000 description 1
- 230000036571 hydration Effects 0.000 description 1
- 238000006703 hydration reaction Methods 0.000 description 1
- 239000005457 ice water Substances 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 238000004215 lattice model Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000003204 osmotic effect Effects 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Railway Tracks (AREA)
- Porous Artificial Stone Or Porous Ceramic Products (AREA)
Abstract
本发明公开了一种无砟轨道冻结与损伤行为计算方法,包括S1、根据动弹性波试验获取已浇筑混凝土试件的多孔骨架体积压缩模量,并计算得到Biot系数;S2、根据CT扫描和图像识别技术,获取混凝土试件孔结构分布参数和整体孔隙率;S3、根据所述孔结构分布参数,计算得到混凝土试件冻结过程中温度‑冻结速率的关系曲线;S4、根据Young‑Laplace方程冰‑未冻水平衡以及冰‑吸附水膜平衡关系,构建平均孔隙压力与温度、孔结构之间的函数关系;S5、基于多孔介质力学、多孔介质内水分迁移达西定律和多孔体系热传导Fourier定律,构建多孔介质应力场‑温度场‑渗流场的耦合控制方程,得到冻融过程中无砟轨道内结构受力、热量传输和水分迁移之间的相互作用关系。
Description
技术领域
本发明属于轨道冻结与损伤行为的技术领域,具体涉及一种无砟轨道冻结与损伤行为计算方法。
背景技术
我国地域辽阔,有相当大的部分地区处于严寒地带,而该区域混凝土结构面临强烈的冻融循环作用。气象统计数据结果表明,兰州地区的年平均冻融循环次数达83.6次,黑龙江哈尔滨冻融循环次数约45次作用,即便江苏南京地区的年平均冻融循环次数也达到19.9次,冻融破坏是我国北方寒冷地区混凝土结构损伤失效的主要原因之一。
无砟轨道结构设计时,虽然考虑了轨道结构的排水及抗冻性,并做了相应的排水及抗冻性设计。但由于施工和养护不当,在雨水丰富地区和排水不畅地段导致有大量的水会沿着表面裂缝或层间离缝渗入轨道结构内部,即使排水较好的区域由于低温环境下轨道板积雪,温度升高融化将导致很大比例的水分沿孔隙及微裂隙渗流进入轨道结构内部。因此,北方寒冷地区部分无砟轨道已出现较为严重的冻结与损伤行为。其中,I型双块式无砟轨道主要由钢轨、扣件、预制轨枕、道床板、混凝土底座或支承层组成,在我国铺设广泛,其中大西、西宝、兰新、宝兰等严寒地带都大量采用了I型双块式为主的轨道结构形式。双块式无砟轨道在桥梁地段采用单元结构,路基和隧道地段采用连续结构。其轨枕为工厂预制,道床板、底座板或支承层采用现场浇筑,对桥梁、路基、隧道等线下基础的适应性较好,但由于为现浇结构,浇筑养护后极易产生干缩裂缝,水分更易渗入轨道内部,因此面临更为严峻的冻结与损伤行为威胁。
Powers提出的静水压理论为混凝土冻融受力的量化分析奠定了基础,此后国内外学者展开了大量研究,陆续提出了渗透压理论(Powers)、结晶压理论(Scherer)、微结冰理论(Setzer)、粘结剥落理论(Valenza)以及基于热力学原理建立的联系宏观结构与微观受力的多孔介质力学理论(Coussy)。在数值计算方面,T.Ueda等人(2009)基于刚体弹簧方法(RBSM),通过引入零强度单元及塑形拉伸弹簧,建立了混凝土中尺度计算模型,分析了冻融引起的应力演化、裂纹扩展、强度退化问题以及冻融后混凝土的受弯性能;Lin Liu,WeiSun等人采用水泥浆体数值微观结构,基于离散化的格构单元法,建立了含量化微观损伤的三维格构模型,分析了饱水状态、外荷载及冻结共同作用下水泥浆体微观结构的受力及微裂纹的分布;Qingli Dai在细观角度采用扩展有限元(XFEM)研究了过冷条件下孔内结晶压力对混凝土内裂纹扩展的影响。
上述研究均通过将经验或理论预测的孔隙压力以外荷载的方式引入到模型中,并未考虑材料应力与孔隙压力的耦合关系。基于水-热-力耦合多孔介质力学物理表达式,采用毛细孔内冰水压力平衡关系,B.Zuber,J.Marchand在微观尺度下研究了低温冻结过程中气孔尺寸及间隔尺寸对水泥基材料变形的影响;段安分析了混凝土试块受低温冻结过程中的温度场、应变场以及孔内压力场的分布规律;曾强研究了冻结速率及气孔内压力边界条件对水泥基材料变形的影响。而混凝土冻结与损伤行为研究在无砟轨道领域中应用尚浅,低温冻结与损伤行为对无砟轨道的影响规律尚不明确。
发明内容
本发明的目的在于针对现有技术中的上述不足,提供一种无砟轨道冻结与损伤行为计算方法,以解决或改善上述的问题。
为达到上述目的,本发明采取的技术方案是:
一种无砟轨道冻结与损伤行为计算方法,其包括:
S1、根据动弹性波试验获取已浇筑混凝土试件的多孔骨架体积压缩模量,并计算得到Biot系数;
S2、根据CT扫描和图像识别技术,获取混凝土试件孔结构分布参数和整体孔隙率;
S3、根据所述孔结构分布参数,计算得到混凝土试件冻结过程中温度-冻结速率的关系曲线;
S4、根据Young-Laplace方程冰-未冻水平衡以及冰-吸附水膜平衡关系,构建平均孔隙压力与温度、孔结构之间的函数关系;
S5、基于多孔介质力学、多孔介质内水分迁移达西定律和多孔体系热传导Friour定律,构建多孔介质应力场-温度场-渗流场的耦合控制方程,得到冻融过程中无砟轨道内结构受力、热量传输和水分迁移之间的相互作用关系。
优选地,步骤S1中计算得到的Biot系数b为:
其中,K0为多孔骨架压缩模量,Km为固体基质压缩模量。
优选地,步骤S3中计算得到混凝土试件冻结过程中温度-冻结速率的关系曲线为:
Vw-i为孔隙内冰的体积含量,其计算表达式为:
Vw-i=φ(Rpeq)-Vads(Rpeq)
吸附层的体积含量Vads为:
孔径Rpeq为:
Rpeq=Req+δ
其中,Req为冰晶体尖端平衡半径;δ为孔隙水膜厚度之和,其表达式为:
其中,θ为摄氏温度;
冰晶体尖端平衡半径Req为:
其中,γls为冰与水表面张力,T0为正常水冰点,Rmax为冰晶体尖端最大半径。
优选地,步骤S4中构建平均孔隙压力与温度、孔结构之间的函数关系为:
P*(θ)=Pw(θ)+X
其中,Pi为孔隙冰压力,Pw为孔隙水压力,γ为冰与水表面张力,Req为冰晶体尖端平衡半径,πw为吸附层水压力,P*(θ)为平均孔隙压力,Pw(θ)为结冰孔,X为未结冰孔;
未结冰孔X的表达式为:
其中,n为整体孔隙率,Rpeq(θ)为结晶的最小孔隙半径,Req(θ)为冰晶体尖端平衡半径。
优选地,步骤S5中构建多孔介质应力场-温度场-渗流场的耦合控制方程分别为:
当不考虑体力的作用下,多孔介质微元的应力场控制方程为:
其中,K0为多孔骨架压缩模量,α0为多孔骨架的体膨胀系数;
忽略对流对热传导的影响,并考虑相变潜热,则多孔介质的温度场控制方程为:
其中,Cm、Cw和Cw分别表示多骨架、水以及冰的比热容,λef、λm、λw、λi分别为多孔体系、骨架、水以及冰的热传导系数,Sw为水饱和度,ρw为水密度,Si为冰饱和度,ρi为冰密度;
混凝土冻结过程中的内水分迁移,即渗流场的控制方程为:
其中,β为多孔弹性储水修正系数,Kw为水的压缩模量,Ki为冰的压缩模量;
其中,αw、αi、α0分别为水、冰、多孔骨架的体膨胀系数。
本发明提供的无砟轨道冻结与损伤行为计算方法,具有以下有益效果:
本发明通过多孔介质内的渗流、传热过程、应力场本构及相互耦合作用,考虑孔隙结构对水分冻结温度的影响,并基于弹性多孔介质力学,研究了冻融过程中无砟轨道内结构受力、热量传输和水分迁移之间的相互作用规律,为无砟轨道混凝土冻结与损伤行为问题的应对及防治工作提供理论指导。
附图说明
图1为孔结构分布曲线。
图2为双块式无砟轨道横截面示意图。
图3为Paste I试验与数值结果。
图4为Paste II试验与数值结果。
图5为水压力空间分布。
图6为支承层孔内水压力演化。
图7为道床板孔内水压力演化。
图8为道床板和支承层结晶压力。
图9为水压力源项组成及演化。
图10为支承层结冰速率。
图11为道床板结冰速率。
图12为道床板和支承层冰饱和度。
图13为温度空间分布。
图14为支承层降温速率。
图15为道床板降温速率。
图16为Tresca应力空间分布规律。
图17为支承层Tresca应力演化规律。
图18为道床板Tresca应力演化规律。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
根据本申请的一个实施例,参考图1,本方案的无砟轨道冻结与损伤行为计算方法,包括:
S1、根据动弹性波试验获取已浇筑混凝土试件的多孔骨架体积压缩模量,并计算得到Biot系数;
S2、根据CT扫描和图像识别技术,获取混凝土试件孔结构分布参数和整体孔隙率;
S3、根据所述孔结构分布参数,计算得到混凝土试件冻结过程中温度-冻结速率的关系曲线;
S4、根据Young-Laplace方程冰-未冻水平衡以及冰-吸附水膜平衡关系,构建平均孔隙压力与温度、孔结构之间的函数关系;
S5、基于多孔介质力学、多孔介质内水分迁移达西定律和多孔体系热传导Friour定律,构建多孔介质应力场-温度场-渗流场的耦合控制方程,得到冻融过程中无砟轨道内结构受力、热量传输和水分迁移之间的相互作用关系。
以下对上述步骤进行详细说明
步骤S1、根据动弹性波试验获取已浇筑混凝土试件的多孔骨架体积压缩模量,并计算得到Biot系数b为:
其中,K0为多孔骨架压缩模量,Km为固体基质压缩模量。
步骤S3、计算得到混凝土试件冻结过程中温度-冻结速率的关系曲线,其具体包括:
根据Gibbs-Thomson方程,孔径与结晶温度的关系表示为:
其中,Tf为孔隙水冰点(K),T0为正常水冰点(273.15K),γls为冰与水表面张力(39×10-3N/m),ρw为水的密度(1000kg/m3),L为水的相变潜热(333.5kJ/kg),Req为温度降至Tf时的冰晶体尖端平衡半径(nm)。
在一般环境下,
将其进行Tylor展开,可简化为:
其中,θ-孔隙水冰点(℃),带入相关数值后可得:
其中,64为计算常数(℃·nm)。
在已结冰的孔隙内,孔隙壁与冰晶之间式中存在一层水膜,孔隙水膜厚度可表示为:
其中,1.93为计算常数℃1/3·nm,因此,温度为θ时,孔径大于Rpeq(nm)的孔隙才会结冰:
Rpeq=Req+δ
假设孔隙中只有水和冰,混凝土完全饱水,则
Si+Sw=1
其中,Si、Sw分别表示冰和水的饱和度,当温度为θ时,假设孔隙为半径r的理想圆柱形,吸附层所占孔隙体积含量为2δ/r,因此结晶孔隙中,吸附层的体积含量为:
孔隙内冰的体积含量Vw-i(m3/m3)表示为:
步骤S4,根据Young-Laplace方程冰-未冻水平衡以及冰-吸附水膜平衡关系,构建平均孔隙压力与温度、孔结构之间的函数关系,其具体步骤包括:
采用高斯函数表示孔结构分布:
孔隙压力,假设在混凝土内孔隙为圆柱形,冰-未冻水接触面为球面,通过Young-Laplace方程冰-未冻水平衡以及冰-吸附水膜平衡关系分别表示为:
施加在混凝土微元孔隙内的孔隙压力可表示为微元内结晶孔隙内吸附水膜压力以及未冻结孔隙内的水压力:
步骤S5、基于多孔介质力学、多孔介质内水分迁移达西定律和多孔体系热传导Friour定律,构建多孔介质应力场-温度场-渗流场的耦合控制方程,得到冻融过程中无砟轨道内结构受力、热量传输和水分迁移之间的相互作用关系,其具体步骤包括:
S5.1、应力场控制方程;
对于各项同性的多孔体系,其应力可分为两部分:多孔体系骨架的有效应力和孔内压力:
σ=σ*-bP*I
σ*=C0:εθ
其中,σ*为骨架有效应力,P*为孔内压力,I为单位张量,εe多孔介质弹性应变,C0为多孔介质骨架的弹性刚度张量,b为Boit系数表征孔隙压力对变形的贡献,对于各项同性材料,一般表示为:
其中,K0和Km分别表示多孔骨架以及固体基质压缩模量。
假设骨架没有塑形应变,考虑到热膨胀效应,多孔骨架的应变可分为温度应变εth以及弹性应变εe:
ε=εe+εth
而骨架的温度应变可表示为:
联立式可得:
由于温度变化是孔隙压力变化的直接原因,令
则:
当不考虑体力的作用下,多孔介质微元的力学平衡方程可表示为:
S5.2、渗流场控制方程;
假设多孔体系内水封迁移服从Darcy渗流定律,由质量守恒及冰、水、骨架本构方程,混凝土冻结过程中的内水分迁移表示为:
其中,k为渗透系数(m2),μ为水的动力粘滞系数(Pa/s);εv为体应变;Kw、Ki分别为水和冰的压缩模量(Pa);αw、ai、α0分别为水、冰以及多孔骨架的体膨胀系数。
S5.3、温度场控制方程;
忽略对流对热传导的影响,并考虑相变潜热,则多孔介质的热传导方程表示为:
其中,Cef、Cm、Cw和Ci分别表示多孔体系、骨架、水以及冰的比热容(J/(kg·K));λef、λm、λm和λi分别表示多孔体系、骨架、水以及冰的热传导系数(W/(m·K));L为水的相变潜热(J/kg)。
根据本申请的一个实施例,本发明基于构建的多孔介质应力场-温度场-渗流场的耦合控制方程,构建I型双块式无砟轨道水-热-力耦合平面模型,模拟大气温度在10h内由10℃等速率降低至-20℃然后在-20℃持续2h过程中,正常状态(h=12W/(m2·K))和做简单保温隔热处理(h=6W/(m2·K))工况下的冻结与损伤行为特征。
参考图1,图中位于上方的折线为Paste I,下方的折线为Paste II,模型中支承层和道床板分别采用Paste I和Paste II的孔结构,其中Paste I的三个孔隙分布峰为327000nm、15.49nm、1.03nm,Paste II的三个孔隙分布峰为11.35nm、6.99nm、3.24nm,因此Paste I孔结构更分散,大孔更多,孔隙率更大。
参考图2,模型中各设置5个数值监测点,下文均表示为点i(i=1~10)。
数值计算时,模型中采用的参数详见下表2-1:
表2-1数值模拟参数
本发明通过有限元方法对建立的混凝土冻结与损伤行为水-热-力耦合理论模型进行求解。计算通过模型验证,得到了整个冻结过程中不同温度状态下的混凝土的应变,并于试验结果相对比。
计算结果并对其分析
1、水压力分布特征
混凝土冻结过程中,未结冰孔隙水迁移滞阻将产生孔隙内部水压力,较大的孔隙内压力将引起混凝土内部微元解体。因此,孔隙水压力是研究混凝土冻结与损伤行为的关键指标。不同工况下无砟轨道内部孔内水压力分布如图5所示。
由图5可知,最大水压力在道床板板角向内部一定深度处,且道床板水压力远大于支承层内部水压力,支承层和道床板各观测点在冻结过程中,水压力的演变规律如图。
由图6和图7可知,由于冻结阻滞效应,轨道板冻结区渗透系数迅速降低,随深度增加水压力有巨大的提升。道床板的水压力增长峰值分别在-5℃及-10℃附近,其分别对应的孔分布峰为11.35nm、6.99nm。支承层的水压力增长峰值在0℃和-5℃附近,其对应的孔分布峰为327000nm、15.49nm。相较于道床板,支承层内水压力随深度增大明显提早,由点1处的0℃提早到点6处的3℃,这并不是由于冻结温度提高,而是支承层渗透系数较小冻结区产生的水压力快速传递到轨道内部引起的。
支承层冻结区域内各点的水压力在h=6W/(m2·K)工况下相较于h=12W/(m2·K)约有0.5Mpa的减小;而道床板冻结区域内各点的水压力在h=6W/(m2·K)工况下相较于h=12W/(m2·K)也约有1Mpa的减小。因此,保温隔热处理可以降低轨道内部水压力,最主要的是有效延缓冻结过程,减少无砟轨道的冻结深度,从而有效降低无砟轨道内的水压力空间分布大小。
由图可知,道床板在-10~15℃区间内结晶压力有增长较快。支承层在0~-5℃区间内结晶压增长较为缓慢且明显小于道床板,-5℃以后开始迅速增大。
由图9可知,在冻结过程中,结冰速率占据水压力源的主导项,尤其是在冻结初期,随着温度逐渐降低,结冰速率所占比重逐渐降低,但仍然占据主导地位。因此,有必要研究支承层和道床板在冻结过程中结冰速率的演化规律。支承层和道床板各观测点在冻结过程中,结冰速率的演变规律如图10及图11所示。
由图10及图11可知,随深度增大,点内结冰速率逐渐减小,结合图6、图7中水压力随冻结深度的增加而变大的关系,揭示出冰阻导致渗透系数减小对水压力的巨大影响。在冻结过程中,在0℃很小温度区间内结冰速率最大。此后支承层在-3℃附近有一个冻结峰,道床板在-6℃和-10℃同样存在冻结峰。在冻结峰内道床板和支承层最大结冰速率分别为0.0035、0.013,结合图6、图7可知道床板水压力更大,由此也反映出渗透系数对水压力的影响大于结冰速率。支承层和道床板在h=6W/(m2·K)工况下冻结区域内的结冰速率约为h=12W/(m2·K)的工况的50%,因此保温隔热处理可有效减小轨道内结冰速度率。
由图12可见支承层内孔隙水在-10℃后结晶速率迅速降低,最终冰饱和度约为0.35,而道床板由于孔结构较为致密,在-10~-15℃区间内冰饱和度仍有显著增大,其最终冰饱和度约为0.55的试件,这是结冰速度更快。
2、温度分布特征
温度是引起混凝土冻结与损伤行为的最直接因素,温度分布规律可以间接判断混凝土的冻结与损伤行为深度。对于不同的表面热传导系数条件下,无砟轨道的最终温度场如图13所示。
由图13可知,在道床板和支承层板角处冻结最深,保温隔热处理后轨道表面的温度略有提高。h=6W/(m2·K)和h=12W/(m2·K)工况下轨道表面温度分别为-10℃和-15℃。其中各监测点的温度变化速率分别如图14及图15所示。
由图14及图15可知,随着深度的增加,轨道内降温速率也越来越小。由于孔隙内水化热的释放,在0℃附近,变温速率有非常大的波动。同时支承层和道床分别在-3℃和-10℃附近区间内降温速率变小,分别对应于图10和图11的结晶峰。保温隔热处理后轨道内部降温速率有显著的减小。
3、应力分布特征
冻结过程中,混凝土孔内压力作用于固体骨架,为分析轨道内的应力分布特征,不同工况下无砟轨道的Tresca应力由图16所示。
由图16可知,由于支承层弹模较低,密实度较小,Boit系数更大,孔内水压力更多的反映到混凝土骨架上,道床板内最大剪应力出现在道床板表层和深度约为5cm处。该结果很好的解释了现场双块式无砟轨道表面冻结与损伤行为剥落情况。而支承层的最大Tresca应力则出现在板角表面和距板角约5cm处。各观测点在不同工况下的Tresca应力演化规律分别如图17及图18所示。
由图17及图18可知,在0℃以下时,随着温度的降低支承层内Tresca应力显著增大,距离表面越深Tresca应力增长越快。h=12W/(m2·K)和h=6W/(m2·K)工况下,支承层最大Tresca应力分别为5.5Mpa和2Mpa,而道床板最大Tresca应力分别为5Mpa和4Mpa,保温隔热处理可减弱冻结与损伤行为对无砟轨道的破坏力度。由于支承层内水压力传递速度较快,Tresca应力也体现出明显的提前性。
根据以上描述,综合可得:
严寒地区无砟轨道受冻结与损伤行为影响较大,长期服役过程中不断出现表面剥蚀的状况。本发明基于多孔介质力学,建立了I型双块式无砟轨道水-热-力耦合模型,分析了冻结作用下无砟轨道应力-应变分布特征、温度以及孔隙压力的分布状态。得出主要结论如下:
(1)结冰速率占据混凝土内冻结水压力源的主导地位,由于水分结冰的阻止效应致使渗透系数降低,距离轨道表面越深的观测点在相同冻结温度下水压力峰值越大。
(2)混凝土的孔结构显著影响轨道内部水分的结冰速率及水压力演化规律。道床板孔径主要集中于3~15nm,冻结峰温度为0℃、-6℃及-11℃,结冰速率较为平缓。而支承层内部大孔结构较多,孔结构分布峰值对应的孔径为300000nm、15nm及1nm,结晶峰温度为0℃及-4℃,冻结速率更快。
(3)支承层点1在-5℃结冰峰结冰速率为0.003kg/(m3·s),水压力为1Mpa。而道床板孔隙率小,渗透系数为支承层的1/30,点6在-12℃结冰峰的结冰速率为0.0007kg/(m3·s),水压力为8.5Mpa。孔隙率及渗透系数对水压力的影响更显著。
(4)虽然道床板水压力更大,但支承层孔隙率更高,Boit系数更大,孔内水压力更多的反映到混凝土骨架上,应力更大,受冻结与损伤行为最大。
(5)保温隔热层处理可有效减小轨道的冻结速率及深度,在h=6W/(m2·K)和h=12W/(m2·K)工况下,轨道表面温度分别为-14℃、-18℃,最大水压力分别为6Mpa、2Mpa。因此可通过无砟轨道两侧堆道砟的方法减轻无砟轨道的冻结与损伤行为。
虽然结合附图对发明的具体实施方式进行了详细地描述,但不应理解为对本专利的保护范围的限定。在权利要求书所描述的范围内,本领域技术人员不经创造性劳动即可做出的各种修改和变形仍属本专利的保护范围。
Claims (5)
1.一种无砟轨道冻结与损伤行为计算方法,其特征在于,包括:
S1、根据动弹性波试验获取已浇筑混凝土试件的多孔骨架体积压缩模量,并计算得到Biot系数;
S2、根据CT扫描和图像识别技术,获取混凝土试件孔结构分布参数和整体孔隙率;
S3、根据所述孔结构分布参数,计算得到混凝土试件冻结过程中温度-冻结速率的关系曲线;
S4、根据Young-Laplace方程冰-未冻水平衡以及冰-吸附水膜平衡关系,构建平均孔隙压力与温度、孔结构之间的函数关系;
S5、基于多孔介质力学、多孔介质内水分迁移达西定律和多孔体系热传导Fourier定律,构建多孔介质应力场-温度场-渗流场的耦合控制方程,得到冻融过程中无砟轨道内结构受力、热量传输和水分迁移之间的相互作用关系。
3.根据权利要求2所述的无砟轨道冻结与损伤行为计算方法,其特征在于,所述步骤S3中计算得到混凝土试件冻结过程中温度-冻结速率的关系曲线为:
Vw-i为孔隙内冰的体积含量,其计算表达式为:
Vw-i=φ(Rpeq)-Vads(Rpeq)
吸附层的体积含量Vads为:
孔径Rpeq为:
Rpeq=Req+δ
其中,Req为冰晶体尖端平衡半径;δ为孔隙水膜厚度之和,其表达式为:
其中,θ为摄氏温度;
冰晶体尖端平衡半径Req为:
其中,γls为冰与水表面张力,T0为正常水冰点,Rmax为冰晶体尖端最大半径,L为水的相变潜热,即单位质量的水分结冰释放的热量。
5.根据权利要求4所述的无砟轨道冻结与损伤行为计算方法,其特征在于,所述步骤S5中构建多孔介质应力场-温度场-渗流场的耦合控制方程分别为:
当不考虑体力的作用下,多孔介质微元的应力场控制方程为:
其中,K0为多孔骨架压缩模量,α0为多孔骨架的体膨胀系数;
忽略对流对热传导的影响,并考虑相变潜热,则多孔介质的温度场控制方程为:
其中,Cm、Cw和Cw分别表示多骨架、水以及冰的比热容,λef、λm、λw、λi分别为多孔体系、骨架、水以及冰的热传导系数,Sw为水饱和度,ρw为水密度,Si为冰饱和度,ρi为冰密度;
混凝土冻结过程中的内水分迁移,即渗流场的控制方程为:
其中,β为多孔弹性储水修正系数,Kw为水的压缩模量,Ki为冰的压缩模量;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910698540.4A CN110390176B (zh) | 2019-07-31 | 2019-07-31 | 一种无砟轨道冻结与损伤行为计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910698540.4A CN110390176B (zh) | 2019-07-31 | 2019-07-31 | 一种无砟轨道冻结与损伤行为计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110390176A CN110390176A (zh) | 2019-10-29 |
CN110390176B true CN110390176B (zh) | 2020-05-19 |
Family
ID=68287818
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910698540.4A Active CN110390176B (zh) | 2019-07-31 | 2019-07-31 | 一种无砟轨道冻结与损伤行为计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110390176B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111125922B (zh) * | 2019-12-28 | 2024-02-20 | 中汽研汽车检验中心(天津)有限公司 | 一种冷空间内下落液滴冻结过程的计算方法 |
CN112487611B (zh) * | 2020-10-27 | 2022-05-17 | 安徽建筑大学 | 上覆压力作用下冻结土体水分迁移模型的构建方法 |
CN112784407B (zh) * | 2021-01-11 | 2022-05-17 | 西南交通大学 | 一种考虑沥青热可逆老化现象的胶结料温度应力计算方法 |
CN113984591B (zh) * | 2021-09-03 | 2023-09-22 | 南京大学 | 一种寒冷区域内多孔介质中lnapl迁移模拟方法 |
CN114428987B (zh) * | 2021-12-24 | 2024-06-14 | 中国水电建设集团十五工程局有限公司 | 多物理场作用下的混凝土重力坝坝体应力分析方法 |
CN114441329A (zh) * | 2022-01-28 | 2022-05-06 | 中国水利水电科学研究院 | 多孔连续介质体的孔隙水压力系数测定方法 |
CN114491781B (zh) * | 2022-04-08 | 2022-07-01 | 西南交通大学 | 一种考虑不确定性的无砟轨道基床损伤预警方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108399308A (zh) * | 2018-03-14 | 2018-08-14 | 武汉理工大学 | 沥青板式无砟轨道轨下基础结构及其准静态设计方法 |
CN108645885A (zh) * | 2018-05-25 | 2018-10-12 | 吉林大学 | 大尺度土体水-热-力-盐四场耦合作用试验系统及方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20140154808A1 (en) * | 2012-12-03 | 2014-06-05 | Gordhanbhai N. Patel | Monitoring system based on etching of metals |
CN105160138B (zh) * | 2015-10-12 | 2018-05-22 | 中铁建大桥工程局集团第一工程有限公司 | 隧道流固耦合稳定性分析方法 |
CN107609304B (zh) * | 2017-09-29 | 2020-10-27 | 中国铁道科学研究院铁道建筑研究所 | 大跨度铁路桥梁的基于phm的故障诊断预测系统及方法 |
CN109063257B (zh) * | 2018-07-02 | 2019-04-26 | 山东科技大学 | 一种煤岩体分区注水渗流-损伤-应力耦合数值模拟方法 |
CN109214121B (zh) * | 2018-10-17 | 2022-01-25 | 哈尔滨工业大学 | 钢筋混凝土结构腐蚀电驱动统一控制方法 |
-
2019
- 2019-07-31 CN CN201910698540.4A patent/CN110390176B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108399308A (zh) * | 2018-03-14 | 2018-08-14 | 武汉理工大学 | 沥青板式无砟轨道轨下基础结构及其准静态设计方法 |
CN108645885A (zh) * | 2018-05-25 | 2018-10-12 | 吉林大学 | 大尺度土体水-热-力-盐四场耦合作用试验系统及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110390176A (zh) | 2019-10-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110390176B (zh) | 一种无砟轨道冻结与损伤行为计算方法 | |
CN109283215B (zh) | 反复冻融损伤后寒区隧道围岩冻胀力的测算方法 | |
Gong et al. | Modeling of water transport in highly saturated concrete with wet surface during freeze/thaw | |
CN108645885A (zh) | 大尺度土体水-热-力-盐四场耦合作用试验系统及方法 | |
CN100554598C (zh) | 混凝土浇筑中温控防裂的施工方法 | |
CN109283086A (zh) | 精准测算反复冻融损伤后寒区隧道岩石弹性模量的方法 | |
Liu et al. | Shotcrete performance-loss due to seepage and temperature coupling in cold-region tunnels | |
Yuan et al. | Study on frost resistance of EPS concrete based on EPS beads wrapping modification | |
Li et al. | Mechanical properties of frozen ballast aggregates with different ice contents and temperatures | |
CN113202480B (zh) | 一种联合承载的双层井壁及其施工方法 | |
CN110219457B (zh) | 一种钢筋模板一体化墙体模板断桥对拉连接杆件 | |
CN107389433B (zh) | 一种用于检测水库结冰对水利建筑冰推力以及冰拔破坏的检测箱及其检测方法 | |
Thomachot et al. | Frost damage of bricks composing a railway tunnel monument in Central Japan: field monitoring and laboratory simulation | |
CN116122307B (zh) | 边坡应急多级加固方法 | |
Dashjamts et al. | Research on consolidation of frozen soils upon thawing | |
Adkins et al. | Freeze-thaw deterioration of concrete pavements | |
CN202689037U (zh) | 一种底板层结构 | |
CN203569726U (zh) | 一种承重型现浇保温层外墙体 | |
Du et al. | Crack control technologies of mass concrete | |
CN219809014U (zh) | 一种寒区隧道自保温装配式二衬结构 | |
Li et al. | Influence of freezing rate on cryo-damage of cementitious material | |
CN216810402U (zh) | 一种混凝土空心砌块墙与清水混凝土外墙锚固结构 | |
Gao et al. | Temperature cracks controlling technology for the mass concrete in solid sections of lower pylon column | |
CN114462251A (zh) | 基于冻融滞后性的分凝冻胀控制方法 | |
Badiane et al. | Monitoring of flexible pavement structures during freezing and thawing |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |