Heterogeneous peridynamic simulation method, device and storage medium for rock material fracture failure

Abstract

The present invention discloses a method, device and storage medium for heterogeneous near-field dynamic simulation of rock material fracture and destruction, which includes establishing a numerical model of the rock material, determining each material point as a void material point or a solid material point, determining whether there is a void bond or a non-void bond between any two solid material points, establishing a motion equation between a first solid material point and all second solid material points, solving the motion equation, and determining the damage value of the first solid material point. The present invention can simulate the heterogeneous structural characteristics of rock material fracture and destruction, and can compare the heterogeneous structural characteristics with the macroscopic mechanical behavior of the rock material, thereby revealing the control effect of the heterogeneous characteristics of the rock on its macroscopic mechanical behavior, thereby obtaining a deeper and more systematic understanding of different lithologies, which is conducive to better controlling the mechanical behavior of rock materials and better serving engineering construction. The present invention is widely used in the field of computer technology.

Description

Technical Field

[0001] The present invention relates to the field of computer technology, and in particular to a method, device and storage medium for simulating the fracture and destruction of rock materials using heterogeneous peridynamics. Background Art

[0002] Numerous experiments and engineering practices have demonstrated that the physical and mechanical properties of rock are fundamental and crucial factors controlling the stability of projects such as slopes and underground caverns. Consequently, in diverse projects such as highways, tunnels, and mines, varying degrees of physical and mechanical property testing and monitoring have been conducted, tailored to the specific rock types and in line with engineering practice. In-depth analysis of their mechanical behavior has been conducted to provide realistic performance indicators for engineering design and a more rational basis for construction control. However, this research has primarily focused on repeated experiments, verification, and analysis of rock's macroscopic mechanical behavior, often treating rock as a homogeneous solid. The study and application of its heterogeneous properties, as well as their correlation with macroscopic properties, lags far behind. In fact, rock heterogeneity often strongly controls its macroscopic mechanical behavior. However, due to a lack of understanding of rock's heterogeneous structural characteristics, our understanding of rock's mechanical behavior remains incomplete and incomplete. Summary of the Invention

[0003] In order to address the current technical issues such as the lack of understanding of the heterogeneous structural characteristics of rocks and therefore insufficient in-depth and systematic understanding of the mechanical behavior of rocks, the present invention aims to provide a method, device and storage medium for heterogeneous peridynamic simulation of rock material fracture failure.

[0004] In one aspect, an embodiment of the present invention includes a method for heterogeneous peridynamic simulation of rock material fracture failure, the method comprising the following steps:

[0005] Establishing a numerical model of rock material; the numerical model includes information of each material point in the rock material;

[0006] Obtain mineral composition information and porosity parameter information of rock materials;

[0007] Assigning material attributes to each of the material points according to the mineral composition information and the porosity parameter information, thereby determining each of the material points as a void material point or a solid material point;

[0008] Obtaining a key attribute between any two of the physical material points, thereby determining whether there is an empty key or a non-empty key between any two of the physical material points;

[0009] For a first physical material point, a motion equation is established between the first physical material point and all second physical material points, the motion equation is solved, and a damage value of the first physical material point is determined; wherein the first physical material point is any one of the physical material points, and the second physical material point is the physical material point located in the near field of the first physical material point and having a non-empty bond with the first physical material point.

[0010] Furthermore, the obtaining of mineral composition information and porosity parameter information of rock materials includes:

[0011] Performing an X-ray diffraction test on the rock material to obtain information on the mineral composition;

[0012] A scanning electron microscope test is performed on the rock material to obtain the porosity parameter information.

[0013] Furthermore, the X-ray diffraction test is performed on the rock material to obtain the mineral composition information, including:

[0014] Performing X-ray diffraction on the rock material to obtain a mineral crystal diffraction pattern of the rock material;

[0015] Performing a phase search on the mineral crystal diffraction pattern to obtain peak characteristics;

[0016] The mineral composition information is determined based on the peak characteristics.

[0017] Furthermore, the rock material is subjected to a scanning electron microscope test to obtain the porosity parameter information, including:

[0018] Performing a scanning electron microscope test on the rock material to obtain a two-dimensional microscopic image of the rock material;

[0019] Using the grayscale value in the two-dimensional microscopic image as a vertical elevation, stretching the image along the vertical elevation to obtain a three-dimensional scanning electron microscope model;

