• 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

半无限空间饱和土一维热固结精确积分解

钮家军, 王秀凯, 凌道盛, 巩师林, 石吉森

钮家军, 王秀凯, 凌道盛, 巩师林, 石吉森. 半无限空间饱和土一维热固结精确积分解[J]. 岩土工程学报, 2021, 43(8): 1426-1433. DOI: 10.11779/CJGE202108007
引用本文: 钮家军, 王秀凯, 凌道盛, 巩师林, 石吉森. 半无限空间饱和土一维热固结精确积分解[J]. 岩土工程学报, 2021, 43(8): 1426-1433. DOI: 10.11779/CJGE202108007
NIU Jia-jun, WANG Xiu-kai, LING Dao-sheng, GONG Shi-lin, SHI Ji-sen. Exact integral solutions for one-dimensional thermal consolidation of semi-infinite saturated soils[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(8): 1426-1433. DOI: 10.11779/CJGE202108007
Citation: NIU Jia-jun, WANG Xiu-kai, LING Dao-sheng, GONG Shi-lin, SHI Ji-sen. Exact integral solutions for one-dimensional thermal consolidation of semi-infinite saturated soils[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(8): 1426-1433. DOI: 10.11779/CJGE202108007

半无限空间饱和土一维热固结精确积分解  English Version

基金项目: 

国家重点研发计划项目 2016YFC0800200

国家自然科学基金项目 51578502

详细信息
    作者简介:

    钮家军(1993— ),男,博士研究生,主要从事多孔介质热–水–力耦合现象的研究。E-mail:niu0918@zju.edu.cn

    通讯作者:

    凌道盛, E-mail:dsling@zju.edu.cn

  • 中图分类号: TU431

Exact integral solutions for one-dimensional thermal consolidation of semi-infinite saturated soils

  • 摘要: 基于Coussy提出的饱和多孔介质热–水–力耦合理论,给出考虑热渗效应和等温热流效应的饱和土热固结耦合控制方程,然后通过正弦和余弦变换方法,结合Hamilton-Cayley定理,给出了两类任意非齐次边界条件下,半无限空间饱和土一维热固结精确积分解的求解方法,避免了Laplace变换处理复杂边界问题时的数值反演。最后,通过算例验证了本文解析解的正确性,并分析了半无限空间表面透水或不透水时,饱和土体在季节性温度荷载作用下的温度、超静孔压和位移响应特性。结果表明:温度、超静孔压和位移受季节性温度荷载影响而呈现周期性波动,表面不透水情形的土体超静孔压和位移响应明显高于相应的表面透水情形的土体响应。
    Abstract: Based on the thermo-hydro-mechanical coupled theory for saturated porous media proposed by Coussy, two coupled governing equations for thermal consolidation of saturated soils are presented, taking thermo-osmosis and thermo- filtration into account. The sine and cosine transforms as well as the Hamilton-Cayley law are used to obtain exact integral solutions for one-dimensional thermal consolidation of semi-infinite saturated soils with two types of arbitrary nonhomogeneous boundary conditions. It helps avoid the errors of numerical inverse Laplace transform when dealing with complex boundary conditions. Finally, case studies are provided to validate the proposed analytical solutions and to investigate the thermal responses of saturated soils with seasonal thermal loads on the permeable or impermeable surface. The results show that temperature change, excess pore water pressure and displacement all fluctuate periodically due to seasonal thermal loads. The responses of the excess pore water pressure and displacement of the saturated soils with an impermeable surface are much higher than those with a permeable surface.
  • 岩体工程多赋存于水环境中,并常处于饱和状态。环境中的水进入岩体孔隙后易发生水岩相互作用,将改变岩石的矿物组成和微细观结构,引起强度、变形等物理力学性能的劣化。为保证水环境中岩体工程的稳定性,研究饱和岩石的水软化特性和机制十分必要。

    目前,许多学者对多种饱和岩石的水软化特性进行了研究。普遍认为,岩石饱和后强度降低,弹性模量减小,泊松比增加,摩擦系数减小。如周翠英等[1]发现泥岩饱和后,其强度下降了约90%。Mann等[2]对3种饱和砂岩进行压缩试验,发现弹性模量减小了6%~19%。熊德国等[3]得到砂岩、砂质泥岩和泥岩的抗压强度水软化系数(饱和岩石抗压强度与干燥岩石抗压强度之比)范围为0.58~0.94。Zhu等[4]得到大理岩饱和后,单轴抗压强度减小了29.18%,泊松比增加13.90%。黄弈茗[5]通过剪切试验得到饱和灰岩的摩擦系数降低了25%。Chen等[6]、Erguler等[7]和Vasarhelyi[8]发现凝灰岩饱和后强度降低了2%~88%。工程上常将单轴抗压强度小于20 MPa的岩石称为软岩,20~40 MPa的称为中硬岩,大于40 MPa的称为硬岩。总结来看,目前软岩的水软化特性研究成果丰富,硬岩的水软化特性研究存在不足。

    为合理解释饱和岩石力学性能的劣化现象,一些学者提出了诸如黏土矿物软化、矿物颗粒摩擦力下降、孔隙水压力等水软化机制[9-13]。然而,受限于试验手段和分析方法,饱和岩石力学性能劣化的物理力学机制尚未达成共识,试验成果多,理论分析成果少。Zhu等[4]综合已有试验成果发现,硬岩的水软化机制可能包括:机制一,黏土矿物软化;机制二,孔隙水压力作用;机制三,摩擦弱化效应;机制四,应力腐蚀。其中,机制一、四分别存在于含黏土矿物和石英矿物的岩石中,且机制一已被公认。机制二、三(统称水岩作用机制,下同)普遍存在于各类饱和岩石中,受到了广泛关注。值得提出的是,上述水软化机制常常同时存在,并相互影响。因此,排除其他针对某一个或几个特定水软化机制,定性和定量描述其存在条件和作用程度成为亟待解决的问题。

    声发射(acoustic emission,AE)是材料加载过程中局部应变能释放所产生的瞬态弹性波,因其可以反映材料内部裂隙发育和破坏特性而被广泛应用于岩石力学中[14-15]。声发射信号的分析方法主要包括参数分析和波形分析。目前一些学者利用声发射特征参数,如声发射事件率、声发射振铃计数、能量等,研究含水岩石的破坏特性[16-18]。虽然声发射参数分析能反映饱和岩石破坏的一些释能特征,但将复杂难懂的波形信号抽象成一系列的参数来加以描述,难免会遗失波形所携带的一些重要的信息。近年来,一些学者开始研究声发射波形信号的频谱特征[14, 19-20]。关于声发射频谱在含水岩石破坏特性研究中的运用,少量学者收获了十分有限的成果。例如,Chugh等[21]、Read等[22]分析了干燥和含水岩石的声发射波形信号的主频差异。张艳博等[23]发现饱和砂岩的声发射信号主频主要分布在35~60 kHz,且主频范围比干燥砂岩集中。然而,基于声发射波形,特别是波形信号的主频统计特征,分析饱和岩石的微观破坏特性进而研究水岩作用机制的成果鲜有报道。

    为分析饱和岩石水岩作用机制的存在条件和作用程度,本文选取不含黏土和硅质矿物的典型硬岩,以排除机制一、四的影响,开展硬岩在干燥、饱和状态下的单轴压缩声发射试验,分析压缩强度、弹性模量、破裂形态等力学软化特性,探讨饱和硬岩破坏的声发射波形信号主频统计特征。结合力学软化特性和声发射主频统计特征,研究硬岩水岩作用机制的存在条件和作用程度。研究成果可丰富饱和硬岩力学软化特性的研究内容,并为量化描述水岩作用机制的研究提供新的视角和途径。

    试验用岩石包括3种:①宝兴大理岩,取样自四川省雅安市宝兴县,②锦屏大理岩,取样自四川省凉山彝族自治州盐源县锦屏一级水电站库岸边坡,③灰岩,取样自四川省凉山彝族自治州雷波县马湖边坡。依据国际岩石力学学会(ISRM)建议方法将岩芯加工成ϕ50 mm×H100 mm的标准圆柱体试样。如图 1所示,采用X射线单晶衍射仪测定3种岩石的矿物成分和各矿物占比,发现宝兴大理岩和灰岩均由单一矿物方解石组成,锦屏大理岩由方解石(12%)、白云石(87%)、氯化钾(0.5%)、氧化铁(0.5%)组成。说明3种岩石均不含黏土矿物和石英矿物,黏土矿物软化和应力腐蚀作用可以忽略。电镜扫描(SEM)图像显示3种岩石的矿物胶结形态、孔隙数量存在较大差异。利用核磁共振岩芯分析系统(NMR)测定试样的孔隙度。试样的物理特性参数如表 1所示。

    图  1  各种硬岩XRD测试结果和SEM图
    Figure  1.  XRD and SEM results of different rocks
    表  1  试样基本物理性质参数
    Table  1.  Basic physical parameters of samples
    岩石类型 天然密度/(kg·m-3) 孔隙度/% 纵波波速/(m·s-1) 矿物成分
    宝兴大理岩 2.69 (±0.04) 6.16 (±0.25) 4262 (±30.73) 方解石
    锦屏大理岩 2.98 (±0.17) 2.68 (±0.13) 4689 (±29.32) 方解石、白云石、氯化钾、氧化铁
    灰岩 2.63 (±0.11) 4.80 (±0.09) 3826 (±18.67) 方解石
    注:括号内数值为试样样本参数的标准差。
    下载: 导出CSV 
    | 显示表格

    为获得干燥、饱和试样,将加工好的试样放入105℃的恒温烘箱内烘48 h,使其达到完全干燥状态。冷却至室温后取一半试样置于干燥箱内保存。采用强制饱和法,将剩余试样置于100 kPa的负压真空抽气机中强制饱和48 h,并静置4 h。随机选取3种岩石的干燥、饱和试样各2个,采用干馏法测定水饱和度,发现其值分别接近0和1,认为获得的干燥、饱和试样满足试验要求。

    试验采用MTS815 Flex Test GT岩石力学试验系统和PCI-2声发射采集系统进行。每个试样布置8个Micro30高精度传感器,可以有效采集试验过程的声发射信号[4, 13]。声发射前置放大器门槛值(增益)为40 dB,采样门槛值为30 dB。试样与传感器间涂抹耦合剂,以保证信号数据的有效采集。试验设备与传感器布置如图 2所示。试验采用LVDT进行轴向位移控制,加载速率为0.05 mm/min。每种岩石包括干燥、饱和试验组各1组,每组4个试样,共计24个。试验过程中,保证力学加载系统与声发射采集系统的时间同步。

    图  2  试验装置和声发射传感器布置图
    Figure  2.  Setup and layout of AE sensors

    3种岩石典型试样的应力应变曲线如图 3所示。可以看出,试样在单轴压缩荷载作用下的应力应变曲线具有明显的阶段性,即加载初期的非线性段(压密段)、线弹性段、峰值应力前的塑性变形段和峰后段。饱和试样的峰值应变(峰值应力对应的应变)比干燥试样的大,表现出应变软化特性。相比较而言,锦屏大理岩力学性能最强,灰岩次之,宝兴大理岩最弱。

    图  3  典型试样应力应变曲线
    Figure  3.  Stress-strain curves of typical rock samples

    干燥、饱和状态下3种岩石的单轴抗压强度σ、弹性模量E及软化系数f图 45所示。其中,软化系数f为饱和试样组的力学参数均值与干燥试样组的力学参数均值的比值,fσfE分别表示单轴抗压强度软化系数和弹性模量软化系数。

    图  4  干燥、饱和状态岩石单轴抗压强度和软化系数
    Figure  4.  Uniaxial compressive strengths and softening coefficients of rocks under dry and saturated conditions
    图  5  干燥、饱和状态岩石弹性模量和软化系数
    Figure  5.  Elastic moduli and softening coefficients of rocks under dry and saturated conditions

    可以看出,3种岩石饱和后,力学性能均出现明显劣化,表现为单轴抗压强度和弹性模量减小。宝兴大理岩、锦屏大理岩、灰岩的单轴抗压强度软化系数分别为0.74,0.94,0.76,弹性模量软化系数分别为0.72,0.92,0.79。相比较而言,宝兴大理岩水软化性能最强,灰岩次之,锦屏大理岩最弱。

    干燥、饱和状态下典型试样的宏观破坏形态与声发射事件定位见图 6,图中声发射事件的颜色对应不同的加载时间。可以看出,声发射事件多发生于试样的宏观破坏面或裂纹处,宏观破裂形态与声发射事件定位能较好吻合。相比较而言,锦屏大理岩破坏时释放的声发射信号比宝兴大理岩和灰岩的多。在压缩荷载作用下,干燥岩石多发生剪切破坏,试样可见明显的单斜面或共轭斜面剪切破坏面。岩石饱和后,平行于试样轴向的张拉裂纹和张拉破裂面增多。

    图  6  岩石宏观破裂形态和声发射事件定位
    Figure  6.  Failure patterns of rocks and location of acoustic emission events

    (1)数据处理

    试验采集了加载全过程岩石释放的声发射波形信号。本文采用声发射信号频谱统计分析方法,即将所有通道采集的海量信号数据进行统计分析。首先,通过MATLAB编程,采用快速傅立叶变换(FFT)获得声发射波形信号的主频参数(图 7)。其次,对获得的波形信号主频进行区段划分。室内岩石的声发射频率位于10~500 kHz范围内[24],且一般低于450 kHz。为便于统计分析,参照李林芮[25]、邓建辉等[26]的研究,采用的主频区段划分方法为:主频0~450 kHz,按每10 kHz划分一个主频段,共45个主频段;主频451~500 kHz,划分为一个主频段,合计46个主频段。

    图  7  典型声发射信号的时域图和频谱图
    Figure  7.  Typical acoustic emission waveforms and their spectra

    (2)声发射主频特征

    干燥、饱和状态下3种岩石典型试样的主频占比统计见图 8。可以发现,在压缩荷载作用下,干燥、饱和岩石的声发射波形主频均集中在2个主频带范围内,称之为高主频带和低主频带。作为邓建辉等[26]关于岩石声发射双主频特性的后续研究,说明双主频现象是岩石破坏的一种固有特性,与岩石类型和含水状态无关。岩石饱和后,高、低主频带范围发生变化。以宝兴大理岩为例,干燥宝兴大理岩的高、低主频带分别集中在230~320 kHz(第24~32主频段)和11~80 kHz(第2~8主频段),饱和宝兴大理岩高、低主频带分别集中在230~270 kHz(第24~27主频段)和0~50 kHz(第1~5主频段)。

    图  8  干燥、饱和状态下岩石的主频占比统计
    Figure  8.  Proportions of dominant frequency of rock samples under dry and saturated conditions

    将占比超过1%的主频段计入样本范围,统计得到干燥、饱和岩石的声发射高、低主频数量占比和能量占比,如图 9所示。干燥、饱和试验组试样主频数量占比均值和能量占比均值如表 2所示。其中能量以声发射信号幅值的平方表征。干燥岩石破坏的声发射以携带高能量的高主频带信号为主,最小高主频信号数量占比均值为69.51%,最小高主频信号能量占比均值为70.70%。3种硬岩饱和后,高主频带波形信号明显减少,低主频带波形信号显著增多。低主频带信号占据主导地位,且携带能量更多。最小低主频信号数量占比均值为58.77%,最小低主频信号能量占比均值为65.01%。

    图  9  岩石高、低主频信号数量占比及能量占比
    Figure  9.  Percentages of high and low dominant frequency bands and corresponding energy of different rocks
    表  2  干燥、饱和岩石的声发射主频参数均值统计
    Table  2.  Average values of dominant frequency parameters of different rocks under dry and saturated conditions
    岩石类型 状态 N/104 低频带 高频带
    PN/% PE/% PN/% PE/%
    宝兴大理岩 干燥 6.89 25.55 29.14 70.29 70.70
    饱和 5.56 73.13 73.16 23.59 26.84
    锦屏大理岩 干燥 2.42 28.63 25.15 70.84 74.35
    饱和 1.31 58.77 65.01 35.18 29.11
    灰岩 干燥 5.79 27.93 21.39 69.51 78.24
    饱和 3.48 68.65 78.52 29.73 20.60
    注:表中NPNPE分别表示波形数、不同主频带信号数量占比、不同主频带信号能量占比。
    下载: 导出CSV 
    | 显示表格

    硬岩的矿物成分决定水软化机制的存在条件。宝兴大理岩、锦屏大理岩和灰岩不含黏土矿物和硅质矿物,可能存在的水软化机制包括孔隙水压力(机制二)和摩擦弱化效应(机制三)。如表 2图 9所示,上述3类硬岩饱和后,低主频信号数量及携带能量显著增加,并占绝对优势。李林芮[25]、邓建辉等[26]提出并验证,岩石的高主频波形信号对应微观剪切破坏,低主频波形信号对应微观拉伸破坏。换言之,本文研究的硬岩饱和后,发生的微观拉伸破坏明显增多,并占据主导地位。而3类岩石只可能存在孔隙水压力和摩擦弱化效应水软化机制,一种比较直观的推论是两种水软化机制中的一种或全部导致微观拉伸破坏数量的增加。

    为论证上述推论,采用ANSYS数值软件对有(无)孔隙压力的单一孔隙周边的应力分布进行模拟分析。将试样孔隙简化为半径为1 mm的圆形孔洞,轴向压力为20 MPa。通过迭代计算获得孔洞周边的应力分布,如图 10所示。对比分析可以得出,有孔隙压力时,孔隙周边的拉应力普遍更大,最大拉应力由19 MPa(无孔隙压力时)增加至21.7 MPa,即孔隙压力导致拉应力增加。因此可以推断,饱和岩石的孔隙水压力作用机制将导致孔隙周边拉应力增加,产生更多微观拉伸破坏。这一推断也从图 6中饱和硬岩的破坏形态得到证实,即饱和硬岩出现更多张拉裂纹和张拉破裂面。而微观拉伸破坏对应声发射低主频波形信号,所以声发射低主频信号数量增加。总结来说,饱和岩石孔隙水压力作用的强弱与声发射低主频信号的多少存在对应关系。

    图  10  单一孔隙周边应力的数值模拟
    Figure  10.  Numerical simulation of stress distribution around a single pore

    3类硬岩饱和后低主频信号占比的增量与单轴抗压强度软化系数的关系如图 11所示。可以得到,低主频信号占比增量与强度软化系数间存在负相关关系,即低主频信号占比增量越大,强度降低越明显。这佐证了孔隙水压力是上述3种硬岩强度降低的主导因素。

    图  11  低频占比增量与单轴抗压强度软化系数关系
    Figure  11.  Relationship between increment of percentage of low dominant frequency and softening coefficient of uniaxial compressive strength

    相比较而言,宝兴大理岩和灰岩的低主频信号增量相差不多,强度软化系数也接近,而锦屏大理岩则表现出巨大差异。笔者认为这与岩石的矿物成分和孔隙结构密切相关。宝兴大理岩和灰岩的矿物成分均为方解石,孔隙度分别为6.16%和4.80%,两者矿物成分和孔隙结构相似,因而软化系数接近,分别为0.74和0.76。锦屏大理岩由多种矿物组成,孔隙度为2.68%。由图 1的SEM图像可以看出,矿物颗粒胶结紧密,未见明显孔隙。可以证实,锦屏大理岩结构致密,未提供孔隙水压力作用的良好条件,发生的微观拉伸破坏比宝兴大理岩和灰岩少,因而声发射低主频信号占比增量更小。综上分析,不含黏土、硅质矿物岩石的矿物成分和孔隙结构决定孔隙水压力作用的强弱。

    大多数岩体工程,如水电站库岸边坡、部分地下工程等,常处于水环境中,易引起岩体物理力学性能劣化,威胁工程的长期稳定和安全。在岩体工程各类分级评估中,包括工程岩体分级、岩体参数评估等,饱和单轴抗压强度作为强制性指标被广泛应用[27]。由于岩石饱和后,强度降低,采用饱和单轴抗压强度进行评估提高了工程的安全系数。但“全盘”要求使用饱和单轴抗压强度使工程稳定性评估有时过于保守,不利于建设的高效和经济。

    不同类型饱和岩石的矿物成分、微结构、孔隙度和连通性等差异巨大,这决定了各种水软化机制的存在条件和作用程度。确定岩石矿物成分是水软化机制分析的首要工作。如前言所述,不含黏土矿物和硅质矿物的岩石可能存在孔隙水压力作用和摩擦弱化效应。本文分析得出,饱和碳酸钙类岩石的水软化机制以孔隙水压力为主导,产生更多的微观拉伸破坏,且孔隙水压力作用强弱受到孔隙度的影响。由于岩石抗压剪不抗拉伸[28],岩体工程分级评估时,应以提高安全系数为原则,即由大孔隙度、碳酸钙类矿物组成的岩体工程处于水环境中时,应严格按饱和单轴抗压强度进行设计评估,以保证工程安全。

    声发射低主频信号增量与强度软化系数的负相关关系可应用于水环境中岩体工程的设计评估。在压缩荷载作用下,饱和岩石的声发射低主频信号增量与孔隙水压力作用强弱相对应,进而可预测发生微观拉伸破坏的多少。普遍认为,微观拉伸破坏与岩石强度间存在相关性。如图 12所示,笔者进行不同加载速率下宝兴大理岩干燥试样的单轴压缩试验,获得低主频信号(拉伸破坏)占比与单轴抗压强度的拟合曲线。发现两者间存在负相关性,即低主频信号越多,岩石强度越低。这与岩石抗压剪不抗拉伸的事实相符,由于岩石内部发生的微观拉伸破坏(低主频信号)多,岩石强度低。说明对应拉伸破坏的低主频信号可用于表征岩石强度变化。本文试验用3种岩石的强度软化系数随声发射低主频信号增量的增加而减小,说明声发射低主频信号可以表征碳酸钙类饱和岩石强度的劣化程度。声发射低主频信号可作为一种新的特征参数应用于含水岩体工程的设计和安全评估。

    图  12  宝兴大理岩强度与低主频信号占比间关系
    Figure  12.  Relationship between uniaxial compressive strength of marble and percentage of low dominant frequency

    对3类硬岩(宝兴大理岩、锦屏大理岩、灰岩)饱和后的力学软化特性和声发射主频统计特征进行了分析,探究了水岩作用机制之孔隙水压力的存在条件和作用程度,取得以下3点结论。

    (1)硬岩饱和后力学软化特性明显,单轴抗压强度和弹性模量减小,出现更多张拉裂纹和张拉破坏面。

    (2)硬岩破坏释放的声发射信号存在双主频现象,且与岩石类型和含水状态无关。硬岩饱和后,声发射高主频带信号明显减少,低主频带信号显著增加。

    (3)水岩作用机制之孔隙水压力是不含黏土、硅质矿物硬岩性能劣化的主导因素,其作用强弱与声发射低主频信号多少存在对应关系,且取决于矿物成分和孔隙结构。

  • 图  1   半无限空间饱和土一维热固结计算模型

    Figure  1.   Model for one-dimensional thermal consolidation of semi-infinite saturated soils

    图  2   表面温度突增时不同时刻温度分布图

    Figure  2.   Distribution of temperature change along depth at different moments subjected to step temperature increase

    图  3   表面温度突增时不同时刻超静孔压分布图

    Figure  3.   Distribution of excess pore pressure along depth at different moments subjected to step temperature increase

    图  4   边界条件(A)不同深度处温度变化过程

    Figure  4.   Evolution of temperature change with time at different depths under boundary condition A

    图  5   边界条件(A)不同深度处超静孔压变化过程

    Figure  5.   Evolution of excess pore pressure with time at different depths under boundary condition A

    图  6   边界条件(A)表面位移变化过程

    Figure  6.   Evolution of surface displacement of soils with time under boundary condition A

    图  7   边界条件(B)不同深度处温度变化过程

    Figure  7.   Evolution of temperature change with time at different depths under boundary condition B

    图  8   边界条件(B)不同深度处超静孔压变化过程

    Figure  8.   Evolution of excess pore pressure with time at different depths under boundary condition B

    图  9   边界条件(B)表面位移变化过程

    Figure  9.   Evolution of surface displacement of soils with time under boundary condition B

    图  10   边界条件(B)不同扩散数表面温度随时间变化

    Figure  10.   Evolution of temperature of surface with time with different coefficients under boundary condition B

    图  11   边界条件(B)不同扩散数表面超静孔压随时间变化

    Figure  11.   Evolution of excess pore pressure of surface with time with different coefficients under boundary condition B

    图  12   边界条件(B)不同扩散数表面位移随时间变化

    Figure  12.   Evolution of displacement of surface with time with different coefficients under boundary condition B

    图  13   边界条件(B)不同扩散数超静孔压随深度变化

    Figure  13.   Distribution of excess pore pressure along depth with different coefficients under boundary condition B

    表  1   饱和土体的材料参数

    Table  1   Material parameters of saturated soils

    孔隙率n泊松比μ土体弹性模量E/MPa土颗粒体积模量Ks/GPa土颗粒密度ρs/(kg·m-3)土颗粒比热容cs/(J·kg-1·K-1)土颗粒热膨胀系数αs/K-1孔隙流体体积模量Kw/GPa孔隙流体密度ρw/(kg·m-3)孔隙流体比热容cw/(J·kg-1·K-1)孔隙流体热膨胀系数αw/K-1土体热传导系数KT/(J·s-1·m-1·K-1)热渗系数Sw/(m2·s-1·K-1)
    0.40.352026008001.5×10-55100042002×10-515×10-11
    下载: 导出CSV
  • [1]

    BOOKER J R, SAVVIDOU C. Consolidation around a spherical heat source[J]. International Journal of Solids and Structures, 1984, 20(11/12): 1079-1090.

    [2]

    BOOKER J R, SAVVIDOU C. Consolidation around a point heat source[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1985, 9(2): 173-184. doi: 10.1002/nag.1610090206

    [3]

    MCTIGUE D F. Thermoelastic response of fluid-saturated porous rock[J]. Journal of Geophysical Research, 1986, 91(B9): 9533-9542. doi: 10.1029/JB091iB09p09533

    [4]

    GIRAUD A, ROUSSET G. Thermoelastic and thermoplastic response of a porous space submitted to a decaying heat source[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1995, 19(7): 475-495. doi: 10.1002/nag.1610190703

    [5]

    SMITH D W, BOOKER J R. Green’s functions for a fully coupled thermoporoelastic material[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1993, 17(3): 139-163. doi: 10.1002/nag.1610170302

    [6]

    SMITH D W, BOOKER J R. Boundary element analysis of linear thermoelastic consolidation[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1996, 20(7): 457-488. doi: 10.1002/(SICI)1096-9853(199607)20:7<457::AID-NAG805>3.0.CO;2-H

    [7]

    BAI M, ABOUSLEIMAN Y. Thermoporoelastic coupling with application to consolidation[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1997, 21(2): 121-132. doi: 10.1002/(SICI)1096-9853(199702)21:2<121::AID-NAG861>3.0.CO;2-W

    [8]

    ZHOU Y, RAJAPAKSE R K N D, GRAHAM J. A coupled thermoporoelastic model with thermo-osmosis and thermal filtration[J]. International Journal of Solids and Structures, 1998, 35(34/35): 4659-4683.

    [9]

    ZHOU Y, RAJAPAKSE R K N D, GRAHAM J. Coupled consolidation of a porous medium with a cylindrical or a spherical cavity[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 1998, 22(6): 449-475. doi: 10.1002/(SICI)1096-9853(199806)22:6<449::AID-NAG930>3.0.CO;2-H

    [10]

    BLOND E, SCHMITT N, HILD F. Response of saturated porous media to cyclic thermal loading[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2003, 27(11): 883-904. doi: 10.1002/nag.301

    [11] 白冰. 循环温度荷载作用下饱和多孔介质热-水-力耦合响应[J]. 工程力学, 2007, 24(5): 87-92. doi: 10.3969/j.issn.1000-4750.2007.05.015

    BAI Bing. Thermo-hydro-mechanical responses of saturated porous media under cyclic thermal loading[J]. Engineering Mechanics, 2007, 24(5): 87-92. (in Chinese) doi: 10.3969/j.issn.1000-4750.2007.05.015

    [12]

    BAI B. Fluctuation responses of saturated porous media subjected to cyclic thermal loading[J]. Computers and Geotechnics, 2006, 33(8): 396-403. doi: 10.1016/j.compgeo.2006.08.005

    [13]

    BAI B. Thermal consolidation of layered porous half-space to variable thermal loading[J]. Applied Mathematics and Mechanics, 2006, 27(11): 1531-1539. doi: 10.1007/s10483-006-1111-1

    [14]

    LIU G B, XIE K H, ZHENG R Y. Model of nonlinear coupled thermo-hydro-elastodynamics response for a saturated poroelastic medium[J]. Science in China Series E: Technological Sciences, 2009, 52(8): 2373-2383. doi: 10.1007/s11431-008-0220-8

    [15]

    LIU G B, XIE K H, ZHENG R Y. Thermo-elastodynamic response of a spherical cavity in saturated poroelastic medium[J]. Applied Mathematical Modelling, 2010, 34(8): 2203-2222. doi: 10.1016/j.apm.2009.10.031

    [16]

    HE L W, JIN Z H. A local thermal nonequilibrium poroelastic theory for fluid saturated porous media[J]. Journal of Thermal Stresses, 2010, 33(8): 799-813. doi: 10.1080/01495739.2010.482358

    [17]

    YANG Y, GUERLEBECK K, SCHANZ T. Thermo-osmosis effect in saturated porous medium[J]. Transport in Porous Media, 2014, 104(2): 253-271. doi: 10.1007/s11242-014-0332-5

    [18]

    AI Z Y, YE Z, ZHAO Z, et al. Time-dependent behavior of axisymmetric thermal consolidation for multilayered transversely isotropic poroelastic material[J]. Applied Mathematical Modelling, 2018, 61: 216-236. doi: 10.1016/j.apm.2018.04.012

    [19]

    AI Z Y, YE Z, LI Y. Quasti-static response of a multilayered transversely isotropic porothermoelastic material subjected to a cylindrical heat source[J]. Computers and Geotechnics, 2019, 107: 201-213. doi: 10.1016/j.compgeo.2018.11.024

    [20]

    AI Z Y, WANG L J. Precise solution to 3D coupled thermohydromechanical problems of layered transversely isotropic saturated porous media[J]. International Journal of Geomechanics, 2018, 18(1): 4017121. doi: 10.1061/(ASCE)GM.1943-5622.0001027

    [21]

    WANG L J, ZHU B, CHEN Y M, CHEN R P, SHI X S. Precise model for predicting excess pore-water pressure of layered soils induced by thermal-mechanical loads[J]. Journal of Engineering Mechanics, 2019, 145(1): 4018114. doi: 10.1061/(ASCE)EM.1943-7889.0001544

    [22] 王路君, 艾智勇. 衰变热源作用下饱和多孔介质热固结问题的扩展精细积分法[J]. 力学学报, 2017, 49(2): 324-334. https://www.cnki.com.cn/Article/CJFDTOTAL-LXXB201702009.htm

    WANG Lu-jun, AI Zhi-yong. EPIM for thermal consolidation problems of saturated porous media subjected to a decaying heat source[J]. Chinese Journal of Theoretical and Applied Mechanics, 2017, 49(2): 324-334. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-LXXB201702009.htm

    [23] 钮家军, 凌道盛, 王秀凯, 等. 饱和单层土体一维热固结精确解[J]. 岩土工程学报, 2019, 41(9): 1715-1723. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201909018.htm

    NIU Jia-jun, LING Dao-sheng, WANG Xiu-kai, et al. Exact solutions for one-dimensional thermal consolidation of single-layer saturated soil[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(9): 1715-1723. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201909018.htm

    [24]

    COUSSY O. Poromechanics[M]. Chichester: John Wiley & Sons Ltd, 2004.

    [25]

    BAI B. Thermal response of saturated porous spherical body containing a cavity under several boundary conditions[J]. Journal of Thermal stresses, 2013, 36(11): 1217-1232. doi: 10.1080/01495739.2013.788389

    [26]

    BAI B, LI T. Irreversible consolidation problem of a saturated porothermoelastic spherical body with a spherical cavity[J]. Applied Mathematical Modelling, 2013, 37(4): 1973-1982. doi: 10.1016/j.apm.2012.05.003

    [27]

    CARNAHAN C L. Thermodynamic coupling of heat and matter flows in near-field regions of nuclear waste repositories[J]. MRS Proceedings, 1983, 26(1): 1023-1030.

    [28]

    HORSEMAN S T, MCEWEN T J. Thermal constraints on disposal of heat-emitting waste in argillaceous rocks[J]. Engineering Geology, 1996, 41(1/2/3/4): 5-16.

图(13)  /  表(1)
计量
  • 文章访问数:  253
  • HTML全文浏览量:  59
  • PDF下载量:  130
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-09-08
  • 网络出版日期:  2022-12-02
  • 刊出日期:  2021-07-31

目录

/

返回文章
返回