Hobbing cutter rock breaking identification method, device and system based on semi-analytical method

Abstract

The invention relates to a hob rock breaking identification method, device and system based on a semi-analytical method, and relates to the technical field of shield cutter rock breaking. The method comprises the following steps: acquiring a cutter model and an initial rock model of the shield tunneling machine; simulating a scene in which the cutter cuts the rock by adopting the cutter model and the initial rock model to obtain a contact area of the cutter model and the initial rock model; dispersing the rock model in the first area into cubic units; assigning material attributes to the cube units, and calculating influence coefficients of the cube units; updating the load value of the cubic unit according to the influence coefficient, and calculating the total stress of the cubic unit according to the load value; determining an equivalent characteristic strain value of the first rock model and a convergence state of the equivalent characteristic strain value; and under the condition that the convergence state meets a preset condition, determining the number of damage units in the first rock model, and determining the damage degree of the first rock model according to the number. According to the scheme, the accuracy of pipeline rock breaking identification is improved.

Description

Technical Field

[0001] The present disclosure relates to the technical field of shield cutter rock breaking, and in particular to a semi-analytical method-based method, device and system for identifying rock breaking by a roller cutter. Background Art

[0002] In related technologies, natural rocks have inherent characteristics such as discontinuity, anisotropy, non-uniformity, and inelasticity. In addition, a large number of cracks are generated inside the rock when the cutter is broken, which often leads to the failure of the continuous medium mechanics assumption to be met. Therefore, it is very difficult to use traditional finite element numerical methods to simulate the cutter rock breaking process, which leads to poor accuracy of the simulation results and difficulty in meeting actual needs. Summary of the Invention

[0003] In order to overcome the problems existing in the related art, the present disclosure provides a method, device and system for identifying rock breaking by a roller cutter based on a semi-analytical method.

[0004] According to a first aspect of an embodiment of the present disclosure, a method for identifying rock breaking by a roller cutter based on a semi-analytical method is provided, comprising:

[0005] Obtain the shield machine's cutter model and initial rock model;

[0006] The tool model and the initial rock model are used to simulate a scene of the tool cutting the rock, and a contact area between the tool model and the initial rock model is obtained;

[0007] Determine a first region based on the contact region, and discretize the rock model in the first region into cubic units; the first region includes the contact region and the neighborhood of the contact region;

[0008] For each cubic unit, assign material attributes to the cubic unit to obtain an updated first rock model, and calculate an influence coefficient of each cubic unit according to the material attributes;

[0009] updating the load value of the cubic unit according to the influence coefficient, and calculating the total stress of the cubic unit according to the load value;

[0010] determining an equivalent characteristic strain value of the first rock model and a convergence state of the equivalent characteristic strain value according to the total stress of each cubic unit;

[0011] When the convergence state satisfies a preset condition, the number of damaged units in the first rock model is determined, and the damage degree of the first rock model is determined according to the number.

[0012] According to a second aspect of an embodiment of the present disclosure, a device for identifying rock breaking by a roller cutter based on a semi-analytical method is provided, comprising:

[0013] An acquisition unit, used to acquire a cutter model and an initial rock model of a shield machine;

[0014] a simulation unit, configured to simulate a scene of a tool cutting a rock using the tool model and the initial rock model, and obtain a contact area between the tool model and the initial rock model;

[0015] Discrete units are used to determine a first area according to the contact area, and discretize the rock model in the first area into cubic units; the first area includes the contact area and the neighborhood of the contact area;

[0016] an updating unit, configured to assign material properties to each cubic unit to obtain an updated first rock model, and calculate an influence coefficient of each cubic unit according to the material properties;

[0017] a calculation unit, configured to update the load value of the cubic unit according to the influence coefficient, and calculate the total stress of the cubic unit according to the load value;

[0018] a first determining unit, configured to determine an equivalent characteristic strain value of the first rock model and a convergence state of the equivalent characteristic strain value according to the total stress of each cubic unit;

[0019] The second determining unit is configured to determine the number of damaged units in the first rock model when the convergence state satisfies a preset condition, and determine the damage degree of the first rock model according to the number.

[0020] According to a third aspect of an embodiment of the present disclosure, an electronic device includes: a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the method described in any one of the first aspects is implemented.

[0021] According to a fourth aspect of an embodiment of the present disclosure, a computer-readable storage medium is provided, on which a computer program is stored. When the computer program is executed by a processor, the method according to any one of the first aspects is implemented.

[0022] According to a fifth aspect of an embodiment of the present disclosure, a computer program product is provided, comprising a computer program, wherein the computer program implements the method as described in any one of the first aspects when executed by a processor.

[0023] The technical solution provided by the embodiments of the present disclosure may include the following beneficial effects: by obtaining the tool model and initial rock model of the shield machine; using the tool model and the initial rock model to simulate the scene of the tool cutting the rock, and obtaining the contact area between the tool model and the initial rock model; determining the first area according to the contact area, and discretizing the rock model in the first area into cubic units; for each cubic unit, assigning material properties to the cubic unit to obtain an updated first rock model, and calculating the influence coefficient of each cubic unit according to the material properties; updating the load value of the cubic unit according to the influence coefficient, and calculating the total stress of the cubic unit according to the load value; determining the equivalent characteristic strain value of the first rock model and the convergence state of the equivalent characteristic strain value according to the total stress of each cubic unit; when the convergence state meets the preset conditions, determining the number of damaged units in the first rock model, and determining the damage degree of the first rock model according to the number, thereby improving the simulation accuracy of the rock model, and then improving the accuracy of pipeline rock breaking identification.

[0024] It is to be understood that the foregoing general description and the following detailed description are exemplary and explanatory only and are not restrictive of the disclosure. BRIEF DESCRIPTION OF THE DRAWINGS

[0025] The accompanying drawings, which are incorporated in and constitute a part of this specification, illustrate embodiments consistent with the invention and, together with the description, serve to explain the principles of the invention.

[0026] Figure 1 The present invention is a flowchart showing a method for identifying rock breaking by a roller cutter based on a semi-analytical method according to an exemplary embodiment.

[0027] Figure 2 The block diagram of a device for identifying rock breaking by a roller cutter based on a semi-analytical method is shown according to an exemplary embodiment.

[0028] Figure 3 The present invention is a block diagram showing an apparatus for a method for identifying rock breaking by a roller cutter based on a semi-analytical method according to an exemplary embodiment. DETAILED DESCRIPTION

[0029] Exemplary embodiments will be described in detail herein, examples of which are illustrated in the accompanying drawings. In the following description, when referring to the drawings, like numbers in different figures represent like or similar elements unless otherwise indicated. The embodiments described in the following exemplary embodiments are not intended to represent all possible embodiments consistent with the present invention. Rather, they are merely examples of apparatus and methods consistent with certain aspects of the present invention, as detailed in the appended claims.

[0030] The terms used in the present disclosure are only for the purpose of describing specific embodiments and are not intended to limit the present disclosure. The singular forms "a", "an" and "the" used in the present disclosure and the appended claims are also intended to include the plural forms unless the context clearly indicates otherwise.

[0031] It should be understood that although the terms first, second, third, etc. may be used to describe various information in the embodiments of the present disclosure, such information should not be limited to these terms. These terms are only used to distinguish information of the same type from each other. For example, without departing from the scope of the embodiments of the present disclosure, the first information may also be referred to as the second information, and similarly, the second information may also be referred to as the first information. Depending on the context, the words "if" and "if" as used herein may be interpreted as "at the time of" or "when" or "in response to determining".

[0032] Furthermore, the various forms of processes shown in the embodiments of this disclosure may be used to reorder, add, or delete steps. For example, the steps described in this application may be performed in parallel, sequentially, or in a different order, as long as the desired results of the technical solutions disclosed in this application can be achieved. This is not a limitation herein.

[0033] In related technologies, natural rocks have inherent characteristics such as discontinuity, anisotropy, non-uniformity, and inelasticity. In addition, a large number of cracks are generated inside the rock when the cutter is broken, which often leads to the failure of the continuous medium mechanics assumption to be met. Therefore, it is very difficult to use traditional finite element numerical methods to simulate the cutter rock breaking process, which leads to poor accuracy of the simulation results and difficulty in meeting actual needs.

[0034] In order to solve the above problems, the present disclosure provides a method, device and system for identifying rock breaking by a roller cutter based on a semi-analytical method, which obtains a tool model and an initial rock model of a shield machine; uses the tool model and the initial rock model to simulate the scene of the tool cutting rock to obtain the contact area between the tool model and the initial rock model; determines a first area based on the contact area, and discretizes the rock model in the first area into cubic units; assigns material properties to each cubic unit to obtain an updated first rock model, and calculates the influence coefficient of each cubic unit based on the material properties; updates the load value of the cubic unit based on the influence coefficient, and calculates the total stress of the cubic unit based on the load value; determines the equivalent characteristic strain value of the first rock model and the convergence state of the equivalent characteristic strain value based on the total stress of each cubic unit; when the convergence state meets the preset conditions, determines the number of damaged units in the first rock model, and determines the damage degree of the first rock model based on the number, thereby improving the simulation accuracy of the rock model and further improving the accuracy of pipeline rock breaking identification.

[0035] Figure 1 FIG. 1 is a flow chart showing a method for identifying rock breaking by a roller cutter based on a semi-analytical method according to an exemplary embodiment. Figure 1 As shown, it should be noted that the semi-analytical method based on the roller cutter rock breaking identification method of the embodiment of the present application is applied to the semi-analytical method based on the roller cutter rock breaking identification device. Figure 1 As shown, the method may include the following steps:

[0036] Step 101: Obtain a cutter model and an initial rock model of a shield machine.

[0037] It should be noted that both the tool model and the initial rock model are three-dimensional models.

[0038] Step 102 : Using the tool model and the initial rock model, simulate the scene of the tool cutting the rock to obtain the contact area between the tool model and the initial rock model.

[0039] Step 103: determine a first region according to the contact region, and discretize the rock model in the first region into cubic units.

[0040] The first area includes the contact area and the neighborhood of the contact area.

[0041] In one embodiment, the area coordinates of the neighborhood may be determined based on the position information of the contact area and a preset distance value.

[0042] Step 104 : assigning material attributes to each cubic unit to obtain an updated first rock model, and calculating the influence coefficient of each cubic unit according to the material attributes.

[0043] In some embodiments of the present application, step 104 may specifically include:

[0044] Assigning different material properties to each cubic unit to obtain an updated first rock model;

[0045] The influence coefficients are calculated for each cubic element based on the material properties; the influence coefficients include the inhomogeneity coefficient.

[0046] In one embodiment, different material properties can be assigned to each element, so that the Young's modulus and tensile / compressive strength of the element follow a Weibull probability distribution to characterize rock heterogeneity. A truncated Mohr-Coulomb criterion is used as the criterion for brittle fracture of the rock mass, reducing the Young's modulus of failed elements to zero. Damaged elements are treated as impurities and processed using the equivalent inclusion method. The same stress value is applied to each element and remains unchanged during the calculation process to account for the influence of confining pressure.

[0047] As an example, the Young's modulus and tensile / compressive strength of an element are calculated using the following formulas:

[0048]

[0049] Where u is a rock mechanical property, such as elastic modulus, tensile strength, etc.; u0 is the mean value of a rock mechanical property; and m is the shape parameter of the probability density function, which is used here as the heterogeneity coefficient.

[0050] As an example, the following formula is used as the criterion for brittle fracture of rock mass:

[0051]

[0052] Among them, σ I is the first principal stress of the element, σ III is the third principal stress of the element; σ C is the compressive strength of rock, σ T is the tensile strength of rock; φ is the internal friction angle of rock.

[0053] Step 105: update the load value of the cube element according to the influence coefficient, and calculate the total stress of the cube element according to the load value.

[0054] In some embodiments of the present application, step 105 may specifically include:

[0055] The conjugate gradient method is used to calculate the first pressure distribution data of the contact area according to the load value;

[0056] The second pressure distribution data of the neighborhood of the contact area is calculated based on the influence coefficient using the Fourier transform method;

[0057] The total stress of the cubic element is calculated based on the first pressure distribution data and the second pressure distribution data.

[0058] In one embodiment, the surface displacement caused by the rock damage zone can be calculated using influence coefficients combined with fast Fourier transforms. Surface deformation affects the contact pressure distribution, thereby changing the subsurface stress and inducing new rock damage. This contact-damage coupling effect is calculated through an iterative method until the change between two consecutive steps is less than a given threshold.

[0059] In one embodiment, the contact pressure distribution may be determined using the following formula:

[0060]

[0061] In the above formula: Δx, Δy are the side lengths of the surface element; α, β are the surface element numbers; Nx, Ny are the number of elements along the x and y directions; p is the contact pressure; W is the total normal force applied to the hob; h is the total gap between the contact surfaces; h 0 is the initial gap between the contact surfaces; u3 and is the surface displacement caused by contact pressure and impurities. The JFNK (Jacoian-Free Newton Krylov) method is used to solve the nonlinear equations of the equivalent inclusion method.

