Strength and structure of Q3paleosol
-
摘要: 古土壤是黄土地层的重要组成部分,但目前对古土壤力学性质的研究相对较少。为研究古土壤的强度及其结构性,以西安某地Q3古土壤为研究对象,通过土力学途径,对原状、重塑、饱和原状古土壤进行无侧限抗压强度试验,得到了相应状态古土壤的无侧限抗压强度,计算得出其初始结构性指标；通过复合幂指数模型,得到了原状、重塑古土壤无侧限抗压强度随含水率变化的表达式,分析了原状、重塑古土壤无侧限抗压强度、扰动灵敏度、浸水灵敏度及构度指标的变化规律。结果表明：古土壤的结构性强,在饱和状态时,其破坏形态仍然为脆性破坏。在密度、湿度、粒度相同条件下,原状古土壤的强度小于重塑古土壤的强度,其扰动灵敏度小于1,与一般黄土的规律截然相反；古土壤的浸水灵敏度、构度指标随着含水率的增加而减小,但在含水率较小或较大时,构度随含水率的变化均不大,且构度指标也可小于1,说明古土壤具有与黄土不同的结构特征。Abstract: Paleosol is an important part of loess stratum, but there are few researches on its mechanical properties.In order to study the strength and structure of paleosol, Q3paleosol in Xi'an is taken as the research object.Through the way of soil mechanics, the unconfined compressive strength test is carried out for undisturbed, remolded and saturated paleosol.The unconfined compressive strength values of the paleosol in the corresponding states are obtained, and the initial structure index is calculated.Based on the composite power exponent model, the expressions for unconfined compressive strength of the undisturbed and remolded paleosol with the water content are obtained, and the rules of the unconfined compressive strength, the disturbance sensitivity, the immersion sensitivity and the structure index of the undisturbed and remolded paleosol are analyzed.The results show that the structure of paleosol is strong, and its failure mode is still brittle when it is saturated.Under the same density, humidity and particle size, the strength of undisturbed paleosol is less than that of remolded paleosol, and its disturbance sensitivity is less than 1, which is contrary to the law of general loess.The immersion sensitivity and structure index of the paleosol decrease with the increase of water content, but when the water content is small or large, the change of structure with water content is not large, and the structure index can be less than 1, which shows that the paleosol has different structure characteristics from loess.
-
0. 引言
位于寒冷地区的水利工程，地基易受极端低温影响而发生冻胀现象，进而造成对水工建筑物结构的破坏。寒区工程超重力试验平台是利用离心力场模拟超重力场，具有克服小比尺模型的应力差异过大，避免了原位试验复杂耗时的优势，且具备热交换、干/湿环境模拟的功能，是开展寒区水利工程模型试验的重要设备。采用该设备进行寒区工程地基冻胀问题的模拟研究时，模型上方的热交换系统可在一定时间内持续为模型提供冷/热边界，利用热交换系统与模型表面间的之间的空气夹层上下表面温差进行热交换。
开展超重力场下寒区模型试验时，试验设计需了解土体与热交换系统间夹层空腔温度变化。而对空气夹层的换热过程可简化为不同尺寸的空腔及土体的耦合对流换热问题。对于空腔结构内部重力导致的对流问题，石峰等[1-2]研究了矩形腔体中Rayleigh-Benard对流结构，Ning等[3]研究了具有水平流的流体对流的周期行为问题，胡宇鹏等[4]研究了腔体内对流的稳定性问题，传统对此类的问题研究以常重力状态下为主[5-7]，对于超重力环境对低温边界下空气夹层-结构对流换热问题研究较少。由于土工模型箱内空气夹层换热过程较为复杂，易受到离心机工作状态的影响，获得超重力对低温对流换热的影响规律是设计试验面临的挑战问题。
本文基于数值模拟开展土体-夹层在超重力条件下对流换热仿真，通过分析土体界面的温度变化以及空腔内换热过程，获得模型表面温度随时间的变化规律，以及空腔内的温度场分布，获得了对换热过程对流、传导相关无量纲数变化的认识，可为开展寒区土体模拟试验提供参考。
1. 夹层传热计算模型
寒区工程超重力试验设备可简化描述为封闭空腔内自然对流传热的物理模型：模型箱内盛放一土体模型，此时模型加速度为Ng（常重力的倍数），方向由离心机主轴指向转臂。热交换板为冷表面，初始温度为T0；模型表面为光滑平面（热面），初始温度为T′0；模型表面与模型箱下表面的距离为L（夹层间隙厚度），该距离内形成一矩形密闭空腔。超重力场下的密闭空腔换热模型简化如图 1所示。
图 1中土体除上表面外均为绝热，空气夹层上表面为冷边界，土体与夹层界面通过对流及热传导进行热交换。可知空腔内流动包含两种运动：沿固壁的边界层运动以及围绕空腔中心的旋转运动，后者具有典型的Rayleigh-Benard（简称RB）对流的特征。RB对流是研究非平衡开放系统中对流运动的稳定性、时空结构、涡环（又称斑图）形成及非线性动力学特性的典型模型之一[1-4]，指在一个封闭的空腔内，其上表面温度恒定，下表面加热，由于底部流体被加热膨胀，相对于上部温度较低的流体其密度较小，下部流体上升，在上升中接触其他流体部分时消耗能量，温度变低密度变大。当底部温度继续升高，某些底部流体的温度足够大，密度足够小，能够上升到顶部而能量不至于全部耗散，从而形成温度差导致腔体内流体运动的流动现象[6]。需要指出的是，存在临界瑞利数（Ra），当空腔内局部Ra高于临界值时，才会出现RB对流[5-7]，如图 3所示，否则仅为热传导状态。
![]()
图 2 空腔界面温度在不同Ng、不同间隙时随时间变化Figure 2. Variation of temperature at interface with time under different Ng and clearances![]()
图 3 20 cm空腔50g状态换热末态时刻温度分布Figure 3. Temperature nephogram of 20 cm thick clearance at overload of 50g and last time step该物理问题与常重力场下的Rayleigh-Benard对流相似，仅重力加速度不同。在Boussinesq近似假设（认为浮力诱导的流体运动中，当温度足够小时，仅在浮力项中考虑密度变化）下，在Ng重力场下的流体力学基本控制方程如式（1）所示：
∇U=0 ，∂U/∂t+(U⋅∇)U=υ∇2U−∇P/ρ0+ρ/ρ0Ngg ，(U⋅∇)T+∂T/∂t=λ/ρ0CP∇2T 。} (1) 式中U为流体速度场；T为温度场；ρ为空气密度；g为重力加速度；υ为空气运动黏度；Cp为空气比热；t为时间；λ为导热系数。由于直接求解物理模型较为复杂，一般采用基于CFD（计算流体动力学）数值方法开展研究。
根据以上分析，模型箱的边界条件可以简化为：矩形方腔左右两侧绝热，上部为空气，下层为土体（密度为1.80 g/cm3，导热系数为2 W/m·℃，热容为4 MJ/m3·℃，热扩散系数为2.5×10-7 m2/s）。问题可简化为在Ng的超重力场下，研究矩形密闭空腔内的温度变化及温度场分布。考虑到空腔温度变化和从土体交换的热量直接相关，土体的热量与体积、质量等参数密切相关，因此建立考虑空气夹层与土体的温度分析模型，以模拟全场温度变化。同时，考虑到模型箱为长方体结构，取二维模型（侧向剖面）以表征真实三维结构，以降低计算量，提高仿真效率。计算模型的边界条件为：土体与夹层界面为换热边界，夹层上表面为-30℃，其它界面均为绝热边界，全场初始温度为10℃；即T0= -30℃；T′0=10℃，侧壁及底面均为绝热边界。根据试验状态取常重力及两种超重力状态，即Ng =1，10和50；根据土体情况取夹层厚度L =1，10和20 cm；针对两种参数组合共9种工况开展分析。
2. 夹层-土体对流换热分析
2.1 模型表面温度随时间变化
给出不同空气夹层厚度时界面温度在不同Ng值及间隙下随时间变化，如图 2所示。
图 2结果温度变化趋势与文献[8]一致。分析表明：随着Ng增大，界面降温加速，高Ng下降温速率大于低Ng工况，以20 cm为例，1g（常重力）、10g和50g状态下前2×105 s的平均降温速率分别为13.5，15，16.5℃/105 s；随着换热时间变长，降温速度随之降低，换热过程减速；降温速率与Ng并非线性关系：1，10，50g状态下200000~400000 s的平均降温速率分别为3.5，3，2℃/105 s；原因在于随着土体温度变低，土体界面与空气温差变小，导致温差驱动的对流换热的热量变小，导致降温过程随之减缓。
其次，在超重力状态下，随着夹层间隙减小，降温速率提升，小间隙下降温速率大于大间隙：以1 cm为例，1g，10g和50g状态下前2×105 s平均降温速率分别为12.5，15.5，17.5℃/105 s。空气与土体界面的换热热流率随间隙增大而减小，随Ng上升而上升：如10cm间隙时，初始阶段1g，10g和50g状态下，提取界面换热热流分别为78，165，325 W/m2；1 cm间隙时，10g和50g状态下界面热流率分别为520，630 W/m2，大于间隙较大时。界面对流换热系数变化趋势也一致。
需要指出的是，由于模型箱总容积一定，当夹层厚度较大时（20 cm），土体质量较小（总热容小），夹层厚度较小时，土体质量大（总热容大），导致通过对流及热传导将土体温度降低所需热量不同。
2.2 温度场分布特性
提取典型时刻全场温度分布及空腔流动速度，取具有代表性的20 cm空腔厚度和50g工况为例，如图 3所示，图 4给出了典型状态空腔内涡环的状态。
![]()
图 4 换热末态10 cm空腔50g状态涡环Figure 4. Vortex ring of 10 cm-clearance and overload 50g at last time step结合数值仿真结果及图 3，4可以看出：①20 cm间隙时，空腔内呈现典型的RB流动特征，形成两个对称的涡环，涡环外部流速较快，内部流速越来越小；出现涡环，超重力使得流动速度变快，1g，10g和50g下涡环中最大速度分别为0.028，0.65，1.02 m/s；在靠近壁面角落部位，会形成小的涡环，随着空腔流场的发展而变化。②10 cm间隙时，不同Ng下均会出现2对涡环；靠近侧边壁面处涡环较细长，中部涡环长宽比趋于1。③1 cm间隙时，1g状态有明显涡环成对出现；10g和50g状态下无明显的涡环发展过程，仅在靠近壁面交角附近出现较小涡环。④局部涡环会随着换热过程而发展。例如，在1 cm工况时，在初始阶段，先在靠近侧边壁面附近形成小的涡环，随之向空腔中部扩展；等到一定程度，涡环发展稳定；随着换热过程继续，涡环消失，仅在壁面附近形成较小的涡流，空腔气体温度随厚度方向变化明显，体现出传热为主的特征。涡环的稳定分布与流场速度、压力及界面换热系数、温度等直接相关。
2.3 空气夹层换热规律
常用Nu数（努塞尔数）和Ra数（瑞利数）描述换热过程中的参数变化[8]，Nu数表征换热系统从流场中吸热的强弱，Ra数表征流场中自然对流的强弱。温度在空气和土体界面上并非完全均匀分布，因此界面上同一时刻对流换热系数、Nu数等并非均匀分布，计算时取界面中点数据。给出土体与空气界面上Nu数随时间的变化如图 5所示。
不同间隙及Ng下，Nu数均随时间降低，在夹层流场启动阶段，Nu数变化较大，表示初始状态流场对流换热较为显著，后续慢慢稳定降低；说明起始阶段（温差最大时刻）换热系统工作效率较高。
3. 结论
通过建立超重力下模型箱土体夹层换热模型，分析了不同厚度及Ng对换热过程的影响规律，研究获得以下结论：
（1）随着超重力Ng值增大，界面降温加速，高Ng值下降温速率大于低Ng值；随着夹层间隙厚度减小，降温速率提升，小间隙下降温速率大于大间隙。
（2）空气与土体界面的换热热量在一定时刻内随间隙增大而减小，随Ng上升而上升，与夹层空腔中涡环发展有关。
（3）超重力状态会使得传热规律发生变化。涡环的稳定分布与流场速度、压力及界面换热系数、温度等直接相关；20 cm间隙时，空腔内呈现典型的RB流动特征，形成两个对称的涡环，超重力使得流动速度变快；10 cm间隙时，不同Ng值下均会出现2对涡环；靠近侧边壁面处涡环较细长，中部涡环长细比小；1 cm间隙时，1g状态有明显涡环成对出现；高Ng值状态下无明显的涡环发展过程，热量传递主要以传热为主。
（4）不同间隙及Ng值下，Nu数均随时间降低，在流场启动阶段，Nu数变化较大；在不同间隙状态下，随着超重力Ng值变大，界面对流换热降低。
-
![]()
图 2 不同含水率Q3古土壤的应力-应变曲线
Figure 2. Stress-strain curves of Q3paleosol with different water contents
表 1 土的物理性质指标
Table 1 Physical properties of soils
初始含水率w0/% 干密度ρd/(g·cm-3) 初始孔隙比e0 塑限wp/% 液限wL/% 17.39 1.64 0.642 22.8 36.8 
下载: 导出CSV
-
[1] 叶万军, 吴云涛, 陈明, 等. 大断面古土壤隧道围岩压力分布规律及支护结构受力特征分析——以银西高铁早胜3号隧道为例[J]. 隧道建设, 2019, 39(3): 355-361. https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD201903004.htm YE Wan-jun, WU Yun-tao, CHEN Ming, et al. Surrounding rock pressure distribution law and supporting structure stress characteristics of large cross-section paleosol tunnel: a case study of Zaosheng No.3 Tunnel on Yinchuan-Xi'an Highspeed Railway[J]. Tunnel Construction, 2019, 39(3): 355-361. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD201903004.htm
[2] 张吉宏. 古土壤地层盾构施工引起的地表沉降分析[J]. 防灾减灾工程学报, 2013, 33(1): 78-83. https://www.cnki.com.cn/Article/CJFDTOTAL-DZXK201301016.htm ZHANG Jihong. Research on ground settlement of shield tunneling in paleosol strata[J]. Jouranal of disaster Prevention and Mitigation Engineering, 2013, 33(1): 78-83. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DZXK201301016.htm
[3] 曹春山, 吴树仁, 潘懋, 等. 古土壤力学特性及其对黄土滑坡的意义[J]. 水文地质工程地质, 2016, 43(5): 127-132. https://www.cnki.com.cn/Article/CJFDTOTAL-SWDG201605019.htm CAO Chun-shan, WU Shu-ren, PAN Mao, et al. Mechanics characteristics of paleosol and its implication to loess landslide[J]. Hydrogeology and Engineering Geology. 2016.43(5): 127-132. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SWDG201605019.htm
[4] 梁莲姬, 孙有斌, CHRISTIAAN J, 等. 黄土中的碳酸盐矿物特征与化学风化[J]. 第四纪研究, 2014, 34(3): 645-653. doi: 10.3969/j.issn.1001-7410.2014.03.18 LIANG Lian-ji, SUN You-bin, CHRISTIAAN J, et al. Characteristics of carbonate minerals in loess and its implication for chemical weathering[J]. Quaternary Sciences, 2014, 34(3): 645-653. (in Chinese) doi: 10.3969/j.issn.1001-7410.2014.03.18
[5] 赵景波. 西安蓝田杨家湾黄土中第一层古土壤(S1)元素含量与环境[J]. 第四纪研究. 2011.31(3): 514-521. doi: 10.3969/j.issn.1001-7410.2011.03.15 ZHAO Jing-bo. Element content and environment of the first layer of paleosol(S1)in yangjiaban loess, lantian, Xi'an[J]. Quaternary Research, 2011, 31(3): 514-521. (in Chinese) doi: 10.3969/j.issn.1001-7410.2011.03.15
[6] 赵景波, 罗小庆, 刘瑞, 等. 关中平原黄土中第1层古土壤发育时的土壤水分研究[J]. 地质学报, 2015, 89(12): 2389-2399. doi: 10.3969/j.issn.0001-5717.2015.12.014 ZHAO Jing-bo, LUO Xiao-qing, LIU Rui, et al. Research in soil moisture during development of S1 paleosol in Guanzhou plain[J]. Acta Geologica Sinica. 2015, 89(12): 2389-2399. (in Chinese) doi: 10.3969/j.issn.0001-5717.2015.12.014
[7] 于永堂, 张继文, 郑建国, 等. 驻波比法测定黄土含水率的标定试验研究[J]. 岩石力学与工程学报, 2015, 34(7): 1462-1469. YU Yong-tang, ZHANG Ji-wen, ZHENG Jian-guo, et al. Calibration of standing-wave ratio sensors for determining water content of loess[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(7): 1462-1469. (in Chinese)
[8] 史贵君, 胡林, 林涛, 等. 关中地区黄土边坡饱和渗透系数变异性和各向异性研究[J]. 灾害学, 2019, 34(增刊1): 213-219. SHI Gui-jun. HU Lin. LIN Tao.et al. Variability and anisotropy of saturated permeability coefficient of loessslopeinGuanzhongRegion[J]. Journalof Catastrophology. 2019.34(S1): 213-219. (in Chinese)
[9] 邓军涛, 张艳, 王娟娟. 黄土古土壤的抗剪强度特性[J]. 水土保持通报, 2015, 35(5): 319-322. https://www.cnki.com.cn/Article/CJFDTOTAL-STTB201505059.htm DENG Jun-tao, ZHANG Yan, WANG Juan-juan. Shear strength characteristics of loess paleosol[J]. Bulletin of Soil and Water Conservation, 2015, 35(5): 319-322. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-STTB201505059.htm
[10] 土工试验方法标准：GB/T 50123—2019[S]. 2019. [11] 王丽琴. 黄土的结构性与湿载变形特性及其评价方法研究[D]. 西安: 西安理工大学, 2017. WANG Li-qin. Study on the Structural and Wet-Load Deformation Characteristics of Loess and Its Evaluation Methods[D]. Xi'an: Xi'an University of Technology, 2017. (in Chinese)
[12] 赵景波. 黄土地层中的CaCO3与环境[J]. 沉积学报, 1993, 11(1): 136-142. https://www.cnki.com.cn/Article/CJFDTOTAL-CJXB199301015.htm ZHAO jing-bo. CaCO3and environment in loess strata[J]. Acta Sedimentaria Sinica, 1993, 11(1): 136-142. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-CJXB199301015.htm
[13] 邵生俊, 郑文, 王正泓, 等. 黄土的构度指标及其试验确定方法[J]. 岩土力学, 2010, 31(1): 15-19. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201001004.htm SHAO Sheng-jun, ZHENG Wen, WANG Zheng-hong et al. Loess structure index and test determination method[J]. Rock and Soil Mechanics, 2010, 31(1): 15-19. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201001004.htm
[14] 王丽琴, 鹿忠刚, 邵生俊. 岩土体复合幂-指数非线性模型[J]. 岩石力学与工程学报, 2017, 36(5): 1269-1278. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201705022.htm WANG Li-qin, LU Zhong-gang, SHAO Sheng-jun. A composite power exponential nonlinear model of rock and soil[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(5): 1269-127. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201705022.htm
计量
- 文章访问数: 174
- HTML全文浏览量: 40
- PDF下载量: 190
