CN109283079A - 一种测定及计算岩石衰减系数及非线性系数的方法 - Google Patents

一种测定及计算岩石衰减系数及非线性系数的方法 Download PDF

Info

Publication number
CN109283079A
CN109283079A CN201811113345.2A CN201811113345A CN109283079A CN 109283079 A CN109283079 A CN 109283079A CN 201811113345 A CN201811113345 A CN 201811113345A CN 109283079 A CN109283079 A CN 109283079A
Authority
CN
China
Prior art keywords
wave
sample
strain
plastic strain
bar
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
CN201811113345.2A
Other languages
English (en)
Other versions
CN109283079B (zh
Inventor
童立红
丁海滨
于洋
徐长节
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
East China Jiaotong University
Original Assignee
East China Jiaotong University
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 East China Jiaotong University filed Critical East China Jiaotong University
Priority to CN201811113345.2A priority Critical patent/CN109283079B/zh
Publication of CN109283079A publication Critical patent/CN109283079A/zh
Application granted granted Critical
Publication of CN109283079B publication Critical patent/CN109283079B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/32—Investigating strength properties of solid materials by application of mechanical stress by applying repeated or pulsating forces
    • G01N3/36—Investigating strength properties of solid materials by application of mechanical stress by applying repeated or pulsating forces generated by pneumatic or hydraulic means
    • 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/0001—Type of application of the stress
    • G01N2203/001—Impulsive
    • 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/003—Generation of the force
    • G01N2203/0042—Pneumatic or hydraulic means
    • G01N2203/0044—Pneumatic means
    • 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/0202—Control of the test
    • G01N2203/0212—Theories, calculations
    • G01N2203/0218—Calculations based on experimental data

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

一种测定及计算岩石衰减系数及非线性系数的方法，包括：将制备好的圆柱形岩石试样放置于霍普金森杆入射杆与透射杆之间，利用应变片记录入射波信号（6）和透射波应变信号（7）；对测得的透射波信号进行频谱分析，得出透射波频谱；由动应力‑应变关系及应变波传播的波动方程，近似得出试样一阶应变及二次谐波应变表达式；计算平均入射波应变幅值和透射波应变幅值比值，结合一阶应变，求解试样衰减系数；计算平均透射波谐波应变幅值与入射应变幅值平方的比值，结合二次谐波，求解试样非线性系数。本发明通过该方法实现颗粒材料的衰减系数及非线性系数的联合测定及计算。

Description

一种测定及计算岩石衰减系数及非线性系数的方法
技术领域
本发明涉及一种本发明涉及一种测定及计算岩石衰减系数及非线性系数的方法，属于土木工程技术领域。
背景技术
岩石和土体等土木工程中常见的地质材料具有特殊的颗粒结构，从而导致其具有非常明显的动力非线性特性。同时，颗粒间的孔隙及材料中存在缺陷(如裂缝、空洞等)对波携带的能量具有吸收效应，从而使波幅随着传播距离的增加而减小，在工程中通常用波的衰减来表述。
目前，对于颗粒材料对波的衰减及应力波作用下颗粒材料表现出的非线性特性已有定性的认识，但如何通过试验结合计算确定颗粒材料的衰减系数及非线性系数仍是一个亟待解决的技术问题。
发明内容
本发明的目的是，为了解决颗粒材料的衰减系数及非线性系数的如何确定问题，本发明提出一种测定及计算岩石衰减系数及非线性系数的方法。
本发明实现的技术方案如下，一种测定及计算岩石衰减系数及非线性系数的方法，包括以下步骤：
步骤1：制备圆柱形岩石试样，将岩石试样放置于霍普金森杆入射杆与透射杆之间，由发射腔产生的高应变波气压推动弹头作用于入射杆，产生的高应变波经入射杆传播至岩石试件，最终传入透射杆。
所述圆柱形岩石试样其直径应与霍普金森杆直径相同，所述圆柱形岩石试样其高度应为直径的1～2倍。
所述发射腔可调节气压大小，以达到产生不同幅值应变波的目的。
所述入射杆和透射杆在距离试件两端面1m处布置应变片，以记录入射波和透射波应变信号。
步骤2：采用快速傅里叶变换对步骤1测得的透射波信号进行频谱分析，得出透射波频谱。
步骤3：根据非线性弹性理论，动应力-应变关系可表示为：
式中：σd为岩石试样动应力；εd为岩石试样动应变；为线性动模量，ρs为岩石试样密度；β为非线性系数；η为岩石试样粘性系数；t为时间。
高应变波在岩石试样中的传播可视为一维波动问题，其波动方程为：
将方程(2)代入方程(1)，可得出非线性平面波方程为：
计算时，假设岩石试样长度远小于高应变波形成距离，同时采用摄动法，可近似得出试样一阶应变及二次谐波应变：
εd1(a，t)＝Re(ε0e-αaejka-jωt)
式中：k＝ω/Cs为冲击波波数，ω＝2πf为角频率，f为基频，可由透射波频谱分析直接得出，Cs为岩石试样中波速可由压电声速装置直接测出或由采集的入射信号和透射信号的时间差算得；α为岩石试样衰减系数；β为试样非线性系数；a为冲击在试样中的传播距离0≤a≤L，L为试样长度，a＝0表示试样最左端，a＝L表示试件最右端；
ε0＝εinT1为试件最左端应变值，εin为入射波应变幅值可由入射杆应变片直接测出，为入射杆-试样间界面透射系数，为试样对波的阻抗比，ρb为入射杆或透射杆的密度；Cb为入射杆或透射杆的波速。
步骤4：引入基频下透射应变εtr 1，入射波应变幅值和透射波应变幅值比值ra，计算岩石衰减系数。
基频下，试样最右端动应变为εd1＝Re(ε0e-αLejkL-jωt)＝Re(εinT1e-αLejkL-jωt)，由此得出透射应变为：
εtr1＝εd1T2＝Re(εinT1T2e-αLejkL-jωt) (5)
式中：εtr1为一阶透射应变幅值，可由透射波频谱直接得到；T2为试样-透射杆界面透射系数，其表达式为，
将T1和T2带入式(5)得出：
通过式(6)可得出入射波应变幅值和透射波应变幅值的比值为：
试验时记录同一试样在不同高应变荷载下入射波应变幅值和透射波应变幅值，计算透射波应变幅值与入射波应变幅值比ra，剔除误差过大的试验结果，求得平均幅值比ra。最后由测得的ra，通过式(7)可计算出衰减系数α。
步骤5：引入透射波的谐波应变εtr2，透射波谐波应变幅值与入射应变幅值平方的比值计算岩石非线性系数。
试样最右端谐波应变由此得出透射杆中谐波应变为：
式(8)中εtr2为透射波信号中的谐波信号，可由透射波频谱直接得到。
通过式(8)可得出比值：
式(9)中衰减系数α可由式(6)得出，其余与参数均为已知值，计算同一试样不同高应变荷载下透射波谐波应变幅值与入射应变幅值平方的比值求出平均比值，最后可通过式(9)求出岩石试样非线性系数β。
本发明的有益效果是，本发明利用应变片记录岩石试样入射波和透射波应变信号，并对测得的透射波信号进行频谱分析，得出透射波频谱；根据非线性弹性理论，由动应力-应变关系及应变波传播的波动方程，采用摄动法，近似得出试样一阶应变及二次谐波应变表达式；引入基频下透射应变，入射波应变幅值和透射波应变幅值比值，计算岩石衰减系数；引入透射波的谐波应变，透射波谐波应变幅值与入射应变幅值平方的比值，计算出岩石非线性系数，从而实现颗粒材料的衰减系数及非线性系数的联合测定及计算方法。
附图说明
图1为测定及计算岩石衰减系数及非线性系数的装置和波形示意图；
图2为本发明测定及计算高应变荷载作用下岩石衰减系数及非线性系数的流程图；
图3为一阶透射波幅值与入射波幅值比；
图4为透射波谐波幅值与入射波幅值平方比；
图中：1为入射波；2为透射波；3、4为应变片；5为岩石试样；6为入射波信号；7为透射波信号；8为围压系统。
具体实施方式
下面结合附图，针对具体实例，对本发明的具体实施方法予以说明。
本实例采用直径为D＝48.4mm，高度为L＝70.5mm的圆柱形砂岩(颗粒粒径0.2～0.5mm)作为试验试样，试验前测得砂岩密度ρs＝2.3g/cm3，波速Cs＝2900m/s，杆件密度ρb＝7.81g/cm3，波速Cb＝5410m/s。本次试验是在无围压条件下进行。
(1)霍普金森杆冲击气压依次设置为0.6MPa、0.7MPa、0.8MPa和0.9MPa。测得不同冲击荷载下入射波应变和透射波应变信号。
(2)对(1)中提取的透射波应变信号进行快速傅里叶变化，得出入射波的基频为f1＝4524Hz，倍频为f2＝9048Hz。
(3)统计(1)中四组冲击荷载下，透射波应变幅值与入射波应变幅值的比值(如图3所示)，求出平均比值为ra≈0.342。
求得阻抗比
求得波数k＝ω/Cs＝2πf1/Cs＝9.8m-1；
利用方程式求得衰减系数α≈4.54m-1。
(4)经快速傅里叶变换得出各组高应变荷载下所得透射波谐波应变幅值εtr2，并计算比值(如图4所示)，得出平均比值为衰减系数采用(3)中计算值α≈4.54m-1，其余参数计算与(3)中类似。
最后利用方程求得高应变荷载作用下岩石衰减系数β≈629。

Claims (2)

1.一种测定及计算岩石衰减系数及非线性系数的方法，其特征在于，所述方法包括以下步骤：
(1)制备圆柱形岩石试样，将岩石试样放置于霍普金森杆的入射杆与透射杆之间，由高压发射腔产生的高应变波气压推动弹头作用于入射杆，高应变波经入射杆传播至岩石试件，最终传入透射杆，利用应变片记录入射波和透射波应变信号；
(2)采用快速傅里叶变换对步骤(1)测得的透射波信号进行频谱分析，得出透射波频谱；
(3)根据非线性弹性理论，由动应力-应变关系及应变波传播的波动方程，采用摄动法，近似得出试样一阶应变及二次谐波应变表达式：
εd1(a,t)＝Re(ε0e-αaejka-jωt)
式中：k＝ω/Cs为冲击波波数；ω＝2πf为角频率；f为基频，可由透射波频谱分析直接得出；Cs为岩石试样中波速，由压电声速装置直接测出或由采集的入射信号和透射信号的时间差算得；α为岩石试样衰减系数；β为试样非线性系数；a为冲击在试样中的传播距离，0≤a≤L，L为试样长度，a＝0表示试样最左端，a＝L表示试件最右端；ε0为试件最左端应变值，ε0＝εinT1，εin为入射波应变幅值，可由入射杆应变片直接测出；为入射杆-试样间界面透射系数；为试样对波的阻抗比；ρb为入射杆或透射杆的密度；Cb为入射杆或透射杆的波速；
(4)引入基频下透射应变εtr1，入射波应变幅值和透射波应变幅值比值ra，计算岩石衰减系数；
通过入射杆-试样和试样-透射杆之间高应变波的传播关系，得出入射波应变幅值和透射波应变幅值的比值；
试验时记录同一试样在不同冲击荷载下入射波应变幅值和透射波应变幅值，计算入射波应变幅值与透射波应变幅值比ra，剔除误差过大的试验结果，求得平均幅值比ra；最后由测得的ra，通过上式计算出衰减系数α；
(5)引入透射波的谐波应变εtr2，透射波谐波应变幅值与入射应变幅值平方的比值计算岩石非线性系数；
通过入射杆-试样和试样-透射杆之间高应变波的传播关系，结合二次谐波得出透射波谐波应变幅值与入射应变幅值平方的比值；
式中，εtr2为透射波信号中的谐波信号，可由透射波频谱直接得到；计算同一试样不同高应变荷载下透射波谐波应变幅值与入射应变幅值平方的比值求出平均比值，最后通过上式求得非线性系数β。
2.根据权利要求1所述的一种测定及计算岩石衰减系数及非线性系数的方法，其特征在于，所述圆柱形岩石试样的直径与霍普金森杆直径相同；所述圆柱形岩石试样的高度为直径的1～2倍。
CN201811113345.2A 2018-09-25 2018-09-25 一种测定及计算岩石衰减系数及非线性系数的方法 Active CN109283079B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811113345.2A CN109283079B (zh) 2018-09-25 2018-09-25 一种测定及计算岩石衰减系数及非线性系数的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811113345.2A CN109283079B (zh) 2018-09-25 2018-09-25 一种测定及计算岩石衰减系数及非线性系数的方法

Publications (2)

Publication Number Publication Date
CN109283079A true CN109283079A (zh) 2019-01-29
CN109283079B CN109283079B (zh) 2020-10-09

Family

ID=65181237

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811113345.2A Active CN109283079B (zh) 2018-09-25 2018-09-25 一种测定及计算岩石衰减系数及非线性系数的方法

Country Status (1)

Country Link
CN (1) CN109283079B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109883599A (zh) * 2019-03-12 2019-06-14 哈尔滨工程大学 一种基于应力波理论的水面爆炸对海底隧道抗暴能影响的试验方法
CN112881241A (zh) * 2021-01-19 2021-06-01 华东交通大学 一种确定颗粒材料模量软化和恢复的方法
CN112903477A (zh) * 2021-01-25 2021-06-04 华东交通大学 一种测定及计算颗粒材料系统剪切强度的方法

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000049458A1 (en) * 1999-02-19 2000-08-24 University Of Southampton Non-linearity compensation in optical devices and transmission systems
JP2007120946A (ja) * 2005-10-24 2007-05-17 National Institute Of Advanced Industrial & Technology 固体材料表層部の物性値測定方法
CN101140159A (zh) * 2007-10-12 2008-03-12 东南大学 微位移光学测试方法及其装置
CN101666780A (zh) * 2009-10-13 2010-03-10 北京工业大学 基于微观形貌的非完好界面接触状态非线性超声评价方法
CN104757999A (zh) * 2015-04-07 2015-07-08 声泰特（成都）科技有限公司 一种基于超声基波和谐波的非线性成像方法和系统
CN105548364A (zh) * 2015-12-10 2016-05-04 北京理工大学 热障涂层结合强度的高阶非线性参数表征方法
CN106688151A (zh) * 2014-09-25 2017-05-17 科磊股份有限公司 使用单片带宽窄化设备的激光组合件及检验系统
CN108175442A (zh) * 2018-01-15 2018-06-19 中国计量大学 一种超声造影剂非线性参数测量方法
CN108287118A (zh) * 2017-12-26 2018-07-17 中国石油天然气股份有限公司 一种岩石低频弹性模量及衰减系数测量方法及系统
CN108490081A (zh) * 2018-04-03 2018-09-04 南昌航空大学 一种基于多参数的合金晶粒尺寸的确定方法及确定系统
CN108509376A (zh) * 2018-03-26 2018-09-07 华东交通大学 一种基于位移模式修正的土压力计算方法

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2000049458A1 (en) * 1999-02-19 2000-08-24 University Of Southampton Non-linearity compensation in optical devices and transmission systems
JP2007120946A (ja) * 2005-10-24 2007-05-17 National Institute Of Advanced Industrial & Technology 固体材料表層部の物性値測定方法
CN101140159A (zh) * 2007-10-12 2008-03-12 东南大学 微位移光学测试方法及其装置
CN101666780A (zh) * 2009-10-13 2010-03-10 北京工业大学 基于微观形貌的非完好界面接触状态非线性超声评价方法
CN106688151A (zh) * 2014-09-25 2017-05-17 科磊股份有限公司 使用单片带宽窄化设备的激光组合件及检验系统
CN104757999A (zh) * 2015-04-07 2015-07-08 声泰特（成都）科技有限公司 一种基于超声基波和谐波的非线性成像方法和系统
CN105548364A (zh) * 2015-12-10 2016-05-04 北京理工大学 热障涂层结合强度的高阶非线性参数表征方法
CN108287118A (zh) * 2017-12-26 2018-07-17 中国石油天然气股份有限公司 一种岩石低频弹性模量及衰减系数测量方法及系统
CN108175442A (zh) * 2018-01-15 2018-06-19 中国计量大学 一种超声造影剂非线性参数测量方法
CN108509376A (zh) * 2018-03-26 2018-09-07 华东交通大学 一种基于位移模式修正的土压力计算方法
CN108490081A (zh) * 2018-04-03 2018-09-04 南昌航空大学 一种基于多参数的合金晶粒尺寸的确定方法及确定系统

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
HYUNJO JEONG ET AL.: "Receiver calibration and the nonlinearity parameter measurement of thick solid samples with diffraction and attenuation corrections", 《ULTRASONICS》 *
刘瑶璐等: "板壳结构中的非线性兰姆波", 《力学进展》 *
刘苗苗等: "基于ＳＨＰＢ波传播系数法试验的衰减弥散修正研究", 《实验力学》 *
张荣繁等: "材料非线性衰减系数的二次谐波测量方法研究", 《声学学报》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109883599A (zh) * 2019-03-12 2019-06-14 哈尔滨工程大学 一种基于应力波理论的水面爆炸对海底隧道抗暴能影响的试验方法
CN112881241A (zh) * 2021-01-19 2021-06-01 华东交通大学 一种确定颗粒材料模量软化和恢复的方法
CN112881241B (zh) * 2021-01-19 2022-10-28 华东交通大学 一种确定颗粒材料模量软化和恢复的方法
CN112903477A (zh) * 2021-01-25 2021-06-04 华东交通大学 一种测定及计算颗粒材料系统剪切强度的方法
CN112903477B (zh) * 2021-01-25 2022-08-26 华东交通大学 一种测定及计算颗粒材料系统剪切强度的方法

Also Published As

Publication number Publication date
CN109283079B (zh) 2020-10-09

Similar Documents

Publication Publication Date Title
CN104142195B (zh) 基于超声波法的钢结构构件内部初始应力检测装置和方法
CN108169330B (zh) 基于非线性超声谐波法的混凝土构件轴向应力无损检测的装置和方法
CN110618198B (zh) 一种保真环境下非接触式测量岩石波速的测试方法
CN107314933B (zh) 动及动静组合载荷下煤岩组合体中煤的力学特性试验方法
Vary Acousto-ultrasonic characterization of fiber reinforced composites
CN109283079A (zh) 一种测定及计算岩石衰减系数及非线性系数的方法
CN101140265B (zh) 采用非线性共振频率的偏移测定混凝土含水量和裂纹的方法
CN105738479B (zh) 弯曲元测试岩土工程材料小应变粘弹性参数的方法及装置
CN106949861B (zh) 一种非线性超声在线监测金属材料应变变化的方法
CN201653999U (zh) 一种智能超声横波纵波综合测试仪
Payan et al. Applying diffuse ultrasound under dynamic loading to improve closed crack characterization in concrete
CN102183585B (zh) 一种岩心选样的方法
CN107271280A (zh) 一种冲击回波‑回弹综合评定混凝土抗压强度方法
CN105842067A (zh) 应力变化与裂缝扩展方向测试装置及方法
CN111077030A (zh) 一种高应变率下混凝土动态力学性能试验装置及方法
Ayala et al. Study of the elastic shear modulus of Bio Bio sand using bender elements in an oedometer
CN110018062A (zh) 一种直剪试验中岩石结构面剪切破坏位置定位方法
CN106680375B (zh) 用于确定材料的弹性模量的空气耦合超声检测方法
CN108593377A (zh) 一种低频测量系统中的样品及制备方法和使用方法
KR100844173B1 (ko) 과실의 물성 측정을 위한 초음파 탐촉자
Peroni et al. High strain-rate compression test on metallic foam using a multiple pulse SHPB apparatus
CN108195725A (zh) 基于shpb的颗粒物质一维应力波衰减效应测定方法
Suwal et al. Upgrading disk transducer to measure elastic waves on coarse-grained granular materials: development and performance revelation
Piriyakul A development of a bender element apparatus
Pengfei et al. A volumetric strain-based method to determine crack initiation stress of low-porosity rocks

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