本发明涉及一种本发明涉及一种测定及计算岩石衰减系数及非线性系数的方法,属于土木工程技术领域。
背景技术:
岩石和土体等土木工程中常见的地质材料具有特殊的颗粒结构,从而导致其具有非常明显的动力非线性特性。同时,颗粒间的孔隙及材料中存在缺陷(如裂缝、空洞等)对波携带的能量具有吸收效应,从而使波幅随着传播距离的增加而减小,在工程中通常用波的衰减来表述。
目前,对于颗粒材料对波的衰减及应力波作用下颗粒材料表现出的非线性特性已有定性的认识,但如何通过试验结合计算确定颗粒材料的衰减系数及非线性系数仍是一个亟待解决的技术问题。
技术实现要素:
本发明的目的是,为了解决颗粒材料的衰减系数及非线性系数的如何确定问题,本发明提出一种测定及计算岩石衰减系数及非线性系数的方法。
本发明实现的技术方案如下,一种测定及计算岩石衰减系数及非线性系数的方法,包括以下步骤:
步骤1:制备圆柱形岩石试样,将岩石试样放置于霍普金森杆入射杆与透射杆之间,由发射腔产生的高应变波气压推动弹头作用于入射杆,产生的高应变波经入射杆传播至岩石试件,最终传入透射杆。
所述圆柱形岩石试样其直径应与霍普金森杆直径相同,所述圆柱形岩石试样其高度应为直径的1~2倍。
所述发射腔可调节气压大小,以达到产生不同幅值应变波的目的。
所述入射杆和透射杆在距离试件两端面1m处布置应变片,以记录入射波和透射波应变信号。
步骤2:采用快速傅里叶变换对步骤1测得的透射波信号进行频谱分析,得出透射波频谱。
步骤3:根据非线性弹性理论,动应力-应变关系可表示为:
式中:σd为岩石试样动应力;εd为岩石试样动应变;
高应变波在岩石试样中的传播可视为一维波动问题,其波动方程为:
将方程(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为入射波应变幅值可由入射杆应变片直接测出,
步骤4:引入基频下透射应变εtr1,入射波应变幅值和透射波应变幅值比值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)得出,其余与参数均为已知值,计算同一试样不同高应变荷载下透射波谐波应变幅值与入射应变幅值平方的比值
本发明的有益效果是,本发明利用应变片记录岩石试样入射波和透射波应变信号,并对测得的透射波信号进行频谱分析,得出透射波频谱;根据非线性弹性理论,由动应力-应变关系及应变波传播的波动方程,采用摄动法,近似得出试样一阶应变及二次谐波应变表达式;引入基频下透射应变,入射波应变幅值和透射波应变幅值比值,计算岩石衰减系数;引入透射波的谐波应变,透射波谐波应变幅值与入射应变幅值平方的比值,计算出岩石非线性系数,从而实现颗粒材料的衰减系数及非线性系数的联合测定及计算方法。
附图说明
图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)经快速傅里叶变换得出各组高应变荷载下所得透射波谐波应变幅值εtr2,并计算比值
最后利用方程