Dynamic crack propagation and arrest investigated with a cracked eccentrically-holed flattened disc of rock
-
摘要: 针对平台圆环构型的优点, 提出偏心圆孔单裂纹平台圆盘(cracked eccentrically holed flattened disc, CEHFD), 该试样具有更长的断裂路径。利用霍普金森压杆加载系统, 径向冲击CEHFD试样, 完成Ⅰ型动态断裂实验。砂岩试样表面粘贴应变片和裂纹扩展计, 用于监测裂纹动态起裂、扩展和止裂的全过程。实验表明, 在整个断裂过程中, 裂纹非匀速扩展, 裂纹扩展速度在裂纹起裂后加速上升, 在裂纹止裂前有明显的减速, 与地震时断层的动态破裂全过程完全吻合。采用实验-数值-解析法得到动态应力强度因子, 其时间历程呈现先增大后减小的趋势。根据断裂过程不同时刻, 得到相应的动态起裂韧度、扩展韧度及止裂韧度。在动态断裂全过程中, 动态扩展韧度为速度的函数, 变化趋势与速度一致, 随着时间先增大后减小; 动态起裂韧度大于动态止裂韧度, 止裂韧度随着裂纹最大扩展速度的增大而降低, 并且有较大的离散性。Abstract: Expanding the advantages of flattened ring geometry, a new test specimen geometry with a longer fracture path was advised to investigate the complete dynamic fracture process.This geometry was cracked eccentrically holed flattened disc (CEHFD) which was more conducive to study the unstable dynamic crack propagation and dynamic crack arrest.In order to study the mode Ⅰ (opening mode) dynamic fracture process of a rock, CEHFD specimens with strain gauges and crack propagation gauge glued on the specimen's surface were diametrically impacted by a split Hopkinson pressure bar.The experiment investigated crack initiation, rapid crack propagation, and crack arrest, all in one specimen.In a single complete fracture event, the crack accelerated after initiating and decelerated before arresting.The process monitored by the laboratory is completely consistent with the whole process of the dynamic rupture of faults in earthquake.A hybrid experimental-numerical-analytical method was used to determine the dynamic stress intensity factor that increased with increasing time.The dynamic fracture toughness was obtained with the fracture time information.Crack propagation velocity was first increased and then decreased with the time, and the dynamic propagation toughness had the same variation tendency with time as it was the function of velocity.The dynamic arrest toughness was smaller than the dynamic initiation toughness, and decreased when the maximum crack propagation speed increased.
-
Brian等[1]认为预测材料的失效是工程中最古老的问题之一, 也是解决得最不完美的问题之一。余寿文[2]把在外载荷及环境作用下, 材料和构件中缺陷的受力与变形视为使裂纹起裂与扩展的推动力, 比喻为“矛”; 把材料抵抗裂纹扩展的能力比喻为“盾”, 认为这一用断裂准则联系着的矛盾两体, 决定了裂纹的力学行为, 从而也就确定了构件和结构的安全及寿命。这一比喻不管对静态断裂, 还是动态断裂而言都是非常贴切的。岩石动态断裂研究涉及土木、矿业、能源、军事和防护等许多工程学科, 以及地质学、地球物理学等科学领域。例如, 地球物理学研究地震断层的高速扩展, 需要了解岩石的动态断裂力学参数及其在地震断裂带上岩性的分布。
动态断裂全过程包括裂纹的动态起裂、快速扩展和止裂过程。Das[3]认为由地震产生的地面加速度, 导致了地震破裂速度的剧烈变化, 从而影响地震灾害的程度, 这是地震学家研究地震断裂动力学的主要动机。在另一篇文章中, Das[4]则明确表示更快的破裂速度将会导致更为严重的地震破坏, 地震带来的灾害程度很大程度上与岩石破裂速度到底有多快相关, 因此在进行抗震结构设计时, 必须加以考虑。Robinson等[5]研究2011年的昆仑山地震, 发现地震破裂起始时速度较为缓慢, 加速到超剪切波速后, 以接近6 km/s的速度扩展近100 km, 在裂纹止裂前, 速度减速扩展。在涉及到裂纹传播速度与裂纹止裂时, 存在着两方面的分歧。一方面认为:裂纹止裂之前观察到了裂纹传播速度的减速[6-7]; 而另一方面则观察到突然快速止裂[8-10]。Ravi-Chandar等[9]详细说明了造成这种分歧的可能原因, 即实验室中实验的试样相对来说是比较小的, 在裂纹扩展过程中应力波的边界反射、加载波与卸载波的相互交叉都会给裂纹尖端造成一个非常复杂的应力历程, 并且会随试样构型而改变, 这就可能造成观测到不同的速度历程。然而, 随着实验技术的发展, 更多的学者倾向于认为裂纹在止裂前存在减速的过程, 但是是否由反射应力波的影响造成, 其影响多大却有待更深入的研究。
在动态断裂研究中, 动态应力强度因子(dynamic stress intensity factor)表征了动态断裂过程中裂纹尖端的应力应变场, 为了研究脆性材料的动态断裂行为, 许多学者将研究的重点放在了动态应力强度因子上。理论分析表明:裂纹尖端的应力分布与裂纹速度有关, 动态应力强度因子也就与速度有关。当裂纹在以某一速度匀速扩展时, 动态应力强度因子可以表示成裂纹速度的函数, Freund[10]给出了这一复杂的表达式。实验研究发现, 动态应力强度因子和裂纹扩展速度之间存在一定的关系, 不同的实验得到不同的结论[11-13]。Fineberg等[14]研究带裂纹板静态加载下的裂纹扩展速度, 认为由于边界应力波反射的影响, 线弹性动态断裂力学有局限性, 在公式中引入加速度的概念。Arakawa等[15]研究静态加载下动态应力强度因子与裂纹速度的关系, 发现裂纹在扩展过程中, 速度有加速上升和减速下降的变化, 并且同一速度对应2个不同的动态应力强度因子。Singh等[16]研究动态加载下双材料界面间的断裂过程, 得到裂纹快速扩展的速度变化和动态应力强度因子的时间历程, 认为裂纹速度的减速现象是由反射应力波造成的。通常, 裂纹速度随着时间变化, 其影响改变了裂纹尖端的动态应力场, 研究裂纹的速度变化是理解动态裂纹扩展的关键。
岩石材料一般为准脆性材料, 通常其抗拉强度比抗压强度要低一个数量级[17]。与金属材料不同, 对岩石进行直接拉伸实验较难, 因此多采用间接拉伸的方式进行实验, 圆盘类试样也多用于测试岩石的动态断裂韧度[18-21]。周妍等[22-23]采用圆孔内单边裂纹平台巴西圆盘(holed single cracked flattened Brazilian disc, HSCFBD)新试样, 分别对其进行了静态和动态断裂韧度的测试, 得到砂岩的静态、动态断裂韧度值。李炼等[24]采用HSCFBD试样首次研究了砂岩的动态断裂全过程, 实验成功检测到试样的破裂历经裂纹的动态起裂、扩展和止裂。
本文中, 在采用HSCFBD试样研究岩石动态断裂的全过程的基础上, 进一步研究裂纹的快速扩展与止裂过程。结合平台圆环类试样的优点, 提出偏心圆孔单裂纹平台圆盘(cracked eccentrically holed flattened disc, CEHFD)这一新型试样, 内圆孔偏离圆心一定距离, 使其断裂带加长, 增加其断裂扩展路径和检测范围, 且可使止裂点远离边界减小其反射应力波相互作用带来的影响。
1. 动态断裂实验
1.1 偏心圆孔单裂纹平台圆盘
CEHFD试样构型如图 1所示。图中R为外圆半径, r为内圆孔半径, 2β为平台加载角度。圆盘圆心为O, 内圆孔圆心为O′, O与O′之间的距离为内圆孔偏离外圆心的距离, 本文定义为偏心距离d; a0为预制裂纹长度; 断裂带长度L=R+d-r。p(t)为动态载荷, 是时间的函数。
实验中所有试样均由青砂岩材料制作, 颗粒细致且紧密, 如图 2所示。经测定该产地青砂岩的泊松比μ=0.21, 弹性模量E=17.67 GPa, 密度ρ=3.055 g/cm3。膨胀波波速cd=2 551.6 m/s, 畸变波波速cs=1 546.0 m/s, Rayleigh波波速cR=1 411.0 m/s。
实验选用试样的尺寸具体为:厚度B=(30±0.2) mm, R=(75±0.2) mm, r=(15±0.2) mm, a0=(8±0.1) mm, 平台角2β=30°, d=(20±0.2) mm, L=R+d-r =80 mm。采用裂缝加工系统预制裂纹, 选用直径为0.15 mm超细精钢砂切割线。使用塞尺测量预制裂纹宽度, 裂缝宽度为0.2 mm, 满足国际岩石力学学会规定的采用岩石试样测量断裂韧度对裂纹尖端宽度尺寸的要求。
1.2 实验加载装置及测量仪器
动态实验采用直径为100 mm的分离式霍普金森压杆(split Hopkinson pressure bar, SHPB)作为冲击加载装置, 其装置如图 3所示。SHPB压杆材料为42CrMo, 弹性模量Eb=210 GPa, 泊松比μb =0.3, 密度ρb=7 850 kg/m3, 一维纵波的理论波速c0 =5 172 m/s。入射杆的长度li=4 500 mm, 透射杆的长度lt=2 500 mm。入射杆上应变片与试样接触左端面的距离l1=1 500 mm, 透射杆上应变片与试样接触右端面的距离l2=1 000 mm。实验均在0.1 MPa加载气压下进行, 产生的压缩应力波由入射杆上应变片和透射杆上应变片测量, 信号经超动态应变仪进行记录后, 转化为电压信号并由瞬态记录仪进行采集, 分别为入射波εi(t)、反射波εr(t)和透射波εt(t)。
为了在动态加载的条件下有效地检测到裂纹扩展的完整过程, 采用裂纹扩展计(crack propagation gauge, CPG)检测装置[23-24], 并进行了改进, 使其既能完整地对整个断裂带进行检测, 又能准确地获得实验数据, 如图 4所示。同时, 在试样的特殊位置设计粘贴应变片(strain gauge, SG), 有利于实验数据和实验现象的分析。CPG具体尺寸及SG粘贴位置见图 4。
将CPG粘贴在裂纹前端, 为了研究应力波边界反射对裂纹动态断裂过程的影响, 在距离加载端面d=5 mm的位置每间隔5 mm粘贴3个应变片SG1、SG2、SG3。在试样另外一面裂纹尖端粘贴2 mm×2 mm的应变片SG4。实验选用CPG共有25根金属丝, 并联总电阻R=1.5 Ω。随着裂纹的扩展, 金属丝依次发生断裂, CPG的并联总电阻就会增加ΔR, CPG两端电压发生变化, 通过监测电压阶跃变化的起跳点即可确定裂纹扩展到该位置的时刻, 已知2根金属丝间距, 可以计算出2根金属丝间的裂纹扩展速度, 从而得到CPG检测范围内裂纹扩展的速度历程。
1.3 动态断裂全过程
图 5所示为动态实验中试样在第1次加载波作用下完全断裂时的入射杆上SGi和CPG电压信号U, 简称为一次断裂, 以CEHFD-2试样为例。SGi记录了加载应力波在入射杆中传播的3个来回, 1个来回即为一次应力波加载。
图 6所示为动态实验中试样在第2次加载波作用下才完全断裂时的入射杆上SGi和CPG的电压信号, 简称为二次断裂。以CEHFD-3试样为例。可知, 在CPG的检测范围内, 试样在第1次加载波的作用下起裂扩展, CPG金属丝并未完全断裂, 表明裂纹在CPG检测范围内的某一个点不再扩展, 即止裂。将这一阶段裂纹的起裂、扩展和止裂定义为阶段Ⅰ。随后, 当第2次加载波到达试样时, 试样再次受到加载波的作用, 裂纹再次起裂扩展, 剩余未断裂的金属丝全部断裂, 超出CPG的检测范围。将这一阶段定义为阶段Ⅱ。阶段Ⅰ的裂纹起裂-扩展-止裂即为完整的动态断裂全过程。
由于实验试样尺寸有限, 需要有较长裂纹扩展距离, 才有可能检测到裂纹的止裂现象, 因而检测止裂的前提条件是有足够的断裂带长度, 断裂带太短, 其裂纹扩展距离短, 导致裂纹尚在扩展阶段, 试样已经完全断裂, 且接近试样边界区域, 应力波和卸载波反射后较为复杂, 即使出现止裂, 若止裂点太靠近这一区域, 复杂应力波也将使止裂难以研究。同时, 较长的断裂带增大了裂纹扩展速度的检测范围, 有利于研究裂纹扩展阶段的速度变化历程。在距离试样透射端5 mm处等距离贴有3个应变片, 即SG1、SG2和SG3, 电压信号如图 7所示。应变片受到应力波作用, 产生应变, 检测到应力波最大值均在裂纹止裂以后, 可知, 裂纹止裂过程, 没有受到反射应力波的影响。止裂前的速度变化也没有受到反射应力波的干扰。
1.3.1 裂纹起裂和止裂时刻
将加载应力波开始作用于试样时的时刻定义为初始时刻t0, 且:
t0=ti+tr2 (1) 根据文献[25-26]中的方法确定入射波和反射波的波头, 分别为ti和tr。
由于CPG第1根金属丝距离裂纹尖端有一段距离, 即金属丝与基底前沿之间的距离l0, 第1根丝断裂对应的时刻t1并不是裂纹的起裂时刻, 而是裂纹扩展l0距离后第1根金属丝断裂的时刻, 如图 4所示。文献[23-24]在求起裂时刻时均假设裂纹起裂后以某一速度匀速扩展, 通过CPG求出这一平均速度va, 起裂时刻tf = t1-l0/va。
但是实验证明, 裂纹扩展的初始阶段并非是匀速的, 采用这一方法得到的起裂时刻有较大偏差, 本文中采用应变片破裂法确定试样的起裂时刻[27-28], 在试样裂纹尖端粘贴应变片SG4, 测得裂尖处的应变波形, 将应变波形对时间进行求导, 得到裂尖处应变率最大的时刻, 图 8所示为CEHFD-2试样和CEHFD-3试样的SG4记录的应变率˙ε信号, 其峰值所对应的时刻作为裂纹通过该应变片所在位置的时刻, 此时刻即为试样的起裂时刻tf[20]。
图 9(a)所示为CEHFD-2试样在阶段Ⅰ的CPG电压信号, 即图 5中该试样裂纹起裂-扩展过程, CPG的25根金属丝全部断裂, 试样一次断裂, 未有止裂。图 9(b)所示为CEHFD-3试样在阶段Ⅰ的CPG电压信号, 即图 6中该试样裂纹起裂-扩展-止裂过程, CEHFD-3试样在第18根丝后止裂。经计算得到CEHFD-3试样的t0=846.3 μs, 其它时间参量均以这一时刻作为起点计算。以t0作为起点, CEHFD-3试样的起裂时刻tf=(929.1-846.3) μs=82.8 μs, CPG断裂第1根丝断裂时刻为t1=103.6 μs, 阶段Ⅰ最后一根丝断裂时刻t18=309.8 μs, 这一时刻后裂纹止裂, 因此假设这一断裂点为裂纹的止裂点, 这一时刻即为止裂时刻ta。
在阶段Ⅱ, CEHFD-3试样裂纹在第2次加载波作用下再次起裂, 如图 10所示, 即图 6中该试样裂纹再起裂后CPG电压信号放大图。将再起裂时CPG第1根金属断裂的时刻规定为试样的再起裂时刻tf′, 直到第25根金属丝断裂, 超出CPG检测范围。以t0为起点, CEHFD-3试样的再起裂时刻t19 =tf′ = (2 643.3-846.3) μs =1 797.0 μs。应力波在SHPB中传播速度按理论一维纵波波速c0进行计算, 应力波在入射杆上传播一个来回需要的时间为2li/c0=1 740.1 μs, 2次起裂时间间隔为(tf-tf′)= (1 797.0-82.8) =1 714.2 μs, 2个时间有所差异, 说明再次起裂所用的起裂时间较短, 这是因为自然裂纹较之人工预制裂纹更为尖锐。动态实验数据见表 1。
表 1 动态实验数据Table 1. Dynamic experimental data试样编号 t0/μs t4/μs tf/μs ta/μs tf′/μs (tf-tf′)/μs CEHFD-1 840.9 916.9 76.0 253.6 1 781.8 1 705.8 CEHFD-2 846.4 928.2 81.8 CEHFD-3 846.3 929.1 82.8 309.8 1 797.0 1 714.2 CEHFD-4 838.0 923.7 85.7 327.5 1 804.2 1 718.5 CEHFD-5 844.1 916.9 72.8 282.8 1 793.0 1 720.2 CEHFD-6 844.0 923.5 79.5 252.6 1 774.6 1 695.1 1.3.2 裂纹扩展速度
CEHFD-2和CEHFD-3试样阶段Ⅰ裂纹在起裂后快速扩展, 其裂纹长度、裂纹扩展速度与时间关系曲线如图 11~12所示。这里做了一个假设, 假设裂纹在2根CPG两根金属丝间的速度是匀速且是直线扩展, 其每2根金属丝间的速度按下式计算:
v=Δl(tn+1−tn) (2) 式中:Δl为2根金属丝间的间距, Δl= l/24=1.67 mm。这里需要注意的是tf与t1之间裂纹扩展距离为l0=3 mm。
可知, 由于一次断裂情况下, 裂纹没有止裂过程, 其a(t)曲线可近似为直线, 如图 11(a)所示, 其裂纹速度增大到最大值后发生振荡, 如图 11(b)所示, 可近似将裂纹扩展阶段的运动看作匀速, 其平均速度为va=527.97 m/s (0.37cR)。二次断裂情况下的裂纹扩展在阶段Ⅰ, 经历了裂纹的起裂-扩展-止裂的过程, 随着时间的增加a(t)曲线并不是均匀变化, 如图 12(a)所示, 这就意味着不可将裂纹作匀速扩展进行处理。可以明显看到, 在裂纹开始起裂后, 裂纹加速上升达到最大速度, 在裂纹止裂前, 速度缓慢下降, 直至速度为零, 裂纹止裂, 如图 12(b)所示。CEHFD-3最大速度为vmax=666.67 m/s (0.47cR), 最小速度为vmin=42.19 m/s (0.03cR)。
2. 动态应力强度因子
实验-数值-解析法是一种将实验记录、数值模拟和解析公式相结合的研究方法[29], 综合考虑材料惯性效应和裂纹扩展速度对动态应力强度因子的影响, 该方法避开了动态应力强度因子难以直接进行测量和计算的问题, 将原本对实验设备的苛刻要求分担到动态实验、数值模拟和近似解析3方面, 使实验中的重点变为动态荷载历程、裂纹起裂止裂时刻和扩展速度等易于检测的信号的记录测量, 而动态应力强度因子等难以直接检测的力学参数则通过有限元等动态数值模拟分析, 并结合解析公式得到。
2.1 实验动态荷载
压缩空气炮驱动炮弹撞击SHPB入射杆, 入射杆上应变片SGi采集的信号, 即入射波εi(t); 弹性压缩波向前传播到达界面后一部分反射回入射杆, 再次通过应变片SGi时采集的信号, 即反射波εr(t)。另一部分弹性压缩波通过试样, 造成试样的高速变形后传播进入透射杆, 透射杆上应变片SGt采集的信号, 即透射波εt(t)。根据SHPB的一维弹性应力波假设, 由入射端的加载波波形可以得到入射端荷载pi(t), 由透射端的加载波波形可以得到透射端荷载pt(t):
{pi(t)=EbAb[εi(t)+εr(t)]pt(t)=EbAbεt(t) (3) 式中:Eb、Ab分别为SHPB弹性模量及横截面面积。如图 13所示得到CEHFD-3试样的pi(t)曲线。其左端动载荷pi(t)最大值pmax=37.98 kN, 时间tmax=44.7 μs。由于本次实验试样为尺寸较大的圆孔圆盘类试样, 应力波传播反射的路径长且次数多, 因而传播时间较长; 试样带有预制裂纹, 岩石动态断裂时会消耗能量, 作用于试样两端的载荷存在显著的时间与空间不均匀性, 即左端的加载载荷pi(t)与右端的加载载荷pt(t)在时间和数量上差异性较大。本文的研究方法为实验-数值-解析法, 不需要满足准静态法研究中的应力均匀性假设, 入射波εi(t)和反射波εr(t)是在近乎一维线弹性的SHPB压杆中传播的, 在荷载叠加方法合理的情况下, pi(t)更接近试样入射端的真实荷载[30], 将pi(t)作为实验的动态加载荷载p(t)。
2.2 数值-动态有限元模拟
高加载率情况下, 快速扩展裂纹尖端附近质点受到裂纹尖端应力集中区影响, 因此, 对于裂纹尖端的运动, 裂纹速度的影响很重要[10]。Freund[10]认为可以将裂纹速度对动态应力强度因子的影响分离出来, 对该时刻静止裂尖的动态应力强度因子时间历程进行普适函数的修正:
KI(d)[t,v(t),a(t)]=k(v)KI(0)[t,a(t)] (4) 式中:KⅠ(d)为动态应力强度因子, KⅠ(0)为静止裂尖的动态应力强度因子, 裂纹扩展速度v和裂纹长度a是时间t的函数, k(v)为普适函数。
采用ANSYS软件动态模拟静止裂纹受外加动荷载时, 其裂尖的动态应力强度因子为KⅠ(0)(t)。首先对经典的“Chen问题”[31]进行分析, 利用位移外推法[32], 计算出一定裂纹长度下的KⅠ(0)(t), 与Chen的有限差分结果基本一致, 验证了采用这一方法的有效性以及高准确率。
CEHFD试样1/2有限元模型如图 14所示, 采用PLANE82平面应变单元, 裂纹尖端应力场和应变场的奇异性采用1/4节点奇异单元, 模型共有3 860个单元和11 837个节点, 选用时间子步步长为0.1 μs。对试样端面进行位移耦合, 使其在加载过程中始终保持一个平面, 以模拟实验过程中SHPB压杆对试样端面的平面加载。起裂时裂纹尺寸a(tf)=a0= 8 mm; 扩展时裂纹尺寸:
a(tn)=a0+l0+(n−1)Δl (5) 可令:a(tn)=an, 止裂时裂纹尺寸a(ta)=aa。
将实验得到的动态载荷时间历程p(t)输入有限元计算程序, 设定其裂纹尺寸为t时刻时的a(t), 将得到载荷p(t)作用下静止裂纹的KⅠ(0)(t)。
2.3 解析-普适函数修正
普适函数是由Freund[10]基于Green函数基本解的概念, 认为在一般荷载作用下的Ⅰ型裂纹, 在不超过瑞雷波速cR时, 以任意速度扩展的动态应力强度因子等于适用于具有该处静止裂尖的动态应力强度因子与瞬时裂纹速度的普适函数的乘积。其近似计算公式为:
k(v)≈(1−v/cR)√1−v/cd (6) 式中:cR为材料的瑞雷波速, cd为材料的膨胀波速。
根据上一节得到静止裂纹的KⅠ(0)(t), 利用式(4)、(6)得到KⅠ(d)(t), 如图 15所示。采用实验-数值-解析法得到CEHFD-3试样整个断裂过程中的动态应力强度因子, 在图中以圆圈表示。实线为KⅠ(d)(t)的四次多项式拟合曲线, 可以看出, KⅠ(d)(t)随着裂纹的扩展呈现先增大后减小的趋势。
3. 结果分析
3.1 动态断裂韧度的确定
不考虑外界温度影响, 动态起裂准则为:
KI(d)(tf)=KIC(D)⋅(K) (7) 式中:KⅠ(d)(tf)表示t = tf时刻的动态应力强度因子, KⅠC(D)为动态起裂韧度, 是抵抗裂纹起裂的材料参数, ˙K=KIC(D)/tf表示动态加载率。当驱动力的KⅠ(d)(tf)超过材料阻力的KⅠC(D)时, 裂纹起裂。CEHFD-3试样的tf=82.8 μs, 根据图 15得到KⅠC(D)=2.65 MPa·m1/2, ˙K=3.20×104 MPa·m1/2/s。
动态扩展韧度与速度有关, 是抵抗动态裂纹扩展的材料参数。动态扩展准则为:
KI(d)(t)=KIC(d)(v) (8) 式中:KⅠ(d)(t)表示tf < t < ta时的动态应力强度因子, KⅠC(d)为动态扩展韧度。
由图 15可知, CEHFD-3试样扩展阶段的KⅠ(d)(t)即为KⅠC(d), KⅠC(d)随着时间的增加先增大后减小。KⅠC(d)与速度有关, 而裂纹扩展速度v(t)为时间的函数, 见图 12(b), 在裂纹快速扩展阶段, CEHFD-3试样裂纹扩展速度在起裂后加速上升, 止裂前减速下降, 速度随着时间先增加后减小。而对于未止裂的CEHFD-2试样, 见图 11(b), 其裂纹扩展速度加速上升后在平均速度va上下波动, 得到的KⅠ(d)(t)也在某一平均值上下波动, CEHFD-2试样的KⅠC(d)值为一常数, KⅠC(d)=3.37 MPa·m1/2。
Freund[10]和Ravi-Chandar[33]提出动态止裂韧度定义:不能维持裂纹继续扩展的最大动态应力强度因子即为动态止裂韧度。动态扩展准则为:
KI(d)(ta)≤KIC(a) (9) 式中:KⅠ(d)(ta)表示t=ta时刻的动态应力强度因子, KⅠC(a)为动态止裂韧度, CEHFD-3试样ta=309.8 μs, 根据图 15得到KⅠC(a)=0.69 MPa·m1/2。
砂岩动态起裂韧度和止裂韧度见表 2, 可知, 动态起裂韧度与加载率有关, 随着加载率的增大而增大; 动态止裂韧度小于动态起裂韧度, 且止裂韧度随着裂纹最大扩展速度的增加而降低, 并且有较大的离散性。
表 2 砂岩的动态起裂和止裂断裂韧度Table 2. Dynamic initiation and propagation toughness of sandstone试样编号 ˙KⅠ/(GPa·m1/2·s-1) KⅠC(d)/(MPa·m1/2) vmax/cR KⅠC(a)/(MPa·m1/2) CEHFD-1 27.4 2.08 0.44 0.72 CEHFD-2 36.1 2.95 0.51 CEHFD-3 32.0 2.65 0.47 0.69 CEHFD-4 41.5 3.56 0.35 1.52 CEHFD-5 24.3 1.77 0.40 1.07 CEHFD-6 26.9 2.14 0.51 0.61 3.2 裂纹止裂过程
假定材料性能值KⅠC或(GⅠC), 即为“盾”; 结构提供的裂纹驱动力K (或G), 即为“矛”。裂纹的力学行为决定于矛盾两体的较量, 当“盾”作用大于“矛”作用时, 显然没有足够的驱动力使裂纹扩展。一般来说, 止裂的量度有两种[10]:一种是从结构的角度考虑。“盾”不变, 而“矛”随着裂纹尺寸的增大其作用减弱, “矛”小于“盾”时产生止裂。另一种是从材料的角度考虑。“盾”随着裂纹尺寸的增大而增大, 当“盾”明显超过“矛”时实现止裂。
通过实验-数值-解析的方法, 得到了CEHFD-3试样随时间变化的动态应力强度因子, 如图 15所示, 并且通过拟合, 得到KⅠ(d)(t)变化趋势曲线, 可视作其裂纹动态扩展的驱动力。KⅠ(d)(t)值变化曲线如图 16所示。假设阻力R为常数, 随着裂纹的扩展, 驱动力的KⅠ(d)(t)先增大后减小。
在裂纹开始扩展之后, 在B点以前KⅠ(d)是增加的, 随后开始减小。到达C点时, KⅠ(d)重新等于阻力R, 如果忽略动能不计, 则该时刻可能发生止裂。此时裂纹扩展长度为a1。C点后有KⅠ(d)小于R, 裂纹仍然能够扩展, 因为在C点裂纹扩展速度不为零, 动能不可忽略, 这部分能量可用于裂纹扩展, 这时观察到的现象就是裂纹速度的降低, 最后在E点动能减小为零, 从而裂纹扩展速度为零, 在E点止裂, 裂纹最终扩展长度为a2。因为瞬时的KⅠ(d)值比KⅠC更小, 因此在速度变为零后, 裂纹止裂, 这一较小的KⅠ(d)值即为动态止裂韧度KⅠC(a)。因此, 不能将裂纹止裂看作是裂纹起裂在时间上的逆过程。Ravi-Chandar指出裂纹止裂是不考虑加载条件的瞬时动态问题, 不能简单的将其看作裂纹起裂的逆过程[33]。
实际当中, 对某些材料, R通常不是常数且与裂纹速度有关, 而快速奔跑裂纹的所有动能都用于裂纹扩展也不太可能。对于许多材料来说, K和R都与裂纹长度和速度有关, 止裂时刻的K并不是一个材料常数[33]。实验结果止裂韧度的离散性也证明了这一点。
4. 结论
(1) 为了更好地研究岩石材料的Ⅰ型动态断裂行为, 特别是裂纹的快速扩展与止裂过程中裂纹尖端的运动情况, 本文提出一种利用岩芯的新型试样CEHFD, 成功地检测到裂纹的动态起裂-扩展-止裂这一完整的动态断裂过程。
(2) 实验表明, 在整个断裂过程中, 裂纹并非匀速扩展, 裂纹扩展速度在裂纹起裂后加速上升, 在裂纹止裂前有明显的减速现象, 裂纹并非突然止裂。在试样透射端附近粘贴应变片, 检测到应力波最大值均出现裂纹止裂以后, 可知, 止裂前的减速过程并不是由反射应力波造成的, 裂纹止裂过程并没有受到反射应力波的影响。
(3) 采用实验-数值-解析法得到动态应力强度因子KⅠ(d)(t), 其时间历程呈现先增大后减小的趋势, 并得到砂岩的动态起裂韧度、扩展韧度及止裂韧度。在完整的动态断裂全过程中, 裂纹扩展速度随着时间先增大后减小, 动态扩展韧度为速度的函数, 变化趋势与速度一致, 随着时间先增大后减小; 动态起裂韧度大于动态止裂韧度, 且止裂韧度随着裂纹最大扩展速度的增大而降低, 并且有较大的离散性。
-
表 1 动态实验数据
Table 1. Dynamic experimental data
试样编号 t0/μs t4/μs tf/μs ta/μs tf′/μs (tf-tf′)/μs CEHFD-1 840.9 916.9 76.0 253.6 1 781.8 1 705.8 CEHFD-2 846.4 928.2 81.8 CEHFD-3 846.3 929.1 82.8 309.8 1 797.0 1 714.2 CEHFD-4 838.0 923.7 85.7 327.5 1 804.2 1 718.5 CEHFD-5 844.1 916.9 72.8 282.8 1 793.0 1 720.2 CEHFD-6 844.0 923.5 79.5 252.6 1 774.6 1 695.1 表 2 砂岩的动态起裂和止裂断裂韧度
Table 2. Dynamic initiation and propagation toughness of sandstone
试样编号 ˙KⅠ/(GPa·m1/2·s-1) KⅠC(d)/(MPa·m1/2) vmax/cR KⅠC(a)/(MPa·m1/2) CEHFD-1 27.4 2.08 0.44 0.72 CEHFD-2 36.1 2.95 0.51 CEHFD-3 32.0 2.65 0.47 0.69 CEHFD-4 41.5 3.56 0.35 1.52 CEHFD-5 24.3 1.77 0.40 1.07 CEHFD-6 26.9 2.14 0.51 0.61 -
[1] BRIAN C, YANG Qingda.In quest of virtual test for structural composite[J].Science, 2006, 314(5802):1102-1107.DOI: 10.1126/science.1131624. [2] 余寿文.断裂力学的历史发展和思考[J].力学与实践, 2015, 39(3):390-394.DOI: 10.6052/1000-0879-15-164.YU Shouwen.The historical development of fracture mechanics and rethinking[J].Mechanics in Engineering, 2015, 39(3):390-394.DOI: 10.6052/1000-0879-15-164. [3] DAS S.Dynamic fracture mechanics in the study of the earthquake rupturing process:Theory and observation[J].Journal of the Mechanics and Physics of Solids, 2003, 51(11/12):1939-1955.DOI: 10.1016/j.jmps.2003.09.025. [4] DAS S.The need to study speed[J].Science, 2007, 317(5840):905-906.DOI: 10.1126/science.1142143. [5] ROBINSON D P, BROUGH C, DAS S.The Mw 7.8, 2001 Kunlunshan earthquake:Extreme rupture speed variability and effect of fault geometry[J].Journal of Geophysical Research, 2006(111):B08303.DOI: 10.1029/2005JB004137. [6] KALTHOFF J F, BEINERT J, WINCKLER S.Experimental analysis of dynamic effects in different crack arrest test specimens[J].ASTM STP 1980, 711:109-127.DOI: 10.1520/STP27443S. [7] KOBAYASHI T, DALLY J W.In fast fracture and crack arrest[J].ASTM STP 1977, 627:257-273. doi: 10.1016-0013-7944(86)90191-8/ [8] SCHARDIN H.Velocity effects in fracture[M].1959: 297-330. [9] RAVI-CHANDAR K, KNAUSS W G.An experimental investigation into dynamic fracture:Ⅰ.Crack initiation and arrest[J].International Journal of Fracture, 1984, 25(4):247-262. doi: 10.1007/BF00963460 [10] FREUND L B.Dynamic fracture mechanics[M].Cambridge:Cambridge University Press, 1990. [11] DALLY J W.Dynamic photo elastic studies of fracture[J].Experimental Mechanics, 1979, 19(10), 349-361. doi: 10.1007/BF02324250 [12] KALTHOFF J F.On some current problems in experimental fracture dynamics[C]//Workshop on Dynamic Fracture.Pasadena: California Institute of Technology, 1983: 11-35. [13] RAVI-CHANDAR K, KNAUSS W G.An experimental investigation into dynamic fracture:Ⅲ.On steady state crack propagation and crack branching[J].International Journal of Fracture, 1984, 26(2):141-154.DOI: 10.1007/BF01157550. [14] FINEBERG J, BOUCHBINDER E.Recent developments in dynamic fracture:some perspectives[J].International Journal of Fracture, 2015, 196(1):33-57.DOI: 10.1007/s10704-015-0038-x. [15] ARAKAWA K, MADA T.Unsteady dynamic crack propagation in a brittle polymer[J].Experimental Mechanics, 2007, 47(5):609-615.DOI: 10.1007/s11340-006-9020-x. [16] SINGH R P, KAVATURU M, SHUKLA A.Initiation, propagation and arrest of an interface crack subjected to controlled stress wave loading[J].International Journal of Fracture, 1997, 83(3):291-304. doi: 10.1023/A:1007358901588 [17] 王自强, 陈少华.高等断裂力学[M].北京:科学出版社, 2009. [18] LEVENT T, CIGDEM K.Effects of geometric factors on mode Ⅰ fracture toughness for modified ring tests[J].International Journal of Rock Mechanics and Mining Sciences, 2012, 51(1):149-161. https://www.deepdyve.com/lp/elsevier/effects-of-geometric-factors-on-mode-i-fracture-toughness-for-modified-xPK0koo9Zh [19] IUNG T, PINEAU A.Dynamic crack propagation and crack arrest investigated with a new specimen geometry:Part Ⅰ:Experimental and numerical calculations[J].Fatigue and Fracture of Engineering Materials and Structures, 1996, 19(11):1357-1367. doi: 10.1111/ffe.1996.19.issue-11 [20] 杨井瑞, 张财贵, 周妍, 等.用CSTBD试样确定砂岩的动态起裂和扩展韧度[J].爆炸与冲击, 2014, 34(3):264-271. http://www.bzycj.cn/CN/abstract/abstract8838.shtmlYANG Jingrui, ZHANG Caigui, ZHOU Yan, et al.Determination of dynamic initiation toughness and propagation toughness of sandstone using CSTBD specimens[J].Explosion and Shock Waves, 2014, 34(3):264-271. http://www.bzycj.cn/CN/abstract/abstract8838.shtml [21] 张盛, 王启智, 谢和平.岩石动态断裂韧度的尺寸效应[J].爆炸与冲击, 2008, 28(6):544-551.DOI: 10.11883/1001-1455(2008)06-0544-08.ZHANG Sheng, WANG Qizhi, XIE Heping.Size effect of rock dynamic fracture toughness[J].Explosion and Shock Waves, 2008, 28(6):544-551.DOI: 10.11883/1001-1455(2008)06-0544-08. [22] 周妍, 张财贵, 杨井瑞, 等.圆孔内单边(或双边)裂纹平台巴西圆盘应力强度因子的全面标定[J].应用数学和力学, 2015, 36(1):16-30. http://d.old.wanfangdata.com.cn/Periodical/yysxhlx201501002ZHOU Yan, ZHANG Caigui, YANG Jingrui, et al.Comprehensive calibration of stress intensity factor for flattened Brazilian disc with holed inner single or double cracks[J].Applied Mathematics and Mechanics, 2015, 36(1):16-30. http://d.old.wanfangdata.com.cn/Periodical/yysxhlx201501002 [23] 周妍, 张财贵, 王启智.用圆孔内单边裂纹平台巴西圆盘和实验-数值-解析法确定砂岩的动态起裂和扩展韧度[J].振动与冲击, 2017, 36(5):37-47. http://d.old.wanfangdata.com.cn/Periodical/zdycj201705007ZHOU Yan, ZHANG Caigui, WANG Qizhi.Determination of dynamic initiation toughness and dynamic propagation toughness of sandstone[J].Journal of Vibration and Shock, 2017, 36(5):37-47. http://d.old.wanfangdata.com.cn/Periodical/zdycj201705007 [24] 李炼, 杨丽萍, 曹富, 等.冲击加载下的砂岩动态断裂全过程的实验和分析[J].煤炭学报, 2016, 41(8):1912-1922. http://d.old.wanfangdata.com.cn/Periodical/mtxb201608006LI Lian, YANG Liping, CAO Fu, et al.Complete dynamic fracture process of sandstone under impact loading:Experiment and analysis[J].Journal of China Coal Society, 2016, 41(8):1912-1922. http://d.old.wanfangdata.com.cn/Periodical/mtxb201608006 [25] 唐志平, 王礼立.SHPB实验的电脑化实验处理系统[J].爆炸与冲击, 1986, 6(4):320-327. http://www.bzycj.cn/CN/abstract/abstract11099.shtmlTANG Zhiping, WANG Lili.A computerized system of data processing used in SHPB experiments[J].Explosion and Shock Waves, 1986, 6(4):320-327. http://www.bzycj.cn/CN/abstract/abstract11099.shtml [26] 李夕兵, 古德生, 赖海辉.冲击荷载下岩石动态应力-应变全图测试中的合理加载波形[J].爆炸与冲击, 1993, 13(2):125-130. http://www.bzycj.cn/CN/abstract/abstract10667.shtmlLI Xibing, GU Desheng, LAI Haihui.On the reasonable loading stress waveforms determined by dynamic stress-strain curves of rocks by SHPB[J].Explosion and Shock Waves, 1993, 13(2):125-130. http://www.bzycj.cn/CN/abstract/abstract10667.shtml [27] GUO W G, LI Y L, LIU Y Y.Analytical and experimental determination of dynamic impact stress intensity factor for 40Cr steel[J].Theoretical and Applied Fracture Mechanics, 1997, 26(1):29-34. doi: 10.1016/S0167-8442(96)00031-6 [28] WEISBROD G, RITTLE D.A method for dynamic fracture toughness determination using short beams[J].International Journal of Fracture, 2000, 104(1):89-103. doi: 10.1023/A:1007673528573 [29] WANG Q Z, YANG J R, ZHANG C G, et al.Determination of dynamic crack initiation and propagation toughness of a rock using a hybrid experimental-numerical approach[J].Journal of Engineering Mechanics, 2016, 142(12):04016097.DOI: 10.1061/(ASCE)EM.1943-7889.0001155. [30] 杨井瑞, 张财贵, 周妍, 等.用SCDC试样测试岩石动态断裂韧度的新方法[J].岩石力学与工程学报, 2015, 34(2):279-292. http://d.old.wanfangdata.com.cn/Periodical/yslxygcxb201502007YANG Jingrui, ZHANG Caigui, ZHOU Yan, et al.A new method for determining dynamic fracture toughness of rock using SCDC specimens[J].Chinese Journal of Rock Mechanics and Engineering, 2015, 34(2):279-292. http://d.old.wanfangdata.com.cn/Periodical/yslxygcxb201502007 [31] CHEN Y M, WILKINS M L.Stress analysis of problems with a three-dimensional, time-dependent computer program[J].International Journal of Fracture, 1976, 12(4):607-617.DOI: 10.1007/BF00034646. [32] CHEN L S, KYANG J H.A modified linear extrapolation formula for determination of stress intensity factors[J].International Journal of Fracture, 1992, 54(1):R3-R8.DOI: 10.1007/BF00040859. [33] RAVI-CHANDAR K K.Dynamic fracture[M].Elsevier, 2004. 期刊类型引用(2)
1. 周磊,朱哲明,董玉清,牛草原,邓帅,姜亚成. 动载荷作用下含倒U型孔洞裂纹试件的动力损伤行为研究. 岩石力学与工程学报. 2021(01): 113-125 . 百度学术
2. 刘宁,胡梦凡,周飞. 基于键基近场动力学理论的单裂纹圆孔板冲击破坏研究. 工程力学. 2020(12): 9-17 . 百度学术
其他类型引用(4)
-