[0020] The total volume and void volume of the three-dimensional scanning electron microscope model are obtained, and the three-dimensional porosity is calculated according to the total volume and the void volume as the porosity parameter information.

[0021] Furthermore, assigning material properties to each of the material points according to the mineral composition information and the porosity parameter information includes:

[0022] Generate a material parameter array according to the mineral composition information and the porosity parameter information;

[0023] Using a shuffling algorithm to shuffle the material parameter array;

[0024] Allocating each of the mineral component composition information and the porosity parameter information in the disrupted material parameter array to a corresponding material point;

[0025] For any of the material points, when the mineral composition information assigned to the material point is greater than a first threshold and the porosity parameter information is less than a second threshold, the material point is determined as a solid material point, otherwise it is determined as a void material point.

[0026] Furthermore, the obtaining of the key attribute between any two of the physical material points, thereby determining whether there is an empty key or a non-empty key between any two of the physical material points, includes:

[0027] For any two of the physical material points, establishing a connection line between the two physical material points;

[0028] When the line connecting two of the entity material points passes through a cell where a void material point is located, it is determined that there is a void bond between the two entity material points; otherwise, it is determined that there is a non-void bond between the two entity material points.

[0029] Furthermore, establishing a motion equation between the first physical material point and all the second physical material points, solving the motion equation, and determining the damage value of the first physical material point includes:

[0030] Establish the following motion equation

[0031]

[0032]

[0033]

[0034]

[0035] Solving the motion equation using an explicit integration algorithm to determine a damage value of the first solid material point;

[0036] in, x represents the first entity material point, represents the second entity material point, Represents the first entity material point x The density, Represents the first entity material point x exist The displacement of time, express The second derivative with respect to time, Represents the second entity material point exist The displacement of time, Represents the first entity material point x Near field, satisfying , is the neighborhood radius, Represents the second entity material point The volume, Represents the first entity material point x exist The force density of the external force or body force at that moment, s represents the bond elongation; represents the bond breaking parameter;

[0037] f Indicates that when the first entity material point x With the second entity material point Bond elongation between s Less than critical elongation When the first physical material point x With the second entity material point The interaction force between Represents the first entity material point x With the second entity material point The relative position in the reference configuration, Represents the first entity material point x With the second entity material point Relative displacement in the reference configuration, critical elongation satisfy ,in represents the energy release rate; c Represents the bond microelastic modulus, satisfying ,in is the elastic modulus of the rock material, is the thickness of the numerical model, is Poisson's ratio.

[0038] Furthermore, solving the motion equation using an explicit integration algorithm to determine the damage value of the first solid material point includes:

[0039] Calculate according to the following formula

[0040] ;

[0041] Obtaining the damage value of the first physical material point; wherein, Represents the first entity material point x exist The damage value at the moment, represents a set consisting of all the second entity material points corresponding to the first entity material points.

[0042] On the other hand, an embodiment of the present invention also includes a computer device including a memory and a processor, the memory being used to store at least one program, and the processor being used to load at least one program to execute a heterogeneous near-field dynamic simulation method for rock material fracture failure in an embodiment.

[0043] On the other hand, an embodiment of the present invention further includes a storage medium storing a program executable by a processor. When the program is executed by the processor, it is used to execute a heterogeneous near-field dynamic simulation method for rock material fracture and destruction in an embodiment.

[0044] The beneficial effects of the present invention are as follows: the heterogeneous peridynamic simulation method for rock material fracture failure in the embodiment can simulate the heterogeneous structural characteristics of rock material fracture failure, and can compare the heterogeneous structural characteristics with the macroscopic mechanical behavior of the rock material, thereby revealing the control effect of the heterogeneous characteristics of the rock on its macroscopic mechanical behavior, thereby obtaining a deeper and more systematic understanding of different rock types, which is conducive to better control of the mechanical behavior of rock materials and better service for engineering construction. BRIEF DESCRIPTION OF THE DRAWINGS

[0045] Figure 1 Schematic diagram of the process of the heterogeneous peridynamic simulation method for rock material fracture failure in the embodiment;

[0046] Figure 2 A schematic diagram of the step of obtaining the bond properties between any two physical material points in the embodiment;

[0047] Figure 3 This is a comparison diagram of the actual failure forms of two types of rocks in the embodiment and the simulation results;

[0048] Figure 4 Comparison of the two rock test and numerical load-displacement curves in the example

[0049] Figure 5 Schematic diagram of the simulation results of the destruction process of rocks with different porosities in the embodiment. DETAILED DESCRIPTION

[0050] Rock heterogeneity controls its macroscopic mechanical properties and fracture behavior. Existing research on rock mechanical properties focuses on treating rock materials as homogeneous or defining material properties based on Weibull random distribution, failing to consider the mineral composition and porosity characteristics of rocks.

[0051] Based on the above principles, this embodiment provides a heterogeneous peridynamic simulation method for rock material fracture failure. Figure 1 The heterogeneous peridynamic simulation method for rock material fracture failure includes the following steps:

[0052] S1. Establish numerical models of rock materials;

[0053] S2. Obtaining mineral composition information and porosity parameter information of rock materials;

[0054] S3. Assigning material properties to each material point based on the mineral composition information and the porosity parameter information, thereby determining each material point as a void material point or a solid material point;

[0055] S4. Obtain the key attribute between any two physical material points, thereby determining whether there is an empty key or a non-empty key between any two physical material points;

[0056] S5. For the first physical material point, establish a motion equation between the first physical material point and all the second physical material points, solve the motion equation, and determine the damage value of the first physical material point.

[0057] The process of steps S1-S5 is as follows Figure 1 Steps S1-S5 may be executed by a computer.

[0058] In step S1, a rock sample to be simulated is obtained, and a numerical model is established based on the sample's geometric dimensions and boundary conditions. The numerical model includes multiple spatial regions, into which each material point in the rock material is evenly discretized. The numerical model then stores information about each material point in the rock material. This information includes information about mineral composition and porosity parameters.

[0059] In step S2, X-ray diffraction and scanning electron microscopy tests are performed on the rock material to obtain mineral composition information and porosity parameter information of the rock material. Specifically, when executing step S2, the following steps may be performed:

[0060] S201. Perform X-ray diffraction (XRD) on the rock material to obtain a mineral crystal diffraction pattern of the rock material;

[0061] S202. Performing a phase search on the mineral crystal diffraction pattern to obtain peak characteristics;

[0062] S203. Determine the mineral composition information based on the peak characteristics;

[0063] S204. Perform a scanning electron microscope (SEM) test on the rock material to obtain a two-dimensional microscopic image of the rock material;

[0064] S205. Using the grayscale value in the two-dimensional microscopic image as the vertical elevation, stretching the image along the vertical elevation to obtain a three-dimensional scanning electron microscope model;

[0065] S206. Obtain the total volume and void volume of the three-dimensional scanning electron microscope model, and calculate the three-dimensional porosity based on the total volume and void volume as porosity parameter information.

[0066] Steps S201-S203 are steps for obtaining mineral composition information of the rock material. During step S202, Jade software can be used to perform a phase search on the mineral crystal diffraction pattern to obtain peak characteristics. The location of the peak indicates the type of mineral component in the rock material, and the size of the peak indicates the proportion of the mineral component in the rock material, thereby obtaining mineral composition information.

[0067] Steps S204-S206 are steps for obtaining porosity parameter information of the rock material. After executing step S204, the two-dimensional microscopic image obtained in step S204 can be preprocessed using the cv2 module of the OpenCV library. This preprocessing includes brightness adjustment, noise reduction, filtering, Gauss transform, and other processes to make the two-dimensional microscopic image clearer. The preprocessed two-dimensional microscopic image is used to execute step S205.

[0068] In step S205, for each pixel in the 2D microscopic image, the plane coordinates of the pixel in the 2D microscopic image remain unchanged, and the vertical elevation of each pixel is set equal to its grayscale value, so that each pixel has a vertical coordinate. The vertical elevation of each pixel is then multiplied by the same coefficient, so that the 2D microscopic image is stretched along the vertical elevation to obtain a 3D scanning electron microscope model.

[0069] In step S206, the total volume and void volume of each unit area in the 3D SEM model obtained in step S205 are calculated. Based on the total volume and void volume of each unit area, the 3D porosity of this unit area is calculated, and the 3D porosity is used as the porosity parameter information. Specifically, the quotient of the void volume and the total volume of the same unit area can be calculated to obtain the 3D porosity of this unit area. Each unit area in the 3D SEM model corresponds to a material point of the rock material, so the 3D porosity of each material point of the rock material can be obtained.

[0070] In this embodiment, when executing step S3, that is, assigning material properties to each material point based on the mineral composition information and the porosity parameter information, the following steps may be specifically performed:

[0071] S301. Generate a material parameter array based on the mineral composition information and porosity parameter information;

[0072] S302. Using a shuffling algorithm to shuffle the material parameter array;

[0073] S303. Assign each mineral component composition information and porosity parameter information in the disrupted material parameter array to a corresponding material point;

[0074] S304. For any material point, when the mineral composition information assigned to the material point is greater than a first threshold and the porosity parameter information is less than a second threshold, the material point is determined as a solid material point; otherwise, it is determined as a void material point.

[0075] By executing steps S201-S206, each material point in the rock material has corresponding mineral composition information and porosity parameter information. In step S301, the corresponding (mineral composition information and porosity parameter information) for each material point can be sorted according to the order of the material points in the numerical model to form a material parameter array. For example, if the order of the material points in the numerical model is material point 1, material point 2, material point 3, etc., then the order in the material parameter array will be (mineral composition information 1, porosity parameter information 1), (mineral composition information 2, porosity parameter information 2), (mineral composition information 3, porosity parameter information 3), etc.

[0076] In step S302, a shuffling algorithm is used to shuffle the order of the arrays in the form of (mineral composition information, porosity parameter information) in the material parameter array. The order of the shuffled material parameter array might be (mineral composition information 3, porosity parameter information 3), (mineral composition information 100, porosity parameter information 100), (mineral composition information 5, porosity parameter information 5), and so on.

[0077] In step S303, the shuffled array of (mineral composition information, porosity parameter information) is sorted according to the order of the material points in the numerical model and assigned to each material point. For example, based on the example in step S302, material point 1 is assigned to (mineral composition information 3, porosity parameter information 3), material point 2 is assigned to (mineral composition information 100, porosity parameter information 100), material point 3 is assigned to (mineral composition information 5, porosity parameter information 5), and so on.

[0078] In step S304, for any material point, the following steps are used to determine whether it is a solid material point or a void material point: if the mineral composition information assigned to this material point is greater than a first threshold (indicating that it contains a large amount of a certain specific mineral component) and the porosity parameter information is less than a second threshold (indicating that the porosity is small), then this material point is determined to be a solid material point; otherwise, this material point is determined to be a void material point.

[0079] Step S304 is performed on all material points, so that the material properties of each material point can be identified, that is, whether a material point is a solid material point or a void material point.

[0080] By executing steps S301 to S304, material properties can be randomly assigned to each material point in the numerical model of the rock material, which is conducive to simulating the heterogeneous peridynamic properties of the rock material fracture failure.

[0081] In this embodiment, when executing step S4, that is, obtaining the key attribute between any two physical material points, thereby determining whether there is an empty key or a non-empty key between any two physical material points, the following steps may be specifically performed:

[0082] S401. For any two physical material points, establish a line between the two physical material points;

[0083] S402. When the line connecting two solid material points passes through a cell where a void material point is located, it is determined that there is a void bond between the two solid material points. Otherwise, it is determined that there is a non-void bond between the two solid material points.

[0084] The principle of steps S401-S402 is as follows Figure 2 As shown. Figure 2 There is interaction between the discrete points in two adjacent solid units (i.e., solid material points), but there is no interaction between the discrete points in a solid unit (i.e., solid material points) and the discrete points in a void unit (i.e., void material points), and the void unit has a blocking effect on the transmission of the interaction (bond). Therefore, whether there is interaction between the discrete points can be determined by judging whether the discrete points are in the void unit (i.e., whether they are void material points) or whether the line between the discrete points passes through the void unit.

[0085] In step S401, for any two physical material points, a line is established between the two physical material points. The line can be virtual and calculated by the coordinates of the two physical material points in the numerical model. In step S402, the position passed by the line calculated in step S401 is determined. If the line passes through a cell where a void material point is located (for example, the distance between the line and a void material point is less than a distance threshold), then it is determined that there is no interaction between the two physical material points in step S401, that is, there is an empty bond between the two physical material points. Otherwise, it is determined that there is an interaction between the two physical material points, that is, there is a non-empty bond between the two physical material points (that is, there is a gap). Figure 2 in the .

