三轴应力-应变应力过程中花岗岩力学特性的试验研究.docxVIP

三轴应力-应变应力过程中花岗岩力学特性的试验研究 1 应变软化本构模型 花岗岩是大陆地壳中最重要的岩石之一。结构均匀，质地坚硬，颜色优美，在中国的土木工程建设中发挥着重要作用。地壳中的岩石均处于三向应力状态，因而研究三向应力状态下岩石强度和变形特性，对于采矿、地下建筑、水利水电等与岩石相关的工程具有重要意义。 岩石在应力达到峰值强度之后，随着变形的继续增加，其强度迅速降到一个较低的水平，这种由于变形引起的岩石材料性能劣化的现象称之为“应变软化”。长期以来，有关岩石应变软化本构模型的研究一直是个热点。日本学者川本眺万最先提出了线形应变软化模型，在岩土工程中应用较多。方德平等应用线性应变软化模型，考虑岩石脆塑性过渡特性，推导了不同围压作用下出现的软–弹、残–软–弹、残–塑–弹三种情况下的围岩位移和应力公式，并进行了有限元计算。K.Dems和Z.Morz提出了应变软化的理想弹脆塑性模型，是对前者的进一步简化，能够反映脆性比较明显的岩土材料的软化现象。 关于应变软化数值模型的建立，沈新善等给出了完整的弹塑性本构积分的数值格式和增量型弹塑性有限元求解算法。从而有了应变空间中应变软化材料与弹塑性硬化材料通用的本构理论；杨超等运用Hoek提出的由主应力圆包络线确定黏聚力和内摩擦角等效值的方法和曲线拟合的方法，研究了围压对软岩峰后软化特性的影响和软岩的宏观物性参数峰后应变软化规律。 在前人的基础上，本文以研究花岗岩峰后力学特性为目的，从花岗岩常规三轴压缩试验出发，得到不同围压下全过程应力–应变曲线。利用轴向应力与轴向应变、侧向应变的关系拟合峰后屈服面并建立与内变量c相关的应变软化模型。 2 试验结果和分析 试验是在中国科学院武汉岩土力学研究所MTS 815.03电伺服岩石力学试验机上进行的，该试验系统配备轴压、围压和孔隙水压3套独立的闭环伺服控制系统，其主要特点是闭环电液伺服控制，整体刚度大。 花岗岩试样取自三峡地下厂房区，其矿物成分为黑云母20%，石英25%，长石45%，透闪石5%。试样按照国际岩石力学学会的建议方法和岩石力学试样规范，加工成直径为50 mm、高度为100 mm的圆柱体标准试样。 试验采用应变控制方式施加轴向压力，加载速率为1.0×10-5s-1。三轴压缩试验中，围压为0～40MPa。图1所示为试验中得到的有代表性花岗岩三轴压缩应力–应变全过程曲线(曲线上的数值为施加的围压)。 测试结果表明，在低围压下，花岗岩峰后应力发生脆性跌落，随着围压的升高，花岗岩的延性特征逐渐明显，峰后的塑性变形能力增强，峰值强度和残余强度都有相应的提高，并且呈现一定的规律性。需要指出的是，对于40 MPa围压下的试验曲线，由于其岩样在端部存在局部的沉积结构面强度缺陷，当缺陷材料达到其承载极限开始屈服时，导致应力–应变关系偏离了通常的形状。 本次试验中共进行了12块花岗岩的应力–应变全过程曲线测试，其中单轴压缩3块，三轴压缩9块，单轴抗压强度为150～177 MPa，残余强度与峰值强度之比在10%左右。如图2所示，与任放和盛谦的测试结果较为接近。 3 花岗岩峰后力学特性的研究 3.1 应变软化模型 为了研究岩石的峰后力学特性，有些学者将岩石的屈服破坏过程以图3的连续的应变软化模型来表征，图3的曲线ABCD表示理想的三线性应变软化模型。 3.2 应变软化参量和屈服应力 本文使用应变软化参量εps作为塑性内状态变量，将三峡花岗岩看作各向同性材料，并假设它的塑性硬化性质也是各向同性的。因而它的屈服面的形式为 式中：σ1为第一主应力；σ2，σ3分别为第二、三主应力且σ2=σ3；εps为应变软化参量，且有 式中：ε1p，ε3p均为塑性主应变分量。 对于常规三轴试验，从试验中得到的是不同围压σ3下差应力σ1-σ3对轴向应变ε1和侧向应变ε3的曲线。由于在现有的塑性力学理论中，大都假设卸载过程是线弹性的，如图3所示，用斜直线(li,i=1,2,3)代表卸载路径，直线的斜率为弹性模量，与横坐标的交于一点，即为岩石的塑性应变。用此方法可以得到不同围压下(σ2=σ3=aj,j=1,2,3)，差应力对不同塑性应变分量ε1p(ε1p=bk,k=1,2,3)和ε3p的曲线。从这些曲线可以得到差应力和应变软化参量εps的关系。 最后可将相应于不同εps值的屈服应力值列成表格。表1所示为三轴压缩试验的4种围压下εps=4.38×10-3时的屈服应力值。 3.3 不同屈服强度的响应面试验结果 假设花岗岩的屈服面形式满足莫尔–库仑屈服准则，则根据表1可以拟合出当εps=4.38×10-3时的后继屈服面，对于不同的εps可以拟合出不同的屈服面，继而得到屈服面的形式为 其中， 图4给出了相应于εps=3.13×10-4，εps=4.38×10-4和εps=9.04×10-4的3条屈服曲