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 PDF

Info

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
Application number
CN201510398349.XA
Other languages
Chinese (zh)
Other versions
CN105004662B (en
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.)
Institute of Mechanics of CAS
Original Assignee
Institute of Mechanics of CAS
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 Institute of Mechanics of CAS filed Critical Institute of Mechanics of CAS
Priority to CN201510398349.XA priority Critical patent/CN105004662B/en
Publication of CN105004662A publication Critical patent/CN105004662A/en
Application granted granted Critical
Publication of CN105004662B publication Critical patent/CN105004662B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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

一种测试岩体结构面接触刚度的方法及装置A method and device for testing the contact stiffness of a rock mass structural surface

技术领域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):

νν == 11 -- cc pp 22 // (( 22 cc pp 22 -- 22 cc sthe s 22 )) -- -- -- (( 11 ))

EE. == 22 ρρ (( 11 ++ νν )) cc sthe s 22 -- -- -- (( 22 ))

其中ν为泊松比、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)

1. test a method for rock mass discontinuity contact stiffness, it is characterized in that: the method comprises the following steps:
(1) choose good structural plane rock mass of appearing as test zone, remove dust and the loose crushing body on test zone surface;
(2) be that symcenter installs vibration transducer in both sides with rock mass discontinuity, each vibration transducer keeps point-blank;
(3) utilize data line to be connected with Acquisition Instrument by vibration transducer, open Acquisition Instrument and make it be in sample states;
(4) carry out hammering in the side of rock mass discontinuity, and produced the compressional wave and shear wave propagated along vibration transducer line direction by adjustment hammering direction;
(5) vibration signal utilizing each vibration transducer perception hammering to produce, and carry out record by Acquisition Instrument;
(6) according to Induction Peried and the distance of each vibration transducer, the longitudinal wave velocity of rockmass, transverse wave speed, compressional wave is calculated by consuming time by structural plane of consuming time, the shear wave of structural plane;
(7) density of the rockmass at test zone place is tested, and calculate elastic modulus and the Poisson ratio of rockmass;
(8) carry out back analysis by numerical method, set up the numerical model similar to test zone, the density of rockmass, elastic modulus and Poisson ratio all get field measurement parameter;
(9) shock load is applied in the side of numerical model, difference consuming time by structural plane of stress wave is obtained by the normal contact stiffness that constantly adjusts structural plane in numerical model and tangential contact stiffness, when stress wave is by the consuming time and step (6) of structural plane consuming time consistent, normal contact stiffness now and tangential contact stiffness are the contact stiffness of test zone structural plane.
2. the method for test rock mass discontinuity contact stiffness according to claim 1, is characterized in that: described structural plane comprises the dry structure face without thickness and the structural plane containing weak intercalated layer.
3. the method for test rock mass discontinuity contact stiffness according to claim 2, it is characterized in that: the vibration transducer quantity of the every side of described structural plane is more than or equal to 2, structural plane is more than or equal to 50cm with the spacing of side senser, and structural plane both sides are less than or equal to 10cm from the nearest sensor of structural plane to the distance of structural plane.
4. the method for test rock mass discontinuity contact stiffness according to claim 3, it is characterized in that: described vibration transducer comprises acceleration transducer, speed pickup, the frequency response of vibration transducer is more than or equal to 1KHz.
5. the method for test rock mass discontinuity contact stiffness according to claim 4, is characterized in that: described vibration transducer is bonding by gypsum, accelerated cement is bonding, and chemical glue mode that is bonding or expansion bolt connection is installed to rock mass.
6. the method for test rock mass discontinuity contact stiffness according to claim 5, is characterized in that: in described step (4), produce compressional wave when hammering direction is vertical with structural plane, produce shear wave when hammering direction is parallel with structural plane.
7. the method for test rock mass discontinuity contact stiffness according to claim 6, it is characterized in that: in described step (7), the Poisson ratio of rockmass is calculated, according to formula (2) calculating elastic modulus according to formula (1):
ν = 1 - c p 2 / ( 2 c p 2 - 2 c s 2 ) - - - ( 1 )
E = 2 ρ ( 1 + v ) c s 2 - - - ( 2 )
Wherein ν is Poisson ratio, E is elastic modulus, ρ is density, c pfor longitudinal wave velocity, c sfor transverse wave speed.
8. the method for test rock mass discontinuity contact stiffness according to claim 7, is characterized in that: the numerical method in described step (8) comprises finite element method, distinct element method, finite volume method, method of finite difference, gridless routing.
9. the method for test rock mass discontinuity contact stiffness according to claim 8, is characterized in that: in the adjustment numerical model in described step (9), the normal contact stiffness of structural plane and the method for tangential contact stiffness comprise dichotomy, line of steepest descent method, method of conjugate gradient, simulated annealing.
10. use a device for the method for test rock mass discontinuity contact stiffness according to claim 9, it is characterized in that: it comprises fullering tool, some vibration transducers, Acquisition Instrument, data processing unit;
Vibration transducer configuration is that symcenter is installed in both sides with rock mass discontinuity, keeps point-blank, the vibration signal that perception hammering produces;
Acquisition Instrument configuration utilizes data line to be connected with vibration transducer, the vibration signal that record hammering produces;
Fullering tool configuration carries out hammering in the side of rock mass discontinuity, and is produced the compressional wave and shear wave propagated along vibration transducer line direction by adjustment hammering direction;
Data processing unit configuration calculates the longitudinal wave velocity of rockmass, transverse wave speed, compressional wave by consuming time by structural plane of consuming time, the shear wave of structural plane; Calculate elastic modulus and the Poisson ratio of rockmass; Carry out back analysis by numerical method, set up the numerical model similar to test zone; The normal contact stiffness of structural plane and tangential contact stiffness in continuous adjustment numerical model, when stress wave by the consuming time and compressional wave of structural plane by consuming time, the shear wave of structural plane by structural plane consuming time consistent time, normal contact stiffness now and tangential contact stiffness are the contact stiffness of test zone structural plane.
CN201510398349.XA 2015-07-08 2015-07-08 A kind of method and device for testing rock mass discontinuity contact stiffness Active CN105004662B (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (6)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
杨德生等: ""剪切波速确定成层土层综合地基刚度系数的研究"", 《工程勘察》 *
王观石等: ""波速和衰减对岩石力学参数和岩体结构变化的敏感性研究》", 《现代矿业》 *
陈卫忠等: ""节理岩体中洞室围岩大变形数值模拟及模型试验研究"", 《岩石力学与工程学报》 *

Cited By (10)

* Cited by examiner, † Cited by third party
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