盾构渣土改良效果评价机器人搅拌探头选型分析.pdf

引用格式:王祥祥,鞠翔宇,荆留杰,等.盾构渣土改良效果评价机器人搅拌探头选型分析J.隧道建设(中英文),2023,43(8):1413.WANG Xiangxiang,JU Xiangyu,JING Liujie.Mixing probe selection of evaluation robot for soil conditioning effect of shieldJ.Tunnel Construction,2023,43(8):1413.收稿日期:2022-07-15;修回日期:2023-05-06基金项目:国家重点研发计划(2018YFB1702504)第一作者简介:王祥祥(1994),男,河南商丘人,2019 年毕业于中国地质大学(北京),地质工程专业,硕士,工程师,现从事隧道掘进机智能感知终端研制工作。E-mail:cugbwxx 。通信作者:鞠翔宇,E-mail:619188157 。盾构渣土改良效果评价机器人搅拌探头选型分析王祥祥1,鞠翔宇1,2,荆留杰1,杨 晨1,游宇嵩1(1.中铁工程装备集团有限公司,河南 郑州 450016;2.山东大学机械工程学院,山东 济南 250061)摘要:为实现土压平衡盾构渣土改良效果的动态监测,须开展渣土改良效果实时评价机器人的研究。搅拌探头作为机器人唯一与渣土接触的部件,其形状对搅拌转矩与贯入阻力检测的灵敏度至关重要。将盾构渣土假定为 Bingham 流体,开展不同形状搅拌探头贯入搅拌不同流变参数渣土的正交仿真试验,观测探头贯入、搅拌渣土过程中贯入阻力与搅拌转矩的变化规律。基于仿真结果,利用层次分析法对搅拌探头进行选型分析,得到可以同时兼顾贯入阻力与搅拌转矩检测灵敏度要求的搅拌探头,并通过室内试验对选型结果进行验证。关键词:土压平衡盾构;渣土改良机器人;搅拌探头;SPH 方法;Bingham 流体DOI:10.3973/j.issn.2096-4498.2023.08.016文章编号:2096-4498(2023)08-1413-12中图分类号:U 45 文献标志码:A开放科学(资源服务)标识码(OSID):MMi ix xi in ng g P Pr ro ob be e S Se el le ec ct ti io on n o of f E Ev va al lu ua at ti io on n R Ro ob bo ot t f fo or r S So oi il l C Co on nd di it ti io on ni in ng g E Ef ff fe ec ct t o of f S Sh hi ie el ld dWANG Xiangxiang1,JU Xiangyu1,2,*,JING Liujie1,YANG Chen1,YOU Yusong1(1.China Railway Engineering Equipment Group Co.,Ltd.,Zhengzhou 450016,Henan,China;2.School of Mechanical Engineering,Shandong University,Jinan 250061,Shandong,China)A Ab bs st tr ra ac ct t:A real-time evaluation robot for soil conditioning effect is investigated to realize dynamic monitoring of the soil conditioning effect of the earth-pressure balance shield.As the only part of the robot that is in contact with the muck,the shape of the stirring probe is of great importance to the sensitivity of the mixing torque and penetration resistance detection.In this study,the shields muck is assumed to be the Bingham fluid,and the orthogonal simulation experiments are conducted to observe the resistance changes when mixing probe penetration and the torque,where the soil is mixed under different rheological parameters.Therefore,based on the simulation results,the selection analysis of the mixing probe is conducted using the analytic hierarchy process method to obtain a mixing probe that could consider both penetration resistance and detection sensitivity requirements of mixing torque.Finally,the selection results are verified by laboratory tests.K Ke ey yw wo or rd ds s:earth-pressure balance shield;soil conditioning evaluation robot;mixing probe;smoothed particle hydrodynamics method;Bingham fluid0 引言近年来,随着我国基础设施建设的快速发展,盾构技术在推动地铁快速建设方面发挥了重要作用1。目前,主流的盾构类型包括土压平衡盾构和泥水平衡盾构。其中,渣土作为土压平衡盾构压力平衡的媒介,一方面可以平衡前方掌子面的水土压力,防止土体失稳;另一方面可以平衡螺机压力,保证设备顺利出渣2。因此,往往需要通过渣土改良使渣土具有良好的流动性、合适的塑性、较低的抗剪强度和黏附强度、较小的渗透系数和一定的压缩性3-7,从而便于排土隧道建设(中英文)第 43 卷和控制土舱压力,降低刀盘结泥饼、刀具磨损、螺机喷涌的风险8。近年来,国内外学者对渣土塑流性的测试方法与评价指标进行了大量的研究。Zumsteg 等9与 Peila等10通过研究发现,黏性渣土中加入泡沫后,其坍落度值明显增加,渣土塑流性增强。闫鑫等11发现当泡沫注入比相同时,砂土含水率低于某值可增加土颗粒间黏聚力,过量水会降低颗粒间黏聚力,导致坍落度增大。叶新宇12拟合大量坍落度试验数据,指出泥质粉砂岩地层盾构渣土的坍落度值与泡沫注入比和含水率呈二次函数关系。王树英等8针对改良后渣土理想坍落度取值范围进行了广泛调研,发现国内外对其的取值范围尚未统一,且差别较大。李培楠等13提出使用水泥胶砂流动度测定仪与稠度仪评价黏性渣土的塑流性,并将试验结果得到的最佳配比应用到上海某地铁盾构隧道施工中,取得了良好的施工效果。Oliveira等14对水泥胶砂流动度仪的试验方法进行了改进,使用渣土在玻璃板上方上下震动 40 次后的直径变化率F40作为渣土塑流性的评价指标,但由于数据量有限,且影响因素较多,作者未给出理想渣土的 F40值。Milligan3与 Maidl 等7用黏稠指数来表征黏土的黏稠状态,进而评价其塑流性,得出稠度指数在 0.4 0.75 区间的渣土具有较好的塑流性,能够较好地满足盾构排土的要求。李培楠等13和 Langmaack15采用砂浆稠度仪测定渣土的塑流性,结果表明随着改良剂掺入比的增加,圆锥的插入深度越大,渣土的塑流性越好。Quebaud 等16与刘大鹏17在室内砂浆搅拌机的基础上增加了功率表,测定搅拌渣土时消耗的能量,并将能量消耗减小量作为渣土塑流性的评价指标。Mori等18研制了一套可以测定在一定压力下渣土抗剪强度的十字剪切试验装置,用于检测不同改良参数下的最大 转 矩 值,以 此 计 算 改 良 土 样 的 抗 剪 强 度。Avunduk 等19提出了一种可手持便捷式扭转剪切试验装置,并通过试验验证其测量结果与经典剪切试验具有很强的相关性,该装置能够快速简便地评价渣土的改良效果。以上研究在一定程度上降低了土压平衡盾构的施工风险,增加了盾构掘进效率,但是不同学者使用的试验仪器和评价指标各不相同,目前缺乏统一的评价标准;此外,由于缺乏盾构智能化方面的深入研究与工程实践,盾构掘进过程中渣土改良状态信息获取滞后,决策控制依赖人为经验致使盾构盲目掘进,导致渣土改良策略与地层条件不适应、掘进效率低等问题一直存在8。因此,有必要研发一款室内渣土改良标准量化试验台,用于改良规律挖掘与评价指标体系的建立;与一款盾构渣土改良效果评价机器人,用于渣土实时取样与检测,并结合掘进参数,指导改良参数的快速调整。机器人与试验台皆通过安装在其上的搅拌装置搅拌渣土,检测搅拌探头插入渣土时的贯入阻力以及探头旋转时的转矩,从而评价渣土的塑流性。搅拌探头作为机器人与试验台中直接与渣土接触的装置,其形状决定了所测贯入阻力与转矩的敏感性。国内外学者对不同搅拌探头的测量效果进行了一定的探索。Li 等20使用十字形搅拌探头测量了气压加载条件下不同贯入深度时改良砂土的搅拌转矩;Mario 等21采用球形搅拌探头(带有球测量系统的流变仪)研究了泡沫改良土的流变特性,得到改良土剪切应力-剪切速率关系;Wan 等22采用圆盘形探头进行了切向黏附试验,探索改良土与圆盘探头之间的黏附强度和转矩变化关系。以上研究为不同学者分别针对单一形状搅拌探头进行的试验,未对不同形状探头的测量效果进行对比分析。Hu 等23使用 CFD 方法对 3 种不同形状探头进行仿真,对比转矩测量效果,但未对探头贯入阻力进行研究。综上所述,本文将渣土看作非牛顿流体,通过平滑粒子 流 体 动 力 学(smoothed particle hydrodynamics,SPH)方法对不同形状搅拌探头贯入渣土以及搅拌过程进行仿真,并在仿真结果的基础上利用层次分析法对不同形状搅拌探头测量的搅拌转矩与贯入阻力敏感性进行分析,以期为机器人与试验台的研发应用提供帮助。1 SPH 理论与渣土流变特性1.1 SPH 基本理论平滑粒子流体动力学是一种无网格或自由网格的数值方法,其不需要定义有限元通常定义的节点和单元,只需用一个点集合代表一个给定物体,SPH 中定义的点通常指粒子或虚拟粒子。SPH 法的核心实为一种插值,宏观物理变量(密度、压力、温度)f(x)能非常方便地借助于一组无序点上的值表示成积分插值计算得到,其方程的构造涉及2 个关键问题24-25:函数的光滑近似逼近、质点的近似逼近。1.1.1 函数的光滑近似逼近应用积分表示函数的核近似法对任意函数和光滑核函数进行逐步积分如下:f(x)=Df(x)W(x-x,h)dx。(1)式中:f(x)为任意空间变量的函数;D 为 x 的积分区间;W(x-x,h)为光滑函数,又被称作插值函数;h为定义光滑函数的影响区域的光滑长度。1.1.2 质点的近似逼近在 SPH 方法中,整个系统是由具有独立的质量、占有独立空间的有限个粒子表示的。应用粒子近似法可让与核近似法相关的连续积分表示式转化为支持域4141第 8 期王祥祥,等:盾构渣土改良效果评价机器人搅拌探头选型分析内所有粒子叠加求和的离散化形式26:f(x)=Ni=1miif(xi)W(x-xi,h)。(2)式中:mi、i分别为粒子的质量和密度;i 为计算质点记号;N 为质点总数。通过式(1)和式(2)便可对任意函数 f(x)在空间坐标内进行 SPH 方法的离散化处理。转化后的不可压缩流体的 Navier-Stokes 公式的形式如下。质量守恒方程为:i=Nj=1mjWijNj=1mjj()Wij。(3)动量守恒方程为:dvidt=Nj=1mji+jijWijxi。(4)能量守恒方程为:deidt=12Nj=1mjPi+PjijvijWijxi+i2iij。(5)式(3)(5)中:v为速度;e为热能;为应力张量;P为压力;为动力黏滞系数;为动剪应变率;和 代表不同坐标方向;i和j代表不同质点,其中i为计算质点记号,j 为邻近质点记号。SPH 方法可以像其他拉格朗日模型那样指定初始条件与边界条件,也允许与其他拉格朗日体的接触相互作用。在变形并非很严重时,SPH 方法的精度比拉格朗日有限元分析要低,但是在诸如液体晃动、波浪工程等大变形仿真的应用上具有很高的精度27。搅拌探头贯入与旋转的过程中,渣土发生极大变形,因此SPH 方法适用于该工况。1.2 渣土流变特性由坍落度试验可知,渣土在自重作用下缓慢流动塌落并最终达到某一稳定高度,表明当剪切应力大于屈服应力时,渣土似流体般发生剪切流动;当剪切应力小于屈服应力时,渣土似固体般不发生变形28,其剪应变率与剪应力之间呈非线性关系,可采用非牛顿流体模型中的 Bingham 流体模型对此进行描述,如式(6)所示。=0+p。(6)式中:0为动切应力(简称动切力或初始屈服应力),Pa;p为塑性黏度,Pas;为剪切应力,Pa;为剪切速率,s-1。2 SPH 数值模拟方案2.1 数值模型搅拌探头选用 S1、P1、B1、B2、C2 5 种型号,搅拌探头的搅拌直径均为 94 mm,如图 1 所示。渣桶底部半径 75 mm,顶部半径 100 mm,渣土在桶内的高度为145 mm,如图 2 所示。每个搅拌探头划分为约 3 000个四面体网格,单元类型均为 C3D10M;渣桶划分为5 858 个四边形网格,单元类型为 S4R;渣土划分为31 764 个六面体网格,单元类型为 C3D8R,计算过程中 C3D8R 单元转化为粒子。图 1 搅拌探头Fig.1 Mixing probe图 2 桶与渣土(单位:mm)Fig.2 Bucket and muck(unit:mm)2.2 边界条件建立几何模型时,以Y 方向作为竖向,以X、Z 方向作为水平方向。在分析过程中始终限制渣桶 X、Y、Z 3 个方向上的位移与旋转自由度,使其完全固定。对于搅拌探头,初始时限制 X、Z 方向上的位移自由度与 X、Y、Z 3 个方向上的旋转自由度,分析时限制 X、Z 方向上的位移与旋转自由度,Y 方向(竖向)的位移设定为-150 mm,旋转速度设定为60 r/min,方向为逆时针旋转。为了模拟机器人和试验台搅拌探头贯入、搅拌的工况,使用幅值曲线控制探头位移与转速的施加。搅拌探头位移与转速施加时机及变化情况如图 3 所示。在 00.5 s,搅拌探头位移沿竖向线性增加,0.5 s 达到最大位移 150 mm 后停止贯入,从而模拟搅拌探头贯入渣土的过程;在 0.51 s,搅拌探头旋转速度由0 rad/s 线性增加至 6.28 rad/s;在 14 s,搅拌探头保持 6.28 rad/s 逆时针匀速旋转,从而模拟渣土搅拌切割的过程。5141隧道建设(中英文)第 43 卷 (a)搅拌探头竖向位移(b)搅拌探头转速图 3 搅拌探头位移与转速施加时机及变化情况Fig.3 Application timing,changes of displacement,and speed of mixing probe2.3 材料参数设置搅拌探头与渣桶为钢材质,由于计算过程中搅拌探头与搅拌桶受力较小,不会发生塑性变形,所以将其视为线弹性材料。使用 Bingham 模型表达渣土的流变特性。为了得到搅拌探头在不同强度渣土中的转矩与贯入阻力的变化规律,参考文献23,29中渣土实测流变参数,将初始屈服应力分别设定为 10、50、150、250、350 Pa,将黏度设定为 100、200、300、400 Pas,并进行正交仿真试验。渣土流变参数见表 1。3 仿真结果分析3.1 搅拌转矩与贯入阻力时变曲线分析以60 r/min 转速的 S1 型搅拌探头,初始屈服应力10 Pa、黏度 300 Pas 的渣土为例,如图 4 所示,搅拌转矩与贯入阻力时变曲线可以划分为 4 个阶段。表 1 渣土流变参数23,29Table 1 Rheological parameters of soil23,29材料名称密度/(t/mm3)初始屈服应力/Pa黏度/(Pas)弹性模量/MPa泊松比渣土1.810-910,50,150,250,350100,200,300,400钢材7.810-9210 0000.3 (a)搅拌转矩随时间变化曲线(b)贯入阻力随时间变化曲线(“-”表示受压,“+”表示受拉)图 4 搅拌转矩与贯入阻力随时间变化曲线Fig.4 Curves of mixing torque and penetration resistance with time3.1.1 t1段:00.12 s搅拌探头匀速向下运动,但还未接触到渣土表面。此时搅拌探头转矩与贯入阻力皆为 0。3.1.2 t2段:0.120.5 s搅拌探头接触到渣土表面后继续匀速向下运动。如图 5(a)所示,在搅拌探头贯入过程中,受到渣土对其竖直向上的阻力 F1y,为了保持匀速向下运动,搅拌探头同时受到方向竖直向下、大小与渣土阻力相同的力,这是一对作用力与反作用力。因此,搅拌探头受压,受力为负,如图 4 所示。6141第 8 期王祥祥,等:盾构渣土改良效果评价机器人搅拌探头选型分析(a)搅拌探头竖直贯入时的渣土作用力(b)搅拌探头旋转时的渣土作用力图 5 搅拌探头受渣土作用力分析Fig.5 Analysis of force exerted on mixing probe by soil根据搅拌探头贯入阻力随时间变化曲线可知,当搅拌探头接触到渣土表面瞬间,所受贯入阻力猛然增大,然后随着贯入深度的增加逐渐增大,直到 0.5 s 瞬时,搅拌探头停止贯入时达到峰值,并迅速变为小于 0的负值。之所以出现小于 0 的负值,是由于 0.5 s 瞬时,搅拌探头还未开始转动,此时搅拌探头受到渣土竖直向上的浮力,为了克服浮力保持静止,搅拌探头同时受到竖直向下的反作用力。根据搅拌探头转矩随时间变化曲线可知,在搅拌探头向下贯入过程中受到渣土对其的反向转矩,并随着贯入深度的增加不断增大。由图 5 可知,这是由于在向下贯入时搅拌探头受到水平向右的分力 F1x,在逆时针旋转时搅拌探头受到水平向左的分力 F2x,2 种状态下搅拌探头所受的水平分力方向相反,因此在向下贯入时产生相对逆时针旋转时的反向转矩。3.1.3 t3段:0.51 s搅拌探头停止贯入,保持竖向的静止,转速开始均匀增加。搅拌探头旋转时,渣土对搅拌探头存在竖直向下的分力 F2y,为了保持竖向静止,搅拌探头同时受到竖直向上的反作用力,此时,搅拌探头受拉,受力为正。随着转速的不断增大,该反作用力不断增大最终保持稳定;随着转速的匀速增加,转矩不断增大,最终保持稳定。3.1.4 t4段:14 s搅拌探头保持匀速旋转,搅拌探头转矩与受力保持稳定。3.2 流变参数对搅拌转矩与贯入阻力影响分析以 B2 型搅拌探头为例,提取不同流变参数条件下搅拌探头的贯入阻力峰值与搅拌转矩平稳值进行分析,如图 6 和图 7 所示。由图可知,随着渣土黏度与初始屈服应力的增大,搅拌转矩与贯入阻力均呈线性增加的趋势。但是与黏度相比,初始屈服应力对搅拌转矩与贯入阻力的影响相对较小。(a)黏度对搅拌转矩的影响(b)初始屈服应力对搅拌转矩的影响图 6 流变参数对搅拌转矩的影响Fig.6 Effect of rheological parameters on mixing torque (a)黏度对贯入阻力的影响(b)初始屈服应力对贯入阻力的影响图 7 流变参数对贯入阻力的影响Fig.7 Effect of rheological parameters on penetration resistance7141隧道建设(中英文)第 43 卷3.3 搅拌探头对搅拌转矩与贯入阻力的影响分析图 8 和图 9 示出初始屈服应力 一定时,不同型号搅拌探头的搅拌转矩与贯入阻力随黏度 的变化曲线。由图可知,对于不同型号的搅拌探头,其搅拌转矩与贯入阻力均随黏度 的增大而线性增大,但是不同型号搅拌探头的搅拌转矩与贯入阻力随黏度变化曲线的斜率存在明显的差别,这说明不同的搅拌探头对于搅拌转矩与贯入阻力的“敏感度”不同,曲线的斜率越大,说明搅拌探头对于所测数据越“敏感”。对不同型号搅拌探头的搅拌转矩与贯入阻力随黏度变化曲线进行线性拟合,得到拟合曲线的斜率如表2 所示。通过对比曲线的斜率不难发现:搅拌转矩敏感度 C2B1B2P1S1,而贯入阻力敏感度 S1P1C2B1B2。(a)=10 Pa(b)=50 Pa (c)=150 Pa(d)=250 Pa(e)=350 Pa图 8 搅拌探头对搅拌转矩的影响Fig.8 Effect of mixing probe on mixing torque8141第 8 期王祥祥,等:盾构渣土改良效果评价机器人搅拌探头选型分析 (a)=10 Pa(b)=50 Pa (c)=150 Pa(d)=250 Pa(e)=350 Pa图 9 搅拌探头对贯入阻力的影响Fig.9 Effect of mixing probe on penetration resistance4 基于层次分析法的搅拌探头选型分析4.1 层次分析法层次分析法(analytic hierarchy process,AHP)是一种定性和定量相结合、系统化、层次化的分析方法。层次分析法的基本步骤包括:1)建立层次结构模型;2)构造成对比较矩阵;3)计算权向量并做一致性检验;4)计算组合权向量。对于搅拌探头选型而言,首先,对 5 个型号的搅拌探头按照搅拌转矩、贯入阻力 2 个选型标准构建层次结构模型;其次,分别构造搅拌转矩与贯入阻力的对比矩阵,计算权重向量并做一致性检验;然后,构造单因素条件下不同搅拌探头敏感性对比矩阵,9141隧道建设(中英文)第 43 卷并分别计算权向量;最后,计算组合权向量并对搅拌探头进行排名。4.2 影响搅拌探头选型的各因素权重计算4.2.1 构造相对重要性矩阵影响搅拌探头选型的主要因素是搅拌转矩与贯入阻力,但是考虑到检测时间、检测数据连续性、抵抗检测环境影响能力等因素,对于搅拌探头来说搅拌转矩的重要性要大于贯入阻力。结合表 3 和表 4,假设搅拌转矩的重要性比贯入阻力略大,通过两两对比,得到相对重要性矩阵 A。表 2 不同型号搅拌探头转矩与贯入阻力随黏度变化曲线斜率Table 2 Slope of torque and penetration resistance versus viscosity curve of different types of mixing probes搅拌探头型号初始屈服应力/Pa转矩随黏度变化曲线斜率斜率均值贯入阻力随黏度变化曲线斜率斜率均值B1B2C2P1S1102.15502.171502.102502.003501.96101.68501.751501.712501.583501.63103.63503.551503.492503.363503.14101.45501.381501.252501.203501.16101.23501.301501.022500.933501.122.0761.6703.4341.2881.1200.210.210.210.210.210.140.140.140.140.150.220.210.210.220.200.380.390.390.390.400.540.540.540.550.540.210.1420.2120.390.542表 3 风险权重评价分值表Table 3 Risk weight evaluation scores风险权重分值 aij定义1i 与 j 重要性一样大3i 比 j 重要性略大5i 比 j 重要性稍大7i 比 j 重要性大得多9i 比 j 重要性大很多2,4,6,8i 比 j 重要性比较结果处于以上结果的中间倒数 j 比 i 重要性比较结果是 i 比 j 重要性比较结果的倒数表 4 搅拌转矩与贯入阻力两两对比结果Table 4 Comparison results of mixing torque and penetration resistancei/j搅拌转矩贯入阻力搅拌转矩13贯入阻力1/31 相对重要性矩阵 A=131/31()。4.2.2 计算矩阵的特征值与特征向量通过计算矩阵 A 的最大特征值 max(A)=2;最大特征 值 对 应 的 特 征 向 量 =(0.948 683 3,0241第 8 期王祥祥,等:盾构渣土改良效果评价机器人搅拌探头选型分析0.316 227 77)T。3)计算随机一致性比率,判断矩阵是否具有满意的一致性。不一致程度的指标CI=max(A)-nn-1。(7)式中 n 为矩阵的阶数。通过计算得出 CI=2-22-1=0,所以矩阵 A 具有满意的一致性。4.2.3 特征向量的归一化将特征向量 归一化后得到 1=(0.75,0.25)T,因此得到搅拌转矩的权重值为 0.75,贯入阻力的权重值为 0.25。4.3 单因素条件下不同搅拌探头敏感性权重计算4.3.1 各搅拌探头搅拌转矩敏感性权重计算转矩随黏度变化曲线斜率是评价搅拌探头转矩敏感度的重要标准。因此,以表 5 中的斜率为两两对比的依据,构造不同搅拌探头的搅拌转矩敏感性相对重要性矩阵,结合 4.2 节中的计算步骤与方法计算归一化的特征向量。表 5 不同搅拌探头转矩随黏度变化曲线斜率Table 5 Slope of torque change curve of different mixing probes with viscosity搅拌探头型号转矩随黏度变化曲线斜率均值B12.076B21.670C23.434P11.288S11.120相对重要性矩阵 B 的最大特征值 max(B)=5,归一化的特征向量 2=(0.217,0.174,0.358,0.134,0.117)T。4.3.2 各搅拌探头贯入阻力敏感性权重计算以表 6 中的斜率为两两对比的依据,构造相对重要性矩阵 C,结合 4.2 节中的计算步骤与方法计算归一化的特征向量。表 6 不同搅拌探头贯入阻力随黏度变化曲线斜率Table 6 Slope of penetration resistance change curve of different mixing probes with viscosity搅拌探头型号贯入阻力随黏度变化曲线斜率均值B10.210B20.142C20.212P10.390S10.542 相对重要性矩阵 C 的最大特征值 max(C)=5,归一化的特征向量 3=(0.140,0.095,0.142,0.261,0.362)T。4.4 选型结果根据 4.2、4.3 节的计算结果可知:搅拌转矩与贯入阻力对影响搅拌探头选型的权向量 1=(0.75,0.25)T;各搅拌探头的搅拌转矩敏感性权向量 2=(0.217,0.174,0.358,0.134,0.117)T;各搅拌探头贯入阻力敏感性权向量 3=(0.140,0.095,0.142,0.261,0.362)T。B1 型搅拌探头的综合得分为 0.750.217+0.250.140=0.198。同理,分别求出所有型号搅拌探头的综合得分如表 7 所示。根据搅拌探头综合评分进行排序:C2B1S1P1B2。因此,C2 型搅拌探头为最佳搅拌探头。表 7 各型号搅拌探头综合得分Table 7 Comprehensive scores of various types of mixing probes搅拌探头型号 搅拌转矩(权重 0.75)贯入阻力(权重 0.25)综合得分B10.2170.1400.198B20.1740.0950.154C20.3580.1420.304P10.1340.2610.166S10.1170.3620.1785 选型结果验证5.1 试验设备及材料搅拌试验是在渣土改良标准量化试验台上完成的(见图10)。该试验台搅拌模块旋转速度为60 r/min,纵向运动模块运动速度为40 mm/s,转矩量程为4 Nm,拉压力量程为 100 N,精度为 0.3%FS,可以完成对渣土的贯入、拉拔与搅拌试验,满足本试验的测试要求。图 10 渣土改良标准量化试验台Fig.10 Standard quantitative test bench for soil conditioning试验用土取自郑州某地铁车站施工现场,取样深度位于盾构区间所在高程范围内。经过测试,该原状土小于 0.075 mm 的细颗粒占总体的 95%以上,塑限1241隧道建设(中英文)第 43 卷含水率为 15%,液限含水率为 25%。5.2 试验方案为了验证选型结果,使用易于加工制作的 C2 型搅拌探头与 B1 型搅拌探头进行试验,渣土的含水量配比与泡沫掺入比如表 8 所示。同时,通过 3.1 节的分析可知:向下贯入时搅拌探头受压,逆时针旋转时搅拌探头受拉,受力状态发生了转变,不利于力的监测。因此,在进行试验时,搅拌探头以顺时针方向旋转,使其在竖向上一直处于受压状态。表 8 试验方案Table 8 Experimental schemes试验编号搅拌探头型号含水率/%泡沫掺入比 FFIR/%N1-1C22070N1-2B12070N2-1C22570N2-2B12570N3-1C23070N3-2B130705.3 试验结果实测搅拌转矩随时间变化曲线如图 11 所示。由图可知,当渣土含水率较低时,转矩在进入渣土初期会出现峰值,随后逐渐趋于平稳;随着渣土含水率的升高,当渣土含水率为 30%时,搅拌探头进入渣土后,转矩立即进入平稳状态。实测贯入阻力时变曲线如图12 所示。由图可知,在搅拌探头向下贯入过程中贯入阻力(探头所受压力)逐渐增大,直到探头停止向下移动,所测贯入阻力达到最大值;当探头开始旋转搅拌渣土时,探头所受压力趋于稳定。图 11 实测搅拌转矩随时间变化曲线Fig.11 Changing curves of measured mixing torque with time在低围压条件下,渣土黏度随着含水率的升高而逐渐降低,且呈近似线性变化的趋势30。因此,取每组试验搅拌转矩平稳值与贯入阻力峰值分别绘制其随含水率的变化曲线,如图 13 和图 14 所示。可以看出,2 种搅拌探头数据变化规律与数值模拟结果(见图 8、图 9)基本一致,可以说明选型结果的合理性。图 12 实测贯入阻力时变曲线(正值表示受压)Fig.12 Changing curves of measured penetration resistance with time(under pressurized condition)图 13 不同搅拌探头搅拌转矩随含水率变化趋势Fig.13 Variation trend of mixing torque of different mixing probes with water content图 14 不同搅拌探头贯入阻力随含水率变化趋势Fig.14 Variation trend of penetration resistance of different mixing probes with water content2241第 8 期王祥祥,等:盾构渣土改良效果评价机器人搅拌探头选型分析6 结论与建议搅拌探头作为机器人和试验台上唯一与渣土发生接触的部件,其对搅拌转矩与贯入阻力检测的准确性与灵敏度至关重要。因此,本文将渣土假定为宾汉姆流体,通过 SPH 方法对不同型号的搅拌探头贯入、搅拌渣土的过程进行动态仿真,运用层次分析法综合考虑不同搅拌探头对贯入阻力与搅拌转矩的检测灵敏度,得到如下结论:1)根据搅拌转矩与贯入阻力时变曲线分析可知,贯入阻力随搅拌探头匀速贯入深度的增加而增加,搅拌转矩随着搅拌探头转速的增加而增加。2)根据贯入阻力与搅拌转矩的影响因素分析可知,贯入阻力与搅拌转矩均随渣土初始屈服应力与塑性黏度的增大而增大;并且,渣土的塑性黏度对其影响更大,对初始剪切屈服应力的影响较小。3)C2 型搅拌探头可以同时兼顾贯入阻力与搅拌转矩灵敏度要求,为本次选型的最佳搅拌探头,可用于粉质黏土地层渣土改良效果评价机器人掘进试验。进一步地,结合本文搅拌探头的选型结果,在标准量化试验台与机器人上使用相同的搅拌探头与渣土容器,使其在相同的检测体系下分别开展室内渣土改良试验与现场搭载试验,用于探索改良规律并建立评价体系;在此基础上,结合机器人现场检测数据与掘进参数评价当前渣土改良状态,并给出改良建议。参考文献(R Re ef fe er re en nc ce es s):1 洪开荣,杜彦良,陈馈,等.中国全断面隧道掘进机发展历程、成就及展望J.隧道建设(中英文),2022,42(5):739.HONG Kairong,DU Yanliang,CHEN Kui,et al.Full-face tunnel boring machines(shields/TBMs)in China:History,achievements,and prospectsJ.Tunnel Construction,2022,42(5):739.2 陈馈,谭顺辉,王江卡,等.盾构施工关键技术M.北京:中国铁道出版社,2020.CHEN Kui,TAN Shunhui,WANG Jiangka,et al.Key technologies in tunnelling by shield machineM.Beijing:China Railway Publishing House,2020.3 MILLIGAN G.Lubrication and soil conditioning in tunnelling,pipe jacking and microtunneling:A state-of-the-art review R.London:Geotechnical Consulting Group,2000.4 MESSERKLINGER S,ZUMSTEG R,PUZRIN A M.A new pressurized vane shear apparatus J.Geotechnical Testing Journal,2011,34(2):112.5 魏康林.土压平衡式盾构施工中“理想状态土体”的探讨J.城市轨道交通研究,2007(1):67.WEI Kanglin.Discussion on the Ideal Soil in the earth pressure balanced shield tunnelling J.Urban Mass Transit,2007(1):67.6 VINAI R,OGGERI C,PEILA D.Soil conditioning of sand for EPB applications:A laboratory research J.Tunnelling and Underground Space Technology,2008,23(3):308.7 MAIDL B,HERRENKNECHT M,MAIDL U,et al.Mechanised shield tunnelling M.2nd ed.Berlin:Ernstand Sohn,2012.8 王树英,刘朋飞,胡钦鑫,等.盾构隧道渣土改良理论与技术研究综述J.中国公路学报,2020,33(5):8.WANG Shuying,LIU Pengfei,HU Qinxin,et al.State-of-the-art on theories and tech