CN110704963B - Rapid method for optimizing tunneling parameters of earth pressure balanced type shield tunneling machine - Google Patents

Rapid method for optimizing tunneling parameters of earth pressure balanced type shield tunneling machine Download PDF

Info

Publication number
CN110704963B
CN110704963B CN201910871045.9A CN201910871045A CN110704963B CN 110704963 B CN110704963 B CN 110704963B CN 201910871045 A CN201910871045 A CN 201910871045A CN 110704963 B CN110704963 B CN 110704963B
Authority
CN
China
Prior art keywords
tunneling
formula
stratum
hob
efficiency
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
CN201910871045.9A
Other languages
Chinese (zh)
Other versions
CN110704963A (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.)
Nanjing Kentop Civil Engineering Technology Co ltd
Nanjing Tech University
Original Assignee
Nanjing Kentop Civil Engineering Technology Co ltd
Nanjing Tech 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 Nanjing Kentop Civil Engineering Technology Co ltd, Nanjing Tech University filed Critical Nanjing Kentop Civil Engineering Technology Co ltd
Priority to CN201910871045.9A priority Critical patent/CN110704963B/en
Publication of CN110704963A publication Critical patent/CN110704963A/en
Application granted granted Critical
Publication of CN110704963B publication Critical patent/CN110704963B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Excavating Of Shafts Or Tunnels (AREA)

Abstract

The invention relates to a rapid method for optimizing tunneling parameters of an earth pressure balance type shield tunneling machine, and belongs to the technical field of underground tunnel engineering construction. For quantitatively optimizing the excavation parameters of the earth pressure balanced type shield machine, the excavation efficiency is taken as an optimized analysis target function, and a prediction equation of the excavation efficiency under a limited sample is obtained by performing symbolic regression and linear regression on actually measured data. And obtaining the tunneling efficiency prediction equation suitable for the earth pressure balance type shield tunneling hard rock by calculating a continuous function of the coefficient of the tunneling efficiency prediction equation under the limited sample along with the change of the saturated uniaxial compressive strength of the tunneling stratum. And then, optimal calculation is carried out in a given definition domain through a step-by-step point taking method, and optimal tunneling parameters of the earth pressure balanced type shield tunneling machine suitable for tunneling hard rock are obtained by combining checking of the tunneling rate. The method is based on engineering investigation data and construction data, has simple and clear calculation process, reasonable method and strong practicability, and can effectively improve the tunneling efficiency of the earth pressure balance type shield tunneling machine.

Description

Rapid method for optimizing tunneling parameters of earth pressure balanced type shield tunneling machine
Technical Field
The invention belongs to the technical field of excavation parameter optimization of a full-face tunnel boring machine, and particularly relates to a rapid method for optimizing excavation parameters of an earth pressure balance type shield machine.
Background
1) A method for predicting the tunneling efficiency of a hard rock tunneling machine (application number: 201610755946.8) the method does not give specific definition of the tunneling efficiency, and the influence factors of the tunneling efficiency are limited to five types of geological features, and the influence of other tunneling parameters on the tunneling efficiency is not considered.
2) TBM tunneling parameter optimization method (application number: 201910231527.8) the method uses the minimum specific energy consumption as the objective function, the integral variables of the work function are displacement and rotation angle, the equation of force and torque can only get the broken line equation through the existing measured value, and the non-tunneling area can not be predicted and optimized in advance.
3) Optimization method of shield tunneling parameters under composite stratum condition (application number: 201610003385.6) and 4) establishment and optimization of a shield tunneling rate model under a composite stratum (lige, fuchsia, guojinbo, et al. establishment and optimization of a shield tunneling rate model under a composite stratum [ J ] modern tunneling technology, 2017 (3)) have limitations in applicability and calculation principle. In the aspect of applicability, the tunneling rate prediction equation is only applicable to one stratum, and the equation is not equally applicable when stratum conditions change. In the aspect of calculation principle, in reference 3), the tunneling rate and the cutter head torque are optimized respectively, but coupling optimization of the tunneling rate and the cutter head torque is not performed, that is, the tunneling parameter group corresponding to the optimal tunneling rate value is different from the tunneling parameter group corresponding to the optimal cutter head torque value, and a final optimization conclusion of comprehensively considering the tunneling rate and the cutter head torque is not given; reference 4) is purely based on the tunneling rate as an optimization target without considering the constraints of equipment loss and load capacity. Further, the constraints of the references 3) and 4) do not give a clear calculation formula, and it is difficult to solve the problem by the K-T method.
5) The method comprises the steps of establishing a shield tunneling rate model and optimizing parameters of a shield tunneling machine under a complex water-rich stratum (Wangqiang. establishing a shield tunneling rate model and optimizing parameters of a shield tunneling machine under a complex water-rich stratum [ J ] water conservancy planning and design, 2019(08):73-78.) by taking a tunneling rate as a single optimized objective function, the restriction of equipment load capacity and loss on tunneling efficiency is not considered, a tunneling rate prediction equation is only suitable for one stratum, and the applicability and principle correctness are limited. In addition, in contrast document 5), optimization calculation is performed by a local optimization method, and the calculation result is greatly affected by the setting of an initial value, resulting in low engineering applicability.
Disclosure of Invention
In order to overcome the defects of the prior art, the invention provides a rapid method for optimizing the tunneling parameters of an earth pressure balanced type shield tunneling machine, which is used for quantitatively optimizing the tunneling parameters of the earth pressure balanced type shield tunneling machine and taking the ratio of the tunneling speed to the cutter head torque as an optimized analysis objective function and specifically comprises the following steps:
step one, establishing a tunneling efficiency statistical sample
The tunneling efficiency β is defined as the tunneling rate v divided by the cutterhead torque T, as formula (1). T directly reflects the operating state of the equipment, and β reflects the relationship between the tunneling rate and the loading capacity of the cutterhead drive motor.
Figure BDA0002202823160000021
And combining the tunneling parameters automatically recorded by the shield tunneling machine during the trial tunneling of the tunnel and the tunneling parameter data of the previous engineering, taking the tunneling speed v, the cutter torque T, the effective thrust F, the soil bin pressure p and the cutter rotating speed n which are actually measured at the same moment as a group of tunneling parameters, and taking β, F, p and n at each recording moment in the trial tunneling stage as a tunneling efficiency statistical sample.
Step two, obtaining a tunneling efficiency prediction equation for the limited sample through symbolic regression and linear regression
Uniaxial compressive strength R to formation saturationcPerforming symbolic regression on the tunneling efficiency β of the shield tunneling machine in different strata respectively, and running the symbolic regression algorithm on the computer with the computation time of less than 1 minute and the fitting precision r2Fitting equations for β above 0.7 are candidate equations.
Counting candidate equations of each stratum, independent variable n2P, np, pF all co-exist in at least one of the candidate equations for each formation, thus in n2P, np and pF are used as independent variables of the tunneling efficiency prediction equation under the limited sample, as shown in formula (2), linear regression is carried out again on the tunneling efficiency statistical sample through the tunneling efficiency prediction equation under the limited sample, and equation coefficients belonging to each stratum are obtained;
β is the predicted value of tunneling efficiency under limited samples.
β*＝a1n2+a2p+a3np+a4pF+a5(2)
Step three tunneling efficiency prediction equation continuation
The equation coefficients in the formula (2) obtained under the condition of the limited samples are discrete equation coefficients for the limited kinds of strata, and the equations corresponding to the discrete equation coefficients need to be extended to obtain the equations suitable for the general strata.
Taking coefficients of prediction equations of the tunneling efficiency in different tunneling strata under the limited samples obtained in the step one as dependent variables, and taking the saturated uniaxial compressive strength R of the tunneling strata ascLinear and nonlinear regression for independent variables to obtain equation coefficients with RcThe continuous function of the change is shown as formula (3), formula (4), formula (5), formula (6) and formula (7). Because the samples are all obtained from hard rocks, the continuation of the tunneling efficiency prediction equation is also in the hard rock range, namely Rc≥30MPa。
a1＝116.588-99.614Rc 0.038,Rc≥30MPa (3)
a2＝46.522-12.532Rc 0.313,Rc≥30MPa (4)
a3＝0.093Rc-4.414,Rc≥30MPa (5)
a4＝1.135×10-6Rc-6.976×10-5,Rc≥30MPa (6)
Figure BDA0002202823160000031
According to the formula (8), the formula (3), the formula (4), the formula (5), the formula (6) and the formula (7) are substituted into the formula (2), so that a tunneling efficiency prediction equation suitable for the soil pressure balance type shield tunneling hard rock is obtained, and β is a tunneling efficiency prediction value of the soil pressure balance type shield tunneling hard rock.
Figure BDA0002202823160000032
Step four, determining an optimized calculation formula and boundary conditions
The optimal calculation formula of the tunneling parameter is shown as a formula (9), pΩ、ΘΩ、ΨΩRespectively reflects the allowable fluctuation range of the soil bin pressure, the effective thrust and the cutter head rotating speed when the stratum is normally tunneled omega, namely [ pmin,pmax]、[Fmin,Fmax]、[nmin,nmax]。
Figure BDA0002202823160000041
Step five-step dot method to calculate β larger value
Influence factors of the optimization problem are quantized, the boundary is clear, linear equal division is adopted to obtain the optimal value grid of the tunneling parameters to be selected, and step-by-step point-taking method gradual operation is carried out on grid nodes.
P is to beΩ、ΘΩ、ΨΩRespectively as emptyAnd defining domains on the X, Y and Z coordinate axes which are mutually orthogonal in the rectangular coordinate system are obtained, and β -star space defining domains are obtained.
Establishing a first step grid, and respectively connecting pΩ、ΘΩ、ΨΩAs m1Dividing equally and making normal planes of coordinate axes of the equally divided points at all equally divided points, wherein the normal planes are intersected to form β -x space definition domain first step grid, and any one first step grid node
Figure BDA0002202823160000042
Satisfies the formulas (10), (11) and (12). Determining tunneling stratum R according to stratum investigation datacThen, the points are put
Figure BDA0002202823160000043
And (3) substituting the formula (8) to obtain β corresponding to each first step grid node in the spatial definition domain of β, arranging β corresponding to each first step grid node in the spatial definition domain of β in descending order of magnitude to obtain a first step grid node β number sequence, and taking the grid nodes corresponding to β which is the first 20% of the maximum in the first step grid nodes β as the first step grid dominant points.
Figure BDA0002202823160000044
Figure BDA0002202823160000045
Figure BDA0002202823160000046
Establishing a second step grid, and taking the dominant point of each first step grid as a body center and the projection lengths of the dominant point in the x axis, the y axis and the z axis as
Figure BDA0002202823160000047
The sides of the cuboid are m2Dividing equally, making normal planes of coordinate axes of all equally divided points, and intersecting all normal planes to form second grid dominant points of the first stepStep grid, any second step grid node
Figure BDA0002202823160000051
Satisfies the formulas (13), (14) and (15). Determining tunneling stratum R according to stratum investigation datacThen, the points are put
Figure BDA0002202823160000052
Figure BDA0002202823160000053
Formula (8) was substituted to obtain β × corresponding to all second-step grid nodes.
Figure BDA0002202823160000054
Figure BDA0002202823160000055
Figure BDA0002202823160000056
β corresponding to the first step grid dominant point and β corresponding to all the second step grid nodes are collected, and the grid nodes corresponding to the first 10% of the maximum β are taken as initial selection optimization parameter points.
And sixthly, checking the tunneling rate.
The positive hob and the side hob are distributed on the cutter head, the positive hob is installed in parallel to the tunneling direction, and a non-zero side hob installation angle exists between the side hob and the tunneling direction. The effective thrust is approximately evenly distributed to each hob, the force parallel to the tunneling direction is F/N, the forces borne by the positive hob and the side hob and parallel to the tunneling direction are approximately evenly F/N, and the number of the positive hobs is NfcWith a number of edge cutters of NlcThe total number of the cutters is N, and the formula (16) shows. According to the balance of forces, the resultant force of the stratum resistance borne by the cutter head is equal to the resultant force of the effective thrust, the torque of the cutter head is equal to the resultant moment of the stratum resistance to the main shaft of the cutter head, and according to the Saint-Venn principle, the sum of the cutter head torques borne by the positive hob is approximate to that of the cutter headIs λ μ (FN)fcR/N) R/2, the cutter head torque can be approximate to the resultant torque of the occlusal force and the frictional resistance generated between the cutter and the rock on the face by the effective thrust, as shown in formula (17), R is the maximum installation radius of the positive hob, and R iswThe mounting radius of the w-th edge hob, thetawThe w-th installation angle of the edge roller cutter is shown, mu is the friction coefficient between steel and the stratum, and lambda is the empirical correction coefficient.
N＝Nfc+Nlc(16)
Figure BDA0002202823160000057
According to the formula (8), the tunneling rate to-be-determined value v is represented by the formula (18).
Figure BDA0002202823160000061
And substituting the cutter head rotating speed n, the soil bin pressure p and the effective thrust F of the primarily selected optimization parameter points into a formula (18) to obtain v corresponding to the points, calculating the v corresponding to each primarily selected optimization parameter point one by one, and taking the tunneling parameter of the primarily selected optimization parameter point when v is the maximum as the optimized tunneling parameter.
Further, in the step one: v, T, F, p and n are respectively in units of mm/min and 106Nm、kN、Bar、r/min。
Further, in the sixth step, λ is an empirical correction coefficient, and λ can be 0.24, 0.21, and 0.2 in a granite formation, an andesite formation, and a limestone formation, respectively.
Drawings
FIG. 1 is a scatter plot of equation coefficients for different formations under a finite sample;
fig. 2 shows the installation angle of the edge roller on the cutter head.
Detailed description of the invention
Aiming at the defects of the prior art, the invention provides a rapid method for optimizing the tunneling parameters of an earth pressure balance type shield tunneling machine, which comprises the following steps:
when the shield tunnels different types of stratums, namely tunnel face stratum distribution is changed, the shield tunneling speed and cutter head torque are changed accordingly. The tunneling speed is larger, the construction period is shorter, but the cutter head torque is also increased, and the cutter head torque is the most main control parameter of the integrity of equipment and reflects the working state of a cutter head main shaft. In order to meet the requirements of construction period and equipment completeness, balance needs to be obtained between the tunneling rate and the cutter head torque, so that the potential of the equipment is excavated to the maximum extent under the condition that the equipment is guaranteed to be complete, and the tunneling rate is effectively improved. Therefore, the shield tunneling parameter control optimization method based on the geological segmentation is provided, and the shield tunneling parameter dynamic optimization control based on the tunneling stratum parameters is realized.
Step one, establishing a statistical sample of tunneling efficiency
The tunneling efficiency β is defined as the tunneling rate v divided by the cutterhead torque T, as formula (1). T directly reflects the operating state of the equipment, and β reflects the relationship between the tunneling rate and the loading capacity of the cutterhead drive motor.
Figure BDA0002202823160000071
Combining tunneling parameters automatically recorded by a shield machine during trial tunneling of a tunnel and tunneling parameter data of a previous project, taking a tunneling rate v, a cutter torque T, an effective thrust F, a soil bin pressure p and a cutter rotating speed n which are actually measured at the same moment as a group of tunneling parameters, and taking β, F, p and n at each recording moment in a trial tunneling stage as a tunneling efficiency statistical sample;
v, T, F, p, n are in mm/min (millimeter per minute), 106Nm, kN, Bar, r/min (circles per minute).
Step two, obtaining a prediction equation of the tunneling efficiency for the limited sample through symbolic regression and linear regression to obtain the saturated uniaxial compressive strength R of the stratumcβ of the shield machine in different strata are respectively subjected to symbolic regression, the operation time of a computer running a symbolic regression algorithm in the background is less than 1 minute, and the fitting precision r2Fitting equations for β above 0.7 are candidate equations.
Counting candidate equations of each stratum, independent variable n2P, np, pF allCo-existing in at least one of the candidate equations of the respective formations, thus in n2P, np and pF are used as independent variables of the tunneling efficiency prediction equation under the limited sample, as shown in formula (2), linear regression is carried out again on the tunneling efficiency statistical sample through the tunneling efficiency prediction equation under the limited sample, and equation coefficients belonging to each stratum are obtained, and are shown in table 1;
β is the predicted value of tunneling efficiency under limited samples.
β*＝a1n2+a2p+a3np+a4pF+a5(2)
TABLE 1 prediction equation coefficient of tunneling efficiency of different saturated uniaxial compressive strength strata
Figure BDA0002202823160000081
Step three tunneling efficiency prediction equation continuation
Table 1 shows equation coefficients obtained by equation (2) under discrete finite sample conditions (the saturated uniaxial compressive strengths of the formation are 42MPa, 54MPa, 62MPa, 90MPa, and 118MPa, respectively), and the equation corresponding to the discrete equation coefficients under the finite sample conditions needs to be extended to obtain the equation coefficients suitable for general formations (the saturated uniaxial compressive strength is different from the sample value), that is, the equation can be applied to the cases when the saturated uniaxial compressive strength of the formation is not equal to 42MPa, 54MPa, 62MPa, 90MPa, and 118 MPa.
Taking coefficients of prediction equations of the tunneling efficiency in different tunneling strata under the limited samples obtained in the step one as dependent variables, and taking the saturated uniaxial compressive strength R of the tunneling strata ascLinear and nonlinear regression for the independent variables, i.e. for the scatter points in FIG. 1, the equation coefficients are obtained as a function of RcThe continuous function of the change is shown as formula (3), formula (4), formula (5), formula (6) and formula (7). Because the samples are all obtained from hard rocks, the continuation of the tunneling efficiency prediction equation is also in the hard rock range, namely Rc≥30MPa。
a1＝116.588-99.614Rc 0.038,Rc≥30MPa (3)
a2＝46.522-12.532Rc 0.313,Rc≥30MPa (4)
a3＝0.093Rc-4.414,Rc≥30MPa (5)
a4＝1.135×10-6Rc-6.976×10-5,Rc≥30MPa (6)
Figure BDA0002202823160000082
According to the formula (8), the formula (3), the formula (4), the formula (5), the formula (6) and the formula (7) are substituted into the formula (2), so that a tunneling efficiency prediction equation suitable for the soil pressure balance type shield tunneling hard rock is obtained, and β is a tunneling efficiency prediction value of the soil pressure balance type shield tunneling hard rock.
Figure BDA0002202823160000091
Step four, determining an optimized calculation formula and boundary conditions
The optimal calculation formula of the tunneling parameter is shown as a formula (9), pΩ、ΘΩ、ΨΩRespectively reflects the allowable fluctuation range of the soil bin pressure, the effective thrust and the cutter head rotating speed when the stratum is normally tunneled omega, namely [ pmin,pmax]、[Fmin,Fmax]、[nmin,nmax]。
Figure BDA0002202823160000092
Taking Shenzhen subway certain tunnel interval as an example, the saturated uniaxial compressive strength R of the tunneling stratumcAnd the pressure is 64.5MPa, the heading efficiency prediction equation corresponding to the stratum is obtained according to the formula (8) and is β ═ 0.115n2+0.346p+1.585np+3.448×10- 6pF +2.739, and according to engineering experience, during normal tunneling, the allowable fluctuation range p of the pressure of the soil bin in the stratumΩIs [0.05,0.4 ]]Allowable fluctuation range psi of cutter head rotation speedΩIs [0.8,2.2 ]]Effective thrust allowable fluctuation range thetaΩIs [6000,15000 ]]。
Step five-step dot method to calculate β larger value
Influence factors of the optimization problem are quantized, the boundary is clear, linear equal division is preferably adopted to obtain a mesh with the optimal value of the tunneling parameter to be selected, and step-by-step point taking method gradual operation is carried out on mesh nodes, so that influence of initial value selection on a calculation result is avoided.
P is to beΩ、ΘΩ、ΨΩAnd β x space domain is obtained as the domain defined on the mutually orthogonal x, y and z coordinate axes of the space rectangular coordinate system.
Establishing a first step grid, and respectively connecting pΩ、ΘΩ、ΨΩAs m1Dividing equally and making normal planes of coordinate axes of the equally divided points at all equally divided points, wherein the normal planes are intersected to form β -x space definition domain first step grid, and any one first step grid node
Figure BDA0002202823160000101
Satisfies the formulas (10), (11) and (12). Determining tunneling stratum R according to stratum investigation datacThen, the points are put
Figure BDA0002202823160000102
And (3) substituting the formula (8) to obtain β corresponding to each first step grid node in the spatial definition domain of β, arranging β corresponding to each first step grid node in the spatial definition domain of β in descending order of magnitude to obtain a first step grid node β number sequence, and taking the grid nodes corresponding to β which is the first 20% of the maximum in the first step grid nodes β as the first step grid dominant points.
Figure BDA0002202823160000103
Figure BDA0002202823160000104
Figure BDA0002202823160000105
Establishing a second step grid, and taking the dominant point of each first step grid as a body center and the projection lengths of the dominant point in the x axis, the y axis and the z axis as
Figure BDA0002202823160000106
The sides of the cuboid are m2Dividing equally, making normal planes of coordinate axes of all equally divided points, intersecting all normal planes to form second step grid of the first step grid dominant point, any one second step grid node
Figure BDA0002202823160000107
Satisfies the formulas (13), (14) and (15). Determining tunneling stratum R according to stratum investigation datacThen, the points are put
Figure BDA0002202823160000108
Figure BDA0002202823160000109
Formula (8) was substituted to obtain β × corresponding to all second-step grid nodes.
Figure BDA00022028231600001010
Figure BDA0002202823160000111
Figure BDA0002202823160000112
β corresponding to the first step grid dominant point and β corresponding to all the second step grid nodes are collected, and the grid nodes corresponding to the first 10% of the maximum β are taken as initial selection optimization parameter points.
Step six tunneling rate check
As shown in figure 2, a positive hob and an edge hob are distributed on the cutterhead, the positive hob is installed in parallel with the tunneling direction, and the edge hob does not exist between the tunneling direction and the edge hobAnd the installation angle of the edge hob is zero. The effective thrust is approximately evenly distributed to each hob, the force parallel to the tunneling direction is F/N, the forces borne by the positive hob and the side hob and parallel to the tunneling direction are approximately evenly F/N, and the number of the positive hobs is NfcWith a number of edge cutters of NlcThe total number of the cutters is N, and the formula (16) shows. According to the balance of forces, the resultant force of the stratum resistance borne by the cutter head is equal to the resultant force of the effective thrust, the torque of the cutter head is equal to the resultant moment of the stratum resistance to the main shaft of the cutter head, and the sum of the torques of the cutter head borne by the positive hob is approximate to lambda mu (FN) according to the Saint-Venn principlefcR/N) R/2, the cutter head torque can be approximate to the resultant torque of the occlusal force and the frictional resistance generated between the cutter and the rock on the face by the effective thrust, as shown in formula (17), R is the maximum installation radius of the positive hob, and R iswThe mounting radius of the w-th edge hob (the distance between the contact point between the hob and the tunnel face and the normal line of the cutterhead plane passing through the center of the cutterhead), thetawThe w-th installation angle of the edge hob is set, mu is the friction coefficient between steel and the stratum, lambda is the empirical correction coefficient, and lambda can be respectively 0.24, 0.21 and 0.2 in the granite stratum, the andesite stratum and the limestone stratum.
N＝Nfc+Nlc(16)
Figure BDA0002202823160000113
According to the formula (8), the tunneling rate to-be-determined value v is represented by the formula (18).
Figure BDA0002202823160000121
And substituting the cutter head rotating speed n, the soil bin pressure p and the effective thrust F of the primarily selected optimization parameter points into a formula (18) to obtain v corresponding to the points, calculating the v corresponding to each primarily selected optimization parameter point one by one, and taking the tunneling parameter of the primarily selected optimization parameter point when v is the maximum as the optimized tunneling parameter.

Claims (3)

1. A rapid method for optimizing the tunneling parameters of an earth pressure balance type shield tunneling machine is characterized by comprising the following steps:
step one, establishing a tunneling efficiency statistical sample
The tunneling efficiency β is defined as the tunneling rate v divided by the cutterhead torque T, as shown in the formula (1), wherein T directly reflects the running state of equipment, and β reflects the relationship between the tunneling rate and the loading capacity of a cutterhead driving motor;
Figure FDA0002473722560000011
combining tunneling parameters automatically recorded by a shield machine during trial tunneling of a tunnel and tunneling parameter data of a previous project, and taking a tunneling speed v, a cutter torque T, an effective thrust F, a soil bin pressure p and a cutter rotating speed n which are actually measured at the same moment as a same group of tunneling parameters;
step two, obtaining a tunneling efficiency prediction equation for the limited sample through symbolic regression and linear regression
Uniaxial compressive strength R to formation saturationcPerforming symbolic regression on the tunneling efficiency β of the shield tunneling machine in different strata respectively, and running the symbolic regression algorithm on the computer with the computation time of less than 1 minute and the fitting precision r2A fitting equation for β above 0.7 as a candidate equation;
counting candidate equations of each stratum, independent variable n2P, np, pF all co-exist in at least one of the candidate equations for each formation, thus in n2P, np and pF are used as independent variables of the tunneling efficiency prediction equation under the limited sample, as shown in formula (2), linear regression is carried out again on the tunneling efficiency statistical sample through the tunneling efficiency prediction equation under the limited sample, and equation coefficients belonging to each stratum are obtained;
β is a predicted value of the tunneling efficiency under a limited sample;
β*＝a1n2+a2p+a3np+a4pF+a5(2)
step three, tunneling efficiency prediction equation continuation
The equation coefficients in the formula (2) obtained under the condition of the limited samples are discrete equation coefficients for the limited kinds of strata, and the equations corresponding to the discrete equation coefficients need to be extended to obtain the equations suitable for the general strata;
taking coefficients of prediction equations of the tunneling efficiency in different tunneling strata under the limited samples obtained in the step one as dependent variables, and taking the saturated uniaxial compressive strength R of the tunneling strata ascLinear and nonlinear regression for independent variables to obtain equation coefficients with RcThe continuous function of the change is shown as formula (3), formula (4), formula (5), formula (6) and formula (7); because the samples are all obtained from hard rocks, the continuation of the tunneling efficiency prediction equation is also in the hard rock range, namely Rc≥30MPa；
a1＝116.588-99.614Rc 0.038,Rc≥30MPa (3)
a2＝46.522-12.532Rc 0.313,Rc≥30MPa (4)
a3＝0.093Rc-4.414,Rc≥30MPa (5)
a4＝1.135×10-6Rc-6.976×10-5,Rc≥30MPa (6)
Figure FDA0002473722560000021
According to the formula (8), the formula (3), the formula (4), the formula (5), the formula (6) and the formula (7) are replaced by the formula (2), so that a tunneling efficiency prediction equation suitable for the soil pressure balance type shield tunneling hard rock is obtained, and β is a tunneling efficiency prediction value of the soil pressure balance type shield tunneling hard rock;
Figure FDA0002473722560000022
step four, determining an optimized calculation formula and boundary conditions
The optimal calculation formula of the tunneling parameter is shown as a formula (9), pΩ、ΘΩ、ΨΩRespectively reflects the pressure of the soil bin, the effective thrust and the cutter when the stratum is normally tunneled at omegaAllowable fluctuation range of disk rotation speed, i.e. [ p ]min,pmax]、[Fmin,Fmax]、[nmin,nmax]；
Figure FDA0002473722560000023
Step five, β is calculated by a step-by-step point-taking method, and β is taken as the initial optimization parameter point corresponding to the top 10% of the maximum grid nodes
Influence factors of the optimization problem are quantized, the boundary is clear, linear equal division is adopted to obtain a mesh with the optimal value of the tunneling parameter to be selected, and step-by-step point-taking method gradual operation is carried out on mesh nodes;
p is to beΩ、ΘΩ、ΨΩRespectively as the definition domains on the orthogonal x, y and z coordinate axes of the space rectangular coordinate system to obtain β × space definition domains;
establishing a first step grid, and respectively connecting pΩ、ΘΩ、ΨΩAs m1Dividing equally and making normal planes of coordinate axes of the equally divided points at all equally divided points, wherein the normal planes are intersected to form β -x space definition domain first step grid, and any one first step grid node (p)i1，Fj1，nk1) Satisfies formula (10), formula (11), formula (12); determining tunneling stratum R according to stratum investigation datacThen, point (p) is addedi1，Fj1，nk1) Substituting formula (8) to obtain β corresponding to each first step grid node in β space definition, arranging β corresponding to each first step grid node in β space definition in descending order to obtain β number sequence of first step grid nodes, and taking β corresponding to the maximum first 20% of β in β number of first step grid nodes as dominant point of first step grid;
Figure FDA0002473722560000031
Figure FDA0002473722560000032
Figure FDA0002473722560000033
establishing a second step grid, and taking the dominant point of each first step grid as a body center and the projection lengths of the dominant point in the x axis, the y axis and the z axis as
Figure FDA0002473722560000034
The sides of the cuboid are m2Dividing equally, making normal planes of coordinate axes of all equally divided points, intersecting all normal planes to form second step grid of the first step grid dominant point, any one second step grid node
Figure FDA0002473722560000035
Satisfies formula (13), formula (14), formula (15); determining tunneling stratum R according to stratum investigation datacThen, the points are put
Figure FDA0002473722560000036
Figure FDA0002473722560000041
Substituting the formula (8) to obtain β corresponding to all second step grid nodes;
Figure FDA0002473722560000042
Figure FDA0002473722560000043
Figure FDA0002473722560000044
β corresponding to the first-step grid dominant point and β corresponding to all second-step grid nodes are collected, and the grid nodes corresponding to the first 10% of the maximum β are taken as initial selection optimization parameter points;
sixthly, checking the tunneling rate
A positive hob and an edge hob are distributed on the cutterhead, the positive hob is installed in parallel to the tunneling direction, and an edge hob installation angle which is not zero exists between the edge hob and the tunneling direction; the effective thrust is approximately evenly distributed to each hob, the force parallel to the tunneling direction is F/N, the forces borne by the positive hob and the side hob and parallel to the tunneling direction are approximately evenly F/N, and the number of the positive hobs is NfcWith a number of edge cutters of NlcThe total number of the cutters is N, and the formula (16) shows; according to the balance of forces, the resultant force of the stratum resistance borne by the cutter head is equal to the resultant force of the effective thrust, the torque of the cutter head is equal to the resultant moment of the stratum resistance to the main shaft of the cutter head, and the sum of the torques of the cutter head borne by the positive hob is approximate to lambda mu (FN) according to the Saint-Venn principlefcR/N) R/2, the cutter head torque can be approximate to the resultant torque of the occlusal force and the frictional resistance generated between the cutter and the rock on the face by the effective thrust, as shown in formula (17), R is the maximum installation radius of the positive hob, and R iswThe mounting radius of the w-th edge hob, thetawSetting the installation angle of the w-th edge hob, wherein mu is the friction coefficient between steel and the stratum, and lambda is the empirical correction coefficient;
N＝Nfc+Nlc(16)
Figure FDA0002473722560000045
according to the formula (8), the tunneling rate to-be-determined value v is shown as the formula (18);
Figure FDA0002473722560000051
and substituting the cutter head rotating speed n, the soil bin pressure p and the effective thrust F of the primarily selected optimization parameter points into a formula (18) to obtain v corresponding to the points, calculating the v corresponding to each primarily selected optimization parameter point one by one, and taking the tunneling parameter of the primarily selected optimization parameter point when v is the maximum as the optimized tunneling parameter.
2. The optimized earth pressure balanced shield as claimed in claim 1The method for quickly tunneling parameters of the mechanism is characterized by comprising the following steps: v, T, F, p and n are respectively in units of mm/min and 106Nm、kN、Bar、r/min。
3. The rapid method for optimizing the tunneling parameters of the earth pressure balance type shield tunneling machine according to claim 1, characterized in that in the sixth step: lambda is an empirical correction coefficient, and can be respectively 0.24, 0.21 and 0.2 in a granite stratum, an andesite stratum and a limestone stratum.
CN201910871045.9A 2019-09-16 2019-09-16 Rapid method for optimizing tunneling parameters of earth pressure balanced type shield tunneling machine Active CN110704963B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910871045.9A CN110704963B (en) 2019-09-16 2019-09-16 Rapid method for optimizing tunneling parameters of earth pressure balanced type shield tunneling machine

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910871045.9A CN110704963B (en) 2019-09-16 2019-09-16 Rapid method for optimizing tunneling parameters of earth pressure balanced type shield tunneling machine

Publications (2)

Publication Number Publication Date
CN110704963A CN110704963A (en) 2020-01-17
CN110704963B true CN110704963B (en) 2020-06-26

Family

ID=69196193

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910871045.9A Active CN110704963B (en) 2019-09-16 2019-09-16 Rapid method for optimizing tunneling parameters of earth pressure balanced type shield tunneling machine

Country Status (1)

Country Link
CN (1) CN110704963B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111946398B (en) * 2020-08-17 2021-06-29 南京工业大学 Composite stratum shield tunneling efficiency field prediction calculation method
CN112196559B (en) * 2020-09-30 2021-08-27 山东大学 TBM operation parameter optimization method based on optimal tunneling speed and optimal cutter consumption
CN113361043B (en) * 2021-06-25 2022-05-06 天津大学 Method and system for predicting specific energy of cutter head of hard rock tunnel boring machine
CN113374488B (en) * 2021-07-28 2022-09-27 中国铁建重工集团股份有限公司 Earth pressure balance shield machine guiding control method and device and readable storage medium
CN117744450B (en) * 2024-02-08 2024-04-16 湖南大学 Simulation method for shield tunneling process

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09184394A (en) * 1995-12-28 1997-07-15 Shimizu Corp Shield machine
CN107180016A (en) * 2017-05-23 2017-09-19 南京工业大学 Hierarchical summation method for predicting hob abrasion loss by using abrasion specific loss index
CN110069893A (en) * 2019-05-09 2019-07-30 中铁工程服务有限公司 A kind of prediction technique of the shield machine boring parameter based on polynomial regression

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105631150A (en) * 2016-01-05 2016-06-01 石家庄铁道大学 Optimization method of shield excavation parameters under condition of compound stratum
CN107248026A (en) * 2017-05-23 2017-10-13 南京工业大学 Quantitative method for predicting shield tunneling parameters by using basic quality indexes of equivalent rock mass
CN109711079B (en) * 2019-01-03 2020-07-31 天津大学 TBM tunneling total thrust determining method and system

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH09184394A (en) * 1995-12-28 1997-07-15 Shimizu Corp Shield machine
CN107180016A (en) * 2017-05-23 2017-09-19 南京工业大学 Hierarchical summation method for predicting hob abrasion loss by using abrasion specific loss index
CN110069893A (en) * 2019-05-09 2019-07-30 中铁工程服务有限公司 A kind of prediction technique of the shield machine boring parameter based on polynomial regression

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
复杂地层盾构掘进速率和刀盘扭矩预测模型及其地层适应性研究;张志奇等;《隧道建设》;20161231;第36卷(第12期);1449-1455 *
长沙地铁下穿湘江土压平衡盾构隧道掘进参数研究;褚东升;《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20130228;C034-315 *

