A method and apparatus for blasting simulation of a phaneritic intrusive rock

Abstract

The application provides a phanerocrystalline intrusive rock blasting simulation method and device, comprising the following steps: constructing an initial rock crystal sample uniaxial compression model according to the mesoscopic physical parameters of a target phanerocrystalline intrusive rock; generating a cohesive force unit on the contact interface between the grains in the uniaxial compression model to obtain a target rock crystal sample uniaxial compression model; performing a uniaxial compression simulation test on the target rock crystal sample uniaxial compression model, obtaining test simulation data, correcting parameters through a back analysis method, obtaining target material parameters of the target phanerocrystalline rock, and constructing a rock crystal medium blasting construction FDEM model; and performing blasting simulation on the target phanerocrystalline intrusive rock according to the rock crystal medium blasting construction FDEM model and blasting parameters. Therefore, the finite-discrete coupling blasting simulation method can not only reduce the cost and time consumption, but also quantitatively represent the dynamic response law and damage mechanism of the rock mass in the blasting process.

Description

A method and apparatus for simulating the blasting of phanerocrystalline intrusive rocks Technical Field

[0001] This application relates to the field of rock simulation blasting technology, and in particular to a blasting simulation method and apparatus for phanerocrystalline intrusive rocks. Background Technology

[0002] In the fields of rock mechanics and blasting engineering, phanerocrystalline intrusive rocks (such as granite, diorite, and gabbro) are renowned for their complex physical and mechanical properties, high hardness, and high compressive strength. These rocks are commonly used in engineering construction, such as mining, tunnel excavation, and hydroelectric power station foundation construction. However, due to the crystalline nature of phanerocrystalline intrusive rocks, their physical characteristics, such as grain size, grain boundary distribution, and internal defects, are complex and highly heterogeneous, making it difficult to accurately predict their mechanical response and failure modes under blasting loads. This presents a significant challenge to blasting construction design and control.

[0003] Specifically, the tightly bound grains within phanerocrystalline intrusive rocks result in extremely high compressive strength, but also increase the energy requirements for blasting operations. Furthermore, due to the distribution and orientation of the crystal structure, phanerocrystalline intrusive rocks exhibit different mechanical properties in different directions, further increasing the complexity of crack propagation during blasting. In addition, during blasting, phanerocrystalline intrusive rocks may experience both intergranular failure (fracture at grain boundaries) and transgranular failure (fracture within the crystal). This makes it difficult to quantitatively characterize the fracture modes of the rock under stress wave action using simple empirical formulas.

[0004] Currently, research on blasting operations for phanerocrystalline intrusive rocks mainly relies on empirical formulas and field tests. However, these methods are not only costly and time-consuming in practical engineering applications, but also struggle to quantify and characterize the dynamic response and failure mechanisms of the rock mass during the blasting process.

[0005] In addition, in recent years, some researchers have attempted to use the finite element method (FEM) to simulate the elastic domain and the discrete element method (DEM) to simulate the plastic and brittle deformation regions, setting contact at the interface to perform FEM-DEM co-simulation of the rock blasting process (e.g., CN115859714B discloses a method for simulating the entire rock blasting process based on FEM-DEM co-simulation). However, since the FEM-DEM coupling method is a technique that combines two numerical methods, the two drastically different assumptions it makes about the continuity and discreteness of the same material can easily lead to problems such as unclear physical meaning and limited model applicability. Summary of the Invention

[0006] In view of this, the purpose of this application is to provide a blasting simulation method and apparatus for phanerocrystalline intrusive rocks. By designing a corresponding finite-discrete coupled blasting simulation method for phanerocrystalline intrusive rocks, not only can the cost and time consumption be reduced, but the dynamic response law and failure mechanism of the rock mass during the blasting process can also be quantitatively characterized.

[0007] This application provides a method for simulating the blasting of phanerocrystalline intrusive rocks, the blasting simulation method comprising:

[0008] Obtain the macroscopic and microscopic physical parameters of the target phanerocrystalline intrusive rock;

[0009] Based on the microscopic physical parameters of the target phanerocrystalline intrusive rock, an initial uniaxial compression model of the rock crystal sample, including mineral composition properties and crystal structure characteristics, was constructed.

[0010] Cohesive elements are generated at the contact interface between grains in the initial uniaxial compression model of the rock crystal sample to obtain the uniaxial compression model of the target rock crystal sample.

[0011] A uniaxial compression simulation test was conducted on the uniaxial compression model of the target rock crystal sample to obtain the experimental simulation data of the target phanerocrystalline intrusive rock;

[0012] Based on the experimental simulation data and the target data in the macroscopic physical parameters, the mineral grain interface damage and fracture mechanical parameters in the uniaxial compression model of the target rock crystal sample are corrected by the inverse analysis method to obtain the target material parameters of the target crystalline material.

[0013] Based on the target material parameters of the target crystalline intrusive rock, an FDEM model for blasting construction in the rock crystal medium is constructed.

[0014] Based on the FDEM model and blasting parameters of the rock crystalline medium blasting construction, a blasting simulation was performed on the target phanerocrystalline intrusive rock.

[0015] Optionally, the step of constructing an initial uniaxial compression model of the rock crystal sample, including mineral composition properties and crystal structure characteristics, based on the microscopic physical parameters of the target phanerocrystalline intrusive rock, includes:

[0016] Based on the microscopic physical parameters of the target phanerocrystalline intrusive rock, mineral-related information and crystal characteristic information in the target phanerocrystalline intrusive rock are determined; wherein, the crystal characteristic information includes: crystal structure, grain size, and crystal orientation distribution;

[0017] Based on the crystal characteristics of the target phanerocrystalline intrusive rock, the Voronoi algorithm is used to determine the grain distribution in the model to be constructed; wherein, each grain is represented as a hexahedral solid unit.

[0018] Based on the mineral-related information of the target phanerocrystalline intrusive rock, attribute values ​​are assigned to each grain to determine the initial uniaxial compression model of the rock crystal sample.

[0019] Optionally, generating cohesive elements at the contact interfaces between grains in the uniaxial compression model of the initial rock crystal sample includes:

[0020] The initial rock crystal sample uniaxial compression model was meshed to determine the contact interface between different grains;

[0021] The corresponding damage criteria and cohesion-related parameters are set for the pre-determined constitutive model, and the model is loaded to generate cohesive units at the contact interface between grains.

[0022] Optionally, the step of conducting a uniaxial compression simulation test on the uniaxial compression model of the target rock crystal sample to obtain the experimental simulation data of the target phanerocrystalline intrusive rock includes:

[0023] An axial load was applied to the uniaxial compression model of the target rock crystal sample, and the stress distribution, crack propagation, and failure mode of the rock were analyzed by finite element simulation to obtain experimental simulation data.

[0024] Optionally, the step of conducting a uniaxial compression simulation test on the uniaxial compression model of the target rock crystal sample to obtain the experimental simulation data of the target phanerocrystalline intrusive rock includes:

[0025] Define the material parameters of the solid elements and cohesive elements in the uniaxial compression model of the target rock crystal sample; wherein the material parameters of the solid elements are set by macroscopic physical parameters.

[0026] Boundary conditions are applied to the uniaxial compression model of the target rock crystal specimen after the material parameters are defined, and simulation processing is performed to generate experimental simulation data.

[0027] Optionally, the step of constructing an FDEM model for blasting construction of rock crystalline media based on the target material parameters of the target phanerocrystalline material includes:

[0028] Based on the target material parameters of the target crystalline material, a macroscopic site model and a microscopic model of the surrounding rock adjacent to the borehole are created respectively.

[0029] The adjacent units and nodes in the macroscopic site model and the microscopic surrounding rock model are merged to form the FDEM model for blasting construction in rock crystalline medium.

[0030] Optionally, the damage criteria include an initial damage criterion and a damage evolution criterion.

[0031] This application embodiment also provides a blasting simulation device for phanerocrystalline intrusive rocks, the blasting simulation device comprising:

[0032] The acquisition module is used to acquire the macroscopic and microscopic physical parameters of the target phanerocrystalline intrusive rock;

[0033] The first construction module is used to construct an initial uniaxial compression model of a rock crystal sample, including mineral composition properties and crystal structure characteristics, based on the microscopic physical parameters of the target phanerocrystalline intrusive rock.

[0034] The generation module is used to generate cohesive elements at the contact interface between grains in the initial uniaxial compression model of the rock crystal sample, so as to obtain the target rock crystal sample uniaxial compression model.

[0035] The test module is used to conduct uniaxial compression simulation tests on the uniaxial compression model of the target rock crystal sample and obtain the test simulation data of the target phanerocrystalline intrusive rock.

[0036] The correction module is used to correct the mineral grain interface damage and fracture mechanical parameters in the uniaxial compression model of the target rock crystal sample based on the experimental simulation data and the target data in the macroscopic physical parameters, and to obtain the target material parameters of the target crystalline material.

[0037] The second construction module is used to construct an FDEM model for blasting construction of rock crystal medium based on the target material parameters of the target crystalline intrusive rock.

[0038] The simulation module is used to simulate the blasting of the target phanerocrystalline intrusive rock based on the FDEM model of the rock crystalline medium blasting construction and the blasting parameters.

[0039] Optionally, when the first construction module is used to construct an initial uniaxial compression model of a rock crystal sample containing mineral composition properties and crystal structure characteristics based on the microscopic physical parameters of the target phanerocrystalline intrusive rock, the first construction module is used to:

[0040] Based on the microscopic physical parameters of the target phanerocrystalline intrusive rock, mineral-related information and crystal characteristic information in the target phanerocrystalline intrusive rock are determined; wherein, the crystal characteristic information includes: crystal structure, grain size, and crystal orientation distribution;

[0041] Based on the crystal characteristics of the target phanerocrystalline intrusive rock, the Voronoi algorithm is used to determine the grain distribution in the model to be constructed; wherein, each grain is represented as a hexahedral solid unit.

[0042] Based on the mineral-related information of the target phanerocrystalline intrusive rock, attribute values ​​are assigned to each grain to determine the initial uniaxial compression model of the rock crystal sample.

[0043] Optionally, when the generation module is used to generate cohesive elements at the contact interfaces between grains in the uniaxial compression model of the initial rock crystal sample, the generation module is used to:

[0044] The initial rock crystal sample uniaxial compression model was meshed to determine the contact interface between different grains;

[0045] The corresponding damage criteria and cohesion-related parameters are set for the pre-determined constitutive model, and the model is loaded to generate cohesive units at the contact interface between grains.

[0046] Optionally, when the test module is used to perform uniaxial compression simulation tests on the uniaxial compression model of the target rock crystal sample to obtain the test simulation data of the target phanerocrystalline intrusive rock, the test module is used for:

[0047] An axial load was applied to the uniaxial compression model of the target rock crystal sample, and the stress distribution, crack propagation, and failure mode of the rock were analyzed by finite element simulation to obtain experimental simulation data.

[0048] Optionally, when the test module is used to perform uniaxial compression simulation tests on the uniaxial compression model of the target rock crystal sample to obtain the test simulation data of the target phanerocrystalline intrusive rock, the test module is used for:

[0049] Define the material parameters of the solid elements and cohesive elements in the uniaxial compression model of the target rock crystal sample; wherein the material parameters of the solid elements are set by macroscopic physical parameters.

[0050] Boundary conditions are applied to the uniaxial compression model of the target rock crystal specimen after the material parameters are defined, and simulation processing is performed to generate experimental simulation data.

[0051] Optionally, when the second construction module is used to construct an FDEM model for blasting construction of rock crystalline media based on the target material parameters of the target phanerocrystalline material, the second construction module is used to:

[0052] Based on the target material parameters of the target crystalline material, a macroscopic site model and a microscopic model of the surrounding rock adjacent to the borehole are created respectively.

[0053] The adjacent units and nodes in the macroscopic site model and the microscopic surrounding rock model are merged to form the FDEM model for blasting construction in rock crystalline medium.

[0054] Optionally, the damage criteria include an initial damage criterion and a damage evolution criterion.

[0055] This application also provides an electronic device, including: a processor, a memory, and a bus. The memory stores machine-readable instructions executable by the processor. When the electronic device is running, the processor communicates with the memory via the bus. When the machine-readable instructions are executed by the processor, the steps of the explosion simulation method described above are performed.

[0056] This application also provides a computer-readable storage medium storing a computer program, which, when executed by a processor, performs the steps of the explosion simulation method described above.

[0057] This application provides a method and apparatus for simulating the blasting of phanerocrystalline intrusive rocks. The blasting simulation method includes: acquiring macroscopic and microscopic physical parameters of the target phanerocrystalline intrusive rock; constructing an initial uniaxial compression model of a rock crystal sample containing mineral composition properties and crystal structure characteristics based on the microscopic physical parameters of the target phanerocrystalline intrusive rock; generating cohesive units at the contact interfaces between grains in the initial uniaxial compression model of the rock crystal sample to obtain a uniaxial compression model of the target rock crystal sample; conducting a uniaxial compression simulation test on the uniaxial compression model of the target rock crystal sample to obtain experimental simulation data of the target phanerocrystalline intrusive rock; correcting the mineral grain interface damage and fracture mechanics parameters in the uniaxial compression model of the target rock crystal sample using a back analysis method based on the experimental simulation data and the target data in the macroscopic physical parameters to obtain the target material parameters of the target phanerocrystalline rock; constructing a rock crystal medium blasting construction FDEM model based on the target material parameters of the target phanerocrystalline intrusive rock; and simulating the blasting of the target phanerocrystalline intrusive rock based on the rock crystal medium blasting construction FDEM model and the blasting parameters.

