CN107507270A - A kind of generation method for being engaged jointed rock mass periodic structure model - Google Patents
A kind of generation method for being engaged jointed rock mass periodic structure model Download PDFInfo
- Publication number
- CN107507270A CN107507270A CN201710650965.9A CN201710650965A CN107507270A CN 107507270 A CN107507270 A CN 107507270A CN 201710650965 A CN201710650965 A CN 201710650965A CN 107507270 A CN107507270 A CN 107507270A
- Authority
- CN
- China
- Prior art keywords
- rock mass
- joint
- voronoi
- jointed rock
- occlusion
- 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.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2210/00—Indexing scheme for image generation or computer graphics
- G06T2210/32—Image data format
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Computer Graphics (AREA)
- Geometry (AREA)
- Software Systems (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The present invention discloses a kind of generation method for being engaged jointed rock mass periodic structure model, comprises the following steps：The actual occlusion jointed rock mass of field investigation, determines its frequency of joints, simulated domain scope etc.；N number of random seed is generated in simulated domain according to frequency of joints, Voronoi structures are generated using random seed, the structure is random joints strata model；Joint seed position is adjusted, generates the jointed rock mass model being consistent with engineering site joint form；Periodic Voronoi diagram shape is generated by jointed rock mass model；Line displacement is entered to the inner boundary of periodicity Voronoi diagram shape according to frequency of joints, finally gives the occlusion jointed rock mass model with periodic structure.The periodic structure model of generation can be used for studying the properties such as mechanical characteristic, the seepage characteristic of the engineering rock mass of different occlusion joint forms, can be well the research service of related rock mass engineering project, can be studied for engineering safety and provide technical guarantee and theory support.
Description
Technical field
The present invention relates to a kind of generation method of rock mass structure model, and in particular to one kind occlusion jointed rock mass is periodically tied
The generation method of structure model, belong to the works such as rock mechanics engineering, Hydraulic and Hydro-Power Engineering, petroleum works, underground engineering, mining engineering
Journey field.
Background technology
The characteristic of rock mass immanent structure is influence and control rock mass deformation, the primary factor destroyed under engineering effect.No
The geological structure of same time, heterogeneity and different scales, joint inside rock mass is caused to show the construction to come in every shape
Characteristic.Common jointed rock mass structure has complete, fragment, occlusion (interlock), disturbance fragmentation, granular media, stratiform etc..
It is extremely common in engineering to be engaged jointed rock mass structure, it is typically characterised by having periodically.With possessing preferable power
The rockmass difference of performance is learned, occlusion jointed rock mass mechanical property in the case of being undisturbed is preferable, but in strong off-load, the water of constructing
Under the operating modes such as reservoir filling, shock effect, mechanical property often significantly weakens, so as to as engineering construction and safe for operation straight
Connect threat.
In consideration of it, the mechanics parameter of occlusion jointed rock mass is always paid attention to by researcher, certainly, from joint rock
The periodicity of body is set out, and occlusion jointed rock mass structural model of the generation with periodic structure, can study such ask well
Topic.
The content of the invention
Goal of the invention：The present invention is for the occlusion jointed rock mass periodic structure generation of Geotechnical Engineering field in the prior art
Problem, there is provided a kind of generation method for being engaged jointed rock mass periodic structure model, this method are calculated based on Voronoi diagram shape
Method, consider the multifactor occlusion jointed rock mass structural model of the generation with periodic structure such as frequency of joints, form.
Technical scheme：A kind of generation method for being engaged jointed rock mass periodic structure model of the present invention, including such as
Lower step：
(1) field investigation is carried out to actual occlusion jointed rock mass, determines that it is engaged joint distribution form, frequency of joints, section
Width and its characterization unit volume are managed, and simulated domain scope R × R is determined according to characterization unit volume；
(2) N number of random seed is generated in simulated domain R × R according to joint distribution form and frequency of joints, wherein, N's
It is worth for the quantity of actual prismatical joint rock mass in the range of simulated domain；Voronoi structures are generated using random seed, should
Voronoi structures are the jointed rock mass model of completely random distribution, are then based on Lloyd relaxed algorithms adjustment joint seed
Position, generation center Voronoi structures, center Voronoi structures are the jointed rock mass being consistent with engineering site joint form
Model；
(3) by the periodic Voronoi diagram shape of above-mentioned center Voronoi structural generations；
(4) line displacement is entered to the inner boundary of periodicity Voronoi diagram shape according to joint width, finally given with the cycle
The occlusion jointed rock mass model of property structure.
In above-mentioned steps (1), carrying out field investigation to actual occlusion jointed rock mass includes the actual occlusion jointed rock mass of investigation
Distribution, occurrence, density, mark be long and width.
In above-mentioned steps (2), the method based on Lloyd relaxed algorithms adjustment joint seed position comprises the steps：
1. according to live joint distribution feature, the random point of respective numbers is generated；
2. according to the distribution of random point, mesh generation is carried out from each point with respect to center position；
3. calculating the position of each grid element center point, random point is adjusted at this；
4. carrying out mesh generation further according to the new situation of random point, repeat the above steps, until random point position changes
Variate completes the adjustment of joint seed position in allowed band.
In above-mentioned steps (3), included by the method for the periodic Voronoi diagram shape of center Voronoi structural generations following
Step：Single structure cell first is chosen from the Voronoi structures of center, unit born of the same parents are from center Voronoi structure position of form center
The square area outwards diffuseed to form, the 1/9 of Voronoi structures centered on its area；X joint is included in said units born of the same parents
Seed, X joint seed translating into 8, periphery region, quantity of formation is 9X joint seed set in R × Zone R domain,
9X joint seed is subjected to Voronoi subdivisions in R × R space, cutting trimming forms periodic Voronoi diagram shape.
Wherein, cutting trimming refer to carry out 9X joint seed in R × R space Voronoi subdivisions formed it is new
After Voronoi diagram shape, the Voronoi diagram shape wherein with single structure cell corresponding position is intercepted out, is periodic
Voronoi diagram shape.
In above-mentioned steps (4), to the line of each Voronoi polygons in the periodicity Voronoi diagram shape of step (3) generation
Bar enters line displacement, the interval between the lines of adjacent Voronoi polygons is consistent with joint width, that is, is formed and actual column section
The joint plane of rock mass same thickness is managed, finally gives the occlusion jointed rock mass model with periodic structure.
Beneficial effect：Compared with prior art, the advantage of the invention is that：The occlusion jointed rock mass of the present invention is periodically tied
The generation method step of structure model is simple, cost is cheap, practical, passes through the occlusion jointed rock mass that the method for the present invention generates
Periodic structure model has filled up in practical study situations such as generating different cycles structure, different joint thickness, different occlusions
Rock mass structure model blank, can so prepare it is different occlusion jointed rock mass periodic structures numerical value samples, so as to
The shear strength test at different three-dimensional joints is enough carried out, obtains changing rule of the shear strength with different joints；And then carry out
Mechanical characteristic, the seepage characteristic of jointed rock mass are engaged, so as to provide technical guarantee and theory support for engineering safety research.
Brief description of the drawings
Fig. 1 is the flow chart of the generation method of the occlusion jointed rock mass periodic structure model of the present invention；
Fig. 2 is using the random joints strata model of random joints seed generation in embodiment；
Fig. 3 is the procedure chart for adjusting joint seed position in embodiment based on Lloyd relaxed algorithms, wherein, Fig. 3 (a) tables
Show random point generation and subdivision, Fig. 3 (b) represents adjustment random point position and again subdivision, and it is random that Fig. 3 (c) expressions continue adjustment
Point position and subdivision, Fig. 3 (d) represent that the distribution adjustment of Voronoi seeds is completed；
Fig. 4 is that the prismatical joint rock for meeting engineering site joint form generated behind the seed position of joint is adjusted in embodiment
Body Model；
Fig. 5 is after the joint seed in embodiment in unit born of the same parents translates to 8, periphery field, in R × Zone R domain generation 9X
Random joints seed；
Fig. 6 is to carry out Voronoi subdivisions in embodiment in R × Zone R domain, and cutting trimming forms periodicity Voronoi diagram shape；
Fig. 7 is the occlusion jointed rock mass periodic structure model ultimately formed in embodiment.
Embodiment
Below exemplified by generating the structural model of the prismatical joint rock mass at some hydropower station dam foundation, with reference to accompanying drawing to the present invention
Technical scheme be described further.
Such as Fig. 1, a kind of generation method of occlusion jointed rock mass periodic structure model of the invention, comprise the following steps：
(1) form of investigation scene occlusion jointed rock mass, statistics are engaged joint distribution, occurrence, mark length, density, width first
The parameters such as degree, it is determined that occlusion jointed rock mass is uniform jointed rock mass；According to statistics achievement, determine that joint distribution form, joint are close
Degree, joint width and characterization unit volume REV, and determine that simulated domain scope is 15m × 15m according to characterization unit volume REV,
The quantity of actual prismatical joint rock mass is 50 in the regional extent；
(2) joint seed is generated：Generated according to joint distribution form and frequency of joints in simulated domain (15m × 15m)
50 random seeds, Voronoi structures are generated using the random seed of generation, the distribution of different Voronoi seeds can be formed not
With the joint of form, the Voronoi structures of generation are the jointed rock mass model of completely random distribution, such as Fig. 3；
(3) on the basis of the jointed rock mass model of random distribution, the position based on Lloyd relaxed algorithms adjustment joint seed
Put, as Fig. 3 (a)~3 (d), method of adjustment are specific as follows：
1. according to live joint distribution feature, the random point of respective numbers is generated；
2. according to the distribution of random point, mesh generation is carried out from each point with respect to center position；
3. calculating the position of each grid element center point, random point is adjusted at this；
4. carrying out mesh generation further according to the new situation of random point, repeat the above steps, until convergence (i.e. random point
Position change value is in allowed band), that is, complete the adjustment of joint seed position.
Center Voronoi structures such as Fig. 4 generated behind adjustment joint seed position, it can be seen that the jointed rock mass model
Middle joint seed is uniformly distributed, and basis is provided for follow-up selected seed.
(4) periodic Voronoi diagram shape is generated：Chosen in the centre position for the center Voronoi diagram shape adjusted single
Structure 5m × 5m of cellular size, wherein including 5 joint seeds.By 5 joint seeds to the field translation of 8, periphery, in 15m
Generate 9 × 5 random joints seeds in × 15m regions, such as Fig. 5, left side is the seed of single structure cell in figure, 5 Voronoi
Seed is translated according to right side arrow sign in figure, new Voronoi seeds is obtained, by new seed in 15m × 15m regions
Carry out Voronoi subdivisions, produce new Voronoi diagram shape, such as Fig. 6, it is periodic pattern；Intercept in new Voronoi diagram shape
With the Voronoi polygons of unit born of the same parents' position correspondence, periodic Voronoi diagram shape is produced, such as square frame area in centre position in Fig. 6
Figure in domain；
(5) according to the actual joint width for being engaged jointed rock mass, the inner boundary of periodicity Voronoi diagram is carried out partially
Move, the joint plane with actual jointed rock mass same thickness is formed, so as to ultimately form the jointed rock mass mould with periodic structure
Type, such as Fig. 7.
In the jointed rock mass structural model finally given, the occurrence scene of can largely simulating of prismatical joint is true
The feature such as rock mass structure, joint width, density is consistent with field actual measurement results., can by numerical simulation means using the model
To obtain related mechanics parameter and the seepage flow characteristics parameter of prismatical joint rock mass at some hydropower station dam foundation, so as to study column at this
The engineering properties at joint.
Claims (6)
1. a kind of generation method for being engaged jointed rock mass periodic structure model, it is characterised in that comprise the following steps：
(1) field investigation is carried out to actual occlusion jointed rock mass, determines that its occlusion joint distribution form, frequency of joints, joint are wide
Degree and its characterization unit volume, and simulated domain scope R × R is determined according to characterization unit volume；
(2) N number of random seed is generated in simulated domain R × R according to joint distribution form and frequency of joints, wherein, N value is
The quantity of actual prismatical joint rock mass in the range of simulated domain；Voronoi structures, Voronoi knots are generated using random seed
Structure is the jointed rock mass model of completely random distribution, is then based on position, the generation of Lloyd relaxed algorithms adjustment joint seed
Center Voronoi structures, center Voronoi structures are the jointed rock mass model being consistent with engineering site joint form；
(3) by the periodic Voronoi diagram shape of above-mentioned center Voronoi structural generations；
(4) line displacement is entered to the inner boundary of periodicity Voronoi diagram shape according to joint width, finally given with periodically knot
The occlusion jointed rock mass model of structure.
2. the generation method of occlusion jointed rock mass periodic structure model according to claim 1, it is characterised in that step
(1) in, it is described actual occlusion jointed rock mass is carried out field investigation include the distribution of the actual occlusion jointed rock mass of investigation, occurrence,
Density, mark length and width.
3. the generation method of occlusion jointed rock mass periodic structure model according to claim 1, it is characterised in that step
(2) in, the method based on Lloyd relaxed algorithms adjustment joint seed position comprises the steps：
1. according to live joint distribution feature, the random point of respective numbers is generated；
2. according to the distribution of random point, mesh generation is carried out from each point with respect to center position；
3. calculating the position of each grid element center point, random point is adjusted at this；
4. carrying out mesh generation further according to the new situation of random point, repeat the above steps, until random point position change value
In allowed band, that is, complete the adjustment of joint seed position.
4. the generation method of occlusion jointed rock mass periodic structure model according to claim 1, it is characterised in that step
(3) in, comprised the steps by the method for the periodic Voronoi diagram shape of center Voronoi structural generations：First from center
Single structure cell is chosen in Voronoi structures, unit born of the same parents outwards diffuse to form from center Voronoi structures position of form center
Area is 1/9 (R × R) square area, includes X joint seed in the unit born of the same parents, by X joint seed to periphery 8
In individual region translate, in R × Zone R domain quantity of formation be 9X joint seed set, by 9X joint seed R × R sky
Interior progress Voronoi subdivisions, cutting trimming form periodic Voronoi diagram shape.
5. the generation method of occlusion jointed rock mass periodic structure model according to claim 4, it is characterised in that described
After cutting trimming refers to that 9X joint seed is carried out to Voronoi subdivisions in R × R space forms new Voronoi diagram shape,
The Voronoi diagram shape wherein with single structure cell corresponding position is intercepted out, as periodic Voronoi diagram shape.
6. the generation method of occlusion jointed rock mass periodic structure model according to claim 1, it is characterised in that step
(4) in, line displacement is entered to the lines of each Voronoi polygons in the periodicity Voronoi diagram shape of step (3) generation, makes phase
Interval between the lines of adjacent Voronoi polygons is consistent with joint width, that is, is formed and actual prismatical joint rock mass same thickness
Joint plane, finally give the occlusion jointed rock mass model with periodic structure.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710650965.9A CN107507270B (en) | 2017-08-02 | 2017-08-02 | A kind of generation method being engaged jointed rock mass periodic structure model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710650965.9A CN107507270B (en) | 2017-08-02 | 2017-08-02 | A kind of generation method being engaged jointed rock mass periodic structure model |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107507270A true CN107507270A (en) | 2017-12-22 |
| CN107507270B CN107507270B (en) | 2019-09-10 |
Family
ID=60690354
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710650965.9A Active CN107507270B (en) | 2017-08-02 | 2017-08-02 | A kind of generation method being engaged jointed rock mass periodic structure model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107507270B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108372571A (en) * | 2018-02-09 | 2018-08-07 | 河海大学 | A kind of preparation method of prismatical joint rock-like materials sample |
| CN109885867A (en) * | 2019-01-08 | 2019-06-14 | 华能澜沧江水电股份有限公司 | Periodic boundary applying method considering rock mass joint |
| CN110362914A (en) * | 2019-07-08 | 2019-10-22 | 四川农业大学 | A kind of irregular random generation method of prismatical joint network model |
| CN111157322A (en) * | 2020-01-19 | 2020-05-15 | 河海大学 | Preparation method of irregular columnar jointed rock mass sample of nested interlayer dislocation zone |
| CN111428359A (en) * | 2020-03-23 | 2020-07-17 | 河海大学 | Anisotropic rock modeling method considering grain occlusion |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104374619A (en) * | 2014-10-21 | 2015-02-25 | 河海大学 | Preparation method of irregular cylindrical joint crack network model rock core test sample |
| US20160147973A1 (en) * | 2014-11-26 | 2016-05-26 | Jeffrey W. Holcomb | Method for the computation of voronoi diagrams |
| CN106339528A (en) * | 2016-08-09 | 2017-01-18 | 鞍钢集团矿业有限公司 | Method for forecasting range of surface movement induced by underground mining of open-pit iron mine end slope |
| CN106547942A (en) * | 2016-09-26 | 2017-03-29 | 昆明理工大学 | A kind of computational methods of the passive lateral rock pressure of flare bedding joint rock mass |
-
2017
- 2017-08-02 CN CN201710650965.9A patent/CN107507270B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104374619A (en) * | 2014-10-21 | 2015-02-25 | 河海大学 | Preparation method of irregular cylindrical joint crack network model rock core test sample |
| US20160147973A1 (en) * | 2014-11-26 | 2016-05-26 | Jeffrey W. Holcomb | Method for the computation of voronoi diagrams |
| CN106339528A (en) * | 2016-08-09 | 2017-01-18 | 鞍钢集团矿业有限公司 | Method for forecasting range of surface movement induced by underground mining of open-pit iron mine end slope |
| CN106547942A (en) * | 2016-09-26 | 2017-03-29 | 昆明理工大学 | A kind of computational methods of the passive lateral rock pressure of flare bedding joint rock mass |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108372571A (en) * | 2018-02-09 | 2018-08-07 | 河海大学 | A kind of preparation method of prismatical joint rock-like materials sample |
| CN109885867A (en) * | 2019-01-08 | 2019-06-14 | 华能澜沧江水电股份有限公司 | Periodic boundary applying method considering rock mass joint |
| CN109885867B (en) * | 2019-01-08 | 2022-10-25 | 华能澜沧江水电股份有限公司 | Periodic boundary applying method considering rock mass joint |
| CN110362914A (en) * | 2019-07-08 | 2019-10-22 | 四川农业大学 | A kind of irregular random generation method of prismatical joint network model |
| CN110362914B (en) * | 2019-07-08 | 2023-01-10 | 四川农业大学 | Random generation method for irregular columnar joint network model |
| CN111157322A (en) * | 2020-01-19 | 2020-05-15 | 河海大学 | Preparation method of irregular columnar jointed rock mass sample of nested interlayer dislocation zone |
| CN111157322B (en) * | 2020-01-19 | 2021-05-11 | 河海大学 | Preparation method of irregular columnar jointed rock mass sample of nested interlayer dislocation zone |
| CN111428359A (en) * | 2020-03-23 | 2020-07-17 | 河海大学 | Anisotropic rock modeling method considering grain occlusion |
| CN111428359B (en) * | 2020-03-23 | 2022-03-08 | 河海大学 | Anisotropic rock modeling method considering grain occlusion |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107507270B (en) | 2019-09-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107507270B (en) | A kind of generation method being engaged jointed rock mass periodic structure model | |
| Rezanejad et al. | Experimental and numerical investigation of the hydrodynamic performance of an oscillating water column wave energy converter | |
| Bae et al. | Lake Sihwa tidal power plant project | |
| Larentis et al. | Gis-based procedures for hydropower potential spotting | |
| CN103544342A (en) | Nuclear power plant breakwater overtopping impact simulation method based on mixed model | |
| CN104460343A (en) | Riverway flood inundation simulation method based on water level monitoring data | |
| CN110046400B (en) | High-precision numerical simulation calculation method for static characteristic analysis of concrete gravity dam based on h-p type finite element method | |
| CN108108563B (en) | Cross-platform hydrodynamics modeling method based on BIM and GIS technology | |
| CN108090283A (en) | The finite element method of tunnel-vehicle coupled vibrations under a kind of DYNAMIC LOADING OF DRIVING TRAIN ON BRIDGES | |
| Oosterbaan et al. | Numerical thermal back-calculation of the kerava solar village underground thermal energy storage | |
| CN108765573A (en) | A kind of analogy method of underground workshop osculum curtain | |
| CN110287634B (en) | Dam abutment deformation simulation method based on volume force application | |
| Beheshti et al. | Seepage analysis of rock-fill dam subjected to water level fluctuation: a case study on Gotvand-Olya dam | |
| CN111455975A (en) | Karst foundation grouting reinforcement construction method based on BIM technology | |
| Su et al. | Deformation-based safety monitoring model for high slope in hydropower project | |
| CN116383917A (en) | Large dynamic deformation numerical simulation method for dykes and dams in seismic liquefaction field | |
| CN107704646A (en) | Modeling method after compact reservoir volume modification | |
| Chu et al. | A process-based approach to sediment transport in the Yangtze Estuary | |
| Xia et al. | Study on coupling of 3 D visualization with numerical simulation for powerhouse excavation of a certain hydro-junction. | |
| Li et al. | Application of new geological modeling technology in secondary development in Daqing oil field | |
| Xie et al. | Dynamic Model and Numerical Simulation of Maximum Turbidity Zone Formation in River Inlet | |
| Zhang et al. | Large scale 3D seepage analysis of whole pumped storage power station project area | |
| Jiang et al. | Research on critical groundwater level under the threshold value of land subsidence in the typical region of Beijing | |
| Longqiang et al. | Application research of 3D geological modeling based on discrete smooth interpolation | |
| Deng et al. | Identification of hydraulic conductivity in aquifer for coupled FEM and adaptive genetic 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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |