CN117669167A - Discrete element numerical simulation method for simulating cutting of rock and soil body in deep sea environment - Google Patents
Discrete element numerical simulation method for simulating cutting of rock and soil body in deep sea environment Download PDFInfo
- Publication number
- CN117669167A CN117669167A CN202311572476.8A CN202311572476A CN117669167A CN 117669167 A CN117669167 A CN 117669167A CN 202311572476 A CN202311572476 A CN 202311572476A CN 117669167 A CN117669167 A CN 117669167A
- Authority
- CN
- China
- Prior art keywords
- rock
- particles
- soil body
- particle
- edge
- 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
- 239000002689 soil Substances 0.000 title claims abstract description 62
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000004088 simulation Methods 0.000 title claims abstract description 26
- 239000011435 rock Substances 0.000 title claims abstract description 7
- 239000002245 particle Substances 0.000 claims abstract description 114
- 230000001133 acceleration Effects 0.000 claims abstract description 8
- 238000003825 pressing Methods 0.000 claims abstract description 7
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 6
- 230000005484 gravity Effects 0.000 claims description 4
- 239000013535 sea water Substances 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012827 research and development Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
A discrete element numerical simulation method for simulating cutting of a rock and soil body in a deep sea environment in the technical field of dredging comprises the following steps: first, a radius r is generated in the simulation area a ～r b As simulated rock-soil mass particles; secondly, the porosity of the rock-soil body particles reaches a target value through scaling operation of the sphere, and certain attribute is given to the rock-soil body particles, so that the rock-soil body particles gradually reach a stress balance state under the action of dead weight; thirdly, identifying particles at the edge of the rock-soil body, and applying pressure perpendicular to the surface of the rock-soil body to the edge particles so as to simulate confining pressure of water on the rock-soil body in a deep sea environment; fourthly, information of the cutting tool is imported, and calculation is carried out through Newton's second law and the relative position and acceleration relation among the tool and the rock-soil body particlesAnd counting the change of the blocked force of the cutting tool with time. The confining pressure is directly applied to soil body edge particles, so that the problem that confining pressure is difficult to solve in cutting numerical simulation in a deep sea environment is solved.
Description
Technical Field
The invention relates to a rock-soil body cutting numerical simulation method in the technical field of dredging, in particular to a discrete element numerical simulation method for simulating rock-soil body cutting in a deep sea environment, which can realize the real-time dynamic display of the topography of a full dredging area and is based on mass conservation and real-time positioning.
Background
With the development of dredging technology, people gradually turn the eyes to deep sea; in order to better exploit deep sea mineral resources, it is often necessary to develop new equipment for deep sea high confining pressure environments. The high confining pressure environment for reproducing deep sea in laboratory is expensive, and researchers put their eyes on numerical simulations. However, it is not easy to accurately reproduce a high confining pressure environment in numerical simulation. How to better simulate the deep sea high confining pressure environment becomes the key of the research and development of the deep sea mining equipment. However, there is no better numerical simulation method in the prior art.
Disclosure of Invention
Aiming at the defects of the prior art, the invention provides a discrete element numerical simulation method for simulating rock-soil body cutting in a deep sea environment, which simulates a deep sea high confining pressure environment by applying vertical stress to rock-soil body edge particles, does not need to be additionally coupled with finite elements, can simulate the rock-soil body cutting and crushing dynamic processes under different water depths, and is convenient to use.
The invention is characterized by the followingThe technical scheme is as follows: the discrete element numerical simulation method for simulating the cutting of the rock and soil body in the deep sea environment is characterized by comprising the following steps of: step one, generating a radius r in a simulation area a ～r b As simulated rock-soil mass particles; secondly, scaling the particles to enable the porosity of the rock-soil particles to reach a target value e, endowing the rock-soil particles with certain elastic modulus, density, friction coefficient and other attributes, and enabling the rock-soil particles to gradually reach a stress balance state under the action of dead weight; step three, recognizing particles at the edge of the rock-soil body, and applying pressure perpendicular to the surface of the rock-soil body to the edge particles so as to simulate confining pressure of water on the rock-soil body in a deep sea environment; and fourthly, importing information such as the position, the speed and the like of the cutting tool, calculating through Newton's second law and the relative position and acceleration relation among the tool, the rock-soil body particles, and counting the change of the blocked force along with time in the movement process of the cutting tool.
In the third step, the edge particles are identified by counting the number of adjacent particles of each rock-soil body particle, and when the number of adjacent particles is lower than a constant n, the rock-soil body particle is considered as the edge particle.
Further, in the third step, the method for judging whether a certain rock-soil body particle is an edge particle is as follows:
step one, respectively identifying coordinate information and radius of each particle center;
step two, calculating the distance between each other particle j and the center of the particle i, wherein the center coordinate information of the particle j is (x) j ，y j ) Radius r j The coordinate information of the center of the particle i is (x i ，y i ) Radius r i The absolute distance between the centers of the particles i and j is
Step three, judging the number n of particles adjacent to the particle i i When d ij ≤1.2*(r i +r j ) The particles are considered to beThe grains i are adjacent;
judging whether the particle i is an edge particle, and when n is i And at less than or equal to e.times.N, the particles are considered to be edge particles. N is a parameter based on dimensions and particle type, n=4 being considered in a general two-dimensional simulation;
wherein r is j Is in units of m, d ij Is in m.
Further, the method of applying pressure to the edge particles perpendicular to the surface of the rock-soil body is as follows:
step one, the coordinate information of the center of the edge particle i is (x i ，y i ) Radius r i Coordinate information (x) identifying adjacent particles of the edge particle i 1 ，y 1 )，(x 2 ，y 2 )，……，(x n ，y n )；
Step two, the x coordinate and the y coordinate of each adjacent particle are averaged,the direction of the pressure is (x) ave -x i ，y ave -y i )；
Step three, calculating the pressure F to which the particles i are subjected i ，F i ＝ρgh•πr i 2 Wherein ρ is the density of seawater in kg/m 3 G is gravity acceleration, g=9.8n/kg; step four, applying pressure to the edge particles i, wherein the pressure is F i The direction is (x ave -x i ，y ave -y i )；
Wherein ρ is the density of seawater in kg/m 3 The method comprises the steps of carrying out a first treatment on the surface of the g is gravity acceleration, g=9.8n/kg; h is the water depth of the position where the particle is located, and the unit is m; f (F) i Is the pressure to which the particles i are subjected, in N.
Compared with the prior art, the invention has the following beneficial effects: the invention is suitable for numerical simulation of discrete elements, can better simulate the effect of a deep sea high confining pressure environment on a rock-soil body, does not need to be coupled with finite elements in the simulation process, and is convenient to use; the method can be used for two-dimensional and three-dimensional numerical simulation and has wide applicability.
Drawings
FIG. 1 is a flow chart of an embodiment of the present invention;
FIG. 2 is a process diagram of numerical simulation of an embodiment of the present invention;
reference numerals in the drawings: 1. a cutting tool 2, edge particles with confining pressure applied, 3, and center particles without confining pressure applied.
Detailed Description
The following describes embodiments of the present invention in detail with reference to the accompanying drawings, and the embodiments and specific operation procedures of the present invention are given by this embodiment on the premise of the technical solution of the present invention, but the protection scope of the present invention is not limited to the following embodiments.
Examples
The invention is shown in the above figures 1 and 2, wherein the first figure is a schematic diagram of an embodiment of the invention, and the second figure is a process diagram of numerical simulation of the invention.
In the implementation process of the invention, the cut soil body is generated by discrete element software, and the soil body radius r a ～r b The method comprises the steps of carrying out a first treatment on the surface of the And (5) giving density to the soil body. The soil body is compacted by only using a radius expansion method (namely, the radius of the soil body is firstly uniformly multiplied by a coefficient smaller than 1 to be contracted until unbalance force among the soil bodies is relatively small, then the radius is divided by the coefficient to be expanded back to the original size), in the process, a certain attribute such as elastic modulus, density, friction coefficient and the like is given to rock-soil body particles, unbalance force and porosity of the soil body are continuously monitored, and after the unbalance force in the soil body is smaller than a certain characteristic value, the soil body is considered to be compacted and stable, and the subsequent simulation research can be carried out.
And identifying the edge particles of the soil body. Step one, respectively identifying coordinate information and radius of each particle center; step two, calculating the distance between each other particle j and the center of the particle i, wherein the center coordinate information of the particle j is (x) j ，y j ) Radius r j The coordinate information of the center of the particle i is (x i ，y i ) Radius r i The absolute distance between the centers of the particles i and j isStep three, judging the number n of particles adjacent to the particle i i When d ij ≤1.2*(r i +r j ) The particle is considered to be adjacent to particle i; judging whether the particle i is an edge particle, and when n is i At +.e.4, the particles are considered to be edge particles.
After identifying the edge particles, confining pressure is applied to the edge particles. Step one, the coordinate information of the center of the edge particle i is (x i ，y i ) Radius r i Coordinate information (x) identifying adjacent particles of the edge particle i 1 ，y 1 )，(x 2 ，y 2 )，……，(x n ，y n ) The method comprises the steps of carrying out a first treatment on the surface of the Step two, the x coordinate and the y coordinate of each adjacent particle are averaged,the direction of the pressure is (x) ave -x i ，y ave -y i ) The method comprises the steps of carrying out a first treatment on the surface of the Step three, calculating the pressure F to which the particles i are subjected i ，F i ＝ρgh•πr i 2 Wherein ρ is the density of seawater in kg/m 3 G is gravity acceleration, g=9.8n/kg; step four, applying pressure to the edge particles i, wherein the pressure is F i The direction is (x ave -x i ，y ave -y i )。
Finally, information such as the position, the speed and the like of the cutting tool is imported, the relative position and the acceleration relation among the tool, the rock-soil body particles are calculated through Newton's second law, and the change of the blocked force with time in the movement process of the cutting tool is counted.
It should be stated that when the cutter cuts the soil body and the soil body breaks, the edge particles are changed, so that the steps of identifying the edge particles and applying confining pressure need to be continuously performed in the calculation.
The foregoing describes a specific mode of operation of the present invention. It is to be understood that the invention is not limited to the particular manner of operation described hereinabove, and that various changes and modifications may be made by one skilled in the art within the scope of the appended claims without affecting the spirit of the invention.
Claims (4)
1. The discrete element numerical simulation method for simulating the cutting of the rock and soil body in the deep sea environment is characterized by comprising the following steps of:
step one, generating a radius r in a simulation area a ～r b As simulated rock-soil mass particles;
secondly, scaling the particles to enable the porosity of the rock-soil particles to reach a target value e, endowing the rock-soil particles with certain properties, and enabling the rock-soil particles to gradually reach a stress balance state under the action of dead weight;
step three, recognizing particles at the edge of the rock-soil body, and applying pressure perpendicular to the surface of the rock-soil body to the edge particles so as to simulate confining pressure of water on the rock-soil body in a deep sea environment;
and fourthly, importing cutting tool information, calculating through Newton's second law and the relative position and acceleration relation among the tools and the rock-soil body particles, and counting the change of the blocked force along with time in the movement process of the cutting tools.
In the third step, the edge particles are identified by counting the number of adjacent particles of each rock-soil body particle, and when the number of adjacent particles is lower than a constant n, the rock-soil body particle is considered as the edge particle.
2. A discrete element numerical simulation method for simulating cutting of a rock and soil body in a deep sea environment according to claim 1, wherein in said step two, properties imparted to the rock and soil body particles include, but are not limited to, modulus of elasticity, density, coefficient of friction; in the fourth step, the information of the introduced cutting tool includes, but is not limited to, position and speed.
3. The discrete element numerical simulation method for simulating cutting of a rock-soil body in a deep sea environment according to claim 1, wherein in the third step, the method for judging whether a certain rock-soil body particle is an edge particle is as follows:
step one, respectively identifying coordinate information and radius of each particle center;
step two, calculating the distance between each other particle j and the center of the particle i, wherein the center coordinate information of the particle j is (x) j ，y j ) Radius r j The coordinate information of the center of the particle i is (x i ，y i ) Radius r i The absolute distance between the centers of the particles i and j is
Step three, judging the number n of particles adjacent to the particle i i When d ij ≤1.2*(r i +r j ) The particle is considered to be adjacent to particle i;
judging whether the particle i is an edge particle, and when n is i And at less than or equal to e.times.N, the particles are considered to be edge particles. N is a parameter based on dimensions and particle type, n=4 being considered in a general two-dimensional simulation;
wherein r is j Is in units of m, d ij Is in m.
4. A discrete element numerical simulation method for simulating cutting of a rock-soil body in a deep sea environment according to claim 1, wherein in the third step, a method of applying a pressure perpendicular to the surface of the rock-soil body to the edge particles is as follows
Step one, the coordinate information of the center of the edge particle i is (x i ，y i ) Radius r i Coordinate information (x) identifying adjacent particles of the edge particle i 1 ，y 1 )，(x 2 ，y 2 )，……，(x n ，y n )；
Step two, the x coordinate and the y coordinate of each adjacent particle are averaged,pressure is thenIs in the direction (x) ave -x i ，y ave -y i )；
Step three, calculating the pressure F to which the particles i are subjected i ，F i ＝ρgh•πr i 2 ；
Step four, applying pressure to the edge particles i, wherein the pressure is F i The direction is (x ave -x i ，y ave -y i )；
Wherein ρ is the density of seawater in kg/m 3 The method comprises the steps of carrying out a first treatment on the surface of the g is gravity acceleration, g=9.8n/kg; h is the water depth of the position where the particle is located, and the unit is m; f (F) i Is the pressure to which the particles i are subjected, in N.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311572476.8A CN117669167A (en) | 2023-11-23 | 2023-11-23 | Discrete element numerical simulation method for simulating cutting of rock and soil body in deep sea environment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311572476.8A CN117669167A (en) | 2023-11-23 | 2023-11-23 | Discrete element numerical simulation method for simulating cutting of rock and soil body in deep sea environment |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117669167A true CN117669167A (en) | 2024-03-08 |
Family
ID=90083679
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311572476.8A Pending CN117669167A (en) | 2023-11-23 | 2023-11-23 | Discrete element numerical simulation method for simulating cutting of rock and soil body in deep sea environment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117669167A (en) |
-
2023
- 2023-11-23 CN CN202311572476.8A patent/CN117669167A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Jiang et al. | New method for characterizing the shear damage of natural rock joint based on 3D engraving and 3D scanning | |
| Di et al. | Macro-mechanical properties of columnar jointed basaltic rock masses | |
| De Barros et al. | Response of a layered viscoelastic half-space to a moving point load | |
| Kulatilake et al. | Estimation of rock mass strength and deformability in 3-D for a 30 m cube at a depth of 485 m at Äspö hard rock laboratory | |
| CN105512400A (en) | Cutting process simulation process for brittle materials | |
| CN114297864B (en) | Cracked loose rock mass slope stability analysis method controlled by steep and gentle dip angles | |
| Forcellini et al. | Response site analyses of 3D homogeneous soil models | |
| CN101009024A (en) | Method for implementing landslide disaster visualization | |
| CN110632653A (en) | Salt structure analysis method, device and system | |
| CN114818427B (en) | Discrete element breakable particle model modeling method based on real particle shape | |
| Barbero et al. | Numerical simulations of compressive tests on bimrock | |
| CN110826276A (en) | Three-dimensional random seepage field simulation method and device based on finite elements | |
| CN115661388A (en) | Three-dimensional fracture network rock mass model characterization method with multi-scale heterogeneity | |
| CN114970235B (en) | Geomechanical heterogeneity-anisotropy modeling method for fractured reservoir | |
| CN113033060B (en) | Optimization method for predicting complex coal seam mining structure | |
| Hong et al. | Comparative study of wave impact loads on circular cylinder by breaking waves | |
| CN117669167A (en) | Discrete element numerical simulation method for simulating cutting of rock and soil body in deep sea environment | |
| Fu et al. | 3D rock mass geometrical modeling with arbitrary discontinuities | |
| JP2023123662A (en) | Liquefaction evaluation model generation device, liquefaction evaluation device, liquefaction evaluation model generation program, liquefaction evaluation program, liquefaction evaluation model generation method, and liquefaction evaluation method | |
| CN108090030B (en) | Processing method of circular single-pile local flow field | |
| Guerra et al. | Multiscale modelling of the slope stability of block-in-matrix materials | |
| CN109033566B (en) | Vehicle load numerical simulation method based on finite element-infinite element model | |
| Lu et al. | Simulating fractures with bonded discrete element method | |
| CN113268874B (en) | LBM-DDA coupling-based accumulation body seepage erosion damage simulation calculation method | |
| CN111080790B (en) | Concave body generation method and device based on convex hull algorithm |
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 |