寒区高等级公路路堑边坡春季浅层滑塌机理研究
|
摘要
我国绝大多数寒冷地区属季节性冻土区,春季路堑边坡浅层滑塌失稳是一种多见的工程病害,一直受到国内外寒冷地区公路工程、岩土工程界的关注。入冬,大地封冻由表及里,土中水分向地表迁移;春季融化亦由表及里,表层融水不能渗入下部未融化的土层中,使上部土体含水量增加、强度减小,使原来的均质边坡成为非均质的边坡;有时冻融就导致浅层滑塌,显然地下水的活动起了关键作用。
本文紧紧抓住冻融过程中地下水分的迁移这个核心问题,以工程应用为目的,土样试验和大型模型试验相结合,深入研究,建立冻融过程土中水分迁移模型,与数值分析结合,探讨路堑边坡冻融失稳的机理。
研究先从各方面条件均能够很好控制的土样试验入手,针对路堑边坡含水量低、地下水埋藏深的特点,采取封闭系统,底部设隔水界面,单向冻结。在土样侧面布置温度传感器,通过顶板的温度控制冻结的速率,实时观测温度沿深度的变化,控制冻结、融化的进程,测试不同深度土中的含水量。原状土试验模拟天然状态,重塑土试验便于区分冻结、融化两个过程的作用,研究土中水分在冻融过程迁移的规律。
试验中同步记录的各深度处的温度,测定冻、融过程中含水量随深度的变化的数据,与初始含水量比较,为进一步建立冻融过程中水分迁移模型奠定了基础。
试验中,选取相同含水量的土样,设定了5种温度模式,测定不同温度梯度下土中水分的迁移量;按五种设定的含水量调制土样,在相同温度梯度条件下,对不同初始含水量的土样,分别测定土中水分的迁移量。对含水量均匀的土样,分别测试完全冻结和冻结后从顶面融化到设定位置的含水量分布,研究融化到设定位置后,冻结面上水分聚集层的位置及其含水量的变化情况。
试验结果表明,冻结并融化到设定位置后,土样未融化部分出现横向和竖向的裂缝、明显的冻融分界线;控制不同顶板温度,不同含水量的试样中各点处温度的变化趋势基本相同;不同初始含水量,冻融后含水量的分布曲线的趋势是相同的,上部土层的含水量增加、下部土层的含水量减少;顶板控温相同,初始含水量大的土样,冻结要消耗的潜热量大,冻结速度慢,迁移的水分多,水分集聚亦多;不同冻结温度和融化温度对试样中水分的迁移影响非常明显,顶板冻结温度越高,水分积聚越明显。重塑的粉质粘土含水量均匀的土样,自上而下冻结后出现含冰量不均匀现象,有的形成冰夹层;冻结后自上而下融化到设定位置后,在冻融交界面之上水分聚集,出现一个含水量最大的薄层。
针对现有的水分迁移机制仍然存在的疑问,冻融过程中水分迁移模型偏于复杂,有的参数值不能直接通过试验获得,对土体工程特性考虑不足,一些建模尚缺乏实验论证,应用的结果往往不够理想的局面,建模不追求对多孔多相介质带相变的固、液、汽、热耦合的统一研究,暂不考虑水分迁移对温度场的反馈影响,强调水气扩散在冻结过程中迁移的重要性,从工程角度出发,宏观上把握水分迁移量与温度梯度的关系,分别考虑冻结、融化过程中水分迁移的机理,借助Fick和Darcy定律进一步改进建立模型。
(1)假定非饱和土中水分在冻结过程中主要以水汽扩散形式的迁移,在此过程中液态水通过转化为水汽参加迁移,未具体考虑转化过程。推导、计算中,设定了冻结锋面以下土层中水汽的浓度方程,建立了浓度-深度随冻结锋面发展的移动坐标系。为考虑温度梯度的影响,在Fick定律公式的右边乘上温度梯度,得到了修正的非饱和土的Fick定律,作为冻结过程的水分迁移模型,发展了计算方法和程序。
本文采用一个综合的扩散项,而不是另增加一个由温度梯度引起的扩散项。从物理上看,表达了温度梯度直接影响扩散率,而不是在原来的水气扩散量的基础上再增加一部分;从数学上看,少了一个参数,便于参数的估计。
(2)假定融化过程中,融土中水分迁移主要是在重力作用下的竖向渗流。用增加温度梯度引起的水分迁移项来修正非饱和土达西定律,作为融化过程中的水分迁移模型。
接着,以试验得到的冻结、融化后含水量沿试样深度的分布曲线为目标函数,以本文所建立的冻、融中水分迁移模型为正演模块,借助遗传算法反演了上述两个模型的参数。分别把反演结果及其平均值代入模型中计算含水量沿土样深度的分布曲线与试验得到的各个曲线对比,论证了模型可靠性及本文思路的可行性。
本文将较大比例尺的室内边坡冻融失稳模型试验设计为土样试验和现场原位观测之间的一个既有一定体积、又可以很好控制各种边界条件的中间环节,弥补土样尺寸小受边界影响大的缺陷。在归纳国内外类似试验的方案、确定模型相似律的基础上,设计、制作了60厘米高的边坡模型以及温度控制条件,在表层密集布设温度、含水量观测传感器,每隔5分钟采集一组数据。试验采用了根据电磁脉冲反射方法研制的时域反射仪TDR,对量值随温度的变化做了标定。经过35天,历时840小时,5个冻融循环,试验显示了水分迁移、冻裂、分凝冰层生成、坡面局部滑塌等现象。通过记录的数据分析,说明温度场控制得比较好,含水量变化比较复杂。对第一个冻融循环的仔细分析,说明了冻结中水分扩散向上迁移,融化过程中向下渗流的规律。依据试验结果阐述了坡面滑塌的主要原因是土层含水量提高引起抗剪强度的降低。
为了研究冻融对土的抗剪强度的影响,定量地揭示含水量提高所引起土层抗剪强度降低的程度,通过制备不同含水量的土样,经过冻融过程,采用三轴不排水剪切试验,测得土的抗剪强度参数。结果表明抗剪强度指标确实下降很大,验证了本文试验中边坡发生滑塌主要是由于浅层土体含水量增加很多,导致土的抗剪强度大幅下降是重要因素的结论。
最后,本文对试验的边坡模型,进行了冻融稳定性数值分析。在有限元分析软件ABAQUS的基础上二次开发,建立了冻融边坡模型的有限元模型,设置了边界条件和单元参数。在垂直坡面3.0米的深度内,按0.3米划分网格为冻融网格,以下为非冻融网格,借助ABAQUS/CAE生成inp形式输入文件。根据模型试验中观测到的温度场,冻结和融化过程中水分的迁移,含水量和三轴试验结果统计关系修正冻融单元的抗剪强度参数,计算了冻结后融化到5个不同深度时模型的变形和内力,结果显现了浅层滑动的可能性。
Most of the cold regions in China are seasonal frozen regions where shallow slide of highway cutting slope is a popular phenomenon in spring. Attentions have been paid to the issue in the areas of road engineering and geotechnical engineering in cold region for many years. In winter, the soil moisture migrates upwards as the ground frozen from the surface to the deep. In spring, the ground thawing starts from the surface to the deep, the water at the shallow layer permeates downwards by gravity and cannot infiltrate into the frozen soil beneath, and then it accumulates in a thin soil layer just above the interface. This results in the increase of the water content and the decrease of the shear strength of the soil layer, and making the original homogeneous slope become inhomogeneous, sometimes which makes the shallow slide. Obviously, the water migration in soil during freeze and thaw process plays a key role.
The soil moisture migration in the freeze-thaw process is taken as the core point in this dissertation; Combining with soil sample and large scale model experiments, two mathematical models are built up separately for the migrations in freeze and thaw processes for engineering application purpose. The instability mechanism of highway cutting slope during freeze-thaw process is studied with also numerical analysis.
The study begins with soil sample experiment, which conditions are well controlled. A close system with watertight boundary at the bottom of the soil sample is adopted, the samples are frozen and thawed in one-way to the practical situation of low water content and deep ground water table in all cutting slope. The freezing rate is handled through the temperature of the top plate, temperature variation with depth is observed regularly time to time from the thermal sensors set at a side of the soil sample, and then the frozen or thaw depth is controlled. Undisturbed soil test is to simulate the natural state, while reshaped soil test is to distinguish the effects of freezing and thawing, the both are exploring the soil moisture migration rules in the two processes.
Temperatures at different depths were measured and recorded simultaneously during the test, and the distribution of water content with depth were measured after the sample frozen or thawed at a given depth, and compared with the original water content. All the data acquired in the tests are building a firm base for the mathematic models.
In the experiments, 5 temperature conditions for samples with the same initial water content were adopted to measure the moisture migration amount under different temperature gradients; while 5 reshaped samples with different given initial water contents were adopted under the same temperature condition also to measure the water migration amounts. The water content distributions were measured when the samples completely frozen or thawed from the top to a given depth afterwards from the initial uniform distribution. The position of the soil layer with the richest water content above the freeze-thaw interface and the variation of the water content at that depth were studied for the thawed cases.
The test results show that some transverse and longitudinal cracks appear in unthawed portion after the sample completely frozen and then thawed to the given depth, a distinguish freeze-thaw interface is also observed; the temperatures at various depths show the similar trend with different top temperatures even with different water contents; the water content curves after freeze-thaw are the similar on the whole, i.e. the water content increases in the upper soil layer and decreases in the lower layer; the more initial water content, the more potential heat consumption in freezing, the slower freezing rate, the more moisture migrated and the more water aggregated above the interface; the effect of temperature condition on moisture migration is quite obvious, the higher freezing temperature at top plate, the more water aggregates above the interface; the inhomogeneous ice content, even segregated ice, appear in reshaped silt-clay samples with initial homogeneous water content after frozen from top to bottom; a thin soil layer with the richest water content appears just above the interface after the sample frozen completely from top and then thawed to the given depth.
There are still queries now to the existing water migration mechanism theories: the moisture migration models are too complicated for engineering practical purpose, values of some parameters in the models could not be obtained from experiment directly, engineering characteristics of soil mass is not taken into account enough, almost of the models are lack of validation, application results of these models are not satisfactory in some cases. From this state of the art, the modeling is not in pursuit of the comprehensive description of real solid, liquid and gas phase coupled with heat in porous media, and the feedback effect of moisture migration on temperature field is not taken into account temporarily in this paper. The importance of moisture diffusion in water migration during freezing is emphasized, and the relationship between water migration amount and temperature gradient is handled in addition, from the engineering point of view. The water migration mechanisms for two processes of freezing and thawing are considered separately, and two models were worked out by modification of the Fick law and the Darcy law.
(1) Assuming the moisture in unsaturated soil migrates in the way of diffusion during freeze process. In this process, liquid water takes part in the migration by transformed into moisture, and the detailed transformation is not considered. The concentration equation of the moisture in the soil beneath the frozen front is presented, and a moving coordinate system is established for the concentration-depth relationship with the freezing front development. In order to take the effect of temperature gradient into account, temperature gradient is multiplied to the right side of the Fick law, thus an improved Fick law for unsaturated soil is obtained as the water migration model of freeze process. A set of calculation procedure and program are developed at the same time.
A comprehensive diffusion term is adopted in this paper, instead of adding a diffusion term caused by temperature gradient. From physics point of view, it describes the fact that temperature gradient influent the diffusion rate directly, not adding an additional amount on the original moisture diffusion. From mathematics point of view, eliminating one parameter is convenient to parameters estimation.
(2) Assuming the water migration in thawed soil is mainly the vertical permeation from the gravity. An additional term of water permeation caused by the temperature gradient is added to the Darcy law for unsaturated soil, as a water migration model for the thaw process.
Furthermore, values of parameters in the two models are inversed by means of the Generic Algorithm, where the objective function is taken as the curve of water content– depth after the frozen and thawed processes respectively obtained from the experiments, and the forward calculation modules is built from the programs for the two models. The water content - depth curves are then calculated from the inversed parameter values and their mean respectively. The resulted curves are compared with those from experiments, thus the reliability of the model and the feasibility of the approach of this paper are validated.
A large scale slope model experiment for instability after frozen-thaw is designed as a link between the soil sample experiment and the in situ experiment in this dissertation, with well-controlled boundary conditions and certain volume to eliminate the effect of boundaries in soil sample test. The designs of similar model tests both at home and abroad are summarized, the similarity law for model is determined. A slope model with 60 centimeters high is designed and made in laboratory, and the temperature condition is preset. Thermal sensors and sensors for water content are buried densely in the upper soil layers of the model. The data are acquired every 5 minutes during the experiment. Time domain reflector (TDR) is adopted in the experiment, which is developed from electromagnetic pulse reflection method, and the measured values are scaled to eliminate the shifting of readings with temperature. During 35 days, 840 hours which consists of 5 freeze-thaw cycles, the experiment shows the phenomenon such as water migration, frozen crack, segregated ice, and local slides on the slope, etc. The analysis on the recorded data shows that the temperature field is well-controlled and the water content variation is complex. The detailed data analysis of the first freeze-thaw cycle indicates rules that the moisture moves upwards by diffusion during the freezing, and permeates downwards by gravity during the thawing. From the results, it can be concluded that the main reason of shallow slide on slope is the significant decrease of shear strength in the thin soil layer just above the freeze-thaw interface resulted from the large increase of water content.
In order to deal with the effect of freeze– thaw on shear strength of soil, to reveal the degree of decrease of soil shear strength from the increase of the water content quantitatively, soil samples with different water contents are reshaped, and are further frozen and then thawed, the shear strength parameters of the soil are tested by means of triaxial undrained shearing experiment. The result shows that the values of shear strength parameters decease certainly to a large extent, thus the above conclusion is validated that the main reason of shallow slides on the model slope is the large decrease of the shear strength of the shallow soil resulted from the increase of water content there.
Finally, numerical analysis of the slope model for stability in the freeze-thaw experiment is carried out. A second development is made by means of ABAQUS software for finite element analysis, a finite element model is constructed for the model in the experiment, and the boundary conditions and element parameters are assigned. In the depth of 3.0 meters vertical from the surface, grids with 0.3 meters distance are freeze-thaw grids whose parameters vary with water content during the test, and others are not. The input data files with“inp”type are generated by ABAQUS/CAE. The shear strengths of freeze– thaw grids are modified from the statistical relationship with water content from the triaxial test mentioned above. The water content data measured in the experiment, the deformations and stresses of elements after frozen and then thawed to 5 depths are calculated respectively. The results show the potential of shallow slide on the slope.
本文紧紧抓住冻融过程中地下水分的迁移这个核心问题,以工程应用为目的,土样试验和大型模型试验相结合,深入研究,建立冻融过程土中水分迁移模型,与数值分析结合,探讨路堑边坡冻融失稳的机理。
研究先从各方面条件均能够很好控制的土样试验入手,针对路堑边坡含水量低、地下水埋藏深的特点,采取封闭系统,底部设隔水界面,单向冻结。在土样侧面布置温度传感器,通过顶板的温度控制冻结的速率,实时观测温度沿深度的变化,控制冻结、融化的进程,测试不同深度土中的含水量。原状土试验模拟天然状态,重塑土试验便于区分冻结、融化两个过程的作用,研究土中水分在冻融过程迁移的规律。
试验中同步记录的各深度处的温度,测定冻、融过程中含水量随深度的变化的数据,与初始含水量比较,为进一步建立冻融过程中水分迁移模型奠定了基础。
试验中,选取相同含水量的土样,设定了5种温度模式,测定不同温度梯度下土中水分的迁移量;按五种设定的含水量调制土样,在相同温度梯度条件下,对不同初始含水量的土样,分别测定土中水分的迁移量。对含水量均匀的土样,分别测试完全冻结和冻结后从顶面融化到设定位置的含水量分布,研究融化到设定位置后,冻结面上水分聚集层的位置及其含水量的变化情况。
试验结果表明,冻结并融化到设定位置后,土样未融化部分出现横向和竖向的裂缝、明显的冻融分界线;控制不同顶板温度,不同含水量的试样中各点处温度的变化趋势基本相同;不同初始含水量,冻融后含水量的分布曲线的趋势是相同的,上部土层的含水量增加、下部土层的含水量减少;顶板控温相同,初始含水量大的土样,冻结要消耗的潜热量大,冻结速度慢,迁移的水分多,水分集聚亦多;不同冻结温度和融化温度对试样中水分的迁移影响非常明显,顶板冻结温度越高,水分积聚越明显。重塑的粉质粘土含水量均匀的土样,自上而下冻结后出现含冰量不均匀现象,有的形成冰夹层;冻结后自上而下融化到设定位置后,在冻融交界面之上水分聚集,出现一个含水量最大的薄层。
针对现有的水分迁移机制仍然存在的疑问,冻融过程中水分迁移模型偏于复杂,有的参数值不能直接通过试验获得,对土体工程特性考虑不足,一些建模尚缺乏实验论证,应用的结果往往不够理想的局面,建模不追求对多孔多相介质带相变的固、液、汽、热耦合的统一研究,暂不考虑水分迁移对温度场的反馈影响,强调水气扩散在冻结过程中迁移的重要性,从工程角度出发,宏观上把握水分迁移量与温度梯度的关系,分别考虑冻结、融化过程中水分迁移的机理,借助Fick和Darcy定律进一步改进建立模型。
(1)假定非饱和土中水分在冻结过程中主要以水汽扩散形式的迁移,在此过程中液态水通过转化为水汽参加迁移,未具体考虑转化过程。推导、计算中,设定了冻结锋面以下土层中水汽的浓度方程,建立了浓度-深度随冻结锋面发展的移动坐标系。为考虑温度梯度的影响,在Fick定律公式的右边乘上温度梯度,得到了修正的非饱和土的Fick定律,作为冻结过程的水分迁移模型,发展了计算方法和程序。
本文采用一个综合的扩散项,而不是另增加一个由温度梯度引起的扩散项。从物理上看,表达了温度梯度直接影响扩散率,而不是在原来的水气扩散量的基础上再增加一部分;从数学上看,少了一个参数,便于参数的估计。
(2)假定融化过程中,融土中水分迁移主要是在重力作用下的竖向渗流。用增加温度梯度引起的水分迁移项来修正非饱和土达西定律,作为融化过程中的水分迁移模型。
接着,以试验得到的冻结、融化后含水量沿试样深度的分布曲线为目标函数,以本文所建立的冻、融中水分迁移模型为正演模块,借助遗传算法反演了上述两个模型的参数。分别把反演结果及其平均值代入模型中计算含水量沿土样深度的分布曲线与试验得到的各个曲线对比,论证了模型可靠性及本文思路的可行性。
本文将较大比例尺的室内边坡冻融失稳模型试验设计为土样试验和现场原位观测之间的一个既有一定体积、又可以很好控制各种边界条件的中间环节,弥补土样尺寸小受边界影响大的缺陷。在归纳国内外类似试验的方案、确定模型相似律的基础上,设计、制作了60厘米高的边坡模型以及温度控制条件,在表层密集布设温度、含水量观测传感器,每隔5分钟采集一组数据。试验采用了根据电磁脉冲反射方法研制的时域反射仪TDR,对量值随温度的变化做了标定。经过35天,历时840小时,5个冻融循环,试验显示了水分迁移、冻裂、分凝冰层生成、坡面局部滑塌等现象。通过记录的数据分析,说明温度场控制得比较好,含水量变化比较复杂。对第一个冻融循环的仔细分析,说明了冻结中水分扩散向上迁移,融化过程中向下渗流的规律。依据试验结果阐述了坡面滑塌的主要原因是土层含水量提高引起抗剪强度的降低。
为了研究冻融对土的抗剪强度的影响,定量地揭示含水量提高所引起土层抗剪强度降低的程度,通过制备不同含水量的土样,经过冻融过程,采用三轴不排水剪切试验,测得土的抗剪强度参数。结果表明抗剪强度指标确实下降很大,验证了本文试验中边坡发生滑塌主要是由于浅层土体含水量增加很多,导致土的抗剪强度大幅下降是重要因素的结论。
最后,本文对试验的边坡模型,进行了冻融稳定性数值分析。在有限元分析软件ABAQUS的基础上二次开发,建立了冻融边坡模型的有限元模型,设置了边界条件和单元参数。在垂直坡面3.0米的深度内,按0.3米划分网格为冻融网格,以下为非冻融网格,借助ABAQUS/CAE生成inp形式输入文件。根据模型试验中观测到的温度场,冻结和融化过程中水分的迁移,含水量和三轴试验结果统计关系修正冻融单元的抗剪强度参数,计算了冻结后融化到5个不同深度时模型的变形和内力,结果显现了浅层滑动的可能性。
Most of the cold regions in China are seasonal frozen regions where shallow slide of highway cutting slope is a popular phenomenon in spring. Attentions have been paid to the issue in the areas of road engineering and geotechnical engineering in cold region for many years. In winter, the soil moisture migrates upwards as the ground frozen from the surface to the deep. In spring, the ground thawing starts from the surface to the deep, the water at the shallow layer permeates downwards by gravity and cannot infiltrate into the frozen soil beneath, and then it accumulates in a thin soil layer just above the interface. This results in the increase of the water content and the decrease of the shear strength of the soil layer, and making the original homogeneous slope become inhomogeneous, sometimes which makes the shallow slide. Obviously, the water migration in soil during freeze and thaw process plays a key role.
The soil moisture migration in the freeze-thaw process is taken as the core point in this dissertation; Combining with soil sample and large scale model experiments, two mathematical models are built up separately for the migrations in freeze and thaw processes for engineering application purpose. The instability mechanism of highway cutting slope during freeze-thaw process is studied with also numerical analysis.
The study begins with soil sample experiment, which conditions are well controlled. A close system with watertight boundary at the bottom of the soil sample is adopted, the samples are frozen and thawed in one-way to the practical situation of low water content and deep ground water table in all cutting slope. The freezing rate is handled through the temperature of the top plate, temperature variation with depth is observed regularly time to time from the thermal sensors set at a side of the soil sample, and then the frozen or thaw depth is controlled. Undisturbed soil test is to simulate the natural state, while reshaped soil test is to distinguish the effects of freezing and thawing, the both are exploring the soil moisture migration rules in the two processes.
Temperatures at different depths were measured and recorded simultaneously during the test, and the distribution of water content with depth were measured after the sample frozen or thawed at a given depth, and compared with the original water content. All the data acquired in the tests are building a firm base for the mathematic models.
In the experiments, 5 temperature conditions for samples with the same initial water content were adopted to measure the moisture migration amount under different temperature gradients; while 5 reshaped samples with different given initial water contents were adopted under the same temperature condition also to measure the water migration amounts. The water content distributions were measured when the samples completely frozen or thawed from the top to a given depth afterwards from the initial uniform distribution. The position of the soil layer with the richest water content above the freeze-thaw interface and the variation of the water content at that depth were studied for the thawed cases.
The test results show that some transverse and longitudinal cracks appear in unthawed portion after the sample completely frozen and then thawed to the given depth, a distinguish freeze-thaw interface is also observed; the temperatures at various depths show the similar trend with different top temperatures even with different water contents; the water content curves after freeze-thaw are the similar on the whole, i.e. the water content increases in the upper soil layer and decreases in the lower layer; the more initial water content, the more potential heat consumption in freezing, the slower freezing rate, the more moisture migrated and the more water aggregated above the interface; the effect of temperature condition on moisture migration is quite obvious, the higher freezing temperature at top plate, the more water aggregates above the interface; the inhomogeneous ice content, even segregated ice, appear in reshaped silt-clay samples with initial homogeneous water content after frozen from top to bottom; a thin soil layer with the richest water content appears just above the interface after the sample frozen completely from top and then thawed to the given depth.
There are still queries now to the existing water migration mechanism theories: the moisture migration models are too complicated for engineering practical purpose, values of some parameters in the models could not be obtained from experiment directly, engineering characteristics of soil mass is not taken into account enough, almost of the models are lack of validation, application results of these models are not satisfactory in some cases. From this state of the art, the modeling is not in pursuit of the comprehensive description of real solid, liquid and gas phase coupled with heat in porous media, and the feedback effect of moisture migration on temperature field is not taken into account temporarily in this paper. The importance of moisture diffusion in water migration during freezing is emphasized, and the relationship between water migration amount and temperature gradient is handled in addition, from the engineering point of view. The water migration mechanisms for two processes of freezing and thawing are considered separately, and two models were worked out by modification of the Fick law and the Darcy law.
(1) Assuming the moisture in unsaturated soil migrates in the way of diffusion during freeze process. In this process, liquid water takes part in the migration by transformed into moisture, and the detailed transformation is not considered. The concentration equation of the moisture in the soil beneath the frozen front is presented, and a moving coordinate system is established for the concentration-depth relationship with the freezing front development. In order to take the effect of temperature gradient into account, temperature gradient is multiplied to the right side of the Fick law, thus an improved Fick law for unsaturated soil is obtained as the water migration model of freeze process. A set of calculation procedure and program are developed at the same time.
A comprehensive diffusion term is adopted in this paper, instead of adding a diffusion term caused by temperature gradient. From physics point of view, it describes the fact that temperature gradient influent the diffusion rate directly, not adding an additional amount on the original moisture diffusion. From mathematics point of view, eliminating one parameter is convenient to parameters estimation.
(2) Assuming the water migration in thawed soil is mainly the vertical permeation from the gravity. An additional term of water permeation caused by the temperature gradient is added to the Darcy law for unsaturated soil, as a water migration model for the thaw process.
Furthermore, values of parameters in the two models are inversed by means of the Generic Algorithm, where the objective function is taken as the curve of water content– depth after the frozen and thawed processes respectively obtained from the experiments, and the forward calculation modules is built from the programs for the two models. The water content - depth curves are then calculated from the inversed parameter values and their mean respectively. The resulted curves are compared with those from experiments, thus the reliability of the model and the feasibility of the approach of this paper are validated.
A large scale slope model experiment for instability after frozen-thaw is designed as a link between the soil sample experiment and the in situ experiment in this dissertation, with well-controlled boundary conditions and certain volume to eliminate the effect of boundaries in soil sample test. The designs of similar model tests both at home and abroad are summarized, the similarity law for model is determined. A slope model with 60 centimeters high is designed and made in laboratory, and the temperature condition is preset. Thermal sensors and sensors for water content are buried densely in the upper soil layers of the model. The data are acquired every 5 minutes during the experiment. Time domain reflector (TDR) is adopted in the experiment, which is developed from electromagnetic pulse reflection method, and the measured values are scaled to eliminate the shifting of readings with temperature. During 35 days, 840 hours which consists of 5 freeze-thaw cycles, the experiment shows the phenomenon such as water migration, frozen crack, segregated ice, and local slides on the slope, etc. The analysis on the recorded data shows that the temperature field is well-controlled and the water content variation is complex. The detailed data analysis of the first freeze-thaw cycle indicates rules that the moisture moves upwards by diffusion during the freezing, and permeates downwards by gravity during the thawing. From the results, it can be concluded that the main reason of shallow slide on slope is the significant decrease of shear strength in the thin soil layer just above the freeze-thaw interface resulted from the large increase of water content.
In order to deal with the effect of freeze– thaw on shear strength of soil, to reveal the degree of decrease of soil shear strength from the increase of the water content quantitatively, soil samples with different water contents are reshaped, and are further frozen and then thawed, the shear strength parameters of the soil are tested by means of triaxial undrained shearing experiment. The result shows that the values of shear strength parameters decease certainly to a large extent, thus the above conclusion is validated that the main reason of shallow slides on the model slope is the large decrease of the shear strength of the shallow soil resulted from the increase of water content there.
Finally, numerical analysis of the slope model for stability in the freeze-thaw experiment is carried out. A second development is made by means of ABAQUS software for finite element analysis, a finite element model is constructed for the model in the experiment, and the boundary conditions and element parameters are assigned. In the depth of 3.0 meters vertical from the surface, grids with 0.3 meters distance are freeze-thaw grids whose parameters vary with water content during the test, and others are not. The input data files with“inp”type are generated by ABAQUS/CAE. The shear strengths of freeze– thaw grids are modified from the statistical relationship with water content from the triaxial test mentioned above. The water content data measured in the experiment, the deformations and stresses of elements after frozen and then thawed to 5 depths are calculated respectively. The results show the potential of shallow slide on the slope.
引文
1张喜发,辛德刚,张冬青等.季节冻土区高速公路路基土中的水分迁移变化.冰川冻土, 2004, 26(4): 454-460
2吴玮江.季节性冻融作用与斜坡整体变形破坏.中国地质灾害与防治学报, 1996, 7(4): 59-64
3许珊珊,高伟.寒区公路路堑边坡抗冻融稳定性分析.低温建筑技术, 2006(2): 10-12
4 L. P. Joseph, Simulation and Field Validation of Heat and Water Flow during Soil Freezing and Thawing. Thesis for Doctoral Degree of Oregon State University, 1987
5王彩霞等.冰冻地区路堑边坡溜方原因分析与稳定性验算方法.东北公路, 2001, 24(1)
6中华人民共和国国家标准,冻土工程地质勘察规范GB50324,中国计划出版社, 2001
7杨润田,林凤桐.多年冻土区水文地质及工程地质学.哈尔滨东北林业大学出版社, 1986
8曹宏章.饱和颗粒土冻结过程中的多场耦合研究.中国科学院研究生院博士学位论文, 2006
9王晓春,张倬元.寒区工程与冻融力学.地学前缘, 2000, 7(增刊): 99-140
10程国栋.中国冻土研究近今进展.地理学报, 1990, 45(2): 220-223
11Г.П.马祖罗夫,冻土物理力学性质.北京煤炭工业出版社, 1980
12Э.Д.叶尔绍夫,工程冻土学,莫斯科大学出版社出版, 1999
13车茜.基于BP神经网络的季节冻土水分迁移机理研究.吉林大学硕士学位论文, 2008
14陈玉超.冻融环境下岩土边坡稳定性研究初探.西安科技大学土木工程学院, 2006: 3-20
15 R. L. Lytton, D. E. Pufahl, C. H. Michalak, H. S. Liang, and B. J. Dempsey, . An Integrated Model of The Climatic Effects on Pavements. Federal Highway Administration Report No. FHWA-RD-1990, 90(33)
16童长江,管枫年.土的冻胀与建筑冻害防治.水利电力出版社, 1985
17 H. A.崔托维奇,张长庆,朱元林译.冻土力学.北京科学出版社, 1985
18吴紫汪,马巍著.冻土强度与蠕变.兰州大学出版社, 1993
19周幼吾等.中国冻土,科学出版社, 2000
20周幼吾等.近40年来东北地区季节冻结与融化层温度变化与气候变暖.第
五届全国冰川冻土学大会论文集(上),甘肃文化出版, 1996: 3-10
21牛富俊,张鲁新,俞祁浩等.青藏高原多年冻土区斜坡类型及典型斜坡稳定性研究.冰川冻土, 2002, 24(5): 608-613
22丁德文,罗学波.冻土热工模型试验的理论基础.科学通报, 1979(8): 360-364
23 B. Y. Liu, M. A. Nearing, L. M. Risse. Slope Gradient Effects on Soil Loss for Slopes . Transactions of the ASAE, 1994, 37(6)
24 B. Y. Liu, M. A. Nearing, P. J. Shi, et al. Slope Length Effects on Soil Loss for Steep Slopes. Soil Sci. Soc. of Am. J., 2000, 64(5)
25铁道部建设司,青藏铁路多年冻土区工程勘察暂行规定, 2001
26铁道第一勘测设计院,西大滩至安多多年冻土区综合地质报告,中国铁道出版社, 2001
27铁道部建设司,青藏铁路多年冻土区工程设计暂行规定, 2002
28铁道第三勘察设计院,冻土工程,中国铁道出版社, 2002
29黄小铭等.青藏高原多年冻土区土冻结融化特征探讨.青藏铁路多年冻土科研成果论文集,中铁西北科学研究院, 2003
30中铁西北科学研究院,青藏铁路风火山多年冻土长期综合观测与试验研究,中国铁道出版社, 2006
31童长江等.青藏公路沿线冻土环境工程地质评价及冻土工程处理.第五届全国冰川冻土学大会论文集(上),甘肃文化出版社, 1996: 365-372
32王国尚等.青藏公路沿线冻土区自然环境及工程环境的变化.第五届全国冰川冻土学大会论文集(上),甘肃文化出版社, 1996: 43-50
33喻文学,武敬民.青藏公路多年冻土区沥青路面路基高度问题.西安公路学院学报,高原多年冻土区道路工程研究报告集,西安公路学院学报编辑部出版, 1983: 25-47
34章金钊.青藏公路沥青路面冻土路基设计原则探讨.第五届全国冰川冻土学大会论文集(上),甘肃文化出版社, 1996: 383-388
35李杨.季节冻土水分迁移模型研究.吉林大学博士学位论文, 2008
36赵安平.季冻区路基土冻胀的微观机理研究.吉林大学博士学位论文, 2008
37李杰.边坡稳定性问题的理论分析及其应用研究.大连理工大学硕士学位论文, 2002
38李雨浓.季冻区高速公路路基冻害影响因素分析.吉林大学硕士学位论文, 2007
39刘鸿绪.对冻结过程力学中若干问题的商榷.黑龙江省低温建筑科学研究所, 1988
40刘鸿绪.建筑基础的冻胀力.第五届全国冰川冻土学大会论文集(下),甘肃文化出版社, 1996
41武鹤,高伟,王国峰等.寒区公路土质路堑边坡滑塌原因及其防治.自然灾害学报, 2006, 15(3): 66-70
42宋高嵩,陈阳东,温福春等.季冻区路基边坡的水稳定性影响分析与防治.公路交通科技, 2005 , 22(6): 65-67
43 A. D. Abrahams, et. al., Hydraulics on a Semiarid Hill Slope. Southern Arizon Earth Surface Processes and Land Forms(21), 1996
44戴惠民,王兴隆.季节冻土区公路桥涵地基土冻胀性的研究.中国公路学报, 1989, 2(4): 18-25
45童长江等,青藏公路沿线冻土环境工程地质评价及冻土工程处理.第五届全国冰川冻土学大会论文集(下),甘肃文化出版社, 1996: 65-72
46雷汉忠.冻融土的工程性质及利用.煤矿设计, 1994(6): 39-44
47陈肖柏,王雅卿,刘建坤等.土的冻结作用与地基.北京科学出版社, 2006: 3-60
48杨平.原状土与冻融土物理力学性能差异性研究.南京林业大学学报, 2001, 25(2): 68-70
49杨平等.人工冻融土物理力学性能研究.冰川冻土, 2002, 24(5): 665-667
50武鹤等.寒区路堑人工土质边坡滑塌原因与稳定技术研究.黑龙江工程学院学报, 2005, 19(2)
51朱选举.冻土中未冻水迁移规律研究.清华大学硕士学位论文, 2000
52郭占荣,荆恩春,聂振龙.冻结期和冻融期土壤水分运移特征分析.水科学进展, 2002, 13(3): 298-303
53梁若筠.基于神经网络方法的冻土水分迁移研究.兰州大学硕士学位论文, 2007
54 R. D. Miller. Soil Freezing in Relation to Pore Water Pressure and Temperature in Permafrost. The North American Contribution to the 2nd International Conference on Permafrost, Yakutsk. National Academy of Sciences, Washington, D. C., 1973
55徐敩祖.国外对冻土中水分迁移课题的研究.冰川冻土, 1982, 4(3): 97-103
56徐敩祖,邓友生.饱水正冻土中水分迁移特征研究.第三届全国冻土学术会议论文选集,科学出版社, 1989: 204-210
57邓友生,徐学祖.非饱和土壤水扩散系数的变化特征.冰川冻土, 1994, 16 (3): 251-258
58徐敩祖,王家澄,张立新.冻土物理学.北京科学出版社, 2001
59 R. Ulusay, C. G?k?eo?lu, H. S?nmez, et al. Causes Mechanism and Environmental Impacts of Instabilities at Himmetoglu Coal Mine and Possible Remedial Measures Sakellariou M. G., A Study of Slope Stability Prediction Using Neural Networks. Geotechnical and Geological Engineering, 2005, 23(4)
60 B. Springer, Intelligent Forecast Procedures for Slope Stability with Evolutionary Artificial Neural Network. Computer Science, 2004, 31(74)
61 S. Iwata. Driving Force for Water Migration in Frozen Clayed Soil. Soil Science and Plant Nutrition, 1980, 26: 215-227
62 H. A. Dinulescu, E. R. G. Ecker. Analysis of the One-dimensional Moisture migration Caused by Temperature Gradients in a Porous Medium. International Journal of Heat and Mass Transfer, 1980, 23(8): 1069-1078.
63 E. Perfect, P. J. Williams. Thermally Induced Water Migration in Frozen Soils. Cold Regions Science and Technology, 1980, (3): 101-109
64 S. Krishnaiah, D. N. Singh. A Methodology to Determine Soil Moisture Movement due to Thermal Gradients. Experimental Thermal and Fluid Science. 2003, (27): 715-721
65 J. M. Konrad, N. R. Morgenstern. Prediction of Frost Heave in the Laboratory during Transient Freezing. Canadian Geotechnical Journal. 1982, 19: 250-259
66原国红.季节冻土水分迁移的机理及数值模拟.吉林大学博士学位论文, 2006
67 C. Dirksen, R. D. Miller. Closed-system Freezing of Unsaturated Soil. Soil Sci. Soc. Am. Proc, 1966, 30(2): 168-173
68王家澄.饱水砂土反复冻融时成冰条件的试验研究.冰川冻土, 1992, 14(2): 101-106
69李萍,徐学祖,王家澄等.不连续分凝冰发育规律的研究.兰州大学学报(自然科学版), 2000, 36(5): 126-133
70 R. L. Berg, G. L. Guymon, T. C. Johnson. Mathematical Model to Correlate Frost Heave of Pavements with Laboratory Predictions. USA Cold Regions Research and Engineering Laboratory, CRREL Report, 1980, 80(10)
71 G. L. Guymon, R. L. Berg, and T. V. Hromadka. Mathematical Model of Frost Heave and Thaw Settlement in Pavements.USA Cold Regions Research and Engineering Laboratory, CRREL Report, 1993, 93(2)
72 S. R. Bigl and S. A. Shoop. Soil Moisture Prediction during Freeze and Thaw Using a Coupled Heat and Moisture Flow Model. USA Cold Regions Research and Engineering Laboratory, CRREL Report, 1994, 111(7)
73 J. M. Konrad, N. R. Morgenstern. The Segregation Potential of a Freezing Soil . Canadian Geotechnical Journal, 1981, 18: 482-491
74 J. R. Philip, D. A.De Veries, Moisture Movement in Porous Material under Temperature Gradient. Trans Am Geophys Union, 1957, 38(2): 222-232
75 R. L. Harlan. Analysis of Coupled Heat-fluid Transport in Partially Frozen Soil . Water Resoures, 1973, 9(5): 1314-1323
76 M. Fremond, M.Mikkola. Thermo-dynamical Modeling of Freezing Soil. Proceedings of 6th International Symposium on Ground Freezing, 1991, 17-24
77徐学祖,邓友生.冻土中水分迁移的试验研究.北京科学出版社, 1991
78徐学祖,王家澄,邓友生.土中水分迁移试验装置的研制.冰川冻土, 1988, 10(2): 189-197
79 J. Wang, G. Cheng. Experimental Study on Conditions for Ice Formation of Saturated Sand in Repeated Freezing and Thawing Circulations. Journal of Glaciology and Geocryology, 1992, 14 (2): 101-106
80高永,胡春元,董智等.土壤冻结过程中水分迁移动向的研究.林业科学, 2000, 36(4): 126-128
81于琳琳.不同人工冻结条件下土的冻胀试验研究.哈尔滨工业大学博士学位论文, 2006: 9-10
82叶乐安,刘春平,劭明安.土壤水、热和溶质耦合运移研究进展.湖南师范大学自然科学学报, 2002, 25(2)
83 G. S. Taylor, J. N. Luthin. A Model for Coupled Heat and Moisture Transfer during Soil Freezing. Can. Geotech, 1978, 15: 548-555
84 G. N. Flerchinger, K. E. Saxton. Simultaneous Heat and Water Model of a Freezing Snow-residue-soil System. ASCE, 1989, 32 (2): 65-78
85 N. Li, B. Chen, F. Chen, X. Xu. The Coupled Heat-moisture-mechanic Model of the Frozen Soil. Cold Regions Science and Technology, 2000, 31 (3): 199-205
86赵建军,董金梅,王沛等.正冻土中的水热耦合模型.天津城市建设学院学报, 2001, 7(1): 47-52
87何平,程国栋,俞祁浩等.饱和冻土中水、热、力耦合模型.冰川冻土, 2001, 23(2): 135-138
88毛雪松,胡长顺,窦明健等.正冻土中水分场和温度场耦合过程的动态观测与分析.冰川冻土, 2003, 25(1): 54-59
89徐学祖,王家澄,邓友生.土中水分迁移试验装置的研制.冰川冻土, 1988, 10(2): 189-197
90徐学祖,邓友生.冻土中水分迁移的试验研究.北京科学出版社, 1991
91中国科学院兰州冰川冻土研究所.冻土的温度水分应力及其相互作用.兰州大学出版社, 1989
92叶佰生,陈肖柏.非饱和土冻结时水热耦合迁移的数值模拟.第四届全国冻土学术会议论文集,北京科学出版社, 1990: 77-89
93李述训,程国栋.冻融土中的水热输送问题.兰州大学出版社, 1995
94雷志栋,尚松浩,杨诗秀等.土壤冻结过程中潜水蒸发规律的模拟研究.水利学报, 1999: 6-9
95郑秀清,樊贵盛.冻融土壤水热迁移数值模型的建立及仿真分析.系统仿真学报, 2001, 13(3): 308-311
96毛雪松.多年冻土地区路基水热力场耦合效应研究.长安大学博士学位论文, 2004
97王铁行,胡长顺.多年冻土地区路基温度场和水分迁移场耦合问题研究.土木工程学报, 2003, 36(12): 93-97
98王铁行,贺再球,赵树德等.非饱和土体气态水迁移试验研究.岩石力学与工程学报, 2005, 24(18): 71-75
99王铁行,赵树德.非饱和土体气态水迁移引起的含水量变化方程.中国公路学报, 2003, 16(2): 18-21
100刘兵.土中水分在冻融过程中的迁移.哈尔滨工业大学硕士学位论文, 2008
101罗小刚,陈湘生,吴成义.冻融对土工参数影响的试验研究.建井技术, 2000, 21(2): 24-26
102杨成松,何平,程国栋等.冻融作用对土体干容重和含水量影响的试验研究.岩石力学与工程学报, 2003, 22(增2): 2695-2699
103杨成松,何平,程国栋等.冻融试验对土中含水量分布的影响.冰川冻土, 2004, 26(增刊): 50-55
104 D. M. Anderson, R. Pusch, E. Penner. Physical and Thermal Properties of Frozen Ground. Geotechnical Engineering for Cold Regions, McGra-Hill, 1978
105 P. J. Williams. Moisture Migration in Frozen Soils. National Academy Press, Washington D. C, 1984
106 S. K. J. Kung, R. S. Steenhuis. Heat and Moisture Transfer in a Partly Frozen Nonheaving Soil. Soil Sci. Soc, Am, 1986, 50: 1114-1122
107 V. J. Lunardini. Heat Transfer with Freezing and Thawing. Elsevier Amesterdam-Oxford-New York-Tokey, 1991
108 S. Li, G. Cheng. Experimental Study on Moisture and Heat Transfer for Lanzhou Yellow Soil during Freezing and Thawing Process. Annual Report on State Key Laboratory of Frozen Soil Engineering, CAREERI, CAS, 1992: 207-216
109孙建乐.非饱和土水分迁移试验研究.西部探矿工程, 2007, (7): 36-38
110盛煜,马巍,侯仲杰.正冻土中水分迁移的迁移势模型.冰川冻土, 1993, 15(1): 140-143
111 S. A. Shoop, S. R. Bigl. Moisture Migration during Freeze and Tthaw of Unsaturated Soils: Modeling and Large Scale Experiments. Cold Regions Science and Technology, 1997, 25: 33-45
112郭力,苗天德,张慧等.饱和正冻土中水热迁移的热力学模型.岩土工程学报, 1998, 20(5): 87-91
113苗天德,郭力,牛永红等.正冻土中水热迁移问题的混合物理论模型.中国科学(D辑), 1999, 29(增1): 8-14
114梁冰,刘晓丽,薛强.非等温入渗条件下土壤中水分运移的解析分析.辽宁工程技术大学学报, 2002, 21(6): 741-744
115毛雪松,王秉纲,胡长顺等.冻土路基水热迁移问题的理论模型及数值模拟.中外公路, 2006a, 26(1): 23-26
116毛雪松,王秉纲,胡长顺.多年冻土路基水分迁移热力学性能分析.路基工程, 2006b, 26(4): 1-3
117 M. Prat. Analysis of Experiments of Moisture Migration Caused by Temperature Differences in Unsaturated Porous Medium by Means of Two-Dimensional Numerical Simulation. International Journal of Heat and Mass Transfer, 1986,29(7): 1033-1039.
118 S. Somasundaram,Y. B. Suh, Win Aung. Analysis of Moisture Migration in Two-dimensional Unsaturated Porous Media with Impermeable Boundaries. International Journal of Heat and Mass Transfer, 1989, (32): 1733-1739
119 L. Lundin. Hydraulic Properties in an Operational Model of Frozen Soil, 1990,118: 289-310.
120张丰帆.季节性冻土构筑物冻胀机理研究及应用.吉林大学硕士学位论文, 2008
121郑秀清,樊贵盛,邢述彦.水分在季节性非饱和冻融土壤中的运动.北京地质出版社, 2002
122杜应吉.地质力学模型试验的研究现状与发展趋势.西北水资源与水工程, 1996, 7(2): 64-67
123靳德武,牛富俊,陈志新等.冻土斜坡模型试验相似分析.地球科学与环境学报, 2004,26(1): 29-31
124朱林楠,李东庆,郭兴民.无外荷载作用下冻土模型试验的相似分析.冰川冻土, 1993, 15(1): 166-169
125李东庆,朱林楠.荷载作用下冻土模型试验相似分析.自然科学进展——国家重点实验室通讯, 1994, 4(3): 323-327
126李之光.相似与模化理论及应用.北京国防工业出版社, 1982: 54-94
127魏先祥,赖远明.相似方法的原理及应用.兰州大学出版社, 2001, 1-38
128徐廷.相似理论与模型试验.中国农业机械出版社, 1982
129杨代泉,沈珠江.非饱和土孔隙气、水、汽、热耦合运动之模拟,岩土工程学报, 2005, 22(3): 357-361
130王绍令,杨梅学,小池俊雄等.时域反射仪在监测青藏高原活动层水分变化过程中的应用.冰川冻土, 2000, 22(1): 78-84
131徐昊.考虑土中水分迁移的路堑边坡冻融失稳数值分析.哈尔滨工业大学硕士学位论文, 2009
132 C. Chen et. Al. Analysis of Time-varying Rainfall Infiltration Induced Landslide, 2005, 48
133杨诗秀,雷志栋,朱强等.土壤冻结条件下水热耦合运移的数值模拟.清华大学学报, 1988, 28(1): 112-120
134 C. L. Shi, J. M. Yu. Numerical Simulation of Water and Heat Transfer during Freezing and Thawing in the Situation of Saturated Soil. China Agriculture Weather, 1998, 19 (4): 21-26
135佟治权等.砂粘土和粘砂土冻融交界面的抗剪强度试验研究报告.青藏铁路多年冻土科研成果论文集,中铁西北科学研究院, 2003
136佟治权.冻融交界面和冻后全融土现场大型直剪试验研究.第二届全国冻土学术会议论文选集,甘肃人民出版社, 1983, 5(1): 304-308
2吴玮江.季节性冻融作用与斜坡整体变形破坏.中国地质灾害与防治学报, 1996, 7(4): 59-64
3许珊珊,高伟.寒区公路路堑边坡抗冻融稳定性分析.低温建筑技术, 2006(2): 10-12
4 L. P. Joseph, Simulation and Field Validation of Heat and Water Flow during Soil Freezing and Thawing. Thesis for Doctoral Degree of Oregon State University, 1987
5王彩霞等.冰冻地区路堑边坡溜方原因分析与稳定性验算方法.东北公路, 2001, 24(1)
6中华人民共和国国家标准,冻土工程地质勘察规范GB50324,中国计划出版社, 2001
7杨润田,林凤桐.多年冻土区水文地质及工程地质学.哈尔滨东北林业大学出版社, 1986
8曹宏章.饱和颗粒土冻结过程中的多场耦合研究.中国科学院研究生院博士学位论文, 2006
9王晓春,张倬元.寒区工程与冻融力学.地学前缘, 2000, 7(增刊): 99-140
10程国栋.中国冻土研究近今进展.地理学报, 1990, 45(2): 220-223
11Г.П.马祖罗夫,冻土物理力学性质.北京煤炭工业出版社, 1980
12Э.Д.叶尔绍夫,工程冻土学,莫斯科大学出版社出版, 1999
13车茜.基于BP神经网络的季节冻土水分迁移机理研究.吉林大学硕士学位论文, 2008
14陈玉超.冻融环境下岩土边坡稳定性研究初探.西安科技大学土木工程学院, 2006: 3-20
15 R. L. Lytton, D. E. Pufahl, C. H. Michalak, H. S. Liang, and B. J. Dempsey, . An Integrated Model of The Climatic Effects on Pavements. Federal Highway Administration Report No. FHWA-RD-1990, 90(33)
16童长江,管枫年.土的冻胀与建筑冻害防治.水利电力出版社, 1985
17 H. A.崔托维奇,张长庆,朱元林译.冻土力学.北京科学出版社, 1985
18吴紫汪,马巍著.冻土强度与蠕变.兰州大学出版社, 1993
19周幼吾等.中国冻土,科学出版社, 2000
20周幼吾等.近40年来东北地区季节冻结与融化层温度变化与气候变暖.第
五届全国冰川冻土学大会论文集(上),甘肃文化出版, 1996: 3-10
21牛富俊,张鲁新,俞祁浩等.青藏高原多年冻土区斜坡类型及典型斜坡稳定性研究.冰川冻土, 2002, 24(5): 608-613
22丁德文,罗学波.冻土热工模型试验的理论基础.科学通报, 1979(8): 360-364
23 B. Y. Liu, M. A. Nearing, L. M. Risse. Slope Gradient Effects on Soil Loss for Slopes . Transactions of the ASAE, 1994, 37(6)
24 B. Y. Liu, M. A. Nearing, P. J. Shi, et al. Slope Length Effects on Soil Loss for Steep Slopes. Soil Sci. Soc. of Am. J., 2000, 64(5)
25铁道部建设司,青藏铁路多年冻土区工程勘察暂行规定, 2001
26铁道第一勘测设计院,西大滩至安多多年冻土区综合地质报告,中国铁道出版社, 2001
27铁道部建设司,青藏铁路多年冻土区工程设计暂行规定, 2002
28铁道第三勘察设计院,冻土工程,中国铁道出版社, 2002
29黄小铭等.青藏高原多年冻土区土冻结融化特征探讨.青藏铁路多年冻土科研成果论文集,中铁西北科学研究院, 2003
30中铁西北科学研究院,青藏铁路风火山多年冻土长期综合观测与试验研究,中国铁道出版社, 2006
31童长江等.青藏公路沿线冻土环境工程地质评价及冻土工程处理.第五届全国冰川冻土学大会论文集(上),甘肃文化出版社, 1996: 365-372
32王国尚等.青藏公路沿线冻土区自然环境及工程环境的变化.第五届全国冰川冻土学大会论文集(上),甘肃文化出版社, 1996: 43-50
33喻文学,武敬民.青藏公路多年冻土区沥青路面路基高度问题.西安公路学院学报,高原多年冻土区道路工程研究报告集,西安公路学院学报编辑部出版, 1983: 25-47
34章金钊.青藏公路沥青路面冻土路基设计原则探讨.第五届全国冰川冻土学大会论文集(上),甘肃文化出版社, 1996: 383-388
35李杨.季节冻土水分迁移模型研究.吉林大学博士学位论文, 2008
36赵安平.季冻区路基土冻胀的微观机理研究.吉林大学博士学位论文, 2008
37李杰.边坡稳定性问题的理论分析及其应用研究.大连理工大学硕士学位论文, 2002
38李雨浓.季冻区高速公路路基冻害影响因素分析.吉林大学硕士学位论文, 2007
39刘鸿绪.对冻结过程力学中若干问题的商榷.黑龙江省低温建筑科学研究所, 1988
40刘鸿绪.建筑基础的冻胀力.第五届全国冰川冻土学大会论文集(下),甘肃文化出版社, 1996
41武鹤,高伟,王国峰等.寒区公路土质路堑边坡滑塌原因及其防治.自然灾害学报, 2006, 15(3): 66-70
42宋高嵩,陈阳东,温福春等.季冻区路基边坡的水稳定性影响分析与防治.公路交通科技, 2005 , 22(6): 65-67
43 A. D. Abrahams, et. al., Hydraulics on a Semiarid Hill Slope. Southern Arizon Earth Surface Processes and Land Forms(21), 1996
44戴惠民,王兴隆.季节冻土区公路桥涵地基土冻胀性的研究.中国公路学报, 1989, 2(4): 18-25
45童长江等,青藏公路沿线冻土环境工程地质评价及冻土工程处理.第五届全国冰川冻土学大会论文集(下),甘肃文化出版社, 1996: 65-72
46雷汉忠.冻融土的工程性质及利用.煤矿设计, 1994(6): 39-44
47陈肖柏,王雅卿,刘建坤等.土的冻结作用与地基.北京科学出版社, 2006: 3-60
48杨平.原状土与冻融土物理力学性能差异性研究.南京林业大学学报, 2001, 25(2): 68-70
49杨平等.人工冻融土物理力学性能研究.冰川冻土, 2002, 24(5): 665-667
50武鹤等.寒区路堑人工土质边坡滑塌原因与稳定技术研究.黑龙江工程学院学报, 2005, 19(2)
51朱选举.冻土中未冻水迁移规律研究.清华大学硕士学位论文, 2000
52郭占荣,荆恩春,聂振龙.冻结期和冻融期土壤水分运移特征分析.水科学进展, 2002, 13(3): 298-303
53梁若筠.基于神经网络方法的冻土水分迁移研究.兰州大学硕士学位论文, 2007
54 R. D. Miller. Soil Freezing in Relation to Pore Water Pressure and Temperature in Permafrost. The North American Contribution to the 2nd International Conference on Permafrost, Yakutsk. National Academy of Sciences, Washington, D. C., 1973
55徐敩祖.国外对冻土中水分迁移课题的研究.冰川冻土, 1982, 4(3): 97-103
56徐敩祖,邓友生.饱水正冻土中水分迁移特征研究.第三届全国冻土学术会议论文选集,科学出版社, 1989: 204-210
57邓友生,徐学祖.非饱和土壤水扩散系数的变化特征.冰川冻土, 1994, 16 (3): 251-258
58徐敩祖,王家澄,张立新.冻土物理学.北京科学出版社, 2001
59 R. Ulusay, C. G?k?eo?lu, H. S?nmez, et al. Causes Mechanism and Environmental Impacts of Instabilities at Himmetoglu Coal Mine and Possible Remedial Measures Sakellariou M. G., A Study of Slope Stability Prediction Using Neural Networks. Geotechnical and Geological Engineering, 2005, 23(4)
60 B. Springer, Intelligent Forecast Procedures for Slope Stability with Evolutionary Artificial Neural Network. Computer Science, 2004, 31(74)
61 S. Iwata. Driving Force for Water Migration in Frozen Clayed Soil. Soil Science and Plant Nutrition, 1980, 26: 215-227
62 H. A. Dinulescu, E. R. G. Ecker. Analysis of the One-dimensional Moisture migration Caused by Temperature Gradients in a Porous Medium. International Journal of Heat and Mass Transfer, 1980, 23(8): 1069-1078.
63 E. Perfect, P. J. Williams. Thermally Induced Water Migration in Frozen Soils. Cold Regions Science and Technology, 1980, (3): 101-109
64 S. Krishnaiah, D. N. Singh. A Methodology to Determine Soil Moisture Movement due to Thermal Gradients. Experimental Thermal and Fluid Science. 2003, (27): 715-721
65 J. M. Konrad, N. R. Morgenstern. Prediction of Frost Heave in the Laboratory during Transient Freezing. Canadian Geotechnical Journal. 1982, 19: 250-259
66原国红.季节冻土水分迁移的机理及数值模拟.吉林大学博士学位论文, 2006
67 C. Dirksen, R. D. Miller. Closed-system Freezing of Unsaturated Soil. Soil Sci. Soc. Am. Proc, 1966, 30(2): 168-173
68王家澄.饱水砂土反复冻融时成冰条件的试验研究.冰川冻土, 1992, 14(2): 101-106
69李萍,徐学祖,王家澄等.不连续分凝冰发育规律的研究.兰州大学学报(自然科学版), 2000, 36(5): 126-133
70 R. L. Berg, G. L. Guymon, T. C. Johnson. Mathematical Model to Correlate Frost Heave of Pavements with Laboratory Predictions. USA Cold Regions Research and Engineering Laboratory, CRREL Report, 1980, 80(10)
71 G. L. Guymon, R. L. Berg, and T. V. Hromadka. Mathematical Model of Frost Heave and Thaw Settlement in Pavements.USA Cold Regions Research and Engineering Laboratory, CRREL Report, 1993, 93(2)
72 S. R. Bigl and S. A. Shoop. Soil Moisture Prediction during Freeze and Thaw Using a Coupled Heat and Moisture Flow Model. USA Cold Regions Research and Engineering Laboratory, CRREL Report, 1994, 111(7)
73 J. M. Konrad, N. R. Morgenstern. The Segregation Potential of a Freezing Soil . Canadian Geotechnical Journal, 1981, 18: 482-491
74 J. R. Philip, D. A.De Veries, Moisture Movement in Porous Material under Temperature Gradient. Trans Am Geophys Union, 1957, 38(2): 222-232
75 R. L. Harlan. Analysis of Coupled Heat-fluid Transport in Partially Frozen Soil . Water Resoures, 1973, 9(5): 1314-1323
76 M. Fremond, M.Mikkola. Thermo-dynamical Modeling of Freezing Soil. Proceedings of 6th International Symposium on Ground Freezing, 1991, 17-24
77徐学祖,邓友生.冻土中水分迁移的试验研究.北京科学出版社, 1991
78徐学祖,王家澄,邓友生.土中水分迁移试验装置的研制.冰川冻土, 1988, 10(2): 189-197
79 J. Wang, G. Cheng. Experimental Study on Conditions for Ice Formation of Saturated Sand in Repeated Freezing and Thawing Circulations. Journal of Glaciology and Geocryology, 1992, 14 (2): 101-106
80高永,胡春元,董智等.土壤冻结过程中水分迁移动向的研究.林业科学, 2000, 36(4): 126-128
81于琳琳.不同人工冻结条件下土的冻胀试验研究.哈尔滨工业大学博士学位论文, 2006: 9-10
82叶乐安,刘春平,劭明安.土壤水、热和溶质耦合运移研究进展.湖南师范大学自然科学学报, 2002, 25(2)
83 G. S. Taylor, J. N. Luthin. A Model for Coupled Heat and Moisture Transfer during Soil Freezing. Can. Geotech, 1978, 15: 548-555
84 G. N. Flerchinger, K. E. Saxton. Simultaneous Heat and Water Model of a Freezing Snow-residue-soil System. ASCE, 1989, 32 (2): 65-78
85 N. Li, B. Chen, F. Chen, X. Xu. The Coupled Heat-moisture-mechanic Model of the Frozen Soil. Cold Regions Science and Technology, 2000, 31 (3): 199-205
86赵建军,董金梅,王沛等.正冻土中的水热耦合模型.天津城市建设学院学报, 2001, 7(1): 47-52
87何平,程国栋,俞祁浩等.饱和冻土中水、热、力耦合模型.冰川冻土, 2001, 23(2): 135-138
88毛雪松,胡长顺,窦明健等.正冻土中水分场和温度场耦合过程的动态观测与分析.冰川冻土, 2003, 25(1): 54-59
89徐学祖,王家澄,邓友生.土中水分迁移试验装置的研制.冰川冻土, 1988, 10(2): 189-197
90徐学祖,邓友生.冻土中水分迁移的试验研究.北京科学出版社, 1991
91中国科学院兰州冰川冻土研究所.冻土的温度水分应力及其相互作用.兰州大学出版社, 1989
92叶佰生,陈肖柏.非饱和土冻结时水热耦合迁移的数值模拟.第四届全国冻土学术会议论文集,北京科学出版社, 1990: 77-89
93李述训,程国栋.冻融土中的水热输送问题.兰州大学出版社, 1995
94雷志栋,尚松浩,杨诗秀等.土壤冻结过程中潜水蒸发规律的模拟研究.水利学报, 1999: 6-9
95郑秀清,樊贵盛.冻融土壤水热迁移数值模型的建立及仿真分析.系统仿真学报, 2001, 13(3): 308-311
96毛雪松.多年冻土地区路基水热力场耦合效应研究.长安大学博士学位论文, 2004
97王铁行,胡长顺.多年冻土地区路基温度场和水分迁移场耦合问题研究.土木工程学报, 2003, 36(12): 93-97
98王铁行,贺再球,赵树德等.非饱和土体气态水迁移试验研究.岩石力学与工程学报, 2005, 24(18): 71-75
99王铁行,赵树德.非饱和土体气态水迁移引起的含水量变化方程.中国公路学报, 2003, 16(2): 18-21
100刘兵.土中水分在冻融过程中的迁移.哈尔滨工业大学硕士学位论文, 2008
101罗小刚,陈湘生,吴成义.冻融对土工参数影响的试验研究.建井技术, 2000, 21(2): 24-26
102杨成松,何平,程国栋等.冻融作用对土体干容重和含水量影响的试验研究.岩石力学与工程学报, 2003, 22(增2): 2695-2699
103杨成松,何平,程国栋等.冻融试验对土中含水量分布的影响.冰川冻土, 2004, 26(增刊): 50-55
104 D. M. Anderson, R. Pusch, E. Penner. Physical and Thermal Properties of Frozen Ground. Geotechnical Engineering for Cold Regions, McGra-Hill, 1978
105 P. J. Williams. Moisture Migration in Frozen Soils. National Academy Press, Washington D. C, 1984
106 S. K. J. Kung, R. S. Steenhuis. Heat and Moisture Transfer in a Partly Frozen Nonheaving Soil. Soil Sci. Soc, Am, 1986, 50: 1114-1122
107 V. J. Lunardini. Heat Transfer with Freezing and Thawing. Elsevier Amesterdam-Oxford-New York-Tokey, 1991
108 S. Li, G. Cheng. Experimental Study on Moisture and Heat Transfer for Lanzhou Yellow Soil during Freezing and Thawing Process. Annual Report on State Key Laboratory of Frozen Soil Engineering, CAREERI, CAS, 1992: 207-216
109孙建乐.非饱和土水分迁移试验研究.西部探矿工程, 2007, (7): 36-38
110盛煜,马巍,侯仲杰.正冻土中水分迁移的迁移势模型.冰川冻土, 1993, 15(1): 140-143
111 S. A. Shoop, S. R. Bigl. Moisture Migration during Freeze and Tthaw of Unsaturated Soils: Modeling and Large Scale Experiments. Cold Regions Science and Technology, 1997, 25: 33-45
112郭力,苗天德,张慧等.饱和正冻土中水热迁移的热力学模型.岩土工程学报, 1998, 20(5): 87-91
113苗天德,郭力,牛永红等.正冻土中水热迁移问题的混合物理论模型.中国科学(D辑), 1999, 29(增1): 8-14
114梁冰,刘晓丽,薛强.非等温入渗条件下土壤中水分运移的解析分析.辽宁工程技术大学学报, 2002, 21(6): 741-744
115毛雪松,王秉纲,胡长顺等.冻土路基水热迁移问题的理论模型及数值模拟.中外公路, 2006a, 26(1): 23-26
116毛雪松,王秉纲,胡长顺.多年冻土路基水分迁移热力学性能分析.路基工程, 2006b, 26(4): 1-3
117 M. Prat. Analysis of Experiments of Moisture Migration Caused by Temperature Differences in Unsaturated Porous Medium by Means of Two-Dimensional Numerical Simulation. International Journal of Heat and Mass Transfer, 1986,29(7): 1033-1039.
118 S. Somasundaram,Y. B. Suh, Win Aung. Analysis of Moisture Migration in Two-dimensional Unsaturated Porous Media with Impermeable Boundaries. International Journal of Heat and Mass Transfer, 1989, (32): 1733-1739
119 L. Lundin. Hydraulic Properties in an Operational Model of Frozen Soil, 1990,118: 289-310.
120张丰帆.季节性冻土构筑物冻胀机理研究及应用.吉林大学硕士学位论文, 2008
121郑秀清,樊贵盛,邢述彦.水分在季节性非饱和冻融土壤中的运动.北京地质出版社, 2002
122杜应吉.地质力学模型试验的研究现状与发展趋势.西北水资源与水工程, 1996, 7(2): 64-67
123靳德武,牛富俊,陈志新等.冻土斜坡模型试验相似分析.地球科学与环境学报, 2004,26(1): 29-31
124朱林楠,李东庆,郭兴民.无外荷载作用下冻土模型试验的相似分析.冰川冻土, 1993, 15(1): 166-169
125李东庆,朱林楠.荷载作用下冻土模型试验相似分析.自然科学进展——国家重点实验室通讯, 1994, 4(3): 323-327
126李之光.相似与模化理论及应用.北京国防工业出版社, 1982: 54-94
127魏先祥,赖远明.相似方法的原理及应用.兰州大学出版社, 2001, 1-38
128徐廷.相似理论与模型试验.中国农业机械出版社, 1982
129杨代泉,沈珠江.非饱和土孔隙气、水、汽、热耦合运动之模拟,岩土工程学报, 2005, 22(3): 357-361
130王绍令,杨梅学,小池俊雄等.时域反射仪在监测青藏高原活动层水分变化过程中的应用.冰川冻土, 2000, 22(1): 78-84
131徐昊.考虑土中水分迁移的路堑边坡冻融失稳数值分析.哈尔滨工业大学硕士学位论文, 2009
132 C. Chen et. Al. Analysis of Time-varying Rainfall Infiltration Induced Landslide, 2005, 48
133杨诗秀,雷志栋,朱强等.土壤冻结条件下水热耦合运移的数值模拟.清华大学学报, 1988, 28(1): 112-120
134 C. L. Shi, J. M. Yu. Numerical Simulation of Water and Heat Transfer during Freezing and Thawing in the Situation of Saturated Soil. China Agriculture Weather, 1998, 19 (4): 21-26
135佟治权等.砂粘土和粘砂土冻融交界面的抗剪强度试验研究报告.青藏铁路多年冻土科研成果论文集,中铁西北科学研究院, 2003
136佟治权.冻融交界面和冻后全融土现场大型直剪试验研究.第二届全国冻土学术会议论文选集,甘肃人民出版社, 1983, 5(1): 304-308