CN113866023A - 一种预测岩石杆中应力波大小的方法 - Google Patents
一种预测岩石杆中应力波大小的方法 Download PDFInfo
- Publication number
- CN113866023A CN113866023A CN202110992728.7A CN202110992728A CN113866023A CN 113866023 A CN113866023 A CN 113866023A CN 202110992728 A CN202110992728 A CN 202110992728A CN 113866023 A CN113866023 A CN 113866023A
- Authority
- CN
- China
- Prior art keywords
- rock
- strain
- rod
- rock rod
- speed
- 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 100
- 238000000034 method Methods 0.000 title claims abstract description 13
- 238000006073 displacement reaction Methods 0.000 claims abstract description 28
- 238000005259 measurement Methods 0.000 claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims abstract description 8
- 238000002474 experimental method Methods 0.000 claims description 3
- 238000012935 Averaging Methods 0.000 claims description 2
- 230000007423 decrease Effects 0.000 claims description 2
- 238000009795 derivation Methods 0.000 claims description 2
- 239000010438 granite Substances 0.000 claims description 2
- 238000012545 processing Methods 0.000 claims description 2
- 238000005507 spraying Methods 0.000 claims description 2
- 230000001502 supplementing effect Effects 0.000 claims description 2
- 230000001902 propagating effect Effects 0.000 claims 1
- 230000007547 defect Effects 0.000 abstract description 2
- 239000000853 adhesive Substances 0.000 description 1
- 230000001070 adhesive effect Effects 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000000691 measurement method Methods 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N3/30—Investigating strength properties of solid materials by application of mechanical stress by applying a single impulsive force, e.g. by falling weight
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P3/00—Measuring linear or angular speed; Measuring differences of linear or angular speeds
- G01P3/36—Devices characterised by the use of optical means, e.g. using infrared, visible, or ultraviolet light
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
- G01S17/02—Systems using the reflection of electromagnetic waves other than radio waves
- G01S17/06—Systems determining position data of a target
- G01S17/08—Systems determining position data of a target for measuring distance only
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/06—Indicating or recording means; Sensing means
- G01N2203/067—Parameter measured for estimating the property
- G01N2203/0676—Force, weight, load, energy, speed or acceleration
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
- G01N2203/02—Details not specific for a particular testing method
- G01N2203/06—Indicating or recording means; Sensing means
- G01N2203/067—Parameter measured for estimating the property
- G01N2203/0682—Spatial dimension, e.g. length, area, angle
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Electromagnetism (AREA)
- Immunology (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Chemical & Material Sciences (AREA)
- Pathology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biochemistry (AREA)
- Computer Networks & Wireless Communication (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Power Engineering (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种预测岩石杆中应力波大小的方法,通过非接触测量得到的端部速度,提供一种精确可靠的计算方法用来确定任意岩石截面的应力、应变、速度和位移参数。通过端部速度就可以确定任意截面位置x处的应力、应变、速度和位移,避免了对多个截面的重复测量。利用激光测距仪与岩石杆右自由端的非接触测量,不需要传统应变片的粘贴,克服了在高温下无法使用应变片的缺点,因此可以在高温状态下得到任意截面的应力、应变、速度和位移,扩展了测量的温度范围。可以通过预测任意截面的应力判断该截面处的应力是否对岩石杆造成破坏。
Description
技术领域
本发明涉及一种预测岩石杆中应力波大小的方法,特别是根据端部速度预测岩石杆中任意一截面应力、应变、速度和位移大小的方法
背景技术
在冲击动力学实验中,如何精确的获得介质内部的应力、应变、速度、位移等参数一直是重要的研究课题。为了研究冲击产生的应力波对岩石造成的应变,通过在岩石杆不同位置上粘贴应变片来获得岩石表面的应变,在得到应变的基础上根据连续方程和运动方程求解得到该截面的速度v和应力σ,在得到速度的基础上通过公式得到截面位移。这种获得截面相关参数的方法受限于应变片的粘贴位置,只能得到岩石杆部分截面的应变。在面对高温岩体时,传统的测量方式由于在高温下应变片和粘合剂的失效,无法得到高温岩体杆的截面应变。随着非接触测量方法的逐渐成熟,较高精度的测量岩体受冲击作用下产生的速度成为了可能。因此利用非接触测量获得的速度进而得到岩石杆任意截面应力、应变、速度和位移成为可能,可以为研究高温岩体在冲击荷载作用下的力学行为提供重要依据。
发明内容
本发明通过非接触测量得到的端部速度,提供一种精确可靠的计算方法用来确定任意岩石截面的应力、应变、速度和位移参数。
本发明的技术方案:
本发明包括对岩石杆密度ρ,端部速度v(L,t)、衰减系数α(ω)和波数k(ω)的测定,进而计算得到频域下端部速度右行波和左行波利用以上结果计算得到频域内任意截面x处的速度应变应力和位移进行傅里叶逆变换得到时域内任意截面x处的应变ε(x,t),速度v(x,t),应力σ(x,t)和位移u(x,t)。
详细描述如下:
得到岩石杆密度ρ:利用卷尺测量得到岩石杆长度L,所采用的岩石杆为圆柱体,可以利用游标卡尺测量得到横截面直径D,进而计算得到横截面积A,其中计算方法为式中π为圆周率,为定值。进而得到岩石杆体积V,其计算方法为V=AL。利用电子秤得到岩石杆质量m。利用公式计算得到岩体密度ρ。
如图2所示建立沿岩石杆轴线方向的x坐标轴,令岩石杆左端受撞击的端部坐标为0,右自由端坐标为L即岩石杆长度。将激光测速仪垂直对准岩石杆右端部,对岩石杆左端施加冲击荷载,得到右自由端在冲击荷载作用下的速度v(L,t),利用公式得到经傅里叶变换后频域下右自由端速度其中ω为傅里叶变换后频域内的谐波频率,t为岩石杆受到冲击后的时间。
利用右自由端速度得到岩石杆任意一截面的应变、速度、应力和位移的公式推导过程如下:
在动态冲击作用下的一维波传播方程可以表示为
和
这里x为截面位置,σ,ε和u分别为应力、应变和轴向位移。
岩石杆在频域下的本构关系可以表示为
波传播系数γ(ω)可以定义为
根据一维波传播方程可以表示为
式的通解表示为
这里为在岩石杆x位置处在频域内的应变,传播系数为γ(ω)=α(ω)+ik(ω),α(ω)为衰减系数,k(ω)为波数;右行波和左行波分别为当x=0时,分别沿x增大和减小方向传播应变的傅里叶变换,可以由测量的右自由端速度确定。根据连续方程可以得到频域下x处的应变为
L处为岩石右端自由面,因此根据应变为0可得
式中传播系数γ(ω)可以通过以下公式确定
可得
对式,,和进行傅里叶逆变换得到岩石杆任意一截面时域下的应变ε(x,t),速度v(x,t),应力σ(x,t)和位移u(x,t)关系。
与现有技术相比较,本发明仅通过端部速度就可以确定任意截面位置x处的应力、应变、速度和位移,避免了对多个截面的重复测量。
本发明利用激光测距仪与岩石杆右自由端的非接触测量,不需要传统应变片的粘贴,克服了在高温下无法使用应变片的缺点,因此可以在高温状态下得到任意截面的应力、应变、速度和位移,扩展了测量的温度范围。
本发明可以通过预测任意截面的应力判断该截面处的应力是否对岩石杆造成破坏。
附图说明
图1为验证方法正确性的实验装置
图2为岩石杆坐标示意图
图中:1空气压缩机;2发射腔室;3入射杆;4测速器;5DIC处理软件;6高速摄影机;7激光测速仪控制器;8激光测速器;9超动态应变仪;10应变片;11岩石杆;12补光灯
具体实施方式
如图1所示,采用现有的冲击装置来实施本发明,下面做进一步的详细说明:
第一步,选择适合的岩石杆试样,本发明基于一维波传播理论预测岩石截面力学参数。因此选择长度1.2m,直径45mm的圆柱形杆,长径比为26.7,本次冲击产生的压缩波长度为80cm,远小于岩石杆直径,满足一维波传播理论假定。
第二步,测量岩石杆密度ρ。首先采用卷尺三次测量岩石杆两端面长度,取平均值,得到岩石杆长度为120.13cm。然后采用游标卡尺测量岩石杆直径,得到直径为4.52cm。利用公式计算得到岩石杆截面面积为16.05cm2。采用电子秤测量得到岩石杆重量为7.453kg,计算得到花岗岩杆密度为3760kg/m3。
第三步,将岩石杆11水平放置在定滑轮上,调整岩石杆位置使中心轴线与入射杆轴线在同一水平线上。将选用的岩石11中间粘贴上应变片10,将应变片与超动态应变仪9相连接采集在冲击荷载的作用下岩石杆中间应变,其中超动态采集频率为100ksps,采用的触发模式为单次触发。将岩石杆中间除应变片位置喷涂散斑,将高速摄影机6拍摄区域对准应变片所测量截面。设定拍摄参数:帧率100kHz,像素256×128,镜头为焦距105mm放大倍率为1:1的微距镜头。采用一盏2000w的LED灯12对散斑区域进行补光。将拍摄结果输入含有DIC软件的电脑5进行实时处理得到岩石杆中间截面的速度和位移。将激光测速仪8垂直对准岩石杆11,采集帧率为100ksps,将测量所得结果输出至电脑7中。
第四步,打开空气压缩机1,使发射腔室2充满高压气体,打开阀门使入射杆3高速离开发射腔室,利用测速仪4测量得到入射杆速度,并同时激发超动态应变仪9,高速摄影机6和激光测速仪8。
Claims (2)
1.一种预测岩石杆中应力波大小的方法,其特征在于:包括对岩石杆密度ρ,端部速度v(L,t)、衰减系数α(ω)和波数k(ω)的测定,进而计算得到频域下端部速度右行波和左行波利用以上结果计算得到频域内任意截面x处的速度应变应力和位移进行傅里叶逆变换得到时域内任意截面x处的应变ε(x,t),速度v(x,t),应力σ(x,t)和位移u(x,t);
得到岩石杆密度ρ:利用卷尺测量得到岩石杆长度L,所采用的岩石杆为圆柱体,可以利用游标卡尺测量得到横截面直径D,进而计算得到横截面积A,其中计算方法为式中π为圆周率,为定值;进而得到岩石杆体积V,其计算方法为V=AL;利用电子秤得到岩石杆质量m;利用公式计算得到岩体密度ρ;
建立沿岩石杆轴线方向的x坐标轴,令岩石杆左端受撞击的端部坐标为0,右自由端坐标为L即岩石杆长度;将激光测速仪垂直对准岩石杆右端部,对岩石杆左端施加冲击荷载,得到右自由端在冲击荷载作用下的速度v(L,t),利用公式得到经傅里叶变换后频域下右自由端速度其中ω为傅里叶变换后频域内的谐波频率,t为岩石杆受到冲击后的时间;
利用右自由端速度得到岩石杆任意一截面的应变、速度、应力和位移的公式推导过程如下:
在动态冲击作用下的一维波传播方程表示为
和
这里x为截面位置,σ,ε和u分别为应力、应变和轴向位移;
岩石杆在频域下的本构关系表示为
波传播系数γ(ω)定义为
根据一维波传播方程表示为
式的通解表示为
这里为在岩石杆x位置处在频域内的应变,传播系数为γ(ω)=α(ω)+ik(ω),α(ω)为衰减系数,k(ω)为波数;右行波和左行波分别为当x=0时,分别沿x增大和减小方向传播应变的傅里叶变换,可以由测量的右自由端速度确定;根据连续方程可以得到频域下x处的应变为
L处为岩石右端自由面,根据应变为0得
式中传播系数γ(ω)通过以下公式确定
可得
对式,,和进行傅里叶逆变换得到岩石杆任意一截面时域下的应变ε(x,t),速度v(x,t),应力σ(x,t)和位移u(x,t)关系;
2.根据权利要求1所述的一种预测岩石杆中应力波大小的方法,其特征在于:第一步,选择适合的岩石杆试样,基于一维波传播理论预测岩石截面力学参数进行分析计算;
第二步,测量岩石杆密度;采用卷尺三次测量岩石杆两端面长度,取平均值,得到岩石杆长度;然后采用游标卡尺测量岩石杆直径,得到直径;利用公式计算得到岩石杆截面面积;采用电子秤测量得到岩石杆重量,计算得到花岗岩杆密度;
第三步,将岩石杆水平放置在定滑轮上,调整岩石杆位置使中心轴线与入射杆轴线在同一水平线上;将选用的岩石中间粘贴上应变片,将应变片与超动态应变仪相连接采集在冲击荷载的作用下岩石杆中间应变,采用的触发模式为单次触发;将岩石杆中间除应变片位置喷涂散斑,将高速摄影机拍摄区域对准应变片所测量截面;采用LED灯对散斑区域进行补光;将拍摄结果输入含有DIC软件的电脑进行实时处理得到岩石杆中间截面的速度和位移;将激光测速仪垂直对准岩石杆,将测量所得结果输出至电脑中;
第四步,打开空气压缩机,使发射腔室充满高压气体,打开阀门使入射杆高速离开发射腔室,利用测速仪测量得到入射杆速度,并同时激发超动态应变仪,高速摄影机和激光测速仪;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110992728.7A CN113866023B (zh) | 2021-08-27 | 2021-08-27 | 一种预测岩石杆中应力波大小的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110992728.7A CN113866023B (zh) | 2021-08-27 | 2021-08-27 | 一种预测岩石杆中应力波大小的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113866023A true CN113866023A (zh) | 2021-12-31 |
CN113866023B CN113866023B (zh) | 2023-11-10 |
Family
ID=78988525
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110992728.7A Active CN113866023B (zh) | 2021-08-27 | 2021-08-27 | 一种预测岩石杆中应力波大小的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113866023B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113865986A (zh) * | 2021-08-27 | 2021-12-31 | 北京工业大学 | 利用高速摄影机和dic技术非接触检测实时高温岩体传播系数装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103245732A (zh) * | 2013-04-10 | 2013-08-14 | 湘潭大学 | 一种基于小波提取的管道注浆质量应力波检测信号处理方法 |
CN107703161A (zh) * | 2017-06-06 | 2018-02-16 | 中冶建筑研究总院有限公司 | 一种冲击应力波检测系统 |
CN109990968A (zh) * | 2019-03-22 | 2019-07-09 | 西北核技术研究所 | 一种基于时间反转法原理的板结构冲击定位方法 |
CN110715865A (zh) * | 2019-10-30 | 2020-01-21 | 中国兵器工业规划研究院 | 脆性材料动态断裂的力学和电磁响应同步测试系统及方法 |
CN112557499A (zh) * | 2020-06-30 | 2021-03-26 | 东南大学 | 一种基于超声波的节理对应力波透反射规律影响的实验方法 |
CN112986012A (zh) * | 2021-02-09 | 2021-06-18 | 北京工业大学 | 一种研究高温下应力波在岩体中传播特性的实验装置 |
CN113091973A (zh) * | 2021-03-04 | 2021-07-09 | 西安交通大学 | 基于内部弹性波非线性特征的激光冲击强化实时监测方法 |
-
2021
- 2021-08-27 CN CN202110992728.7A patent/CN113866023B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103245732A (zh) * | 2013-04-10 | 2013-08-14 | 湘潭大学 | 一种基于小波提取的管道注浆质量应力波检测信号处理方法 |
CN107703161A (zh) * | 2017-06-06 | 2018-02-16 | 中冶建筑研究总院有限公司 | 一种冲击应力波检测系统 |
CN109990968A (zh) * | 2019-03-22 | 2019-07-09 | 西北核技术研究所 | 一种基于时间反转法原理的板结构冲击定位方法 |
CN110715865A (zh) * | 2019-10-30 | 2020-01-21 | 中国兵器工业规划研究院 | 脆性材料动态断裂的力学和电磁响应同步测试系统及方法 |
CN112557499A (zh) * | 2020-06-30 | 2021-03-26 | 东南大学 | 一种基于超声波的节理对应力波透反射规律影响的实验方法 |
CN112986012A (zh) * | 2021-02-09 | 2021-06-18 | 北京工业大学 | 一种研究高温下应力波在岩体中传播特性的实验装置 |
CN113091973A (zh) * | 2021-03-04 | 2021-07-09 | 西安交通大学 | 基于内部弹性波非线性特征的激光冲击强化实时监测方法 |
Non-Patent Citations (3)
Title |
---|
张振;王永刚;: "基于激光干涉测试技术的分离式Hopkinson压杆实验测试系统", 爆炸与冲击, no. 05 * |
杨仁树;李炜煜;李永亮;方士正;朱晔;: "3种岩石动态拉伸力学性能试验与对比分析", 煤炭学报, no. 09 * |
金解放;王熙博;袁伟;: "轴向静载荷作用下变截面红砂岩应力波衰减特性试验研究", 矿业研究与开发, no. 05 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113865986A (zh) * | 2021-08-27 | 2021-12-31 | 北京工业大学 | 利用高速摄影机和dic技术非接触检测实时高温岩体传播系数装置 |
CN113865986B (zh) * | 2021-08-27 | 2023-12-29 | 北京工业大学 | 利用高速摄影机和dic技术非接触检测实时高温岩体传播系数装置 |
Also Published As
Publication number | Publication date |
---|---|
CN113866023B (zh) | 2023-11-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108375501B (zh) | 一种基于分离式霍普金森压杆实验技术的数据处理方法 | |
Hayashi et al. | Wave structure analysis of guided waves in a bar with an arbitrary cross-section | |
Zhang et al. | A review of dynamic experimental techniques and mechanical behaviour of rock materials | |
Marais et al. | Material testing at high strain rate using the split Hopkinson pressure bar | |
Marineau et al. | Investigation of Mach 10 boundary layer stability of sharp cones at angle-of-attack, part 1: experiments | |
CN106949861B (zh) | 一种非线性超声在线监测金属材料应变变化的方法 | |
Louar et al. | Explosive driven shock tube loading of aluminium plates: Experimental study | |
Wilson et al. | Fastener failure detection using a surface acoustic wave strain sensor | |
CN102564856A (zh) | 基于数字图像相关的塑性多缺陷材料m积分测量方法 | |
Poozesh et al. | Full field inspection of a utility scale wind turbine blade using digital image correlation | |
Panowicz et al. | Tensile split Hopkinson bar technique: numerical analysis of the problem of wave disturbance and specimen geometry selection | |
Marineau et al. | Sharp cone boundary-layer transition and stability at Mach 14 | |
CN113866023A (zh) | 一种预测岩石杆中应力波大小的方法 | |
Chen et al. | Strain monitoring on damaged composite laminates using digital image correlation | |
Zhu et al. | The contact characteristics analysis for rod fastening rotors using ultrasonic guided waves | |
Jin et al. | The yielding behavior and plastic deformation of oxygen-free copper under biaxial quasi-static and dynamic loadings | |
CN113688544A (zh) | 一种主被动结合的复合材料损伤定量辨识方法 | |
Moćko et al. | Application of FEM in the assessment of phenomena associated with dynamic investigations on a miniaturised DICT testing stand | |
Lodygowski et al. | Dry sliding friction experiments at elevated velocities | |
Zhenkun et al. | Noncontact optical measurement of CTOA and CTOD for interface crack in DCB test | |
CN104913988A (zh) | 基于霍普金森原理的混凝土轴心抗拉强度测量方法 | |
Capponi et al. | Thermoelasticity and ArUco marker-based model validation of polymer structure: application to the San Giorgio’s bridge inspection robot | |
Zanarini | Full Field ESPI Vibration Measurements to Predict Fatigue Behavior | |
House et al. | Film data reduction from Taylor impact tests | |
Marinho et al. | Comparison between direct and indirect measurement methods for bulge tests |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |