基于细观均匀化的膨润土有效热传导特性研究

基于均匀化方法和椭球夹杂问题基本解给出了非饱和岩土体的有效热传导特性模型。该模型可以考虑夹杂形态、体积分数和空间分布及夹杂间相互作用对有效热传导特性的影响,反映了由于非均匀夹杂引起有效热传导张量的各向异性特性。讨论了夹杂宽高比以及夹杂与基质热传导系数比对有效热传导张量各向异性系数的影响。将岩土体看作固体基质和孔隙夹杂构成的非均匀材料,探讨了岩土体有效热传导系数随孔隙形态、孔隙率和饱和度的变化规律。最后,应用上述模型对高庙子膨润土(GMZ01)有效热传导系数进行预测并与其他模型预测结果进行对比分析。研究表明,本文模型对GMZ01膨润土有效热传导系数具有较好的预测能力,但更准确的预测需根据膨润土的孔隙结构采用多层次均匀化方法。研究成果对于高放核废料深地质处置库缓冲材料的热-水-力耦合特性具有一定参考价值。     

基于微观结构重塑的非饱和冻土导热系数预测

为预测非饱和冻土的导热性能,基于土体微观结构,提出了非饱和冻土特征结构识别算法和多元素生成算法,并将该算法与传统有限单元法组合,建立非饱和冻土导热系数蒙特卡洛预测模型。通过土体SEM电镜图像,采用逆向四参数增长识别法识别土体中各组分含量、大小以及各方向分布概率;改进传统的四参数随机增长法,提出了考虑土、水、冰和气的多元素生成算法;基于生成的非饱和冻土模型,通过蒙特卡洛方法获得非饱和冻土导热系数,并与规范中冻土导热系数进行对比,验证了蒙特卡洛法预测模型的合理性（平均误差<4%）;通过多因素分析研究孔隙率、颗粒大小、土体导热性、饱和度以及结冰率对非饱和冻土导热性影响,各因素与导热系数的相关系数依次为：-0.352、-0.098、0.641、0.520和0.060,影响大小为：土颗粒导热性>饱和度>孔隙率>土颗粒大小>结冰率。各影响因素对非饱和冻土导热系数影响可以归纳为对热通量形成“热链”密度、宽度、连通性、热流承载力以及对“热桥”通量的影响。     

非饱和土热力学参数确定的探讨

研究非饱和土热力学的力学性质,必须建立能量守恒方程,表征能量的两个基本物理力学参数（体积比热系数和热传导系数）的确定是求解能量守恒方程的重要因素。对于不同颗粒组成的土体,对其干燥颗粒有效热力学参数进行了分析,对固体、液体和气体三相共存的热力学参数进行了研究,并对计算结果进行了对比分析,在此基础上,提出了较适合计算热传导系数的方法,并将此方法和体积平均法、自洽法及Hashin-Strikman方法进行了比较。经与试验结果对比表明：提出的方法能较好地描述非饱和土的热力学参数特性。     

非饱和粉质粘土固结压缩特性及体变试验研究

为进一步研究与基质吸力相关联的非饱和土固结压缩特性,扩展非饱和土固结理论在工程实际中的适用性,通过分析非饱和粉质粘土水土特征曲线变化规律,对非饱和土固结变形机理进行研究。试验结果表明：非饱和土的最终沉降量仅与土骨架的压缩模量有关。对于不同饱和度的非饱和土而言,固结速度随初始饱和度的增加而减小,饱和土固结过程所需要的时间比非饱和土固结过程所需要的时间短。由于孔隙流体的压缩性导致高饱和土体的瞬时沉降比低饱和土体的瞬时沉降小,但高饱和土体的后期固结沉降受饱和度和吸力的影响,比低饱和土体要大得多。      

非饱和重塑黄土基质吸力变化特征与物理性质的关系

本文以兰州重塑黄土为对象,通过室内试验分析6组非饱和重塑黄土土-水特征曲线分段特征值与主要物性指标的统计关系,探讨了非饱和重塑黄土物理性质对其基质吸力变化特征的影响。结果表明,对于非饱和重塑黄土,分段特征含水量、特征饱和度、特征基质吸力与干密度关系密切,与颗粒级配关系较弱。随干密度增加,过渡区特征饱和度、残余饱和度线性增加,饱和含水量、过渡区特征含水量线性降低,残余基质吸力及过渡区特征基质吸力分别为线性及非线性增加。过渡区特征饱和度、残余饱和度与土体中粘粒含量线性正相关,与粉粒含量负相关。饱和含水量与粉粒含量线性正相关,残余基质吸力与其线性负相关。重塑黄土的微观结构显示这些关系受制于土体的孔隙结构和颗粒间的相互关系。     

考虑残余含气量的非饱和土的水力耦合本构模型

大量的非饱和土干湿循环试验表明,当土体处于吸湿过程直至吸力降低为0 kPa时,土体并不能达到完全饱和状态,还存在一定的残余气体。在高饱和度时,由于残余气体以封闭气泡的形式分布在土体中,土体呈现较大的压缩性,使其与饱和土的性质不同。在这种状态下,现有的非饱和土本构模型预测到的土饱和度为100%,与试验结果存在一定的偏差。为了使本构模型在高饱和度状态时具有较高的精度,对非饱和土的毛细滞回和塑性变形耦合本构模型进行了修正,使其能够考虑残余含气量的影响。通过预测与实测结果比较,证明了新模型能够有效地模拟残余含气量对非饱和土力学特性的影响。     

一个高饱和度非饱和土的本构模型

高饱和度的非饱和土中由于气体处于封闭状态,其内部气压的变化必将对土体的行为产生影响。首先,对高饱和度非饱和土特性进行探讨和研究,随后,在已有非饱和土模型框架基础上,采用广义有效应力原理,建立一个适用于高饱和度条件下的非饱和土的弹塑性本构模型。模型中引入气相耗散的影响,在硬化方程中考虑封闭气体压力改变的影响。最后,利用已有的试验结果来对模型进行验证,并将模型预测结果与前人模型进行对比,表明模型预测可以很好地预测土体的行为,尤其是在高饱和度条件下其结果比其他模型更加接近实际情况。     

非饱和土中圆柱形衬砌隧道的瞬态动力响应

在以往关于圆柱形衬砌隧道的瞬态动力响应中,衬砌周围土体大多假定为弹性介质或饱和介质。然而,自然界中的土体大多为非饱和介质。考虑土体与衬砌结构的动力相互作用及动荷载引起的附加质量密度的影响,研究了瞬态荷载作用下非饱和土中无限长深埋圆柱形衬砌隧道的动力响应。基于多孔介质混合物理论和连续介质力学理论,建立了非饱和土中圆柱形衬砌隧道受到瞬态荷载作用时衬砌及周围土体的控制方程,利用Durbin数值反演法得到了衬砌及土体在时间域的动力响应。数值分析了饱和度对瞬态荷载下径向位移、径向应力、环向应力和孔隙水压力的影响。结果表明：饱和度对衬砌及周围土体的瞬态响应影响显著;饱和度对径向位移沿径向的衰减影响较小,对环向应力和孔隙压力沿径向的衰减影响较大。     

PS加固对非饱和遗址土的渗透特性影响研究

高模数硅酸钾（简称PS）对西北干旱区的土遗址具有明显的加固效果,其加固后的土体也具有较好的渗透性。采用VJ-T非饱和土渗透仪和常规压力板仪测试了PS加固前后饱和样的渗透系数和非饱和土的土-水特征曲线(SWCC),并利用相关模型预测了非饱和土的渗透系数。结果显示,PS加固以后,在各个含水条件下的渗透性都有一定程度的提高,7%PS加固样表现非常明显。其原因是PS与土颗粒进行了非常复杂的物理化学变化,在反应中颗粒的几何形状发生变化,大量棱角弱化或消失,土颗粒磨圆度变好。同时,孔隙壁变得光滑、平直,孔隙弯曲因子减小,有效孔隙增大,这些因素导致了非饱和黏土的渗透性能有较大幅度提高。电镜测试分析了PS加固以后黏土颗粒及孔隙的变化,较好地验证了渗透性试验的结果     

