CN109766621B - TBM cutter head torque determination method and system - Google Patents

TBM cutter head torque determination method and system Download PDF

Info

Publication number
CN109766621B
CN109766621B CN201910004148.5A CN201910004148A CN109766621B CN 109766621 B CN109766621 B CN 109766621B CN 201910004148 A CN201910004148 A CN 201910004148A CN 109766621 B CN109766621 B CN 109766621B
Authority
CN
China
Prior art keywords
rock
torque
dimensionless
unit
coefficient
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
Application number
CN201910004148.5A
Other languages
Chinese (zh)
Other versions
CN109766621A (en
Inventor
张茜
周思阳
亢一澜
侯振德
蔡宗熙
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN201910004148.5A priority Critical patent/CN109766621B/en
Publication of CN109766621A publication Critical patent/CN109766621A/en
Application granted granted Critical
Publication of CN109766621B publication Critical patent/CN109766621B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Excavating Of Shafts Or Tunnels (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention discloses a TBM cutter torque determining method and system. The method comprises the following steps: determining parameters influencing the cutter head torque; determining a dimensionless quantity and a dimensionless cutterhead torque in the TBM tunneling process according to the parameters; obtaining a coefficient identification result of the dimensionless quantity by adopting an elastic network algorithm according to the dimensionless quantity; multiplying each dimensionless quantity coefficient identification result by the corresponding dimensionless quantity to obtain a cutter head torque calculation model; and determining the TBM cutter torque according to the cutter torque calculation model. By adopting the method or the system, the cutter torque in the TBM equipment tunneling process can be quickly and accurately calculated.

Description

TBM cutter head torque determination method and system
Technical Field
The invention relates to the field of TBM cutter torque calculation, in particular to a TBM cutter torque determination method and a TBM cutter torque determination system.
Background
A full-face Tunnel Boring Machine (TBM) is factory-like flow line tunnel construction equipment integrated by systems of machinery, electricity, liquid, light, gas and the like, and is mainly used for tunnel excavation of rock stratums. The torque of the cutter head is an important performance parameter and is used for maintaining the cutter head to continuously rotate to crush front rocks. At present, the research of predicting the TBM cutter head torque based on engineering data is mainly divided into a mechanical model based on stress analysis and prediction based on engineering data analysis. Among such studies based on engineering data analysis, there is a calculation method based on dimensional analysis, which has high universality. However, for predicting the cutter head torque of the heading machine in hard rock, the research based on dimension analysis does not eliminate irrelevant influence factors, so that the model is complicated and the requirements of hard rock tunnel engineering are not facilitated.
Disclosure of Invention
The invention aims to provide a TBM cutter torque determining method and a TBM cutter torque determining system, which can screen out main parameters influencing the TBM cutter torque and help TBM technical managers to adjust the cutter torque in time.
In order to achieve the purpose, the invention provides the following scheme:
a TBM cutterhead torque determination method, comprising:
determining parameters influencing the cutter head torque, wherein the parameters comprise rock saturated compressive strength, tunneling speed, cutter head rotating speed, rock firmness coefficient, rock tensile strength, rock friction coefficient, rock cohesion and rock longitudinal wave speed;
determining a dimensionless quantity and a dimensionless cutterhead torque in the TBM tunneling process according to the parameters;
obtaining a coefficient identification result of the dimensionless quantity by adopting an elastic network algorithm according to the dimensionless quantity;
multiplying each dimensionless quantity coefficient identification result by the corresponding dimensionless quantity to obtain a cutter head torque calculation model;
and determining the TBM cutter torque according to the cutter torque calculation model.
Optionally, the determining a dimensionless amount and a dimensionless cutterhead torque in the TBM tunneling process according to the parameters specifically includes:
according to the formula pi1=η、
Figure BDA0001934799480000029
π3=f、
Figure BDA0001934799480000022
Respectively determining dimensionless quantity pi in the process of TBM tunneling1,π2,π3,π4,π5
According to the formula
Figure BDA0001934799480000023
Determining non-dimensionalized cutterhead torque in TBM tunneling process
Figure BDA0001934799480000024
Wherein, pi1,π2,π3,π4,π5Five dimensionless quantities, η rock strength coefficient, sigmatTensile strength of rock, unit: kpa; sigmacIs the saturated compressive strength of rock, unit: kpa; f is the rock friction coefficient; c is rock cohesion in units: kpa; vpIs the rock longitudinal wave velocity, unit:m/s; v is the tunneling speed, unit: m/s; omega is the rotating speed of the cutter head, unit: r/min;
Figure BDA0001934799480000025
the cutter torque is dimensionless; t is TBM cutter torque, unit: kNm.
Optionally, obtaining the coefficient identification result of the dimensionless quantity by using an elastic network algorithm according to the dimensionless quantity specifically includes:
combining said dimensionless quantity pi1,π2,π3,π4,π5As independent variable, dimensionless cutter torque
Figure BDA0001934799480000026
Using an elastic network algorithm to adopt a formula as a dependent variable
Figure BDA0001934799480000027
Obtaining a dimensionless coefficient;
wherein n is the sample size; x is an independent variable, including1,π2,π3,π4,π5α is an identification coefficient including α12345(ii) a y is a dependent variable and y is a non-dimensionalized cutter torque
Figure BDA0001934799480000028
w is a penalty weight; r is the ratio occupied by the L1 norm regularization.
Optionally, the step of multiplying each non-dimensional coefficient identification result by the corresponding non-dimensional coefficient to obtain a cutter head torque calculation model specifically includes:
multiplying each dimensionless coefficient identification result by the corresponding dimensionless coefficient to obtain a cutter head torque calculation model T- α1ησcv3ω-32σtv3ω-33cv3ω-34cv3ω-35Vpσcv2ω-3
Wherein, α12345Coefficient of dimensionless quantity η rock strength coefficient sigmatTensile strength of rock, unit: kpa; sigmacIs the saturated compressive strength of rock, unit: kpa; f is the rock friction coefficient; c is rock cohesion in units: kpa; vpIs the rock longitudinal wave velocity, unit: m/s; v is the tunneling speed, unit: m/s; omega is the rotating speed of the cutter head, unit: r/min; t is TBM cutter torque, unit: kNm.
A TBM cutterhead torque determination system, comprising:
the parameter determination module is used for determining parameters influencing the cutter head torque, wherein the parameters comprise rock saturated compressive strength, tunneling speed, cutter head rotating speed, rock firmness coefficient, rock tensile strength, rock friction coefficient, rock cohesion and rock longitudinal wave speed;
the dimensionless parameter determining module is used for determining dimensionless quantity and dimensionless cutterhead torque in the TBM tunneling process according to the parameters;
the coefficient identification result determining module is used for obtaining the coefficient identification result of the dimensionless quantity by adopting an elastic network algorithm according to the dimensionless quantity;
the model establishing module is used for multiplying each dimensionless quantity coefficient identification result with the corresponding dimensionless quantity to obtain a cutter torque calculation model;
and the cutter head torque calculation module is used for determining the TBM cutter head torque according to the cutter head torque calculation model.
Optionally, the dimensionless parameter determining module specifically includes:
a dimensionless quantity parameter determining unit for determining the quantity of the measured quantity according to the formula pi1=η、
Figure BDA00019347994800000311
π3=f、
Figure BDA00019347994800000310
Figure BDA0001934799480000033
Respectively determining dimensionless quantity pi in the process of TBM tunneling1,π2,π3,π4,π5
A dimensionless cutter torque determination unit for determining the torque of the cutter according to a formula
Figure BDA0001934799480000034
Determining non-dimensionalized cutterhead torque in TBM tunneling process
Figure BDA0001934799480000035
Wherein, pi1,π2,π3,π4,π5Five dimensionless quantities, η rock strength coefficient, sigmatTensile strength of rock, unit: kpa; sigmacIs the saturated compressive strength of rock, unit: kpa; f is the rock friction coefficient; c is rock cohesion in units: kpa; vpIs the rock longitudinal wave velocity, unit: m/s; v is the tunneling speed, unit: m/s; omega is the rotating speed of the cutter head, unit: r/min;
Figure BDA0001934799480000036
the cutter torque is dimensionless; t is TBM cutter torque, unit: kNm.
Optionally, the coefficient identification result determining module specifically includes:
a coefficient recognition result determination unit for determining the dimensionless quantity pi1,π2,π3,π4,π5As independent variable, dimensionless cutter torque
Figure BDA0001934799480000037
Using an elastic network algorithm to adopt a formula as a dependent variable
Figure BDA0001934799480000038
Obtaining a dimensionless coefficient;
wherein n is the sample size;x is an independent variable, including1,π2,π3,π4,π5α is an identification coefficient including α12345(ii) a y is a dependent variable and y is a non-dimensionalized cutter torque
Figure BDA0001934799480000039
w is a penalty weight; r is the ratio occupied by the L1 norm regularization.
Optionally, the model building module specifically includes:
the model establishing unit is used for multiplying each non-dimensional coefficient identification result by the corresponding non-dimensional coefficient to obtain a cutter head torque calculation model T which is α1ησcv3ω-32σtv3ω-33cv3ω-34cv3ω-35Vpσcv2ω-3
Wherein, α12345Coefficient of dimensionless quantity η rock strength coefficient sigmatTensile strength of rock, unit: kpa; sigmacIs the saturated compressive strength of rock, unit: kpa; f is the rock friction coefficient; c is rock cohesion in units: kpa; vpIs the rock longitudinal wave velocity, unit: m/s; v is the tunneling speed, unit: m/s; omega is the rotating speed of the cutter head, unit: r/min; t is TBM cutter torque, unit: kNm.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects: the invention provides a TBM cutter head torque determination method, which comprises the following steps: determining parameters influencing the cutter head torque, wherein the parameters comprise rock saturated compressive strength, tunneling speed, cutter head rotating speed, rock firmness coefficient, rock tensile strength, rock friction coefficient, rock cohesion and rock longitudinal wave speed; determining a dimensionless quantity and a dimensionless cutterhead torque in the TBM tunneling process according to the parameters; obtaining a coefficient identification result of the dimensionless quantity by adopting an elastic network algorithm according to the dimensionless quantity; multiplying each dimensionless quantity coefficient identification result by the corresponding dimensionless quantity to obtain a cutter head torque calculation model; and determining the TBM cutter torque according to the cutter torque calculation model. By adopting the method or the system, the cutter torque in the TBM equipment tunneling process can be quickly and accurately calculated. On the basis of a basic mechanical rule, the elastic network algorithm is fused, the algorithm is a weighted algorithm between a lasso algorithm and a ridge regression algorithm, and a model obtained based on the algorithm can inherit the sparsity of the lasso algorithm and the stability of the ridge regression, so that effective decision support is provided for TBM technical managers to adjust the cutter torque in time.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
FIG. 1 is a flow chart of a TBM cutterhead torque determination method in accordance with an embodiment of the present invention;
figure 2 is a block diagram of a TBM cutterhead torque determination system in accordance with an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The invention aims to provide a TBM cutter torque determining method and a TBM cutter torque determining system, which can screen out main parameters influencing the TBM cutter torque and help TBM technical managers to adjust the cutter torque in time.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
FIG. 1 is a flow chart of a TBM cutterhead torque determination method according to an embodiment of the invention. As shown in fig. 1, a TBM cutterhead torque determination method includes:
step 101: determining parameters influencing the cutter head torque, wherein the parameters comprise rock saturated compressive strength, tunneling speed, cutter head rotating speed, rock firmness coefficient, rock tensile strength, rock friction coefficient, rock cohesion and rock longitudinal wave speed;
step 102: determining a dimensionless quantity and a dimensionless cutterhead torque in the TBM tunneling process according to the parameters;
step 103: obtaining a coefficient identification result of the dimensionless quantity by adopting an elastic network algorithm according to the dimensionless quantity;
step 104: multiplying each dimensionless quantity coefficient identification result by the corresponding dimensionless quantity to obtain a cutter head torque calculation model;
step 105: and determining the TBM cutter torque according to the cutter torque calculation model.
Step 102, specifically comprising:
according to the formula pi1=η、
Figure BDA0001934799480000055
π3=f、
Figure BDA0001934799480000052
Respectively determining dimensionless quantity pi in the process of TBM tunneling1,π2,π3,π4,π5
According to the formula
Figure BDA0001934799480000053
Determining non-dimensionalized cutterhead torque in TBM tunneling process
Figure BDA0001934799480000054
Wherein, pi1,π2,π3,π4,π5Five dimensionless quantities, η rock strength coefficient, sigmatTensile strength of rock, unit: kpa; sigmacIs the saturated compressive strength of rock, unit: kpa; f is the rock friction coefficient; c is rock cohesion in units: kpa; vpIs the rock longitudinal wave velocity, unit: m/s; v is the tunneling speed, unit: m/s; omega is the rotating speed of the cutter head, unit: r/min;
Figure BDA0001934799480000061
the cutter torque is dimensionless; t is TBM cutter torque, unit: kNm.
Step 103, specifically comprising:
combining said dimensionless quantity pi1,π2,π3,π4,π5As independent variable, dimensionless cutter torque
Figure BDA0001934799480000062
Using an elastic network algorithm to adopt a formula as a dependent variable
Figure BDA0001934799480000063
Obtaining a dimensionless coefficient;
wherein n is the sample size; x is an independent variable, including1,π2,π3,π4,π5α is an identification coefficient including α12345(ii) a y is a dependent variable and y is a non-dimensionalized cutter torque
Figure BDA0001934799480000065
w is a penalty weight; r is the ratio occupied by the L1 norm regularization.
L1 norm regularization adds L1 norm to the cost function, so that the learning result meets the sparsification requirement, and people can conveniently extract features.
Step 104, specifically comprising:
each will beMultiplying the non-dimensional quantity coefficient identification result by the corresponding non-dimensional quantity to obtain a cutter head torque calculation model T- α1ησcv3ω-32σtv3ω-33cv3ω-34cv3ω-35Vpσcv2ω-3
Wherein T is cutterhead torque with unit of kNm α12345Coefficient of dimensionless quantity η rock strength coefficient sigmatTensile strength of rock, unit: kpa; sigmacIs the saturated compressive strength of rock, unit: kpa; f is the rock friction coefficient; c is rock cohesion in units: kpa; vpIs the rock longitudinal wave velocity, unit: m/s; v is the tunneling speed, unit: m/s; omega is the rotating speed of the cutter head, unit: r/min;
Figure BDA0001934799480000064
the cutter torque is dimensionless; t is TBM cutter torque, unit: kNm.
On the basis of a basic mechanical rule, the elastic network algorithm is fused, the algorithm is a weighted algorithm between a lasso algorithm and a ridge regression algorithm, and a model obtained based on the algorithm can inherit the sparsity of the lasso algorithm and the stability of the ridge regression, so that effective decision support is provided for TBM technical managers to adjust the cutter torque in time.
Figure 2 is a block diagram of a TBM cutterhead torque determination system in accordance with an embodiment of the present invention. As shown in fig. 2, a TBM cutterhead torque determination system includes:
the parameter determining module 201 is used for determining parameters influencing the cutter head torque, wherein the parameters comprise rock saturated compressive strength, tunneling speed, cutter head rotating speed, rock firmness coefficient, rock tensile strength, rock friction coefficient, rock cohesion and rock longitudinal wave speed;
the dimensionless parameter determining module 202 is used for determining dimensionless quantity and dimensionless cutterhead torque in the TBM tunneling process according to the parameters;
a coefficient identification result determining module 203, configured to obtain a coefficient identification result of the dimensionless quantity by using an elastic network algorithm according to the dimensionless quantity;
the model establishing module 204 is used for multiplying each non-dimensional coefficient identification result by the corresponding non-dimensional coefficient to obtain a cutter head torque calculation model;
and the cutterhead torque calculation module 205 is used for determining the TBM cutterhead torque according to the cutterhead torque calculation model.
The dimensionless parameter determining module 202 specifically includes:
a dimensionless quantity parameter determining unit for determining the quantity of the measured quantity according to the formula pi1=η、
Figure BDA00019347994800000711
π3=f、
Figure BDA00019347994800000710
Figure BDA0001934799480000073
Respectively determining dimensionless quantity pi in the process of TBM tunneling1,π2,π3,π4,π5
A dimensionless cutter torque determination unit for determining the torque of the cutter according to a formula
Figure BDA0001934799480000074
Determining non-dimensionalized cutterhead torque in TBM tunneling process
Figure BDA0001934799480000075
Wherein, pi1,π2,π3,π4,π5Five dimensionless quantities, η rock strength coefficient, sigmatTensile strength of rock, unit: kpa; sigmacIs the saturated compressive strength of rock, unit: kpa; f is the rock friction coefficient; c is rock cohesion in units: kpa; vpIs the rock longitudinal wave velocity, unit: m/s; v is the tunneling speed, unit: m/s; omega is the rotating speed of the cutter head, unit: r/min;
Figure BDA0001934799480000076
the cutter torque is dimensionless; t is TBM cutter torque, unit: kNm.
The coefficient identification result determining module 203 specifically includes:
a coefficient recognition result determination unit for determining the dimensionless quantity pi1,π2,π3,π4,π5As independent variable, dimensionless cutter torque
Figure BDA0001934799480000077
Using an elastic network algorithm to adopt a formula as a dependent variable
Figure BDA0001934799480000078
Obtaining a dimensionless coefficient;
wherein n is the sample size; x is an independent variable, including1,π2,π3,π4,π5α is an identification coefficient including α12345(ii) a y is a dependent variable and y is a non-dimensionalized cutter torque
Figure BDA0001934799480000079
w is a penalty weight; r is the ratio occupied by the L1 norm regularization.
The model building module 204 specifically includes:
the model establishing unit is used for multiplying each non-dimensional coefficient identification result by the corresponding non-dimensional coefficient to obtain a cutter head torque calculation model T which is α1ησcv3ω-32σtv3ω-33cv3ω-34cv3ω-35Vpσcv2ω-3
Wherein, α12345Coefficient of dimensionless quantity η rock strength coefficient sigmatTensile strength of rock, unit: kpa; sigmacIs the saturated compressive strength of rock, unit: kpa; f is the rock friction coefficient; c is rock cohesion in units: kpa; vpIs the rock longitudinal wave velocity, unit: m/s; v is the tunneling speed, unit: m/s; omega is the rotating speed of the cutter head, unit: r/min; t is TBM cutter torque, unit: kNm.
Specific example 1:
the method for calculating the cutter head torque according to the present invention will be described in detail with reference to the following embodiments. The method comprises the following specific steps:
1. the parameters influencing the cutter head torque T are listed, η is a rock strength coefficient, and sigmatTensile strength of rock, unit: kpa; sigmacIs the saturated compressive strength of rock, unit: kpa; f is the rock friction coefficient; c is rock cohesion in units: kpa; vpIs the rock longitudinal wave velocity, unit: m/s; v is the tunneling speed, unit: m/s; omega is the rotating speed of the cutter head, unit: r/min.
2. The dimensionless and dimensionless cutterhead torque magnitudes are calculated according to the following equations (1) - (6), respectively. Wherein, pi1,π2,π3,π4,π5Five dimensionless quantities, η rock strength coefficient, sigmatTensile strength of rock, unit: kpa; sigmacIs the saturated compressive strength of rock, unit: kpa; f is the rock friction coefficient; c is rock cohesion in units: kpa; vpIs the rock longitudinal wave velocity, unit: m/s; v is the tunneling speed, unit: m/s; omega is the rotating speed of the cutter head, unit: r/min;
Figure BDA0001934799480000081
the cutter torque is dimensionless; t is TBM cutter torque, unit: kNm.
π1=η (1)
Figure BDA0001934799480000082
π3=f (3)
Figure BDA0001934799480000083
Figure BDA0001934799480000084
Figure BDA0001934799480000085
3. Summarizing the values calculated according to the formulas 1 to 6, and dividing pi1,π2,π3,π4,π5Using the five dimensionless quantities as independent variables and the dimensionless cutter torque as dependent variables, and using elastic network algorithm to obtain the dimensionless coefficient α1~α5The identification calculation formula (2) is as follows:
α1=0,α2=0,α3=0,α4=0,α5=2395.68.
4. identifying α the coefficient of each dimensionless quantity1~α5With five dimensionless quantities pi1~π5Multiplying respectively to obtain a calculation model of the TBM cutter head torque T: as shown in the following formula:
T=2395.68Vpσcv2w-3(7)
engineering data as shown in the following table:
σc v ω η σt f c Vp
20530 0.000167 1.2 17 1644 0.26 170 3120
the magnitude of the TBM cutterhead torque T available in the belt-in (7) is as follows:
T=2476.63(kNm)。
the embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other. For the system disclosed by the embodiment, the description is relatively simple because the system corresponds to the method disclosed by the embodiment, and the relevant points can be referred to the method part for description.
The principles and embodiments of the present invention have been described herein using specific examples, which are provided only to help understand the method and the core concept of the present invention; meanwhile, for a person skilled in the art, according to the idea of the present invention, the specific embodiments and the application range may be changed. In view of the above, the present disclosure should not be construed as limiting the invention.

Claims (6)

1. A TBM cutterhead torque determination method, comprising:
determining parameters influencing the cutter head torque, wherein the parameters comprise rock saturated compressive strength, tunneling speed, cutter head rotating speed, rock firmness coefficient, rock tensile strength, rock friction coefficient, rock cohesion and rock longitudinal wave speed;
determining a dimensionless quantity and a dimensionless cutterhead torque in the TBM tunneling process according to the parameters;
obtaining a coefficient identification result of the dimensionless quantity by adopting an elastic network algorithm according to the dimensionless quantity;
multiplying each dimensionless quantity coefficient identification result by the corresponding dimensionless quantity to obtain a cutter head torque calculation model;
determining TBM cutter torque according to the cutter torque calculation model;
the method for determining the dimensionless quantity and the dimensionless cutterhead torque in the TBM tunneling process according to the parameters specifically comprises the following steps:
according to the formula pi1=η、
Figure FDA0002542989160000011
π3=f、
Figure FDA0002542989160000012
Respectively determining dimensionless quantity pi in the process of TBM tunneling1,π2,π3,π4,π5
According to the formula
Figure FDA0002542989160000013
Determining non-dimensionalized cutterhead torque in TBM tunneling process
Figure FDA0002542989160000014
Wherein, pi1,π2,π3,π4,π5Five dimensionless quantities η being rockStone firmness factor; sigmatTensile strength of rock, unit: kpa; sigmacIs the saturated compressive strength of rock, unit: kpa; f is the rock friction coefficient; c is rock cohesion in units: kpa; vpIs the rock longitudinal wave velocity, unit: m/s; v is the tunneling speed, unit: m/s; omega is the rotating speed of the cutter head, unit: r/min;
Figure FDA0002542989160000017
the cutter torque is dimensionless; t is cutter head torque, unit: kNm.
2. The method for determining the torque of the cutterhead of the TBM according to claim 1, wherein the obtaining of the coefficient identification result of the dimensionless quantity by using an elastic network algorithm according to the dimensionless quantity specifically comprises:
combining said dimensionless quantity pi1,π2,π3,π4,π5As independent variable, dimensionless cutter torque
Figure FDA0002542989160000015
Using an elastic network algorithm to adopt a formula as a dependent variable
Figure FDA0002542989160000016
Obtaining a dimensionless coefficient;
wherein n is the sample size; x is an independent variable, including1,π2,π3,π4,π5α is an identification coefficient including α12345(ii) a y is a dependent variable and y is a non-dimensionalized cutter torque
Figure FDA0002542989160000023
w is a penalty weight; r is the ratio occupied by the L1 norm regularization.
3. The method for determining the torque of the cutterhead of the TBM according to claim 2, wherein the step of multiplying each non-dimensional coefficient identification result by the corresponding non-dimensional coefficient to obtain a cutterhead torque calculation model specifically comprises:
multiplying each dimensionless coefficient identification result by the corresponding dimensionless coefficient to obtain a cutter head torque calculation model T- α1ησcv3ω-32σtv3ω-33cv3ω-34cv3ω-35Vpσcv2ω-3
Wherein, α12345Coefficient of dimensionless quantity η rock strength coefficient sigmatTensile strength of rock, unit: kpa; sigmacIs the saturated compressive strength of rock, unit: kpa; f is the rock friction coefficient; c is rock cohesion in units: kpa; vpIs the rock longitudinal wave velocity, unit: m/s; v is the tunneling speed, unit: m/s; omega is the rotating speed of the cutter head, unit: r/min; t is cutter head torque, unit: kNm.
4. A TBM cutterhead torque determination system, comprising:
the parameter determination module is used for determining parameters influencing the cutter head torque, wherein the parameters comprise rock saturated compressive strength, tunneling speed, cutter head rotating speed, rock firmness coefficient, rock tensile strength, rock friction coefficient, rock cohesion and rock longitudinal wave speed;
the dimensionless parameter determining module is used for determining dimensionless quantity and dimensionless cutterhead torque in the TBM tunneling process according to the parameters;
the coefficient identification result determining module is used for obtaining the coefficient identification result of the dimensionless quantity by adopting an elastic network algorithm according to the dimensionless quantity;
the model establishing module is used for multiplying each dimensionless quantity coefficient identification result with the corresponding dimensionless quantity to obtain a cutter torque calculation model;
the cutter head torque calculation module is used for determining TBM cutter head torque according to the cutter head torque calculation model;
the dimensionless parameter determining module specifically includes:
a dimensionless quantity parameter determining unit for determining the quantity of the measured quantity according to the formula pi1=η、
Figure FDA0002542989160000021
π3=f、
Figure FDA0002542989160000022
Respectively determining dimensionless quantity pi in the process of TBM tunneling1,π2,π3,π4,π5
A dimensionless cutter torque determination unit for determining the torque of the cutter according to a formula
Figure FDA0002542989160000031
Determining non-dimensionalized cutterhead torque in TBM tunneling process
Figure FDA0002542989160000032
Wherein, pi1,π2,π3,π4,π5Five dimensionless quantities, η rock strength coefficient, sigmatTensile strength of rock, unit: kpa; sigmacIs the saturated compressive strength of rock, unit: kpa; f is the rock friction coefficient; c is rock cohesion in units: kpa; vpIs the rock longitudinal wave velocity, unit: m/s; v is the tunneling speed, unit: m/s; omega is the rotating speed of the cutter head, unit: r/min;
Figure FDA0002542989160000035
the cutter torque is dimensionless; t is cutter head torque, unit: kNm.
5. The TBM cutterhead torque determination system of claim 4, wherein the coefficient identification result determination module specifically includes:
a coefficient recognition result determination unit for determining the dimensionless quantity pi1,π2,π3,π4,π5As independent variable, dimensionless cutter torque
Figure FDA0002542989160000033
Using an elastic network algorithm to adopt a formula as a dependent variable
Figure FDA0002542989160000034
Obtaining a dimensionless coefficient;
wherein n is the sample size; x is an independent variable, including1,π2,π3,π4,π5α is an identification coefficient including α12345(ii) a y is a dependent variable and y is a non-dimensionalized cutter torque
Figure FDA0002542989160000036
w is a penalty weight; r is the ratio occupied by the L1 norm regularization.
6. The TBM cutterhead torque determination system of claim 5, wherein the model building module specifically includes:
the model establishing unit is used for multiplying each non-dimensional coefficient identification result by the corresponding non-dimensional coefficient to obtain a cutter head torque calculation model T which is α1ησcv3ω-32σtv3ω-33cv3ω-34cv3ω-35Vpσcv2ω-3
Wherein, α12345Coefficient of dimensionless quantity η rock strength coefficient sigmatTensile strength of rock, unit: kpa; sigmacIs the saturated compressive strength of rock, unit: kpa; f is the rock friction coefficient; c is in rockForce concentration, unit: kpa; vpIs the rock longitudinal wave velocity, unit: m/s; v is the tunneling speed, unit: m/s; omega is the rotating speed of the cutter head, unit: r/min; t is cutter head torque, unit: kNm.
CN201910004148.5A 2019-01-03 2019-01-03 TBM cutter head torque determination method and system Active CN109766621B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910004148.5A CN109766621B (en) 2019-01-03 2019-01-03 TBM cutter head torque determination method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910004148.5A CN109766621B (en) 2019-01-03 2019-01-03 TBM cutter head torque determination method and system

Publications (2)

Publication Number Publication Date
CN109766621A CN109766621A (en) 2019-05-17
CN109766621B true CN109766621B (en) 2020-08-21

Family

ID=66452619

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910004148.5A Active CN109766621B (en) 2019-01-03 2019-01-03 TBM cutter head torque determination method and system

Country Status (1)

Country Link
CN (1) CN109766621B (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102147826A (en) * 2011-03-15 2011-08-10 天津大学 Method for calculating optimal driving speed of tunneling machine under different geologies
CN204436401U (en) * 2014-12-11 2015-07-01 徐工集团凯宫重工南京有限公司 A kind of shield machine cutter head torque input structure
CN106351667A (en) * 2016-08-26 2017-01-25 中国铁建重工集团有限公司 Cutter head torque abnormality monitoring system for TBM

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5019452B2 (en) * 2007-09-05 2012-09-05 いわみ中央農業協同組合 Ripe grape harvest method and packaged grapes
CN102129499B (en) * 2011-03-15 2012-07-11 天津大学 Calculation method of drive loads of cutterhead of tunnel borer
CN108388689B (en) * 2018-01-16 2021-08-31 天津大学 Shield cutter head torque calculation method based on dimensional analysis

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102147826A (en) * 2011-03-15 2011-08-10 天津大学 Method for calculating optimal driving speed of tunneling machine under different geologies
CN204436401U (en) * 2014-12-11 2015-07-01 徐工集团凯宫重工南京有限公司 A kind of shield machine cutter head torque input structure
CN106351667A (en) * 2016-08-26 2017-01-25 中国铁建重工集团有限公司 Cutter head torque abnormality monitoring system for TBM

Also Published As

Publication number Publication date
CN109766621A (en) 2019-05-17

Similar Documents

Publication Publication Date Title
CN106777528B (en) The holographic forecast method of mine air-required volume
CN110685651B (en) Yield splitting method and system for multilayer commingled production gas well
CN108388689B (en) Shield cutter head torque calculation method based on dimensional analysis
CN105139388A (en) Method and apparatus for building facade damage detection in oblique aerial image
CN112950771B (en) Detection method and system for surrounding rock structure and related components
CN112946783B (en) Hydrate saturation determination method, device and equipment
CN111140244B (en) Intelligent support grade recommendation method for hard rock heading machine
CN112610277B (en) Method for rapidly predicting geological disaster risk of water-rich weak surrounding rock tunnel
WO2019091134A1 (en) Method for predicting water-rich levels of sandstone in coal seam roof
CN108280289B (en) Rock burst danger level prediction method based on local weighted C4.5 algorithm
CN112632852A (en) Karst area subway tunnel shield tunneling speed prediction method and device
CN109711079B (en) TBM tunneling total thrust determining method and system
CN113361043B (en) Method and system for predicting specific energy of cutter head of hard rock tunnel boring machine
CN109766621B (en) TBM cutter head torque determination method and system
CN111986041B (en) Oil reservoir development method and device
CN109736819B (en) Method and system for determining total shield tunneling thrust
CN112115537B (en) TBM (Tunnel boring machine) digability grading method
CN116663903A (en) Drilling and blasting method construction risk evaluation method based on combined weighting-nonlinear FAHP
CN116402339A (en) Method, system, equipment and medium for evaluating shield tunnel construction risk level
CN114036696B (en) Cutterhead torque prediction method and system based on neural network model fine adjustment
CN113901939B (en) Rock burst danger level prediction method based on fuzzy correction, storage medium and equipment
CN113836617B (en) Method and system for analogically predicting front fissure of tunneling surface based on time sequence model
CN109488321B (en) A kind of Cutter Head Torque in Shield Tunneling determines method and system
CN111222227B (en) Mine curve determining method and system
CN113933297A (en) Tunnel surrounding rock grading method and device, electronic equipment and medium

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