水热力耦合下公路隧道结构各向异性冻胀力分布特征研究

寒区围岩的冻胀会对衬砌结构带来附加荷载,威胁着衬砌结构的耐久性与安全性。为了研究季节冻土区隧道冻胀力的分布特征,基于能量守恒与质量守恒基本原理,建立了考虑正交各向异性冻胀变形的正冻围岩水-热-力耦合模型;结合青沙山隧道冻胀力的监测结果,验证了隧道冻胀力模型的可靠性。进而,针对东天山特长隧道的抗冻需求,构建了相应的隧道冻胀数值模型,对东天山隧道的温度场、水分场与冻胀力分布特征进行了研究,并分析了最低气温、初始地层含水率、已冻与未冻围岩模量比、冻胀变形各向异性系数对冻胀力分布的影响。东天山隧道冻胀的数值仿真结果表明:隧道的冻结深度分布并不均匀,在拱脚处最小,仰拱中心处最大,两处冻结深度相差48 cm。冻结范围的差异分布主要为拱脚处几何曲率较大,冷量需辐射至更广的空间所致。同时,由于拱脚处冻深最小,且拱脚的几何曲率较顶拱与仰拱更大,导致衬砌的拱脚处弯曲折叠最大,von Mises应力最大。单个冻融周期内,隧道冻深范围内围岩的含水率可分为冻结、融化、滞水、滞水消散共4个阶段。随着冻融循环次数的增加,拱顶与拱侧的滞水阶段体积含水率分别升高了10.46%与4.21%,拱脚与衬砌底部的滞水阶段体积含水率略有降低。围岩的冻胀对衬砌同时产生了法向应力与切向应力,其中顶拱、仰拱主要体现为压应力,拱脚的压应力较小且局部体现为拉应力。不同最低气温、初始地层含水率、冻结与未冻结围岩模量比、冻胀变形正交各向异性系数下的冻胀力分布模式相同。其中衬砌外围的法向应力整体呈“蘑菇形”。最低气温的降低、冻胀变形正交各向异性系数的增大,分别通过增大冻结范围与促进冻胀应变方向集中化,最终导致隧道冻胀力整体数值的增大,两者对隧道冻胀力影响显著。冻结与未冻结围岩模量比、初始地层含水率对隧道冻胀力存在一定程度的影响:冻结与未冻结围岩模量比与隧道冻胀力负相关,初始地层含水率与隧道冻胀力正相关。冻胀变形各向异性系数对冻胀力的大小与分布影响最大。同时,经过多年的冻融循环,不同初始地层含水率下衬砌外围的水分场差距缩小,导致初始地层含水率对第20年冻胀力的大小与分布影响最小。整体而言,隧道冻胀力的分布主要是温度场与衬砌几何之间博弈的结果:其中拱脚处最小的冻结深度导致该处冻胀力较小,但衬砌几何受力后顶拱与仰拱变形导致的拱脚弯曲折叠并向外侧围岩挤压,会增大拱脚处的压应力。因此,对于衬砌厚度较小的隧道而言,拱脚向外侧围岩挤压的效应更为明显,这导致拱脚处冻胀力较大。综合来看,寒区隧道的冻胀力分布需同时考虑温度条件、水分条件、围岩冻胀变形的各向异性、衬砌几何等因素的影响。