[0086] In this embodiment, the goal of executing step S5 is to calculate the displacement and damage value of any physical material point at any time within a certain period of time. Figure 1 , we can take the i-th physical material point among all the physical material points and record it as the first physical material point x(i) , and set the current time to By executing step S5 once, the first entity material point is calculated x(i) exist Displacement of time and damage value , after traversing all i and all , the displacement and damage value of any solid material point at any time can be obtained.

[0087] Reference Figure 1 , for the i-th physical material point, that is, the first physical material point x(i) , in its near field Select other physical points and traverse the near field All the physical material points within will be combined with the first physical material point x(i) The physical material points with non-empty bonds between them are screened out and recorded as the second physical material points .

[0088] Reference Figure 1 , for the first solid material point x(i) Its near field All second entity material points within , establish the motion equations (1)-(4) between them:

[0089] (1)

[0090] (2)

[0091] (3)

[0092] (4)

[0093] In the equations of motion (1)-(4), the default first entity material point x is the i-th physical point, i.e. x(i) , so the symbol i can be ignored. The meaning of the symbols in the equation of motion is: s represents the bond elongation, x represents the first solid material point, represents the second entity material point, Represents the first solid material point x The density, Represents the first solid material point x exist The displacement of time, express The second derivative with respect to time, Represents the second entity material point exist The displacement of time, Represents the first solid material point x Near field, satisfying , is the neighborhood radius, Represents the second entity material point The volume, Represents the first solid material point x exist The force density of the external force or body force at that moment, Represents the first solid material point x With the second entity material point The relative position in the reference configuration, Represents the first solid material point x With the second entity material point The relative displacement in the reference configuration, c Represents the bond microelastic modulus, satisfying ,in is the elastic modulus of rock material, is the thickness of the numerical model, is Poisson's ratio.

[0094] In the equation of motion (2), f Represents the first solid material point x With the second entity material point The interaction force between represents the bond breaking parameter, . Reference Figure 2 , for any second entity material point (jth solid material point (j) ), the second entity material point can be calculated by the motion equation (4) (j) With the first solid material point x(i) Bond elongation between s , the bond elongation s and critical elongation For comparison, if s ≤ , then the second entity material point can be calculated according to the motion equation (2) (j) With the first solid material point x(i) The interaction force between them is the bond force, otherwise the second entity material point will be interrupted (j) With the first solid material point x(i) The bond between the two entities, such as the second entity point (j) Corresponding f Set to 0. Traverse all j to get any second entity material point Corresponding to the equation of motion (2) f The value of .

[0095] In the equation of motion (2), the critical elongation used is It can be determined by equation (5):

[0096] (5)

[0097] In equation (5), Represents the energy release rate.

[0098] In determining the first solid material point x(i) and all corresponding second entity material points After that, the values ​​of the parameters in the motion equations (1)-(4) can be determined, and the motion equations (1)-(4) can be solved by explicit integration algorithm to obtain the first solid material point x(i) exist Displacement of time , according to the solution results of the motion equations (1)-(4), solve equation (6):

[0099] (6)

[0100] Thus, the first entity material point is obtained x(i) exist Damage value at the moment .

[0101] In equation (6), Represented by the first entity material point x(i) All corresponding second entity material points A collection of components.

[0102] In this embodiment, refer to Figure 1 After executing a round of step S5, the first entity material point can be obtained x(i) exist Displacement of time and damage value .Keep The first physical point is obtained by setting i=i+1 and executing step S5 again. x(i+1) exist The displacement and damage value at the moment, until any first solid material point is obtained Displacement and damage values ​​at the moment.

[0103] In this embodiment, refer to Figure 1 , after traversing all the first entity material points After the displacement and damage values ​​at the moment, set = +1, and set i back to the initial value (for example, 1), and execute step S5 again to obtain the value of any first entity material point in +1 moment, until all moments within a specific time period are traversed, and finally the displacement and damage values ​​of any physical material point at any moment can be obtained, completing the heterogeneous peridynamic simulation of rock material fracture and destruction.

[0104] The heterogeneous peridynamic simulation method for rock material fracture failure in this embodiment can simulate the heterogeneous structural characteristics of rock material fracture failure, and can compare the heterogeneous structural characteristics with the macroscopic mechanical behavior of the rock material, thereby revealing the control effect of the heterogeneous characteristics of the rock on its macroscopic mechanical behavior, thereby obtaining a deeper and more systematic understanding of different rock types, which is conducive to better control of the mechanical behavior of rock materials and better serving engineering construction.

[0105] A computer program for executing the heterogeneous peridynamic simulation method for rock material fracture failure in this embodiment can be written and written into a computer device or storage medium. When the computer program is read out and run, the heterogeneous peridynamic simulation method for rock material fracture failure in this embodiment is executed, thereby achieving the same technical effect as the heterogeneous peridynamic simulation method for rock material fracture failure in the embodiment.

[0106] Implementation 1:

[0107] The heterogeneous peridynamic simulation method for rock material fracture failure described in the embodiment can be verified through a rock tensile strength test. Fine-grained sandstone and sandy mudstone were used as rock specimens for the tensile strength test. Rock specimens with an aspect ratio of 1:2 were used, and cylindrical specimens measuring 25 mm x 50 mm were used. The test procedure involved placing the specimen between the upper and lower bearing plates of a material testing machine. The machine was then operated and loaded at a rate of 0.1 mm / min until the specimen failed.

[0108] At the same time, XRD and SEM tests were carried out on the two rocks respectively to obtain the rock mineral composition and proportion as well as porosity parameters.

[0109] According to the mechanical properties test results of the two rocks, a disk model with a diameter of 50 mm was established, and the material point spacing was , the total number of material points is 30965, the near field range , time step Δ t = 0.01s. For the sandy mudstone model, the porosity was set to 20%, and the mineral composition included quartz (74.3%), nickel chlorite (13.7%), low-sodium feldspar (1.8%), kaolinite (7.8%), and mica (2.4%). For the fine-grained sandstone model, the porosity was set to 25%, and the mineral composition included quartz (81.7%), nickel chlorite (6.5%), low-sodium feldspar (5.8%), kaolinite (5.2%), and mica (0.8%). A 5×10 -4 mm / s displacement boundary condition. Comparison of the actual failure modes of two rocks with the simulation results. Figure 3 The comparison of the two rock test and numerical load-displacement curves is shown in Figure 4 shown.

[0110] according to Figure 3 and Figure 4 The verification results shown show that the heterogeneous peridynamic simulation method for rock material fracture failure in the embodiment has a good simulation effect close to the actual failure morphology of rock.

[0111] Implementation 2:

[0112] The heterogeneous peridynamic simulation method for rock material fracture failure in the embodiment can be used to analyze the effect of porosity on the rock failure process.

[0113] Create a rock model with mineral composition including quartz (75%), nickel chlorite (14%), low sodium feldspar (2%), kaolinite (7%), and mica (2%). , the total number of material points is 30965, the near field range , time step Δ t = 0.01s. Apply 5×10-4 mm / s displacement boundary condition. For structural models with different porosity (0, 5%, 15%, 20%, 30%), the calculation results of 5%, 30%, 50%, 70%, and 100% of the peak load were taken for observation. The fracture development process of the model is as follows: Figure 5 shown.

[0114] according to Figure 5 It can be seen from the analysis results shown that the heterogeneous peridynamic simulation method for rock material fracture failure in the embodiment has a good simulation effect on the rock failure process.

[0115] It should be noted that, unless otherwise specified, when a feature is referred to as being "fixed" or "connected" to another feature, it may be directly fixed or connected to the other feature, or it may be indirectly fixed or connected to the other feature. In addition, the descriptions of up, down, left, right, etc. used in this disclosure are only relative to the relative positional relationships of the components of the present disclosure in the accompanying drawings. The singular forms of "a", "" and "the" used in this disclosure are also intended to include the plural forms, unless the context clearly indicates otherwise. In addition, unless otherwise defined, all technical and scientific terms used in this embodiment have the same meaning as those generally understood by those skilled in the art. The terms used in the specification of this embodiment are only for describing specific embodiments and are not intended to limit the invention. The term "and / or" used in this embodiment includes any combination of one or more related listed items.

[0116] It should be understood that, although the terms first, second, third, etc. may be used to describe various elements in the present disclosure, these elements should not be limited to these terms. These terms are only used to distinguish elements of the same type from each other. For example, without departing from the scope of the present disclosure, the first element may also be referred to as the second element, and similarly, the second element may also be referred to as the first element. The use of any and all examples or exemplary language ("for example", "such as", etc.) provided in the present embodiment is only intended to better illustrate embodiments of the present invention, and unless otherwise required, will not impose limitations on the scope of the present invention.

[0117] It should be appreciated that embodiments of the present invention can be implemented or practiced by computer hardware, a combination of hardware and software, or by computer instructions stored in a non-transitory computer-readable memory. The methods can be implemented in a computer program using standard programming techniques - including a non-transitory computer-readable storage medium configured with a computer program, wherein the storage medium so configured causes the computer to operate in a specific and predefined manner - according to the methods and figures described in the specific embodiments. Each program can be implemented in a high-level procedural or object-oriented programming language to communicate with the computer system. However, if desired, the program can be implemented in assembly or machine language. In any case, the language can be a compiled or interpreted language. In addition, the program can be run on a programmed application-specific integrated circuit for this purpose.

[0118] Furthermore, the operations of the processes described in this embodiment may be performed in any suitable order, unless otherwise indicated in this embodiment or otherwise clearly contradicted by the context. The processes described in this embodiment (or variations and / or combinations thereof) may be executed under the control of one or more computer systems configured with executable instructions, and may be implemented as code (e.g., executable instructions, one or more computer programs, or one or more applications) that is executed collectively on one or more processors, by hardware, or a combination thereof. A computer program includes multiple instructions that can be executed by one or more processors.

[0119] Furthermore, the method can be implemented in any type of computing platform that is operably connected to a suitable computer, including but not limited to a personal computer, a minicomputer, a mainframe, a workstation, a network or distributed computing environment, a separate or integrated computer platform, or in communication with a charged particle tool or other imaging device, etc. Various aspects of the present invention can be implemented as machine-readable code stored on a non-transitory storage medium or device, whether removable or integrated into a computing platform, such as a hard disk, an optical read and / or write storage medium, RAM, ROM, etc., so that it can be read by a programmable computer, and when the storage medium or device is read by the computer, it can be used to configure and operate the computer to perform the process described herein. In addition, the machine-readable code, or portions thereof, can be transmitted over a wired or wireless network. When such media includes instructions or programs that implement the above steps in conjunction with a microprocessor or other data processor, the invention of this embodiment includes these and other different types of non-transitory computer-readable storage media. When programmed according to the methods and techniques of the present invention, the present invention also includes the computer itself.

[0120] The computer program can be applied to input data to perform the functions of the present embodiment, thereby converting the input data to generate output data that is stored in a non-volatile memory. The output information can also be applied to one or more output devices such as a display. In a preferred embodiment of the present invention, the converted data represents a physical and tangible object, including a specific visual depiction of the physical and tangible object produced on the display.

[0121] The above are merely preferred embodiments of the present invention. The present invention is not limited to the aforementioned embodiments. As long as the technical effects of the present invention are achieved by the same means, any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention shall be included within the scope of protection of the present invention. Within the scope of protection of the present invention, various modifications and variations of the technical solutions and / or implementation methods may be made.

Claims

1. A heterogeneous peridynamic simulation method for rock material fracture failure, characterized by: The heterogeneous peridynamic simulation method for rock material fracture failure includes: Establishing a numerical model of rock material; the numerical model includes information of each material point in the rock material; Obtain mineral composition information and porosity parameter information of rock materials; Assigning material attributes to each of the material points according to the mineral composition information and the porosity parameter information, thereby determining each of the material points as a void material point or a solid material point; Obtaining a key attribute between any two of the physical material points, thereby determining whether there is an empty key or a non-empty key between any two of the physical material points; For a first physical material point, establishing a motion equation between the first physical material point and all second physical material points, solving the motion equation, and determining a damage value of the first physical material point; wherein the first physical material point is any one of the physical material points, and the second physical material points are the physical material points located in the near field of the first physical material point and having a non-empty bond with the first physical material point; Assigning material properties to each of the material points according to the mineral composition information and the porosity parameter information includes: Generate a material parameter array according to the mineral composition information and the porosity parameter information; Using a shuffling algorithm to shuffle the material parameter array; Allocating each of the mineral component composition information and the porosity parameter information in the disrupted material parameter array to a corresponding material point; For any of the material points, when the mineral composition information assigned to the material point is greater than a first threshold and the porosity parameter information is less than a second threshold, the material point is determined as a solid material point, otherwise it is determined as a void material point.

2. The heterogeneous peridynamic simulation method for rock material fracture failure according to claim 1 is characterized in that: The obtaining of mineral composition information and porosity parameter information of rock materials includes: Performing an X-ray diffraction test on the rock material to obtain information on the mineral composition; A scanning electron microscope test is performed on the rock material to obtain the porosity parameter information.

3. The heterogeneous peridynamic simulation method for rock material fracture failure according to claim 2 is characterized in that: The X-ray diffraction test is performed on the rock material to obtain the mineral composition information, including: Performing X-ray diffraction on the rock material to obtain a mineral crystal diffraction pattern of the rock material; Performing a phase search on the mineral crystal diffraction pattern to obtain peak characteristics; The mineral composition information is determined based on the peak characteristics.

4. The heterogeneous peridynamic simulation method for rock material fracture failure according to claim 2, characterized in that: The scanning electron microscope test is performed on the rock material to obtain the porosity parameter information, including: Performing a scanning electron microscope test on the rock material to obtain a two-dimensional microscopic image of the rock material; Using the grayscale value in the two-dimensional microscopic image as a vertical elevation, stretching the image along the vertical elevation to obtain a three-dimensional scanning electron microscope model; The total volume and void volume of the three-dimensional scanning electron microscope model are obtained, and the three-dimensional porosity is calculated according to the total volume and the void volume as the porosity parameter information.

5. The heterogeneous peridynamic simulation method for rock material fracture failure according to claim 1 is characterized in that: The acquiring of the key attribute between any two of the physical material points, thereby determining whether there is an empty key or a non-empty key between any two of the physical material points, includes: For any two of the physical material points, establishing a connection line between the two physical material points; When the line connecting two of the entity material points passes through a cell where a void material point is located, it is determined that there is a void bond between the two entity material points; otherwise, it is determined that there is a non-void bond between the two entity material points.

6. The heterogeneous peridynamic simulation method for rock material fracture failure according to any one of claims 1 to 5, characterized in that: The step of establishing a motion equation between the first physical material point and all the second physical material points, solving the motion equation, and determining the damage value of the first physical material point includes: Establish the following motion equation Solving the motion equation using an explicit integration algorithm to determine a damage value of the first solid material point; in, x represents the first entity material point, represents the second entity material point, Represents the first entity material point x The density, Represents the first entity material point x exist The displacement of time, express The second derivative with respect to time, Represents the second entity material point exist The displacement of time, Represents the first entity material point x Near field, satisfying , is the neighborhood radius, Represents the second entity material point The volume, Represents the first entity material point x exist The force density of the external force or body force at that moment, s represents the bond elongation; represents the bond breaking parameter; f Indicates that when the first entity material point x With the second entity material point Bond elongation between s Less than critical elongation When the first physical material point x With the second entity material point The interaction force between Represents the first entity material point x With the second entity material point The relative position in the reference configuration, Represents the first entity material point x With the second entity material point Relative displacement in the reference configuration, critical elongation satisfy ,in represents the energy release rate; c Represents the bond microelastic modulus, satisfying ,in is the elastic modulus of the rock material, is the thickness of the numerical model, is Poisson's ratio.

7. The heterogeneous peridynamic simulation method for rock material fracture failure according to claim 6, characterized in that: Solving the motion equation using an explicit integration algorithm to determine the damage value of the first solid material point includes: Calculate according to the following formula ； Obtaining the damage value of the first physical material point; wherein, Represents the first entity material point x exist The damage value at the moment, represents a set consisting of all the second entity material points corresponding to the first entity material points.

8. A computer device, characterized in that: The method comprises a memory and a processor, wherein the memory is used to store at least one program, and the processor is used to load at least one program to execute the heterogeneous peridynamic simulation method for rock material fracture and failure according to any one of claims 1 to 7.

9. A computer-readable storage medium storing a program executable by a processor, characterized in that: The program executable by the processor is used to execute the heterogeneous peridynamic simulation method for rock material fracture and failure as described in any one of claims 1 to 7 when executed by the processor.

Images