CN107563087A - Structural plane roughness coefficient statistical method under optional sampling spacing condition - Google Patents
Structural plane roughness coefficient statistical method under optional sampling spacing condition Download PDFInfo
- Publication number
- CN107563087A CN107563087A CN201710822964.8A CN201710822964A CN107563087A CN 107563087 A CN107563087 A CN 107563087A CN 201710822964 A CN201710822964 A CN 201710822964A CN 107563087 A CN107563087 A CN 107563087A
- Authority
- CN
- China
- Prior art keywords
- value
- roughness
- structural plane
- jrc
- fractal dimension
- 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
Abstract
A kind of structural plane roughness coefficient statistical method under optional sampling spacing condition, comprises the following steps:1) according to optional test direction, the two-dimensional silhouette hatching that n bars are parallel to each other is drawn out along experiment direction, the drawing for drawing contour curve is scanned and is converted to TIF forms;Matlab programs are write according to yardstick method, wherein i-th two dimensional cross-section line is chosen, calculates corresponding fractal dimension D;2) the fractal dimension D value for the two dimensional cross-section curve being calculated is substituted into formula (2), calculates roughness value characteristic value JRC corresponding to i-th two dimensional cross-section linei;3) and then to other 1 curve of n, also according to step 2) and step 3), corresponding roughness value characteristic value JRC is calculatedi, structural plane is finally counted along experiment direction, roughness value average value when sampling interval is Δ x.The present invention more accurately can describe rock structural plane roughness using fractal dimension.
Description
Technical field
The present invention relates to the roughness value statistical method under a kind of spacing condition based on optional sampling, is adopted suitable for basis
Sample spacing estimates the occasion of structural plane roughness coefficient with 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.
After Turk and Carr research, fractal dimension is considered as a conjunction for quantifying rock texture surface roughness
Suitable parameter.Lee etc. calculates fractal dimension D using five 2,4,6,8,10mm sampling intervals, and Bae analyzes sampling interval and is
1st, 2,4,8,16,32, resulting fractal dimension D value during 64mm, and Zhu Yuxue, Li Yanrong etc. be not in the literature to being provided
Sampler body spacing.The fractal dimension D being calculated under more sampling interval is analyzed in studying above, but is not provided suitably
Sampling interval.There has been no influence of researcher's system research sampling interval to fractal dimension D value in analyzing and researching above, this is also
The result of calculation of different researchers is different or even where the reason for mutual conflict.
Change of the seminar for the rock mass discontinuity two-dimensional silhouette fractal dimension of curve under the conditions of 91 sampling interval sections
Law analysis and research show that, when sampling interval r is 0.3,0.6,0.9,1.2mm, the fractal dimension D value calculated has more
Representativeness, now, formula JRC=1126.95D-1127.50 can preferably description scheme surface roughness JRC and fractal dimension Ds
Between relation.
The content of the invention
In order to overcome existing structure surface roughness coefficient fractal evaluation method can not description scheme surface roughness exactly
Deficiency, the present invention provide the structural plane roughness coefficient statistical method under a kind of spacing condition based on optional sampling, can be more
Rock structural plane roughness accurately is described using fractal dimension.
The technical solution adopted for the present invention to solve the technical problems is:
A kind of structural plane roughness coefficient statistical method under optional sampling spacing condition, the statistical method include following
Step:
1) according to optional test direction, according to certain interval, draw out that n bars are parallel to each other along experiment direction two
Profile cross section line is tieed up, is scanned the drawing for drawing contour curve using large scale scanner, and be converted to TIF forms;
Matlab programs are write according to yardstick method, choose wherein i-th two dimensional cross-section line, wherein i is less than or equal to n, sets
Sampling interval is r, and C is undetermined constant, and corresponding fractal dimension D is calculated according to equation below;
LogN (r)=- Dlogr+C (1)
2) the fractal dimension D value for the two dimensional cross-section curve being calculated is substituted into formula (2), calculates i-th two dimension and cut open
Roughness value characteristic value JRC corresponding to upper threadi;
JRC=1126.95D-1127.50 (2)
3) and then to other n-1 bars curves, also according to step 2) and step 3), corresponding roughness value is calculated
Characteristic value JRCi, structural plane is finally counted along experiment direction, roughness value average value when sampling interval is Δ xWherein:
The present invention technical concept be:Based under different sampling interval sections, the FRACTAL DIMENSION of rock mass discontinuity contour curve
Several changing rules, propose a kind of statistical method of the rock structural plane roughness coefficient under optional sampling spacing.
Beneficial effects of the present invention are mainly manifested in:The fractal dimension calculated is more representative, can be more accurate
Rock structural plane roughness is described using fractal dimension.
Brief description of the drawings
Fig. 1 is the schematic diagram of Barton10 bar nominal contour curves.
Fig. 2 is Barton10 bar standard contour lines JRC calculated value under optional sampling spacing condition.
Embodiment
The invention will be further described below in conjunction with the accompanying drawings.
Referring to Figures 1 and 2, the structural plane roughness coefficient statistical method under a kind of optional sampling spacing condition, the system
Meter method comprises the following steps:
1) according to optional test direction, according to certain interval, draw out that n bars are parallel to each other along experiment direction two
Profile cross section line is tieed up, is scanned the drawing for drawing contour curve using large scale scanner, and be converted to TIF forms;
Matlab programs are write according to yardstick method, choose wherein i-th two dimensional cross-section line, wherein i is less than or equal to n, sets
Sampling interval is r, and C is undetermined constant, and corresponding fractal dimension D is calculated according to equation below;
LogN (r)=- Dlogr+C (1)
2) the fractal dimension D value for the two dimensional cross-section curve being calculated is substituted into formula (2), calculates i-th two dimension and cut open
Roughness value characteristic value JRC corresponding to upper threadi;
JRC=1126.95D-1127.50 (2)
3) and then to other n-1 bars curves, also according to step 2) and step 3), corresponding roughness value is calculated
Characteristic value JRCi, structural plane is finally counted along experiment direction, roughness value average value when sampling interval is Δ xWherein:
Example 1:Analyzed from Barton 10 nominal contour curve comparisons.A kind of knot under optional sampling spacing condition
Structure surface roughness coefficients statistics method, the statistical method comprise the following steps:
1) in order to verify the feasibility of proposed structural plane roughness coefficient statistical method, Barton 10 standards are chosen
Contour curve (size is 10cm, sees Fig. 1) is analyzed.To set sampling interval r be 0.3,0.6,0.9,1.2mm, press
Calculate the fractal dimension D of ten two dimensional cross-section lines respectively according to formula (1):
2) the fractal dimension D value of be calculated ten contour curves is substituted into formula (2), calculates selected two dimension
Roughness value characteristic value JRC corresponding to hatchingi, wherein i=1~10, numbering is represented respectively as 1~10 different roughness
Curve.
Table 1 is r=0.3, when 0.6,0.9,1.2mm, the roughness value feature of Barton 10 nominal contour curves
Value JRC.
Numbering | JRC inverse values | JRC calculated values | Relative error/% |
1 | 0.4 | 0.39 | 2.78 |
2 | 2.8 | 2.72 | 2.77 |
3 | 5.8 | 5.68 | 2.12 |
4 | 6.7 | 6.89 | 2.88 |
5 | 9.5 | 9.26 | 2.48 |
6 | 10.8 | 11.13 | 3.04 |
7 | 12.8 | 13.13 | 2.56 |
8 | 14.5 | 14.15 | 2.40 |
9 | 16.7 | 17.12 | 2.50 |
10 | 18.7 | 19.25 | 2.96 |
Table 1
3) by the calculated value of JRC in table 1, according to formula δ=[JRCCalculated value-JRCInverse value)/JRCInverse value] × 100% calculates
Relative error corresponding to the roughness value characteristic value JRC under optional sampling spacing condition is obtained, as shown in Figure 2.It can be found that
The JRC values that relational expression between the JRC and D that newly propose is estimated and the goodness of fit for testing inverse value are higher, and relative error exists
Within 5%.Therefore, the result that new formula is calculated can preferably meet to require, JRC values are more representative.
Example 2:The protolith structural plane of selection comes from Zhejiang Province Changshan County, Zhejiang Province, belongs to calcareous slate, has tabular knot
It structure, can be flaked along foliation direction, the rock structural face of better quality can be obtained, used for research.
Structural plane roughness coefficient statistical method under a kind of optional sampling spacing condition of the present embodiment, including following step
Suddenly:
1) the calcareous slate structural plane (planar dimension is 1100mm × 1100mm) for choosing Changshan County, Zhejiang Province is protolith examination
Sample, along direction initialization, it is that interval is oriented measurement to structural plane with 15 °, obtains 24 two-dimensional silhouette curves, swept using large-scale
Retouching instrument has structural plane contour line to be scanned drafting, and is converted to TIF forms.
2) the sampling interval r of setting structure face two-dimensional silhouette curve be 0.3,0.6,0.9,1.2mm, according to formula (1) point
Not Ji Suan 24 two-dimensional silhouette curves fractal dimension D.
3) the fractal dimension D value for the two-dimensional silhouette curve being calculated is substituted into formula (2), two-dimensional silhouette selected by calculating
Roughness value characteristic value JRC corresponding to curve1..., JRC24;Then roughness system of the structural plane along experiment direction is counted
Number average value,The as roughness value of structural plane.
Claims (1)
- A kind of 1. structural plane roughness coefficient statistical method under optional sampling spacing condition, it is characterised in that:The statistics side Method comprises the following steps:1) according to optional test direction, according to certain interval, draw out the two dimension that n bars are parallel to each other along experiment direction and take turns Wide hatching, the drawing for drawing contour curve is scanned using large scale scanner, and is converted to TIF forms;Matlab programs are write according to yardstick method, choose wherein i-th two dimensional cross-section line, wherein i is less than or equal to n, sets sampling Spacing is r, and C is undetermined constant, and corresponding fractal dimension D is calculated according to equation below;LogN (r)=- Dlogr+C (1)2) the fractal dimension D value for the two dimensional cross-section curve being calculated is substituted into formula (2), calculates i-th two dimensional cross-section line Corresponding roughness value characteristic value JRCi;JRC=1126.95D-1127.50 (2)3) and then to other n-1 bars curves, also according to step 2) and step 3), corresponding roughness value feature is calculated Value JRCi, structural plane is finally counted along experiment direction, roughness value average value when sampling interval is Δ xIts In:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710822964.8A CN107563087B (en) | 2017-09-13 | 2017-09-13 | Structural surface roughness coefficient statistical method under optimal sampling interval condition |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710822964.8A CN107563087B (en) | 2017-09-13 | 2017-09-13 | Structural surface roughness coefficient statistical method under optimal sampling interval condition |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107563087A true CN107563087A (en) | 2018-01-09 |
CN107563087B CN107563087B (en) | 2020-12-15 |
Family
ID=60980743
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710822964.8A Active CN107563087B (en) | 2017-09-13 | 2017-09-13 | Structural surface roughness coefficient statistical method under optimal sampling interval condition |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107563087B (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109099880A (en) * | 2018-07-12 | 2018-12-28 | 杜时贵 | Rock structural plane roughness coefficient universe searches for measurement method |
CN111512434A (en) * | 2018-02-19 | 2020-08-07 | 富士电机株式会社 | Semiconductor module and method for manufacturing the same |
CN113378909A (en) * | 2021-06-07 | 2021-09-10 | 武汉科技大学 | Grading characterization method, device and medium for roughness coefficient of rock joint surface |
CN114152221A (en) * | 2021-10-27 | 2022-03-08 | 北京工业大学 | Fractal dimension-based combined material contact surface roughness determination method |
CN115930847A (en) * | 2022-09-30 | 2023-04-07 | 中国科学院武汉岩土力学研究所 | Quantitative determination method for roughness evaluation index of three-dimensional structure surface |
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 |
CN1654954A (en) * | 2005-01-12 | 2005-08-17 | 杜时贵 | Method for taking values of mechanical and hydraulic properties of rock mass structural plane |
CN103558094A (en) * | 2013-09-23 | 2014-02-05 | 绍兴文理学院 | Method for representatively sampling subsize rock model structural surface sample sampled based on layering probability |
CN105758361A (en) * | 2016-02-01 | 2016-07-13 | 绍兴文理学院 | Quantitative evaluation method of structural surface roughness coefficient anisotropy |
-
2017
- 2017-09-13 CN CN201710822964.8A patent/CN107563087B/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 |
CN1654954A (en) * | 2005-01-12 | 2005-08-17 | 杜时贵 | Method for taking values of mechanical and hydraulic properties of rock mass structural plane |
CN103558094A (en) * | 2013-09-23 | 2014-02-05 | 绍兴文理学院 | Method for representatively sampling subsize rock model structural surface sample sampled based on layering probability |
CN105758361A (en) * | 2016-02-01 | 2016-07-13 | 绍兴文理学院 | Quantitative evaluation method of structural surface roughness coefficient anisotropy |
Non-Patent Citations (2)
Title |
---|
杜时贵等: "岩体结构面起伏幅度尺寸效应的试验研究", 《工程地质学报》 * |
罗战友等: "岩石结构面粗糙度系数尺寸效应的推拉试验研究", 《岩土力学》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111512434A (en) * | 2018-02-19 | 2020-08-07 | 富士电机株式会社 | Semiconductor module and method for manufacturing the same |
CN111512434B (en) * | 2018-02-19 | 2023-08-08 | 富士电机株式会社 | Semiconductor module and method for manufacturing the same |
US11749581B2 (en) | 2018-02-19 | 2023-09-05 | Fuji Electric Co., Ltd. | Semiconductor module and method for manufacturing same |
CN109099880A (en) * | 2018-07-12 | 2018-12-28 | 杜时贵 | Rock structural plane roughness coefficient universe searches for measurement method |
CN109099880B (en) * | 2018-07-12 | 2020-08-11 | 宁波大学 | Rock mass structural plane roughness coefficient global search measuring method |
CN113378909A (en) * | 2021-06-07 | 2021-09-10 | 武汉科技大学 | Grading characterization method, device and medium for roughness coefficient of rock joint surface |
CN113378909B (en) * | 2021-06-07 | 2022-08-12 | 武汉科技大学 | Grading characterization method, device and medium for roughness coefficient of rock joint surface |
CN114152221A (en) * | 2021-10-27 | 2022-03-08 | 北京工业大学 | Fractal dimension-based combined material contact surface roughness determination method |
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 |
Also Published As
Publication number | Publication date |
---|---|
CN107563087B (en) | 2020-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107563087A (en) | Structural plane roughness coefficient statistical method under optional sampling spacing condition | |
CN104834806B (en) | Structural plane roughness coefficient dimensional effect sample representation evaluation method | |
CN107656902A (en) | Structural plane roughness coefficient statistical method under different sampling intervals | |
CN101050958A (en) | Measuring method for large size structure surface contour curve roughness coefficient | |
CN106482674B (en) | The approximate expression method of structural plane roughness dimensional effect based on middle intelligence number function | |
CN104833333B (en) | Structural plane roughness coefficient dimensional effect sample characterization unit determines method | |
CN108534843B (en) | Single gate flow calibration method and device for open channel water delivery channel | |
CN105606463B (en) | A kind of rock mass structural plane shearing strength integrated evaluating method based on middle intelligence function | |
CN106570338B (en) | Method for determining sampling precision of contour line in roughness size effect of rock mass structural plane | |
CN107036561B (en) | The anisotropic approximate expression method of structural plane roughness based on middle intelligence number function | |
CN115060452B (en) | Panoramic error detection method applied to large wind tunnel spray pipe molded surface | |
CN103902798B (en) | Data preprocessing method | |
CN104881564B (en) | The construction method of structural plane roughness coefficient dimensional effect probability density estimation | |
CN107843215B (en) | Based on the roughness value fractal evaluation model building method under optional sampling spacing condition | |
CN109543236B (en) | Method for determining roughness statistics sample number of rock mass structural plane | |
WO2022021967A1 (en) | Method for high-precision measurement of degree of bending and degree of twisting of connecting rod | |
CN113868926A (en) | Method for constructing spatial distribution model of water quality parameters of culture pond | |
CN104765013B (en) | Calibrating three-axle magnetic sensor method | |
CN1815184A (en) | Non-homogenous structural shear-strength sampling method | |
CN206235284U (en) | A kind of dedicated slot model | |
CN113291207A (en) | Dynamic measurement method of rigid subway overhead line system | |
CN110516281B (en) | Evaluation method for quantitatively representing anisotropic size effect of rock mass structural plane | |
CN103675218B (en) | A kind of papermaking-method reconstituted tobaccos slurry uniformity coefficient characterizing method | |
Wang et al. | Retracted: One-Variable Linear Regression Mathematical Model of Color Reading and Material Concentration Identification | |
CN206291836U (en) | A kind of special angular templet |
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 |