[0058] Thus, by comprehensively acquiring and analyzing macroscopic and microscopic physical parameters, this application constructs a crystal sample model that more realistically reflects the internal structural characteristics of the target phanerocrystalline intrusive rock, thereby improving the accuracy of the simulation results. Furthermore, through uniaxial compression simulation experiments and back analysis methods, the grain interface damage and fracture mechanical parameters are corrected, making the final target material parameters closer to the mechanical properties of real rocks, ensuring the accuracy of blasting simulation. Constructing an FDEM (Finite Discrete Element Method) model of the rock crystalline medium can effectively handle the complex coupling problem of rock mass fracture and particle movement, improving the construction efficiency and simulation performance of the blasting construction model. Combining blasting parameters with simulation analysis can intuitively demonstrate the dynamic response law of the rock mass, thus providing an optimization basis for actual blasting design and reducing resource waste and unnecessary cost input in actual experiments. Moreover, the FDEM (Finite-Discrete Coupled Method) in this scheme utilizes a cohesive element model to achieve an automatic transition from continuous to discontinuous deformation of brittle materials such as rock and concrete, thereby realizing the natural spatial and temporal coupling of FEM and DEM. Furthermore, through microstructural response analysis at the crystal level, it is possible to quantitatively describe the evolution of microcracks, interface failure mechanisms, and crack propagation patterns in rock masses during blasting, providing a scientific basis for studying the failure behavior of phanerocrystalline intrusive rocks.

[0059] Therefore, this method achieves full-process characterization of phanerocrystalline intrusive rocks from microstructure to macroscopic failure behavior with low cost and resource consumption, and provides an efficient and scientific technical means for rock blasting engineering.

[0060] To make the above-mentioned objectives, features and advantages of this application more apparent and understandable, preferred embodiments are described below in detail with reference to the accompanying drawings. Attached Figure Description

[0061] To more clearly illustrate the technical solutions of the embodiments of this application, the accompanying drawings used in the embodiments will be briefly introduced below. It should be understood that the following drawings only show some embodiments of this application and should not be regarded as a limitation of the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort.

[0062] Figure 1 is a flowchart of a blasting simulation method for phanerocrystalline intrusive rocks provided in an embodiment of this application;

[0063] Figure 2 is a schematic diagram of a uniaxial compression model of an initial rock crystal sample provided in this application;

[0064] Figure 3 is a schematic diagram of a uniaxial compression model of a target rock crystal sample provided in this application;

[0065] Figure 4 shows the results of a uniaxial compression test on a rock provided in this application;

[0066] Figure 5 is a rock blasting FDEM model diagram provided in this application;

[0067] Figure 6 is a schematic diagram of a rock blasting result provided in this application;

[0068] Figure 7 is a schematic diagram of the structure of a blasting simulation device for phanerocrystalline intrusive rocks provided in an embodiment of this application;

[0069] Figure 8 is a schematic diagram of the structure of an electronic device provided in an embodiment of this application. Detailed Implementation

[0070] To make the objectives, technical solutions, and advantages of the embodiments of this application clearer, the technical solutions of the embodiments of this application will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this application, and not all embodiments. The components of the embodiments of this application described and shown in the accompanying drawings can generally be arranged and designed in various different configurations. Therefore, the following detailed description of the embodiments of this application provided in the accompanying drawings is not intended to limit the scope of the claimed application, but merely represents selected embodiments of this application. Based on the embodiments of this application, every other embodiment obtained by those skilled in the art without inventive effort falls within the scope of protection of this application.

[0071] In the fields of rock mechanics and blasting engineering, phanerocrystalline intrusive rocks (such as granite, diorite, and gabbro) are renowned for their complex physical and mechanical properties, high hardness, and high compressive strength. These rocks are commonly used in engineering construction, such as mining, tunnel excavation, and hydroelectric power station foundation construction. However, due to the crystalline nature of phanerocrystalline intrusive rocks, their physical characteristics, such as grain size, grain boundary distribution, and internal defects, are complex and highly heterogeneous, making it difficult to accurately predict their mechanical response and failure modes under blasting loads. This presents a significant challenge to blasting construction design and control.

[0072] Specifically, the tightly bound grains within phanerocrystalline intrusive rocks result in extremely high compressive strength, but also increase the energy requirements for blasting operations. Furthermore, due to the distribution and orientation of the crystal structure, phanerocrystalline intrusive rocks exhibit different mechanical properties in different directions, further increasing the complexity of crack propagation during blasting. In addition, during blasting, phanerocrystalline intrusive rocks may experience both intergranular failure (fracture at grain boundaries) and transgranular failure (fracture within the crystal). This makes it difficult to quantitatively characterize the fracture modes of the rock under stress wave action using simple empirical formulas.

[0073] Currently, research on blasting operations for phanerocrystalline intrusive rocks mainly relies on empirical formulas and field tests. However, these methods are not only costly and time-consuming in practical engineering applications, but also struggle to quantify and characterize the dynamic response and failure mechanisms of the rock mass during the blasting process.

[0074] In recent years, with the improvement of computer computing efficiency and the development of numerical simulation technology, numerical simulation methods based on finite-discrete coupling (FDEM) have gradually shown unique advantages in the study of the mechanical behavior of rock media during construction. The FDEM method combines the accuracy of the finite element method in continuous media simulation with the physical similarity of the discrete element method in fracture and discontinuous media simulation. By inserting a layer of cohesive elements with a thickness of 0 at the interface of the solid model elements through the pre-set crack edges (faces), it can reveal the damage and fracture behavior of continuous media under external forces from a physical perspective. Therefore, it is suitable for simulating the blasting process of engineering rock masses and exploring its mechanical mechanisms.

[0075] Currently, there are still many challenges in using continuous-discrete coupling technology to simulate the explosive mechanics behavior and blasting process of phanerocrystalline intrusive rocks with obvious crystal characteristics and significant anisotropy. For example, how to simulate the intergranular and transgranular fracture characteristics of rocks during blasting, and how to consider the influence of heterogeneous factors such as rock grain size, mineral composition and content on blasting effect.

[0076] Based on this, the present application provides a blasting simulation method and apparatus for phanerocrystalline intrusive rocks. By designing a finite-discrete coupled blasting simulation method, not only can the cost and time consumption be reduced, but the dynamic response law and failure mechanism of the rock mass during the blasting process can also be quantitatively characterized.

[0077] Please refer to Figure 1, which is a flowchart of a blasting simulation method for phanerocrystalline intrusive rocks provided in an embodiment of this application. As shown in Figure 1, the blasting simulation method provided in this embodiment includes:

[0078] S101. Obtain the macroscopic and microscopic physical parameters of the target phanerocrystalline intrusive rock.

[0079] S102. Based on the microscopic physical parameters of the target phanerocrystalline intrusive rock, construct an initial uniaxial compression model of the rock crystal sample, including mineral composition properties and crystal structure characteristics.

[0080] S103. Generate cohesive elements at the contact interface between grains in the initial uniaxial compression model of the rock crystal sample to obtain the uniaxial compression model of the target rock crystal sample.

[0081] S104. Conduct a uniaxial compression simulation test on the uniaxial compression model of the target rock crystal sample to obtain the experimental simulation data of the target phanerocrystalline intrusive rock.

[0082] S105. Based on the experimental simulation data and the target data in the macroscopic physical parameters, the mineral grain interface damage and fracture mechanical parameters in the uniaxial compression model of the target rock crystal sample are corrected by the inverse analysis method to obtain the target material parameters of the target crystalline material.

[0083] S106. Based on the target material parameters of the target crystalline intrusive rock, construct an FDEM model for rock crystal medium blasting construction.

[0084] S107. Based on the FDEM model of rock crystal medium blasting construction and blasting parameters, blasting simulation is performed on the target phanerocrystalline intrusive rock.

[0085] Regarding step S101, the phanerocrystalline intrusive rock refers to rocks with a distinct crystalline structure formed after magma intrudes into the Earth's crust and cools, such as granite, diorite, and gabbro.

[0086] The macroscopic physical parameters refer to the mechanical and physical properties of the rock at the overall scale, used to describe the overall performance of the rock. For example, the macroscopic parameters may include rock density, compressive strength, tensile strength, elastic modulus, Poisson's ratio, etc.

[0087] The microphysical parameters refer to the mineral composition, structure, and interface characteristics of rocks at the microscale, used to describe the detailed internal properties of rocks. For example, the microphysical parameters may include mineral composition, crystal structure, grain size, and crystal orientation distribution.

[0088] The macroscopic physical parameters can be obtained by collecting macroscopic characteristics of phanerocrystalline intrusive rocks through experiments or existing data. The microscopic physical parameters can be obtained by means of microscopy, X-ray diffraction (XRD), etc.

[0089] For step S102, for example, this step may include establishing an initial model with mineral grains, crystal boundaries and mineral distribution using numerical simulation tools (such as DEM or FEM) based on microscopic physical parameters, i.e., determining the initial uniaxial compression model of the rock crystal sample.

[0090] Here, the mineral composition properties may include the mechanical properties and proportions of various minerals (such as quartz, feldspar, mica, etc.) in the rock. The crystal structure characteristics describe the internal geometry, arrangement, and grain boundary properties of the mineral crystals.

[0091] A uniaxial compression model is a model that simulates the mechanical response of a rock when it is loaded in a single direction. It is used to study the compression behavior and fracture mechanism of rocks.

[0092] In one embodiment provided in this application, the step of constructing an initial uniaxial compression model of a rock crystal sample, including mineral composition properties and crystal structure characteristics, based on the microscopic physical parameters of the target phanerocrystalline intrusive rock, includes:

[0093] S1021. Based on the microscopic physical parameters of the target phanerocrystalline intrusive rock, determine the mineral-related information and crystal characteristic information of the target phanerocrystalline intrusive rock.

[0094] S1022. Based on the crystal characteristics information of the target phanerocrystalline intrusive rock, the Voronoi algorithm is used to determine the distribution state of the grains in the model to be constructed.

[0095] S1023. Based on the mineral-related information of the target phanerocrystalline intrusive rock, assign attribute values ​​to each grain to determine the initial uniaxial compression model of the rock crystal sample.

[0096] For step S1021, the crystal characteristic information includes: crystal structure, grain size, and crystal orientation distribution.

[0097] For example, the crystal structure can be determined by microscopy and crystal diffraction data, which can be used to determine the crystal lattice structure (such as monoclinic, orthorhombic, etc.) of crystals in rocks.

[0098] The grain size can be determined using particle analysis software to statistically analyze the average diameter and distribution range of grains in the rock (e.g., 5–15 mm).

[0099] The spatial distribution of crystal orientation can be analyzed using a pole figure or an orientation distribution function (ODF).

[0100] For step S1022, the representation of each grain in the constructed model is a hexahedral solid unit.

[0101] In this step, when using the Voronoi algorithm to create an initial uniaxial compression model of the rock crystal sample, an example may include: randomly generating seed points based on the grain size distribution of the target rock, with each point representing the center position of a grain; using the Voronoi algorithm to partition the space, generating a random grain distribution to form a model that conforms to the characteristics of real crystals; and assigning each grain to a hexahedral solid element, defining boundary regions between grains.

[0102] Here, grains are represented using hexahedral solid units. The grain morphology distribution is determined using the Voronoi spatial partitioning algorithm, which identifies the spatial location of discrete data points based on the crystal mineral type, grain quantity, and size, thereby constructing a Delaunay triangulation to establish a Thiessen polygon polycrystalline model that most closely approximates the actual morphology of rock grains. This method can directly generate a standard Poisson-Voronoi partitioning model and a grain growth model with wider gradation and higher sphericity using the grain partitioning module of the open-source software Neper.

[0103] Regarding step S1023, the rules followed when assigning attribute values ​​include: randomly assigning the mineral type of the grains according to the proportion of mineral composition, and assigning different physical properties to each mineral grain according to its mechanical parameters (such as elastic modulus and Poisson's ratio).

[0104] For example, please refer to Figure 2, which is a schematic diagram of a uniaxial compression model of an initial rock crystal sample provided in this application. As shown in Figure 2, the uniaxial test sample model is 50 mm high and 25 mm in diameter, and consists of 1382 grains. Among them, feldspar accounts for about 65%, with an average grain size of 3.5 mm; quartz accounts for about 35%, with an average grain size of 2.5 mm.

[0105] Regarding step S103, the intergranular contact interface refers to the boundary region where mineral grains in a rock come into contact with each other, which is the main site for crack initiation.

[0106] Cohesive element: In numerical simulation, this element is used to describe the strength of material interfaces and can simulate the tensile, shear, and failure behavior of interfaces.

[0107] In one embodiment provided in this application, generating cohesive elements at the contact interfaces between grains in the uniaxial compression model of the initial rock crystal sample includes:

[0108] S1031. Mesh the initial rock crystal sample uniaxial compression model to determine the contact interface between different grains.

[0109] S1032. Set the corresponding damage criteria and cohesion-related parameters for the pre-determined constitutive model, and load the model to generate cohesion units at the contact interface between grains.

[0110] For step S1031, in this step, each grain in the initial rock crystal sample model is divided into finite element meshes, the complex grain geometry is discretized into a finite number of mesh elements, the mesh division results are analyzed, and the contact area between each two adjacent grains is identified.

[0111] The interface shape is usually a complex curved surface or a plane. Based on the node distribution of the contact points or contact surfaces, the position and area of ​​the contact interface are recorded.

[0112] It should be noted that meshing discretizes complex geometries into simple element meshes (such as hexahedrons or tetrahedrons) to facilitate numerical calculations in finite element analysis. In this application, the meshing is divided into hexahedral contact interfaces, which refer to the boundary regions between adjacent grains and are crucial for the transfer of mechanical forces.

[0113] In step S1032, a suitable constitutive model for the phanerocrystalline intrusive rock (such as a traction-separation bilinear constitutive model) is selected to describe the stress-strain relationship of the material under stress. Damage criteria are defined, for example, based on shear strength or tensile failure strength, and interface failure conditions are set (such as maximum stress or energy release rate exceeding a threshold). In the identified contact interface region, cohesive elements are inserted, and physical parameters are assigned to the cohesive elements, including cohesive strength, interfacial elastic modulus, and damage parameters. The stiffness type of the cohesive elements is specified as traction type. For example, it may include three fracture criteria: opening (Type I), slip (Type II), and tearing (Type III).

[0114] Among them, the constitutive model is a mathematical model that describes the stress and strain relationship of a material under the action of external forces, and it is the foundation of rock mechanics analysis.

[0115] Damage criteria are conditions used to determine when damage or failure occurs in a material or interface, such as the maximum tensile stress criterion.

[0116] Cohesive elements are a type of finite element used to simulate the cohesive mechanical behavior at contact interfaces, representing the shearing, stretching, and failure processes at the interface. This embodiment can generate cohesive elements at crystal interfaces and within the crystal.

[0117] For example, the damage criteria may include an initial damage criterion and a damage evolution criterion.

[0118] The initial damage criterion can be selected from four damage criteria: Maxe Damage, Maxs Damage, Quads Damage, and Quade Damage. The damage evolution criterion can be selected from either displacement-based or energy-based damage evolution criteria. This scheme can choose either the Maxs initial damage criterion or the displacement-based damage evolution criterion.

[0119] For example, please refer to Figure 3, which is a schematic diagram of a uniaxial compression model of a target rock crystal sample provided in this application. As shown in Figure 3, the grains are divided into 306,266 C3D4 solid elements and 111,295 COH3D6 cohesive elements. Material parameters are assigned to the grains and grain boundary elements that characterize different minerals.

[0120] Continuing with step S103, in one optional implementation, when creating the uniaxial compression model of the target rock crystal sample, the process includes: first, using Neper software to generate a polycrystalline model containing grain boundary cohesive units (initial uniaxial compression model of the rock crystal sample); then, using the global cohesive unit plugin to generate an intragranular cohesive unit model (cohesive unit). The strength of the grain boundary cohesive units is slightly lower than that of the intragranular cohesive units to match the physical and mechanical process of phanerocrystalline rocks first undergoing intergranular fracture under external load, followed by transgranular fracture.

[0121] For step S104, this step includes: applying uniaxial compression loading to the uniaxial compression model of the target rock crystal sample, recording the rock change data during the loading process, and obtaining experimental simulation data. For example, the experimental simulation data may include data such as stress-strain curves, crack distribution, and failure modes.

[0122] Among them, the uniaxial compression simulation test is a numerical simulation method to simulate the mechanical behavior of rock under compression in a single direction and observe its failure mode and mechanical response.

[0123] Experimental simulation data is the output data generated during the numerical simulation process, including information such as stress, strain, and crack propagation.

[0124] In one embodiment provided in this application, the step of performing a uniaxial compression simulation test on the uniaxial compression model of the target rock crystal sample to obtain the test simulation data of the target phanerocrystalline intrusive rock includes: applying an axial load to the uniaxial compression model of the target rock crystal sample, and analyzing the stress distribution, crack propagation and failure mode of the rock through finite element simulation to obtain the test simulation data.

[0125] Here, a vertical axial load (usually applied by displacement control or force control) is applied to the top or bottom of the uniaxial compression model of the target rock crystal specimen. The load can be applied gradually (quasi-static loading) or at a certain speed (dynamic loading) to simulate the actual uniaxial compression test conditions. Then, the finite element method (FEM) is used to perform numerical analysis on the loading process to obtain test simulation data.

[0126] The experimental simulation data may specifically include the following key data:

[0127] Stress-strain curve: Records the stress-strain relationship of the model during loading.

[0128] Crack data: including the number, location, propagation path, and length of cracks.

[0129] Failure mode: The overall failure mode of the rock model, such as slip along the interface or grain fragmentation.

[0130] It should be noted that these simulation data are used to further correct model parameters (such as interface mechanical parameters) or verify the accuracy of the model.

[0131] In another embodiment provided in this application, the step of performing a uniaxial compression simulation test on the uniaxial compression model of the target rock crystal sample to obtain the experimental simulation data of the target phanerocrystalline intrusive rock includes: defining the material parameters of the solid elements and cohesive elements in the uniaxial compression model of the target rock crystal sample; wherein the material parameters of the solid elements are set by macroscopic physical parameters; applying boundary conditions to the uniaxial compression model of the target rock crystal sample after the material parameters are defined and performing simulation processing to generate experimental simulation data.

[0132] First, based on the macroscopic physical parameters of the target rock (such as density, elastic modulus, Poisson's ratio, etc.), appropriate mechanical properties are assigned to the solid elements to ensure that the overall mechanical behavior of the model conforms to the properties of real rocks. Simultaneously, the material parameters of the cohesive elements (including cohesive strength, fracture toughness, interfacial elastic modulus, etc.) are determined. Then, constraints and loading conditions are applied to the target rock crystal sample model. Finally, the compression process is simulated using finite element methods or explicit dynamic methods (such as discrete element method or explicit finite difference method), and relevant data are recorded and experimental simulation data is acquired during the simulation.

[0133] For example, constraint methods may include: bottom fixed constraint (preventing displacement or rotation of the bottom of the model and providing symmetry conditions).

[0134] The loading conditions may include: displacement-controlled loading (gradually applying displacement, typically increasing at a rate of several millimeters per second), force-controlled loading (gradually applying force, typically increasing at a rate of several megapascals per second), etc.

[0135] The experimental simulation data may include: stress-strain relationship, crack propagation path, damage distribution, etc.

[0136] For example, please refer to Figure 4, which shows the uniaxial compression test results of a rock according to this application. As shown in Figure 4, the test simulation data is a stress-strain curve, which illustrates the deformation of the uniaxial compression specimen. This test simulation data is used to compare and analyze the stress-strain curves of numerical simulation and laboratory tests. The parameters of the cohesive element constitutive model are then corrected in reverse using the laboratory test results as a guideline, ensuring that the numerical model results match the laboratory test results.

[0137] For step S105, for example, this step includes: adjusting the mineral grain interface parameters (such as strength and cohesion) in the model using back analysis methods (such as optimization algorithms) based on the differences between the simulation test data and the actual experimental data. This adjustment is iteratively repeated until the model simulation results are highly consistent with the experimental results, thus obtaining accurate target material parameters.

[0138] Here, the calibration specifically involves calibrating the parameters of the cohesion constitutive model. For example, the method for calibrating the parameters of the cohesion constitutive model is as follows: First, select initial material parameters based on literature and calculate to obtain the initial stress-strain relationship curve of the rock sample under compression. Then, compare and analyze the calculated results with the experimental results, and, considering the sensitivity and influence of the material parameters on the results, gradually increase or decrease the parameters accordingly. Finally, ensure that the calculated stress-strain curve is similar in shape and amplitude to the experimental curve. In this example, the simulation results can reflect both intergranular fracture and transgranular fracture during the sample compression failure process; furthermore, it can simulate the full stress-strain curve and identify the post-peak failure strength of the rock.

[0139] Among them, the inverse analysis method is a method that adjusts the model parameters to make the numerical simulation results gradually approach the experimental results, and is often used for parameter calibration.

[0140] The target material parameters are the final, corrected parameters used to characterize the mechanical and physical properties of the target rock.

[0141] For step S106, for example, this step may include: using the FDEM (Finite Discrete Element Method) tool, and combining the target material parameters to establish a blasting model containing a mineral crystalline medium. Here, the model may incorporate numerical descriptions of the rock fracture process, such as stress wave propagation and crack propagation.

[0142] Among them, the rock crystal medium is the microstructure medium composed of mineral crystals and grains in the rock.

[0143] The FDEM model is a finite discrete element method model that can simultaneously simulate the continuous medium behavior and fracture process of rocks.

[0144] In one embodiment provided in this application, constructing an FDEM model for blasting construction of rock crystal medium based on the target material parameters of the target phanerocrystalline substance includes:

[0145] Based on the target material parameters of the target crystalline material, a macroscopic site model and a microscopic model of the surrounding rock adjacent to the borehole are created respectively; adjacent units and nodes in the macroscopic site model and the microscopic model of the surrounding rock are merged to form a rock crystalline medium blasting construction FDEM model.

[0146] This step may include: First, determining the target material parameters of the target phanerocrystalline intrusive rock (for example, the rock's elastic modulus, Poisson's ratio, density, tensile strength, compressive strength, cohesive strength, and fracture toughness, etc.), then selecting a suitable numerical calculation method (such as finite element method, finite difference method, or boundary element method, etc.) to build the model and create the macroscopic site model, and simultaneously selecting a suitable microscopic simulation method (such as discrete element method, microscopic finite element method, etc.) to build the model and create the surrounding rock microscopic model; then, assembling the two into a whole and merging adjacent elements and nodes to obtain the FDEM model for blasting construction of rock crystalline media.

[0147] It should be noted that the FDEM model for blasting operations in rock crystalline media is a macro-micro coupled model. This model includes not only crystal geometry information but also grain boundaries and cohesive units within the crystal generated by Neper software. Furthermore, this model can be used to simulate the blasting response of a site, and is particularly suitable for exploring the destructive and damage effects of blasting loads on the surrounding rock.

[0148] For example, please refer to Figure 5, which is an FDEM model diagram of rock blasting provided in this application. As shown in Figure 5, a is the macroscopic model (macroscopic site model), and b is the microscopic model (mesoscopic model of surrounding rock). As shown in Figure 5, the global model is 100m long, 50m wide, and 1m thick, composed of 57,125 C3D6 elements and 57,077 COH3D6 cohesive elements; the sub-model has a planar size of 20cm × 20cm, divided into 61,022 C3D6 elements by 3,756 grains, of which 2,441 are feldspar mineral grains and 1,315 are quartz grains. Cohesive elements are used to bond the grains and the grain boundaries.

[0149] For step S107, blasting parameters (such as explosive charge, hole spacing, and delay) are introduced into the FDEM model to simulate the blasting process. The dynamic response of the rock mass during blasting is then obtained, which may include results such as crack distribution, fracture range, and flyrock trajectory.

[0150] For example, this step may include: setting blasting parameters such as explosive type, charge quantity, borehole layout, and detonation sequence according to blasting design requirements; simulating the processes of rock fragmentation, crack propagation, stress wave propagation, and the effects of explosive gases under blasting loads. The blasting effect can be evaluated by analyzing the simulation results, including the extent, morphology, and distribution characteristics of the pulverized and fractured zones, as well as the impact of the blast on the surrounding rock.

[0151] In this way, the blasting test of phanerocrystalline intrusive bodies can be simulated, and the internal stress state, rock fracture and fragmentation development process of the rock mass can be analyzed.

[0152] For example, please refer to Figure 6, which is a schematic diagram of a rock blasting result provided in this application. As shown in Figure 6, a is a stress cloud diagram, and b is a fracture distribution diagram. In this way, the blasting test of phanerocrystalline intrusive bodies can be simulated, and the internal stress state of the rock mass, rock fracture, and fragmentation development process can be analyzed simultaneously.

[0153] Thus, by comprehensively acquiring and analyzing macroscopic and microscopic physical parameters, this application constructs a crystal sample model that more realistically reflects the internal structural characteristics of the target phanerocrystalline intrusive rock, thereby improving the accuracy of the simulation results. Furthermore, through uniaxial compression simulation experiments and back analysis methods, the grain interface damage and fracture mechanical parameters are corrected, making the final target material parameters closer to the mechanical properties of real rocks, ensuring the accuracy of blasting simulation. Constructing an FDEM (Finite Discrete Element Method) model of the rock crystalline medium can effectively handle the complex coupling problem of rock mass fracture and particle movement, improving the construction efficiency and simulation performance of the blasting construction model. Analyzing blasting parameters during the simulation process can intuitively demonstrate the dynamic response law of the rock mass, thus providing an optimization basis for actual blasting design and reducing resource waste and unnecessary cost input in actual experiments. In addition, through microstructural response analysis at the crystal level, the evolution of microcracks, interface failure mechanisms, and crack propagation laws of the rock mass during blasting can be quantitatively described, providing a scientific basis for studying the failure behavior of phanerocrystalline intrusive rocks.

[0154] Therefore, this method achieves full-process characterization of phanerocrystalline intrusive rocks from microstructure to macroscopic failure behavior with low cost and resource consumption, and provides an efficient and scientific technical means for rock blasting engineering.

[0155] Based on the same inventive concept, this application also provides a blasting simulation device corresponding to the blasting simulation method. Since the principle of the device in this application is similar to that of the blasting simulation method described above, the implementation of the device can refer to the implementation of the method, and the repeated parts will not be described again.

[0156] Please refer to Figure 7, which is a schematic diagram of the structure of a blasting simulation device for phanerocrystalline intrusive rocks provided in an embodiment of this application. As shown in Figure 7, the blasting simulation device 700 includes:

[0157] The acquisition module 710 is used to acquire the macroscopic and microscopic physical parameters of the target phanerocrystalline intrusive rock;

[0158] The first construction module 720 is used to construct an initial uniaxial compression model of a rock crystal sample, including mineral composition properties and crystal structure characteristics, based on the microscopic physical parameters of the target phanerocrystalline intrusive rock.

[0159] The generation module 730 is used to generate cohesive elements at the contact interface between grains in the initial uniaxial compression model of the rock crystal sample to obtain the target uniaxial compression model of the rock crystal sample.

[0160] The test module 740 is used to conduct a uniaxial compression simulation test on the uniaxial compression model of the target rock crystal sample to obtain the test simulation data of the target phanerocrystalline intrusive rock;

[0161] The correction module 750 is used to correct the mineral grain interface damage and fracture mechanical parameters in the uniaxial compression model of the target rock crystal sample based on the experimental simulation data and the target data in the macroscopic physical parameters, and to obtain the target material parameters of the target crystalline material.

[0162] The second construction module 760 is used to construct an FDEM model for blasting construction of rock crystal medium based on the target material parameters of the target crystalline intrusive rock.

[0163] The simulation module 770 is used to simulate the blasting of the target phanerocrystalline intrusive rock based on the FDEM model of the rock crystalline medium blasting construction and the blasting parameters.

[0164] Optionally, when the first construction module 720 is used to construct an initial uniaxial compression model of a rock crystal sample containing mineral composition properties and crystal structure characteristics based on the microscopic physical parameters of the target phanerocrystalline intrusive rock, the first construction module 720 is used to:

[0165] Based on the microscopic physical parameters of the target phanerocrystalline intrusive rock, mineral-related information and crystal characteristic information in the target phanerocrystalline intrusive rock are determined; wherein, the crystal characteristic information includes: crystal structure, grain size, and crystal orientation distribution;

[0166] Based on the crystal characteristics of the target phanerocrystalline intrusive rock, the Voronoi algorithm is used to determine the grain distribution in the model to be constructed; wherein, each grain is represented as a hexahedral solid unit.

[0167] Based on the mineral-related information of the target phanerocrystalline intrusive rock, attribute values ​​are assigned to each grain to determine the initial uniaxial compression model of the rock crystal sample.

[0168] Optionally, when the generation module 730 is used to generate cohesive elements at the contact interfaces between grains in the uniaxial compression model of the initial rock crystal sample, the generation module 730 is used to:

[0169] The initial rock crystal sample uniaxial compression model was meshed to determine the contact interface between different grains;

[0170] The corresponding damage criteria and cohesion-related parameters are set for the pre-determined constitutive model, and the model is loaded to generate cohesive units at the contact interface between grains.

[0171] Optionally, when the test module 740 is used to perform a uniaxial compression simulation test on the uniaxial compression model of the target rock crystal sample to obtain the test simulation data of the target phanerocrystalline intrusive rock, the test module 740 is used for:

[0172] An axial load was applied to the uniaxial compression model of the target rock crystal sample, and the stress distribution, crack propagation, and failure mode of the rock were analyzed by finite element simulation to obtain experimental simulation data.

[0173] Optionally, when the test module 740 is used to perform a uniaxial compression simulation test on the uniaxial compression model of the target rock crystal sample to obtain the test simulation data of the target phanerocrystalline intrusive rock, the test module 740 is used for:

[0174] Define the material parameters of the solid elements and cohesive elements in the uniaxial compression model of the target rock crystal sample; wherein the material parameters of the solid elements are set by macroscopic physical parameters.

[0175] Boundary conditions are applied to the uniaxial compression model of the target rock crystal specimen after the material parameters are defined, and simulation processing is performed to generate experimental simulation data.

[0176] Optionally, when the second construction module 760 is used to construct an FDEM model for blasting construction of rock crystal medium based on the target material parameters of the target crystalline substance, the second construction module 760 is used to:

[0177] Based on the target material parameters of the target crystalline material, a macroscopic site model and a microscopic model of the surrounding rock adjacent to the borehole are created respectively.

[0178] The adjacent units and nodes in the macroscopic site model and the microscopic surrounding rock model are merged to form the FDEM model for blasting construction in rock crystalline medium.

[0179] Optionally, the damage criteria include an initial damage criterion and a damage evolution criterion.

[0180] Please refer to Figure 8, which is a schematic diagram of the structure of an electronic device provided in an embodiment of this application. As shown in Figure 8, the electronic device 800 includes a processor 810, a memory 820, and a bus 830.

[0181] The memory 820 stores machine-readable instructions that can be executed by the processor 810. When the electronic device 800 is running, the processor 810 and the memory 820 communicate through the bus 830. When the machine-readable instructions are executed by the processor 810, the steps in the method embodiments shown in Figures 1 to 6 above can be performed. For specific implementation methods, please refer to the method embodiments, which will not be repeated here.

[0182] This application also provides a computer-readable storage medium storing a computer program. When the computer program is run by a processor, it can execute the steps in the method embodiments shown in Figures 1 to 6 above. For specific implementation details, please refer to the method embodiments, which will not be repeated here.

[0183] Those skilled in the art will understand that, for the sake of convenience and brevity, the specific working processes of the systems, devices, and units described above can be referred to the corresponding processes in the foregoing method embodiments, and will not be repeated here.

[0184] In the several embodiments provided in this application, it should be understood that the disclosed systems, apparatuses, and methods can be implemented in other ways. The apparatus embodiments described above are merely illustrative. For example, the division of units is only a logical functional division, and in actual implementation, there may be other division methods. Furthermore, multiple units or components may be combined or integrated into another system, or some features may be ignored or not executed. Additionally, the shown or discussed mutual couplings, direct couplings, or communication connections may be through some communication interfaces; indirect couplings or communication connections between devices or units may be electrical, mechanical, or other forms.

[0185] The units described as separate components may or may not be physically separate. The components shown as units may or may not be physical units; that is, they may be located in one place or distributed across multiple network units. Some or all of the units can be selected to achieve the purpose of this embodiment according to actual needs.

[0186] In addition, the functional units in the various embodiments of this application can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit.

[0187] If the aforementioned functions are implemented as software functional units and sold or used as independent products, they can be stored in a processor-executable, non-volatile, computer-readable storage medium. Based on this understanding, the technical solution of this application, in essence, or the part that contributes to the prior art, or a portion of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) to execute all or part of the steps of the methods described in the various embodiments of this application. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0188] Finally, it should be noted that the above-described embodiments are merely specific implementations of this application, used to illustrate the technical solutions of this application, and not to limit them. The scope of protection of this application is not limited thereto. Although this application has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that any person skilled in the art can still modify or easily conceive of changes to the technical solutions described in the foregoing embodiments, or make equivalent substitutions for some of the technical features, within the scope of the technology disclosed in this application. Such modifications, changes, or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of this application, and should all be covered within the scope of protection of this application. Therefore, the scope of protection of this application should be determined by the scope of the claims.

Claims

1. A method for simulating the blasting of phanerocrystalline intrusive rocks, characterized in that, The blasting simulation method includes: acquiring the macroscopic and microscopic physical parameters of the target phanerocrystalline intrusive rock; constructing an initial uniaxial compression model of the rock crystal sample, including mineral composition properties and crystal structure characteristics, based on the microscopic physical parameters of the target phanerocrystalline intrusive rock; generating cohesive elements at the contact interfaces between grains in the initial uniaxial compression model of the rock crystal sample to obtain the uniaxial compression model of the target rock crystal sample; conducting uniaxial compression simulation tests on the uniaxial compression model of the target rock crystal sample to obtain the experimental simulation data of the target phanerocrystalline intrusive rock; correcting the mineral grain interface damage and fracture mechanical parameters in the uniaxial compression model of the target rock crystal sample using a back analysis method based on the experimental simulation data and the target data in the macroscopic physical parameters to obtain the target material parameters of the target phanerocrystalline rock; constructing a rock crystal medium blasting construction FDEM model based on the target material parameters of the target phanerocrystalline intrusive rock; conducting blasting simulation on the target phanerocrystalline intrusive rock based on the rock crystal medium blasting construction FDEM model and blasting parameters; and constructing a rock crystal medium blasting construction FDEM model based on the target material parameters of the target phanerocrystalline intrusive rock. The physical parameters are used to construct an initial uniaxial compression model of a rock crystal sample, including mineral composition attributes and crystal structure characteristics. This includes: determining mineral-related information and crystal characteristic information in the target phanerocrystalline intrusive rock based on its microscopic physical parameters; wherein the crystal characteristic information includes crystal structure, grain size, and crystal orientation distribution; determining the grain distribution state in the model to be constructed using the Voronoi algorithm based on the crystal characteristic information of the target phanerocrystalline intrusive rock; wherein each grain is represented as a hexahedral solid element; assigning attribute values ​​to each grain based on the mineral-related information of the target phanerocrystalline intrusive rock to determine the initial uniaxial compression model of the rock crystal sample; generating cohesive elements at the contact interfaces between grains in the initial uniaxial compression model of the rock crystal sample, including: meshing the initial uniaxial compression model of the rock crystal sample to determine the contact interfaces between different grains; setting corresponding damage criteria and corresponding cohesive related parameters for the pre-determined constitutive model, and loading the model to generate cohesive elements at the contact interfaces between grains.

2. The blasting simulation method according to claim 1, characterized in that, The step of conducting a uniaxial compression simulation test on the uniaxial compression model of the target rock crystal sample to obtain the test simulation data of the target phanerocrystalline intrusive rock includes: applying an axial load to the uniaxial compression model of the target rock crystal sample, and analyzing the stress distribution, crack propagation and failure mode of the rock through finite element simulation to obtain the test simulation data.

3. The blasting simulation method according to claim 1, characterized in that, The step of conducting a uniaxial compression simulation test on the uniaxial compression model of the target rock crystal sample to obtain the experimental simulation data of the target phanerocrystalline intrusive rock includes: defining the material parameters of the solid elements and cohesive elements in the uniaxial compression model of the target rock crystal sample; wherein the material parameters of the solid elements are set through macroscopic physical parameters; applying boundary conditions to the uniaxial compression model of the target rock crystal sample after the material parameters are defined and performing simulation processing to generate experimental simulation data.

4. The blasting simulation method according to claim 1, characterized in that, The step of constructing a rock crystal medium blasting construction FDEM model based on the target material parameters of the target crystalline substance includes: creating a macroscopic site model and a surrounding rock microscopic model adjacent to the borehole based on the target material parameters of the target crystalline substance; merging adjacent units and nodes in the macroscopic site model and the surrounding rock microscopic model to form a rock crystal medium blasting construction FDEM model.