[0062] Step 106: determining the equivalent characteristic strain value of the first rock model and the convergence state of the equivalent characteristic strain value according to the total stress of each cubic unit.

[0063] Step 107: When the convergence state satisfies a preset condition, the number of damaged units in the first rock model is determined, and the damage degree of the first rock model is determined according to the number.

[0064] It can be understood that after the rock model is updated in the above manner, the rock model has a higher accuracy, and in this case, the number of damaged units determined has a higher accuracy.

[0065] In one embodiment, when the equivalent characteristic strain value converges, the number of damaged units in the first rock model is determined, and the damage degree of the first rock model is determined according to the number.

[0066] As an example, the above quantity may be compared with a plurality of preset threshold intervals to determine the threshold interval to which the quantity belongs and the damage level corresponding to the threshold interval.

[0067] In some embodiments of the present application, the method may further include: when the convergence state does not meet the preset conditions, updating the displacement state information of the tool model, returning to execute the step of determining the first area according to the contact area based on the displacement state information, and discretizing the rock model in the first area into cubic units until the convergence state meets the preset conditions.

[0068] In some embodiments of the present application, before determining the number of damaged units in the first rock model, the method further includes:

[0069] determining whether the total stress is greater than a preset tensile strength value;

[0070] When the total stress is greater than a preset tensile strength value, determining that the damage cause of the first region is tensile strength;

[0071] If the total stress is less than or equal to the preset tensile strength value, determining whether the total stress is greater than the preset shear strength;

[0072] When the total stress is greater than the preset shear strength, it is determined that the damage cause of the first region is shear.

[0073] According to the semi-analytical method for identifying rock breaking by a roller cutter proposed in an embodiment of the present application, a cutter model and an initial rock model of a shield machine are obtained; the cutter model and the initial rock model are used to simulate the scene of the cutter cutting rock to obtain the contact area between the cutter model and the initial rock model; a first area is determined based on the contact area, and the rock model in the first area is discretized into cubic units; for each cubic unit, material properties are assigned to the cubic unit to obtain an updated first rock model, and the influence coefficient of each cubic unit is calculated based on the material properties; the load value of the cubic unit is updated based on the influence coefficient, and the total stress of the cubic unit is calculated based on the load value; the equivalent characteristic strain value of the first rock model and the convergence state of the equivalent characteristic strain value are determined based on the total stress of each cubic unit; when the convergence state meets the preset conditions, the number of damaged units in the first rock model is determined, and the degree of damage of the first rock model is determined based on the number, thereby improving the simulation accuracy of the rock model and further improving the accuracy of pipeline rock breaking identification.

[0074] Figure 2 This is a block diagram of a device for identifying rock breaking by a roller cutter based on a semi-analytical method according to an exemplary embodiment. Figure 2 The device includes an acquisition unit 201, a simulation unit 202, a discrete unit 203, an update unit 204, a calculation unit 205, a first determination unit 206 and a second determination unit 207.

[0075] The acquisition unit 201 is used to acquire the cutter model and initial rock model of the shield machine;

[0076] The simulation unit 202 is used to simulate the scene of the tool cutting the rock using the tool model and the initial rock model to obtain the contact area between the tool model and the initial rock model;

[0077] Discrete units 203 are used to determine a first region according to the contact region, and discretize the rock model in the first region into cubic units; the first region includes the contact region and the neighborhood of the contact region;

[0078] An updating unit 204 is configured to assign material attributes to each cubic unit to obtain an updated first rock model, and calculate an influence coefficient of each cubic unit based on the material attributes;

[0079] A calculation unit 205 is used to update the load value of the cubic unit according to the influence coefficient and calculate the total stress of the cubic unit according to the load value;

[0080] a first determining unit 206 for determining an equivalent characteristic strain value of the first rock model and a convergence state of the equivalent characteristic strain value according to the total stress of each cubic unit;

[0081] The second determining unit 207 is configured to determine the number of damaged units in the first rock model when the convergence state satisfies a preset condition, and determine the damage degree of the first rock model according to the number.

[0082] In some embodiments of the present application, the updating unit 204 may be specifically configured to:

[0083] Assigning different material properties to each cubic unit to obtain an updated first rock model;

[0084] The influence coefficients are calculated for each cubic element based on the material properties; the influence coefficients include the inhomogeneity coefficient.

[0085] In some embodiments of the present application, the calculation unit 205 may be specifically configured to:

[0086] The conjugate gradient method is used to calculate the first pressure distribution data of the contact area according to the load value;

[0087] The second pressure distribution data of the neighborhood of the contact area is calculated based on the influence coefficient using the Fourier transform method;

[0088] The total stress of the cubic element is calculated based on the first pressure distribution data and the second pressure distribution data.

[0089] In some embodiments of the present application, the apparatus may further include:

[0090] The updating unit is further used to update the displacement state information of the tool model when the convergence state does not meet the preset conditions, and return to execute the steps of determining the first area according to the contact area and discretizing the rock model in the first area into cubic units based on the displacement state information until the convergence state meets the preset conditions.

[0091] In some embodiments of the present application, the apparatus may further include:

[0092] The determination unit is further used to determine whether the total stress is greater than a preset tensile strength value;

[0093] The determining unit is further configured to determine that the damage cause of the first region is tension when the total stress is greater than a preset tensile strength value;

[0094] The determination unit is further configured to determine whether the total stress is greater than a preset shear strength when the total stress is less than or equal to a preset tensile strength value;

[0095] The determination unit is further configured to determine that the damage cause of the first region is shear when the total stress is greater than a preset shear strength.

[0096] Regarding the apparatus in the above embodiment, the specific manner in which each module performs operations has been described in detail in the embodiment of the method, and will not be elaborated here.

[0097] According to the semi-analytical method-based rock cutter identification device proposed in an embodiment of the present application, the cutter model and initial rock model of the shield machine are obtained; the cutter model and the initial rock model are used to simulate the scene of the cutter cutting the rock to obtain the contact area between the cutter model and the initial rock model; the first area is determined according to the contact area, and the rock model in the first area is discretized into cubic units; for each cubic unit, material properties are assigned to the cubic unit to obtain an updated first rock model, and the influence coefficient of each cubic unit is calculated according to the material properties; the load value of the cubic unit is updated according to the influence coefficient, and the total stress of the cubic unit is calculated according to the load value; the equivalent characteristic strain value of the first rock model and the convergence state of the equivalent characteristic strain value are determined according to the total stress of each cubic unit; when the convergence state meets the preset conditions, the number of damaged units in the first rock model is determined, and the damage degree of the first rock model is determined according to the number, thereby improving the simulation accuracy of the rock model and further improving the accuracy of pipeline rock breaking identification.

[0098] Figure 3 This is a block diagram of an apparatus for a semi-analytical method for identifying rock breakage using a roller cutter, according to an exemplary embodiment. For example, apparatus 300 may be an electronic device, such as a mobile phone, a computer, a digital broadcast terminal, a messaging device, a game console, a tablet device, a medical device, a fitness device, a personal digital assistant, or the like.

[0099] Reference Figure 3 , apparatus 300 may include one or more of the following components: a processing component 302 , a memory 304 , a power component 306 , a multimedia component 308 , an audio component 310 , an input / output (I / O) interface 312 , a sensor component 314 , and a communication component 316 .

[0100] The processing component 302 generally controls the overall operation of the device 300, such as operations associated with display, phone calls, data communications, camera operation, and recording operations. The processing component 302 may include one or more processors 320 to execute instructions to perform all or part of the steps of the above-described method. In addition, the processing component 302 may include one or more modules to facilitate interaction between the processing component 302 and other components. For example, the processing component 302 may include a multimedia module to facilitate interaction between the multimedia component 308 and the processing component 302.

[0101] The memory 304 is configured to store various types of data to support operations on the device 300. Examples of such data include instructions for any application or method operating on the device 300, contact data, phone book data, messages, pictures, videos, etc. The memory 304 can be implemented by any type of volatile or non-volatile storage device, or a combination thereof, such as static random access memory (SRAM), electrically erasable programmable read-only memory (EEPROM), erasable programmable read-only memory (EPROM), programmable read-only memory (PROM), read-only memory (ROM), magnetic memory, flash memory, magnetic disk, or optical disk.

[0102] The power component 306 provides power to the various components of the device 300. The power component 306 may include a power management system, one or more power supplies, and other components associated with generating, managing, and distributing power to the device 300.

