Blasting simulation method and device for crystal intrusive rock

Abstract

The invention provides a blasting simulation method and device for crystal intruding rock, and the method comprises the steps: building an initial rock crystal sample uniaxial compression model according to the mesoscopic physical parameters of target crystal intruding rock; generating a cohesive force unit for a contact interface between crystal grains in the uniaxial compression model to obtain a uniaxial compression model of the target rock crystal sample; performing a uniaxial compression simulation test on the uniaxial compression model of the target rock crystal sample to obtain test simulation data, correcting parameters through a back analysis method to obtain target material parameters of a target crystal, and constructing a rock crystal medium blasting construction FDEM model; and according to the rock crystal medium blasting construction FDEM model and the blasting parameters, blasting simulation is carried out on the target crystal intruding rock. Therefore, according to the scheme, by designing the finite-discrete coupling blasting simulation method, the cost can be reduced, the time consumption can be reduced, and the dynamic response rule and the damage mechanism of the rock mass in the blasting process can be quantitatively represented.

Description

Technical Field

[0001] This application relates to the technical field of rock simulation blasting, and in particular to a blasting simulation method and device for phaneritic intrusive rocks. Background Art

[0002] In the fields of rock mechanics and blasting engineering, phaneritic intrusive rocks (such as granite, diorite, gabbro, etc.) are known for their complex physical and mechanical properties, high hardness, and high compressive strength. These rocks are commonly found in engineering construction, such as mine exploitation, tunnel excavation, and hydropower station foundation construction. However, due to the crystal properties within phaneritic 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 the mechanical response and failure mode under blasting loads. This poses a huge challenge to blasting construction design and control.

[0003] Specifically, because the grains within phaneritic intrusive rocks are tightly bound, it results in extremely high compressive capacity, but at the same time increases the energy requirements in blasting construction. And due to the distribution and directionality of the crystal structure, phaneritic intrusive rocks exhibit different mechanical properties in different directions, further increasing the complexity of crack propagation during blasting. In addition, during the blasting process, phaneritic intrusive rocks may undergo intergranular failure (fracture at grain boundaries) or transgranular failure (fracture within crystals). This makes it difficult to quantitatively characterize the fracture mode of the rock under stress wave action through simple empirical formulas.

[0004] Currently, the research on blasting construction for phaneritic intrusive rocks mainly relies on empirical formulas and on-site tests. However, these methods not only have high costs and long time consumption in practical engineering applications, but also are difficult to quantitatively characterize the dynamic response law and failure mechanism of rock masses during blasting.

[0005] In addition, in recent years, some people have tried to use the finite element method to simulate the elastic domain, the discrete element method to simulate the plastic and brittle deformation regions, and set contacts at the interface to conduct FEM-DEM joint simulation of the rock blasting process (for example, CN115859714B discloses a method for simulating the whole process of rock blasting based on FEM-DEM joint simulation). However, since the finite element and discrete element coupling (FEM-DEM Coupling) method is a technology that combines two numerical methods of finite element and discrete element, the completely different assumptions of continuity and discreteness made for the same material by this method are likely to cause problems such as unclear physical meaning and limited model applicability. Summary of the Invention

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

[0007] An embodiment of the present application provides a blasting simulation method for phanerocrystalline intrusive rock, and the blasting simulation method includes:

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

[0009] According to the mesoscopic physical parameters of the target phanerocrystalline intrusive rock, construct an initial uniaxial compression model of the rock crystal specimen including mineral composition attributes and crystal structure characteristics;

[0010] Generate cohesive force units at the contact interfaces between grains in the initial uniaxial compression model of the rock crystal specimen to obtain a target uniaxial compression model of the rock crystal specimen;

[0011] Conduct a uniaxial compression simulation test on the target uniaxial compression model of the rock crystal specimen to obtain the test simulation data of the target phanerocrystalline intrusive rock;

[0012] According to the test simulation data and the target data in the macroscopic physical parameters, through the inverse analysis method, correct the mineral grain interface damage and fracture mechanics parameters in the target uniaxial compression model of the rock crystal specimen to obtain the target material parameters of the target phanerocrystalline;

[0013] Construct a rock crystal medium blasting construction FDEM model according to the target material parameters of the target phanerocrystalline intrusive rock;

[0014] Conduct a blasting simulation on the target phanerocrystalline intrusive rock according to the rock crystal medium blasting construction FDEM model and the blasting parameters.

[0015] Optionally, the constructing an initial uniaxial compression model of the rock crystal specimen including mineral composition attributes and crystal structure characteristics according to the mesoscopic physical parameters of the target phanerocrystalline intrusive rock includes:

[0016] According to the mesoscopic physical parameters of the target phanerocrystalline intrusive rock, determine the mineral-related information and crystal characteristic information in the target phanerocrystalline intrusive rock; wherein, the crystal characteristic information includes: crystal structure, grain size, and crystal orientation distribution;

[0017] According to 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 representation form of each grain is a hexahedral solid element;

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

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

[0020] Performing mesh division on the uniaxial compression model of the initial rock crystal specimen to determine the contact interfaces between different grains;

[0021] Setting corresponding damage criteria and corresponding cohesive - force - related parameters for a pre - determined constitutive model, and performing model loading to realize the generation of cohesive force units at the contact interfaces between grains.

[0022] Optionally, performing a uniaxial compression simulation test on the uniaxial compression model of the target rock crystal specimen to obtain the test simulation data of the target phanerocrystalline intrusive rock, including:

[0023] Applying an axial load to the uniaxial compression model of the target rock crystal specimen, and analyzing the stress distribution, crack propagation, and failure mode of the rock through finite - element simulation to obtain the test simulation data.

[0024] Optionally, performing a uniaxial compression simulation test on the uniaxial compression model of the target rock crystal specimen to obtain the test simulation data of the target phanerocrystalline intrusive rock, including:

[0025] Defining the material parameters of the solid elements and cohesive force units in the uniaxial compression model of the target rock crystal specimen; where the material parameters of the solid elements are set through macroscopic physical parameters;

[0026] Applying boundary conditions to the uniaxial compression model of the target rock crystal specimen after the definition of material parameters and performing simulation processing to generate the test simulation data.

[0027] Optionally, constructing an FDEM model for rock crystal medium blasting construction according to the target material parameters of the target phanerocrystalline, including:

[0028] According to the target material parameters of the target phanerocrystalline, creating a macroscopic site model and a mesoscopic model of the surrounding rock adjacent to the borehole respectively;

[0029] Merging adjacent units and nodes in the macroscopic site model and the mesoscopic model of the surrounding rock to form an FDEM model for rock crystal medium blasting construction.

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

[0031] The embodiment of the present application also provides a blasting simulation device for phanerocrystalline intrusive rock, and the blasting simulation device includes:

[0032] An acquisition module, configured to acquire the macroscopic physical parameters and mesoscopic physical parameters of the target phanerocrystalline intrusive rock;

[0033] A first construction module, configured to construct an initial uniaxial compression model of a rock crystal specimen including mineral composition attributes and crystal structure characteristics according to the mesoscopic physical parameters of the target phanerocrystalline intrusive rock;

[0034] A generation module, configured to generate cohesive force units for the contact interfaces between grains in the initial uniaxial compression model of the rock crystal specimen to obtain a target uniaxial compression model of the rock crystal specimen;

[0035] A test module, configured to perform a uniaxial compression simulation test on the target uniaxial compression model of the rock crystal specimen to obtain test simulation data of the target phanerocrystalline intrusive rock;

[0036] A correction module, configured to correct the mineral grain interface damage and fracture mechanics parameters in the target uniaxial compression model of the rock crystal specimen by an inverse analysis method according to the test simulation data and the target data in the macroscopic physical parameters to obtain the target material parameters of the target phanerocrystalline;

[0037] A second construction module, configured to construct a blasting construction FDEM model of a rock crystal medium according to the target material parameters of the target phanerocrystalline intrusive rock;

[0038] A simulation module, configured to perform a blasting simulation on the target phanerocrystalline intrusive rock according to the blasting construction FDEM model of the rock crystal medium and the blasting parameters.

[0039] Optionally, when the first construction module is configured to construct an initial uniaxial compression model of a rock crystal specimen including mineral composition attributes and crystal structure characteristics according to the mesoscopic physical parameters of the target phanerocrystalline intrusive rock, the first construction module is configured to:

[0040] Determine the mineral-related information and crystal characteristic information in the target phanerocrystalline intrusive rock according to the mesoscopic physical parameters of the target phanerocrystalline intrusive rock; wherein, the crystal characteristic information includes: crystal structure, grain size, and crystal orientation distribution;

[0041] Determine the distribution state of grains in the model to be constructed by using the Voronoi algorithm according to the crystal characteristic information of the target phanerocrystalline intrusive rock; wherein, the representation form of each grain is a hexahedral solid element;

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

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

[0044] Perform mesh division on the uniaxial compression model of the initial rock crystal specimen to determine the contact interfaces between different grains;

[0045] Set the corresponding damage criterion and the corresponding cohesive - force - related parameters for the constitutive model determined in advance, and perform model loading to realize the generation of cohesive force units at the contact interfaces between grains.

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

[0047] Apply an axial load to the uniaxial compression model of the target rock crystal specimen, and analyze the stress distribution, crack propagation, and failure mode of the rock through finite - element simulation to obtain the test simulation data.

[0048] Optionally, when the test module is used to perform a uniaxial compression simulation test on the uniaxial compression model of the target rock crystal specimen 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 force units in the uniaxial compression model of the target rock crystal specimen; among them, the material parameters of the solid elements are set through macroscopic physical parameters;

[0050] Apply boundary conditions to the uniaxial compression model of the target rock crystal specimen after the definition of material parameters and perform simulation processing to generate test simulation data.

[0051] Optionally, when the second construction module is used to construct the FDEM model for rock crystal medium blasting construction according to the target material parameters of the target phanerocrystalline, the second construction module is used for:

[0052] According to the target material parameters of the target phanerocrystalline, create a macroscopic site model and a mesoscopic surrounding - rock model adjacent to the borehole respectively;

[0053] Merge the adjacent units and nodes in the macroscopic site model and the mesoscopic surrounding - rock model to form the FDEM model for rock crystal medium blasting construction.

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

[0055] An embodiment of the present 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 runs, the processor communicates with the memory through the bus. When the machine-readable instructions are executed by the processor, the steps of the blasting simulation method as described above are executed.

[0056] An embodiment of the present application also provides a computer-readable storage medium. A computer program is stored on the computer-readable storage medium. When the computer program is run by a processor, the steps of the blasting simulation method as described above are executed.

[0057] A blasting simulation method and device for phanerocrystalline intrusive rock provided by an embodiment of the present application. The blasting simulation method includes: obtaining the macroscopic physical parameters and mesoscopic physical parameters of the target phanerocrystalline intrusive rock; constructing an initial uniaxial compression model of a rock crystal specimen including mineral composition attributes and crystal structure characteristics according to the mesoscopic physical parameters of the target phanerocrystalline intrusive rock; generating cohesive force units at the contact interfaces between grains in the initial uniaxial compression model of the rock crystal specimen to obtain a target uniaxial compression model of the rock crystal specimen; performing a uniaxial compression simulation test on the target uniaxial compression model of the rock crystal specimen to obtain the test simulation data of the target phanerocrystalline intrusive rock; correcting the mineral grain interface damage and fracture mechanics parameters in the target uniaxial compression model of the rock crystal specimen by an inverse analysis method according to the test simulation data and the target data in the macroscopic physical parameters to obtain the target material parameters of the target phanerocrystalline; constructing a blasting construction FDEM model of a rock crystal medium according to the target material parameters of the target phanerocrystalline intrusive rock; and performing blasting simulation on the target phanerocrystalline intrusive rock according to the blasting construction FDEM model of the rock crystal medium and the blasting parameters.

[0058] In this way, through the comprehensive acquisition and analysis of macroscopic physical parameters and mesoscopic physical parameters, the crystal specimen model constructed by this application can more realistically reflect the internal structural characteristics of the target phanerocrystalline intrusive rock, thereby improving the accuracy of the simulation results. And through the uniaxial compression simulation test and the back-analysis method, the grain boundary damage and fracture mechanics parameters are corrected, so that the finally obtained target material parameters are closer to the mechanical properties of real rocks, providing guarantee for the accuracy of blasting simulation. Constructing the FDEM (Finite Discrete Element Method) model of rock crystal medium can effectively handle the complex coupling problems of rock mass fracture and particle movement, improving the construction efficiency and simulation performance of the blasting construction model. Combining the analysis of blasting parameters during the simulation process can intuitively display the dynamic response law of the rock mass, thereby providing an optimization basis for actual blasting design and reducing the waste of resources and unnecessary cost investment in actual tests. And the FDEM (Finite-Discrete Coupling Method) in this solution realizes the automatic transition of brittle materials such as rocks and concrete from continuous to discontinuous deformation by using the cohesive element model, thus realizing the natural coupling of FEM and DEM in space and time. In addition, through the microscopic structural response analysis at the crystal level, the microcrack evolution, interface failure mechanism and crack penetration law of the rock mass during blasting can be quantitatively described, providing a scientific basis for studying the failure behavior of phanerocrystalline intrusive rocks.

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

[0060] To make the above objects, features and advantages of this application more obvious and understandable, the following specifically gives preferred embodiments and, in conjunction with the accompanying drawings, the detailed description is as follows. Brief Description of the Drawings

[0061] In order to more clearly illustrate the technical solutions of the embodiments of this application, the following will briefly introduce the drawings required to be used in the embodiments. It should be understood that the following drawings only show some embodiments of this application, and therefore should not be regarded as limiting the scope. For those of ordinary skill in the art, without creative efforts, other related drawings can also be obtained based on these drawings.

[0062] Figure 1 It is a flow chart of a blasting simulation method for a phanerocrystalline intrusive rock provided by an embodiment of this application;

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

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

[0065] Figure 4 A diagram of the uniaxial compression test results of rocks provided for this application;

[0066] Figure 5 A diagram of the rock blasting FDEM model provided for this application;

[0067] Figure 6 A schematic diagram of the rock blasting results provided for this application;

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

[0069] Figure 8 A schematic diagram of the structure of an electronic device provided in the embodiments of this application. Detailed implementation manners

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

[0071] In the fields of rock mechanics and blasting engineering, phanerocrystalline intrusive rocks (such as granite, diorite, gabbro, etc.) are well-known for their complex physical and mechanical properties, high hardness, and high compressive strength. These rocks are commonly found in engineering construction, such as mine exploitation, tunnel excavation, and hydropower station foundation construction. However, due to the crystal properties within 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 the mechanical response and failure mode under blasting loads. This poses a huge challenge to blasting construction design and control.

[0072] Specifically, due to the tight binding of grains within phaneritic intrusive rocks, they exhibit extremely high compressive strength, but this also increases the energy requirements during blasting operations. Additionally, due to the distribution and orientation of the crystal structure, phaneritic intrusive rocks exhibit different mechanical properties in different directions, further increasing the complexity of crack propagation during blasting. Moreover, during the blasting process, phaneritic intrusive rocks may experience intergranular failure (fracture at grain boundaries) or transgranular failure (fracture within the crystal). This makes it difficult to quantitatively characterize the fracture mode of the rock under the action of stress waves through simple empirical formulas.

[0073] Currently, research on blasting operations for phaneritic intrusive rocks mainly relies on empirical formulas and on-site tests. However, these methods are not only costly and time-consuming in practical engineering applications but also difficult to quantitatively characterize the dynamic response laws and failure mechanisms of rock masses during blasting.

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

[0075] Currently, there are still many challenges in using continuous-discrete coupling technology to simulate the explosion mechanics behavior and blasting construction process of phaneritic 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 heterogeneity factors such as rock grain size, mineral composition, and content on blasting effects.

[0076] Based on this, the embodiments of this application provide a blasting simulation method and device for phaneritic intrusive rocks. By designing a finite-discrete coupling blasting simulation method, it is not only possible to reduce costs and time consumption but also to quantitatively characterize the dynamic response laws and failure mechanisms of rock masses during blasting.

[0077] Please refer to Figure 1 , Figure 1 which is a flowchart of a blasting simulation method for phaneritic intrusive rocks provided by the embodiments of this application. As Figure 1 shown in

[0078] S101. Obtain the macroscopic physical parameters and mesoscopic physical parameters of the target phaneritic intrusive rock.

[0079] S102. According to the mesoscopic physical parameters of the target phaneritic intrusive rock, construct an initial uniaxial compression model of the rock crystal specimen that includes mineral composition attributes and crystal structure characteristics.

[0080] S103. Generate cohesive units for the contact interfaces between grains in the initial uniaxial compression model of the rock crystal specimen to obtain the target uniaxial compression model of the rock crystal specimen.

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

[0082] S105. According to the test simulation data and the target data in the macroscopic physical parameters, use the inverse analysis method to correct the mineral grain interface damage and fracture mechanics parameters in the target uniaxial compression model of the rock crystal specimen to obtain the target material parameters of the target phaneritic.

[0083] S106. According to the target material parameters of the target phaneritic intrusive rock, construct a FDEM model for rock crystal medium blasting construction.

[0084] S107. According to the FDEM model for rock crystal medium blasting construction and the blasting parameters, conduct a blasting simulation on the target phaneritic intrusive rock.

[0085] Regarding step S101, the phaneritic intrusive rock refers to a rock with an obvious crystal structure formed after magma intrudes into the earth's crust and cools, such as granite, diorite, gabbro.

[0086] The macroscopic physical parameters refer to the mechanical and physical properties at the overall scale of the rock, which are 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 mesoscopic physical parameters refer to the mineral composition, structure, and interface characteristics of the rock at the microscopic scale, which are used to describe the detailed characteristics inside the rock. For example, the mesoscopic physical parameters may include mineral composition, crystal structure, grain size, crystal orientation distribution, etc.

[0088] Among them, the macroscopic physical parameters can be obtained by experiments or existing data to collect the macroscopic characteristics of the phaneritic intrusive rock. The mesoscopic physical parameters can be obtained by means such as microscopes, X-ray diffraction (XRD), etc.

[0089] For step S102, by way of example, this step may include establishing an initial model with mineral grains, crystal boundaries, and mineral distribution based on mesoscopic physical parameters using a numerical simulation tool (such as DEM or FEM), that is, determining the uniaxial compression model of the initial rock crystal specimen.

[0090] Here, the mineral composition attributes may include the mechanical properties and proportional distribution of various minerals (such as quartz, feldspar, mica, etc.) in the rock. The crystal structure characteristics describe the geometric morphology, arrangement pattern, and grain boundary characteristics inside the mineral crystal.

[0091] The uniaxial compression model refers to simulating the mechanical response of a rock under unidirectional loading, and is used to study the compression behavior and fracture mechanism of the rock. <{

[0092] In an implementation provided by the present application, the constructing an uniaxial compression model of the initial rock crystal specimen including mineral composition attributes and crystal structure characteristics according to the mesoscopic physical parameters of the target phaneritic intrusive rock includes:

[0093] S1021. Determine the mineral-related information and crystal characteristic information in the target phaneritic intrusive rock according to the mesoscopic physical parameters of the target phaneritic intrusive rock.

[0094] S1022. According to the crystal characteristic information of the target phaneritic intrusive rock, use the Voronoi algorithm to determine the distribution state of the grains in the model to be constructed.

[0095] S1023. Assign attributes to each grain according to the mineral-related information of the target phaneritic intrusive rock to determine the uniaxial compression model of the initial rock crystal specimen.

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

[0097] By way of example, the crystal structure can be used to determine the lattice structure (such as monoclinic system, orthorhombic system, etc.) of the crystals in the rock through microscopy and crystal diffraction data.

[0098] The grain size can use particle analysis software to statistically analyze the average diameter and distribution range of the grains in the rock (such as 5–15 mm).

[0099] The crystal orientation distribution can analyze the spatial distribution of crystal orientations through a pole figure or a crystal orientation distribution function (ODF).

[0100] For step S1022, the characterization form of each grain in the constructed model is a hexahedral solid element.

[0101] In this step, when using the Voronoi algorithm to establish the uniaxial compression model of the initial rock crystal specimen, for example, it may include: based on the grain size distribution of the target rock, randomly generating seed points, with each point representing the center position of a grain. Using the Voronoi algorithm to divide the space, generating a random grain distribution, and forming a model that conforms to the characteristics of real crystals. And each grain is assigned as a hexahedral solid element, and a boundary region is defined between the grains.

[0102] Here, the grains are characterized by hexahedral solid elements. The grain morphology distribution uses the Voronoi space division algorithm to determine the spatial positions of discrete data points according to the crystal mineral type, the number and size of grains, so as to construct a Delaunay triangulation network to establish a Thiessen polygon polycrystal model that is closest to the actual morphology of rock grains. This method can be determined according to the grain dissection module of the open-source software Neper to directly generate a standard Poisson-Voronoi dissection model and a grain growth model with a wider grading and higher sphericity.

[0103] Regarding step S1023, in this step, when performing attribute assignment, the rules followed include: randomly assigning the mineral types of grains according to the proportion of mineral components, and each mineral grain is assigned different physical properties according to its mechanical parameters (such as elastic modulus, Poisson's ratio).

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

[0105] Regarding step S103, the grain-to-grain contact interface: refers to the boundary region where mineral grains in the rock contact each other, and is the main location for crack initiation.

[0106] Cohesive force unit: a unit used to describe the interface strength of materials in numerical simulation, which can simulate the tensile, shear, and failure behaviors of the interface.

[0107] In an embodiment provided by this application, generating a cohesive force unit for the grain-to-grain contact interface in the uniaxial compression model of the initial rock crystal specimen includes:

[0108] S1031. Performing mesh dissection on the uniaxial compression model of the initial rock crystal specimen to determine the grain-to-grain contact interfaces.

[0109] S1032. Setting corresponding damage criteria and corresponding cohesion-related parameters for a predetermined constitutive model, and loading the model to generate cohesion units at the contact interface between grains.

[0110] Regarding step S1031, in this step, finite element meshing is performed on each grain in the initial rock crystal sample model, the complex grain geometry is discretized into a finite number of mesh units, the meshing 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 surface, the position and area of the contact interface are recorded.

[0112] It should be noted that meshing discretizes complex geometric bodies into simple unit grids (such as hexahedrons and tetrahedrons) to facilitate numerical calculations in finite element analysis. In this application, the meshing is divided into hexahedrons. The contact interface refers to the boundary area between adjacent grains and is an important part of mechanical transmission.

[0113] For step S1032, a constitutive model suitable for phaneritic intrusive rock is selected in this step (such as the traction-separation bilinear constitutive model) to describe the stress-strain relationship of the material under stress. And define the damage criterion, for example, based on shear strength or tensile failure strength, set the interface failure condition (such as the maximum stress or energy release rate exceeds the threshold). In the identified contact interface area, insert a cohesion unit and assign physical parameters to the cohesion unit, including cohesion strength, interface elastic modulus and damage parameters. The stiffness type of the cohesion unit is specified as traction type. For example, it can 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 relationship between stress and strain of a material under external forces and is the basis of rock mechanics analysis.

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

[0116] The cohesion unit is a finite element unit used to simulate the cohesive mechanical behavior on the contact interface and can represent the shear, stretch and failure processes of the interface. In this embodiment, the cohesion unit can be generated at the crystal interface and inside the crystal.

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

[0118] The initial damage criterion can be selected from any one of the following four damage criteria: the maximum nominal strain criterion (Maxe Damage), the maximum nominal stress criterion (MaxsDamage), the quadratic nominal strain criterion (Quads Damage), and the quadratic nominal stress criterion (Quade Damage). The damage evolution criterion can be selected from the damage evolution criteria based on displacement or energy. This solution can select the Maxs initial damage criterion and the damage evolution criterion based on displacement.

[0119] For examples, please refer to Figure 3 , Figure 3 which is a schematic diagram of a uniaxial compression model of a target rock crystal specimen provided for this application. As Figure 3 shown, Figure 3 the grain shown is divided into 306,266 C3D4 solid elements and 111,295 COH3D6 cohesive elements (cohesion elements), and material parameters are assigned to the grain and grain boundary elements representing different minerals respectively.

[0120] Continuing with step S103, in an alternative implementation manner when creating the uniaxial compression model of the target rock crystal specimen, the implementation process includes: First, use Neper software to generate a polycrystalline model (the initial uniaxial compression model of the rock crystal specimen) containing grain boundary cohesion elements; then, use the global cohesion element plug-in to generate an intragranular cohesion element model (cohesion elements). Among them, the strength of the grain boundary cohesion elements is slightly lower than that of the intragranular cohesion elements to conform to the physical and mechanical process that the phaneritic rock first undergoes intergranular fracture and then transgranular fracture under external loads.

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

[0122] Among them, the uniaxial compression simulation test is to simulate the mechanical behavior of the rock under unidirectional compression through numerical simulation means, and observe its failure mode and mechanical response.

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

[0124] In an implementation provided by the present application, performing a uniaxial compression simulation test on the uniaxial compression model of the target rock crystal specimen to obtain test simulation data of the target phaneritic intrusive rock, including: applying an axial load to the uniaxial compression model of the target rock crystal specimen, and analyzing the stress distribution, crack propagation, and failure mode of the rock through finite element simulation to obtain test simulation data.

[0125] Here, a vertical axial load is applied to the top or bottom of the uniaxial compression model of the target rock crystal specimen (usually loaded in a displacement-controlled or force-controlled manner), and the application of the load can be gradually increased (quasi-static loading) or loaded at a certain speed (dynamic loading) to simulate the actual uniaxial compression test conditions. Then, the finite element method (FEM) is used to numerically analyze the loading process to obtain test simulation data.

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

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

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

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

[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 implementation provided by the present application, performing a uniaxial compression simulation test on the uniaxial compression model of the target rock crystal specimen to obtain test simulation data of the target phaneritic intrusive rock, including: defining the material parameters of the solid elements and cohesion elements in the uniaxial compression model of the target rock crystal specimen; where 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 specimen after the definition of the material parameters and performing simulation processing to generate test simulation data.

[0132] Here, first, according to 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. At the same time, determine the material parameters of the cohesion elements (the material parameters include cohesion strength, fracture toughness, interface elastic modulus, etc.); then, apply constraints and loading conditions to the target rock crystal specimen model; finally, use the finite element or explicit dynamics method (such as the discrete element method or the explicit finite difference method) to simulate the compression process, and record relevant data and obtain test simulation data during the simulation process.

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

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

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

[0136] Exemplarily, please refer to Figure 4 , Figure 4 , which is a uniaxial compression test result diagram of rock provided by this application. As Figure 4 shown Figure 4 The test simulation data shown is a stress-strain curve and shows 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 test, and with the laboratory test results as the criterion, the constitutive model parameters of the cohesive unit are corrected backward so that the results of the numerical model are consistent with the laboratory test results.

[0137] For step S105, exemplarily, this step includes: according to the difference between the simulation test data and the actual experimental data, using the inverse analysis method (such as optimization algorithm) to adjust the mineral grain interface parameters (such as strength, cohesion) in the model. And continuously iterate and adjust until the model simulation results are highly consistent with the experimental results to obtain accurate target material parameters.

[0138] Here, when performing the correction, it may specifically be to correct the constitutive model parameters of cohesion. Exemplarily, the method for correcting the constitutive model parameters of cohesion is specifically: First, select initial material parameters according to the literature for calculation to obtain the initial stress-strain relationship curve of rock sample compression; then, compare and analyze the calculation results with the test results, and in combination with the sensitivity degree and influence effect of the material parameters on the results, gradually increase or decrease the parameters step by step; finally, make the calculated stress-strain curve similar in shape and close in amplitude to the test curve. In this example, the simulation results can not only reflect the intergranular fracture during the compression failure process of the specimen, but also reflect the transgranular fracture; in addition, 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 of gradually approaching the experimental results by adjusting the model parameters, and is often used for parameter calibration.

[0140] The target material parameters are the mechanical and physical characteristic parameters used to characterize the target rock after final correction.

[0141] For step S106, by way of example, this step may include: using an FDEM (Finite Discrete Element Method) tool to establish a blasting model containing mineral crystal media in combination with target material parameters. Here, a numerical description of the rock fracture process, such as stress wave propagation, crack propagation and other phenomena, may be introduced into the model.

[0142] Among them, the rock crystal medium is a microscopic structural 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 an implementation manner provided by the present application, constructing an FDEM model for rock crystal medium blasting construction according to the target material parameters of the target phanerocrystalline includes:

[0145] According to the target material parameters of the target phanerocrystalline, respectively create a macroscopic site model and a mesoscopic surrounding rock model adjacent to the borehole; merge adjacent units and nodes in the macroscopic site model and the mesoscopic surrounding rock model to form an FDEM model for rock crystal medium blasting construction.

[0146] This step may include: First, determine the target material parameters of the target phanerocrystalline intrusive rock (by way of example, may include elastic modulus, Poisson's ratio, density, tensile strength, compressive strength, cohesive strength and fracture toughness of the rock, etc.), then select a suitable numerical calculation method (such as finite element, finite difference or boundary element, etc.) for model construction to create the macroscopic site model, and at the same time select a suitable mesoscopic simulation method (such as discrete element method, mesoscopic finite element method, etc.) for model construction to create the mesoscopic surrounding rock model; then, assemble the two into a whole and merge adjacent units and nodes to obtain an FDEM model for rock crystal medium blasting construction.

[0147] It should be noted that the FDEM model for rock crystal medium blasting construction is a macro-mesoscopic coupling model. This model contains both crystal geometric structure information and grain boundaries and Cohesive units inside the crystal generated by Neper software. At the same time, this model can be used to simulate the blasting response of the site, and is particularly suitable for exploring the damage and injury effects of blasting loads on the surrounding rock.

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

[0149] Regarding step S107, in this step, blasting parameters (such as explosive amount, hole spacing, delay time, etc.) are introduced into the FDEM model to simulate the blasting process. And the dynamic response law of the rock mass during blasting is obtained. For example, the results may include crack distribution, fragmentation range, and flying rock trajectory, etc.

[0150] For example, this step may include: setting blasting parameters such as explosive type, charge amount, hole layout, initiation sequence, etc. according to the blasting design requirements; simulating the processes of rock fragmentation, crack propagation, stress wave propagation, and the action of explosive-generated gases under the action of blasting loads. The blasting effect can be evaluated through the analysis of the simulation results, including the range, shape, and distribution characteristics of the crushed zone and the fissure zone, as well as the influence of blasting on the surrounding rocks.

[0151] In this way, the blasting test of phaneritic intrusions can be simulated, and at the same time, the internal stress state of the rock mass, the process of rock fracture and fragmentation can be analyzed.

[0152] For example, please refer to Figure 6 , Figure 6 which is a schematic diagram of the rock blasting result provided by this application. As Figure 6 shown, a is the stress nephogram and b is the fissure distribution diagram. In this way, the blasting test of phaneritic intrusions can be simulated, and at the same time, the internal stress state of the rock mass, the process of rock fracture and fragmentation can be analyzed.

[0153] In this way, through the comprehensive acquisition and analysis of macroscopic physical parameters and mesoscopic physical parameters, the crystal specimen model constructed in this application can more realistically reflect the internal structural characteristics of the target phanerocrystalline intrusive rock, thereby improving the accuracy of simulation results. And through the uniaxial compression simulation test and the back-analysis method, the grain boundary damage and fracture mechanics parameters are corrected, so that the finally obtained target material parameters are closer to the mechanical properties of real rocks, providing guarantee for the accuracy of blasting simulation. Constructing the FDEM (Finite Discrete Element Method) model of rock crystal medium can effectively handle the complex coupling problems of rock mass fracture and particle movement, and improve the construction efficiency and simulation performance of the blasting construction model. Combining the blasting parameters for analysis during the simulation process can intuitively display the dynamic response law of the rock mass, thereby providing an optimization basis for actual blasting design and reducing the waste of resources and unnecessary cost investment in actual tests. In addition, through the microscopic structure response analysis at the crystal level, the microcrack evolution, interface failure mechanism and crack penetration law of the rock mass during the blasting process can be quantitatively described, providing a scientific basis for studying the failure behavior of phanerocrystalline intrusive rocks.

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

[0155] Based on the same inventive concept, the embodiments of this application also provide a blasting simulation device corresponding to the blasting simulation method. Since the principle of solving problems by the device in the embodiments of this application is similar to the above-mentioned blasting simulation method in the embodiments of this application, 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 , Figure 7 which is a schematic structural diagram of a blasting simulation device for phanerocrystalline intrusive rocks provided by the embodiments of this application. As Figure 7 shown in, the blasting simulation device 700 includes:

[0157] An acquisition module 710, configured to acquire the macroscopic physical parameters and mesoscopic physical parameters of the target phanerocrystalline intrusive rock;

[0158] A first construction module 720, configured to construct an initial uniaxial compression model of a rock crystal specimen including mineral composition attributes and crystal structure characteristics according to the mesoscopic physical parameters of the target phanerocrystalline intrusive rock;

[0159] A generation module 730, configured to generate cohesive force units for the contact interfaces between grains in the initial uniaxial compression model of the rock crystal specimen to obtain a target uniaxial compression model of the rock crystal specimen;

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

[0161] The calibration module 750 is used to calibrate the mineral grain interface damage and fracture mechanics parameters in the uniaxial compression model of the target rock crystal specimen by means of back analysis according to the test simulation data and the target data in the macroscopic physical parameters, so as to obtain the target material parameters of the target phaneritic;

[0162] The second construction module 760 is used to construct a FDEM model for rock crystal medium blasting construction according to the target material parameters of the target phaneritic intrusive rock;

[0163] The simulation module 770 is used to conduct a blasting simulation on the target phaneritic intrusive rock according to the FDEM model for rock crystal 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 specimen including mineral composition attributes and crystal structure characteristics according to the mesoscopic physical parameters of the target phaneritic intrusive rock, the first construction module 720 is used for:

[0165] According to the mesoscopic physical parameters of the target phaneritic intrusive rock, determine the mineral-related information and crystal characteristic information in the target phaneritic intrusive rock; wherein, the crystal characteristic information includes: crystal structure, grain size, and crystal orientation distribution;

[0166] According to the crystal characteristic information of the target phaneritic intrusive rock, use the Voronoi algorithm to determine the distribution state of the grains in the model to be constructed; wherein, the characterization form of each grain is a hexahedral solid element;

[0167] According to the mineral-related information of the target phaneritic intrusive rock, assign attributes to each grain to determine the initial uniaxial compression model of the rock crystal specimen.

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

[0169] Perform mesh division on the initial uniaxial compression model of the rock crystal specimen to determine the contact interfaces between different grains;

[0170] Set the corresponding damage criterion and the corresponding cohesive force-related parameters for the constitutive model determined in advance, and perform model loading to realize the generation of cohesive force units at the contact interfaces between the 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 specimen to obtain test simulation data of the target phanerocrystalline intrusive rock, the test module 740 is configured to:

[0172] Apply an axial load to the uniaxial compression model of the target rock crystal specimen, and analyze the stress distribution, crack propagation, and failure mode of the rock through finite element simulation to obtain test 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 specimen to obtain test simulation data of the target phanerocrystalline intrusive rock, the test module 740 is configured to:

[0174] Define the material parameters of the solid elements and cohesion elements in the uniaxial compression model of the target rock crystal specimen; wherein the material parameters of the solid elements are set through macroscopic physical parameters;

[0175] Apply boundary conditions to the uniaxial compression model of the target rock crystal specimen after the material parameters are defined and perform simulation processing to generate test simulation data.

[0176] Optionally, when the second construction module 760 is used to construct a rock crystal medium blasting construction FDEM model according to the target material parameters of the target phanerocrystalline, the second construction module 760 is configured to:

[0177] Create a macroscopic site model and a mesoscopic surrounding rock model adjacent to the borehole according to the target material parameters of the target phanerocrystalline;

[0178] Merge adjacent units and nodes in the macroscopic site model and the mesoscopic surrounding rock model to form a rock crystal medium blasting construction FDEM model.

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

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

[0181] The memory 820 stores machine-readable instructions executable by the processor 810. When the electronic device 800 runs, the processor 810 communicates with the memory 820 through the bus 830. When the machine-readable instructions are executed by the processor 810, they can execute as described above Figures 1 to 6The steps in the method embodiments shown can be specifically implemented by referring to the method embodiments, which will not be elaborated here.

[0182] The embodiments of the present application also provide a computer-readable storage medium, on which a computer program is stored. When the computer program is run by a processor, it can execute the steps in the method embodiments shown above Figures 1 to 6 The steps in the method embodiments shown can be specifically implemented by referring to the method embodiments, which will not be elaborated here.

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

[0184] In the several embodiments provided by the present application, it should be understood that the disclosed systems, devices, and methods can be implemented in other ways. The device embodiments described above are only illustrative. For example, the division of the units is only a logical function division, and there can be other division methods in actual implementation. For another example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the displayed or discussed couplings or direct couplings or communication connections to each other can be through some communication interfaces. The indirect couplings or communication connections of the devices or units can be in electrical, mechanical, or other forms.

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

[0186] In addition, in each embodiment of the present application, the functional units can be integrated into one processing unit, or each unit can exist physically alone, or two or more units can be integrated into one unit.

[0187] When the above-mentioned functions are implemented in the form of software functional units and sold or used as independent products, they can be stored in a non-volatile computer-readable storage medium executable by a processor. Based on such understanding, the technical solution of the present application, in essence, or the part that contributes to the prior art, or a part of this 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 for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods described in various embodiments of the present application. The aforementioned storage medium includes: various media that can store program codes, such as USB flash drives, mobile hard disks, read-only memories (ROM), random access memories (RAM), magnetic disks, or optical discs.

[0188] Finally, it should be noted that the above-mentioned embodiments are only specific implementation manners of the present application, used to illustrate the technical solutions of the present application, rather than limiting them. The protection scope of the present application is not limited thereto. Although the present application has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that: any person skilled in the art within the technical scope disclosed in the present application can still modify the technical solutions recorded in the foregoing embodiments, or can easily think of changes, or perform equivalent replacements on some of the technical features; and these modifications, changes, or replacements do not make the essence of the corresponding technical solutions deviate from the spirit and scope of the technical solutions of the embodiments of the present application, and should all be covered by the protection scope of the present application. Therefore, the protection scope of the present application should be subject to the protection scope of the claims.

Claims

1. A blasting simulation method for phaneritic intrusive rock, characterized in that, The blasting simulation method includes: Obtaining the macroscopic physical parameters and mesoscopic physical parameters of the target phaneritic intrusive rock; Constructing an initial uniaxial compression model of a rock crystal specimen containing mineral composition attributes and crystal structure characteristics according to the mesoscopic physical parameters of the target phaneritic intrusive rock; Generating cohesive force units at the contact interfaces between grains in the initial uniaxial compression model of the rock crystal specimen to obtain a target uniaxial compression model of the rock crystal specimen; Conducting a uniaxial compression simulation test on the target uniaxial compression model of the rock crystal specimen to obtain the test simulation data of the target phaneritic intrusive rock; According to the test simulation data and the target data in the macroscopic physical parameters, correcting the mineral grain interface damage and fracture mechanics parameters in the target uniaxial compression model of the rock crystal specimen through an inverse analysis method to obtain the target material parameters of the target phaneritic; Constructing a blasting construction FDEM model of a rock crystal medium according to the target material parameters of the target phaneritic intrusive rock; Conducting blasting simulation on the target phaneritic intrusive rock according to the blasting construction FDEM model of the rock crystal medium and the blasting parameters.

2. The blasting simulation method according to claim 1, wherein The constructing an initial uniaxial compression model of a rock crystal specimen containing mineral composition attributes and crystal structure characteristics according to the mesoscopic physical parameters of the target phaneritic intrusive rock includes: Determining the mineral-related information and crystal characteristic information in the target phaneritic intrusive rock according to the mesoscopic physical parameters of the target phaneritic intrusive rock; wherein, the crystal characteristic information includes: crystal structure, grain size, and crystal orientation distribution; Determining the distribution state of grains in the model to be constructed by using the Voronoi algorithm according to the crystal characteristic information of the target phaneritic intrusive rock; wherein, the characterization form of each grain is a hexahedral solid element; Assigning attributes to each grain according to the mineral-related information of the target phaneritic intrusive rock to determine the initial uniaxial compression model of the rock crystal specimen.

3. The blasting simulation method according to claim 1, wherein, The generating cohesive force units at the contact interfaces between grains in the initial uniaxial compression model of the rock crystal specimen includes: Performing mesh division on the initial uniaxial compression model of the rock crystal specimen to determine the contact interfaces between different grains; Setting corresponding damage criteria and corresponding cohesive force-related parameters for the constitutive model determined in advance, and performing model loading to realize the generation of cohesive force units at the contact interfaces between grains.

4. The blasting simulation method according to claim 1, wherein The conducting a uniaxial compression simulation test on the target uniaxial compression model of the rock crystal specimen to obtain the test simulation data of the target phaneritic intrusive rock includes: Applying an axial load to the target uniaxial compression model of the rock crystal specimen, and analyzing the stress distribution, crack propagation, and failure mode of the rock through finite element simulation to obtain the test simulation data.

5. The blasting simulation method according to claim 1, characterized in that The conducting a uniaxial compression simulation test on the target uniaxial compression model of the rock crystal specimen to obtain the test simulation data of the target phaneritic intrusive rock includes: Defining the material parameters of the solid elements and cohesive force units in the target uniaxial compression model of the rock crystal specimen; wherein the material parameters of the solid elements are set through macroscopic physical parameters; After defining the material parameters, boundary conditions are applied to the uniaxial compression model of the target rock crystal specimen, and simulation processing is carried out to generate test simulation data.

6. The blasting simulation method according to claim 1, characterized in that Constructing a rock crystal medium blasting construction FDEM model according to the target material parameters of the target phaneritic includes: Creating a macroscopic site model and a mesoscopic surrounding rock model adjacent to the borehole respectively according to the target material parameters of the target phaneritic; Merging adjacent units and nodes in the macroscopic site model and the mesoscopic surrounding rock model to form a rock crystal medium blasting construction FDEM model.

7. The blasting simulation method according to claim 3, wherein The damage criterion includes an initial damage criterion and a damage evolution criterion.

8. A blasting simulation device for phaneritic intrusive rock, characterized in that, The blasting simulation device includes: An acquisition module for acquiring macroscopic physical parameters and mesoscopic physical parameters of the target phaneritic intrusive rock; A first construction module for constructing an initial uniaxial compression model of a rock crystal specimen including mineral composition attributes and crystal structure characteristics according to the mesoscopic physical parameters of the target phaneritic intrusive rock; A generation module for generating cohesive force units at the contact interfaces between grains in the initial uniaxial compression model of the rock crystal specimen to obtain a target uniaxial compression model of the rock crystal specimen; A test module for performing a uniaxial compression simulation test on the target uniaxial compression model of the rock crystal specimen to obtain test simulation data of the target phaneritic intrusive rock; A correction module for correcting the mineral grain interface damage and fracture mechanics parameters in the target uniaxial compression model of the rock crystal specimen by an inverse analysis method according to the test simulation data and the target data in the macroscopic physical parameters to obtain the target material parameters of the target phaneritic; A second construction module for constructing a rock crystal medium blasting construction FDEM model according to the target material parameters of the target phaneritic intrusive rock; A simulation module for performing a blasting simulation on the target phaneritic intrusive rock according to the rock crystal medium blasting construction FDEM model and the blasting parameters.

9. An electronic device, characterized in that, Including: A processor, a memory, and a bus. The memory stores machine-readable instructions executable by the processor. When the electronic device runs, the processor communicates with the memory through the bus. When the machine-readable instructions are run by the processor, the steps of the blasting simulation method according to any one of claims 1 to 7 are executed.

10. A computer-readable storage medium, characterized in that, A computer program is stored on the computer-readable storage medium. When the computer program is run by the processor, the steps of the blasting simulation method according to any one of claims 1 to 7 are executed.

Images