CN106570338B - 岩体结构面粗糙度尺寸效应中轮廓线采样精度确定方法 - Google Patents

岩体结构面粗糙度尺寸效应中轮廓线采样精度确定方法 Download PDF

Info

Publication number
CN106570338B
CN106570338B CN201610998652.8A CN201610998652A CN106570338B CN 106570338 B CN106570338 B CN 106570338B CN 201610998652 A CN201610998652 A CN 201610998652A CN 106570338 B CN106570338 B CN 106570338B
Authority
CN
China
Prior art keywords
structural surface
contour line
fourier series
size
line
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.)
Active
Application number
CN201610998652.8A
Other languages
English (en)
Other versions
CN106570338A (zh
Inventor
杜时贵
雍睿
杨小聪
夏才初
刘文连
代永新
马成荣
黄曼
叶军
李国平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Shaoxing
Original Assignee
University of Shaoxing
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Shaoxing filed Critical University of Shaoxing
Priority to CN201610998652.8A priority Critical patent/CN106570338B/zh
Publication of CN106570338A publication Critical patent/CN106570338A/zh
Application granted granted Critical
Publication of CN106570338B publication Critical patent/CN106570338B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
    • G16Z99/00—Subject matter not provided for in other main groups of this subclass

Landscapes

  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • Image Analysis (AREA)

Abstract

一种结构面粗糙度尺寸效应研究中轮廓线采样精度确定方法，包括以下步骤：(1)提取系列尺寸结构面轮廓线坐标数据；(2)采用傅里叶级数对结构面轮廓线进行逼近；(3)计算不同阶次条件下的近似拟合值；(4)设置傅里叶级数逼近的阈值；(5)计算不同阶次条件下，傅里叶级数拟合值与结构面实测值的均方差；(6)计算系列尺寸结构面轮廓曲线的傅里叶级数的最低阶次；(7)分析最低阶次平均值与结构面长度的关系；(8)依据傅里叶级数的频率关系，得到最大采样间距；(9)推测给定尺寸L的结构面轮廓线的建议采样精度。本发明能定量确定不同尺寸条件下结构面轮廓线采样精度。

Description

岩体结构面粗糙度尺寸效应中轮廓线采样精度确定方法
技术领域
本发明属于岩石力学工程领域，涉及一种结构面粗糙度尺寸效应研究所涉及的轮廓线采样精度确定方法。
背景技术
岩体结构面粗糙度直接影响着岩体强度、变形以及渗流特征。国内外学者开展了大量关于岩体结构面粗糙度评价方法的研究。为实现结构面表面粗糙性质的定量化评价，学者们一般通过结构面表面轮廓的坐标数据，计算表征粗糙度的各项几何参数，进而建立它们与结构面粗糙度系数的数理统计关系。粗糙度尺寸效应是结构面的固有属性之一，越来越多的国内外学者对其开展了研究。然而，人们分析结构面粗糙度尺寸效应时，往往会采用不同的采样精度，并且没有对采样精度的选用进行分析和解释，这使得粗糙度评价结果存在很大的随意性。另一方面，人们在研究小尺寸结构面粗糙度时，往往采用了非常高的采样精度。倘若采样相同的采样精度研究大尺寸结构面，会产生非常多的数据，影响了粗糙度的计算效果与效率。Fardin等(2004)在结构面表面粗糙度尺寸效应研究中提出了以 0.2mm为10cm结构面轮廓线的采样精度，其他尺寸结构面的采样精度通过等比例扩大得到。这种方法虽然从一定程度上给出了结构面粗糙度尺寸效应研究中轮廓线采样精度的建议方法，但该建议方法依然属于经验方法，缺乏理论依据。因此，为了研究结构面粗糙度尺寸效应，提出一种确定系列尺度结构面轮廓线采样精度的定量方法显得至关重要。
发明内容
为了克服现有技术的无法确定系列尺度结构面轮廓线采样精度的不足，本发明提供了一种系列尺寸结构面采样精度的确定方法，为进一步的结构面粗糙度尺寸效应研究提供了保障。
本发明解决其技术问题所采用的技术方案是：
一种岩体结构面粗糙度尺寸效应研究中轮廓线采样精度确定方法，包括以下步骤：
(1)选定所需要分析的天然结构面，等间距布置大尺寸结构面测线，然后采用高精度三维激光扫描仪进行测量，并提取各条结构面轮廓线上的坐标数据；
(2)对获得的大尺寸结构面进行分割，以轮廓线左端作为起点，提取不同尺寸条件下的轮廓曲线。
(3)采用傅里叶级数对不同尺寸条件下的结构面轮廓线上各点坐标进行逼近,
Figure GDA0001848821010000021
其中，L为结构面测线长度，y是结构面轮廓线上各点的高度，x是轮廓线上各点的水平坐标，a0、an与bn为傅里叶级数的系数，n为阶次；
(4)不断调整n的大小，计算不同阶次条件下傅里叶级数的近似拟合值
Figure GDA0001848821010000022
(5)在结构面轮廓线纵坐标最大值与最小值之间作一条水平直线，该水平线将结构面轮廓线划分为水平线以上部分和水平线以下部分，不断调整水平线的位置，直至水平线上下两部分曲线面积相等。计算结构面轮廓线上各点到该水平线的距离的均方差；
(6)设置傅里叶级数逼近值与结构面轮廓线的拟合系数，譬如5％，将轮廓线上各点到水平线距离的均方差乘以拟合系数作为傅里叶级数逼近的阈值；
(7)按如下公式，分别计算不同阶次条件下相同水平坐标条件下傅里叶级数的近似拟合值
Figure GDA0001848821010000031
与实测结构面轮廓线坐标Yi的均方差：
Figure GDA0001848821010000032
其中，i为结构面轮廓线上坐标点编号；m为结构面轮廓线上的点数；
(8)分析En随阶次n增大的变化规律，判断傅里叶级数逼近的阈值的大小关系；当阶次n增大到某一值时，En小于傅里叶级数逼近的阈值，此时，阶次n为傅里叶级数逼近所要求的最低阶次；
(9)计算得到每条大尺寸结构面轮廓曲线在各个尺寸条件下傅里叶级数的最低阶次；
(10)分析得到各个尺寸条件下傅里叶级数的最低阶次平均值与结构面轮廓曲线长度的关系f(L)；
(11)依据香农采样定理，只有当采样频率大于2倍以上图形的最高频率时，图像的基本信息才能得到准确表达。依据傅里叶级数的频率关系，得到最大采样间距的表达式：
Figure GDA0001848821010000033
(12)根据尺寸效应研究中结构面轮廓线尺寸L，推测建议采样精度。
本发明的有益效果主要表现在：定量确定系列尺度结构面轮廓线采样精度。
附图说明
图1为实测结构面起伏轮廓曲线图。
图2为最小阶次n的平均值与结构面尺寸系数的关系曲线图。
具体实施方式
下面结合附图对本发明作进一步描述。
参照图1和图2，一种岩体结构面粗糙度尺寸效应研究中轮廓线采样精度确定方法，包括以下步骤：
(1)选定所需要分析的天然结构面，等间距布置大尺寸结构面测线，然后采用高精度三维激光扫描仪进行测量，并提取各条结构面轮廓线上的坐标数据；
(2)对获得的大尺寸结构面进行分割，以轮廓线左端作为起点，提取不同尺寸条件下的轮廓曲线。
(3)采用傅里叶级数对不同尺寸条件下的结构面轮廓线上各点坐标进行逼近,
Figure GDA0001848821010000041
其中，L为结构面测线长度，y是结构面轮廓线上各点的高度，x是轮廓线上各点的水平坐标，a0、an与bn为傅里叶级数的系数，n为阶次；
(4)不断调整n的大小，计算不同阶次条件下傅里叶级数的近似拟合值
Figure GDA0001848821010000042
(5)在结构面轮廓线纵坐标最大值与最小值之间作一条水平直线，该水平线将结构面轮廓线划分为水平线以上部分和水平线以下部分，不断调整水平线的位置，直至水平线上下两部分曲线面积相等。计算结构面轮廓线上各点到该水平线的距离的均方差；
(6)设置傅里叶级数逼近值与结构面轮廓线的拟合系数，譬如5％，将轮廓线上各点到水平线距离的均方差乘以拟合系数作为傅里叶级数逼近的阈值；
(7)按如下公式，分别计算不同阶次条件下相同水平坐标条件下傅里叶级数的近似拟合值
Figure GDA0001848821010000043
与实测结构面轮廓线坐标Yi的均方差：
Figure GDA0001848821010000051
其中，i为结构面轮廓线上坐标点编号；m为结构面轮廓线上的点数。
(8)分析En随阶次n增大的变化规律，判断傅里叶级数逼近的阈值的大小关系。当阶次n增大到某一值时，En小于傅里叶级数逼近的阈值，此时阶次n为傅里叶级数逼近所要求的最低阶次。
(9)计算得到每条大尺寸结构面轮廓曲线在各个尺寸条件下傅里叶级数的最低阶次。
(10)分析得到各个尺寸条件下傅里叶级数的最低阶次平均值与结构面轮廓曲线长度的关系f(L)。
(11)依据香农采样定理，只有当采样频率大于2倍以上图形的最高频率时，图像的基本信息才能得到准确表达。依据傅里叶级数的频率关系，得到最大采样间距的表达式：
Figure GDA0001848821010000052
(12)根据尺寸效应研究中结构面轮廓线尺寸L，推测建议采样精度。
实例：一种结构面粗糙度尺寸效应研究中轮廓线采样精度确定方法，包括以下步骤：
(1)选定浙江省常山县板岩作为所需要分析的天然结构面，等间距布置3 条典型大尺寸结构面测线，然后采用高精度三维激光扫描仪进程测量，各结构面轮廓线的坐标数据如图1所示；
(2)对获得的1000mm长的大尺寸结构面进行分割，以轮廓线左端作为起点，分别提取100mm～1000mm尺寸条件下的轮廓曲线。
(3)采用傅里叶级数对所得到的结构面轮廓线上各点坐标进行逼近，不断调整n的大小，计算不同阶次条件下傅里叶级数的近似拟合值
Figure GDA0001848821010000053
(4)在结构面轮廓线纵坐标最大值与最小值之间作一条水平直线，不断调整水平线的位置，直至水平线上下两部分曲线面积相等。计算结构面轮廓线上各点到该水平线的距离的均方差等于Emax；将轮廓线上各点到水平线距离的均方差乘以拟合系数5％得到傅里叶级数逼近的阈值 5％Emax；
(5)计算不同阶次条件下相同水平坐标条件下傅里叶级数的近似拟合值
Figure GDA0001848821010000061
与实测结构面轮廓线坐标Yi的均方差En。
(6)计算当En等于傅里叶级数逼近的阈值5％Emax时的阶次，并将其作为傅里叶级数逼近所要求的最低阶次n。三条大尺寸结构面轮廓曲线的最低阶次计算结果如表1-3所示。表1为各个尺寸条件下结构面轮廓曲线1 最低阶次的计算结果。
Figure GDA0001848821010000062
表1
表2为各个尺寸条件下结构面轮廓曲线2最低阶次的计算结果：
Figure GDA0001848821010000063
Figure GDA0001848821010000071
表2
表3为各个尺寸条件下结构面轮廓曲线3最低阶次的计算结果：
Figure GDA0001848821010000072
表3
(7)分析得到各个尺寸条件下傅里叶级数的最低阶次平均值与结构面轮廓曲线长度的关系f(L)，如图2所示。
(8)依据香农采样定理，只有当采样频率大于2倍以上图形的最高频率时，图像的基本信息才能得到准确表达。依据傅里叶级数的频率关系，得到最大采样间距的表达式：
Figure GDA0001848821010000081
其中，L10＝10cm。
(9)根据尺寸效应研究中结构面轮廓线尺寸L，譬如L＝60cm，带入上述公式计算得到采样精度建议值应小于2.716mm。

Claims (1)

1.一种岩体结构面粗糙度尺寸效应中轮廓线采样精度确定方法，其特征在于：所述确定方法包括以下步骤：
(1)选定所需要分析的天然结构面，等间距布置大尺寸结构面测线，然后采用三维激光扫描仪进行测量，并提取各条结构面轮廓线上的坐标数据；
(2)对获得的大尺寸结构面进行分割，以轮廓线左端作为起点，提取不同尺寸条件下的轮廓曲线；
(3)采用傅里叶级数对不同尺寸条件下的结构面轮廓线上各点坐标进行逼近,
Figure FDA0002319324380000011
其中，L为结构面轮廓曲线长度，y是结构面轮廓线上各点的高度，x是轮廓线上各点的水平坐标，a0、an与bn为傅里叶级数的系数，n为阶次；
(4)不断调整n的大小，计算不同阶次条件下傅里叶级数的近似拟合值
Figure FDA0002319324380000015
(5)在结构面轮廓线纵坐标最大值与最小值之间作一条水平直线，该水平线将结构面轮廓线划分为水平线以上部分和水平线以下部分，不断调整水平线的位置，直至水平线上下两部分曲线面积相等；计算结构面轮廓线上各点到该水平线的距离的均方差；
(6)设置傅里叶级数逼近值与结构面轮廓线的拟合系数，将轮廓线上各点到水平线距离的均方差乘以拟合系数作为傅里叶级数逼近的阈值；
(7)按如下公式，分别计算不同阶次条件下相同水平坐标条件下傅里叶级数的近似拟合值
Figure FDA0002319324380000012
与实测结构面轮廓线坐标Yi的均方差：
Figure FDA0002319324380000013
其中，i为结构面轮廓线上坐标点编号；m为结构面轮廓线上的坐标点总数；
(8)分析En随阶次n增大的变化规律，判断傅里叶级数逼近的阈值的大小关系；当阶次n增大到临界值时，En小于傅里叶级数逼近的阈值，此时阶次n为傅里叶级数逼近所要求的最低阶次；
(9)计算得到每条大尺寸结构面轮廓曲线在各个尺寸条件下傅里叶级数的最低阶次；
(10)分析得到各个尺寸条件下傅里叶级数的最低阶次平均值与结构面轮廓曲线长度的关系f(L)；
(11)依据傅里叶级数的频率关系，得到最大采样间距的表达式：
Figure FDA0002319324380000014
(12)根据尺寸效应研究中结构面轮廓线尺寸L，推测建议采样精度。
CN201610998652.8A 2016-11-14 2016-11-14 岩体结构面粗糙度尺寸效应中轮廓线采样精度确定方法 Active CN106570338B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610998652.8A CN106570338B (zh) 2016-11-14 2016-11-14 岩体结构面粗糙度尺寸效应中轮廓线采样精度确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610998652.8A CN106570338B (zh) 2016-11-14 2016-11-14 岩体结构面粗糙度尺寸效应中轮廓线采样精度确定方法

Publications (2)

Publication Number Publication Date
CN106570338A CN106570338A (zh) 2017-04-19
CN106570338B true CN106570338B (zh) 2020-05-08

Family

ID=58542631

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610998652.8A Active CN106570338B (zh) 2016-11-14 2016-11-14 岩体结构面粗糙度尺寸效应中轮廓线采样精度确定方法

Country Status (1)

Country Link
CN (1) CN106570338B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107577884B (zh) * 2017-09-13 2020-10-27 绍兴文理学院 结构面粗糙度系数统计参数z2与采样间距关系模型的构建方法
CN109099880B (zh) * 2018-07-12 2020-08-11 宁波大学 岩体结构面粗糙度系数全域搜索测量方法
CN109029334B (zh) * 2018-07-12 2020-06-02 宁波大学 岩体结构面粗糙度系数尺寸效应全域搜索测量方法
CN109084708B (zh) * 2018-07-25 2020-04-21 深圳大学 颗粒二维轮廓表面整体粗糙度计算方法
CN116342974A (zh) * 2023-03-20 2023-06-27 北京医准智能科技有限公司 模型训练方法、膝关节分割方法、装置及设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102141545A (zh) * 2010-11-27 2011-08-03 江西理工大学 基于爆破地震波时空衰减规律的岩体力学参数测试方法
CN103926156A (zh) * 2014-04-10 2014-07-16 内蒙古科技大学 一种三维岩石结构面剪切力学行为多重分形评估方法
CN104697436A (zh) * 2015-02-15 2015-06-10 中国科学院西安光学精密机械研究所 一种基于傅里叶级数的圆感应同步器误差模型分析方法
CN104881564A (zh) * 2015-03-09 2015-09-02 绍兴文理学院 结构面粗糙度系数尺寸效应概率密度函数模型的构建方法
CN105205035A (zh) * 2015-10-13 2015-12-30 哈尔滨工程大学 一种非均匀弹性约束边界条件矩形板结构面内振动分析方法
CN105737768A (zh) * 2016-02-01 2016-07-06 绍兴文理学院 基于Cosine相似度量的结构面粗糙度系数评价方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6915225B2 (en) * 2003-05-15 2005-07-05 Northrop Grumman Corporation Method, apparatus and system for digital data resampling utilizing fourier series based interpolation
ATE556646T1 (de) * 2004-01-13 2012-05-15 Univ Toledo Noninvasive, doppelbrechungskompensierte messpolarimeter
DE102010021672A1 (de) * 2010-05-27 2011-12-01 Deutsches Zentrum für Luft- und Raumfahrt e.V. Verfahren zur Ermittlung einer Wirbelgeometrieänderung

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102141545A (zh) * 2010-11-27 2011-08-03 江西理工大学 基于爆破地震波时空衰减规律的岩体力学参数测试方法
CN103926156A (zh) * 2014-04-10 2014-07-16 内蒙古科技大学 一种三维岩石结构面剪切力学行为多重分形评估方法
CN104697436A (zh) * 2015-02-15 2015-06-10 中国科学院西安光学精密机械研究所 一种基于傅里叶级数的圆感应同步器误差模型分析方法
CN104881564A (zh) * 2015-03-09 2015-09-02 绍兴文理学院 结构面粗糙度系数尺寸效应概率密度函数模型的构建方法
CN105205035A (zh) * 2015-10-13 2015-12-30 哈尔滨工程大学 一种非均匀弹性约束边界条件矩形板结构面内振动分析方法
CN105737768A (zh) * 2016-02-01 2016-07-06 绍兴文理学院 基于Cosine相似度量的结构面粗糙度系数评价方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
APPLICATION OF DISCRETE FOURIER SERIES TO THE FINITE ELEMENT STRESS ANALYSIS OF AXI-SYMMETRIC SOLIDS;J.Y.LA1等;《INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN ENGINEERING》;19991231;第619-647页 *
不同加载方式下岩石断裂面形貌的多尺度分析;潘学哉;《中国博士学位论文全文数据库（工程科技Ⅱ辑）》;20140815(第08期);第C038-16页 *

Also Published As

Publication number Publication date
CN106570338A (zh) 2017-04-19

Similar Documents

Publication Publication Date Title
CN106570338B (zh) 岩体结构面粗糙度尺寸效应中轮廓线采样精度确定方法
CN107656902B (zh) 不同采样间距下的结构面粗糙度系数统计方法
CN105066912B (zh) 酸刻蚀物理模拟实验中岩板表面扫描数据的步长标定方法
CN105571511B (zh) 一种船舶外板成型精度在线检测方法
CN104776817A (zh) 滑动式瓷砖表面平整度检测装置及其方法
CN107563087B (zh) 最优采样间距条件下的结构面粗糙度系数统计方法
CN105783842B (zh) 一种大螺距外螺纹加工表面形貌分布特性的检测方法
CN104260074A (zh) 划线检测装置
CN106482674A (zh) 基于中智数函数的结构面粗糙度尺寸效应的近似表达方法
CN104089599B (zh) 用于接触式测头测量中提取二维轮廓的准形态学滤波方法
CN107330142A (zh) 一种基于stl模型的在机检测的测点法矢估算方法
Frącz et al. Aspects of verification and optimization of sheet metal numerical simulations process using the photogrammetric system
CN107036561A (zh) 基于中智数函数的结构面粗糙度各向异性的近似表达方法
CN206291836U (zh) 一种专用角度样板
CN110826280A (zh) 基于有限元模拟的筒形件拉深凸耳改善的工艺优化方法
CN204036447U (zh) 划线检测装置
CN202853515U (zh) 用于测量u型钢板桩锁口尺寸的游标卡尺
CN109615657A (zh) 基于点云数据的含螺纹目标物位姿的计算方法
CN205482655U (zh) 用于测量倒角的游标卡尺
CN106546189A (zh) 一种系列尺寸岩体结构面起伏幅度各向异性的评价方法
CN106441155B (zh) 一种结构面轮廓线采样精度的确定方法
CN105352402B (zh) 斜面内环槽宽度测量量具
CN209295833U (zh) 一种汽车零部件曲面位置度的检测工装
CN106441032A (zh) 一种超薄浮法玻璃量化弯曲度值的检测方法
CN107895061B (zh) 共轨喷油器平面密封性评价用粗糙度参数的选定方法

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: He Zhihai

Inventor before: Yong Rui

Inventor before: Li Bo

Inventor before: Ye Jun

Inventor before: Wu Faquan

Inventor before: Xia Caichu

Inventor before: Liang Qifeng

Inventor before: Huang Man

Inventor before: Ma Chengrong

GR01 Patent grant
GR01 Patent grant