5. The blasting simulation method according to claim 1, characterized in that, The damage criteria include the initial damage criterion and the damage evolution criterion.

6. A blasting simulation device for phanerocrystalline intrusive rocks, characterized in that, The blasting simulation device includes: an acquisition module for acquiring the macroscopic and microscopic physical parameters of the target phanerocrystalline intrusive rock; a first construction module for constructing an initial uniaxial compression model of the rock crystal sample, including mineral composition properties and crystal structure characteristics, based on the microscopic physical parameters of the target phanerocrystalline intrusive rock; a generation module for generating cohesive elements at the contact interfaces between grains in the initial uniaxial compression model of the rock crystal sample, thereby obtaining the uniaxial compression model of the target rock crystal sample; and an experimental module for conducting uniaxial compression simulation tests on the uniaxial compression model of the target rock crystal sample to obtain the macroscopic and microscopic physical parameters of the target phanerocrystalline intrusive rock. The system includes: a simulation module for rock intrusion test data; a calibration module for correcting the mineral grain interface damage and fracture mechanics parameters in the uniaxial compression model of the target rock crystal sample using back analysis based on the simulation test data and the target data in the macroscopic physical parameters, to obtain the target material parameters of the target crystalline rock; a second construction module for constructing a rock crystal medium blasting FDEM model based on the target material parameters of the target crystalline intrusive rock; and a simulation module for blasting simulation of the target crystalline intrusive rock based on the rock crystal medium blasting FDEM model and blasting parameters. When constructing an initial uniaxial compression model of a rock crystal sample containing mineral composition and crystal structure characteristics based on the microscopic physical parameters of the target phanerocrystalline intrusive rock, the first construction module is used to: determine mineral-related information and crystal characteristic information in the target phanerocrystalline intrusive rock based on the microscopic physical parameters of the target phanerocrystalline intrusive rock; wherein, the crystal characteristic information includes: crystal structure, grain size, and crystal orientation distribution; and, based on the crystal characteristic information of the target phanerocrystalline intrusive rock, use the Voronoi algorithm to determine the distribution state of the grains in the model to be constructed; wherein, the characteristic shape of each grain is... The formula is a hexahedral solid element; based on the mineral information of the target phanerocrystalline intrusive rock, attribute values ​​are assigned to each grain to determine the initial uniaxial compression model of the rock crystal sample; when the generation module is used to generate cohesive elements at the contact interface between grains in the initial uniaxial compression model of the rock crystal sample, the generation module is used to: mesh the initial uniaxial compression model of the rock crystal sample to determine the contact interface between different grains; set the corresponding damage criteria and corresponding cohesion-related parameters for the pre-determined constitutive model, and load the model to realize the generation of cohesive elements at the contact interface between grains.

7. An electronic device, characterized in that, include: The device includes a processor, a memory, and a bus. The memory stores machine-readable instructions executable by the processor. When the electronic device is running, the processor communicates with the memory via the bus, and the machine-readable instructions are executed by the processor to perform the steps of the blasting simulation method as described in any one of claims 1 to 5.

8. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores a computer program that, when executed by a processor, performs the steps of the blasting simulation method as described in any one of claims 1 to 5.

Images