基于分形理论的SWCC边界曲线滞后效应模型研究

微观土颗粒及孔隙分布的非均匀性及由此引起的瓶颈效应是造成非饱和土土-水特征曲线(SWCC)滞后效应的主要原因。引入分形理论,考虑非饱和土孔径及渗流路径的微观分形特性,提出了一个用于描述水在非饱和土中渗流的毛细管模型。模型中将非饱和土孔隙简化为一系列具有不同孔径大小的毛细弯管,其孔径大小及弯曲程度假定服从分形规律。在此基础上,推导得了非饱和土的吸湿与脱湿过程的饱和度S_e~-水头高度h来描述土-水特征曲线滞后效应的特征方程以及饱和度S_e~-相对水力传导系数Kr特征方程。与室内观测结果及已有研究的对比表明,该模型相比以往方法,可更好地模拟非饱和土土-水特征曲线的滞后效应。对非饱和土吸湿与脱湿过程滞后效应的本质进行了对比分析,揭示了滞后效应产生的根本原因在于土体中流通孔隙大小的非均匀性。     

非饱和（冻）土导热系数预估模型研究

岩土材料的导热系数是岩土工程温度场分析及建筑热工计算中的重要参数。研究旨在建立一个基于归一化导热系数概念和以土的干燥和饱和状态导热系数为基准的非饱和土导热系数的通用预估模型。通过对文献中328组实测数据的分析发现,将同类土在不同密实度条件下的各种导热系数-含水率曲线簇进行归一化处理后,可以得到惟一的归一化导热系数kr与饱和度Sr（归一化含水率）关系,1/kr与1/Sr呈相关性较好的线性关系,而每支1/kr-1/Sr直线均通过坐标（1,1）点的斜率由土质类型决定。据此提出了一个集成土质类型、密实度（孔隙率）和含水率（饱和度）等因素综合影响的融土和冻土导热系数通用预估模型,并给出了导热系数预估模型中土质参数的取值范围,以及融土和冻土处于完全干燥状态和饱和状态的确定方法。对预估模型进行验证结果表明,所提出的非饱和土导热系数预估模型具有较好准确性。     

Improving a thermal conductivity model of unsaturated soils based on multivariate distribution analysis

Soil thermal conductivity (k) is a key parameter for the design of energy geo-structures, and it depends on many soil properties such as saturation degree, porosity, mineralogical composition, soil type and others. Capturing these diversified influencing factors in a soil thermal conductivity model is a challenging task for engineers due to the nonlinear dependencies. In this study, a multivariate distribution approach was utilized to improve an existing soil thermal conductivity model, Cote and Konrad model, by quantitatively considering the impacts of dry density (ρ_d), porosity (n), saturation degree (S_r), quartz content (m_q), sand content (m_s) and clay content (m_c) on thermal conductivity of unsaturated soils. A large database containing these seven soil parameters was compiled from the literature to support the multivariate analysis. Simplified bivariate and multivariate correlations for improving the Cote and Konrad model were derived analytically and numerically to consider different influencing factors. By incorporating these simplified correlations, the predicted k values were more concentrated around the measured values with the coefficient of determination (R²) increased from 0.83 to 0.95. It is concluded that the developed correlations with the information of different soil properties provide an efficient, rational and simple way to predict soil thermal conductivity more accurately. Moreover, the quartz content is a more important factor than the porosity that shall be considered in the establishment of thermal conductivity models for unsaturated soils with high quartz content.

     

Simulation of liquefaction of unsaturated soil using critical state soil model

Several researchers have reported that the mean effective stress of unsaturated soils having a relatively high degree of saturation gradually decreases under fully undrained cyclic loading conditions, and such soils can be finally liquefied like saturated soils. This paper describes a series of simulations of fully undrained cyclic loading on unsaturated soils, conducted using an elastoplastic model for unsaturated soils. This model is a critical state soil model formulated using effective stress tensor for unsaturated soils, which incorporates the following concepts: (a) the volumetric movement of the state boundary surface containing the critical state line owing to the variation in the degree of saturation; (b) the soil water characteristic curve considering the effects of specific volume and hydraulic hysteresis; and (c) the subloading surface concept for considering the effect of density. Void air is assumed to be an ideal gas obeying Boyle's law. The proposed model is validated through comparisons with past results. The simulation results show that the proposed model properly describes the fully undrained cyclic behavior of unsaturated soils, such as liquefaction, compression, and an increase in the degree of saturation. Finally, the effects of the degree of saturation, void ratio, and confining pressure on the cyclic strength of unsaturated soils are described by the simulation results. The liquefaction resistance of unsaturated soils increases as the degree of saturation and the void ratio decrease, and as the confining pressure increases. Furthermore, the degree of saturation has a greater effect on the liquefaction resistance than the confining pressure and void ratio. Copyright © 2017 John Wiley & Sons, Ltd.     

基于细观模拟的软土导热系数数值预测模型

热环境控制是城市地下空间安全运行的有力保障,科学预测地下围岩的导热性能是地下空间环控系统热负荷评估的基础。软土是地下空间开发中一类常见土体,现有模型主要适合于预测中低含水率范围内土体导热系数的变化,而对于高含水率软土,合适的导热系数预测模型较少。基于细观模拟,本文提出了一种能有效预测软土导热系数的数值模型。该模型除了能够反映含水率、干密度等常规因素影响外,还可考虑矿物组成以及粒径分布等的影响。最后,通过与苏通GIL管廊工程中20个软土样的实测导热系数进行对比以验证数值模型的可靠性。首先借助激光粒度分析仪和X射线衍射分析测试了矿物组成和粒度分布,代入模型进行数字建模并通过细观导热模拟得到导热系数预测值。导热系数模拟预测值与实测值对比结果显示:模拟预测值基本在实测值±20%范围内,验证了本文模型的可靠性,表明了该模型在预测高含水率软土导热系数方面的潜力。此外,该模型还可以直观地展示土内各处局部热流的分布特征,这为深入认识土体导热行为的机理奠定了基础。本文研究可为软土以及土体导热系数的预测评价提供新的思路和方法。     

Correlation of unsaturated soil and dielectric property for monitoring of subsurface characteristics: development of unsaturated dielectric mixing models and its application

The behavior due to rainfall infiltrating the ground plays a role in landslides, groundwater recharge and various other ground responses. Most of these geotechnical behaviors have a correlation between soil pore space and soil volumetric water content in the unsaturated and saturated soil porous media. Therefore, the soil porosity associated with soil pores and the distribution of volumetric water content are significantly important hydrological characteristics. In the case of shallow slope failure such as landslide, the infiltration activity due to the connectivity of soil pore spaces in a porous media is induced. Slope failure may be attributed to the effect of a wetting front with the slope due to liquid infiltration, which changes the volumetric water content, soil matric suction and shear strength of the slope. This study was performed with an unsaturated injection test using a frequency domain reflectometry (FDR) dielectric device which measures the dielectric constant of unsaturated soil and the study then proposed the unsaturated dielectric mixing models to calculate soil porosity and effective porosity of unsaturated soils. From the experimental results the ratio of effective porosity to porosity of soils are measured in a range of 70–85%. These experimental results show a decrease of about 5–10% for unsaturated soil compared to the ratio of effective porosity to porosity of saturated soil. The infiltration passages of tracer material are restricted within the pore connectivity in the unsaturated soil which is caused by dead-pores in the soil. Using the FDR device and the unsaturated dielectric mixing models, we can consider the acquisition of physical properties to detect the infiltration activity, the response of the dielectric constant along with the injected tracer and hydrological parameters for the unsaturated soil porous media.     

A macroscopic model for predicting the relative hydraulic permeability of unsaturated soils