[0103] The multimedia component 308 includes a screen that provides an output interface between the device 300 and the user. In some embodiments, the screen may include a liquid crystal display (LCD) and a touch panel (TP). If the screen includes a touch panel, the screen can be implemented as a touch screen to receive input signals from the user. The touch panel includes one or more touch sensors to sense touches, slides, and gestures on the touch panel. The touch sensor can not only sense the boundaries of the touch or slide action, but also detect the duration and pressure associated with the touch or slide operation. In some embodiments, the multimedia component 308 includes a front camera and / or a rear camera. When the device 300 is in an operating mode, such as a shooting mode or a video mode, the front camera and / or the rear camera can receive external multimedia data. Each front camera and rear camera can be a fixed optical lens system or have a focal length and optical zoom capability.

[0104] The audio component 310 is configured to output and / or input audio signals. For example, the audio component 310 includes a microphone (MIC) that is configured to receive external audio signals when the device 300 is in an operating mode, such as a call mode, a recording mode, and a voice recognition mode. The received audio signals may be further stored in the memory 304 or transmitted via the communication component 316. In some embodiments, the audio component 310 further includes a speaker for outputting audio signals.

[0105] I / O interface 312 provides an interface between processing component 302 and peripheral interface modules, such as a keyboard, click wheel, buttons, etc. These buttons may include but are not limited to: a home button, volume buttons, a start button, and a lock button.

[0106] The sensor assembly 314 includes one or more sensors for providing various aspects of the status assessment of the device 300. For example, the sensor assembly 314 can detect the open / closed state of the device 300, the relative positioning of components, such as the display and keypad of the device 300. The sensor assembly 314 can also detect changes in the position of the device 300 or a component of the device 300, the presence or absence of user contact with the device 300, the orientation or acceleration / deceleration of the device 300, and temperature changes of the device 300. The sensor assembly 314 may include a proximity sensor configured to detect the presence of nearby objects without any physical contact. The sensor assembly 314 may also include an optical sensor, such as a CMOS or CCD image sensor, for use in imaging applications. In some embodiments, the sensor assembly 314 may also include an accelerometer, a gyroscope sensor, a magnetic sensor, a pressure sensor, or a temperature sensor.

[0107] The communication component 316 is configured to facilitate wired or wireless communication between the device 300 and other devices. The device 300 can access a wireless network based on a communication standard, such as WiFi, 2G or 3G, or a combination thereof. In an exemplary embodiment, the communication component 316 receives a broadcast signal or broadcast-related information from an external broadcast management system via a broadcast channel. In an exemplary embodiment, the communication component 316 also includes a near field communication (NFC) module to facilitate short-range communication. For example, the NFC module can be implemented based on radio frequency identification (RFID) technology, infrared data association (IrDA) technology, ultra-wideband (UWB) technology, Bluetooth (BT) technology and other technologies.

[0108] In an exemplary embodiment, the apparatus 300 may be implemented by one or more application-specific integrated circuits (ASICs), digital signal processors (DSPs), digital signal processing devices (DSPDs), programmable logic devices (PLDs), field programmable gate arrays (FPGAs), controllers, microcontrollers, microprocessors, or other electronic components to perform the above-described method.

[0109] In an exemplary embodiment, a non-transitory computer-readable storage medium including instructions is also provided, such as a memory 304 including instructions, which can be executed by the processor 320 of the apparatus 300 to perform the above method. For example, the non-transitory computer-readable storage medium can be a ROM, a random access memory (RAM), a CD-ROM, a magnetic tape, a floppy disk, an optical data storage device, etc.

[0110] In an exemplary embodiment, a computer program product is also provided, comprising a computer program, which implements the above method when executed by the processor 320 of the apparatus 300 .

[0111] Other embodiments of the present invention will readily occur to those skilled in the art after considering the specification and practicing the invention disclosed herein. This application is intended to cover any variations, uses, or adaptations of the present invention that follow the general principles of the invention and include common knowledge or customary techniques in the art not disclosed herein. The description and examples are to be considered as exemplary only, with the true scope and spirit of the invention being indicated by the following claims.

[0112] It should be understood that the present invention is not limited to the exact construction described above and shown in the drawings, and that various modifications and changes may be made without departing from the scope thereof. The scope of the present invention is limited only by the appended claims.

Claims

