复合地层下盾构施工模态系统刀盘荷载及振动特性研究王鲁琦1,2， 刘海宁1,2,*， 李 苗2， 周建军1， 刘汉东2 (1. 盾构及掘进技术国家重点实验室， 河南 郑州 450001；2. 华北水利水电大学岩土工程与水工结构研究院， 河南 郑州 450045) 摘要：基于苏埃通道盾构工程和盾构及掘进技术国家重点实验室的盾构施工模态系统的刀盘实体，利用JHC本构模型对岩体进行仿真，通过显示动力学软件LS-DYNA对不同模式下的滚刀群切削复合地层进行数值模拟，并结合滚刀、刀盘的力学分析及傅立叶变换函数，得到滚刀和刀盘的三向荷载-时间历程曲线及荷载的频谱图。通过分析荷载曲线及频谱图可知： 1)在荷载特性方面，切削过程中滚刀和刀盘荷载均具有突变性、随机性和冲击性；当滚刀运动经过软硬岩层交界面时，滚刀切点荷载发生突变，有可能发生剧烈震荡。2)在振动特性方面，刀盘荷载的频率主要集中在0～5 Hz，属于典型的低频振动，并与滚刀的频率分布一致； 软岩及硬岩仅对滚刀荷载的幅值产生影响，对其频率分布影响有限。结合刀盘在软硬不均地层条件下静力学的分析结果，确定刀盘在不均匀受力状态下的最大应力、应变位置，为苏埃通道滚刀破岩物理模型试验刀盘测试系统的构建提供了依据。 关键词：复合地层； 盾构； 模型试验； 刀盘荷载； 刀盘振动； LS-DYNA 0 引言当盾构施工穿越软硬不均地层时，刀盘刀具将承受冲击荷载，极易导致刀具产生严重磨损和异常损坏，对施工安全和工程进展造成严重的影响[1-3]。盾构造价昂贵且移动困难，在施工现场进行刀盘刀具荷载、振动的监测以及相关试验具有较大的难度，且复杂的现场施工条件对于实时数据的监测也会产生较大的影响； 同时，由于施工现场条件的局限性，无法对盾构的各项参数进行试验性调整； 因而，在室内开展盾构施工模型试验具有较高的实际意义和理论价值[4-5]。 针对隧道下穿软硬不均地层所采取的处治措施，国内外还没有形成统一的认知和定论，更没有形成统一思路。竺维彬等[6]讨论了复合地层(孤石、基岩侵入体等)高强度岩体的分类、特性以及盾构施工存在的风险，提出采用地面引孔下药施爆的控制钻爆法，严格控制破碎炸块体量以保证盾构的正常掘进。李明阳等[7]基于Terzaghi松动土压力模型与CSM模型，对广州地铁3号线大石—汉溪区间不同风化程度下混合岩地层的土压平衡盾构掘进过程进行了模拟分析，验证了综合考虑软土和硬岩的复合地层掘进参数模拟分析方法的有效性。已有研究多偏重于施工技术及掘进参数的选定，从刀盘荷载及振动特性角度进行研究和分析还比较少。 为了探究盾构穿越软硬不均地层时刀盘刀具的荷载及振动特性，结合汕头市苏埃跨海通道的实际工况，通过盾构及掘进技术国家重点实验室的盾构掘进模态综合试验台进行模型试验研究。该盾构掘进模态综合试验台主要由机械结构、液压泵站、掘进装置、旋转装置、螺旋输送装置和控制系统等构成，可在垂直和水平状态下，通过设置不同的刀具材质、刀间距、破岩切削速度和进给量进行试验，进而统计分析出不同刀具参数下的破岩效率和刀具寿命，为特定的工程配置相应材质的刀具并提出建议性刀具布置方案； 但是，该系统对于振动特性和刀盘刀具荷载的实时监测功能有所欠缺，无法为苏埃通道滚刀破岩物理模型试验提供必要的数据支撑。为此，在苏埃通道模型材料数据的基础上，通过对该盾构施工模态系统进行数值分析，得到其刀盘刀具的荷载、振动特性，进而为盾构施工模态试验台刀盘测试系统的构建提供理论依据。 1 软硬不均地层岩体本构模型的建立为了得到盾构施工模态系统中刀盘刀具的实时荷载历程曲线，本文采用有限元软件LS-DYNA进行数值模拟。盾构破岩是一个复杂的动态过程，是具有几何、材料以及接触等多种非线性特性的复合系统[8]。LS-DYNA软件采用的显式算法能够有效降低处理大规模非线性动力学问题时出现计算不收敛的概率，在接触问题上，其采用的对称罚函数法适用于滚刀破岩仿真时滚刀与岩体之间离散方程中时间积分的计算[8-11]。 在分析汕头市苏埃过海隧道工程地质条件的基础上，优选出一段软硬不均地层作为室内掘进物理试验的研究对象，所选典型断面的软、硬岩层之比为1∶1。根据所选断面的工程地质特征，概化该段隧道的工程地质力学模型，确定微风化花岗岩和软弱地层(中粗砂)的物理力学参数，见表1。 表1 典型岩土层的物理力学参数 Table 1 Physico-mechanical parameters of typical strata  围岩类别抗压强度Rc/MPa黏聚力c/MPa内摩擦角φ/(°)弹性模量E/MPa泊松比μ密度ρ/(kg/m3)微风化花岗岩2102162．57．1×1050．202760中粗砂28220．252000 根据Froude比例定律的要求，进行长度、应力、力和时间相似比尺设计[12-13]。 1)长度比尺的确定。汕头苏埃过海隧道施工拟采用开挖直径(lP)为14.93 m的复合式泥水盾构，盾构施工模态试验系统的切削直径(lM)为2.28 m，则长度比尺
(1) 2)应力比尺 kσ=kρ·kg·kl=kγ·kl。 (2)式中： Kρ为密度比尺； Kg为加速度比尺； Kγ为重度比尺。 本次物理模型试验主要模拟软硬不均地层盾构掘进中的受力、变形性能及其破岩效率，则岩土体的力学性能是重要的模拟参数，而对其密度模拟相似要求不高。其中，kγ可为任意值，所以，原型和模型的应力比尺、几何比尺可独立选取。根据试验精度要求和试验工作量的大小，初步确定本次试验的应力比尺kσ=1/5。 3)力比尺
(3) 4)扭矩比尺
(4) 5)时间比尺
(5) 根据相似要求，可得盾构地层物理力学参数如表2所示，其中软岩的相似材料采用淮南矿区动载模型试验的塑性混凝土。 表2 盾构地层物理力学模拟参数 Table 2 Physico-mechanical parameters of shield boring strata  围岩类别抗压强度Rc/MPa黏聚力c/MPa内摩擦角φ/(°)弹性模量E/MPa泊松比μ密度ρ/(kg/m3)硬岩相似材料424．262．51．42×1050．202400软岩相似材料0．8284．40．252000 采用LS-DYNA软件中的JHC动态损伤模型对岩体进行仿真。该本构模型在考虑损伤力学对岩石破坏演化影响的基础上，引进了损伤因子，并基于岩石损伤演化过程对岩体进行力学概化，能够有效地模拟滚刀破碎侵彻岩体时岩石的断裂破坏过程[14-16]。 JHC本构模型共涉及21个参数，可分为3类： 1)强度参数，包括抗压强度fc、抗拉强度T、粘结强度A、压力硬化参数B、压力硬化指数N、应变率系数C、标准化最大强度Smax、参考应变率 由岩石常规试验数据可以得到参数fc、T、E、μ、plock和ρ的取值。根据国内学者的相关研究[17-18]，C值与岩石体的初始抗压强度无关，一般取0.007； A、B、N、Smax的参数值只与岩石体静态抗压强度有关，参考文献进行取值； D2=1.0，εfmin=0.01， μlock=0.1； 其他参数 D1=0.01/(1/6+T)； (6)
(7) pcrush=fc/3。 (8) 根据上述分析，由盾构地层物理模拟参数得到软、硬岩的JHC模型参数如表3和表4所示。 2 滚刀切削模型的建立盾构施工模态系统采用17 in(432 mm)标准滚刀，滚刀由刀轴、刀体及刀圈构成。本文模拟的重点在于获得滚刀刀圈破岩过程中的荷载及振动特性，因而，将滚刀简化为中空薄圈，只保留刀圈及刀体[8,19]，其中，外侧为刀圈，内侧为刀体，如图1所示。 刀圈和刀体均设定为线弹性、各向同性，密度为 7 800 kg/m3，弹性模量并为2.10×1011 Pa，泊松比为0.3，且刀体设定为刚体，并约束其除切削方向以外的位移。 表3 软岩JHC本构模型参数 Table 3 Parameters of JHC model of soft rock  剪切模量G/Pa粘结强度A压力硬化参数B应变率系数C抗压强度fc/Pa抗拉强度T/Pa参考应变率ε·0/s-1最小塑性应变εfmin压碎压力pcrush/Pa压碎体积应变μcrush1．76×1060．791．60．0078．00×1058．00×10410．012．67×1050．911压实压力plock/Pa压密体积应变μlock损伤系数D2压力系数k1/Pa压力系数k2/Pa压力系数k3/Pa压力硬化指数N标准化最大强度Smax损伤系数D1失效类型fs8．00×1080．118．50×1010-1．71×10112．08×10110．6170．37495-1 表4 硬岩JHC本构模型参数 Table 4 Parameters of JHC model of hard rock  剪切模量G/Pa粘结强度A压力硬化参数B应变率系数C抗压强度fc/Pa抗拉强度T/Pa参考应变率ε·0/s-1最小塑性应变εfmin压碎压力pcrush/Pa压碎体积应变μcrush5．92×10100．791．60．0074．20×1074．20×10610．011．40×1071．775×10-4压实压力plock/Pa压密体积应变μlock损伤系数D2压力系数k1/Pa压力系数k2/Pa压力系数k3/Pa压力硬化指数N标准化最大强度Smax损伤系数D1失效类型fs8．00×1060．118．50×1010-1．71×10112．08×10110．6170．37495-1 盾构及掘进技术国家重点实验室的盾构施工模态系统内置滚刀14把，包括6把中心滚刀和8把正滚刀。中心滚刀为6把单刃滚刀楔装到一个刀轴上，正滚刀分布在中心滚刀周围，一共形成14条滚刀轨迹线。刀盘CAD图及ANSYS三维示意图如图2所示。 为了简化仿真模型，进而缩短仿真时间，将该试验台的刀盘切削模型简化为2个基本模型： 1)中心滚刀群切削模型，采用6把滚刀切削圆盘状岩盘； 2)正滚刀切削模型，采用7#、8#、9#、10#正滚刀进行模拟，这4把正滚刀切削半径较小，相对位置具有代表性，所切削岩盘取圆环状，在保证滚刀轨迹均能够实时监测的同时，可减少模型计算量。滚刀切削模型如图3所示。  (a) 标准滚刀 (b) 标准滚刀简化模型 图1 标准滚刀结构示意图 Fig. 1 Structure diagrams of standard disc cutter (a) CAD图 (b) 三维示意图 图2 刀盘示意图 Fig. 2 Sketch diagrams of cutterhead (a) 中心滚刀 (b) 正滚刀 图3 滚刀切削模型 Fig. 3 Cutting models of disc cutter 以中心滚刀切削模型为例，简述其有限元仿真的关键点。在剖分网格方面，刀圈和刀体为自由剖分，岩体为扫略(sweep)剖分，网格大小均设定为20 mm。在接触方面，刀圈和刀体设定为自动单面(single surface)接触，保证刀圈和刀体之间不发生分离； 刀圈和岩体设定为面面侵蚀(eroding surface)接触，在岩体单元失效删除后，该接触能够保证滚刀刀圈与岩体单元继续进行侵蚀模拟。此外，通过无反射边界条件的设定，将模型中的岩体仿真为无限大，其底面和侧面各施加30 MPa和15 MPa的围压，使之符合现实工况中的岩体状态。 根据典型软硬不均地层盾构隧道的工程地质特征和地质力学模型，确定盾构的推力和转速。初步确定刀盘转速为5 r/min，并将该速度以角速度0.63 rad/s的形式施加到岩盘上； 掘进速度为5 mm/s，施加到刀体上。中心滚刀刀轴线的初始位置与软硬岩层分界线相交，以便研究滚刀经过软硬不均地层交界面时滚刀荷载的变化，其初始状态的有限元模型如图4所示。 图4 软硬不均地层中心滚刀初始状态 Fig. 4 Initial state of center disc cutter in heterogeneous ground 3 滚刀的荷载及振动特性利用LS-DYNA软件的求解器对有限元仿真模型进行求解，可以得到三轴坐标系下滚刀与岩体接触点的荷载-时间历程曲线。中心滚刀及正滚刀在运动过程中的受力分析如图5所示(图中θ为极角，Fx、Fy、Fz为全局坐标系下的荷载，Fr、Fs、Fv分别为滚刀的滚动力、侧向力和垂直力)。 (a) 中心滚刀 (b) 正滚刀 图5 滚刀受力分析[11] Fig. 5 Force analysis of disc cutter[11] 根据滚刀的受力分析可知，中心滚刀的三向荷载与三轴坐标荷载之间的关系： Fr=-Fxsin θ-Fzcos θ； (9) Fs=Fzsin θ-Fxcos θ； (10) Fv=Fy。 (11) 正滚刀的三向荷载与三轴坐标荷载之间的关系： Fr=Fxsin θ+Fzcos θ； (12) Fs=-Fzsin θ+Fxcos θ； (13) Fv=Fy。 (14) 在仿真模型计算结果的基础上，结合滚刀三向荷载计算公式(9)—(14)，参考工程的实测荷载谱对数据进行修正，得到盾构施工模态系统在软硬不均地层条件下中心滚刀及正滚刀的荷载-时间历程曲线，如图6和图7所示。 (a) 侧向力 (b) 垂直力 (c) 滚动力 图6 中心滚刀三向荷载-时间历程曲线 Fig. 6 Time-dependent three-dimensional load curves of center disc cutter 通过MATLAB软件调用傅立叶计算函数对滚刀的三向荷载-时间历程数据进行处理，得到滚刀切削过程的频谱图。滚刀切削硬岩的三向荷载频谱图如图8所示，软岩的荷载频谱图除幅值不同以外，其频率分布与硬岩相同。 对荷载-时间历程曲线及荷载频谱图进行数据统计，见表5。 通过对滚刀三向荷载-时间历程曲线进行分析，得出软硬不均地层条件下滚刀在切削过程中的荷载特性。 (a) 侧向力 (b) 垂直力 (c) 滚动力 图7 正滚刀三向荷载-时间历程曲线 Fig. 7 Time-dependent three-dimensional load curves of front disc cutter (a) 中心滚刀 (b) 正滚刀 图8 滚刀三向荷载频谱图 Fig. 8 Spectrographs of three-dimensional load of disc cutter 1)滚刀在切削过程中，滚刀切点处的三向荷载均具有突变性、随机性和冲击性。 2)当滚刀仅在软岩或硬岩中切削时，滚刀的三向荷载-时间历程曲线呈周期性变化，且该周期与滚刀转速有关，中心滚刀以π/2为一个周期(3 s)，正滚刀以π/4为一个周期(1.5 s)。 3)滚刀从软岩到硬岩以及从硬岩到软岩的运动过程中，滚刀切点荷载会发生突变，在垂直力方向存在最大荷载差，此时滚刀将会产生较大程度的振动，有可能会加速滚刀的磨损； 滚刀荷载曲线在软岩和硬岩中产生的周期性变化相似，由于荷载值跨度较大，在图6和图7中未能清晰显示。 表5 软硬不均地层滚刀荷载及频谱统计 Table 5 Statistics of disc cutter loads and frequency spectrum in heterogeneous ground 岩石类别中心滚刀侧向力均值/N垂直力均值/N滚动力均值/N频率/Hz正滚刀侧向力均值/N垂直力均值/N滚动力均值/N频率/Hz软岩2．50×1022．48×1039．05×1020～52．50×1022．48×1039．05×1020～5，120～130硬岩1．25×1041．24×1051．87×1040～51．24×1041．25×1051．87×1040～5，120～130 4)虽然正滚刀和中心滚刀的时间历程曲线振动形式不同，但是两者的侧向力、垂直力及滚动力的均值接近。 对滚刀三向荷载频谱图进行分析，中心滚刀频率集中在0～5 Hz，正滚刀频率主要分布在0～5 Hz，少部分分布在120～130 Hz； 滚刀在软硬不均地层交界线处运动时，荷载幅值在25 N和0.5 N产生交替，发生剧烈变化，对滚刀产生的振动影响最大。 4 刀盘荷载及振动特性忽略滚刀荷载传递过程中的能量损失，并且不考虑滚刀切削过程中惯性力和惯性力矩对刀盘的影响，在修正滚刀三向荷载数据的基础上，将极限荷载设定为滚刀垂直力，并使滚刀的滚动力均值、垂直力均值及侧向力均值满足Fr=0.15Fv、Fs=0.1Fv，可推导得到刀盘的荷载状态[20]。 盾构开挖方向上刀盘的轴向力
(15) 沿刀盘面水平轴线的刀盘径向力 Fx=Fxi+Fxj。 (16) 沿刀盘面竖直轴线的刀盘径向力 Fy=Fyi+Fyj。 (17) 由垂直力对盘面产生的刀盘倾覆力矩
(18)
(19) 由中心滚刀及正滚刀绕刀盘中心轴线的力矩之和产生的刀盘扭矩
(20) 式(15)—(20)中： n为中心滚刀数量，即n=6； m为正滚刀数量，即m=8； Fvi、Fxi、Fyi、Fri分别为第i把中心滚刀的垂直力、横向不平衡力、纵向不平衡力和滚动力，N； Fvj、Fxj、Fyj、Frj分别为第j把正滚刀的垂直力、横向不平衡力、纵向不平衡力和滚动力，N； lix、liy分别为中心滚动垂直力到x轴、y轴的距离，m； ljx、ljy分别为正滚刀垂直力到x轴、y轴的距离，m； di、dj分别为中心滚刀、正滚刀到刀盘中心的距离，m。 由式(15)—(20)，结合修正后的三向荷载数据，得到盾构施工模态系统的刀盘荷载-时间历程曲线，该曲线为滚刀在均质地层(仅为硬岩条件)下的刀盘荷载状态，如图9所示。 采用MATLAB软件对刀盘荷载数据进行傅立叶变换，得到刀盘频谱图，如图10所示。
(a) 刀盘横向不平衡力
(b) 刀盘纵向不平衡力
(c) 刀盘轴向力
(d) 刀盘纵向倾覆力矩
(e) 刀盘横向倾覆力矩
(f) 刀盘扭矩 图9 刀盘荷载-时间历程曲线 Fig. 9 Time-dependent load curves of cutterhead (a) 刀盘横向不平衡力 (b) 刀盘纵向不平衡力 (c) 刀盘轴向力 (d) 刀盘纵向倾覆力矩 (e) 刀盘横向倾覆力矩 (f) 刀盘扭矩 图10 刀盘荷载频谱图 Fig. 10 Spectrographs of load of cutterhead 由图9和图10分析可知： 刀盘荷载同样具有突变性、随机性和冲击性的特点，刀盘荷载的频率主要集中在0～5 Hz，属于低频振动，与滚刀的主要频率分布一致； 刀盘荷载在120 Hz和200 Hz处有少量分布，极少量分布在300 Hz和400 Hz处； 对比滚刀和刀盘的荷载频谱图可知，刀盘荷载的频率基本由滚刀荷载频率组成。 盾构施工模态系统在进行软硬不均地层的掘进物理模型试验时，滚刀处于自转及公转状态，切削过程中受力点及滚刀空间位置呈螺旋线变化，很难利用传感器或应变片对滚刀荷载和频率进行实时监测。通过在刀盘上布设加速度传感器，对刀盘的三轴加速度进行监测，可以推断得到滚刀实时的荷载、振动特性，为滚刀破岩机制的研究提供数据支撑。通过有限元仿真模拟可知，盾构刀盘在掘进过程中的振动频率主要集中在0～5 Hz、120 Hz，为加速度传感器的选用提供了依据。 此外，通过ANSYS的静力学分析软件，对刀盘在软硬不均地层条件下的应力、应变分布进行了分析，确定了刀盘在不均匀受力状态下的最大应力、应变位置，得到传感器、应变片的最佳布置设点，为刀盘测试系统的构建提供依据，如图11所示。
(a) 刀盘应力状态 (b) 刀盘应变状态 (c) 传感器及应变片拟布设位置 图11 刀盘测试系统静力学分析 Fig. 11 Static analysis of cutterhead test system 5 结论与讨论基于苏埃跨海隧道的工程实例及盾构及掘进国家重点实验室的盾构施工模态系统的刀盘实体，通过LS-DYNA软件对滚刀切削不均匀地层进行数值仿真，得到不均匀地层条件下滚刀及刀盘的荷载及振动特性： 1)在切削过程中，中心滚刀及正滚刀切点处的三向荷载、刀盘荷载均具有突变性、随机性和冲击性，其中，正滚刀与中心滚刀的时间历程曲线振动形式不同，但三向荷载均值接近。 2)当滚刀仅在软岩或硬岩中切削时，滚刀及刀盘的荷载-时间历程曲线呈周期性变化，且该周期与刀盘转速有关。 3)滚刀切削过程中，软岩及硬岩等地层条件对滚刀荷载的频率分布影响有限，但对其幅值影响较大。滚刀在软硬岩的交界面运动时，滚刀切点荷载会发生突变，其中，垂直力方向上存在最大荷载差； 荷载频谱图的幅值也产生剧烈变化，此时滚刀可能产生较大程度的振动，会加速滚刀的磨损。 4)刀盘荷载的频率全部由滚刀荷载频率组成，主要集中在0～5 Hz，属于典型的低频振动，与滚刀的主要频率分布具有一致性，因此，可以在刀盘上布设传感器和应变片，对刀盘的荷载及振动特性进行实时监测，进而推导得到滚刀的实时状态。 基于上述结论，结合静力学分析结果所确定的刀盘不均匀受力状态下的最大应力、应变位置，为苏埃通道滚刀破岩物理模型试验刀盘荷载、振动实时监测系统的构建提供了依据。 随着后续工作的开展，基于刀盘实时监测系统及相似材料的研究成果，通过物理模型试验，可获得表征刀盘振动及荷载特性的实测数据。结合本文的理论推算值可对刀盘系统的本构模型进行修正，并进行数据演算结果的反推，确定该工况下的最优掘进参数，进而指导苏埃通道盾构的实际施工。 参考文献(References)： [1] 杨书江.盾构在硬岩及软硬不均地层施工技术研究[D].上海: 上海交通大学，2006.(YANG Shujiang. Research on tunneling technology in hard and mix ground by shield machine[D].Shanghai: Shanghai Jiaotong University, 2006. (in Chinese)) [2] 谭忠盛，洪开荣，万姜林，等.软硬不均地层盾构的研究及掘进技术[J].岩石力学与工程学报,2006，25(增刊2)： 3945-3952.(TAN Zhongsheng，HONG Kairong，WAN Jianglin，et al. Study of composite shield and construction technique in complex uneven strata[J]. Chinese Journal of Rock Mechanics and Engineering, 2006， 25(S2)： 3945-3952. (in Chinese)) [3] 宋天田，周顺华. 复合地层条件下盾构刀盘设计研究[J].地下空间与工程学报, 2007，3(3)： 479-482.(SONG Tiantian, ZHOU Shunhua. Design analysis of cutterhead of shield machine in mixed strata[J]. Chinese Journal of Underground Space and Engineering, 2007，3(3)： 479-482. (in Chinese)) [4] 施虎，杨华勇，龚国芳，等.盾构掘进机关键技术及模拟试验台现状与展望[J].浙江大学学报(工学版)，2013，47(5)： 741-749.(SHI Hu, YANG Huayong, GONG Guofang, et al. Key technologies of shield tunneling machine and present status and prospect of test rigs for tunneling simulation[J].Journal of Zhejiang University (Engineering Science), 2013， 47(5)： 741-749. (in Chinese)) [5] 姚乐.复合地层土压平衡盾构掘进适应性评价[D].北京： 北京交通大学，2015.(YAO Le. Research on adaptability evaluation for EPB shield tunneling in composite strata[D].Beijing: Beijing Jiaotong University, 2015. (in Chinese)) [6] 竺维彬，黄威然，孟庆彪，等. 盾构工程孤石及基岩侵入体爆破技术研究[J].现代隧道技术，2011，48(5)： 12-17.(ZHU Weibin, HUANG Weiran, MENG Qingbiao, et al. Research on controlled blasting technology for boulders and invading bedrock in shield construction[J]. Modern Tunnelling Technology, 2011， 48(5)： 12-17. (in Chinese)) [7] 李明阳，杨海涛，邹高明. 复合地层土压平衡盾构掘进参数模拟分析研究[J].隧道建设, 2012， 32(3)： 287-295.(LI Mingyang, YANG Haitao, ZOU Gaoming. Simulation analysis of boring parameters of EPB shield in complex strata[J]. Tunnel Construction, 2012, 32(3)： 287-295. (in Chinese)) [8] 吴起星.复合地层中盾构机滚刀破岩力学分析[D].广州： 暨南大学, 2011.(WU Qixing. Mechanical analysis of rock breaking by disc cutter during shield tunneling in mixed face ground conditions[D].Guangzhou: Jinan University, 2011.(in Chinese)) [9] 张魁，夏毅敏，谭青，等.不同围压条件下TBM刀具破岩模式的数值研究[J].岩土工程学报, 2010, 32(11)： 1780-1787.(ZHANG Kui，XIA Yimin，TAN Qing，et al．Numerical study of modes of breaking rock by TBM cutter under different confining pressures[J]. Chinese Journal of Geotechnical Engineering, 2010, 32(11)： 1780-1787. (in Chinese)) [10] 谭青，屠昌锋，暨智勇，等.盘形滚刀切割岩石的数值仿真[J].现代制造工程， 2009(2): 62-64.(TAN Qing，TU Changfeng，JI Zhiyong，et al．Numerical simulation of rock cutting by disc cutter[J].Modern Manufacturing Engineering， 2009(2): 62-64. (in Chinese)) [11] 杨静.TBM刀盘系统振动特性分析与掘进现场测试研究[D].大连： 大连理工大学，2014.(YANG Jing. Vibration characteristics analysis of TBM cutter-head system and research on tunneling testing[D].Dalian： Dalian University of Technology, 2014. (in Chinese)) [12] 王伟，李小春.无临空面深层岩体爆破增渗模型试验相似律研究及相似条件的实现[J].岩土力学, 2009, 30(7)： 1892-1898.(WANG Wei, LI Xiaochun. Study of similarity law of model test of blasting enhanced permeability in deep rock mass without free surface and its implementation[J].Rock and Soil Mechanics, 2009, 30(7): 1892-1898. (in Chinese)) [13] 袁宗盼，陈新民，袁媛，等. 地质力学模型相似材料配比的正交试验研究[J]. 防灾减灾工程学报， 2014, 34(2)： 197-202.(YUAN Zongpan, CHEN Xinmin, YUAN Yuan, et al. Study of ratio of similar material of geomechanical model based on orthogonal test[J].Journal of Disaster Prevention and Mitigation Engineering, 2014, 34(2)： 197-202. (in Chinese)) [14] TU Zhenguo，LU Yong. Evaluation of typical concrete material models used in hydrocodes for high dynamic response simulations[J]. International Journal of Impact Engineering， 2009, 36(1): 132-146. [15] 高红，郑颖人，冯夏庭.材料屈服与破坏的探索[J].岩石力学与工程学报，2006，25(12)： 2515-2522.(GAO Hong， ZHENG Yingren， FENG Xiating．Exploration on yield and failure of materials[J].Chinese Journal of Rock Mechanics and Engineering， 2006， 25(12): 2515-2522. (in Chinese)) [16] 张东风，陈星明，肖正学.脆性岩石中爆炸扩腔作用的数值模拟及分析[J].爆破，2014,31(1)： 60-64.(ZHANG Dongfeng, CHEN Xingming, XIAO Zhengxue. Numerical simulation of bulling blasting in brittle rock[J].Blasting, 2014, 31(1)： 60-64. (in Chinese)) [17] 张凤国，李恩征.大应变、高应变率及高压强条件下混凝土的计算模型[J].爆炸与冲击, 2002, 22(3): 198-202.(ZHANG Fengguo, LI Enzheng. A computational model for concrete subjected to large strains, high strain rates, and high pressures[J].Explosion and Shock Waves, 2002, 22(3): 198-202. (in Chinese)) [18] 熊益波，胡永乐，徐进，等.混凝土Johnson-Holmquist模型极限面参数确定[J].兵工学报, 2010, 31(6): 746-751.(XIONG Yibo, HU Yongle, XU Jin, et al. Determining failure surface parameters of the Johnson-Honlmquist concrete constitutive model[J].Acta Armamentarii, 2010, 31(6): 746-751. (in Chinese)) [19] 霍军周，杨静，孙伟，等.不同模式下 TBM 刀群三维回转切削仿真与优化设计[J].哈尔滨工程大学学报, 2014, 35(11)： 1403-1408.(HUO Junzhou，YANG Jing，SUN Wei，et al. Simulation and optimization design of three-dimensional rotating cutting action of TBM cutter group with different modes[J].Journal of Harbin Engineering University, 2014, 35(11)： 1403-1408. (in Chinese)) [20] ZHANG Z X, KOU S Q, Lindqvist P A. In-situ measurements of cutter forces on boring machine at sp hard rock laboratory[J].Rock Mechanics and Rock Engineering, 2003, 36(1): 63-83. Study of Vibration and Loads Characteristics of Modal Cutterhead of Shield Boring in Composite Ground WANG Luqi1, 2, LIU Haining1, 2, *, LI Miao2, ZHOU Jianjun1, LIU Handong2 (1.State Key Laboratory of Shield Machine and Boring Technology, Zhengzhou 450001, Henan, China; 2. Institute of Geotechnical Engineering and Hydraulic Structure Engineering,North China University of Water Resources and Electric Power, Zhengzhou 450045, Henan, China) Abstract：The rock is simulated by JHC constitutive model, the rock broken by shield cutter group under different modes are numerically stimulated by dynamic analysis software LS-DYNA; and then the time-dependent three-dimensional load and load spectrogram of disc cutter and cutterhead of shield are derived by mechanical analysis of disc cutter and cutterhead and Fourier transform, by taking Su’ai shield tunnel and cutterhead of TBM under modal system for examples. The analytical results show that: 1) The loads of disc cutter and the cutterhead show mutagenicity，randomness and impact during rock breaking; the load of disc cutter varies suddenly and the disc cutter vibrates seriously when the disc cutter cutting the interface of hard rock and soft rock. 2) The loading frequency of the cutterhead, 0-5 Hz and belongs to low frequency, coincides with that of disc cutter; the soft rock and hard rock only affect the amplitude of disc cutter load, and affect the distribution of the frequency little. Furthermore, the positions of maximum stress and strain of cutterhead under heterogeneous stress condition are determined by statics analysis of cutterhead in heterogeneous strata. The results can provide reference for cutterhead test system of physical modeling experiment of rock broken by disc cutter of Su’ai Tunnel. Keywords：composite strata; shield; model test; load on cutterhead; cutterhead vibration; LS-DYNA 收稿日期：2016-03-21; 修回日期：2016-12-15 基金项目：国家高技术研究发展计划(“863”计划)专项(2012AA041802)； 盾构及掘进技术国家重点实验室开放基金(SMBT2012K002)； 福建省地质灾害重点实验室开放基金(FJKLGP2012K004)； 华北水利水电大学研究生教育创新计划基金(YK2015-12) 第一作者简介：王鲁琦(1993—)，男，山东菏泽人，华北水利水电大学地质工程专业在读硕士，研究方向为工程地质及岩土工程。E-mail： 377030692@qq.com。*通讯作者： 刘海宁， E-mail： lhngeo@yahoo.com.cn。 DOI：10.3973/j.issn.1672-741X.2017.01.017 中图分类号：U 455.3 文献标志码：A 文章编号：1672-741X(2017)01-0103-10 |
|
|
。
。
。
。
；