A macroscopic model for predicting the relative hydraulic permeability of unsaturated soils is proposed. In this model, pores in unsaturated soils are considered to be parallel flow tubes. The water flow in the pores is assumed to take place in the water film on the inside wall of the flow tubes. The viscosity of pore water is considered to be different from the viscosity of pure water and variable with the variation of degree of saturation. The values of tortuosity factor and pore shape factor of unsaturated soils are estimated theoretically. The theoretical model is verified using experimental data for 32 different soils. For application in engineering practice, the value of viscosity of pore water in different soils is proposed.     

Investigation of the effect of specific interfacial area on strength of unsaturated granular materials by X-ray tomography

This paper studies the effect of interfacial areas (air–water interfaces and solid–water interfaces) on material strength of unsaturated granular materials. High-resolution X-ray computed tomography technique is employed to measure the interfacial areas in wet glass bead samples. The scanned 3D images are trinarized into three phases and meshed into representative volume elements (RVEs). An appropriate RVE size is selected to represent adequate local information. Due to the local heterogeneity of the material, the discretized RVEs of the scanned samples actually cover a very large range of degree of saturation and porosity. The data of RVEs present the relationship between the specific interfacial areas and degree of saturation and gives boundaries where the interfacial area of a whole sample should fall in. In parallel, suction-controlled direct shear tests have been carried out on glass beads and the material strength has been corroborated with two effective stress definitions related to the specific air–water interfacial areas and fraction of wetted solid surface, respectively. The comparisons show that the specific air–water interfacial area reaches the peak at about 25% of saturation and contributes significantly to the material strength (up to 60% of the total capillary strength). The wetted solid surface obtained from X-ray CT is also used to estimate Bishop’s coefficient χ based on the second type of effective stress definition, which shows a good agreement with the measured value. This work emphasizes the importance to include interface terms in effective stress formulations of unsaturated soils. It also suggests that the X-ray CT technique and RVE-based multiscale analysis are very valuable in the studies of multiphase geomaterials.

     

Desiccation shrinkage of non‐clayey soils: a numerical study

A mesoscale model of desiccation of soil based on the evolution of the pore system idealized as bimodal is numerically examined. A simplified evolution of the model reveals a series of characteristics that qualitatively agree with the observed macroscopic experimental findings. The principal mechanism is deemed to be driven by the surface evaporation and water outflow generating a pore pressure gradient resulting in the shrinkage mainly of the largest pores. The amount of shrinkage is a function of (negative) pore pressure and is controlled by the compressibility of the solid matrix. The numerical model includes also the ensuing partial saturation stage initiated by the air entry simulated as a scenario with a moving phase interface inside the pore. The proposed model can be extended beyond the two‐mode porosity soils, to include the multi‐modal porosity, or its statistical representation.Copyright © 2012 John Wiley & Sons, Ltd.     

A hierarchical thermo‐hydro‐plastic constitutive model for unsaturated soils and its numerical implementation

Unsaturated soils are highly heterogeneous 3‐phase porous media. Variations of temperature, the degree of saturation, and density have dramatic impacts on the hydro‐mechanical behavior of unsaturated soils. To model all these features, we present a thermo‐hydro‐plastic model in which the hydro‐mechanical hardening and thermal softening are incorporated in a hierarchical fashion for unsaturated soils. This novel constitutive model can capture heterogeneities in density, suction, the degree of saturation, and temperature. Specifically, this constitutive model has 2 ingredients: (1) it has a “mesoscale” mechanical state variable—porosity and 3 environmental state variables—suction, the degree of saturation, and temperature; (2) both temperature and mechanical effects on water retention properties are taken into account. The return mapping algorithm is applied to implement this model at Gauss point assuming an infinitesimal strain. At each time step, the return mapping is conducted only in principal elastic strain space, assuming no return mapping in suction and temperature. The numerical results obtained by this constitutive model are compared with the experimental results. It shows that the proposed model can simulate the thermo‐hydro‐mechanical behavior of unsaturated soils with satisfaction. We also conduct shear band analysis of an unsaturated soil specimen under plane strain condition to demonstrate the impact of temperature variation on shear banding triggered by initial material heterogeneities.