1. A method for identifying rock breaking by a roller cutter based on a semi-analytical method, characterized in that: include: Obtain the shield machine's cutter model and initial rock model; The tool model and the initial rock model are used to simulate a scene of the tool cutting the rock, and a contact area between the tool model and the initial rock model is obtained; Determine a first region based on the contact region, and discretize the rock model in the first region into cubic units; the first region includes the contact region and the neighborhood of the contact region; For each cubic unit, assign material attributes to the cubic unit to obtain an updated first rock model, and calculate an influence coefficient of each cubic unit according to the material attributes; updating the load value of the cubic unit according to the influence coefficient, and calculating the total stress of the cubic unit according to the load value; determining an equivalent characteristic strain value of the first rock model and a convergence state of the equivalent characteristic strain value according to the total stress of each cubic unit; When the convergence state satisfies a preset condition, the number of damaged units in the first rock model is determined, and the damage degree of the first rock model is determined according to the number.

2. The semi-analytical method for identifying rock breaking by a roller cutter according to claim 1, characterized in that: Assigning material attributes to each cubic unit to obtain an updated first rock model, and calculating an influence coefficient of each cubic unit according to the material attributes, including: Assigning different material properties to each cubic unit to obtain an updated first rock model; An influence coefficient of each cubic element is calculated based on the material properties; the influence coefficient includes a non-uniformity coefficient.

3. The semi-analytical method for identifying rock breaking by a roller cutter according to claim 1, characterized in that: The updating of the load value of the cubic unit according to the influence coefficient and the calculation of the total stress of the cubic unit according to the load value include: Calculating first pressure distribution data of the contact area according to the load value using a conjugate gradient method; Using a Fourier transform method, calculating second pressure distribution data of a neighborhood of the contact area according to the influence coefficient; The total stress of the cubic unit is calculated according to the first pressure distribution data and the second pressure distribution data.

4. The semi-analytical method for identifying rock breaking by a roller cutter according to claim 1, characterized in that: Also includes: When the convergence state does not meet the preset conditions, the displacement state information of the tool model is updated, and based on the displacement state information, the step of determining the first area according to the contact area and discretizing the rock model in the first area into cubic units is returned to be executed until the convergence state meets the preset conditions.

5. The method for identifying rock breaking by a roller cutter based on a semi-analytical method according to claim 1, characterized in that: Before determining the number of damaged units in the first rock model, the method further includes: determining whether the total stress is greater than a preset tensile strength value; When the total stress is greater than a preset tensile strength value, determining that the damage cause of the first region is tensile strength; When the total stress is less than or equal to a preset tensile strength value, determining whether the total stress is greater than a preset shear strength; When the total stress is greater than the preset shear strength, it is determined that the damage cause of the first region is shear.

6. A device for identifying rock breaking by a roller cutter based on a semi-analytical method, characterized in that: include: An acquisition unit, used to acquire a cutter model and an initial rock model of a shield machine; a simulation unit, configured to simulate a scene of a tool cutting a rock using the tool model and the initial rock model, and obtain a contact area between the tool model and the initial rock model; Discrete units are used to determine a first area according to the contact area, and discretize the rock model in the first area into cubic units; the first area includes the contact area and the neighborhood of the contact area; an updating unit, configured to assign material properties to each cubic unit to obtain an updated first rock model, and calculate an influence coefficient of each cubic unit according to the material properties; a calculation unit, configured to update the load value of the cubic unit according to the influence coefficient, and calculate the total stress of the cubic unit according to the load value; a first determining unit, configured to determine an equivalent characteristic strain value of the first rock model and a convergence state of the equivalent characteristic strain value according to the total stress of each cubic unit; The second determining unit is configured to determine the number of damaged units in the first rock model when the convergence state satisfies a preset condition, and determine the damage degree of the first rock model according to the number.

7. The semi-analytical method-based rock cutter identification device according to claim 1, characterized in that: The updating unit is specifically configured to: Assigning different material properties to each cubic unit to obtain an updated first rock model; An influence coefficient of each cubic element is calculated based on the material properties; the influence coefficient includes a non-uniformity coefficient.

8. An electronic device, characterized in that: include: A memory, a processor, and a computer program stored in the memory and executable on the processor, wherein when the processor executes the computer program, the method according to any one of claims 1 to 5 is implemented.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the method according to any one of claims 1 to 5 is implemented.

10. A computer program product comprising a computer program, characterized in that The computer program implements the method according to any one of claims 1 to 5 when executed by a processor.

Images