CN108388689B - 一种基于量纲分析的盾构刀盘扭矩计算方法 - Google Patents
一种基于量纲分析的盾构刀盘扭矩计算方法 Download PDFInfo
- Publication number
- CN108388689B CN108388689B CN201810034645.5A CN201810034645A CN108388689B CN 108388689 B CN108388689 B CN 108388689B CN 201810034645 A CN201810034645 A CN 201810034645A CN 108388689 B CN108388689 B CN 108388689B
- Authority
- CN
- China
- Prior art keywords
- unit
- torque
- cutter head
- main control
- cutter
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Excavating Of Shafts Or Tunnels (AREA)
Abstract
一种基于量纲分析的盾构刀盘扭矩计算方法,包括:确定影响刀盘扭矩的参数;根据所述的影响刀盘扭矩的参数分别计算主控项及无量纲化刀盘扭矩大小;将计算得到的主控项及无量纲化刀盘扭矩进行汇总;利用利用线性回归方法,得到主控项的系数识别;根据主控项系数识别结果,得到刀盘扭矩的计算模型。本发明其模型预测结果可为装备的智能化设计提供理论依据,为不同地质的安全施工提供参考。在满足模型两侧量纲一致的前提下,综合考虑刀盘直径、掘进速度、刀盘转速等影响因素,因而预估准确可靠。并且能根据掘进地层的条件以及操作参数的变化调整扭矩参数,为掘进机装备动力系统的设计以及施工过程中扭矩参数的实时调整提供可靠的理论依据。
Description
技术领域
本发明涉及一种盾构刀盘扭矩的计算方法。特别是涉及一种基于量纲分析的预估盾构刀盘扭矩的计算方法。
背景技术
盾构是一种隧道掘进的专用工程机械,现已广泛用于地铁、铁路、公路、市政、水电等隧道工程。盾构的刀盘扭矩用于维持刀盘不断旋转切削前方土体,是贯穿装备运行始终的核心参量,也是装备能效控制与安全状态监控的基础参量,与装备尺寸和构型、地质力学参数、掘进控制参数等密切相关,影响因素十分复杂。预估刀盘扭矩是装备参数设计与施工调控的重要理论基础,也是盾构快速、安全、高效掘进的基本保证。
目前刀盘扭矩的已有理论建模方法主要围绕刀盘各个部件的几何特征、受载状态与土体本构,通过受力分析建立刀盘扭矩的预测模型。为满足理想的理论建模条件,预估过程中会存在假定与近似等问题。对于土体这类复杂材料,由于其物理力学性质深受土体成分、含水量、地质构造等因素影响,目前土体强度理论都只能适用于特定加载条件下某类土的物理力学性质,还没有准确的本构关系来描述其在实际工程中的力学性质。比如在计算刀盘摩擦扭矩这一重要分量时,刀盘与周边土体的摩擦系数需人为假定。这些假定条件都会导致刀盘扭矩的预估不准。在施工过程中刀盘扭矩的设定不当将直接造成刀具切削效率低、刀盘被卡事故等问题。
近年来,随着计算机技术的不断发展,基于数据的各类挖掘算法可为海量工程数据反演识别分析及建模提供有力的技术条件。然而,现阶段数据挖掘手段如人工神经网络、支持向量机等算法,其工程应用更多仅停留在数据识别层面,不能解释目标参数与影响因素间的因果关系。其“输入-输出”的黑箱结构难以建立表征问题的显式模型,直接导致这些数据分析方法得到的结果经验性过强,无法提取同类问题的共性规律,难以为盾构刀盘扭矩提供普适的建模方法。
因此,建立一种准确的、普适性高的刀盘扭矩预测方法,已成为本行业之急需。
发明内容
本发明所要解决的技术问题是,通过量纲分析以及统计工具,提供一种为实现隧道掘进装备参数控制依据的盾构刀盘扭矩的计算方法。
本发明所采用的技术方案是:一种基于量纲分析的盾构刀盘扭矩计算方法,包括如下步骤:
1)确定影响刀盘扭矩的参数;
2)根据所述的影响刀盘扭矩的参数分别计算主控项及无量纲化刀盘扭矩大小;
3)将计算得到的主控项及无量纲化刀盘扭矩进行汇总;
4)利用利用线性回归方法,得到主控项的系数识别;
5)根据主控项系数识别结果,得到刀盘扭矩的计算模型。
步骤1)所述的影响刀盘扭矩的主要参数包括:刀盘扭矩T,单位:kNm;土体剪切模量G,单位:kpa;土体容重γ,单位:kN/m3;隧道埋深H,单位:m;刀盘宽度W,单位:m;掘进速度v,单位:m/s;静止土压力系数K0;土舱压力pe,单位:kpa;盾构直径D,单位:m;刀盘转速ω,单位:r/min。
步骤2)所述的主控项有3个:J1、J2、J3,计算公式分别如下:
式中,K0为静止土压力系数;γ,单位:kN/m3,为土体容重;G,单位:kpa,为土体剪切模量;D,单位:m,为盾构直径;pe,单位:kpa,为土舱压力;v,单位:m/s,为掘进速度;ω,单位:r/min,为刀盘转速;W,单位:m,为刀盘厚度;H,单位:m为隧道埋深;T,单位:kNm,为刀盘扭矩。
矩阵J中第1列代表主控项J1、矩阵J中第2列代表主控项J2、矩阵J中第3列代表主控项J3;矩阵J与向量的第1~n行依次对应工程数据中第1~n个样本依据3个主控项J1、J2、J3和无量纲化刀盘扭矩的计算公式所计算的值。
步骤5)包括:
将系数识别结果β1、β2、β3与三个主控项J1、J2、J3分别相乘,得到刀盘扭矩T的计算模型如下:
其中,K0为静止土压力系数;γ,单位:kN/m3,为土体容重;G,单位:kpa,为土体剪切模量;D,单位:m,为盾构直径;pe,单位:kpa,为土舱压力;v,单位:m/s,为掘进速度;ω,单位:r/min,为刀盘转速;W,单位:m,为刀盘厚度;H,单位:m为隧道埋深。
本发明的一种基于量纲分析的盾构刀盘扭矩计算方法,其模型预测结果可为装备的智能化设计提供理论依据,为不同地质的安全施工提供参考。在满足模型两侧量纲一致的前提下,综合考虑刀盘直径、掘进速度、刀盘转速等影响因素,因而预估准确可靠。并且能根据掘进地层的条件以及操作参数的变化调整扭矩参数,为掘进机装备动力系统的设计以及施工过程中扭矩参数的实时调整提供可靠的理论依据。
具体实施方式
下面结合实施例对本发明的一种基于量纲分析的盾构刀盘扭矩计算方法做出详细说明。
本发明的一种基于量纲分析的盾构刀盘扭矩计算方法,包括如下步骤:
1)确定影响刀盘扭矩的参数;包括:
刀盘扭矩T,单位:kNm;土体剪切模量G,单位:kpa;土体容重γ,单位:kN/m3;隧道埋深H,单位:m;刀盘宽度W,单位:m;掘进速度v,单位:m/s;静止土压力系数K0;土舱压力pe,单位:kpa;盾构直径D,单位:m;刀盘转速ω,单位:r/min。
2)根据所述的影响刀盘扭矩的参数分别计算主控项及无量纲化刀盘扭矩大小;
所述的主控项有3个:J1、J2、J3,计算公式分别如下:
式中,K0为静止土压力系数;γ,单位:kN/m3,为土体容重;G,单位:kpa,为土体剪切模量;D,单位:m,为盾构直径;pe,单位:kpa,为土舱压力;v,单位:m/s,为掘进速度;ω,单位:r/min,为刀盘转速;W,单位:m,为刀盘厚度;H,单位:m为隧道埋深;T,单位:kNm,为刀盘扭矩。
3)将计算得到的主控项及无量纲化刀盘扭矩进行汇总;
4)利用线性回归方法,得到主控项的系数识别;包括:
矩阵J中第1列代表主控项J1、矩阵J中第2列代表主控项J2、矩阵J中第3列代表主控项J3;矩阵J与向量的第1~n行依次对应工程数据中第1~n个样本依据3个主控项J1、J2、J3和无量纲化刀盘扭矩的计算公式所计算的值。
5)根据主控项系数识别结果,得到刀盘扭矩的计算模型。包括:
将系数识别结果β1、β2、β3与三个主控项J1、J2、J3分别相乘,得到刀盘扭矩T的计算模型如下:
其中,K0为静止土压力系数;γ,单位:kN/m3,为土体容重;G,单位:kpa,为土体剪切模量;D,单位:m,为盾构直径;pe,单位:kpa,为土舱压力;v,单位:m/s,为掘进速度;ω,单位:r/min,为刀盘转速;W,单位:m,为刀盘厚度;H,单位:m,为隧道埋深。
下面给出具体实例:
1、列出影响刀盘扭矩T的主要参数及具体数值
T | G | γ | H | W | v | K<sub>0</sub> | p<sub>e</sub> | D | ω |
1082.35 | 12000 | 20.00 | 12.50 | 1.00 | 32.57 | 0.49 | 122.94 | 6.60 | 0.81 |
1011.76 | 23310 | 26.06 | 13.20 | 1.00 | 54.29 | 0.39 | 171.18 | 6.60 | 1.32 |
1200.00 | 16310 | 18.18 | 11.60 | 1.00 | 32.86 | 0.35 | 152.35 | 6.60 | 1.03 |
800.00 | 12740 | 14.59 | 10.20 | 1.00 | 40.29 | 0.43 | 110.06 | 6.60 | 1.19 |
941.18 | 19390 | 19.55 | 9.90 | 1.00 | 45.43 | 0.31 | 154.71 | 6.60 | 1.22 |
上表中,T(kNm)为刀盘扭矩;G(kpa)为土体剪切模量;γ(kN/m3)为土体容重;H(m)为隧道埋深;W(m)为刀盘宽度;v(m/s)为掘进速度;K0为静止土压力系数;pe(kpa)为土舱压力;D(m)为盾构直径;ω(r/min)为刀盘转速。
2、将上表中的各参数代入如下1-4式,分别计算主控项及无量纲化刀盘扭矩大小。其中,等式左边的J1、J2、J3为三个主控项,为无量纲化刀盘扭矩。等式右边,K0为静止土压力系数;γ(kN/m3)为土体容重;G(kpa)为土体剪切模量;D(m)为盾构直径;pe(kpa)为土舱压力;v(m/s)为掘进速度;ω(r/min)为刀盘转速;W(m)为刀盘厚度;H(m)为隧道埋深;T,单位:kNm,为刀盘扭矩。
3、将按1-4式计算得到的值汇总,如下表所示:
矩阵J中第1列代表主控项J1、矩阵J中第2列代表主控项J2、矩阵J中第3列代表主控项J3;矩阵J与向量的第1~n行依次对应工程数据中第1~n个样本依据3个主控项有J1、J2、J3和无量纲化刀盘扭矩的计算公式所计算的值。从而得到主控项的系数识别如下:
β1=9.09×10-5β2=5.47×10-4β3=9.72×10-2
5、将计算得到的系数识别结果β1、β2、β3的值与三个主控项J1、J2、J3分别相乘,得到刀盘扭矩T的计算模型,如式(7)所示:
将如下表所示的工程数据:
G | γ | H | W | v | K<sub>0</sub> | p<sub>e</sub> | D | ω |
22140 | 18.80 | 10.6 | 1.00 | 22.88 | 0.45 | 165.6 | 8.00 | 1.08 |
带入式(7)可得盾构刀盘扭矩T的大小如下:
T=1069.74(kNm)。
Claims (4)
1.一种基于量纲分析的盾构刀盘扭矩计算方法,其特征在于,包括如下步骤:
1)确定影响刀盘扭矩的参数;
2)根据所述的影响刀盘扭矩的参数分别计算主控项及无量纲化刀盘扭矩大小;
所述的主控项有3个:J1、J2、J3,计算公式分别如下:
式中,K0为静止土压力系数;γ,单位:kN/m3,为土体容重;G,单位:kpa,为土体剪切模量;D,单位:m,为盾构直径;pe,单位:kpa,为土舱压力;v,单位:m/s,为掘进速度;ω,单位:r/min,为刀盘转速;W,单位:m,为刀盘厚度;H,单位:m为隧道埋深;T,单位:kNm,为刀盘扭矩
3)将计算得到的主控项及无量纲化刀盘扭矩进行汇总;
4)利用线性回归方法,得到主控项的系数识别;
5)根据主控项系数识别结果,得到刀盘扭矩的计算模型。
2.根据权利要求1所述的一种基于量纲分析的盾构刀盘扭矩计算方法,其特征在于,步骤1)所述的影响刀盘扭矩的主要参数包括:刀盘扭矩T,单位:kNm;土体剪切模量G,单位:kpa;土体容重γ,单位:kN/m3;隧道埋深H,单位:m;刀盘宽度W,单位:m;掘进速度v,单位:m/s;静止土压力系数K0;土舱压力pe,单位:kpa;盾构直径D,单位:m;刀盘转速ω,单位:r/min。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810034645.5A CN108388689B (zh) | 2018-01-16 | 2018-01-16 | 一种基于量纲分析的盾构刀盘扭矩计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810034645.5A CN108388689B (zh) | 2018-01-16 | 2018-01-16 | 一种基于量纲分析的盾构刀盘扭矩计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108388689A CN108388689A (zh) | 2018-08-10 |
CN108388689B true CN108388689B (zh) | 2021-08-31 |
Family
ID=63076257
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810034645.5A Active CN108388689B (zh) | 2018-01-16 | 2018-01-16 | 一种基于量纲分析的盾构刀盘扭矩计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108388689B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109488321B (zh) * | 2019-01-03 | 2019-11-29 | 天津大学 | 一种盾构刀盘扭矩确定方法及系统 |
CN109766621B (zh) * | 2019-01-03 | 2020-08-21 | 天津大学 | 一种tbm刀盘扭矩确定方法及系统 |
CN110362899B (zh) * | 2019-07-01 | 2023-07-25 | 中铁隧道局集团有限公司 | 用于超大直径盾构过基岩突起段的刀盘扭矩计算方法 |
CN112065421B (zh) * | 2020-10-10 | 2022-04-22 | 中国铁建重工集团股份有限公司 | 一种掘进机刀盘自动定位方法 |
CN112329172B (zh) * | 2020-11-05 | 2022-12-27 | 上海交通大学 | 基于并联式神经网络的盾构机刀盘扭矩预测方法及系统 |
CN113435055B (zh) * | 2021-07-08 | 2022-11-22 | 上海交通大学 | 盾构刀盘扭矩领域自适应迁移预测方法和系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4151642A (en) * | 1976-10-04 | 1979-05-01 | General Dynamics Corporation | Device for pre-setting tools for numerical controlled machine |
CN105631150A (zh) * | 2016-01-05 | 2016-06-01 | 石家庄铁道大学 | 复合地层情况下盾构掘进参数的优化方法 |
CN105736008A (zh) * | 2016-02-01 | 2016-07-06 | 中铁建大桥工程局集团第三工程有限公司 | 富水砂、卵地层多平衡综合顶管机刀盘刀具的配置方法 |
-
2018
- 2018-01-16 CN CN201810034645.5A patent/CN108388689B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4151642A (en) * | 1976-10-04 | 1979-05-01 | General Dynamics Corporation | Device for pre-setting tools for numerical controlled machine |
CN105631150A (zh) * | 2016-01-05 | 2016-06-01 | 石家庄铁道大学 | 复合地层情况下盾构掘进参数的优化方法 |
CN105736008A (zh) * | 2016-02-01 | 2016-07-06 | 中铁建大桥工程局集团第三工程有限公司 | 富水砂、卵地层多平衡综合顶管机刀盘刀具的配置方法 |
Non-Patent Citations (3)
Title |
---|
Disc cutter wear prediction for a hard rock TBM cutterhead based on energy analysis;Zhang qian等;《Tunnelling and Underground Space Technology》;20150831;第50卷;第324-333页 * |
复杂地层盾构掘进速率和刀盘扭矩预测模型及其地层适应性研究;张志奇等;《隧道建设》;20161230;第36卷(第12期);第1449-1455页 * |
盾构施工振动振源的影响因素研究;王鑫等;《地震工程学报》;20140930;第36卷(第3期);第592-598页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108388689A (zh) | 2018-08-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108388689B (zh) | 一种基于量纲分析的盾构刀盘扭矩计算方法 | |
CN108446418B (zh) | 一种基于量纲分析的盾构掘进总推力计算方法 | |
Obregon et al. | Probabilistic approach for open pit bench slope stability analysis–A mine case study | |
CN109685378B (zh) | 一种基于数据挖掘的tbm施工围岩可掘性分级方法 | |
Vlachopoulos et al. | Appropriate uses and practical limitations of 2D numerical analysis of tunnels and tunnel support response | |
Zhao et al. | Effects of jointed rock mass and mixed ground conditions on the cutting efficiency and cutter wear of tunnel boring machine | |
Zhang et al. | Modeling of the thrust and torque acting on shield machines during tunneling | |
WO2019063024A1 (zh) | 硬岩tbm掘进控制参数智能决策方法及系统 | |
Liu et al. | Optimal earth pressure balance control for shield tunneling based on LS-SVM and PSO | |
CN109918822A (zh) | 基于随机森林回归的盾构掘进参数偏差计算方法及系统 | |
CN103870677A (zh) | 一种掘进机的掘进参数设定方法 | |
Zhang et al. | Mechanical characterization of the load distribution on the cutterhead–ground interface of shield tunneling machines | |
Yan et al. | Identification of geological characteristics from construction parameters during shield tunnelling | |
CN110630274B (zh) | 一种盾构隧道复合地层组段实施系统及实施方法 | |
CN109854259A (zh) | 一种盾构机施工掘进参数最佳取值范围的获取方法与系统 | |
CN103197609A (zh) | 数控加工动态特征建模方法 | |
CN111144001A (zh) | 一种基于bp神经网络的矿山井筒工程tbm控制方法 | |
CN106777568B (zh) | 全断面岩石掘进机刀盘掘进载荷的确定方法 | |
Qiu et al. | TBM tunnel surrounding rock classification method and real-time identification model based on tunneling performance | |
Hu et al. | Improved SVR method for predicting the cutting force of a TBM cutter using linear cutting machine test data | |
Song et al. | A cutter head torque forecast model based on multivariate nonlinear regression for EPB shield tunneling | |
Gu et al. | The application of intuitionistic fuzzy set-TOPSIS model on the level assessment of the surrounding rocks | |
CN104834986B (zh) | 基于全局综合敏感性分析的隧洞施工进度动态控制方法 | |
CN114036696B (zh) | 一种基于神经网络模型微调的刀盘扭矩预测方法及系统 | |
Liu et al. | Optimization Control of Energy Consumption in Tunneling System of Earth Pressure Balance Shield Tunneling Machine. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |