本发明涉及冻土物理学，热力学和流体动力学领域，特别涉及一种上覆压力作用下冻结土体水分迁移模型的构建方法。
背景技术：
：上覆压力对土(岩)体冻胀有一定的抑制作用，目前主要反映在降低冻结温度和影响水分重分布(驱动力)两方面：1)上覆压力对土体冻结温度的影响随着上覆压力的增加，导致矿物颗粒之间的接触应力增大，降低的了土体的冻结温度。《冻土未冻水含量与压力关系的实验研究》(张立新等，《冰川冻土》，第124-127页，1998，20(2))公开了：冻结温度随压力增大呈线性降低，斜率接近于-0.075℃/mpa，因此上覆压力影响了土体中水分的相态转换，在同一土体中，对冻土施加的上覆压力越大，未冻水含量增加的越多，且可用广义clapeyron方程的固-液相平衡方程来确定上覆压力与冻结温度的降低量。此外，《phenomenonandmechanismoffrostheaving》(kaplar[j].highwayresearchrecord，第1-13页1970，304，)认为这种外力作用，实质上是对土体系统做功，而土体克服上覆压力所做的功等于活性薄膜处激发的自由能，因此压力作用下的冻胀是能量的一种转换形式。2)上覆压力对土体水分重分布(驱动力)的影响上覆压力会影响水分迁移的驱动力，进而减少冻土中水分向冻结锋面的迁移量。关于荷载对这种驱动力影响，目前有两种意见：一种是从能量的角度解释了上覆压力对水分迁移的抑制机制，例如：《thethermodynamicsoffrostdamagetoporoussolids》(everett，[j].physicalandinorganicchemistry，第1541–1551页，1961，57)认为水分迁移的驱动力是土颗粒中冰-水之间的毛细吸力，当上覆压力增加到等于冰-水界面所产生的界面能量时，冻结面就不吸水，土体的冻胀则停止，这时上覆压力被称为“中断压力”，且该“中断压力”与温度条件无关。一种是从压力的角度解释了上覆压力对水分迁移的抑制机制，例如：《theunfrozenwaterandtheapparentspecificheatcapacityoffrozensoils》(anderson，[j].prceedingsof2thinternationalconferenceonpermafrost，第512-521页，1973)认为水分迁移的驱动力是一种负压力，由于冻土中未冻水比自由水具有较低的孔隙水压，且孔隙水压是温度的函数，因此当自由水压为零时，对应的孔隙水压为负值；当施加上覆压力时，孔隙水压由负转正、逐渐增大时，此时迁移驱动力逐渐消失，水分迁移趋近停止。但是，根据《压力作用下冻结粉质粘土的水分迁移特征》(明锋，2015)主导的冻胀试验发现，当上覆压力为1mpa时，土体的冻胀率接近于零，水分迁移趋向停止。但此时未冻水的冻结温度只降低-0.075℃，对未冻水含量变化的影响很小，几乎可以忽略不计。显然，上覆压力导致冻结温度降低不是抑制水分迁移的主要原因。此外，根据《studyonthefreezingcharacteristicsofsiltyclayunderhighloadingconditions》(guan，[j].coldregionsscienceandtechnology，第26-31页，2015，110)主导的高应力水平条件下的一维冻胀加载试验，发现上覆压力为6mpa时，仍检测到明显的水分再分配现象，但此时的上覆压力远远大于毛细吸力或者孔隙水负压；显然，上覆压力导致驱动力降低也不是抑制水分迁移的主要原因。综上可知，上述两种观点都不是上覆压力抑制冻胀的主要原因。目前关于上覆压力对冻胀影响的研究还处于理论探索阶段，且尚未能有试验进行验证。因此，探究上覆压力对水分迁移的抑制机理具有重要的理论价值和现实意义。本发明依托安徽省自然科学基金资助项目(2008085me165)，对上覆压力作用下冻结土体水分迁移特征进行研究。技术实现要素：本发明的目的在于提供一种上覆压力作用下冻结土体水分迁移模型的构建方法，首次从上覆压力控制冰的分凝角度，对水分迁移的抑制机理展开分析。由于上覆压力对分凝-冻结温度产生直接影响，而分凝-冻结温度又直接控制冰透镜体暖端的表面吸附力和冻结缘渗透系数，因此，可认为上覆压力通过控制冰的分凝，从而抑制水分迁移的速度。为了实现上述目的，本发明提供如下技术方案：一种上覆压力作用下冻结土体水分迁移模型的构建方法，构建方法包括如下步骤：步骤1：获取平衡状态下的理论冰压方程；步骤2：获取冰透镜体暖端的分凝-冻结温度方程；步骤3：根据上覆压力、步骤1和步骤2中得到的理论冰压方程、分凝-冻结温度方程获取冰透镜体暖端的迁移驱动力方程；步骤4：获得冰透镜体暖端的表面吸附力方程；步骤5：获取冻结缘区的总渗透系数方程；步骤6：获取冰透镜体暖端的水分迁移速度显式方程，完成上覆压力作用下冻结土体水分迁移模型的构建。进一步地，在上述的上覆压力作用下冻结土体水分迁移模型的构建方法中，在所述步骤1中，根据广义clapeyron方程，可得平衡状态下的理论冰压方程：式(5)中：ps0为平衡状态下薄膜水受到的理论冰压；l为水的相变潜热，可取3.34×109cm2/s2；vs为冰的比容；ta为纯冰的绝对冻结温度，可取273.15k；t为冰-水界面的温度，t为变量。进一步地，在上述的上覆压力作用下冻结土体水分迁移模型的构建方法中，在所述步骤2中，具体包括：步骤2.1：假定冰透镜体形成之前实际冰压等于理论冰压；步骤2.2：当实际冰压等于上覆压力与分离压力之和时，则会产生新的分凝冰，获取分凝冰形成条件方程；步骤2.3：根据分凝冰形成条件方程获取冰透镜体暖端的分凝-冻结温度方程；实际冰压方程为：结合分凝冰形成机制可知，当实际冰压等于上覆压力与分离压力之和时，则会产生新的分凝冰，即：ps＝pob psep(3)分凝冰形成条件方程的获取：式(2)和(3)中：ps0为平衡状态下薄膜水受到的理论冰压；ps为实际冰压；pob为上覆压力；psep为分离压力；结合公式(2)、(3)，可得分凝冰形成条件方程：变换公式(3)可得分凝-冻结温度方程：式(4)和(5)中：ts为冰透镜体暖端的分凝-冻结温度。进一步地，在上述的上覆压力作用下冻结土体水分迁移模型的构建方法中，在所述步骤3中，具体包括：根据程桦提出的薄膜水迁移驱动力模型，可知冰透镜体暖端的迁移驱动力方程等于冰透镜体暖端的理论冰压与上覆压力之差：式(6)中：pld为迁移驱动力，ts代表t在冰透镜体暖端的分凝-冻结温度。进一步地，在上述的上覆压力作用下冻结土体水分迁移模型的构建方法中，在所述步骤4中，具体包括：步骤4.1：获取压-吸变量热力学换算系数λ，λ＝-vs/vl＝-1.09，vs为冰的比容，vl是水的比容，λ取-1.09；步骤4.2：根据迁移驱动力方程和热力学换算系数λ，获取冰透镜体暖端的表面吸附力方程；式(7)中：plb为表面吸附力，pld为迁移驱动力，λ为压-吸变量的热力学换算系数。进一步地，在上述的上覆压力作用下冻结土体水分迁移模型的构建方法中，在所述步骤5中：根据托马斯统计的试验结果可知，未冻结区、冻结缘区和已冻结区之间的渗透系数不同，因此设定未冻结区的渗透系数kuf，kuf也为饱和土体在常温条件下的渗透系数，冻结缘区的总渗透系数kff，已冻结区的渗透系数kf，托马斯推荐的冻结缘渗透系数方程为：kff＝kuf[1-(ts-tf)]-β(8)根据所述步骤2得到的分凝-冻结温度方程、公式(8)获取冻结缘区的总渗透系数方程，冻结缘区的总渗透系数方程为：式(8)和(9)中：β为幂，无量纲；tf为冻结锋面的冻结温度，可取0℃；kuf为饱和土体在常温条件下的渗透系数；kff为冻结缘区的总渗透系数。进一步地，在上述的上覆压力作用下冻结土体水分迁移模型的构建方法中，在所述步骤6中，根据达西定律，结合冰透镜体暖端的表面吸附力plb、冻结缘区的总渗透系数kff，获取冰透镜体暖端的水分迁移速度显式方程，完成压力作用下冻结土体水分迁移模型的构建，具体包括：假设冻结缘区内的水分迁移速度符合达西定律，结合公式(7)、(9)可得水分迁移速度显式方程：式(10)中，lt为水分迁移的渗透路径、ρw为水的密度、g为重力加速度、vff为冰透镜体暖端的水分迁移速度。进一步地，在上述的上覆压力作用下冻结土体水分迁移模型的构建方法中，利用冻结土体的试样在上覆压力作用下进行水分迁移模型的构建，所述冻结土体的试样包括未冻结区、冻结缘区和已冻结区，所述未冻结区的一端与所述冻结缘区的一端连接，所述冻结缘区的另一端与所述已冻结区的一端连接，所述未冻结区、所述冻结缘区和所述已冻结区依次排列设置，在所述已冻结区内并靠近所述冻结缘区的一侧形成冰透镜体，在所述未冻结区内靠近所述冻结缘区的一侧形成冻结锋面。进一步地，在上述的上覆压力作用下冻结土体水分迁移模型的构建方法中，所述未冻结区的另一端为热端，所述已冻结区的另一端为冷端，所述冰透镜体靠近冻结缘区的一端为所述冰透镜体的暖端。进一步地，在上述的上覆压力作用下冻结土体水分迁移模型的构建方法中，所述冻结缘区内设置有若干土颗粒，每个所述土颗粒外围均包裹有一层未冻水膜，水分可通过所述未冻水膜从所述未冻结区向所述冰透镜体的暖端迁移。分析可知，本发明公开一种上覆压力作用下冻结土体水分迁移模型的构建方法，首次从上覆压力控制冰的分凝角度，对水分迁移的抑制机理展开分析。由于上覆压力对分凝-冻结温度产生直接影响，而分凝-冻结温度又直接控制冰透镜体暖端的表面吸附力和冻结缘渗透系数，因此，可认为上覆压力通过控制冰的分凝，从而抑制水分迁移的速度，达到抑制水分迁移的目的。附图说明构成本申请的一部分的说明书附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。其中：图1为压力作用下冻结土体的试样示意图。图2为压力作用下水分迁移模型的示意图。图3为水分迁移的计算速度与试验速度对比图。附图标记说明：1未冻结区、2冻结缘区、3冰透镜体、4已冻结区、5冻结锋面、6土颗粒、h1未冻结区的高度、h2冻结缘区的厚度、h3已冻结区的高度、h4渗透路径长度、h5试样的初始高度。具体实施方式下面将参考附图并结合实施例来详细说明本发明。各个示例通过本发明的解释的方式提供而非限制本发明。实际上，本领域的技术人员将清楚，在不脱离本发明的范围或精神的情况下，可在本发明中进行修改和变型。例如，示为或描述为一个实施例的一部分的特征可用于另一个实施例，以产生又一个实施例。因此，所期望的是，本发明包含归入所附权利要求及其等同物的范围内的此类修改和变型。在本发明的描述中，术语“纵向”、“横向”、“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明而不是要求本发明必须以特定的方位构造和操作，因此不能理解为对本发明的限制。本发明中使用的术语“相连”、“连接”、“设置”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接；可以是直接相连，也可以通过中间部件间接相连；可以是有线电连接、无线电连接，也可以是无线通信信号连接，对于本领域的普通技术人员而言，可以根据具体情况理解上述术语的具体含义。所附附图中示出了本发明的一个或多个示例。详细描述使用了数字和字母标记来指代附图中的特征。附图和描述中的相似或类似标记的已经用于指代本发明的相似或类似的部分。如本文所用的那样，用语“第一”、“第二”、“第三”以及“第四”可互换地使用，以将一个构件与另一个区分开，且不旨在表示单独构件的位置或重要性。如图1至图3所示，根据本发明的实施例，提供了一种上覆压力作用下冻结土体水分迁移模型的构建方法。首先，准备冻结土体的试样利用冻结土体的试样在上覆压力作用下进行水分迁移模型的构建。冻结土体的试样包括未冻结区1、冻结缘区2和已冻结区4，未冻结区1的一端与冻结缘区2的一端连接，冻结缘区2的另一端与已冻结区4的一端连接，未冻结区1、冻结缘区2和已冻结区4依次排列设置。通过未冻结区1、冻结缘区2和已冻结区4模拟自然界冻土的基本结构，在已冻结区4内并靠近冻结缘区2的一侧形成冰透镜体3，在未冻结区1内靠近冻结缘区2的一侧形成冻结锋面5。未冻结区1的另一端为热端，已冻结区4的另一端为冷端。冰透镜体3靠近冻结缘区2的一端为冰透镜体3的暖端。冻结缘区2内设置有若干土颗粒6，每个土颗粒6外围均包裹有一层厚度约为50nm的未冻水膜(也可称为薄膜水)，水分可通过未冻水膜从未冻结区1向冰透镜体3的暖端迁移。向已冻结区4的另一端施加朝向冻结缘区2的上覆压力，向未冻结区1的另一端施加朝向冻结缘区2的外部水压力，由于未冻结区1、冻结缘区2和已冻结区4之间的冻结状态不同，因此设定未冻结区1的水分迁移速度为vuf，冻结缘区2的水分迁移速度为vff，已冻结区4的水分迁移速度为vf。施加的上覆压力增大时，冰透镜体3的形成位置会向上移动，即上覆压力越大，冰透镜体3的位置会更靠近冷端。上覆压力i(pobi)，上覆压力ii(pobii)分别代表上覆压力的两种情况，假设上覆压力i为50kpa，上覆压力ii为200kpa。在上覆压力ii作用下形成的冰透镜体3的位置位于在上覆压力i作用下形成的冰透镜体3的上方。本发明讨论了上覆压力增大导致迁移驱动力和冻结缘总渗透系数的变化。如图2所示，上覆压力i(pobi)作用下：冰透镜体3暖端的分凝-冻结温度为ts1、冻结缘区2的总渗透系数为kff1、冰透镜体暖端的水分迁移速度为vff1；上覆压力ii(pobii)作用下：冰透镜体3暖端的分凝-冻结温度为ts2、冻结缘区2的总渗透系数为kff2、冰透镜体暖端的水分迁移速度为vff2。未冻结区1的高度为h1，冻结缘区2的厚度为h2，已冻结区4的高度(包含冰透镜体3的厚度)为h3，渗透路径长度(包含未冻结区1的高度与冻结缘区2的厚度)为h4，h4＝h1 h2，试样的初始高度为h5，h5＝h1 h2 h3。选择试样为高度为11cm，直径为10cm的标准圆柱体。且在试验开始前，对饱和粉土的土体样品施加1mpa的压力并固结24h，以保证土体在不同上覆压力条件下物理力学性质的恒定。试验期间，试样放入环境箱内，冷端(顶端)循环液的温度tc设置为-10℃，热端(底端)循环液的温度tw设置为 1℃，环境箱温度也设置为 1℃，且温度的波动范围控制在±0.1℃。与此同时，沿试样高度从下而上在2.5cm、5.0cm、7.5cm和10.0cm四个层位点(测点)均匀布置4个温度、水分集成传感器，其传感器型号为mec10，含水量可选量程为：0～50％和0～100％；分辨率在0～50％内为0.03％，在50～100％内为1％；精度在0～50％内为2％，在50～100％内为2％。温度量程为-40～80℃，分辨率为0.1℃，精度为±0.5℃；以监测冻结过程中温度和未冻水含量的变化。其次，利用上述冻结土体的试样，构建上覆压力作用下冻结土体水分迁移模型，构建方法包括如下步骤：步骤1：获取平衡状态下的理论冰压方程沿着试样的长度方向设置x轴，试样为圆柱体时，试样的长度方向也是试样的轴线方向，根据《冻土毛细-薄膜水分迁移机制及其试验验证[j].岩土工程学报，2020，42(10)：1790-1799》，程桦(chenghua)，提出的广义克拉珀龙(clapeyron)力学模型，给出薄膜水沿x轴的迁移驱动力pld(迁移驱动力是指驱使薄膜水从热端向冷端迁移的压强差)、实际冰压、上覆压力i、上覆压力ii、理论冰压(理论冰压为土体内部温度对应的平衡状态条件下的冰压力，与负温度值成正比)，根据广义克拉珀龙(clapeyron)方程获取平衡状态下的理论冰压方程。如图1所示，理论冰压的获取：根据广义clapeyron方程，可得平衡状态下的理论冰压方程：式(1)中：ps0为平衡状态下薄膜水受到的理论冰压；l为水的相变潜热，可取3.34×109cm2/s2；vs为冰的比容；ta为纯冰的绝对冻结温度，可取273.15k；t为冰-水界面的温度(摄氏度)，t为变量。步骤2：获取冰透镜体暖端的分凝-冻结温度方程根据理论冰压方程、上覆压力、分离压力获取冰透镜体暖端的分凝-冻结温度方程。分离压力一般为20-100kpa，可通过巴西劈裂试验进行测量。具体包括：步骤2.1：根据《internalstressesinfrozenground》(williams，[j].canadiangeotechnicaljournal.1985，22:413-416)，假定冰透镜体形成之前实际冰压等于理论冰压；步骤2.2：当实际冰压等于上覆压力与分离压力之和时，则会产生新的分凝冰，获取分凝冰形成条件方程；步骤2.3：根据分凝冰形成条件方程获取冰透镜体暖端的分凝-冻结温度方程。如图1所示，由理论冰压方程(1)可知，理论冰压ps0随温度降低会逐渐增大，根据williams(1985)可知，在新的冰透镜体形成之前实际冰压等于理论冰压，实际冰压方程为：结合分凝冰形成机制可知，当实际冰压等于上覆压力与分离压力之和时，则会产生新的分凝冰，即：ps＝pob psep(3)分凝冰形成条件方程的获取：式(3)中：ps为实际冰压；pob为上覆压力；psep为分离压力。结合公式(2)、(3)，可得分凝冰形成条件方程：变换公式(4)可得分凝-冻结温度方程：式(4)和(5)中：ts为冰透镜体暖端的分凝-冻结温度。步骤3：根据步骤1和步骤2中得到的理论冰压方程、上覆压力、分凝-冻结温度方程获取冰透镜体暖端的迁移驱动力方程。根据冰透镜体暖端的理论冰压与上覆压力之差获取冰透镜体暖端的迁移驱动力方程具体包括：根据程桦(chenghua)提出的薄膜水迁移驱动力模型，可知冰透镜体暖端的迁移驱动力方程等于冰透镜体暖端的理论冰压与上覆压力之差：式(6)中：pld为迁移驱动力，ts代表t在冰透镜体暖端的分凝-冻结温度。步骤4：获得冰透镜体暖端的表面吸附力方程根据冰透镜体暖端的迁移驱动力方程和程桦(chenghua)的压-吸变量代换，获得冰透镜体暖端的表面吸附力方程。具体包括：步骤4.1：获取压-吸变量热力学换算系数λ，λ＝-vs/vl＝-1.09，vs为冰的比容，vl是水的比容，λ取-1.09；步骤4.2：根据迁移驱动力方程和热力学换算系数λ，获取冰透镜体暖端的表面吸附力方程；根据迁移驱动力方程，并根据程桦(chenghua)提出的压-吸变量代换，可得冰透镜体暖端的表面吸附力方程：式(7)中：plb为表面吸附力，pld为迁移驱动力，λ为压-吸变量的热力学换算系数。步骤5：获取冻结缘区的总渗透系数方程根据托马斯(thomas)统计的试验结果可知，未冻结区1、冻结缘区2和已冻结区4之间的渗透系数不同，因此设定未冻结区1的渗透系数kuf，kuf也为饱和土体在常温条件下的渗透系数，冻结缘区2的总渗透系数kff，已冻结区4的渗透系数kf，kf趋向于零。托马斯(thomas)推荐的冻结缘渗透系数方程为：kff＝kuf[1-(ts-tf)]-β(8)根据步骤2得到的分凝-冻结温度方程、公式(8)获取冻结缘区的总渗透系数方程。具体包括：如图1所示，通过托马斯(thomas)推荐的冻结缘渗透系数方程，可给出冻结缘区的渗透系数与分凝-冻结温度之间的关系式、即冻结缘区的总渗透系数方程为：式(8)和(9)中：β为幂，无量纲；tf为冻结锋面的冻结温度，可取0℃；kuf为饱和土体在常温条件下的渗透系数；kff为冻结缘区的总渗透系数。步骤6：获取冰透镜体暖端的水分迁移速度显式方程根据步骤4得到的冰透镜体暖端的表面吸附力方程、步骤5得到的冻结缘区的总渗透系数方程，在达西定律的基础上，获取冰透镜体暖端的水分迁移速度显式方程，完成压力作用下冻结土体水分迁移模型的构建。以根据压力作用下冻结土体水分迁移模型全面合理地解释附加荷载对土体冻胀的抑制现象。根据达西定律，结合冰透镜体3暖端的表面吸附力plb、冻结缘区的总渗透系数kff，获取冰透镜体3暖端的水分迁移速度显式方程，完成压力作用下冻结土体水分迁移模型的构建，具体包括：根据冰透镜体暖端的表面吸附力方程、冻结缘区的总渗透系数方程，结合达西定律，获取冰透镜体暖端的水分迁移速度显式方程。水分迁移速度显式方程的获取：假设冻结缘区2内的水分迁移速度符合达西定律，结合公式(7)、(9)可得：式(10)中，lt为水分迁移的渗透路径、ρw为水的密度、g为重力加速度、vff为冰透镜体暖端的水分迁移速度。本发明涉及的参数如下：t为冰-水界面的温度，t为变量，单位为℃；ta为纯冰的冻结温度(绝对温度)，单位为k，且ta＝273.15k；ps0为平衡状态下薄膜水受到的理论冰压，单位为mpa；ps为实际冰压，单位为mpa；pob为上覆压力，单位为mpa；psep为分离压力，单位为mpa；pld为迁移驱动力，单位为mpa；plb为表面吸附力，单位为mpa；λ为压-吸变量的热力学换算系数，取-1.09；kuf为饱和土体在常温条件下的渗透系数，单位为cm/s；kff为冻结缘区的总渗透系数，单位为cm/s；tf为冻结锋面的冻结温度，单位为℃；ts为冰透镜体暖端的分凝-冻结温度，单位为℃；g为重力加速度，单位为m/s2；vff为冰透镜体暖端的水分迁移速度，单位为cm/s；vs为冰的比容，单位为s2/m；vl是水的比容，单位为s2/m；l为纯水的相变潜热，单位为cm2/s2；lt为水分迁移的渗透路径，单位为cm；β为幂，无量纲；ρw为水的密度，单位为g/cm3。本发明还提出了一种对上述的上覆压力作用下冻结土体水分迁移模型的检验方法，本试样为高度为11cm，直径为10cm的标准圆柱体。且在试验开始前，对饱和粉土样品施加1mpa的压力并固结24h，以保证土体在不同上覆压力条件下物理力学性质的恒定。试验期间，冷端(顶端)循环液的温度设置为-10℃，热端(底端)循环液的温度设置为 1℃，环境箱温度也设置为 1℃，且温度的波动范围控制在±0.1℃。由于对试样热端施加温度tw为 1℃，对冷端施加温度tc为-10℃，因此沿着x轴可获得冻结锋面温度tf为0℃--0.03℃。与此同时，沿试样高度从下而上在2.5cm、5.0cm、7.5cm和10.0cm四个层位点(测点)均匀布置4个温度、水分集成传感器，其传感器型号为mec10，含水量可选量程为：0～50％和0～100％；分辨率在0～50％内为0.03％，在50～100％内为1％；精度在0～50％内为2％，在50～100％内为2％。温度量程为-40～80℃，分辨率为0.1℃，精度为±0.5℃；以监测冻结过程中温度和未冻水含量的变化。由于是开放系统条件下的冻胀试验，因此对试样底部进行补水，并保证试验过程中补给水头的稳定。试验过程中分别施加50kpa、100kpa、200kpa、400kpa、1000kpa的上覆荷载，以检验不同上覆荷载作用下，土样的冻胀率。具体检验过程如下：以冻结的土体是饱和粉土为例，根据上覆压力不同(0.05mpa、0.1mpa、0.2mpa、0.4mpa、1.0mpa)，设计了开放系统条件下一维冻胀加载试验，对压力作用下冻结土体的水分迁移模型正确性进行验证，首先引入主要参数及其取值如表1所示：表1计算模型主要参数取值参数取值参数取值l/(cm2/s2)334.56×107g/(m/s2)9.81ρw/(g/cm3)1.0ta/(k)273.15β8tf/(℃)0此外，在冻胀试验开始前，测得饱和粉土的分离压力psep为25kpa，饱和粉土的渗透系数kuf＝10-8cm/s，在试验过程中，测得渗透路径lt长度为3cm。其他主要参数见表2所示：表2不同荷载条件下水分迁移速度为便于分析，将表2中的水分迁移速度理论计算值和一维冻胀加载试验测得的试验值进行对比，对比图如图3所示。根据图3和表2可知，水分迁移速度的理论计算值(速度-计算值vff)和一维冻胀加载试验值(速度-试验值vff)的平均相对误差在10％以内，且呈现规律基本一致。虽然上覆压力pob为0.4mpa时，速度-计算值vff为0.9708×10-7cm/s，速度-试验值vff为0.7449×10-7cm/s，相对误差超过10％，但是由于冻胀试验的控制因素多，波动较大，所以对整体变化规律影响不大，影响可忽略。因此，压力作用下冻结土体的水分迁移模型可准确预测不同附加荷载对土体冻胀的抑制情况。本发明实施例提供的技术方案带来的有益效果是：首次从上覆压力控制冰晶分凝的角度，全面合理的阐明了上覆压力对水分迁移的抑制机理和控制因素。研究结果表明，随着上覆压力增大，冰透镜体暖端的表面吸附力几乎保持不变，但冻结缘区的总渗透系数随之呈高阶幂函数形式减小；鉴于水分迁移速度符合达西定律，据此，给出了压力作用下水分迁移速度显式方程。该模型简化了参数变量，根据冻结缘区2内的水分迁移速度可快速预测饱和土体在附加荷载作用下的冻胀率和冻胀量，对完善和发展现有冻胀理论，指导工程实践等方面，如：预测地基、路基变形，防控建筑物、构筑物冻胀破坏等方面具有重要的理论意义和应用价值。以上所述仅为本发明的优选实施例，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。当前第1页1 2 3
技术特征：1.一种上覆压力作用下冻结土体水分迁移模型的构建方法，其特征在于，构建方法包括如下步骤：
步骤1：获取平衡状态下的理论冰压方程；
步骤2：获取冰透镜体暖端的分凝-冻结温度方程；
步骤3：根据上覆压力、步骤1和步骤2中得到的理论冰压方程、分凝-冻结温度方程获取冰透镜体暖端的迁移驱动力方程；
步骤4：获得冰透镜体暖端的表面吸附力方程；
步骤5：获取冻结缘区的总渗透系数方程；
步骤6：获取冰透镜体暖端的水分迁移速度显式方程，完成上覆压力作用下冻结土体水分迁移模型的构建。
2.根据权利要求1所述的上覆压力作用下冻结土体水分迁移模型的构建方法，其特征在于，在所述步骤1中，
根据广义clapeyron方程，可得平衡状态下的理论冰压方程：
式(2)中：ps0为平衡状态下薄膜水受到的理论冰压；l为水的相变潜热，可取3.34×109cm2/s2；vs为冰的比容；ta为纯冰的绝对冻结温度，可取273.15k；t为冰-水界面的温度，t为变量。
3.根据权利要求1所述的上覆压力作用下冻结土体水分迁移模型的构建方法，其特征在于，在所述步骤2中，具体包括：
步骤2.1：假定冰透镜体形成之前实际冰压等于理论冰压；
步骤2.2：当实际冰压等于上覆压力与分离压力之和时，则会产生新的分凝冰，获取分凝冰形成条件方程；
步骤2.3：根据分凝冰形成条件方程获取冰透镜体暖端的分凝-冻结温度方程；
实际冰压方程为：
结合分凝冰形成机制可知，当实际冰压等于上覆压力与分离压力之和时，则会产生新的分凝冰，即：
ps＝pob psep(3)
分凝冰形成条件方程的获取：式(2)和(3)中：ps0为平衡状态下薄膜水受到的理论冰压；ps为实际冰压；pob为上覆压力；psep为分离压力；
结合公式(2)、(3)，可得分凝冰形成条件方程：
变换公式(3)可得分凝-冻结温度方程：
式(4)和(5)中：ts为冰透镜体暖端的分凝-冻结温度。
4.根据权利要求1所述的上覆压力作用下冻结土体水分迁移模型的构建方法，其特征在于，在所述步骤3中，具体包括：
根据程桦提出的薄膜水迁移驱动力模型，可知冰透镜体暖端的迁移驱动力方程等于冰透镜体暖端的理论冰压与上覆压力之差：
式(6)中：pld为迁移驱动力，ts代表t在冰透镜体暖端的分凝-冻结温度。
5.根据权利要求1所述的上覆压力作用下冻结土体水分迁移模型的构建方法，其特征在于，在所述步骤4中，具体包括：
步骤4.1：获取压-吸变量热力学换算系数λ，λ＝-vs/vl＝-1.09，vs为冰的比容，vl是水的比容，λ取-1.09；
步骤4.2：根据迁移驱动力方程和热力学换算系数λ，获取冰透镜体暖端的表面吸附力方程；
式(7)中：plb为表面吸附力，pld为迁移驱动力，λ为压-吸变量的热力学换算系数。
6.根据权利要求1所述的上覆压力作用下冻结土体水分迁移模型的构建方法，其特征在于，在所述步骤5中：
根据托马斯统计的试验结果可知，未冻结区、冻结缘区和已冻结区之间的渗透系数不同，因此设定未冻结区的渗透系数kuf，kuf也为饱和土体在常温条件下的渗透系数，冻结缘区的总渗透系数kff，已冻结区的渗透系数kf，
托马斯推荐的冻结缘渗透系数方程为：
kff＝kuf[1-(ts-tf)]-β(8)
根据所述步骤2得到的分凝-冻结温度方程、公式(8)获取冻结缘区的总渗透系数方程，冻结缘区的总渗透系数方程为：
式(8)和(9)中：β为幂，无量纲；tf为冻结锋面的冻结温度，可取0℃；kuf为饱和土体在常温条件下的渗透系数；kff为冻结缘区的总渗透系数。
7.根据权利要求1所述的上覆压力作用下冻结土体水分迁移模型的构建方法，其特征在于，在所述步骤6中，
根据达西定律，结合冰透镜体暖端的表面吸附力plb、冻结缘区的总渗透系数kff，获取冰透镜体暖端的水分迁移速度显式方程，完成压力作用下冻结土体水分迁移模型的构建，具体包括：
假设冻结缘区内的水分迁移速度符合达西定律，结合公式(7)、(9)可得水分迁移速度显式方程：
式(10)中，lt为水分迁移的渗透路径、ρw为水的密度、g为重力加速度、vff为冰透镜体暖端的水分迁移速度。
8.根据权利要求1所述的上覆压力作用下冻结土体水分迁移模型的构建方法，其特征在于，利用冻结土体的试样在上覆压力作用下进行水分迁移模型的构建，
所述冻结土体的试样包括未冻结区、冻结缘区和已冻结区，所述未冻结区的一端与所述冻结缘区的一端连接，所述冻结缘区的另一端与所述已冻结区的一端连接，所述未冻结区、所述冻结缘区和所述已冻结区依次排列设置，在所述已冻结区内并靠近所述冻结缘区的一侧形成冰透镜体，在所述未冻结区内靠近所述冻结缘区的一侧形成冻结锋面。
9.根据权利要求8所述的上覆压力作用下冻结土体水分迁移模型的构建方法，其特征在于，
所述未冻结区的另一端为热端，所述已冻结区的另一端为冷端，所述冰透镜体靠近冻结缘区的一端为所述冰透镜体的暖端。
10.根据权利要求9所述的上覆压力作用下冻结土体水分迁移模型的构建方法，其特征在于，
所述冻结缘区内设置有若干土颗粒，每个所述土颗粒外围均包裹有一层未冻水膜，水分可通过所述未冻水膜从所述未冻结区向所述冰透镜体的暖端迁移。
技术总结本发明提供一种上覆压力作用下冻结土体水分迁移模型的构建方法，构建方法包括如下步骤：步骤1：获取平衡状态下的理论冰压方程；步骤2：获取冰透镜体暖端的分凝‑冻结温度方程；步骤3：根据上覆压力、步骤1和步骤2中得到的理论冰压方程、分凝‑冻结温度方程获取冰透镜体暖端的迁移驱动力方程；步骤4：获得冰透镜体暖端的表面吸附力方程；步骤5：获取冻结缘区的总渗透系数方程；步骤6：获取冰透镜体暖端的水分迁移速度显式方程，完成上覆压力作用下冻结土体水分迁移模型的构建。根据冻结缘区内的水分迁移速度可快速预测饱和土体在附加荷载作用下的冻胀率和冻胀量，该发明完善和发展了现有冻胀理论，可指导工程实践。
技术研发人员：马茂艳;陈汉青;曹广勇;林键;彭世龙;彭曙光
受保护的技术使用者：安徽建筑大学
技术研发日：2020.11.09
技术公布日：2021.03.12
