CN105004662A - Method for testing contact rigidity of rock discontinuity structural plane, and apparatus thereof - Google Patents
Method for testing contact rigidity of rock discontinuity structural plane, and apparatus thereof Download PDFInfo
- Publication number
- CN105004662A CN105004662A CN201510398349.XA CN201510398349A CN105004662A CN 105004662 A CN105004662 A CN 105004662A CN 201510398349 A CN201510398349 A CN 201510398349A CN 105004662 A CN105004662 A CN 105004662A
- Authority
- CN
- China
- Prior art keywords
- contact stiffness
- rock mass
- structural plane
- structural surface
- test
- 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
- 239000011435 rock Substances 0.000 title claims abstract description 77
- 238000012360 testing method Methods 0.000 title claims abstract description 64
- 238000000034 method Methods 0.000 title claims abstract description 54
- 230000001133 acceleration Effects 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 7
- 239000000428 dust Substances 0.000 claims description 4
- 238000011439 discrete element method Methods 0.000 claims description 3
- 239000010440 gypsum Substances 0.000 claims description 3
- 229910052602 gypsum Inorganic materials 0.000 claims description 3
- 239000004568 cement Substances 0.000 claims description 2
- 230000004044 response Effects 0.000 claims description 2
- 238000002922 simulated annealing Methods 0.000 claims description 2
- 238000002945 steepest descent method Methods 0.000 claims description 2
- 239000000126 substance Substances 0.000 claims description 2
- 230000008447 perception Effects 0.000 claims 2
- 230000000644 propagated effect Effects 0.000 claims 2
- 239000003292 glue Substances 0.000 claims 1
- 230000006698 induction Effects 0.000 claims 1
- 238000005259 measurement Methods 0.000 claims 1
- 230000035939 shock Effects 0.000 claims 1
- 230000001902 propagating effect Effects 0.000 abstract description 11
- 238000005070 sampling Methods 0.000 abstract description 6
- 238000002474 experimental method Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 3
- 238000012956 testing procedure Methods 0.000 description 3
- 238000009412 basement excavation Methods 0.000 description 2
- 238000002939 conjugate gradient method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000010438 granite Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000010008 shearing Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- 235000019738 Limestone Nutrition 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 238000004026 adhesive bonding Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000003776 cleavage reaction Methods 0.000 description 1
- 238000012669 compression test Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 239000011229 interlayer Substances 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 230000007017 scission Effects 0.000 description 1
Landscapes
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
本发明公开了一种测试岩体结构面接触刚度的方法,包括步骤:(1)选取测试区域;(2)安装振动传感器;(3)打开采集仪使其处于采样状态;(4)在岩体结构面的一侧进行锤击,并通过调整锤击方向产生沿着振动传感器连线方向传播的纵波及横波;(5)通过采集仪进行记录;(6)计算纵波波速、横波波速、纵波通过结构面的耗时、横波通过结构面的耗时;(7)计算完整岩体的弹性模量及泊松比;(8)借助数值方法进行反演分析,建立与测试区域相似的数值模型;(9)在数值模型的一侧施加冲击载荷,通过不断调整数值模型中结构面的法向接触刚度和切向接触刚度,获得测试区域结构面的接触刚度。还提供了使用该方法的装置。
The invention discloses a method for testing the contact stiffness of a rock mass structural surface, which comprises the steps of: (1) selecting a test area; (2) installing a vibration sensor; (3) turning on the acquisition instrument so that it is in a sampling state; One side of the structural surface of the body is hammered, and by adjusting the hammering direction, the longitudinal wave and shear wave propagating along the direction of the vibration sensor connection are generated; (5) record through the acquisition instrument; (6) calculate the longitudinal wave velocity, shear wave velocity, longitudinal wave Time-consuming through the structural surface and shear wave passing through the structural surface; (7) Calculate the elastic modulus and Poisson's ratio of the complete rock mass; (8) Perform inversion analysis with the help of numerical methods, and establish a numerical model similar to the test area (9) Apply an impact load on one side of the numerical model, and obtain the contact stiffness of the structural surface in the test area by continuously adjusting the normal contact stiffness and tangential contact stiffness of the structural surface in the numerical model. Devices using the method are also provided.
Description
技术领域technical field
本发明属于地质学和力学的技术领域,具体地涉及一种测试岩体结构面接触刚度的方法,以及使用该方法的装置。The invention belongs to the technical field of geology and mechanics, and in particular relates to a method for testing the contact stiffness of a rock mass structural surface and a device using the method.
背景技术Background technique
地质体中存在大量断层、节理、劈理等结构面,上述结构面的力学性质将直接影响地质体的稳定程度及潜在的失稳破坏模式。结构面的力学性质包括结构面法向接触刚度、切向接触刚度、粘聚力、内摩擦角及抗拉强度等。在上述结构面的力学性质中,结构面的接触刚度不仅影响地质体的静态应力应变关系,更直接影响地质体的动态力学行为,如应力波在地质体中的传播规律、地质体在动载荷作用下的破坏演化规律等。There are a large number of structural planes such as faults, joints, and cleavages in the geological body. The mechanical properties of the above structural planes will directly affect the stability of the geological body and the potential failure mode of instability. The mechanical properties of the structural surface include the normal contact stiffness, tangential contact stiffness, cohesion, internal friction angle and tensile strength of the structural surface. Among the above-mentioned mechanical properties of the structural surface, the contact stiffness of the structural surface not only affects the static stress-strain relationship of the geological body, but also directly affects the dynamic mechanical behavior of the geological body, such as the propagation law of the stress wave in the geological body, the dynamic load of the geological body The law of damage evolution under the action of
目前,对于结构面粘聚力、内摩擦角及抗拉强度的测试及实验技术已经较为成熟,如原位剪切实验,室内直剪实验、室内三轴实验及室内巴西劈裂实验等。但关于结构面接触刚度的测试方法及分析技术较少,一般采用单轴压缩实验完成法向接触刚度的测试,采用直剪实验完成切向接触刚度的测试。法向接触刚度测试时,首先测试完整岩块的法向刚度,而后测试含结构面岩体的法向刚度(测试时结构面保持水平状态),最后根据相关公式计算出结构面的法向接触刚度。切向接触刚度测试时,将结构面作为直剪面,通过对结构面上下的完整部分施加切向载荷,观察结构面处的切向位移随切向力的变化规律,进而计算出结构面的切向接触刚度。At present, the testing and experimental techniques for structural surface cohesion, internal friction angle and tensile strength are relatively mature, such as in-situ shearing experiment, indoor direct shearing experiment, indoor triaxial experiment and indoor Brazilian splitting experiment, etc. However, there are few test methods and analysis techniques for the contact stiffness of the structural surface. Generally, the uniaxial compression test is used to complete the test of the normal contact stiffness, and the direct shear test is used to complete the test of the tangential contact stiffness. When testing the normal contact stiffness, first test the normal stiffness of the complete rock block, then test the normal stiffness of the rock mass with structural plane (the structural plane remains horizontal during the test), and finally calculate the normal contact of the structural plane according to the relevant formula stiffness. When testing the tangential contact stiffness, the structural surface is used as a direct shear surface, and by applying a tangential load to the complete part above and below the structural surface, observing the change law of the tangential displacement at the structural surface with the tangential force, and then calculating the structural surface Tangential contact stiffness.
上述测试结构面接触刚度的方法只能在实验室完成,因此需要进行现场取样及试件加工,工序较为复杂、对结构面扰动较大、且无法反映原岩应力对结构面刚度的影响。此外,受实验设备的限制,待测试的试样及结构面尺寸一般在分米量级,无法测量大型结构面的接触刚度。最后,上述方法测试的刚度值为结构面的静态刚度,当把此测试的刚度值用于动态问题分析时,将存在较大的误差。The above method of testing the contact stiffness of the structural surface can only be done in the laboratory, so on-site sampling and specimen processing are required, the process is relatively complicated, the structural surface is greatly disturbed, and the influence of the original rock stress on the structural surface stiffness cannot be reflected. In addition, due to the limitation of experimental equipment, the size of the sample and structural surface to be tested is generally on the order of decimeters, and the contact stiffness of large structural surfaces cannot be measured. Finally, the stiffness value tested by the above method is the static stiffness of the structural surface. When the stiffness value of this test is used for dynamic problem analysis, there will be a large error.
发明内容Contents of the invention
本发明的技术解决问题是:克服现有技术的不足,提供一种测试岩体结构面接触刚度的方法,其能够在野外直接对结构面岩体实施测量分析,测试精度高,测试工序简单,能反映原岩应力对结构面刚度的影响,且可以测量大型结构面的接触刚度。The technical solution of the present invention is to overcome the deficiencies of the prior art and provide a method for testing the contact stiffness of the rock mass structural surface, which can directly measure and analyze the structural surface rock mass in the field, with high testing accuracy and simple testing procedures. It can reflect the influence of original rock stress on structural surface stiffness, and can measure the contact stiffness of large structural surfaces.
本发明的技术解决方案是:这种测试岩体结构面接触刚度的方法,该方法包括以下步骤:The technical solution of the present invention is: the method for this test rock mass structural surface contact stiffness, this method comprises the following steps:
(1)选取露头较好的结构面岩体作为测试区域,清除测试区域表面的粉尘及松散破碎体;(1) Select the rock mass with better structural surface as the test area, and remove the dust and loose broken bodies on the surface of the test area;
(2)以岩体结构面为对称中心在两侧安装振动传感器,各振动传感器保持在一条直线上;(2) Vibration sensors are installed on both sides with the rock mass structural surface as the symmetrical center, and each vibration sensor is kept on a straight line;
(3)利用数据线将振动传感器与采集仪相连,打开采集仪使其处于采样状态;(3) Utilize the data line to connect the vibration sensor with the acquisition instrument, open the acquisition instrument to make it in the sampling state;
(4)在岩体结构面的一侧进行锤击,并通过调整锤击方向产生沿着振动传感器连线方向传播的纵波及横波;(4) Hammering is performed on one side of the structural surface of the rock mass, and by adjusting the direction of the hammering, the longitudinal wave and the transverse wave propagating along the direction of the vibration sensor connection line are generated;
(5)利用各振动传感器感知锤击产生的振动信号,并通过采集仪进行记录;(5) Utilize each vibration sensor to perceive the vibration signal generated by hammering, and record it through the acquisition instrument;
(6)根据各振动传感器的起振时间和距离,计算出完整岩体的纵波波速、横波波速、纵波通过结构面的耗时、横波通过结构面的耗时;(6) According to the start-up time and distance of each vibration sensor, calculate the longitudinal wave velocity, shear wave velocity, time consumption of longitudinal wave passing through the structural surface, and time consumption of shear wave passing through the structural surface of the complete rock mass;
(7)对测试区域所在的完整岩体的密度进行测试,并计算出完整岩体的弹性模量及泊松比;(7) Test the density of the complete rock mass where the test area is located, and calculate the elastic modulus and Poisson's ratio of the complete rock mass;
(8)借助数值方法进行反演分析,建立与测试区域相似的数值模型,完整岩体的密度、弹性模量及泊松比均取现场实测参数;(8) Inversion analysis is carried out by means of numerical methods, and a numerical model similar to that of the test area is established. The density, elastic modulus and Poisson's ratio of the complete rock mass are all measured on site;
(9)在数值模型的一侧施加冲击载荷,通过不断调整数值模型中结构面的法向接触刚度和切向接触刚度获得应力波通过结构面的不同耗时,当应力波通过结构面的耗时与步骤(6)的耗时一致时,此时的法向接触刚度和切向接触刚度为测试区域结构面的接触刚度。(9) Apply an impact load on one side of the numerical model, and obtain the different time consumption of the stress wave passing through the structural surface by continuously adjusting the normal contact stiffness and tangential contact stiffness of the structural surface in the numerical model. When the time consumption is the same as that of step (6), the normal contact stiffness and tangential contact stiffness at this time are the contact stiffness of the structural surface in the test area.
还提供了一种使用这种测试岩体结构面接触刚度的方法的装置,其包括锤击工具、若干振动传感器、采集仪、数据处理单元;Also provided is a device using the method for testing the contact stiffness of a rock mass structural surface, which includes a hammering tool, several vibration sensors, an acquisition instrument, and a data processing unit;
振动传感器配置来以岩体结构面为对称中心在两侧安装,保持在一条直线上,感知锤击产生的振动信号;The vibration sensor is configured to install on both sides with the rock mass structure plane as the symmetrical center, keep in a straight line, and sense the vibration signal generated by the hammering;
采集仪配置来利用数据线与振动传感器相连,记录锤击产生的振动信号;The acquisition instrument is configured to use the data line to connect with the vibration sensor to record the vibration signal generated by the hammer;
锤击工具配置来在岩体结构面的一侧进行锤击,并通过调整锤击方向产生沿着振动传感器连线方向传播的纵波及横波;The hammering tool is configured to hammer on one side of the structural surface of the rock mass, and by adjusting the hammering direction, the longitudinal wave and the transverse wave propagating along the direction of the vibration sensor connection line are generated;
数据处理单元配置来计算完整岩体的纵波波速、横波波速、纵波通过结构面的耗时、横波通过结构面的耗时;计算完整岩体的弹性模量及泊松比;借助数值方法进行反演分析,建立与测试区域相似的数值模型;不断调整数值模型中结构面的法向接触刚度和切向接触刚度,当应力波通过结构面的耗时与纵波通过结构面的耗时、横波通过结构面的耗时一致时,此时的法向接触刚度和切向接触刚度为测试区域结构面的接触刚度。The data processing unit is configured to calculate the longitudinal wave velocity of the complete rock mass, the shear wave velocity, the time consumption of the longitudinal wave passing through the structural surface, and the time consumption of the shear wave passing through the structural surface; calculate the elastic modulus and Poisson's ratio of the complete rock mass; use the numerical method to inverse Perform analysis and establish a numerical model similar to the test area; constantly adjust the normal contact stiffness and tangential contact stiffness of the structural surface in the numerical model. When the time consumption of the structural surface is consistent, the normal contact stiffness and tangential contact stiffness at this time are the contact stiffness of the structural surface in the test area.
本发明通过在结构面一侧分别按照垂直于结构面及平行于结构面两种方向锤击岩体,使其分别产生沿着传感器连线方向传播的纵波及横波,借助采集仪捕捉各传感器感知到的起振时间,根据各传感器间的距离计算完整岩体的纵波波速、横波波速、纵波通过结构面的耗时及横波通过结构面的耗时,根据岩体的密度、纵波波速及横波波速计算完整岩体的弹性模量及泊松比,利用数值计算方法对纵波及横波通过结构面的耗时进行分析计算,通过不断调整数值模型中结构面的法向接触刚度及切向接触刚度,使得数值计算获得的纵波耗时及横波耗时与已计算的耗时一致,此法向接触刚度及切向接触刚度即为测试区域结构面的接触刚度,因此这种测试岩体结构面接触刚度的方法能够在野外直接对结构面岩体实施测量分析,测试精度高,测试工序简单,能反映原岩应力对结构面刚度的影响,且可以测量大型结构面的接触刚度。The present invention hammers the rock mass in two directions perpendicular to the structural surface and parallel to the structural surface on one side of the structural surface, so that the longitudinal waves and transverse waves propagating along the direction of the sensor connection line are respectively generated, and the sensors sensed by the acquisition instrument are captured. According to the distance between the sensors, calculate the longitudinal wave velocity, shear wave velocity, time consumption of longitudinal wave passing through the structural surface and time consumption of shear wave passing through the structural surface of the complete rock mass, according to the density of rock mass, longitudinal wave velocity and shear wave velocity Calculate the elastic modulus and Poisson's ratio of the complete rock mass, and use the numerical calculation method to analyze and calculate the time-consuming time for longitudinal waves and shear waves to pass through the structural surface. By continuously adjusting the normal contact stiffness and tangential contact stiffness of the structural surface in the numerical model, The longitudinal wave time consumption and shear wave time consumption obtained by numerical calculation are consistent with the calculated time consumption. The normal contact stiffness and tangential contact stiffness are the contact stiffness of the structural surface in the test area. Therefore, the contact stiffness of the test rock mass structural surface The method can directly measure and analyze structural surface rock mass in the field, has high test accuracy, simple testing procedure, can reflect the influence of original rock stress on structural surface stiffness, and can measure the contact stiffness of large structural surfaces.
附图说明Description of drawings
图1示出了根据本发明的单结构面锤击实验系统布设图。Fig. 1 shows the layout diagram of the single-structure surface hammering experiment system according to the present invention.
图2示出了根据本发明的双结构面锤击实验系统布设图。Fig. 2 shows the layout diagram of the double-structure surface hammering experiment system according to the present invention.
图3示出了根据本发明的测试岩体结构面接触刚度的方法的流程图。Fig. 3 shows a flowchart of a method for testing the contact stiffness of a rock mass structural surface according to the present invention.
具体实施方式Detailed ways
从图3中可以看出,这种测试岩体结构面接触刚度的方法,该方法包括以下步骤:As can be seen from Fig. 3, the method for testing the contact stiffness of the rock mass structural surface comprises the following steps:
(1)选取露头较好的结构面岩体作为测试区域,清除(利用毛刷等)测试区域表面的粉尘及松散破碎体;(1) Select the rock mass with better structural surface as the test area, and remove (using a brush, etc.) the dust and loose broken bodies on the surface of the test area;
(2)以岩体结构面为对称中心在两侧安装振动传感器,各振动传感器保持在一条直线上;(2) Vibration sensors are installed on both sides with the rock mass structural surface as the symmetrical center, and each vibration sensor is kept on a straight line;
(3)利用数据线将振动传感器与采集仪相连,打开采集仪使其处于采样状态;(3) Utilize the data line to connect the vibration sensor with the acquisition instrument, open the acquisition instrument to make it in the sampling state;
(4)在岩体结构面的一侧进行锤击,并通过调整锤击方向产生沿着振动传感器连线方向传播的纵波及横波;(4) Hammering is performed on one side of the structural surface of the rock mass, and by adjusting the direction of the hammering, the longitudinal wave and the transverse wave propagating along the direction of the vibration sensor connection line are generated;
(5)利用各振动传感器感知锤击产生的振动信号,并通过采集仪进行记录;(5) Utilize each vibration sensor to perceive the vibration signal generated by hammering, and record it through the acquisition instrument;
(6)根据各振动传感器的起振时间和距离,计算出完整岩体的纵波波速、横波波速、纵波通过结构面的耗时、横波通过结构面的耗时;(6) According to the start-up time and distance of each vibration sensor, calculate the longitudinal wave velocity, shear wave velocity, time consumption of longitudinal wave passing through the structural surface, and time consumption of shear wave passing through the structural surface of the complete rock mass;
(7)对测试区域所在的完整岩体的密度进行测试,并计算出完整岩体的弹性模量及泊松比;(7) Test the density of the complete rock mass where the test area is located, and calculate the elastic modulus and Poisson's ratio of the complete rock mass;
(8)借助数值方法进行反演分析,建立与测试区域相似的数值模型,完整岩体的密度、弹性模量及泊松比均取现场实测参数;(8) Inversion analysis is carried out by means of numerical methods, and a numerical model similar to that of the test area is established. The density, elastic modulus and Poisson's ratio of the complete rock mass are all measured on site;
(9)在数值模型的一侧施加冲击载荷,通过不断调整数值模型中结构面的法向接触刚度和切向接触刚度获得应力波通过结构面的不同耗时,当应力波通过结构面的耗时与步骤(6)的耗时一致时,此时的法向接触刚度和切向接触刚度为测试区域结构面的接触刚度。(9) Apply an impact load on one side of the numerical model, and obtain the different time consumption of the stress wave passing through the structural surface by continuously adjusting the normal contact stiffness and tangential contact stiffness of the structural surface in the numerical model. When the time consumption is the same as that of step (6), the normal contact stiffness and tangential contact stiffness at this time are the contact stiffness of the structural surface in the test area.
另外,所述结构面包括无厚度的干性结构面和含软弱夹层的结构面。In addition, the structural surface includes a dry structural surface without thickness and a structural surface with a weak interlayer.
另外,所述结构面每侧的振动传感器数量大于等于2支,结构面同侧传感器的间距大于等于50cm,结构面两侧离结构面最近的传感器到结构面的距离小于等于10cm。In addition, the number of vibration sensors on each side of the structural surface is greater than or equal to 2, the distance between the sensors on the same side of the structural surface is greater than or equal to 50cm, and the distance from the sensors closest to the structural surface on both sides of the structural surface to the structural surface is less than or equal to 10cm.
另外,所述振动传感器包括加速度传感器、速度传感器,振动传感器的频响大于等于1KHz。In addition, the vibration sensor includes an acceleration sensor and a speed sensor, and the frequency response of the vibration sensor is greater than or equal to 1KHz.
另外,所述振动传感器通过石膏粘接、速凝水泥粘接,化学胶水粘接或膨胀螺栓连接的方式安装到岩体。In addition, the vibration sensor is installed to the rock mass through gypsum bonding, quick-setting cement bonding, chemical glue bonding or expansion bolt connection.
另外,所述步骤(4)中,锤击方向与结构面垂直时产生纵波,锤击方向与结构面平行时产生横波。In addition, in the step (4), longitudinal waves are generated when the hammering direction is perpendicular to the structural surface, and transverse waves are generated when the hammering direction is parallel to the structural surface.
另外,所述步骤(7)中,根据公式(1)计算完整岩体的泊松比,根据公式(2)计算弹性模量:In addition, in the step (7), the Poisson's ratio of the complete rock mass is calculated according to the formula (1), and the modulus of elasticity is calculated according to the formula (2):
其中ν为泊松比、E为弹性模量、ρ为密度、cp为纵波波速、cs为横波Among them, ν is Poisson's ratio, E is elastic modulus, ρ is density, c p is longitudinal wave velocity, c s is transverse wave
波速。wave speed.
另外,所述步骤(8)中的数值方法包括有限元法、离散元法、有限体积法、有限差分法、无网格法。In addition, the numerical methods in the step (8) include finite element method, discrete element method, finite volume method, finite difference method, and meshless method.
另外,所述步骤(9)中的调整数值模型中结构面的法向接触刚度和切向接触刚度的方法包括二分法、最速下降线法、共轭梯度法、模拟退火算法。In addition, the method for adjusting the normal contact stiffness and tangential contact stiffness of the structural surface in the numerical model in the step (9) includes the dichotomy method, the steepest descent method, the conjugate gradient method, and the simulated annealing algorithm.
还提供了一种使用这种测试岩体结构面接触刚度的方法的装置,其包括锤击工具、若干振动传感器、采集仪、数据处理单元;Also provided is a device using the method for testing the contact stiffness of a rock mass structural surface, which includes a hammering tool, several vibration sensors, an acquisition instrument, and a data processing unit;
振动传感器配置来以岩体结构面为对称中心在两侧安装,保持在一条直线上,感知锤击产生的振动信号;The vibration sensor is configured to install on both sides with the rock mass structure plane as the symmetrical center, keep in a straight line, and sense the vibration signal generated by the hammering;
采集仪配置来利用数据线与振动传感器相连,记录锤击产生的振动信号;The acquisition instrument is configured to use the data line to connect with the vibration sensor to record the vibration signal generated by the hammer;
锤击工具配置来在岩体结构面的一侧进行锤击,并通过调整锤击方向产生沿着振动传感器连线方向传播的纵波及横波;The hammering tool is configured to hammer on one side of the structural surface of the rock mass, and by adjusting the hammering direction, the longitudinal wave and the transverse wave propagating along the direction of the vibration sensor connection line are generated;
数据处理单元配置来计算完整岩体的纵波波速、横波波速、纵波通过结构面的耗时、横波通过结构面的耗时;计算完整岩体的弹性模量及泊松比;借助数值方法进行反演分析,建立与测试区域相似的数值模型;不断调整数值模型中结构面的法向接触刚度和切向接触刚度,当应力波通过结构面的耗时与纵波通过结构面的耗时、横波通过结构面的耗时一致时,此时的法向接触刚度和切向接触刚度为测试区域结构面的接触刚度。The data processing unit is configured to calculate the longitudinal wave velocity of the complete rock mass, the shear wave velocity, the time consumption of the longitudinal wave passing through the structural surface, and the time consumption of the shear wave passing through the structural surface; calculate the elastic modulus and Poisson's ratio of the complete rock mass; use the numerical method to inverse Perform analysis and establish a numerical model similar to the test area; constantly adjust the normal contact stiffness and tangential contact stiffness of the structural surface in the numerical model. When the time consumption of the structural surface is consistent, the normal contact stiffness and tangential contact stiffness at this time are the contact stiffness of the structural surface in the test area.
本发明通过在结构面一侧分别按照垂直于结构面及平行于结构面两种方向锤击岩体,使其分别产生沿着传感器连线方向传播的纵波及横波,借助采集仪捕捉各传感器感知到的起振时间,根据各传感器间的距离计算完整岩体的纵波波速、横波波速、纵波通过结构面的耗时及横波通过结构面的耗时,根据岩体的密度、纵波波速及横波波速计算完整岩体的弹性模量及泊松比,利用数值计算方法对纵波及横波通过结构面的耗时进行分析计算,通过不断调整数值模型中结构面的法向接触刚度及切向接触刚度,使得数值计算获得的纵波耗时及横波耗时与已计算的耗时一致,此法向接触刚度及切向接触刚度即为测试区域结构面的接触刚度,因此这种测试岩体结构面接触刚度的方法能够在野外直接对结构面岩体实施测量分析,测试精度高,测试工序简单,能反映原岩应力对结构面刚度的影响,且可以测量大型结构面的接触刚度。The present invention hammers the rock mass in two directions perpendicular to the structural surface and parallel to the structural surface on one side of the structural surface, so that the longitudinal waves and transverse waves propagating along the direction of the sensor connection line are respectively generated, and the sensors sensed by the acquisition instrument are captured. According to the distance between the sensors, calculate the longitudinal wave velocity, shear wave velocity, time consumption of longitudinal wave passing through the structural surface and time consumption of shear wave passing through the structural surface of the complete rock mass, according to the density of rock mass, longitudinal wave velocity and shear wave velocity Calculate the elastic modulus and Poisson's ratio of the complete rock mass, and use the numerical calculation method to analyze and calculate the time-consuming time for longitudinal waves and shear waves to pass through the structural surface. By continuously adjusting the normal contact stiffness and tangential contact stiffness of the structural surface in the numerical model, The longitudinal wave time consumption and shear wave time consumption obtained by numerical calculation are consistent with the calculated time consumption. The normal contact stiffness and tangential contact stiffness are the contact stiffness of the structural surface in the test area. Therefore, the contact stiffness of the test rock mass structural surface The method can directly measure and analyze structural surface rock mass in the field, has high test accuracy, simple testing procedure, can reflect the influence of original rock stress on structural surface stiffness, and can measure the contact stiffness of large structural surfaces.
以下给出本发明的两个具体实施例。Two specific examples of the present invention are given below.
实施例1:Example 1:
对某露天矿灰岩结构面的法向接触刚度及切向接触刚度进行测试分析,实验过程如图1所示。选取结构面2露头较好的区域,利用毛刷对结构面2及两侧完整岩体1进行清理,清除表面浮尘及松散破碎体。利用石膏,将振动加速度传感器3至6按一定间距固定在结构面2两侧的完整岩体1上,传感器布设时需保证传感器3至6在一条直线上,且垂直于结构面2。加速度传感器4和5到结构面2的距离均为5cm,加速度传感器3和6到结构面的距离均为50cm。利用数据线7将加速度传感器3至6连接至采集仪8,打开采集仪8使其处于采样状态。选取锤击点9及锤击点11两个锤击点,在锤击点9按照锤击方向10进行锤击,产生向结构面2方向传播的纵波;在锤击点11按照锤击方向12进行锤击,产生向结构面2方向传播的横波。利用加速度传感器3至6对振动信号进行拾取,并利用采集仪8对振动信号进行记录。根据传感器3至6的起振时间及传感器3至6间的距离,计算出完整岩体1的纵波波速为3917m/s,横波波速为2165m/s,纵波通过结构面2的耗时为0.21ms,横波通过结构面的耗时为0.65ms。采用天平、量筒等对测试区域岩体1的密度进行测试,为2500kg/m3,由此计算出该结构面2两侧完成岩体1的弹性模量为30GPa,泊松比为0.28。采用可变形块体离散元方法进行数值反演,参数调整策略选用二分法,经过计算分析,结构面2的法向接触刚度为10.1GPa/m,切向接触刚度为3.5GPa/m。The normal contact stiffness and tangential contact stiffness of a limestone structural plane in an open-pit mine were tested and analyzed. The experimental process is shown in Figure 1. Select an area with better outcropping on the structural surface 2, and use a brush to clean the structural surface 2 and the complete rock mass 1 on both sides to remove floating dust and loose broken bodies on the surface. Using gypsum, fix the vibration acceleration sensors 3 to 6 on the complete rock mass 1 on both sides of the structural surface 2 at a certain distance. The distances from acceleration sensors 4 and 5 to the structural surface 2 are both 5 cm, and the distances from acceleration sensors 3 and 6 to the structural surface are both 50 cm. Use the data line 7 to connect the acceleration sensors 3 to 6 to the acquisition instrument 8, and turn on the acquisition instrument 8 to make it in the sampling state. Select two hammering points, hammering point 9 and hammering point 11, and perform hammering according to hammering direction 10 at hammering point 9 to generate longitudinal waves propagating in the direction of structural surface 2; at hammering point 11, follow hammering direction 12 Hammering is performed to generate shear waves propagating in the direction of structural surface 2. Acceleration sensors 3 to 6 are used to pick up the vibration signals, and the acquisition instrument 8 is used to record the vibration signals. According to the start-up time of sensors 3 to 6 and the distance between sensors 3 to 6, it is calculated that the longitudinal wave velocity of the complete rock mass 1 is 3917m/s, the shear wave velocity is 2165m/s, and the time taken for the longitudinal wave to pass through the structural surface 2 is 0.21ms , the shear wave takes 0.65ms to pass through the structural surface. The density of the rock mass 1 in the test area was tested with a balance and a measuring cylinder, and it was 2500kg/m 3 . From this, the elastic modulus of the completed rock mass 1 on both sides of the structural surface 2 was calculated to be 30GPa, and the Poisson's ratio was 0.28. The deformable block discrete element method is used for numerical inversion, and the dichotomy method is selected as the parameter adjustment strategy. After calculation and analysis, the normal contact stiffness of structural surface 2 is 10.1GPa/m, and the tangential contact stiffness is 3.5GPa/m.
实施例2:Example 2:
对某边坡开挖揭露出的花岗岩结构面进行法向接触刚度及切向接触刚度的测试分析,实验过程如图2所示。对开挖揭露出的两个结构面2及周围完整岩体1进行清理,清除表层杂质。利用膨胀螺栓将振动速度传感器3、4、5、6、13、14安装至结构面2周围完整岩体1的对应位置;并确保传感器3、4、5、6、13、14处于一条直线上,且垂直于结构面2。传感器4、传感器5到左侧结构面2的距离均为10cm,传感器6、传感器13到右侧结构面2的距离也为10cm,传感器3与传感器4之间、传感器5与传感器6之间、传感器13与传感器14之间的距离为5m。利用数据线7将传感器3、4、5、6、13、14与采集仪8相连,打开采集仪8使其处于采样状态。在锤击点9按照锤击方向10进行锤击,产生向两个结构面2传播的纵波;在锤击点11按照锤击方向12进行锤击,产生向两个结构面2传播的横波。根据传感器3、4、5、6、13、14记录的振动起跳时间及各传感器的间距计算获得,纵波通过右侧结构面2的耗时为0.33ms,纵波通过左侧结构面2的耗时为0.24ms,横波通过右侧结构面2的耗时为0.78ms,横波通过左侧结构面2的耗时为0.65ms,完整岩块的纵波波速为4307m/s,横波波速为2527m/s。根据现场测试获得的岩体密度2600kg/m3,计算出该花岗岩的弹性模量为41GPa、泊松比为0.24。利用含结构面模型的有限体积法进行数值分析,参数调整方法选为共轭梯度法,经过计算分析,右侧结构面2、左侧结构面2的法向接触刚度分别为9.4GPa/m、15GPa/m,右侧结构面2、左侧结构面2的切向接触刚度分别为2.8GPa/m及6.4GPa/m。The normal contact stiffness and tangential contact stiffness were tested and analyzed for the granite structural surface exposed by the excavation of a certain slope. The experimental process is shown in Figure 2. Clean up the two structural surfaces 2 exposed by the excavation and the surrounding complete rock mass 1 to remove surface impurities. Use expansion bolts to install the vibration velocity sensors 3, 4, 5, 6, 13, 14 to the corresponding positions of the complete rock mass 1 around the structural surface 2; and ensure that the sensors 3, 4, 5, 6, 13, 14 are on a straight line , and perpendicular to the structural plane 2. The distances from sensor 4 and sensor 5 to the left structural surface 2 are both 10 cm, and the distances from sensor 6 and sensor 13 to the right structural surface 2 are also 10 cm. Between sensor 3 and sensor 4, between sensor 5 and sensor 6, The distance between sensor 13 and sensor 14 is 5 m. The sensors 3, 4, 5, 6, 13, 14 are connected to the acquisition instrument 8 by the data line 7, and the acquisition instrument 8 is turned on to make it in the sampling state. Hammering at the hammering point 9 in the hammering direction 10 generates longitudinal waves propagating to the two structural surfaces 2 ; hammering at the hammering point 11 in the hammering direction 12 generates transverse waves propagating to the two structural surfaces 2 . Calculated based on the vibration take-off time recorded by sensors 3, 4, 5, 6, 13, and 14 and the distance between the sensors, the time taken for the longitudinal wave to pass through the structural surface 2 on the right is 0.33ms, and the time for the longitudinal wave to pass through the structural surface 2 on the left It takes 0.24ms for the shear wave to pass through the structural surface 2 on the right, 0.78ms, and 0.65ms for the shear wave to pass through the structural surface 2 on the left. According to the rock mass density of 2600kg/m 3 obtained from the field test, the elastic modulus of the granite is calculated to be 41GPa, and the Poisson's ratio is 0.24. The numerical analysis was carried out by using the finite volume method with structural surface model, and the parameter adjustment method was selected as the conjugate gradient method. After calculation and analysis, the normal contact stiffness of structural surface 2 on the right and structural surface 2 on the left were 9.4GPa/m, 15GPa/m, the tangential contact stiffness of the right structural surface 2 and the left structural surface 2 are 2.8GPa/m and 6.4GPa/m respectively.
以上所述,仅是本发明的较佳实施例,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属本发明技术方案的保护范围。The above are only preferred embodiments of the present invention, and are not intended to limit the present invention in any form. Any simple modifications, equivalent changes and modifications made to the above embodiments according to the technical essence of the present invention are still within the scope of this invention. The protection scope of the technical solution of the invention.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510398349.XA CN105004662B (en) | 2015-07-08 | 2015-07-08 | A kind of method and device for testing rock mass discontinuity contact stiffness |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510398349.XA CN105004662B (en) | 2015-07-08 | 2015-07-08 | A kind of method and device for testing rock mass discontinuity contact stiffness |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105004662A true CN105004662A (en) | 2015-10-28 |
CN105004662B CN105004662B (en) | 2017-10-20 |
Family
ID=54377415
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510398349.XA Active CN105004662B (en) | 2015-07-08 | 2015-07-08 | A kind of method and device for testing rock mass discontinuity contact stiffness |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105004662B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106248324A (en) * | 2016-07-08 | 2016-12-21 | 中国电建集团华东勘测设计研究院有限公司 | The determination method of weak structural face intensive parameter of growing up |
CN106290576A (en) * | 2016-08-12 | 2017-01-04 | 中国电建集团成都勘测设计研究院有限公司 | A kind of Engineering Geology Classification method of compound alterated rocks/band |
CN107817177A (en) * | 2017-09-14 | 2018-03-20 | 吉林大学 | A kind of material blast Deformation Experiments method |
CN109084889A (en) * | 2018-08-30 | 2018-12-25 | 华南农业大学 | A kind of vibration signal detection method |
CN110940587A (en) * | 2019-12-17 | 2020-03-31 | 北京工业大学 | A method for calculating the contact stiffness of bed-foundation interface based on multi-scale theory |
CN112729757A (en) * | 2020-12-04 | 2021-04-30 | 中国航空工业集团公司沈阳空气动力研究所 | Evaluation optimization method for connection end design of chip hinge moment balance |
CN114330065A (en) * | 2021-12-30 | 2022-04-12 | 中国科学院力学研究所 | Numerical algorithm for simulating self-adaptive crack propagation |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4637252A (en) * | 1982-09-23 | 1987-01-20 | The Uniroyal Goodrich Tire Company | Test instrument |
US20030233876A1 (en) * | 2000-06-22 | 2003-12-25 | The Boeing Company | Apparatus for in-situ nondestructive measurement of young's modulus of plate structures |
CN1645101A (en) * | 2005-01-12 | 2005-07-27 | 杜时贵 | Determining method for mechanics quality stable threshold of rock structural face |
US20130121961A1 (en) * | 2011-11-11 | 2013-05-16 | The Johns Hopkins University | Treatment modalities to prevent or treat loss of cardiovascular function in aging humans |
CN103792289A (en) * | 2014-02-14 | 2014-05-14 | 江西理工大学 | Full-waveform information test method for joint rigidity of rock mass |
CN104007462A (en) * | 2014-04-16 | 2014-08-27 | 孙赞东 | Fracture prediction method based on attenuation anisotropy |
-
2015
- 2015-07-08 CN CN201510398349.XA patent/CN105004662B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4637252A (en) * | 1982-09-23 | 1987-01-20 | The Uniroyal Goodrich Tire Company | Test instrument |
US20030233876A1 (en) * | 2000-06-22 | 2003-12-25 | The Boeing Company | Apparatus for in-situ nondestructive measurement of young's modulus of plate structures |
CN1645101A (en) * | 2005-01-12 | 2005-07-27 | 杜时贵 | Determining method for mechanics quality stable threshold of rock structural face |
US20130121961A1 (en) * | 2011-11-11 | 2013-05-16 | The Johns Hopkins University | Treatment modalities to prevent or treat loss of cardiovascular function in aging humans |
CN103792289A (en) * | 2014-02-14 | 2014-05-14 | 江西理工大学 | Full-waveform information test method for joint rigidity of rock mass |
CN104007462A (en) * | 2014-04-16 | 2014-08-27 | 孙赞东 | Fracture prediction method based on attenuation anisotropy |
Non-Patent Citations (3)
Title |
---|
杨德生等: ""剪切波速确定成层土层综合地基刚度系数的研究"", 《工程勘察》 * |
王观石等: ""波速和衰减对岩石力学参数和岩体结构变化的敏感性研究》", 《现代矿业》 * |
陈卫忠等: ""节理岩体中洞室围岩大变形数值模拟及模型试验研究"", 《岩石力学与工程学报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106248324A (en) * | 2016-07-08 | 2016-12-21 | 中国电建集团华东勘测设计研究院有限公司 | The determination method of weak structural face intensive parameter of growing up |
CN106248324B (en) * | 2016-07-08 | 2018-08-31 | 中国电建集团华东勘测设计研究院有限公司 | The determination method for weak structural face intensive parameter of growing up |
CN106290576A (en) * | 2016-08-12 | 2017-01-04 | 中国电建集团成都勘测设计研究院有限公司 | A kind of Engineering Geology Classification method of compound alterated rocks/band |
CN106290576B (en) * | 2016-08-12 | 2018-12-21 | 中国电建集团成都勘测设计研究院有限公司 | A kind of Engineering Geology Classification method of compound alterated rocks/band |
CN107817177A (en) * | 2017-09-14 | 2018-03-20 | 吉林大学 | A kind of material blast Deformation Experiments method |
CN107817177B (en) * | 2017-09-14 | 2020-04-21 | 吉林大学 | An experimental method for material explosion deformation |
CN109084889A (en) * | 2018-08-30 | 2018-12-25 | 华南农业大学 | A kind of vibration signal detection method |
CN110940587A (en) * | 2019-12-17 | 2020-03-31 | 北京工业大学 | A method for calculating the contact stiffness of bed-foundation interface based on multi-scale theory |
CN112729757A (en) * | 2020-12-04 | 2021-04-30 | 中国航空工业集团公司沈阳空气动力研究所 | Evaluation optimization method for connection end design of chip hinge moment balance |
CN114330065A (en) * | 2021-12-30 | 2022-04-12 | 中国科学院力学研究所 | Numerical algorithm for simulating self-adaptive crack propagation |
Also Published As
Publication number | Publication date |
---|---|
CN105004662B (en) | 2017-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105004662B (en) | A kind of method and device for testing rock mass discontinuity contact stiffness | |
Owolabi et al. | Crack detection in beams using changes in frequencies and amplitudes of frequency response functions | |
Agioutantis et al. | Potential of acoustic emissions from three point bending tests as rock failure precursors | |
Gao et al. | Estimation of effect of voids on frequency response of mountain tunnel lining based on microtremor method | |
Iturrioz et al. | Experimental analysis and truss-like discrete element model simulation of concrete specimens under uniaxial compression | |
Shirole et al. | Damage monitoring in rock specimens with pre-existing flaws by non-linear ultrasonic waves and digital image correlation | |
Lee et al. | Evolution of the shear wave velocity during shaking modeled in centrifuge shaking table tests | |
Sang et al. | Quantifying fatigue-damage and failure-precursors using ultrasonic coda wave interferometry | |
Wang et al. | Structural damage detection using local damage factor | |
Gao et al. | Voids delineation behind tunnel lining based on the vibration intensity of microtremors | |
Moayedi et al. | A new real-time monitoring technique in calculation of the py curve of single thin steel piles considering the influence of driven energy and using strain gauge sensors | |
Rathje et al. | Development of an in situ dynamic liquefaction test | |
CN119026403B (en) | Pile foundation load influence assessment method, medium and system by pile sinking soil compaction effect | |
Sebastian et al. | Resonant column tests and nonlinear elasticity in simulated rocks | |
JP2000065809A (en) | Measuring method for damage degree of ground | |
Wang et al. | Dynamic response of gravity dam model with crack and damage detection | |
Jiang et al. | The nature frequency identification of tunnel lining based on the microtremor method | |
Voznesenskii et al. | On the evaluation of rock integrity around mine workings with anchorage by the shock-spectral method | |
Kwasniewski et al. | Self-excited acoustical system for stress measurement in mass rocks | |
Wang et al. | Relationship between local damage and structural dynamic behavior | |
Shokri et al. | Crack source location by acoustic emission monitoring method in RC strips during in-situ load test | |
Liao et al. | Monitoring of a prestressed bridge model by ultrasonic measurement and vibration recordings | |
Saha et al. | Dynamic behaviour of rock joint during shear wave propagation: The influence of joint orientation | |
GR1010172B (en) | Micro-seismic arrangement and analysis method with use of ultrasounds | |
RU2114426C1 (en) | Method of quality inspection of construction structures |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |