CN107036561A - The anisotropic approximate expression method of structural plane roughness based on middle intelligence number function - Google Patents
The anisotropic approximate expression method of structural plane roughness based on middle intelligence number function Download PDFInfo
- Publication number
- CN107036561A CN107036561A CN201610998576.0A CN201610998576A CN107036561A CN 107036561 A CN107036561 A CN 107036561A CN 201610998576 A CN201610998576 A CN 201610998576A CN 107036561 A CN107036561 A CN 107036561A
- Authority
- CN
- China
- Prior art keywords
- structural plane
- coefficient
- plane roughness
- roughness coefficient
- directions
- 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
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/30—Measuring arrangements characterised by the use of optical techniques for measuring roughness or irregularity of surfaces
- G01B11/303—Measuring arrangements characterised by the use of optical techniques for measuring roughness or irregularity of surfaces using photoelectric detection means
-
- 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
Abstract
A kind of anisotropic approximate expression method of structural plane roughness based on middle intelligence number function, obtains structural plane three-dimensional point cloud coordinate data；Different dimensional structure facial contour lines and its coordinate data, each structural plane roughness coefficient for surveying section of measurement are uniformly extracted in different measurement directions；Under different size conditions, statistical analysis is carried out respectively to the structural plane roughness coefficient of all directions, the assembly average and standard deviation for obtaining structural plane roughness coefficient is calculated；Assembly average and standard deviation sum are regard as lower limit as higher limit, difference；The anisotropy figure under different size conditions is drawn respectively；Different directions higher limit and lower limit are fitted using elliptic equation, polar middle intelligence number function representation of structural plane roughness coefficient is obtained；Calculate and obtain interval range under different angle conditions；Obtain the interval range of structural plane all directions roughness value under different size conditions.The energy anisotropic uncertainty of approximate expression structural plane roughness of the invention.
Description
Technical field
The invention belongs to field of engineering technology, it is related to a kind of anisotropic middle intelligence number function representation side of structural plane roughness
Method, it is specifically of the invention from the anisotropic essence of structural plane roughness coefficient, do not know to ask with reference to the expression of middle intelligence number
The advantage of topic, it is proposed that the middle intelligence number function representation method of different directions structural plane roughness.
Background technology
The mechanical property of rock mass discontinuity depends primarily on the coarse relief feature on structural plane surface, how accurate comprehensive
The roughness properties in description scheme face are the element tasks of structural face shear strength research.But, due to structure surface development formation ground
The complexity of matter condition, it has anisotropism, dimensional effect and anisotropic build-in attribute, and this significantly increases structural plane
The difficulty of the quantitative assessment of coarse property.Therefore, further investigation rock structural plane roughness anisotropy and dimensional effect are ground
Study carefully and have great importance.
2007, Du Shigui etc. disclosed a kind of obtaining value method for anisotropic rock mass structural plane shearing, was survey with 10 °
Amount interval, the measurement result according to structural plane roughness coefficient on different directions has obtained structural plane anisotropy under different sizes
Figure, research finds that structural plane roughness anisotropy rule is not fully identical under different size conditions, but exists each other certain
Contact.
2007, Du Shigui disclosed a kind of potential glide direction shearing strength of rock mass discontinuity and determines method, the invention
In using 10cm as survey segment length, it is proposed that structural plane JRC anisotropy figure determine method.
2014, Chen Shijiang, Guo Lingfei disclosed a kind of rock structural face pattern anisotropy evaluation method, by obtaining
Rock structural face three-dimensional appearance data, according to experimental variations function formula and dimensional analysis, become journey with parameter and built with base station
Characterize the parameter of structural plane anisotropic character.
2016, Du Shigui etc. disclosed a kind of structural plane roughness coefficient anisotropy evaluation method, specifically should
Invention carries out level than changing to the roughness value of structural plane all directions, and level can be used to represent that structural plane is thick than the coefficient of conversion
The degree of irregularity of roughness coefficient all directions distribution, and then it is anisotropic quantitative to realize rock structural plane roughness coefficient
Evaluate.
However, not accounting for the Statistical Distribution of the structural plane coefficient of roughness in all directions in the studies above, greatly
It all only have studied the anisotropic properties of single size structural plane roughness, it is impossible to expression structure surface roughness index anisotropy
The uncertain regularity of distribution, the nondeterministic function equation of all directions structural plane roughness is not set up.
The content of the invention
In order to overcome in the prior art without appropriate method expression structure surface roughness it is anisotropic it is probabilistic not
Foot, the present invention provides intelligence number function in utilization, it is proposed that a kind of probabilistic approximate expression of structural plane roughness anisotropy
Method.
The technical solution adopted for the present invention to solve the technical problems is：
A kind of anisotropic approximate expression method of structural plane roughness based on middle intelligence number function, comprises the following steps：
(1) large-scale structure interview sample is scanned using three-dimensional laser instrument, obtains structural plane three-dimensional point cloud number of coordinates
According to；
(2) three dimensional point cloud according to coarse scale structures face, different dimensional structures are uniformly extracted in different measurement directions
Facial contour line and its coordinate data, measure each structural plane roughness coefficient for surveying section accordingly；
(3) under different size conditions, statistical analysis is carried out respectively to the structural plane roughness coefficient of all directions, calculated
Obtain the assembly average of structural plane roughness coefficientWith standard deviation sigma；
(4) assembly average of structural plane rugosity coefficient and standard deviation sum are analyzed as structural plane roughness coefficient
Higher limit JRCup；The difference of the assembly average of structural plane rugosity coefficient and standard deviation is analyzed as structural plane roughness coefficient
Lower limit JRCdown；
(5) under different size conditions, the anisotropy of difference rendering architecture surface roughness coefficient higher limit and lower limit
Figure；Different directions higher limit and lower limit are fitted using elliptic equation, and obtain the match value of elliptic equation；
(6) according to the anisotropic upper limit match value of structural plane roughness coefficient and lower limit match value, structural plane is obtained
Polar middle intelligence number function representation of roughness value；
(7) the middle intelligence number function representation according to structural plane roughness coefficient, calculates and obtains structural plane under different angle conditions
The interval range of roughness value；
(8) according to step (3)~(7), the interval of structural plane all directions roughness value under different size conditions is obtained
Scope.
Further, in the step (5), each point in elliptic equation interval range is expressed according to polar form：
Wherein, n1For the major axis of the elliptic equation fitting of structural plane roughness coefficient lower limit under a certain size；n2For a certain chi
The short axle of the elliptic equation fitting of very little lower structural plane roughness coefficient lower limit；For structural plane roughness coefficient under a certain size
2 times of the difference of the major axis of the elliptic equation fitting of upper and lower limit；For structural plane roughness coefficient upper and lower limit under a certain size
2 times of difference of short axle of elliptic equation fitting, I is indeterminacy section scope [0,0.5]；θ is measurement direction.
In the step (6), middle intelligence number function representation is：
Beneficial effects of the present invention are mainly manifested in：Approximate expression structural plane roughness anisotropy is uncertain.
Brief description of the drawings
Fig. 1 is structural plane Point Cloud Data from Three Dimension Laser Scanning and cross-sectional data extraction figure.
Fig. 2 is the upper and lower limitation of 10cm structural plane roughness coefficients and elliptic equation fitting result schematic diagram.
Embodiment
The invention will be further described below in conjunction with the accompanying drawings.
Referring to Figures 1 and 2, a kind of anisotropic approximate expression method of structural plane roughness based on middle intelligence number function,
Comprise the following steps：
(1) large-scale structure interview sample is scanned using three-dimensional laser instrument, obtains structural plane three-dimensional point cloud number of coordinates
According to；
(2) three dimensional point cloud according to coarse scale structures face, different dimensional structures are uniformly extracted in different measurement directions
Facial contour line and its coordinate data, measure each structural plane roughness coefficient for surveying section accordingly；
(3) under different size conditions, statistical analysis is carried out respectively to the structural plane roughness coefficient of all directions, calculated
Obtain the assembly average of structural plane roughness coefficientWith standard deviation sigma；
(4) assembly average of structural plane rugosity coefficient and standard deviation sum are analyzed as structural plane roughness coefficient
Higher limit JRCup；The difference of the assembly average of structural plane rugosity coefficient and standard deviation is analyzed as structural plane roughness coefficient
Lower limit JRCdown；
(5) under different size conditions, the anisotropy of difference rendering architecture surface roughness coefficient higher limit and lower limit
Figure；Different directions higher limit and lower limit are fitted using elliptic equation, and obtain the match value of elliptic equation；
(6) according to the anisotropic upper limit match value of structural plane roughness coefficient and lower limit match value, structural plane is obtained
Polar middle intelligence number function representation of roughness value；
(7) the middle intelligence number function representation according to structural plane roughness coefficient, calculates and obtains structural plane under different angle conditions
The interval range of roughness value；
(8) according to step (3)~(7), the interval of structural plane all directions roughness value under different size conditions is obtained
Scope.
Example：A kind of probabilistic approximate expression method of structural plane roughness anisotropy, comprises the following steps：
(1) field condition selectes Changshan County, Zhejiang Province slate structural plane (1m × 1m), uses high-precision three-dimensional laser scanning
Instrument is scanned, and obtains complicated topography three-dimensional coordinate data, and its cloud data is as shown in Figure 1；
(2) three dimensional point cloud according to coarse scale structures face, different dimensional structures are uniformly extracted in different measurement directions
Illustrated in facial contour line and its coordinate data, Fig. 1 measurement direction for 0 ° with 180 ° when, the arrangement of Extracting contour, survey
Measure the structural plane roughness coefficient that structural plane contour line under each size condition surveys section；
(3) under different size conditions, statistical analysis is carried out respectively to the structural plane roughness coefficient of all directions, calculated
Obtain the assembly average of structural plane roughness coefficientWith standard deviation sigma；
(4) assembly average of structural plane rugosity coefficient and standard deviation sum are analyzed as structural plane roughness coefficient
Higher limit；The difference of the assembly average of structural plane rugosity coefficient and standard deviation is analyzed down as structural plane roughness coefficient
Limit value；Upper limit value and lower limit value of the size for the structural plane roughness coefficient JRC in 10cm all directions is illustrated in Fig. 2；
(5) the anisotropy figure of structural plane roughness coefficient higher limit and lower limit under different size conditions is drawn；Using
Elliptic equation is fitted to different directions higher limit and lower limit, and obtains the match value of elliptic equation；Size is in Fig. 2
The major and minor axis of 10cm structural plane roughness coefficients JRC higher limit elliptic equation fitting result is respectively 13.07 and 9.92, and under
The major and minor axis of limit value elliptic equation fitting result is respectively 8.65 and 5.81.Table 1 is structural plane roughness under the conditions of different scale
The fitting result of upper and lower limit.
Table 1
(6) according to the anisotropic upper limit match value of structural plane roughness coefficient and lower limit match value, by elliptic equation
Each point in interval range is expressed according to polar form：
Wherein, n1For the major axis of the elliptic equation fitting of structural plane roughness coefficient lower limit under a certain size；n2For a certain chi
The short axle of the elliptic equation fitting of very little lower structural plane roughness coefficient lower limit；For structural plane roughness coefficient under a certain size
2 times of the difference of the major axis of the elliptic equation fitting of upper and lower limit；For structural plane roughness coefficient upper and lower limit under a certain size
2 times of difference of short axle of elliptic equation fitting, I is indeterminacy section scope [0,0.5]；θ is measurement direction.
And then obtain polar middle intelligence number function representation of structural plane roughness coefficient：
(7) by taking 10cm structural fece samples as an example, n1=8.65,n2=5.81, μn2=8.23, bring measurement angle into
Degree, you can when calculating the interval range, such as 0 ° of direction of structural plane roughness coefficient under the different angle conditions of acquisition is most of thick
Roughness coefficient is in interval range [8.66,13.07]；During 45 ° of directions, most of roughness value be in interval range [7.37,
11.61]。
(8) according to step (3)~(7), the interval of structural plane all directions roughness value under different size conditions is obtained
Scope.Table 2, table 3 are the forecast interval scope of structural plane roughness coefficient under the conditions of different sizes, different directions.Table 2 is
The forecast interval scope of 10cm-50cm structural plane roughness coefficients.
Table 2
Table 3 is the forecast interval scope of 60cm-100cm structural plane roughness coefficients.
Table 3.
Claims (2)
1. a kind of anisotropic approximate expression method of structural plane roughness based on middle intelligence number function, it is characterised in that：It is described
Approximate expression method comprises the following steps：
(1) large-scale structure interview sample is scanned using three-dimensional laser instrument, obtains structural plane three-dimensional point cloud coordinate data；
(2) three dimensional point cloud according to coarse scale structures face, different dimensional structure face wheels are uniformly extracted in different measurement directions
Profile and its coordinate data, measure each structural plane roughness coefficient for surveying section accordingly；
(3) under different size conditions, statistical analysis is carried out respectively to the structural plane roughness coefficient of all directions, calculating is obtained
The assembly average of structural plane roughness coefficientWith standard deviation sigma；
(4) upper limit for analyzing the assembly average of structural plane rugosity coefficient and standard deviation sum as structural plane roughness coefficient
Value JRCup；The difference of the assembly average of structural plane rugosity coefficient and standard deviation is analyzed down as structural plane roughness coefficient
Limit value JRCdown；
(5) under different size conditions, the anisotropy figure of difference rendering architecture surface roughness coefficient higher limit and lower limit；Adopt
Different directions higher limit and lower limit are fitted with elliptic equation, and obtain the match value of elliptic equation；
(6) according to the anisotropic upper limit match value of structural plane roughness coefficient and lower limit match value, structural plane is obtained coarse
Spend polar middle intelligence number function representation of coefficient；
(7) the middle intelligence number function representation according to structural plane roughness coefficient, calculates structural plane under the different angle conditions of acquisition coarse
Spend the interval range of coefficient；
(8) according to step (3)~(7), the interval range of structural plane all directions roughness value under different size conditions is obtained.
2. the structural plane roughness anisotropic approximate expression method as claimed in claim 1 based on middle intelligence number function, its
It is characterised by：In the step (5), each point in elliptic equation interval range is expressed according to polar form：
Wherein, n1For the major axis of the elliptic equation fitting of structural plane roughness coefficient lower limit under a certain size；n2For under a certain size
The short axle of the elliptic equation fitting of structural plane roughness coefficient lower limit；It is upper and lower for structural plane roughness coefficient under a certain size
2 times of the difference of the major axis of the elliptic equation fitting of limit；For the ellipse of structural plane roughness coefficient upper and lower limit under a certain size
2 times of the difference of the short axle of equation model, I is indeterminacy section scope [0,0.5]；θ is measurement direction；
In the step (6), middle intelligence number function representation is：
1
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610998576.0A CN107036561B (en) | 2016-11-14 | 2016-11-14 | The anisotropic approximate expression method of structural plane roughness based on middle intelligence number function |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610998576.0A CN107036561B (en) | 2016-11-14 | 2016-11-14 | The anisotropic approximate expression method of structural plane roughness based on middle intelligence number function |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107036561A true CN107036561A (en) | 2017-08-11 |
| CN107036561B CN107036561B (en) | 2019-05-31 |
Family
ID=59530349
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610998576.0A Active CN107036561B (en) | 2016-11-14 | 2016-11-14 | The anisotropic approximate expression method of structural plane roughness based on middle intelligence number function |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107036561B (en) |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108445188A (en) * | 2018-04-16 | 2018-08-24 | 浙江大学 | The following slope stability expression of harshness coefficient size effect of rock structural face based on middle intelligence interval function |
| CN110263393A (en) * | 2019-06-04 | 2019-09-20 | 绍兴文理学院 | Rock mass discontinuity anisotropy evaluation method based on anisotropy index of variability AVC |
| CN110348051A (en) * | 2019-06-04 | 2019-10-18 | 绍兴文理学院 | Rock structural face anisotropic character evaluation method based on two-dimensional appearance parameter |
| CN110362857A (en) * | 2019-06-04 | 2019-10-22 | 绍兴文理学院 | The theoretical model analytic method of rock structural face anisotropic elastic solid feature |
| CN110415283A (en) * | 2019-07-03 | 2019-11-05 | 绍兴文理学院 | Parse the fractal evaluation method of rock mass discontinuity anisotropic dimensions effect character |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1645050A (en) * | 2004-12-15 | 2005-07-27 | 金华职业技术学院 | Making method for harshness coefficient table of rock structural face |
| CN103886613A (en) * | 2014-04-10 | 2014-06-25 | 内蒙古科技大学 | Rock structure surface shape anisotropism evaluation method |
| CN104881564A (en) * | 2015-03-09 | 2015-09-02 | 绍兴文理学院 | Establishing method of joint roughness coefficient size effect probability density function model |
-
2016
- 2016-11-14 CN CN201610998576.0A patent/CN107036561B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN1645050A (en) * | 2004-12-15 | 2005-07-27 | 金华职业技术学院 | Making method for harshness coefficient table of rock structural face |
| CN103886613A (en) * | 2014-04-10 | 2014-06-25 | 内蒙古科技大学 | Rock structure surface shape anisotropism evaluation method |
| CN104881564A (en) * | 2015-03-09 | 2015-09-02 | 绍兴文理学院 | Establishing method of joint roughness coefficient size effect probability density function model |
Non-Patent Citations (1)
| Title |
|---|
| 黄江 等: "三维激光扫描技术在高边坡危岩体调查中的应用与讨论", 《长江科学院院报》 * |
Cited By (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108445188A (en) * | 2018-04-16 | 2018-08-24 | 浙江大学 | The following slope stability expression of harshness coefficient size effect of rock structural face based on middle intelligence interval function |
| CN110263393A (en) * | 2019-06-04 | 2019-09-20 | 绍兴文理学院 | Rock mass discontinuity anisotropy evaluation method based on anisotropy index of variability AVC |
| CN110348051A (en) * | 2019-06-04 | 2019-10-18 | 绍兴文理学院 | Rock structural face anisotropic character evaluation method based on two-dimensional appearance parameter |
| CN110362857A (en) * | 2019-06-04 | 2019-10-22 | 绍兴文理学院 | The theoretical model analytic method of rock structural face anisotropic elastic solid feature |
| CN110263393B (en) * | 2019-06-04 | 2023-05-16 | 绍兴文理学院 | Rock mass structural plane anisotropy evaluation method based on anisotropic mutation index AVC |
| CN110348051B (en) * | 2019-06-04 | 2023-05-16 | 绍兴文理学院 | Rock structural surface anisotropy characteristic evaluation method based on two-dimensional morphological parameters |
| CN110362857B (en) * | 2019-06-04 | 2023-05-16 | 绍兴文理学院 | Theoretical model analysis method for anisotropic distribution characteristics of rock structural surface |
| CN110415283A (en) * | 2019-07-03 | 2019-11-05 | 绍兴文理学院 | Parse the fractal evaluation method of rock mass discontinuity anisotropic dimensions effect character |
| CN110415283B (en) * | 2019-07-03 | 2022-04-01 | 绍兴文理学院 | Fractal evaluation method for analyzing anisotropic size effect characteristics of rock mass structural plane |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107036561B (en) | 2019-05-31 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107036561B (en) | The anisotropic approximate expression method of structural plane roughness based on middle intelligence number function | |
| CN106482674B (en) | The approximate expression method of structural plane roughness dimensional effect based on middle intelligence number function | |
| CN104834806B (en) | Structural plane roughness coefficient dimensional effect sample representation evaluation method | |
| CN113587835A (en) | Method for checking and accepting bridge engineering quality by using three-dimensional laser scanning technology | |
| CN106644206B (en) | A kind of pressure from surrounding rock calculation method for shallow tunnel | |
| CN103605135B (en) | A kind of road feature extraction method based on section subdivision | |
| CN110780307B (en) | Method for obtaining road cross section based on storage battery car-mounted laser point cloud mobile measurement system | |
| CN104613904B (en) | Comprehensive determination method for coefficient of roughness of rock mass structure surface | |
| Lu et al. | New study of etching patterns of acid-fracture surfaces and relevant conductivity | |
| CN108445188A (en) | The following slope stability expression of harshness coefficient size effect of rock structural face based on middle intelligence interval function | |
| CN103775071B (en) | Adopt the measuring method that coal and rock crack develops | |
| CN103114848A (en) | Formation fracture space reconstruction method based on rock core measure | |
| CN108871241B (en) | A kind of method that degree of roughness changes after the coarse crack acid etch of quantitative assessment | |
| CN106769276B (en) | Three-dimensional structure face aliquot part choosing method based on Dice similarity measure | |
| CN105606463B (en) | A kind of rock mass structural plane shearing strength integrated evaluating method based on middle intelligence function | |
| CN110457795A (en) | High-rise charming appearance and behaviour displacement state appraisal procedure neural network based | |
| CN105603956B (en) | The determination method of strong rammer groundwork ultimate bearing capacity | |
| CN107563087A (en) | Structural plane roughness coefficient statistical method under optional sampling spacing condition | |
| CN106570338A (en) | Method for determining contour sampling precision in structural plane roughness size effect research | |
| CN104881564B (en) | The construction method of structural plane roughness coefficient dimensional effect probability density estimation | |
| CN108460422B (en) | Submarine geomorphy kind identification method based on depth distribution feature | |
| CN109543236A (en) | Method is determined based on the rock structural plane roughness statistical sample number of variation lines several levels score analysis | |
| CN111199109B (en) | Coupling method for dividing rock mass homogeneous region based on box-counting dimension and cluster analysis | |
| CN106846481B (en) | Geological profile generation method | |
| CN107340159B (en) | Method for selecting representative sample of shear strength of three-dimensional structural surface |
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 | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Du Shigui Inventor after: Li Guoping Inventor after: Yong Rui Inventor after: Yang Xiaocong Inventor after: Xia Caichu Inventor after: Liu Wenlian Inventor after: Dai Yongxin Inventor after: Ma Chengrong Inventor after: Huang Man Inventor after: Ye Jun Inventor before: Du Shigui Inventor before: Fu Xi Inventor before: Yong Rui Inventor before: Ye Jun Inventor before: Li Bo Inventor before: Xia Caichu Inventor before: Zhang Guozhu Inventor before: Huang Man Inventor before: Ma Chengrong Inventor before: He Zhihai |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |