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

基于非傅里叶定律的近场动力学热力耦合模型及花岗岩热损伤破裂模拟

周露明, 朱珍德, 谢兴华, 吕茂淋

周露明, 朱珍德, 谢兴华, 吕茂淋. 基于非傅里叶定律的近场动力学热力耦合模型及花岗岩热损伤破裂模拟[J]. 岩土工程学报, 2024, 46(11): 2391-2400. DOI: 10.11779/CJGE20230769
引用本文: 周露明, 朱珍德, 谢兴华, 吕茂淋. 基于非傅里叶定律的近场动力学热力耦合模型及花岗岩热损伤破裂模拟[J]. 岩土工程学报, 2024, 46(11): 2391-2400. DOI: 10.11779/CJGE20230769
ZHOU Luming, ZHU Zhende, XIE Xinghua, LÜ Maolin. Non-Fourier law-based peridynamic thermo-mechanical coupling model and simulation of thermal damage and fracture in granite[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(11): 2391-2400. DOI: 10.11779/CJGE20230769
Citation: ZHOU Luming, ZHU Zhende, XIE Xinghua, LÜ Maolin. Non-Fourier law-based peridynamic thermo-mechanical coupling model and simulation of thermal damage and fracture in granite[J]. Chinese Journal of Geotechnical Engineering, 2024, 46(11): 2391-2400. DOI: 10.11779/CJGE20230769

基于非傅里叶定律的近场动力学热力耦合模型及花岗岩热损伤破裂模拟  English Version

基金项目: 

国家自然科学基金项目 51878249

国家自然科学基金项目 41831278

详细信息
    作者简介:

    周露明(1994—),男,博士研究生,主要从事岩石力学与地下能源工程方面的研究工作。E-mail: 190204030004@hhu.edu.cn

  • 中图分类号: TU45

Non-Fourier law-based peridynamic thermo-mechanical coupling model and simulation of thermal damage and fracture in granite

  • 摘要: 研究岩石的热损伤破裂特征对于地热开采等工程具有重要意义。在传统常规态基近场动力学理论框架内,通过引入双相滞后(DPL)模型,提出了基于非傅里叶热传导定律推导得到的热力耦合模型。通过平板瞬态热传导问题及Lac du Bonnet(LdB)花岗岩的热损伤破裂试验对模型进行了验证。分析了温度梯度弛豫时间和热流弛豫时间对岩石热损伤破裂特征的影响。结果表明:近场动力学热力耦合模型模拟结果很好地反映了LdB花岗岩的热损伤破裂特征及温度分布的不连续性;温度梯度弛豫时间对热传导起促进作用,试件热损伤破裂程度随温度梯度弛豫时间的增加而增大,随热流弛豫时间的增加而减小。研究为深入理解岩石的热损伤破裂行为、优化地热能开采工程提供了有益的探索。
    Abstract: It is of great significance to study the thermal damage and fracture characteristics of rocks for deep rock projects, such as geothermal exploitation. Within the framework of the traditional classical ordinary state-based peridynamic theory, a thermo-mechanical coupling model based on the non-Fourier heat conduction law is proposed by introducing a dual-phase-lag model. The model is validated through the transient heat conduction problems in the plate and thermal damage and fracture tests on LdB granite. It is found that the simulated results accurately reflect the thermal damage and fracture characteristics, as well as the discontinuity in temperature distribution, of LdB granite. Furthermore, through numerical analysis, it is observed that the temperature gradient relaxation time promotes heat conduction, leading to an increase in the degree of the thermal damage and fracture with an increase in the temperature gradient relaxation time, while it decreases with an increase in the heat flux relaxation time. This research provides valuable insights for a better understanding of the thermal damage and fracture behavior of rocks and offers beneficial explorations for optimizing the geothermal energy extraction projects.
  • 随着中国经济社会的快速发展,城市规模和城市化水平不断提高,人口过度集中对交通与环境的压力不断增大,城市地铁因其高效便捷,对地面占用小的特点,逐渐成为一剂解决城市交通拥堵与环境污染的“灵丹妙药”。在这种背景下,城市隧道的施工技术便显得尤为重要。盾构法施工对环境影响小,不受地形、地貌、江河水域等条件限制,以及施工安全、快速等诸多优点[1],成为城市隧道施工的主要方法之一。但是,由于盾构机造价高昂,以现阶段中国工程造价水准,还不具备一项工程完全折旧的能力。中国的许多盾构机都是根据工程通过改造后反复使用[2]。为防止盾尾空隙中的泥水、土砂和注浆浆液从盾尾间隙渗漏到盾构机内部,严重时出现涌水涌砂甚至盾构机被埋的重大安全事故,必须要保障盾尾密封的效果[3]

    虽然盾构技术的发展已经较为成熟但是由于地质条件的复杂多变以及管理水平上的差异,盾尾刷密封损坏引起整个盾尾密封系统失效的情况时有发生,处置不当会造成严重的工程安全事故。2018年2月佛山市轨道交通2号线一期工程土建一标段盾构区间右线工地因盾尾密封系统被外部水土压力击穿出现严重渗漏,形成的透水涌砂通道从而使隧道周边地基被掏空,引发盾构机姿态出现突变造成管片隧道结构处于扭、拉、压和弯不利受力和大变形而坍塌,大量泥砂在丰富的地下水作用下形成了快速泥砂流涌入隧道,最终造成12人死亡,直接经济损失约5323.8万元[4]

    盾尾密封问题已经引起了工程界和学术界的不断重视。Chao等[5]根据一些工程事故实例中的结构与力学特征分析了多层盾尾密封系统在不同荷载类型下由于外部水土压力与油脂腔内压力不平衡以及盾构机姿态变化或施工不当导致的盾尾间隙增大进而导致盾尾密封失效的情况。关于盾尾刷的设计选型与性能试验方面,李奕等[6]总结出传统盾尾刷存在的一系列缺陷。同时结合各自所研究的实际工程对象,李奕等[6]、霍志光[7]对盾尾刷进行创新设计,通过一系列措施,对保护板、钢丝和压紧板结构及特性的优化,提高了盾尾刷的抗压性能。根据盾尾密封的实际失效形式,耐磨性与弹性是评估盾尾刷的两个关键性能。沈桂丽等[8]运用液压原理、PLC及运算放大控制原理,通过对动力参数和控制系统的详细研究,设计出一台套用于检测盾尾刷耐磨性和弹性的试验平台,通过测试耐磨性和弹性,评价盾尾刷的质量。Wei[9]设计一套可模拟不同盾尾间隙工况条件下的盾尾刷加载模拟试验装置, 并对盾尾刷在不同的盾尾间隙下的工作状态(加载变形和破坏模式)进行了试验研究。Ye等[10]探究了盾尾刷对管片的荷载分布,并对不同钢板厚度的钢板盾尾刷以及钢丝盾尾刷在盾尾间隙240 mm情况下所产生的盾尾刷贴合力进行试验。钟小春等[11]建立盾尾刷-油脂腔密封系统密封单元试验装置,探讨了盾尾密封性能的影响因素。

    以上总结中可以看出,目前主要针对盾尾刷设计和选材工程界开展一些实用性的研究工作,并设计了相关了的盾尾刷试验平台。但对于盾尾刷试验及成果分析总结,以及其性能评价方法缺乏系统的研究。本文基于中铁隧道集团设备公司开发的盾尾刷及盾尾刷试验平台开展了盾尾刷动态与静态的盾尾刷压缩试验,引入盾尾刷贴合力与弹性系数这对力学指标以及压紧板残余高度、残余角度等盾尾刷结构指标作为盾尾刷力学性能的评价指标,结合盾尾刷密封油脂逃逸试验,建立了盾尾密封性能与盾尾刷力学性能之间的关系,系统地探讨了盾尾刷性能评价方法。

    盾尾刷是保障盾尾密封的关键构件,焊接在盾壳上,与拼装管片共同形成2~3道环形空腔,并向腔中泵入油脂,确保盾尾密封效果及安全。

    试验中所用的盾尾刷样品的构造及细部尺寸见图 1,其主要由保护板(前板)、钢丝和压紧板(后板)3部分组成,前板由4块100 mm×160 mm×5 mm不锈钢钢板分两层构成,后板由2块150 mm×105 mm×5 mm和2块150 mm×205 mm×5 mm不锈钢钢板分两层错缝搭接构成,搭接宽度为50 mm。压紧板与水平底座钢板初始夹角55°,中间夹有4层15 mm厚度钢丝,呈阶梯状布置,相邻层长度相差20 mm。

    图  1  盾尾刷构造及细部尺寸
    Figure  1.  Structure and detailed dimensions of shield tail brush

    盾尾刷压缩试验平台由试验台、支架、固定盾尾刷的试验基座、液压油缸、压力传感器及显示器、操作台等主要部分组成,装置如图 2

    图  2  盾尾刷压缩试验平台示意图
    Figure  2.  Platform of compression tests on shield tail brush

    (1)将液压油缸升起,把待进行试验的盾尾刷固定在试验台上,确保安装位置准确并稳固。

    (2)通过操作台调节油缸行程上下限,控制单次压缩油缸往复时间,调节完毕即可开始试验。

    (3)通过直尺和量角器测量盾尾刷压紧板残余高度及其与底座平板残余角度。通过压力传感显示器,记录盾尾刷压缩过程中不同盾尾间隙下的接触压力数据,该压力可以反映盾尾刷与管片环外壁的贴合紧密程度,将此接触压力定义为盾尾刷贴合力。

    (4)记录盾尾刷经过1次、50次、100次、200次、400次、800次、1000次压缩时盾尾刷贴合力以及压缩后的压紧板残余高度及其与基座平板的残余角度等数据。

    盾尾刷静态压缩试验同样在盾尾刷压缩试验平台上进行。前3个步骤与动态压缩试验完全相同,静态压缩试验第4步为将盾尾刷压至设定的盾尾间隙并保持,记录0~48 h盾尾刷贴合力的数据变化。

    工程上常常将油脂损失量作为评价盾尾密封性能好坏的指标。为了将盾尾刷力学性能与油脂损失量联系起来,开展了盾尾刷密封系统油脂逃逸试验。油脂逃逸压力是指在一定密封油脂注入压力下,密封油脂开始从盾尾刷与管片之间的缝隙或者盾尾刷之间的缝隙发生被动逃逸时的油脂压强值。超过油脂逃逸压力时油脂的损耗将明显加大,因此研究油脂逃逸压力及油脂逃逸速率对于盾尾刷密封系统意义重大。很显然,油脂发生逃逸是油脂压力超过外部土水压力或者盾尾空隙浆液压力的结果,因此下文中的油脂压力考虑为油脂腔中的实际油脂压力与外部压力之差。

    本试验采用逐级增压的油脂压力施加方式,起始压力为0.3 MPa,压力梯度为0.3 MPa,共4级:0.3,0.6,0.9,1.2 MPa。测量在50,75,90 mm3种盾尾间隙在使用锥入度236.6的油脂时的逃逸压力、逃逸速率的变化规律,具体试验装置示意图及照片见图 3

    图  3  盾尾密封系统示意图及实物图
    Figure  3.  Schematic diagram and physical image of shield tail sealing system

    在实际盾构施工中,受盾构姿态变化、管片拼装等影响,盾尾间隙一般情况下在50~110 mm发生变化,不同盾尾间隙下盾尾刷的贴合力的不同,进而会影响到盾尾密封系统的密封性。

    为得到不同盾尾间隙状态下盾尾刷弹性性能的变化情况,本文采用盾尾刷动态压缩试验,设计4种盾尾间隙110,90,75,50 mm,测试其弹性性能变化情况。需要说明的是,在实际工程中盾尾间隙是盾尾与管片之间的间隙即盾尾内径与管片外径之差。在本试验中盾尾间隙是盾尾刷压缩产生最大变形时盾尾刷此时的残余高度,如图 4所示。

    图  4  盾尾刷压缩试验
    Figure  4.  Compression tests on shield tail brush

    (1)下压次数的影响

    不同盾尾间隙下盾尾刷的贴合力随下压次数变化情况试验结果如图 5。通过对比四种不同盾尾间隙下盾尾刷反复压缩后贴合力变化情况发现,随着下压次数的增多,不同盾尾间隙下的盾尾刷贴合力均有所减小,盾尾刷弹性性能均有不同程度的降低,且前期(0~200次动态压缩)盾尾刷贴合力损失最为集中,当下压次数达到400次以上时,贴合力变化不明显,超过600次以后基本稳定。

    图  5  贴合力随下压次数变化情况
    Figure  5.  Variation of adhesion force of shield tail brush with number of compression

    累计压缩至1000次时,盾尾间隙90 mm时,贴合力为396 N左右,约为初始值的65%;盾尾间隙75 mm时,贴合力为1000 N左右,约为初始值的77%;盾尾间隙位于最小值50 mm时,贴合力为1360 N左右,为初始值的90%。在常见的盾尾间隙60~90 mm,贴合力的损失率大约为10%~35%。

    (2)盾尾间隙的影响

    为进一步分析动态压缩试验中盾尾刷弹性性能的变化情况,将累计压缩1,50,200,400,800,1000次时的盾尾刷贴合力随压缩量(盾尾间隙)变化的情况加以统计,结果如图 6所示。

    图  6  盾尾间隙50~120 mm时贴合力变化情况
    Figure  6.  Variation of adhesion force of shield tail brush with shield tail gap of 50~120 mm

    图 6可以看出,贴合力在盾尾间隙90 mm出现了突变,随着盾尾间隙进一步减少,贴合力快速增加,当压缩至极限50 mm,贴合力增长又呈现了趋缓的趋势。这主要是因为盾尾刷由保护板、钢丝束和压紧板复合而成的构件。当盾尾间隙较大时,主要由保护板承担压力荷载;盾尾间隙进一步减少时,保护板、钢丝束和压紧板开始协同工作;进一步压缩至极限时,出现了明显了塑性,贴合力的增加趋势趋缓。

    盾构掘进过程中,盾尾拼装管片环会发生一定量的上浮[12],盾尾刷被控制在60~100 mm盾尾间隙随着盾构机姿态发生缓慢变化,可以认为处于相对静止压缩状态。盾尾刷的动态压缩试验能否反映这种工作受力特性,还需开展盾尾刷静态压缩试验。将盾尾刷在盾尾间隙为50,75 mm下静态压缩48 h,分析其保护板与管片贴合力随时间变化情况,试验结果如图 7中红色曲线所示。试验发现,盾尾刷在静压状态下,其保护板与管片的贴合力随时间的增加而逐渐减小,减小趋势由快变慢,当静压40 h以上时,贴合力变化已较为缓慢,趋于一个稳定值。

    图  7  动态试验与静态试验结果比较
    Figure  7.  Comparison between dynamic and static test results

    由于盾尾刷的静压试验时间成本高而动态试验则只需10 min即可完成,能否采用动态试验来对盾尾刷的力学特性性能评判,下面通过两者试验结果对比来说明。本小节以盾尾间隙50,75 mm为例进行比较,如图 7所示。

    从中可以看出,盾尾间隙50 mm情况下,盾尾刷初始弹性性能相近的情况下,盾尾刷贴合力的损失均集中于各自试验的前半段,其中对于动态试验,前500次压缩造成的贴合力的损失占总损失的80%,对于静态试验来说,前24 h静压造成的贴合力的损失占总损失的86%,而当试验结束时,动态压缩1000次与静压48 h后盾尾刷的残余盾尾刷贴合力分别为1362.2,1366 N,盾尾间隙75 mm情况下则分别为980,1000 N,盾尾刷贴合力基本一致。因此,可以采用动态压缩试验代替静态压缩试验开展盾尾刷力学性能试验研究,提高试验效率。

    盾尾密封系统油脂逃逸试验记录了在盾尾间隙50,75,90 mm时,对于锥入度236.6的油脂,油脂逃逸速度与油脂压力关系,如图 8。油脂逃逸速率随着油脂压力的逐级增加而不断上升,当盾尾间隙为50 mm(盾尾刷贴合力较大)时,油脂逃逸速率与油脂压力基本呈线性增加趋势,而盾尾间隙为75,90 mm(盾尾刷贴合力较小)时,油脂逃逸速率随油脂压力升高呈非线性增大趋势,说明此时盾尾刷贴合力不足以平衡油脂腔内压力,引起密封失效。

    图  8  油脂逃逸速度与油脂压力关系
    Figure  8.  Relationship between grease escape speed and grease pressure

    盾尾密封系统油脂逃逸试验中采用的前后各并排两个的盾尾刷,其宽为400 mm,将试验测得的油脂逃逸速率换算成每环油脂逃逸量(对于常规盾构直径6.28 m),得出不同油脂压力下,每环油脂逃逸量与盾尾间隙的关系,如图 9

    图  9  每环油脂逃逸量与盾尾间隙关系
    Figure  9.  Relationship between grease escape per ring and shield tail gap

    实际工程中很多盾尾刷的损坏属于非正常磨损引起的,结合实际工程经验,总结出3种盾尾刷破坏模式。

    (1)盾尾刷无显著磨损下油脂压力不足抵抗外部水土压力引起的密封失效。

    (2)盾尾刷非正常损坏引起油脂逃逸导致油脂压力不足以抵抗外部水土压力引起的密封失效。

    (3)盾尾刷被盾尾空隙中的砂浆侵入导致盾尾刷弹性丧失引起的密封失效。

    后两种破坏模式的工程实际图片如图 10所示。本试验针对的是盾尾刷未损坏以及未硬化的情况。

    图  10  两种典型的盾尾刷破坏方式
    Figure  10.  Two typical damage patterns of shield tail brush

    实际工程中,由于盾构机姿态和管片姿态处于动态变化过程中,其盾尾间隙是也往往是变化的,其范围主要为通常为60~90 mm,极限时可达到50~110 mm。因此盾尾刷弹性变形能力应按照盾尾间隙50~110 mm来讨论。

    盾尾刷的弹性系数定义为单位弹性变形下所需要施加的贴合力,即盾尾刷的弹性系数为其贴合力与弹性变形的比值。动态压缩800次后盾尾刷的弹性系数与盾尾间隙的关系如图 11所示。

    图  11  弹性系数与盾尾间隙关系
    Figure  11.  Relationship between elastic coefficient and shield tail gap

    从中可以看出,弹性系数与盾尾间隙的变化趋势整体可以分为盾尾间隙50~90 mm与90~110 mm两段线性变化,盾尾刷的弹性系数随盾尾间隙的减少而增大,其变化速率会在盾尾间隙90 mm前后发生突变。这是由于盾尾刷作为组合材料,当盾尾间隙较大,即变形量较小时,主要是由保护板进行受力;而当盾尾间隙逐渐较小,盾尾刷保护板、压紧板、钢丝刷等各个部分均逐渐参与受力,因此弹性系数也逐渐增大。

    根据初始盾尾刷高度H0、盾尾间隙S、弹性变形占比α和弹性系数K可以按照下列公式计算盾尾刷的贴合力F

    F=Kα(H0S) (1)

    压紧板高度和压紧板与底座角度的变化情况也是评价盾尾刷弹性性能损失的重要参考因素,与盾尾刷贴合力相比是更为直观具象的指标。试验中,不同盾尾间隙下累计下压至1000次过程中盾尾刷压紧板高度和其与底座角度值分别见图 1213所示。

    图  12  压紧板残余高度变化情况
    Figure  12.  Variation of residual height of pressure plate
    图  13  压紧板残余角度变化情况
    Figure  13.  Variation of residual angle of pressure plate

    通过对不同工况下的盾尾刷反复动态压缩可以发现,盾尾刷压紧板的高度和其与底座所呈角度有明显减小趋势,尤其以前100次压缩过程中减小趋势最为明显,经过1000次压缩后压紧板高度与压紧板和底座所呈角度趋于均趋于稳定,其值称为压紧板残余高度Hr、压紧板残余角度ψr,其中最不利状况下(最小盾尾间隙50 mm,反复压缩1000次),压紧板高度从初始192 mm下降至135 mm,其与底座所呈角度也从初始55.5°下降至35.5°。

    盾尾刷单次压缩所产生的塑性变形量为前一次残余量与后一次残余量之差,进而通过压紧板高度计算得出盾尾刷的塑性变化量与压缩次数关系,如图 14。可以看出盾尾刷的塑性变形基本集中于前200次内,前200次压缩所造成的塑性变形占总塑性变形的85%。选取各个盾尾间隙下压1000次后的总塑性变形量与盾尾间隙作图,如图 15。可以看出总塑性变形量与盾尾间隙基本呈线性降低关系,各个盾尾间隙下弹性变形与塑性变形占总变形比例如图 16所示。弹性变形对总变形量的占比基本随着盾尾间隙的减少几乎线性减少,而塑性变形的占比情况则刚好相反。当盾尾间隙处于110 mm左右,弹性变形占比约为90%,此时盾尾刷的压缩主要为弹性变形;当盾尾间隙压缩至50 mm时,弹性变形占比仅为60%。

    图  14  盾尾刷塑性变形量与压缩次数关系
    Figure  14.  Relationship between plastic deformation of shield tail brush and number of compressions
    图  15  总塑性变形量与盾尾间隙关系
    Figure  15.  Relationship between total plastic deformation and shield tail gap
    图  16  弹性变形与塑性变形占比
    Figure  16.  Ratios of elastic to plastic deformation

    总结盾尾刷动静态压缩试验,对于试验盾尾刷,不同盾尾间隙下的各个评价指标的对应关系见表 1

    表  1  盾尾刷力学性能指标
    Table  1.  Mechanical parameters of shield tail brush
    盾尾间隙/mm 盾尾刷贴合力/N 压紧板残余高度/mm 压紧板残余角度/(°) 弹性系数/(N·mm-1)
    50 1362.2 135 35.6 16
    75 997.2 159 44.5 11.9
    90 392 174 48 4.7
    110 294 183 51.4 4.02
    下载: 导出CSV 
    | 显示表格

    根据压紧板残余高度Hr、盾尾间隙S、弹性系数K可以按照下列公式计算盾尾刷的贴合力F

    F=K(HrS) (2)

    以标准盾尾间隙75 mm为例,此时Hr=159 mm,S=75 mm,K=11.9 N/mm,代入式(2),计算得出盾尾刷贴合力F=999.6 N,与动态压缩试验中数据一致。因此,可以根据弹性系数采用式(2)计算计算获得盾尾刷贴合力。工程中一般要求最终盾尾刷的压紧板残余角度大于35°,这说明本试验盾尾刷性能优良。

    盾尾密封能性主要依赖于盾尾油脂压力来抵抗盾尾空隙中的砂浆、泥水压力;而盾尾油脂压力的压力保持则依赖于盾尾刷与管片外壁之间的贴合力大小。盾尾密封工作性能如图 17所示,为了保证盾尾密封系统在盾构机和管片姿态不断动态变化中的密封性,需要确保盾尾刷贴合力不低于某一允许值,并且要维持稳定。工程上一般会使用油脂允许消耗量作为盾尾密封性的评价指标。一般来说,对于直径为6.28 m的常规盾构机掘进一环一般需耗时2 h,消耗油脂约为50~80 L[13]

    图  17  盾尾密封工作性能流程图
    Figure  17.  Flow chart of service behavior of shield tail sealing

    参考第3节所得出的每环油脂逃逸量与盾尾间隙的关系,同时考虑盾尾间隙又与盾尾刷贴合力、压紧板残余高度等评价指标密切相关,可以通过盾尾间隙将盾尾刷贴合力和每环油脂逃逸量建立关系。以油脂压力0.6 MPa为例,将每环油脂逃逸量60 L作为允许消耗量,其对应盾尾间隙约为83 mm,该盾尾间隙对应的盾尾刷贴合力约为710 N。由此做出盾尾密封系统的安全区域,如图 18中阴影所示。安全区域含义为当盾尾间隙位于50~83 mm,此时盾尾刷贴合力能够有效减少油脂的逃逸(保证每环掘进过程中油脂逃逸量小于允许值),且保持住油脂腔内压力,更好地平衡盾尾空隙中的泥水、注浆压力,形成良好的盾尾密封。因此可将710 N作为油脂压力0.6 MPa时保障盾尾密封性的盾尾刷贴合力允许值。

    图  18  油脂逃逸量与盾尾刷评价指标关系(0.6 MPa)
    Figure  18.  Relationship between grease consumption and evaluation parameter of shield tail brush(0.6 MPa)

    通过盾尾密封性能试验建立了油脂逃逸速率与盾尾间隙之间的关系,以及盾尾刷力学试验建立了盾尾刷的贴合力、压紧板残余高度、压紧板残余角度、盾尾刷弹性系数与盾尾间隙的关系,结合以上两方面的试验成果,从而可以只通过盾尾刷力学性能对盾尾密封系统的密封性能进行评价,避免开展盾尾密封逃逸试验。

    (1)盾尾间隙50 mm情况下,1000次动态压缩与48静态压缩后盾尾刷贴合力分别为1362.2,1366 N,盾尾间隙75 mm情况下则分别为980, 1000 N,基本一致。因此,为提高试验测试效率,可采用动态压缩1000次后来评价盾尾刷力学性能。

    (2)盾尾刷贴合力受压缩次数显著影响,其中随着下压次数的增多,盾尾刷贴合力不断减小,减小速率由快变慢。前期(0~200次动态压缩)盾尾刷贴合力损失最为集中,当累计下压次数达到600次以后,盾尾刷的各项弹性性能参数减小已非常缓慢,趋于稳定,总的贴合力损失约占初始贴合力的10%~35%。

    (3)盾尾刷贴合力受盾尾间隙显著影响。当盾尾间隙压缩在90 mm以上范围内反复压缩盾尾刷时,贴合力随盾尾间隙的压缩变化小;当盾尾间隙进一步压缩至75 mm时,由于盾尾刷各部分参与受力,盾尾刷贴合力快速增加,进一步压缩至50 mm时,盾尾刷出现了明显了塑性,贴合力的增加趋势趋缓。

    (4)盾尾刷的弹性变形占比随着盾尾间隙减少呈线性降低,相反其塑性变形占比随着盾尾间隙减少呈线性增加。当盾尾间隙110 mm,弹性和塑性变形占比分别为90%,10%;当盾尾间隙50 mm,弹性和塑性变形占比分别为60%,40%。

    (5)盾尾刷弹性系数与盾尾间隙的变化趋势整体可以分为盾尾间隙50~90 mm,90~110 mm两段。总的趋势是盾尾刷的弹性系数随盾尾间隙的减少而增大。盾尾刷模量在盾尾间隙90 mm前后发生斜率突变与其贴合力在90 mm突变是一致的。

    (6)形成一套关于密封盾尾刷性能的评价方法。首先通过油脂逃逸试验及每环管片掘进中的油脂允许损耗量判断出盾尾密封所需要最小贴合力,其次结合盾尾刷的弹性变形占比、盾尾间隙允许变化范围和盾尾刷弹性系数等评价指标确定盾尾刷贴合力,最后对比获得盾尾刷密封安全区域,以此来判断所选用的盾尾刷是否满足盾尾密封性能的要求。

  • 图  1   近场动力学相互作用的物质点

    Figure  1.   Material points of peridynamic interactions

    图  2   近场动力学变形示意图

    Figure  2.   Schematic diagram of peridynamic deformation

    图  3   PD模型离散示意图

    Figure  3.   Diagrammatic sketch of PD model discretization

    图  4   均质方板的几何模型和边界条件

    Figure  4.   Geometric modeling and boundary conditions of homogeneous square plate

    图  5   不同时刻m收敛结果

    Figure  5.   Results of m-convergence at different time

    图  6   不同时刻δ收敛结果

    Figure  6.   Results of δ-convergence at different time

    图  7   LdB花岗岩热损伤破裂模型

    Figure  7.   Thermal damage fracture model for LdB granite

    图  8   LdB花岗岩破裂声发射结果[25]

    Figure  8.   Acoustic emission results of LdB granite

    图  9   损伤破裂和温度分布演化过程的PD模拟结果

    Figure  9.   Simulated results of damage fracture and evolution process of temperature distribution by PD

    图  10   文献[3132]的花岗岩热破裂声发射结果

    Figure  10.   Acoustic emission results by Ishida et al

    图  11   不同温度梯度弛豫时间条件下试件损伤率变化

    Figure  11.   Variation of damage rate of specimen under different temperature gradient relaxation time conditions

    图  12   不同温度梯度弛豫时间时x = 0处的温度分布

    Figure  12.   Distribution of temperature at x = 0 under different temperature gradient relaxation time conditions

    图  13   不同热流弛豫时间条件下试件损伤率变化

    Figure  13.   Variation of damage rate of specimen under different heat flux relaxation time conditions

    图  14   不同热流弛豫时间时x = 0处的温度分布

    Figure  14.   Distribution of temperature at x = 0 under different heat flux relaxation time conditions

    表  1   瞬态热传导离散参数值

    Table  1   Discrete parameter values for transient heat conduction

    条件 物质点分布 点间距/m 近场半径/m 非局部比
    m
    收敛
    250×250 0.040 0.060 1.5
    400×400 0.025 0.060 2.4
    500×500 0.020 0.060 3.0
    δ
    收敛
    250×250 0.040 0.120 3.0
    400×400 0.025 0.075 3.0
    500×500 0.020 0.060 3.0
    下载: 导出CSV
  • [1]

    ZHANG J Y, SHEN Y J, YANG G S, et al. Inconsistency of changes in uniaxial compressive strength and P-wave velocity of sandstone after temperature treatments[J]. Journal of Rock Mechanics and Geotechnical Engineering, 2021, 13(1): 143-153. doi: 10.1016/j.jrmge.2020.05.008

    [2] 郤保平, 吴阳春, 王帅, 等. 热冲击作用下花岗岩力学特性及其随冷却温度演变规律试验研究[J]. 岩土力学, 2020, 41(增刊1): 83-94.

    XI Baoping, WU Yangchun, WANG Shuai, et al. Experimental study on mechanical properties of granite under thermal shock and its evolution with cooling temperature[J]. Rock and Soil Mechanics, 2020, 41(S1): 83-94. (in Chinese)

    [3]

    WANG G Y, YANG D, LIU S W, et al. Experimental study on the anisotropic mechanical properties of oil shales under real-time high-temperature conditions[J]. Rock Mechanics and Rock Engineering, 2021, 54(12): 6565-6583. doi: 10.1007/s00603-021-02624-7

    [4]

    YANG S Q, HUANG Y H, TIAN W L, et al. Effect of high temperature on deformation failure behavior of granite specimen containing a single fissure under uniaxial compression[J]. Rock Mechanics and Rock Engineering, 2019, 52(7): 2087-2107. doi: 10.1007/s00603-018-1725-5

    [5]

    ZHANG F, ZHAO J J, HU D W, et al. Laboratory investigation on physical and mechanical properties of granite after heating and water-cooling treatment[J]. Rock Mechanics and Rock Engineering, 2018, 51(3): 677-694. doi: 10.1007/s00603-017-1350-8

    [6] 唐世斌, 唐春安, 朱万成, 等. 热应力作用下的岩石破裂过程分析[J]. 岩石力学与工程学报, 2006, 25(10): 2071-2078. doi: 10.3321/j.issn:1000-6915.2006.10.019

    TANG Shibin, TANG Chun'an, ZHU Wancheng, et al. Numerical investigation on rock failure process induced by thermal stress[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(10): 2071-2078. (in Chinese) doi: 10.3321/j.issn:1000-6915.2006.10.019

    [7]

    JIAO Y Y, ZHANG X L, ZHANG H Q, et al. A coupled thermo-mechanical discontinuum model for simulating rock cracking induced by temperature stresses[J]. Computers and Geotechnics, 2015, 67: 142-149. doi: 10.1016/j.compgeo.2015.03.009

    [8] 严成增. FDEM-TM方法模拟岩石热破裂[J]. 岩土工程学报, 2018, 40(7): 1198-1204. doi: 10.11779/CJGE201807005

    YAN Chengzeng. Simulating thermal cracking of rock using FDEM-TM method[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(7): 1198-1204. (in Chinese) doi: 10.11779/CJGE201807005

    [9]

    YAO J, ZHANG X, SUN Z X, et al. Numerical simulation of the heat extraction in 3D-EGS with thermal-hydraulic- mechanical coupling method based on discrete fractures model[J]. Geothermics, 2018, 74: 19-34. doi: 10.1016/j.geothermics.2017.12.005

    [10]

    SILLING S A. Reformulation of elasticity theory for discontinuities and long-range forces[J]. Journal of the Mechanics and Physics of Solids, 2000, 48(1): 175-209. doi: 10.1016/S0022-5096(99)00029-0

    [11]

    SILLING S A, ASKARI E. A meshfree method based on the peridynamic model of solid mechanics[J]. Computers and Structures, 2005, 83(17/18): 1526-1535.

    [12]

    GAO C L, ZHOU Z Q, LI L P, et al. Strength reduction model for jointed rock masses and peridynamics simulation of uniaxial compression testing[J]. Geomechanics and Geophysics for Geo-Energy and Geo-Resources, 2021, 7(2): 34. doi: 10.1007/s40948-021-00232-x

    [13] 马鹏飞, 李树忱, 袁超, 等. 基于SED准则的近场动力学及岩石类材料裂纹扩展模拟[J]. 岩土工程学报, 2021, 43(6): 1109-1117. doi: 10.11779/CJGE202106014

    MA Pengfei, LI Shuchen, YUAN Chao, et al. Simulations of crack propagation in rock-like materials by peridynamics based on SED criterion[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(6): 1109-1117. (in Chinese) doi: 10.11779/CJGE202106014

    [14] 马鹏飞, 李树忱, 周慧颖, 等. 岩石材料裂纹扩展的改进近场动力学方法模拟[J]. 岩土力学, 2019, 40(10): 4111-4119.

    MA Pengfei, LI Shuchen, ZHOU Huiying, et al. Simulations of crack propagation in rock-like materials using modified peridynamic method[J]. Rock and Soil Mechanics, 2019, 40(10): 4111-4119. (in Chinese)

    [15]

    ZHANG H, QIAO P Z. An extended state-based peridynamic model for damage growth prediction of bimaterial structures under thermomechanical loading[J]. Engineering Fracture Mechanics, 2018, 189: 81-97. doi: 10.1016/j.engfracmech.2017.09.023

    [16]

    WANG Y T, ZHOU X P, KOU M M. An improved coupled thermo-mechanic bond-based peridynamic model for cracking behaviors in brittle solids subjected to thermal shocks[J]. European Journal of Mechanics-A, 2019, 73: 282-305. doi: 10.1016/j.euromechsol.2018.09.007

    [17]

    WANG Y T, ZHOU X P. Peridynamic simulation of thermal failure behaviors in rocks subjected to heating from boreholes[J]. International Journal of Rock Mechanics and Mining Sciences, 2019, 117: 31-48. doi: 10.1016/j.ijrmms.2019.03.007

    [18]

    YANG Z, YANG S Q, CHEN M. Peridynamic simulation on fracture mechanical behavior of granite containing a single fissure after thermal cycling treatment[J]. Computers and Geotechnics, 2020, 120: 103414. doi: 10.1016/j.compgeo.2019.103414

    [19]

    YANG Z, YANG S Q, TIAN W L. Peridynamic simulation of fracture mechanical behaviour of granite specimen under real-time temperature and post-temperature treatments[J]. International Journal of Rock Mechanics and Mining Sciences, 2021, 138: 104573. doi: 10.1016/j.ijrmms.2020.104573

    [20]

    CATTANEO C. A form of heat conduction equation which eliminates the paradox of instantaneous propagation[J]. Compte Rendus, 1958, 247: 431-433.

    [21]

    VERNOTTE P. Paradoxes in the continuous theory of the heat conduction[J]. Compte Rendus, 1958, 246: 3154-3155.

    [22]

    TZOU D Y. A unified field approach for heat conduction from macro- to micro-scales[J]. Journal of Heat Transfer, 1995, 117(1): 8-16. doi: 10.1115/1.2822329

    [23]

    AL-NIMR M A, NAJI M, ABDALLAH R I. Thermal behavior of a multi-layered thin slab carrying periodic signals under the effect of the dual-phase-lag heat conduction model[J]. International Journal of Thermophysics, 2004, 25(3): 949-966. doi: 10.1023/B:IJOT.0000034247.32646.d4

    [24]

    WANG L J, XU J F, WANG J X. A peridynamic framework and simulation of non-Fourier and nonlocal heat conduction[J]. International Journal of Heat and Mass Transfer, 2018, 118: 1284-1292. doi: 10.1016/j.ijheatmasstransfer.2017.11.074

    [25]

    JANSEN D P, CARLSON S R, YOUNG R P, et al. Ultrasonic imaging and acoustic emission monitoring of thermally induced microcracks in Lac du Bonnet granite[J]. Journal of Geophysical Research: Solid Earth, 1993, 98(B12): 22231-22243. doi: 10.1029/93JB01816

    [26]

    MADENCI E, OTERKUS E. Peridynamic Theory and Its Applications[M]. New York: Springer, 2014.

    [27]

    HOSSEINI-TEHRANI P, ESLAMI M R. BEM analysis of thermal and mechanical shock in a two-dimensional finite domain considering coupled thermoelasticity[J]. Engineering Analysis with Boundary Elements, 2000, 24(3): 249-257. doi: 10.1016/S0955-7997(99)00063-6

    [28]

    WANG Y T, ZHOU X P, XU X. Numerical simulation of propagation and coalescence of flaws in rock materials under compressive loads using the extended non-ordinary state-based peridynamics[J]. Engineering Fracture Mechanics, 2016, 163: 248-273. doi: 10.1016/j.engfracmech.2016.06.013

    [29]

    WANNE T S, YOUNG R P. Bonded-particle modeling of thermally fractured granite[J]. International Journal of Rock Mechanics and Mining Sciences, 2008, 45(5): 789-799. doi: 10.1016/j.ijrmms.2007.09.004

    [30]

    GOFFREDO DE PORTU. Introduction to Mechanical Behaviour of Ceramics[M]. Faenza: CNR, 1992.

    [31]

    ISHIDA T, KINOSHITA N, WAKABAYASHI N. Acoustic emission monitoring during thermal cracking of a granite block heated in a center hole[C]//Proceedings of the 3rd Asian Rock Mechanics Symposium. Kyoto, 2004.

    [32]

    HUANG X, TANG S B, TANG C A, et al. Numerical simulation of cracking behavior in artificially designed rock models subjected to heating from a central borehole[J]. International Journal of Rock Mechanics and Mining Sciences, 2017, 98: 191-202. doi: 10.1016/j.ijrmms.2017.07.016

    [33] 胡学功, 刘登瀛, 周定伟. 细陶瓷棒的脉冲热流波动传播行为[J]. 应用科学学报, 2002, 20(3): 305-308. doi: 10.3969/j.issn.0255-8297.2002.03.019

    HU Xuegong, LIU Dengying, ZHOU Dingwei. The behavior of pulse heat flux wave propagation in a fine round ceramic bar[J]. Journal of Applied Sciences, 2002, 20(3): 305-308. (in Chinese) doi: 10.3969/j.issn.0255-8297.2002.03.019

    [34]

    KAMINSKI W. Hyperbolic heat conduction equation for materials with a nonhomogeneous inner structure[J]. Journal of Heat Transfer, 1990, 112(3): 555-560. doi: 10.1115/1.2910422

    [35]

    GUO S L, WANG B L. Thermal shock fracture of a cylinder with a penny-shaped crack based on hyperbolic heat conduction[J]. International Journal of Heat and Mass Transfer, 2015, 91: 235-245. doi: 10.1016/j.ijheatmasstransfer.2015.07.081

    [36]

    GUO S L, WANG B L. Thermal shock cracking behavior of a cylinder specimen with an internal penny-shaped crack based on non-fourier heat conduction[J]. International Journal of Thermophysics, 2016, 37(2): 17. doi: 10.1007/s10765-015-2029-6

  • 期刊类型引用(6)

    1. 刘耀徽,恩和得力海,黄亚平,张耘获,杨坪. 时移跨孔地震全波形反演方法在水泥搅拌桩检测中的应用研究. 地球物理学进展. 2024(05): 2078-2089 . 百度学术
    2. 张丁宁. 上海滨海平原区土的热物理指标与静力触探Ps值的关系研究. 城市道桥与防洪. 2022(07): 197-199+25 . 百度学术
    3. 贾剑青,赵阳阳,贾超,辛成平,张帮鑫. 湿陷性黄土地基水泥土搅拌桩加固效果研究. 铁道工程学报. 2022(07): 18-24 . 百度学术
    4. 何杨闽. 物联网全过程管控水泥搅拌桩关键技术研究. 珠江水运. 2021(11): 36-38 . 百度学术
    5. 王建梅. 高铁软弱斜坡段路基处理技术研究. 施工技术. 2020(07): 49-53 . 百度学术
    6. 赵洋洋,杨昌民. 基于RBF神经网络对基桩完整性的预测. 施工技术. 2020(S1): 69-72 . 百度学术

    其他类型引用(1)

图(14)  /  表(1)
计量
  • 文章访问数:  319
  • HTML全文浏览量:  46
  • PDF下载量:  65
  • 被引次数: 7
出版历程
  • 收稿日期:  2023-08-13
  • 网络出版日期:  2024-03-24
  • 刊出日期:  2024-10-31

目录

/

返回文章
返回