CN107843215A - Based on the roughness value fractal evaluation model building method under optional sampling spacing condition - Google Patents
Based on the roughness value fractal evaluation model building method under optional sampling spacing condition Download PDFInfo
- Publication number
- CN107843215A CN107843215A CN201710823597.3A CN201710823597A CN107843215A CN 107843215 A CN107843215 A CN 107843215A CN 201710823597 A CN201710823597 A CN 201710823597A CN 107843215 A CN107843215 A CN 107843215A
- Authority
- CN
- China
- Prior art keywords
- fractal dimension
- sampling interval
- value
- sampling
- spacing
- 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
Abstract
A kind of roughness value fractal evaluation model building method under spacing condition based on optional sampling, comprises the following steps：1) the high pixel picture of Barton m bar standard contour lines is chosen, extracts the data message of m bar canonical profile lines in picture, the profile of m bar canonical profile lines is obtained and carries out post-processing；2) program is write in matlab according to yardstick method, by the structural plane contour line information obtained in (1) import, different sampling interval r is set, obtain m bar canonical profile lines corresponding to N (r) values, calculate fractal dimension D；3) dispersion degree of fractal dimension D value under 91 sampling interval sections is analyzed；4) changing rule of fractal dimension D value under different sampling interval intervals is explored, draws optional sampling spacing r, builds the function model of structural plane roughness coefficient under optional sampling spacing.The present invention can more precisely estimate structural plane roughness coefficient JRC according to optional sampling spacing.
Description
Technical field
Consider that sampling interval evaluates the model structure influenceed to structural plane roughness coefficient fractal dimension the present invention relates to a kind of
Construction method, suitable for estimating the occasion of structural plane roughness coefficient according to sampling interval and fractal dimension D value.
Background technology
Rock mass discontinuity controls the mechanical properties such as the deformation and failure of rock mass, and the mechanical property of rock mass discontinuity and its
Surface topography is closely related.Characteristic parameters of the JRC as description scheme face surface topography, by the extensive concern of scholars.Close
Mainly have in rock structural plane roughness coefficient JRC evaluation method following several：Empirical estimation method, statistical parameter method, straight flange
Method and fractal dimension method etc..A kind of quantitative method of the fractal dimension method as evaluation rock mass discontinuity surface roughness, by
The extensive concern of scholar is arrived.
Computational methods on fractal dimension D have many kinds, as yardstick method, package topology, variogram method, spectroscopic methodology, from
Affine fractal method, h-L methods etc., wherein in computational methods on two-dimensional fractal dimension D, one of the most commonly used method is yardstick
Method.Lee etc. calculates fractal dimension D using five 2,4,6,8,10mm sampling intervals, Bae analyze sampling interval for 1,2,4,
8th, 16,32, resulting fractal dimension D value during 64mm, and Zhu Yuxue, Li Yanrong etc. do not provide specific sampling in the literature
Spacing.More than study in analyze the fractal dimension D being calculated under more sampling interval, but do not provide between suitable sampling
Away from section.
Week, wound soldier when proposing that critical yardstick takes 1~3mm, can be obtained using different yardstick r measurement joint hatchings
" appropriate fractal dimension ".The method based on linear scale such as Kulatilake proposes the new general of suitably sized scope
Read, Jang etc. applies to this method, obtains preferable effect.More than analyze and research in not yet network analysis sampling interval
Influence to fractal dimension, this be also different researchers result of calculation is different in addition where the reason for mutual conflict.
The content of the invention
In order to overcome existing structural plane roughness coefficient fractal dimension evaluation method can not consider sampling interval pair
Its deficiency influenceed, the present invention provide the structural plane roughness coefficient fractal evaluation mould under a kind of spacing condition based on optional sampling
The construction method of type, structural plane roughness coefficient JRC can more precisely be estimated according to optional sampling spacing.
The technical solution adopted for the present invention to solve the technical problems is：
A kind of roughness value fractal evaluation model building method under spacing condition based on optional sampling, methods described bag
Include following steps：
1) the high pixel picture of Barton m bar standard contour lines is chosen, is carried using " image " function of matlab softwares
The data message of m bars canonical profile line in picture is taken, the profile of m bar canonical profile lines is obtained and carries out post-processing, save as
Jpg forms；
2) program is write in matlab according to yardstick method, the structural plane contour line information obtained in (1) is imported
In matlab softwares, different sampling interval r is set, obtains N (r) values corresponding to m bar canonical profile lines, C is undetermined constant, is pressed
Corresponding fractal dimension D is calculated according to equation below；
LogN (r)=- Dlogr+C (1)
3) according to formula CV=σ/μ, wherein, σ is standard deviation, and μ is average value, and CV is the coefficient of variation, to sampling interval
When section is R=a~bmm, 0.1mm≤a≤9.1mm, 1.0mm≤b≤10.0mm, the dispersion degree of fractal dimension D value are carried out
Analysis, draws：1. the coefficient of variation of fractal dimension D value gradually increases with the increase of sampling interval section and hatching roughness
Greatly；2. as sampling interval section R=0.1mm-1.2mm, the coefficient of variation CV of profile line fractal dimension D values is respectively less than 0.05,
Dispersion degree is relatively low.
4) explore under different sampling interval intervals, the changing rule of fractal dimension D value, draw optional sampling spacing r, and then
The function model of structural plane roughness coefficient under optional sampling spacing can be built：
JRC=1126.95D-1127.50 (2)
Beneficial effects of the present invention are mainly manifested in：(1) it can be considered that sampling interval is to structural plane fractal dimension evaluation side
The influence of method, so that the roughness value characteristic value JRC calculated is more representative.(2) can relatively accurately count
Calculate under optional sampling spacing, JRC values corresponding to standard contour line line, and new JRC formula have preferable applicability.
Brief description of the drawings
Fig. 1 is the schematic diagram of Barton10 bar standard contour lines.
Fig. 2 is the coefficient of variation CV changing rules of Article 7 hatching.
Fig. 3 is the fractal dimension D value contrast of hatching under different sampling intervals.
Embodiment
The invention will be further described below in conjunction with the accompanying drawings.
1~Fig. 3 of reference picture, the structural plane roughness coefficient fractal evaluation model under a kind of spacing condition based on optional sampling
According to method, comprise the following steps：
1) the high pixel picture of Barton m bar standard contour lines is chosen, is carried using " image " function of matlab softwares
The data message of m bars canonical profile line in picture is taken, the profile of m bar canonical profile lines is obtained and carries out post-processing, save as
Jpg forms；
2) program is write in matlab according to yardstick method, the structural plane contour line information obtained in (1) is imported
In matlab softwares, different sampling interval r is set, obtains N (r) values corresponding to m bar canonical profile lines, C is undetermined constant, is pressed
Corresponding fractal dimension D is calculated according to equation below；
LogN (r)=- Dlogr+C (1)
3) according to formula CV=σ/μ, wherein, σ is standard deviation, and μ is average value, and CV is the coefficient of variation, to sampling interval
When section is R=a~bmm, 0.1mm≤a≤9.1mm, 1.0mm≤b≤10.0mm, the dispersion degree of fractal dimension D value are carried out
Analysis.Draw：1. the coefficient of variation of fractal dimension D value gradually increases with the increase of sampling interval section and hatching roughness
Greatly；2. as sampling interval section R=0.1mm-1.2mm, the coefficient of variation CV of profile line fractal dimension D values is respectively less than 0.05,
Dispersion degree is relatively low.
4) explore under different sampling interval intervals, the changing rule of fractal dimension D value, draw optional sampling spacing r, and then
Build the function model of structural plane roughness coefficient under optional sampling spacing：
JRC=1126.95D-1127.50 (2).
It is as follows as research object, embodiment that the present embodiment chooses ten standard contour lines of Barton：
1) the high pixel photo of Barton ten (taking m=10) bar nominal contour curves is chosen respectively, utilizes matlab
The data message of ten canonical profile lines in " image " function extraction picture of software, obtain the profile of ten canonical profile lines
And post-processing is carried out, jpg forms are saved as, as shown in Figure 1；
2) program is write in matlab according to yardstick method, the structural plane contour line information obtained in (1) is imported
In matlab softwares, set sampling interval section be respectively 0.1mm-1.0mm, 0.2mm-1.1mm ..., 9.1mm-10mm, obtain
To N (r) values corresponding to ten canonical profile lines, corresponding fractal dimension D value is calculated according to equation below；
LogN (r)=- Dlogr+C (1)
3) 91 are sampled according to formula CV=σ/μ (wherein, σ is standard deviation, and μ is average value, and CV is the coefficient of variation)
Under Space Interval, the dispersion degree of fractal dimension D value is analyzed, as shown in Figure 2.Draw：1. the variation lines of fractal dimension D value
Number gradually increases with the increase of sampling interval section and hatching roughness；2. in 0.1mm-1.2mm sampling intervals section
Interior, the coefficient of variation CV of profile line fractal dimension D values is respectively less than 0.05, and dispersion degree is relatively low.
4) it is respectively r to set sampling interval1=0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0,1.1,
1.2mm；r2=0.2,0.4,0.6,0.8,1.0,1.2mm；r3=0.3,0.6,0.9,1.2mm；r4=0.6,1.2mm, calculate
To corresponding fractal dimension D, comparing result is shown in Fig. 3, so as to draw optional sampling spacing r=0.3,0.6,0.9,1.2mm；
And then build under optional sampling spacing condition, the model of roughness value JRC and fractal dimension is as follows：
JRC=1126.95D-1127.50 R=0.998 (2)
Wherein, R is fitting parameter.
Claims (1)
- A kind of 1. roughness value fractal evaluation model building method under spacing condition based on optional sampling, it is characterised in that： Methods described comprises the following steps：1) the high pixel picture of Barton m bar standard contour lines is chosen, figure is extracted using " image " function of matlab softwares The data message of m bars canonical profile line in piece, obtain the profile of m bar canonical profile lines and carry out post-processing, save as jpg lattice Formula；2) program is write in matlab according to yardstick method, it is soft that the structural plane contour line information obtained in (1) is imported into matlab In part, different sampling interval r is set, obtains N (r) values corresponding to m bar canonical profile lines, C is undetermined constant, according to following public affairs Formula calculates corresponding fractal dimension D；LogN (r)=- Dlogr+C (1)3) according to formula CV=σ/μ, wherein, σ is standard deviation, and μ is average value, and CV is the coefficient of variation, to sampling interval section For R=a~bmm when, 0.1mm≤a≤9.1mm, 1.0mm≤b≤10.0mm, the dispersion degree of fractal dimension D value are analyzed, Draw：1. the coefficient of variation of fractal dimension D value gradually increases with the increase of sampling interval section and hatching roughness；② As sampling interval section R=0.1mm-1.2mm, the coefficient of variation CV of profile line fractal dimension D values is respectively less than 0.05, discrete journey Spend relatively low；4) explore under different sampling interval intervals, the changing rule of fractal dimension D value, draw optional sampling spacing r, and then build The function model of structural plane roughness coefficient under optional sampling spacing：JRC=1126.95D-1127.50 (2).
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710823597.3A CN107843215B (en) | 2017-09-13 | 2017-09-13 | Based on the roughness value fractal evaluation model building method under optional sampling spacing condition |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710823597.3A CN107843215B (en) | 2017-09-13 | 2017-09-13 | Based on the roughness value fractal evaluation model building method under optional sampling spacing condition |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107843215A true CN107843215A (en) | 2018-03-27 |
| CN107843215B CN107843215B (en) | 2019-10-11 |
Family
ID=61683279
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710823597.3A Active CN107843215B (en) | 2017-09-13 | 2017-09-13 | Based on the roughness value fractal evaluation model building method under optional sampling spacing condition |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107843215B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109460603A (en) * | 2018-11-07 | 2019-03-12 | 绍兴文理学院 | JRC parameter equation modification method based on adaptive correction function |
| CN115930847A (en) * | 2022-09-30 | 2023-04-07 | 中国科学院武汉岩土力学研究所 | Quantitative determination method for roughness evaluation index of three-dimensional structure surface |
| CN116008139A (en) * | 2023-03-27 | 2023-04-25 | 华中科技大学 | Evaluation method and evaluation system for fractal dimension of particles in dispersion system |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5612700A (en) * | 1995-05-17 | 1997-03-18 | Fastman, Inc. | System for extracting targets from radar signatures |
| CN1645049A (en) * | 2004-12-15 | 2005-07-27 | 金华职业技术学院 | Method for determining dimensional size effects of toughness coefficients of typical stone structural surface |
| CN103644866A (en) * | 2013-12-10 | 2014-03-19 | 中国地质大学（武汉） | Rock mass structure surface roughness evaluation method overcoming size effect |
| CN104239711A (en) * | 2014-09-09 | 2014-12-24 | 同济大学 | Method for determining joint roughness |
-
2017
- 2017-09-13 CN CN201710823597.3A patent/CN107843215B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5612700A (en) * | 1995-05-17 | 1997-03-18 | Fastman, Inc. | System for extracting targets from radar signatures |
| CN1645049A (en) * | 2004-12-15 | 2005-07-27 | 金华职业技术学院 | Method for determining dimensional size effects of toughness coefficients of typical stone structural surface |
| CN103644866A (en) * | 2013-12-10 | 2014-03-19 | 中国地质大学（武汉） | Rock mass structure surface roughness evaluation method overcoming size effect |
| CN104239711A (en) * | 2014-09-09 | 2014-12-24 | 同济大学 | Method for determining joint roughness |
Non-Patent Citations (2)
| Title |
|---|
| 孙辅庭等: "Barton 标准剖面 JRC 与独立于离散间距的统计参数关系研究", 《岩石力学与工程学报》 * |
| 温韬等: "基于分形理论的岩体结构面粗糙度影响因素研究", 《人民长江》 * |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109460603A (en) * | 2018-11-07 | 2019-03-12 | 绍兴文理学院 | JRC parameter equation modification method based on adaptive correction function |
| CN109460603B (en) * | 2018-11-07 | 2023-05-16 | 绍兴文理学院 | JRC parameter formula correction method based on self-adaptive correction function |
| CN115930847A (en) * | 2022-09-30 | 2023-04-07 | 中国科学院武汉岩土力学研究所 | Quantitative determination method for roughness evaluation index of three-dimensional structure surface |
| CN115930847B (en) * | 2022-09-30 | 2023-09-22 | 中国科学院武汉岩土力学研究所 | Quantitative determination method for roughness evaluation index of three-dimensional structural surface |
| CN116008139A (en) * | 2023-03-27 | 2023-04-25 | 华中科技大学 | Evaluation method and evaluation system for fractal dimension of particles in dispersion system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107843215B (en) | 2019-10-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107843215A (en) | Based on the roughness value fractal evaluation model building method under optional sampling spacing condition | |
| CN103916896B (en) | Anomaly detection method based on multi-dimensional Epanechnikov kernel density estimation | |
| CN107038292A (en) | A kind of many output of wind electric field correlation modeling methods based on adaptive multivariable nonparametric probability | |
| Park | Testing exponentiality based on the Kullback-Leibler information with the type II censored data | |
| CN101763638B (en) | Method for classifying cerebral white matter fiber tracts in diffusion tensor nuclear magnetic resonance image | |
| Rieder | On convergence rates of inexact Newton regularizations | |
| CN107563087A (en) | Structural plane roughness coefficient statistical method under optional sampling spacing condition | |
| CN103177450A (en) | Image scene segmentation and layering joint solution method based on component set sampling | |
| CN106952339A (en) | A kind of Points Sample method theoretical based on optimal transmission | |
| CN107656902A (en) | Structural plane roughness coefficient statistical method under different sampling intervals | |
| CN103237320B (en) | Wireless sensor network quantizes based on mixing the method for tracking target that Kalman is merged | |
| CN113237885B (en) | Building performance evaluation method based on structural health monitoring data | |
| CN114201924B (en) | Solar irradiance prediction method and system based on transfer learning | |
| CN104933418A (en) | Population size counting method of double cameras | |
| CN107895366B (en) | Imaging method and system for color evaluation and computer readable storage device | |
| CN107036561A (en) | The anisotropic approximate expression method of structural plane roughness based on middle intelligence number function | |
| CN109738723A (en) | A kind of electric energy meter three-phase automatic identification method | |
| CN106525466A (en) | Robust filtering method and system for a key part of motor train unit braking system | |
| CN108491927A (en) | A kind of data processing method and device based on neural network | |
| CN113326744B (en) | Method and system for detecting on-orbit state abnormity of spacecraft | |
| CN102332042B (en) | Modeling method for quartz flexible accelerometer starting model | |
| CN104881564B (en) | The construction method of structural plane roughness coefficient dimensional effect probability density estimation | |
| CN111210051B (en) | User electricity consumption behavior prediction method and system | |
| CN109857975A (en) | A kind of electric energy meter platform area and three-phase automatic identification method | |
| EP1297742A3 (en) | A method of collecting measurement data during automatically milking an animal |
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 |