温度作用下花岗岩断裂行为损伤力学分析.docVIP

温度作用下花岗岩断裂行为损伤力学分析 　　由于岩石是含有微裂隙、微孔洞等初始缺陷的天然材料,当受热或承受一定载荷后,不可避免地会在其中产生大量细观裂纹,并随着温度的升高或载荷的增大而逐渐扩展,在一定程度上表现出材料受力性能的劣化,说明温度、载荷对岩石造成了损伤[1-4]。目前,国内外学者对温度影响下岩石损伤特性做了大量的研究工作。刘泉声等[5]从花岗岩弹性模量随温度的变化规律入手,提出了热损伤的概念,在此基础上导出了热损伤演化方程和一维热力耦合弹脆性损伤本构方程。徐燕萍等[6]推导了岩石热弹塑性力学特性本构方程。付文生等[7]提出采用两个损伤变量来描述各向同性材料的损伤变化规律。王利等[8]定义了弹性应变表示的一维损伤变量及其本构模型,利用双剪强度理论将其推广至三维模型。但从力学意义上讲,损伤变量选取的一个重要原则是容易与宏观力学量建立联系并且易于测量。声发射是材料在外力作用下,其组成微粒产生错位或裂纹等释放出能量而产生的弹性波。声发射检测直接反映了试样内部的损伤,因而声发射累积数可用于表征试样的损伤度[9-11]。因此本文基于岩样在不同温度作用下的力学实验和声发射实验,讨论了由于载荷作用引起的机械损伤和由于温度作用引起的热损伤,建立了热力耦合损伤本构方程,同时研究了花岗岩热损伤开裂机理。 　　1　实验 　　1.1　方法 　　实验所采用的岩样为产于某矿区的花岗岩,按岩石力学性质试验规范加工成25 mm×50 mm的圆柱体。 　　温度自室温开始,以2℃/s的速率匀速加温到50,100,200,300,400,500,600,700,800,900,1 000,1 100,1 200℃,恒温20 min后自然冷却至常温。实验采用位移控制方式,以0.001 5 mm/s的加载速率一次加载至破坏,实验过程中测量岩样载荷、变形及声发射参数等。 　　1.2　设备 　　温度的施加采用高温炉MTS652.02实现,加载设备采用MTS815电液伺服材料试验系统,声发射监测采用由美国声学物理公司生产的PCI-2声发射仪。 　　2　结果与分析 　　2.1　温度作用下花岗岩力学性质实验结果 　　由于每组岩样的应力-应变曲线具有相近的分布形式,大致均经历了压密、弹性、弹塑性变形和破坏4个阶段。该文列出若干具有代表性的应力-应变曲线,如图1所示[12]。实验所得花岗岩的弹性模量随温度的变化规律如图2所示。 　　随着温度的升高,花岗岩的脆性减弱,延性增强,800℃之前,尽管岩样脆性降低,强度变小,但是变化不大,温度超过800℃,岩样强度突然降低,呈现出明显的塑性行为和峰后行为,1 200℃岩样基本失去承载能力。花岗岩各温度段弹性模 量具有较大的离散性,这是由于岩石自身的结构差异所致。按最小二乘法拟合的回归曲线表现为三次抛物线型,弹性模量随温度升高总体呈下降趋势。 　　2.2不同温度作用下花岗岩的变形与声发射率 　　选取声发射振铃计数率对岩样的声发射特性进行描述,以此来反映岩样内部破裂事件发生的频率。绘制不同温度作用后花岗岩的应力-应变全过程曲线与声发射曲线如图3所示。 　　在各个温度点,花岗岩应力-应变曲线与声发射曲线吻合较好,声发射总是发生在加载曲线的转折点,转折点意味着能量的突变,这与声发射是岩体内部的能量突然释放相一致。800℃前声发射信号呈现出明显的脆性破坏特征,承受800℃高温后,岩样呈现塑性破坏特征,峰值强度之后,试件的残余塑性变形仍释放出较密集的声发射信号,与岩样在800℃后宏观力学性能突然降低相一致。 　　3　损伤分析 　　3.1　由载荷作用产生的机械损伤 　　由于岩石微缺陷的演化是一种随机变化,可以认为微元的强度分布服从Weibull分布,其分布密度函数为[13]:(ε)=m/αεm-1exp(-εm/α),此处(ε)是岩样在加载过程中微元损伤率的一种度量。从本质上看,由于声发射活动规律是一种统计的规律,因此,它必然与材料内部缺陷的统计分布规律一致。若整个截面全破坏时的声发射积累为Ωm,则岩样受压应变增至时的声发射积累为:Ω=Ωm∫ε0(x)dx,积分后得: 　　Ω/Ωm= 1-exp(-εm/α) (1)损伤参量D是材料损伤程度的度量,与微元破坏的概率密度之间存在着如下关系: 　　D=∫ε0(x)dx= 1-exp(-εm/α) (2) 　　比较式(1)和式(2),得到由于载荷作用而产生的机械损伤为:D=Ω/Ωm。 　　根据声发射实验数据绘制岩样在不同温度下的机械损伤曲线,如图4所示。 　　可以看出,该文提出的损伤变量D能够较好地反映花岗岩的微观损伤机制。在初始压密阶段,花岗岩内部的原生孔隙被压密,损伤变量为0;随后损伤变量小幅度增大,这一阶段的变形主要是由固体颗粒本身的弹性变形引起的;当应力大约达到极限强度的60%时,即σ=0