CN117291083A - Rock discrete element model building method, rock mechanics multi-scale computing method, system and medium - Google Patents
Rock discrete element model building method, rock mechanics multi-scale computing method, system and medium Download PDFInfo
- Publication number
- CN117291083A CN117291083A CN202311165711.XA CN202311165711A CN117291083A CN 117291083 A CN117291083 A CN 117291083A CN 202311165711 A CN202311165711 A CN 202311165711A CN 117291083 A CN117291083 A CN 117291083A
- Authority
- CN
- China
- Prior art keywords
- rock
- parameters
- model
- discrete
- mesoscopic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 239000011435 rock Substances 0.000 title claims abstract description 140
- 238000000034 method Methods 0.000 title claims abstract description 69
- 238000004364 calculation method Methods 0.000 title claims abstract description 6
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 56
- 239000011707 mineral Substances 0.000 claims abstract description 56
- 238000012360 testing method Methods 0.000 claims abstract description 33
- 239000013078 crystal Substances 0.000 claims abstract description 21
- 238000004088 simulation Methods 0.000 claims abstract description 17
- 239000002245 particle Substances 0.000 claims description 36
- 230000006835 compression Effects 0.000 claims description 12
- 238000007906 compression Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 12
- 239000004568 cement Substances 0.000 claims description 11
- 238000013401 experimental design Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 8
- 238000013461 design Methods 0.000 claims description 8
- 239000000463 material Substances 0.000 claims description 6
- 238000005457 optimization Methods 0.000 claims description 6
- 241000613130 Tima Species 0.000 claims description 5
- 238000012216 screening Methods 0.000 claims description 4
- 238000002474 experimental method Methods 0.000 claims description 3
- 238000011084 recovery Methods 0.000 claims description 3
- 238000005096 rolling process Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 238000011160 research Methods 0.000 description 5
- 238000005553 drilling Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 241000282414 Homo sapiens Species 0.000 description 2
- 238000003491 array Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000005065 mining Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 201000004569 Blindness Diseases 0.000 description 1
- 238000009411 base construction Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000265 homogenisation Methods 0.000 description 1
- 238000000329 molecular dynamics simulation Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/25—Design optimisation, verification or simulation using particle-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention belongs to the technical field of rock mechanics, and provides a rock discrete element model building method, a rock mechanics multi-scale computing method, a system and a medium, wherein the method comprises the following steps: based on the rock sample, acquiring microscopic mechanical parameters of the rock; establishing an accurate mineral crystal model of the rock, and acquiring stress-strain curves, damage forms and macroscopic mechanical parameters of the rock; establishing a discrete meta model; determining and adjusting mesoscopic parameters, and carrying out a uniaxial simulation test based on a discrete meta-model; and (3) comparing the stress-strain curve and the damage form in the step (S4) with the stress-strain curve and the damage form in the step (S2), if the error of the comparison result is within the allowable range, determining the optimal combination of the mesoscopic parameters, and establishing a rock discrete element model, otherwise, adjusting the mesoscopic parameters and performing a uniaxial simulation test until the error of the comparison result is within the allowable range. Compared with the prior art, the invention has the advantages of lower use cost and ensured reliability of the result.
Description
Technical Field
The invention belongs to the technical field of rock mechanics, and particularly relates to a rock discrete element model building method, a rock mechanics multi-scale computing system and a rock mechanics medium.
Background
The exploration of outer space has important significance for researching the origin of solar systems, mining resource exploitation, space observation stations and the like, and the exploration activities of human beings on the outer space such as planetary scientific research, mining resource exploitation, space observation stations, space base construction and the like can not leave the support of the related theory of planetary geotechnical mechanics and engineering machinery; therefore, the research on the mechanical properties of the planet rocks has important significance for the space construction task to be smoothly carried out.
At present, the research on macroscopic mechanical properties and engineering response of outer space rocks is less, simulation such as outer space drilling, collision and impact is helpful for people to know the properties of the outer space rocks, and discrete element models of the outer space rocks need to be built.
In the prior art, in the process of establishing discrete elements of rock, the mesoscopic parameters of the rock are required to be calibrated in advance, and the process of calibrating the mesoscopic parameters is as follows: and adjusting the calibrated parameters in the discrete element model established by the discrete element to perform macroscopic test simulation (uniaxial compression and the like), comparing the result obtained by the discrete element simulation with the result obtained by the macroscopic mechanical test (uniaxial compression and the like) obtained by a laboratory, and continuously adjusting the calibrated parameters according to the parameter difference between the discrete element model and the result until the macroscopic parameters of the discrete element and the macroscopic mechanical parameters obtained by the laboratory are within the error allowable range, wherein the obtained microscopic parameters can establish the discrete element model of the outer space rock. However, outer space merle is expensive and rare, small in size and random, and difficult to process into rock samples required by traditional macroscopic mechanical tests, and these problems result in the fact that existing discrete meta-model technical solutions are not suitable for space rocks, and certain difficulties are brought to the establishment of discrete meta-models of space rocks.
Disclosure of Invention
The invention aims to provide a rock discrete meta-model building method, which aims to solve the technical problems that outer space meteorites in the prior art are expensive, rare, small in size and random, and are difficult to process into rock samples required by traditional macroscopic mechanical tests, the existing discrete meta-model technical scheme is not suitable for space rocks, and a certain difficulty is brought to building the discrete meta-model of the space rocks.
In order to achieve the above purpose, the technical scheme adopted by the invention is as follows:
a rock discrete element model building method comprises the following steps:
step S1: acquiring micromechanic parameters of the rock based on the rock sample, the micromechanic parameters including young's modulus of various minerals and mineral cements in the rock;
step S2: according to the micromechanics parameters, an accurate mineral crystal model of the rock is established, boundary conditions are added for uniaxial compression, and stress-strain curves, damage forms and macroscopic mechanical parameters of the rock are obtained, wherein the macroscopic mechanical parameters comprise elastic modulus E and Poisson's ratio v;
step S3: the elastic modulus E and the Poisson ratio v are adopted as intrinsic parameters, and an initial discrete element model is established;
step S4: determining and adjusting mesoscopic parameters, and carrying out a uniaxial simulation test based on a discrete element model to obtain a stress-strain curve and a damage form in discrete elements;
step S5: comparing the stress-strain curve and the damage form in the step S4 with the stress-strain curve and the damage form in the step S2, if the error of the comparison result is within the allowable range, executing the step S7, otherwise, executing the step S6;
step S6: repeating the steps S4 to S5 until the error of the comparison result is within the allowable range;
step S7: determining a preferred combination of the mesoscopic parameters;
step S8: and establishing a rock discrete element model according to the mesoscopic parameters determined in the step S7.
In one embodiment, the specific method of step S1 is as follows:
s1.1, according to a rock sample, acquiring a mineral image of the rock through an integrated mineral analyzer TIMA, acquiring a volume porosity image of the rock through a scanning electron microscope SEM, and respectively testing Young modulus of various minerals and mineral cements in the rock sample by utilizing a nanoindentation and an atomic force microscope AFM;
step S1.2, expanding the mineral image and the volume porosity image obtained in the step S1.1 from cm multiplied by cm to m multiplied by m through image cross-scale expansion;
s1.3, capturing a specific size based on the expanded image to obtain a captured image, and establishing an accurate mineral crystal model in the step S2.
In one embodiment, the specific method of step S1.2 is as follows: dividing a cm multiplied by cm picture into a plurality of parts according to mineral grains, then carrying out rotation and mirror image processing on each part, and finally splicing the processed parts into an m multiplied by m picture, wherein the using times of each part in the splicing process are equal.
In one embodiment, the specific method of step S2 is as follows:
s2.1, establishing a two-dimensional rectangular model with the same pixel size as the intercepted picture in finite element software, and dividing square grids with the unit of 1;
s2.2, endowing the rock micromechanics parameters into a two-dimensional rectangular model through material positions corresponding to each pixel of the intercepted picture, and inserting mineral cement between different minerals to obtain an accurate mineral crystal model;
and S2.3, based on an accurate mineral crystal model, adding boundary conditions to perform uniaxial compression to obtain stress-strain curves, damage forms and macroscopic mechanical parameters of the rock.
In one embodiment, the specific method of step S3 is as follows: in discrete meta-software, adopting elastic modulus E and Poisson ratio v as intrinsic parameters, creating a particle material, setting parameters of a particle container, and generating a discrete meta-model in the particle container; the cohesive contact model Hertz-Mindlin was selected during the creation of the discrete metamodel to bind the rock particles as a whole.
In one embodiment, the mesoscopic parameter in step S4 is determined and adjusted by the method of orthogonal experimental design and the method of CCD experimental design, and the specific steps are as follows:
s4.1, determining mesoscopic parameters to be calibrated, wherein the mesoscopic parameters comprise intrinsic parameters, particle contact parameters and particle combination parameters;
s4.2, screening out microscopic parameters which have obvious influence on each macroscopic parameter by adopting an orthogonal test design method;
and S4.3, estimating the nonlinear relation between the test index and the factors by adopting CCD experimental design, and adjusting the microscopic parameters according to the orthogonal experiment and CCD experimental design results.
In one embodiment, the intrinsic parameters in step S4.1 include elastic modulus, poisson' S ratio, which are provided by the macro-mechanical parameters; the particle contact parameters include a collision recovery coefficient, a static friction coefficient, and a rolling friction coefficient, the particle binding parameters include a normal stiffness, a tangential stiffness, a normal strength, a tangential strength, and a contact radius, and the particle contact parameters and the particle binding parameters are determined by referencing a quantitative relationship between macro-micro parameters of the rock.
In order to achieve the above purpose, the invention also provides a rock mechanics multi-scale calculation method, which applies the rock discrete element model established by the rock discrete element model establishment method.
In order to achieve the above object, the present invention further provides a system for implementing the method for building a discrete rock meta-model, comprising:
micromechanics parameter module: acquiring micromechanic parameters of the rock based on the rock sample, the micromechanic parameters including young's modulus of various minerals and mineral cements in the rock;
macroscopic mechanical parameter module: according to the micromechanics parameters, an accurate mineral crystal model of the rock is established, boundary conditions are added for uniaxial compression, and stress-strain curves, damage forms and macroscopic mechanical parameters of the rock are obtained, wherein the macroscopic mechanical parameters comprise elastic modulus E and Poisson's ratio v;
an initial discrete meta-model module: the elastic modulus E and the Poisson ratio v are adopted as intrinsic parameters, and a discrete element model is established;
and a mesoscopic parameter module: determining and adjusting mesoscopic parameters, and carrying out a uniaxial simulation test based on a discrete element model to obtain a stress-strain curve and a damage form in discrete elements;
comparison module: comparing the stress-strain curve and the damage form in the mesoscopic parameter module with the stress-strain curve and the damage form in the macroscopic mechanical parameter module, if the error of the comparison result is within the allowable range, determining the preferred combination of the mesoscopic parameters, otherwise, adjusting the mesoscopic parameters;
and the mesoscopic parameter adjustment module: adjusting the mesoscopic parameters and performing a uniaxial simulation test until the error of the comparison result is within an allowable range;
the mesoscopic parameter optimization module: determining a preferred combination of the mesoscopic parameters;
rock discrete element model module: and establishing a discrete meta-model of the rock according to the mesoscopic parameters determined by the mesoscopic parameter optimization module.
To achieve the above object, the present invention also provides a computer-readable storage medium having stored thereon a computer program to be executed by a processor to implement the rock discrete meta-model building method.
Compared with the prior art, the invention has the following beneficial effects:
(1) According to the invention, SEM, TIMA, AFM, nanoindentation, finite element and discrete element are comprehensively utilized, the rock macro-mechanical parameters are obtained through an accurate mineral crystal model by using the rock micro-mechanical parameters obtained through SEM, TIMA, AFM and nanoindentation, and the obtained rock macro-mechanical parameters are combined with the discrete element to build a discrete element model of the rock; compared with the traditional rock macroscopic mechanical test and discrete element combination for establishing a discrete element model of the rock, the method provided by the invention can acquire the macroscopic mechanical property of the rock based on a small amount of planetary rock, has lower use cost and ensures the reliability of the result;
(2) The method provided by the invention can simulate the problems of drilling, collision, impact and the like of outer space exploration through the discrete element model of the rock established by coupling the accurate mineral crystal model and the discrete element, and provides a method for human to know the mechanical properties of outer space planets.
Drawings
FIG. 1 is a schematic flow chart of the present invention-example 1.
Fig. 2 is a schematic block diagram of embodiment 2 of the present invention.
Detailed Description
The present invention will be further described in detail with reference to examples so as to enable those skilled in the art to more clearly understand and understand the present invention. It should be understood that the following specific embodiments are only for explaining the present invention, and it is convenient to understand that the technical solutions provided by the present invention are not limited to the technical solutions provided by the following embodiments, and the technical solutions provided by the embodiments should not limit the protection scope of the present invention.
It should be noted that the illustrations provided in the following embodiments merely illustrate the basic concept of the present invention by way of illustration, so that only the components related to the present invention are shown in the drawings and are not drawn according to the number, shape and size of the components in actual implementation, the form, number and proportion of the components in actual implementation may be arbitrarily changed, and the layout of the components may be more complicated.
Example 1
As shown in fig. 1, the present embodiment provides a method for building a discrete rock meta-model, and the design principle of the method is as follows: according to the method, the rock discrete element model is built by coupling the accurate mineral crystal model and the discrete elements, and compared with the discrete element model built by combining the traditional rock macroscopic mechanical test and the discrete elements, the method has the advantages that the use cost is low, and meanwhile, the reliability of a result is guaranteed.
In this embodiment, a rock discrete meta-model building method specifically includes the following steps:
1. based on rock sample, obtaining micromechanics parameters of rock
Firstly, according to a rock sample, acquiring a mineral image of the rock through an integrated mineral analyzer TIMA, acquiring a volume porosity image of the rock through a scanning electron microscope SEM, and respectively testing Young modulus of various minerals and mineral cements (intergranular phases) in the rock by utilizing a nanoindentation and an atomic force microscope AFM; the mineral image comprises the microstructure, composition and shape of the rock, and the Young modulus is the micromechanics parameter of the rock sample;
then, expanding the step mineral image and the volumetric porosity image from cm×cm to m×m by image cross-scale expansion; the specific process of expansion is as follows: dividing a cm multiplied by cm picture into a plurality of parts according to mineral grains, then carrying out rotation and mirror image processing on each part, and finally splicing the processed parts into an m multiplied by m picture, wherein the using times of each part in the splicing process are equal.
Finally, capturing a specific size based on the expanded image to obtain a captured picture for establishing an accurate mineral crystal model; wherein, the specific size is a manually set size.
2. Obtaining stress-strain curve, damage form and macroscopic mechanical parameters of rock
Firstly, a two-dimensional rectangular model with the same pixel size as a cut picture is established in finite element software, and square grids with the unit of 1 are divided; secondly, endowing the rock micromechanics parameters into a two-dimensional rectangular model through material positions corresponding to each pixel of the intercepted picture, and inserting mineral cement between different minerals to obtain an accurate mineral crystal model; based on an accurate mineral crystal model, adding boundary conditions to perform uniaxial compression to obtain stress-strain curves, damage forms and macroscopic mechanical parameters of the rock; the precise mineral crystal model can be used for acquiring the macroscopic mechanical parameters of the rock by using methods such as homogenization, molecular dynamics and the like, and the means belong to conventional substitution for a person skilled in the art.
3. Establishing an initial discrete meta-model
In discrete meta-software, adopting elastic modulus E and Poisson ratio v as intrinsic parameters of a discrete meta-model, creating a particle material, setting parameters of a particle container, and generating the discrete meta-model in the particle container; selecting a cohesive contact model Hertz-Mindlin in the process of establishing the discrete element model to bond rock particles into a whole; wherein the discrete meta-software is not particularly limited.
4. Single-axis simulation test is carried out based on discrete element model to obtain stress-strain curve and damage form in discrete element
Determining a mesoscopic parameter by using an orthogonal test design and a CCD (charge coupled device) test design method, and adjusting the mesoscopic parameter to perform a uniaxial simulation test in a discrete element to obtain a stress-strain curve and a damage form in the discrete element; the method for screening the mesoscopic parameters and calibrating the mesoscopic parameters by the orthogonal test design and the CCD test design is as follows:
firstly, determining mesoscopic parameters to be calibrated, wherein the mesoscopic parameters comprise intrinsic parameters, particle contact parameters and particle combination parameters; wherein the intrinsic parameters include elastic modulus, poisson's ratio, which are provided by macroscopic mechanical parameters; the particle contact parameters comprise a collision recovery coefficient, a static friction coefficient and a rolling friction coefficient, the particle combination parameters comprise normal rigidity, tangential rigidity, normal strength, tangential strength and contact radius, and the particle contact parameters and the particle combination parameters are determined by referring to the quantitative relation between macro-micro parameters of the rock; the purpose of confirming the microscopic parameters is to ensure the accuracy of the discrete meta-model of the rock and avoid the blindness of research;
then screening out microscopic parameters which have obvious influence on each macroscopic parameter by adopting an orthogonal test design method; wherein the mesoscopic parameters include one or more of the intrinsic parameters, particle contact parameters, and particle binding parameters described above;
finally, estimating the nonlinear relation between the test index and the factors by adopting a CCD experimental design, and adjusting the mesoscopic parameters according to the orthogonal experiment and the CCD experimental design result; wherein the adjusted mesoscopic parameter is the mesoscopic parameter screened in the previous step.
5. Comparing the stress-strain curve and the damage form in the fourth step with the stress-strain curve and the damage form in the second step
If the error of the comparison result is within the allowable range, executing the seventh step, otherwise, executing the sixth step, wherein the allowable range of the error is set manually, and the person skilled in the art can adjust the error according to the actual research requirement; it should be noted that the model size and boundary conditions during uniaxial compression of discrete elements are the same as those of uniaxial compression of an accurate mineral crystal model.
6. After adjusting the microscopic parameters, performing discrete element simulation until the error of the comparison result is within an allowable range; and repeating the fourth to fifth steps, wherein the adjusted mesoscopic parameters mainly comprise one or more of particle contact parameters and particle combination parameters, after the mesoscopic parameters are adjusted, performing a uniaxial simulation test based on a discrete element model to obtain a stress-strain curve and a damage form in the discrete element, and comparing the stress-strain curve and the damage form in the second step until the error of the comparison result is within an allowable range.
7. Determining preferred combinations of mesoscopic parameters
The preferred combination of the mesoscopic parameters is that the stress-strain curve and the damage form in the fourth step are compared with the stress-strain curve and the damage form in the second step, and when the comparison result is within the allowable error range, the mesoscopic parameters are determined, wherein the mesoscopic parameters comprise one or more of intrinsic parameters, particle contact parameters and particle combination parameters.
8. And establishing a rock discrete meta-model according to the mesoscopic parameters determined in the step seven.
The method can realize the establishment of a final rock discrete element model by coupling an accurate mineral crystal model with discrete elements, and establishes the rock discrete element model, compared with the traditional rock macroscopic mechanical test and discrete element combination to establish the rock discrete element model, the method has the advantages that the use cost is lower, the reliability of results is ensured, and the problems of drilling, collision, impact and the like of outer space exploration can be simulated through the established rock discrete element model.
Example 2
As shown in fig. 2, the present embodiment provides a system for implementing the rock discrete meta-model building method provided in embodiment 1, which specifically includes:
micromechanics parameter module: acquiring micromechanic parameters of the rock based on the rock sample, the micromechanic parameters including young's modulus of various minerals and mineral cements in the rock;
macroscopic mechanical parameter module: according to the micromechanics parameters, an accurate mineral crystal model of the rock is established, boundary conditions are added for uniaxial compression, and stress-strain curves, damage forms and macroscopic mechanical parameters of the rock are obtained, wherein the macroscopic mechanical parameters comprise elastic modulus E and Poisson's ratio v;
an initial discrete meta-model module: the elastic modulus E and the Poisson ratio v are adopted as intrinsic parameters, and a discrete element model is established;
and a mesoscopic parameter module: determining and adjusting mesoscopic parameters, and carrying out a uniaxial simulation test based on a discrete element model to obtain a stress-strain curve and a damage form in discrete elements;
comparison module: comparing the stress-strain curve and the damage form in the mesoscopic parameter module with the stress-strain curve and the damage form in the macroscopic mechanical parameter module, if the error of the comparison result is within the allowable range, determining the preferred combination of the mesoscopic parameters, otherwise, adjusting the mesoscopic parameters;
and the mesoscopic parameter adjustment module: adjusting the mesoscopic parameters and performing a uniaxial simulation test until the error of the comparison result is within an allowable range;
the mesoscopic parameter optimization module: determining a preferred combination of the mesoscopic parameters;
rock discrete element model module: and establishing a discrete meta-model of the rock according to the mesoscopic parameters determined by the mesoscopic parameter optimization module.
It should be noted that the structure and/or principle of each module corresponds to the steps in the rock discrete meta-model building method described in embodiment 1 one by one, so that the description is omitted here.
It should be noted that, it should be understood that the division of each module of the above system is merely a division of a logic function, and may be fully or partially integrated into a physical entity in actual implementation, or may be physically separated, and the modules may be fully implemented in a form of software called by a processing element; or can be realized in hardware; the method can also be realized in a form of calling software by a processing element, and the method can be realized in a form of hardware by a part of modules. For example, a module may be a processing element that is set up separately, may be implemented in a chip of an apparatus, may be stored in a memory of the apparatus in the form of program codes, may be called by a processing element of the apparatus and perform functions of a module, and may be implemented similarly. In addition, all or part of the modules can be integrated together or can be independently implemented. The processing element described herein may be an integrated circuit having signal processing capabilities. In implementation, each step of the above method or each module above may be implemented by an integrated logic circuit of hardware in a processor element or an instruction in a software form.
For example, the modules above may be one or more integrated circuits configured to implement the methods above, such as: one or more specific integrated circuits, or one or more microprocessors, or one or more field programmable gate arrays, etc. For another example, when a module above is implemented in the form of a processing element scheduler code, the processing element may be a general purpose processor, such as a central processing unit or other processor that may invoke the program code. For another example, the modules may be integrated together and implemented in a system-on-chip form.
Example 3
The embodiment provides a rock mechanics multiscale computing method, which is applied to the rock discrete element model established by the rock discrete element model establishing method provided by the embodiment 1, and can simulate the problems of drilling, collision, impact and the like of outer space exploration, thereby providing a method for people to know the mechanical properties of outer space planets.
Example 4
The present embodiment provides a computer-readable storage medium having stored thereon a computer program that is executed by a processor to implement the rock discrete element model building method provided in embodiment 1. Those of ordinary skill in the art will appreciate that: all or part of the steps of implementing the method provided in embodiment 1 may be implemented by hardware associated with a computer program, where the computer program may be stored in a computer readable storage medium, and when executed, the program performs steps including the method provided in embodiment 1; and the storage medium includes: various media that can store program code, such as ROM, RAM, magnetic or optical disks.
Example 5
The embodiment provides a rock discrete meta-model building device, which comprises: a processor and a memory; the memory is used for storing a computer program; the processor is connected to the memory, and is configured to execute the computer program stored in the memory, so that the discrete rock meta-model building device performs the discrete rock meta-model building method provided in embodiment 1.
Specifically, the memory includes: various media capable of storing program codes, such as ROM, RAM, magnetic disk, U-disk, memory card, or optical disk.
Preferably, the processor may be a general-purpose processor, including a central processor, a network processor, etc.; but also digital signal processors, application specific integrated circuits, field programmable gate arrays or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components.
The above embodiments are merely illustrative of the principles of the present invention and its effectiveness, and are not intended to limit the invention. Modifications and variations may be made to the above-described embodiments by those skilled in the art without departing from the spirit and scope of the invention. Accordingly, it is intended that all equivalent modifications and variations of the invention be covered by the claims, which are within the ordinary skill of the art, be within the spirit and scope of the present disclosure.
Claims (10)
1. The rock discrete element model building method is characterized by comprising the following steps of:
step S1: acquiring micromechanic parameters of the rock based on the rock sample, the micromechanic parameters including young's modulus of various minerals and mineral cements in the rock;
step S2: according to the micromechanics parameters, an accurate mineral crystal model of the rock is established, boundary conditions are added for uniaxial compression, and stress-strain curves, damage forms and macroscopic mechanical parameters of the rock are obtained, wherein the macroscopic mechanical parameters comprise elastic modulus E and Poisson's ratio v;
step S3: the elastic modulus E and the Poisson ratio v are adopted as intrinsic parameters, and an initial discrete element model is established;
step S4: determining and adjusting mesoscopic parameters, and carrying out a uniaxial simulation test based on a discrete element model to obtain a stress-strain curve and a damage form in discrete elements;
step S5: comparing the stress-strain curve and the damage form in the step S4 with the stress-strain curve and the damage form in the step S2, if the error of the comparison result is within the allowable range, executing the step S7, otherwise, executing the step S6;
step S6: repeating the steps S4 to S5 until the error of the comparison result is within the allowable range;
step S7: determining a preferred combination of the mesoscopic parameters;
step S8: and establishing a rock discrete element model according to the mesoscopic parameters determined in the step S7.
2. The method for building a discrete rock meta-model according to claim 1, wherein the specific method in the step S1 is as follows:
s1.1, according to a rock sample, acquiring a mineral image of the rock through an integrated mineral analyzer TIMA, acquiring a volume porosity image of the rock through a scanning electron microscope SEM, and respectively testing Young modulus of various minerals and mineral cements in the rock sample by utilizing a nanoindentation and an atomic force microscope AFM;
step S1.2, expanding the mineral image and the volume porosity image obtained in the step S1.1 from cm multiplied by cm to m multiplied by m through image cross-scale expansion;
s1.3, capturing a specific size based on the expanded image to obtain a captured image, and establishing an accurate mineral crystal model in the step S2.
3. The method for building a discrete element model of rock according to claim 2, wherein the specific method of step S1.2 is as follows: dividing a cm multiplied by cm picture into a plurality of parts according to mineral grains, then carrying out rotation and mirror image processing on each part, and finally splicing the processed parts into an m multiplied by m picture, wherein the using times of each part in the splicing process are equal.
4. A method for building a discrete rock meta-model according to claim 3, wherein the specific method in step S2 is as follows:
s2.1, establishing a two-dimensional rectangular model with the same pixel size as the intercepted picture in finite element software, and dividing square grids with the unit of 1;
s2.2, endowing the rock micromechanics parameters into a two-dimensional rectangular model through material positions corresponding to each pixel of the intercepted picture, and inserting mineral cement between different minerals to obtain an accurate mineral crystal model;
and S2.3, based on an accurate mineral crystal model, adding boundary conditions to perform uniaxial compression to obtain stress-strain curves, damage forms and macroscopic mechanical parameters of the rock.
5. The method for building a discrete element model of rock according to claim 4, wherein the specific method of step S3 is as follows: in discrete meta-software, adopting elastic modulus E and Poisson ratio v as intrinsic parameters, creating a particle material, setting parameters of a particle container, and generating a discrete meta-model in the particle container; the cohesive contact model Hertz-Mindlin was selected during the creation of the discrete metamodel to bind the rock particles as a whole.
6. The method for building a discrete rock meta-model according to claim 5, wherein the mesoscopic parameters in the step S4 are determined and adjusted by the method of orthogonal experimental design and the method of CCD experimental design, and the specific steps are as follows:
s4.1, determining mesoscopic parameters to be calibrated, wherein the mesoscopic parameters comprise intrinsic parameters, particle contact parameters and particle combination parameters;
s4.2, screening out microscopic parameters which have obvious influence on each macroscopic parameter by adopting an orthogonal test design method;
and S4.3, estimating the nonlinear relation between the test index and the factors by adopting CCD experimental design, and adjusting the microscopic parameters according to the orthogonal experiment and CCD experimental design results.
7. The rock discrete meta-model building method according to claim 6, wherein the intrinsic parameters in step S4.1 include elastic modulus, poisson' S ratio, which are provided by the macro-mechanical parameters; the particle contact parameters include a collision recovery coefficient, a static friction coefficient, and a rolling friction coefficient, the particle binding parameters include a normal stiffness, a tangential stiffness, a normal strength, a tangential strength, and a contact radius, and the particle contact parameters and the particle binding parameters are determined by referencing a quantitative relationship between macro-micro parameters of the rock.
8. A rock mechanics multiscale computing method, characterized by applying a rock discrete meta model built by the rock discrete meta model building method according to any one of claims 1-7.
9. A system for implementing the rock discrete meta-model building method according to any one of claims 1 to 7, comprising:
micromechanics parameter module: acquiring micromechanic parameters of the rock based on the rock sample, the micromechanic parameters including young's modulus of various minerals and mineral cements in the rock;
macroscopic mechanical parameter module: according to the micromechanics parameters, an accurate mineral crystal model of the rock is established, boundary conditions are added for uniaxial compression, and stress-strain curves, damage forms and macroscopic mechanical parameters of the rock are obtained, wherein the macroscopic mechanical parameters comprise elastic modulus E and Poisson's ratio v;
an initial discrete meta-model module: the elastic modulus E and the Poisson ratio v are adopted as intrinsic parameters, and a discrete element model is established;
and a mesoscopic parameter module: determining and adjusting mesoscopic parameters, and carrying out a uniaxial simulation test based on a discrete element model to obtain a stress-strain curve and a damage form in discrete elements;
comparison module: comparing the stress-strain curve and the damage form in the mesoscopic parameter module with the stress-strain curve and the damage form in the macroscopic mechanical parameter module, if the error of the comparison result is within the allowable range, determining the preferred combination of the mesoscopic parameters, otherwise, adjusting the mesoscopic parameters;
and the mesoscopic parameter adjustment module: adjusting the mesoscopic parameters and performing a uniaxial simulation test until the error of the comparison result is within an allowable range;
the mesoscopic parameter optimization module: determining a preferred combination of the mesoscopic parameters;
rock discrete element model module: and establishing a discrete meta-model of the rock according to the mesoscopic parameters determined by the mesoscopic parameter optimization module.
10. A computer readable storage medium having stored thereon a computer program, wherein the computer program is executed by a processor to implement the rock discrete meta model building method according to any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311165711.XA CN117291083A (en) | 2023-09-11 | 2023-09-11 | Rock discrete element model building method, rock mechanics multi-scale computing method, system and medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311165711.XA CN117291083A (en) | 2023-09-11 | 2023-09-11 | Rock discrete element model building method, rock mechanics multi-scale computing method, system and medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117291083A true CN117291083A (en) | 2023-12-26 |
Family
ID=89238054
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311165711.XA Pending CN117291083A (en) | 2023-09-11 | 2023-09-11 | Rock discrete element model building method, rock mechanics multi-scale computing method, system and medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117291083A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117610395A (en) * | 2024-01-24 | 2024-02-27 | 西安交通大学 | Characterization method, device, equipment and medium for compression hardening memory effect of crystalline rock |
| CN117744412A (en) * | 2024-02-19 | 2024-03-22 | 西安交通大学 | Rock nonlinear mechanical simulation method, system, equipment and medium |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109030202A (en) * | 2018-06-19 | 2018-12-18 | 湘潭大学 | A kind of method of quick determining rock fragile materials discrete element analysis parameter |
| CN110765572A (en) * | 2019-09-12 | 2020-02-07 | 中国科学院武汉岩土力学研究所 | Continuous discontinuous numerical simulation method for single triaxial test of almond-shaped basalt |
| CN111159855A (en) * | 2019-12-12 | 2020-05-15 | 湘潭大学 | Simulation calculation method for crushing sepiolite in stirring mill |
| CN112417715A (en) * | 2020-10-26 | 2021-02-26 | 山东大学 | Rock mass fracture simulation method and system under true triaxial servo loading state |
| CN116338136A (en) * | 2022-10-24 | 2023-06-27 | 武汉大学深圳研究院 | Drilling cuttings mechanical parameter measurement method based on accurate mineral crystal simulation model |
-
2023
- 2023-09-11 CN CN202311165711.XA patent/CN117291083A/en active Pending
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109030202A (en) * | 2018-06-19 | 2018-12-18 | 湘潭大学 | A kind of method of quick determining rock fragile materials discrete element analysis parameter |
| CN110765572A (en) * | 2019-09-12 | 2020-02-07 | 中国科学院武汉岩土力学研究所 | Continuous discontinuous numerical simulation method for single triaxial test of almond-shaped basalt |
| CN111159855A (en) * | 2019-12-12 | 2020-05-15 | 湘潭大学 | Simulation calculation method for crushing sepiolite in stirring mill |
| CN112417715A (en) * | 2020-10-26 | 2021-02-26 | 山东大学 | Rock mass fracture simulation method and system under true triaxial servo loading state |
| CN116338136A (en) * | 2022-10-24 | 2023-06-27 | 武汉大学深圳研究院 | Drilling cuttings mechanical parameter measurement method based on accurate mineral crystal simulation model |
Non-Patent Citations (3)
| Title |
|---|
| SUN, LEI等: "Simulating the entire progressive failure process of rock slopes using the combined finite-discrete element method", 《COMPUTERS & GEOTECHNICS》, vol. 141, 31 January 2022 (2022-01-31), pages 104557 * |
| XUHAI TANG等: "Determining Young\'s modulus of granite using accurate grain-based modeling with microscale rock mechanical experiments", 《INTERNATIONAL JOURNAL OF ROCK MECHANICS AND MINING SCIENCES》, vol. 157, 18 July 2022 (2022-07-18), pages 105167, XP087150315, DOI: 10.1016/j.ijrmms.2022.105167 * |
| 刘曼曼: "考虑粗骨料破碎的混凝土力学特性细观模拟", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》, no. 02, 15 February 2023 (2023-02-15), pages 11 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117610395A (en) * | 2024-01-24 | 2024-02-27 | 西安交通大学 | Characterization method, device, equipment and medium for compression hardening memory effect of crystalline rock |
| CN117610395B (en) * | 2024-01-24 | 2024-04-16 | 西安交通大学 | Characterization method, device, equipment and medium for compression hardening memory effect of crystalline rock |
| CN117744412A (en) * | 2024-02-19 | 2024-03-22 | 西安交通大学 | Rock nonlinear mechanical simulation method, system, equipment and medium |
| CN117744412B (en) * | 2024-02-19 | 2024-05-10 | 西安交通大学 | Rock nonlinear mechanical simulation method, system, equipment and medium |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN117291083A (en) | Rock discrete element model building method, rock mechanics multi-scale computing method, system and medium | |
| Wang et al. | Three-dimensional numerical study on the failure characteristics of intermittent fissures under compressive-shear loads | |
| Nikolic et al. | Rock mechanics model capable of representing initial heterogeneities and full set of 3D failure mechanisms | |
| Mahabadi et al. | Influence of microscale heterogeneity and microstructure on the tensile behavior of crystalline rocks | |
| Tejchman et al. | Shearing of a narrow granular layer with polar quantities | |
| Collop et al. | Use of the distinct element method to model the deformation behavior of an idealized asphalt mixture | |
| Kumar et al. | Macroscopic model with anisotropy based on micro–macro information | |
| Vallejos et al. | Calibration and verification of two bonded-particle models for simulation of intact rock behavior | |
| Bažant et al. | Microplane model M5 with kinematic and static constraints for concrete fracture and anelasticity. I: Theory | |
| Yahya et al. | A unified representation of the plasticity, creep and relaxation behavior of rocksalt | |
| Sun et al. | Simulating the entire progressive failure process of rock slopes using the combined finite-discrete element method | |
| Bai et al. | Numerical investigation of the mechanical and damage behaviors of veined gneiss during true-triaxial stress path loading by simulation of in situ conditions | |
| He et al. | Development of new three-dimensional coal mass strength criterion | |
| Grgic | Constitutive modelling of the elastic–plastic, viscoplastic and damage behaviour of hard porous rocks within the unified theory of inelastic flow | |
| Sugiura et al. | Toward understanding the origin of asteroid geometries-variety in shapes produced by equal-mass impacts | |
| Kosteski et al. | A lattice discrete element method to model the falling-weight impact test of PMMA specimens | |
| CN112131633A (en) | Fluid-solid coupling simulation method and system based on coarse graining calculation theory | |
| Oboudi et al. | Description of inherent and induced anisotropy in granular media with particles of high sphericity | |
| Chen et al. | A spatial decomposition parallel algorithm for a concurrent atomistic-continuum simulator and its preliminary applications | |
| Navas et al. | Optimal transportation meshfree method in geotechnical engineering problems under large deformation regime | |
| You et al. | Three-dimensional microstructural modeling framework for dense-graded asphalt concrete using a coupled viscoelastic, viscoplastic, and viscodamage model | |
| Sett et al. | Soil uncertainty and its influence on simulated G/G max and damping behavior | |
| Balzani et al. | Construction of statistically similar representative volume elements | |
| Piazolo et al. | Numerical simulations of microstructures using the Elle platform: a modern research and teaching tool | |
| Dattola et al. | A distinct element method numerical investigation of compaction processes in highly porous cemented granular materials |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |