CN-116738537-B - Shield hob structure optimization method, system, electronic equipment and storage medium

CN116738537BCN 116738537 BCN116738537 BCN 116738537BCN-116738537-B

Abstract

The invention discloses a shield hob structure optimization method, a shield hob structure optimization system, electronic equipment and a storage medium, and relates to the field of shield construction. The method comprises the steps of constructing a rock HJC constitutive model through physical parameters and constitutive parameters of the rock, selecting sensitive factors for hob structure optimization, constructing a hob model, carrying out hob breaking numerical simulation based on finite element software LS-DYNA to obtain counter force and crushed rock volume when a hob is cut, constructing a rock breaking specific energy prediction model and a hob abrasion prediction model by adopting a XGBoost algorithm, carrying out multi-objective optimization by taking the rock breaking specific energy prediction model and the hob abrasion prediction model as objective optimization functions and the sensitive factors as decision variables by using a NSGAII method to obtain an optimal solution, enabling the optimal solution to be the optimal hob parameters, and determining the optimal hob structure based on the optimal hob parameters. The method realizes the prediction and multi-objective optimization of the rock breaking efficiency and abrasion of the hob in a low-cost mode, and provides guidance for hob selection and tunneling parameter selection before shield construction.

Inventors

YANG WEIMIN
SONG XUAN
WANG MEIXIA
ZHANG ZHIYUAN
BAI YIFAN
WANG SENWEI
TIAN CONG
LIU LANG

Assignees

山东大学

Dates

Publication Date: 20260505
Application Date: 20230605

Claims (8)

1. The shield hob structure optimization method is characterized by comprising the following steps: Constructing a rock HJC constitutive model through physical parameters of the rock and constitutive parameters, wherein the physical parameters comprise rock density, compressive strength, elastic modulus, poisson ratio, shear modulus and bulk modulus, and the constitutive parameters comprise strain rate effect parameters, limit surface parameters, pressure parameters and damage parameters; selecting sensitive factors for hob structure optimization, and constructing a hob model, wherein the sensitive factors comprise a blade fillet, a blade width and a cutter penetration; the hob model and the rock HJC constitutive model are imported into finite element software LS-DYNA, hob rock breaking numerical simulation is conducted, and counter force and rock crushing volume when hob cutting are obtained, wherein the counter force when hob cutting comprises rolling force, normal force and lateral force; Based on the sensitive factors, the counterforce applied by the hob during cutting and the crushed rock volume, constructing a rock breaking specific energy prediction model and a hob abrasion prediction model by adopting an XG Boost algorithm; Taking the rock breaking specific energy prediction model and the hob abrasion prediction model as target optimization functions, taking the sensitive factors as decision variables, and performing multi-target optimization by using a NSGAII method to obtain an optimal solution, wherein the optimal solution is an optimal hob parameter; And determining an optimal hob structure based on the optimal hob parameters.
2. The method for optimizing the structure of the shield hob according to claim 1, wherein the rock HJC constitutive model comprises an intensity model, a damage evolution equation and a state equation; The expression of the intensity model is as follows: Wherein, the In order to characterize the equivalent stress, In order to characterize the pressure, the pressure is, To characterize strain rate, A, B, N and All are limit surface parameters, C is a strain rate effect parameter, and D is a damage parameter; The expression of the damage evolution equation is as follows: Wherein, the 、 Calculating an equivalent plastic strain delta and a plastic volume strain delta for the cell within the cycle; For the equivalent plastic strain and equivalent plastic volume strain at the current calculation step, The tensile strength is normalized for the material and, 、 Is a material damage parameter; is the minimum plastic strain at which the material breaks; The state equation comprises a state equation of an elastic compression stage, a state equation of a compaction deformation stage and a state equation of a deformation stage after compaction; the expression of the state equation of the elastic compression stage is as follows: ; The expression of the state equation of the compaction deformation stage is as follows: The expression of the state equation of the deformation stage after compaction is as follows: ; Wherein, the Is the water purifying pressure, K is the bulk modulus, mu is the bulk strain, T is the maximum tensile water purifying pressure of the material, Is the elastic limit water purification pressure, In order to compact the limit net water pressure, For a volumetric strain corresponding to the elastic limit, To limit the water purifying pressure The corresponding volume strain is used to determine, For the corresponding purified water pressure of the volume deformation before unloading, F is an unloading proportionality coefficient, Is the corrected volume strain; 、、 Is a pressure constant.
3. The shield hob structure optimization method according to claim 1, wherein the expression of the rock breaking specific energy prediction model is as follows: The expression of the hob abrasion prediction model is as follows: Wherein, the In order to break the rock in a specific energy, For the abrasion of the hob cutter, 、 Respectively rolling force and normal force, 、 V is the penetration degree of the hob, the rolling distance and the rock breaking volume, I is the energy abrasion rate, Is the coefficient of friction.
4. The shield hob configuration optimization method according to claim 1, characterized in that before constructing a rock breaking specific energy prediction model and a hob abrasion prediction model by using an XG Boost algorithm based on the sensitivity factor, the reaction force applied by the hob during cutting and the crushed rock volume, the method further comprises: and adopting a noise reduction algorithm to reduce the noise of the counter force applied to the hob during cutting.
5. A shield cutter structure optimization system, comprising: The rock HJC constitutive model construction module is used for constructing a rock HJC constitutive model through physical parameters of the rock and constitutive parameters, wherein the physical parameters comprise rock density, compressive strength, elastic modulus, poisson ratio, shear modulus and bulk modulus, and the constitutive parameters comprise strain rate effect parameters, limit surface parameters, pressure parameters and damage parameters; The hob model construction module is used for selecting sensitive factors for hob structure optimization and constructing a hob model, wherein the sensitive factors comprise a blade fillet, a blade width and a cutter penetration; The hob breaking numerical simulation module is used for guiding the hob model and the rock HJC constitutive model into finite element software LS-DYNA, and carrying out hob breaking numerical simulation to obtain counter force and crushed rock volume when the hob is cut, wherein the counter force comprises rolling force, normal force and lateral force when the hob is cut; the prediction model construction module is used for constructing a rock breaking specific energy prediction model and a hob abrasion prediction model by adopting an XG Boost algorithm based on the sensitive factors, the counter force applied by the hob during cutting and the crushed rock volume; the optimizing module is used for carrying out multi-objective optimization by taking the rock breaking specific energy prediction model and the hob abrasion prediction model as objective optimizing functions and the sensitive factors as decision variables and using a NSGAII method to obtain an optimal solution, wherein the optimal solution is an optimal hob parameter; And determining an optimal hob structure based on the optimal hob parameters.
6. The shield tunneling cutter structure optimization system of claim 5, further comprising: And the noise reduction module is used for carrying out noise reduction treatment on counter force applied to the hob during cutting by adopting a noise reduction algorithm.
7. An electronic device comprising a memory for storing a computer program and a processor that runs the computer program to cause the electronic device to perform the shield hob configuration optimization method according to any one of the claims 1-4.
8. A computer-readable storage medium, characterized in that it stores a computer program which, when executed by a processor, implements the shield hob structure optimization method according to any one of the claims 1-4.