Also Published As

Publication number Publication date
CN110704963A (en) 2020-01-17

Similar Documents

Publication Publication Date Title
CN110704963B (en) Rapid method for optimizing tunneling parameters of earth pressure balanced type shield tunneling machine
Wang et al. Disc cutter wear prediction for a hard rock TBM cutterhead based on energy analysis
Nelson TBM performance analysis with reference to rock properties
CN109307493B (en) Real-time detection method for abrasion condition of constant-section disc cutter of hard rock tunnel boring machine
Torabi et al. Study of the influence of geotechnical parameters on the TBM performance in Tehran–Shomal highway project using ANN and SPSS
Zhang et al. Theoretical prediction of wear of disc cutters in tunnel boring machine and its application
Frough et al. Application of RMR for estimating rock-mass–related TBM utilization and performance parameters: a case study
CN110378069A (en) Method for predicting accumulated mass loss rate of cutter of tunneling machine
Gao et al. Evolution of local mine stiffness with mining process: insight from physical and numerical modeling
Koukoutas et al. Settlements due to single and twin tube urban EPB shield tunnelling
CN109117556B (en) Shield propulsion distance prediction method based on shield cutter head cutter partition cutting performance
CN108984817A (en) A kind of TBM tool abrasion real time evaluating method
Ocak et al. Performance prediction of impact hammer using ensemble machine learning techniques
Bilgin et al. The effect of rock weathering and transition zones on the performance of an EPB-TBM in complex geology near Istanbul, Turkey
She et al. Prediction model for disc cutter wear during hard rock breaking based on plastic removal abrasiveness mechanism
CN113435056B (en) Shield utilization rate prediction and operation parameter optimization method and system based on SVR and PSO
Tang et al. Deviation correction strategy for the earth pressure balance shield based on shield-soil interactions
Frough et al. Study of the correlation between RMR and TBM downtimes
Jain et al. Estimation of the performance of the tunnel boring machine (TBM) using uniaxial compressive strength and rock mass rating classification (RMR)–A case study from the Deccan traps, India
Samadi et al. Developing GEP technique for prediction of EPB-TBM performance in limestone strata
CN115774912A (en) Mechanism and data combined driving TBM hob abrasion prediction method
Stavropoulou et al. Characterization of rock masses based on geostatistical joint mapping and rock boring operations
Sun et al. Evaluation of TBM cutterhead vibration under complicated condition
CN113158561A (en) TBM operation parameter optimization method and system suitable for various rock mass conditions
Bitetti et al. On the effects of Line 6 tunnel excavation in Naples

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
CP02 Change in the address of a patent holder

Address after: 210000 rooms 10C and 10d, building 02, No.77, alfalfa Garden Street, Qinhuai District, Nanjing City, Jiangsu Province

Patentee after: NANJING KENTOP CIVIL ENGINEERING TECHNOLOGY Co.,Ltd.

Patentee after: NANJING TECH University

Address before: 210000 105, Jia Dong Cun, Yuhuatai oil District, Nanjing, Jiangsu.

Patentee before: NANJING KENTOP CIVIL ENGINEERING TECHNOLOGY Co.,Ltd.

Patentee before: NANJING TECH University

CP02 Change in the address of a patent holder
CP02 Change in the address of a patent holder

Address after: 210000 105, Jia Dong Cun, Yuhuatai oil District, Nanjing, Jiangsu.

Patentee after: NANJING KENTOP CIVIL ENGINEERING TECHNOLOGY Co.,Ltd.

Patentee after: NANJING TECH University

Address before: 210000 rooms 10C and 10d, building 02, No.77, alfalfa Garden Street, Qinhuai District, Nanjing City, Jiangsu Province

Patentee before: NANJING KENTOP CIVIL ENGINEERING TECHNOLOGY Co.,Ltd.

Patentee before: NANJING TECH University

CP02 Change in the address of a patent holder