Description

Shield hob structure optimization method, system, electronic equipment and storage medium Technical Field The invention relates to the technical field of shield construction, in particular to a shield hob structure optimization method, a system, electronic equipment and a storage medium. Background The shield method has become a main method for urban tunnel construction due to the advantages of high mechanization degree, high construction efficiency, small environmental disturbance and the like. The disc-type hob is used as a direct rock breaking tool of the shield tunneling machine, and has important influence on the efficiency and safety of shield tunnel construction. The factors such as the cutter width, the penetration, the cutter edge angle and the like of the hob obviously influence the rock breaking efficiency, and simultaneously have great influence on the abrasion of the hob. Most of traditional research methods are based on related preconditions due to the limitation of test amount, and considered factors are limited, so that partial rules are often obtained, and the actual requirements of engineering are difficult to meet. Meanwhile, in the aspect of optimization research, system consideration on hob rock breaking efficiency and hob abrasion is absent. In application, excessive pursuit of rock breaking efficiency may cause the increase of hob abrasion, frequent tool changing adversely affects the construction progress, and the construction cost is increased. Disclosure of Invention The invention aims to provide a shield hob structure optimization method, a system, electronic equipment and a storage medium, which are used for solving the problem that the prior art lacks consideration of hob rock breaking efficiency and hob abrasion. In order to achieve the above object, the present invention provides the following solutions: a shield hob structure optimization method comprises the following steps: Constructing a rock HJC constitutive model through physical parameters of the rock and constitutive parameters, wherein the physical parameters comprise rock density, compressive strength, elastic modulus, poisson ratio, shear modulus and bulk modulus, and the constitutive parameters comprise strain rate effect parameters, limit surface parameters, pressure parameters and damage parameters; selecting sensitive factors for hob structure optimization, and constructing a hob model, wherein the sensitive factors comprise a blade fillet, a blade width and a cutter penetration; the hob model and the rock HJC constitutive model are imported into finite element software LS-DYNA, hob rock breaking numerical simulation is conducted, and counter force and rock crushing volume when hob cutting are obtained, wherein the counter force when hob cutting comprises rolling force, normal force and lateral force; Based on the sensitive factors, the counterforce applied by the hob during cutting and the crushed rock volume, constructing a rock breaking specific energy prediction model and a hob abrasion prediction model by adopting an XG Boost algorithm; Taking the rock breaking specific energy prediction model and the hob abrasion prediction model as target optimization functions, taking the sensitive factors as decision variables, and performing multi-target optimization by using a NSGAII method to obtain an optimal solution, wherein the optimal solution is an optimal hob parameter; And determining an optimal hob structure based on the optimal hob parameters. Optionally, the rock HJC constitutive model comprises an intensity model, a damage evolution equation and a state equation; The expression of the intensity model is as follows: Wherein σ * is the characterization equivalent stress, p * is the characterization pressure, For characterizing strain rate, A, B, N and S MAX are both limit surface parameters, C is strain rate effect parameter, and D is damage parameter; The expression of the damage evolution equation is as follows: wherein Δε p、Δμp is the equivalent plastic strain increase and plastic volume strain increase of a cell within a calculated cycle; For the equivalent plastic strain and equivalent plastic volume strain under the current calculation step, T * is the normalized tensile strength of the material, D 1、D2 is the damage parameter of the material, EF MIN is the minimum plastic strain when the material is damaged; The state equation comprises a state equation of an elastic compression stage, a state equation of a compaction deformation stage and a state equation of a deformation stage after compaction; the expression of the state equation of the compression stage is as follows: p=Kμ;-T(1-D)≤p≤pc; The expression of the state equation of the compaction deformation stage is as follows: p=p0-[(1-F)K+FK1](μ0-μ)) The expression of the state equation of the deformation stage after compaction is as follows: Wherein p is the water purifying pressure, K is the bulk modulus, mu is the bulk strain, T is the maximum t