Research on Damage Model of Rock Under Freeze-Thaw Cycles Based on Maximum Tensile Strain Criterion
-
摘要:
为反映寒区遭受冻融的岩石在外力作用下的应力-应变全过程，基于统计损伤力学理论，假设岩石微元体强度服从Weibull分布，结合最大拉应变破坏准则，建立了考虑冻融和荷载耦合作用的寒区冻融岩石损伤模型；推导了模型参数的理论解，利用已有的试验结果和已有损伤本构对模型进行了验证，探讨了不同冻融次数下岩石的总损伤曲线演化规律，并对模型参数进行了分析；采用数值模拟方法计算了冻融循环对隧道工程稳定性的影响. 结果表明：本文所建立的模型可较好地重现岩石的应力-应变全过程，并能反映岩石的峰后强度；不同冻融次数下岩石的总损伤演变曲线呈“S”形，总损伤演变曲线可分为3个阶段，即初始损伤阶段、加速阶段和完全损伤阶段；Weibull分布参数$ m $和$ {f}_{0} $分别代表岩石的脆性特征和塑性特征；随着冻融次数的增加，隧道围岩的垂直位移、最大主应力及衬砌最大主应力均逐渐增加，但隧道围岩的应力增幅较小，经历10次冻融循环后，隧道拱顶的围岩下沉量和底部的围岩隆起量分别增加17.87%和19.24%，隧道围岩最大拉应力和压应力极值分别增加了2.70%和2.01%，隧道衬砌最大拉应力和压应力极值分别增加了21.52%和17.87%.
Abstract:In order to reflect the whole process of stress-strain of rock subjected to freeze thawing in the cold region under external force, this paper uses the theory of statistical damage mechanics and assumes that the strength of rock microelements obeys Weibull distribution. Then, according to the maximum tensile strain failure criterion, a damage model of rock subjected to freeze thawing in the cold region is established by considering the coupling effect of freeze thawing and load. The theoretical solution of the model parameters is deduced, and the damage model is verified by the test results and damage constitutive model of the predecessors. The evolution law of the total damage curve of the rock under different freeze-thaw cycles is discussed, and the model parameters are analyzed. Finally, the influence of the freeze-thaw cycle on the stability of tunnel engineering is calculated by the numerical simulation method. The results show that the model established in this paper can well reproduce the stress-strain curve of the rock and reflect the post-peak strength of the rock. The total damage evolution curve of the rock under different freeze-thaw cycles is S-shaped, and it can be divided into three stages: initial damage stage, accelerated stage, and complete damage stage. Weibull distribution parameters
m and $ {f}_{0} $ represent the brittle and plastic characteristics of rock, respectively. Additionally, with the increase in freeze-thaw cycles, the vertical displacement and the maximum principal stress of the tunnel surrounding rock increase. Meanwhile, the maximum principal stress of tunnel lining also increases. However, the stress of the tunnel surrounding rock increases slightly. After 10 freeze-thaw cycles, the settlement of surrounding rock at the top of the tunnel arch and the uplift of surrounding rock at the bottom increase by 17.87% and 19.24%, respectively; the maximum tensile stress and compressive stress extremes of the tunnel surrounding rock increase by 2.70% and 2.01%, respectively, while those of the tunnel lining increase by 21.52% and 17.87%, respectively.-
Key words:
- rock /
- freeze-thaw cycle /
- failure criteria /
- damage model /
- parameter analysis
-
据相关数据统计，永久性冻土区和季节性冻土区约占我国国土面积的3/4，当环境温度低于0 ℃时寒区岩体内水分冻结，温度高于0 ℃时则消融. 随着季节变换及昼夜更替环境温度在0 ℃上下震荡，岩体内水分发生周期性的冻结和消融，即冻融循环[1-2]. 随着“一带一路”倡议的实施，大量的工程建设在西部寒区兴起，越来越多的岩土工程遭遇到冻融问题，其中岩石在冻融循环作用下的损伤模型一直是研究者关注的问题，对寒区岩石的损伤劣化机理及岩土工程稳定控制具有重要的理论意义[3-6].
目前，对冻融循环引起的岩石损伤劣化机理已经有了较为一致的认识：岩石矿物颗粒在低温作用下体积收缩，而孔隙和裂隙内的水冻结成冰产生约9.08%的体积膨胀，由于不同矿物颗粒缩胀率有所差异，导致颗粒间的缩胀变形不协调，矿物颗粒与微孔隙之间产生巨大的冻胀力，对岩石造成一定的损伤[7]；温度升高时，岩石内孔隙、裂隙中的冰体融化，随着冻胀力的释放及水分的迁移，进一步促进岩石的损伤；在长期的冻融循环作用下岩石内部微裂纹不断萌生和发展，造成岩石承载力下降[8-9].
迄今为止，国内外学者从宏观、微观及细观等多方面对冻融岩石的损伤劣化机理进行了大量研究：在冻融岩石损伤模型方面，张慧梅等[10]基于岩石内部缺陷分布的随机性，运用损伤力学理论建立了温度-荷载耦合作用下岩石的损伤模型，并通过试验验证模型的合理性；Huang等[11]假设岩石微元体强度服从Weibull分布，以弹性模型作为损伤变量结合最大拉应变准则建立了冻融岩石的损伤模型，对模型进行了验证，并应用于寒区隧道的稳定性分析；Jia等[12]将冻融作用视为疲劳损伤，建立了冻融岩石的疲劳损伤模型，揭示了反复冻融下岩石力学强度衰减机理；为实现对冻融岩石变形全过程的预测，袁超等[13]选用Drucker-Prager准则，建立了能反映冻融岩石各变形阶段特征的损伤模型；张峰瑞等[14]对经历不同冻融循环次数的花岗岩开展细观特征分析和剪切蠕变试验，基于试验结果提出了冻融岩石损伤黏性元件，构建了花岗岩冻融剪切蠕变本构模型.
总体看来，对岩石的冻融损伤模型研究较少，模型多基于莫尔-库伦准则及Drucker-Prager准则，不能准确反映岩石的张拉破坏特性和峰后强度特征. 众所周知，岩石的抗拉强度远低于抗压强度，岩石微孔隙内水结成冰形成的冻胀力可视为拉应力，岩石在冻融和受荷作用下易发生张拉破坏. 因此，本文基于最大拉应变破坏准则和统计损伤理论建立了考虑冻融和荷载耦合作用的岩石损伤模型，推导了模型参数的理论解，利用已有研究结果对本文建立的损伤模型进行了验证，探讨了不同冻融次数下岩石的总损伤曲线演化规律，对模型参数进行了分析，并采用数值模拟方法计算了冻融作用对隧道工程稳定性的影响规律.
1. 损伤本构模型建立
1.1 损伤变量及本构方程
基于Lemaitre应变等价准则[15]，岩石在荷载作用下的有效应力为
σi1=σi1−D， (1) 式中：${\sigma }_{{{\rm{i}}}}$为名义应力；$D$为受荷损伤变量.
天然岩石中含有大量微孔隙和裂隙，为了表征岩石微元体的强度引入Weibull分布函数，则岩石微元体的强度概率密度函数可表示为
P(ˉf)=mf0(ˉff0)m−1exp(−(ˉff0)m)， (2) 式中：$\bar{f}\mathrm{为}$微元体强度变量；$ {f}_{0} $、$ m\mathrm{为} $分布函数参数.
受荷损伤变量D为失效单元数${N}_{{\rm{f}}}$与总单元数$ {N}$的比值，如式（3）.
D=NfN. (3) 区间[0，$\bar{f}$]范围内失效单元数可由式（4）计算.
Nf(ˉf)=∫ˉf0NP(y)dy=N(1−exp(−(ˉff0)m))， (4) 则可得到受荷损伤变量$ \mathrm{为} $
D=1−exp(−(ˉff0)m). (5) 参照张慧梅等[16]提出的冻融损伤模型，将冻融作用下承受荷载的岩石划分为冻融损伤部分、受荷损伤部分和未损伤部分，则
σi=σi1A1+σr(An+A2)， (6) A=A1+An+A2， (7) 式中：${\sigma }_{{\rm{r}}}$为损伤部分承受的残余应力；$A$为受荷面积；${A}_{1}$为未损伤部分面积；${A}_{{\rm{n}}}$为冻融损伤部分面积；${A}_{2}$为受荷损伤部分面积.
岩石的冻融损伤变量、受荷损伤变量和总损伤变量如式（8）～（10）.
Dn=AnA， (8) D=A2A−An， (9) Dm=An+A2A. (10) 联立式（8）～（10）可得
Dm=D+Dn−DDn. (11) 联立式（1）和式（6）、（7）、（11）可得
σi=σi1A1+σr(An+A2)A=σi1(1−Dm)+σrDm. (12) 假定岩石只在主应力方向发生损伤，且未损伤部分服从广义胡克定律，结合岩石变形协调关系得到冻融作用下承受三轴压缩荷载岩石在主应力方向的损伤模型为
σ1=(Enε1+2μnσ3)（1−Dm）+σrDm， (13) 式中：$ {E}_{n} $为岩石冻融循环n次后的弹性模量；$ {\mu }_{n} $为岩石冻融循环n次后的泊松比；$ {\sigma }_{3} $为围压；${\sigma }_{1} $为主应力方向应力；ε1为主应力方向应变.
由宏观损伤力学可知，岩石的冻融损伤变量为
Dn=1−EnE0， (14) 式中：${E}_{0}$为未受冻融作用岩石的弹性模量.
将式（5）、（14）代入式（11）、（13）得
σ1=EnE0(Enε1+2μnσ3−σr)exp(−(ˉff0)m)+σr. (15) 假设岩石微元体强度服从最大拉应变准则，如式（16）.
ˉf=f(ε)=ε1. (16) 将式（16）代入式（15）得
σ1=EnE0(Enε1+2μnσ3−σr)exp(−(ε1f0)m)+σr. (17) 1.2 模型参数确定
参数$ {f}_{0} $、$ m $可用峰值应力法确定[11]，峰值点处应力、应变为
σ1c=σ1, (18) ε1c=ε1, (19) 式中：${\sigma }_{1} $为轴向应力.
峰值点处应力对应变的偏微分等于0，则
∂σ1∂ε1|(σ1c，ε1c)=0. (20) 将式（18）、（19）代入式（17）得
σ1c=EnE0(Enε1c+2μnσ3−σr)exp(−(ε1cf0)m)+σr. (21) 将式（21）代入式（20）得
E2nE0exp(−(ε1cf0)m)[1−m(1+2μnσ3−σrEnε1c)(ε1cf0)m=0. (22) 由式（22）可得
m(1+2μnσ3−σrEnε1c)(ε1cf0)m=1， (23) 则
(ε1cf0)m=Enε1cm(Enε1c+2μnσ3−σr). (24) 将式（24）代入式（21）得
σ1c=EnE0(Enε1c+2μnσ3−σr)×exp(−Enε1cm(Enε1c+2μnσ3−σr))+σr. (25) 由式（25）及Weibull分布函数参数定义可得
m=|Enε1c(Enε1c+2μnσ3−σr)lnEn(Enε1c+2μnσ3−σr)E0(σ1c−σr)|. (26) 由式（24）得
f0=ε1c/[Enε1cm(Enε1c+2μnσ3−σr)]1m. (27) 最后，引入修正系数k对式（17）进行修正得
σ1=EnE0（Enε1+2μnσ3−kσr）exp(−(ε1f0)m)+kσr， (28) $ k $的取值为
k={ε1ε1c,0⩽ε1<ε1c，1,ε1⩾ε1c. (29) 式（28）及式（29）即为本文所建立的基于最大拉应变破坏准则并考虑冻融和荷载耦合作用的寒区冻融岩石损伤模型.
2. 模型验证
2.1 实例一：低孔隙率硬质花岗岩
Tan等[17]对采自西藏嘎隆拉山高速公路隧道中的花岗岩开展了饱和状态下，温度变化范围为−40～+ 40 ℃，循环次数高达150次的冻融循环试验. 文献[17]所采用的花岗岩孔隙率为0.67%，平均单轴抗压强度为135.73 MPa，属于致密的硬质岩石. 为验证本文所建立损伤本构模型的适用范围及合理性，利用文献中围压为10 MPa，冻融循环次数为0、50、100、150次的试验曲线与本文损伤模型计算的理论曲线进行对比分析. 同时，基于Huang等[11]建立的损伤本构模型计算了低孔隙率硬质花岗岩的应力-应变曲线，将其作为对比曲线之一.
表1给出了不同冻融次数下低孔隙率硬质花岗岩的力学参数. 表1中，N为冻融次数. 由于冻融作用对硬质岩石的泊松比影响程度较低，因此，本文在计算过程中，假设不同冻融次数下低孔隙率硬质花岗岩的泊松比数值不变，均为0.15. 试验与理论模型的对比结果如图1所示，图中：理论模型1代表本文所建立的损伤本构模型；理论模型2代表文献[11]中的损伤本构模型.
表 1 低孔隙率硬质花岗岩力学参数Table 1. Mechanical parameters of hard granite with low porosityN/次 围压/
MPa峰值应力/
MPa峰值点
应变残余强度/
MPa弹性模量/
GPa0 10 185.14 0.0052 40 36.71 50 160.02 0.0068 140 26.67 100 123.29 0.0058 90 22.99 150 106.24 0.0064 40 21.12 
图 1 低孔隙率硬质花岗岩试验与理论曲线对比（${\sigma }_{3} $ = 10 MPa）Figure 1. Comparison of experimental and theoretical curves of hard granite with low-porosity (${\sigma }_{3} $ = 10 MPa)从图1中可以看出：当花岗岩冻经历0、50、100次冻融循环作用后，试验数据曲线与理论模型1、理论模型2在峰值强度前均具有较好的对应关系；不同的是，在峰后阶段，理论模型1可较好反映出岩石的残余强度，而理论模型2的应力在峰值强度后跌落速率较快，无法体现岩石的残余强度特征；当冻融循环次数达到150次时，理论模型2在峰值强度前与试验数据吻合度较理论模型1高，此时，虽然理论模型1的峰值强度与试验数据曲线之间存在一定的差距，但曲线整体趋势较为一致. 可见，本文所建立的冻融岩石损伤本构模型可以适用于低孔隙率硬质岩石，且能较好地体现岩石的残余强度特征.
2.2 实例二：高孔隙率软质红砂岩
张慧梅等[16]对未受冻融循环作用及受冻融循环后的红砂岩开展了三轴力学试验，分析了不同冻融循环次数和围压下红砂岩的力学特性. 文献[16]所采用的红砂岩孔隙度为14.26%，单轴抗压强度约为4 MPa，属于高孔隙率软岩[18-19]. 为验证本文所建立的损伤本构模型在高孔隙率软岩中的适用性，利用文献[16]中的试验数据曲线与本文计算的理论曲线进行对比分析. 类似的，此处也加入了与文献[11]中的损伤本构模型的对比. 表2给出了不同冻融次数和围压下高孔隙率软质红砂岩的力学参数. 试验与理论模型的对比结果如图2所示.
表 2 高孔隙率软质红砂岩力学参数Table 2. Mechanical parameters of soft red sandstone with high porosityN/次 围压/MPa 峰值应力/MPa 峰值点应变 残余强度/MPa 弹性模量/GPa 泊松比 0 6 24.866 0.016 22.408 1.649 0.254 5 4 19.132 0.013 14.822 1.452 0.257 10 2 12.701 0.011 7.020 1.156 0.262 由图2可知：对于经历不同冻融作用的高孔隙率软质红砂岩，在峰值应力前，试验数据与理论模型1、理论模型2的吻合性均较好；在峰值强度后阶段，理论模型2应力快速跌落，与试验数据中岩石具有一定的残余强度形变规律不符；理论模型1则可较好地反映高孔隙率软质红砂岩的峰后形变特征. 以上现象进一步证明了理论模型1的合理性，也证明了本文所建立的损伤本构可适用于孔隙率较高的软岩.
由图2、3可知：本文所建立的冻融岩石损伤本构模型具有较广泛的适用性，不仅适用于低孔隙率硬岩，也适用于高孔隙率软岩；同时，本文所建立的模型还可较好地反映出不同类别岩石的峰后形变特征.
图 2 高孔隙率软质红砂岩试验与理论曲线对比Figure 2. Comparison of experimental and theoretical curves of soft red sandstone with high porosity3. 损伤变量演化规律
由式（5）、（11）、（14）及文献[16-17]中的试验数据计算得到两类岩石在冻融循环下的总损伤演变曲线，如图3和图4所示.
图 3 低孔隙率硬质花岗岩总损伤演变曲线Figure 3. Total damage evolution curves of hard granite with low-porosity
图 4 高孔隙率软质红砂岩总损伤演变曲线Figure 4. Total damage evolution curves of soft red sandstone with high porosity从图3中可以看出：不同冻融循环次数下花岗岩的总损伤演变曲线均呈“S”形；未受冻融影响的花岗岩初始阶段总损伤近似为0，受冻融循环作用的花岗岩在初始阶段总损伤为冻融损伤量，且初始阶段的总损伤量随冻融次数的增加而增加；随着外力的加载，岩石内部微裂纹和缺陷不断发展扩张，受荷损伤开始快速增加，花岗岩的总损伤曲线也呈现出加速损伤的趋势，岩石逐渐进入屈服阶段；当岩石强度逐渐丧失后，总损伤曲线斜率开始变缓，并最终趋于稳定值1，此时岩石已完全破坏，总损伤达到最大值. 由图4可知：红砂岩的总损伤曲线也呈现出近似“S”形的趋势，其演变规律与花岗岩类似，不再展开描述.
基于以上分析，将低孔隙率硬质花岗岩和高孔隙率软质红砂岩总损伤演变曲线分为初始损伤阶段、加速阶段和完全损伤阶段（图5）：在初始损伤阶段，岩石内的总损伤量较低，且随着应变的增加，总损伤量近乎不变；在加速损伤阶段，总损伤量随应变的增加而急剧增长，岩石内产生大量的裂纹；在完全损伤阶段，总损伤的增长趋势逐渐变缓并最终达到稳定，岩石在此阶段发生宏观破坏.
图 5 冻融受荷岩石总损伤演变规律Figure 5. The evolution of total damage of rock under loading and freeze-thaw action4. 模型参数分析
4.1 参数m和$ {\mathit{f}}_{0} $对应力-应变曲线的影响
为确定本文所建立的损伤模型中参数$ m $和$ {f}_{0} $的物理意义及对应力-应变曲线的影响，以冻融循环50次，围压为10 MPa的低孔隙率硬质花岗岩以及冻融循环5次，围压为4 MPa的高孔隙率软质红砂岩为例进行参数分析，如图6、7所示.
图 7 参数$ {f}_{0} $对应力-应变曲线的影响Figure 7. Influence of parameter $ {f}_{0} $ on the stress-strain curve从图6中可以看出：随着$ m $的增加，低孔隙率硬质花岗岩和高孔隙率软质红砂岩的峰值强度均有所增加；在峰值应力后，应力跌落速率随$ m $的增大而增大，说明岩石的脆性特征逐渐增强. 从图7中可以看出：随着$ {f}_{0} $的增加，两类岩石的应力-应变曲线峰值点升高且向右移动，达到峰值应力所需的应变量也随之增加；岩石的塑性变形特征随着$ {f}_{0} $的增加而增加. 由此可以推断出，参数$ m $和$ {f}_{0} $的物理意义分别为岩石的脆性变形特征和塑性变形特征.
4.2 参数$ m $和$ {f}_{0} $对总损伤曲线的影响
由式（5）、（11）和式（14）计算得到总损伤量，并分析参数$ m $和$ {f}_{0} $对总损伤量的影响，如图8、9所示. 由图8可知：随着参数$ m $的增加，2类岩石的初始损伤阶段皆有所延长，即，岩石进入加速损伤阶段的临界应变值变大；加速损伤阶段的斜率随$ m $的增加而增加，且岩石进入完全损伤阶段的临界应变值随$ m $的增加而减小，说明了岩石提前进入完全损伤阶段. 由图9可知：随着参数$ {f}_{0} $的增加，2类岩石的初始损伤阶段皆有所延长，这与参数$ m $对总损伤曲线的影响规律一致；不同的是，随着$ {f}_{0} $的增加，2类岩石加速损伤阶段的斜率均逐渐减小，进入完全损伤阶段的临界应变值增大，即岩石将逐渐滞后进入完全损伤阶段. 以上现象进一步表明参数$ m $可反映岩石的脆性形变特征，而参数$ {f}_{0} $则代表岩石的塑性形变特征.
图 9 参数$ {f}_{0} $对总损伤量的影响Figure 9. Influence of parameter $ {f}_{0} $ on the total damage5. 工程算例分析
为分析冻融循环对大型工程稳定性的影响，本节将以隧道工程为算例进行数值模拟分析. 图10给出了工程算例的计算模型. 如图所示，该隧道断面为三心拱形，拱顶埋深24.5 m，隧道围岩为软质红砂岩. 该隧道采用全断面开挖法，开挖循环进尺为2 m，计算开挖总长度为20 m. 隧道支护方式为锚喷支护，喷射混凝土等级为C20，锚喷衬砌厚度0.15 m，锚杆截面面积为490.63 mm2，长度为4 m，锚杆间弧长1.8 m.
图 10 隧道工程算例计算模型（单位：m）Figure 10. Calculation model of tunnel engineering example （unit: m）基于文献[16]中所给出不同冻融次数下的红砂岩力学参数，在大型岩土工程有限元软件MIDAS GTS中分别计算了冻融循环次数为0、5、10次时，隧道开挖后的位移场和应力场分布. 不同工况下隧道围岩的力学参数见表3，在计算过程中假设围岩的容重和内摩擦不随冻融次数而改变. 计算得到的围岩竖向位移场、围岩最大主应力矢量场以及衬砌结构最大主应力场分别如图11～13所示.
表 3 不同冻融次数下隧道围岩力学参数Table 3. Mechanical parameters of surrounding rock for tunnel under different freeze-thaw cyclesN/次 弹性模量/GPa 黏聚力/MPa 泊松比 容重/（kN·m−3） 内摩擦角/（°） 0 1.387 1.91 0.258 22 30 5 1.295 1.62 0.259 10 1.156 1.41 0.262 
图 12 不同冻融次数下隧道围岩最大主应力矢量场Figure 12. Maximum principal stress vector field of surrounding rock for tunnel under different freeze-thaw cycles从图11中可以看出：隧道开挖后拱顶和底部的围岩垂直位移量较大；当围岩经历0、5、10次冻融循环后，隧道拱顶的最大下沉量分别为3.75 、3.99 、4.42 mm，经历10次冻融循环后，隧道拱顶的围岩下沉量增加了17.87%；类似地，当围岩经历0、5、10次冻融循环后，隧道底部的最大隆起量分别为4.73 、5.05 、5.64 mm，经历10次冻融循环后，隧道底部的围岩隆起量增加了19.24%.
从图12中可以看出：隧道开挖后拱顶与底部分布拉应力积聚区，隧道两侧围岩分布压应力积聚区；当围岩经历0、5、10次冻融循环后，隧道围岩的最大拉应力极值分别为0.297 、0.300、0.305 MPa，经历10次冻融循环后，隧道围岩最大拉应力极值增加了2.7%；当围岩经历0、5、10次冻融循环后，隧道围岩的最大压应力极值分别为0.448、0.450、0.457 MPa，经历10次冻融循环后，隧道围岩的最大压应力极值增加了2.01%. 隧道围岩的最大主应力随冻融循环变化较小，与文献[11]中得出的结论一致，这是由于围岩内主要是自重应力，自重应力在冻融循环下基本保持不变. 因此，围岩的最大主应力变化量很小，而由于岩石力学性质的改变，围岩位移随冻融循环的变化量则较大.
此外，隧道衬砌结构的应力分布也受到冻融作用的影响，从图13中可以看出：隧道衬砌拱顶和底部主要为拉应力积聚区，衬砌两侧分布压应力积聚区，随着冻融循环次数的增加，隧道衬砌拉、压应力极值也逐渐增大；当围岩经历0、5、10次冻融循环后，隧道衬砌最大拉应力极值分别为4.88、5.25、5.93 MPa，经历10次冻融循环后，隧衬砌最大拉应力极值增加了21.52%；当围岩经历0、5、10次冻融循环后，隧道衬砌最大压应力极值分别为3.75、3.99、4.42 MPa，经历10次冻融循环后，隧衬砌最大压应力极值增加了17.87%. 冻融作用后隧道围岩的力学性质劣化，导致围岩压力传递至衬砌结构，引起衬砌应力的增加. 由此可知，虽然冻融作用对隧道围岩应力场影响不大，但随着围岩的软化也可引起隧道结构的破坏，影响隧道的整体稳定性. 因此，在寒区进行岩土工程建设时，应特别注意对冻融循环对岩土性质的影响并采取适当的防治措施.
图 11 不同冻融次数下隧道围岩竖向位移Figure 11. Vertical displacement of surrounding rock for tunnel under different freeze-thaw cycles
图 13 不同冻融次数下隧道衬砌最大主应力Figure 13. Maximum principal stress in tunnel lining under different freeze-thaw cycles6. 结 论
本文基于统计损伤力学理论，假设岩石微元体强度服从Weibull分布，结合最大拉应变破坏准则建立了考虑冻融和荷载耦合作用的寒区冻融岩石损伤本构模型；基于试验数据及已有的损伤本构模型验证了所建立的损伤本构模型的正确性，并对总损伤演变过程及模型参数进行了分析；最后，在数值模拟软件中分析了冻融循环对隧道工程稳定性的影响，得到的主要结论如下：
1） 以低孔隙率硬质花岗岩和高孔隙率软质红砂岩为验证对象，结合试验数据及已有的损伤本构模型确定所建立模型的合理性和适用范围，结果表明，本文所建立的模型可较好地重现2类岩石的应力-应变全过程曲线，且能较好地反映岩石的峰后强度.
2） 不同冻融循环次数下低孔隙率硬质花岗岩和高孔隙率软质红砂岩的总损伤演变曲线均呈“S”形，且总损伤演变曲线可分为初始损伤阶段、加速阶段和完全损伤阶段.
3） Weibull分布参数$ m $和$ {f}_{0} $分别代表岩石的脆性特征和塑性特征；随着$ m $增加，2类岩石的应力峰值逐渐增加，峰后应力跌落速率逐渐增大，总损伤曲线加速损伤阶段的斜率也逐渐增加，岩石提前进入完全损伤阶段；随着$ {f}_{0} $的增加，岩石达到峰值应力所对应的应变逐渐增加，总损伤曲线加速损伤阶段斜率逐渐减小，岩石滞后进入完全损伤阶段.
4） 采用数值模拟方法，计算了冻融循环对隧道围岩位移、应力及隧道衬砌应力的影响，随着冻融次数的增加，隧道围岩的垂直位移、最大主应力及衬砌最大主应力均逐渐增加；冻融循环对隧道围岩的应力分布影响较小，但由于围岩物理力学性质的劣化，围岩压力传递至衬砌可导致隧道结构产生损伤.
-
图 1 低孔隙率硬质花岗岩试验与理论曲线对比（${\sigma }_{3} $ = 10 MPa）
Figure 1. Comparison of experimental and theoretical curves of hard granite with low-porosity (${\sigma }_{3} $ = 10 MPa)
图 2 高孔隙率软质红砂岩试验与理论曲线对比
Figure 2. Comparison of experimental and theoretical curves of soft red sandstone with high porosity
图 3 低孔隙率硬质花岗岩总损伤演变曲线
Figure 3. Total damage evolution curves of hard granite with low-porosity
图 4 高孔隙率软质红砂岩总损伤演变曲线
Figure 4. Total damage evolution curves of soft red sandstone with high porosity
图 5 冻融受荷岩石总损伤演变规律
Figure 5. The evolution of total damage of rock under loading and freeze-thaw action
图 7 参数$ {f}_{0} $对应力-应变曲线的影响
Figure 7. Influence of parameter $ {f}_{0} $ on the stress-strain curve
图 9 参数$ {f}_{0} $对总损伤量的影响
Figure 9. Influence of parameter $ {f}_{0} $ on the total damage
图 10 隧道工程算例计算模型（单位：m）
Figure 10. Calculation model of tunnel engineering example （unit: m）
图 12 不同冻融次数下隧道围岩最大主应力矢量场
Figure 12. Maximum principal stress vector field of surrounding rock for tunnel under different freeze-thaw cycles
图 11 不同冻融次数下隧道围岩竖向位移
Figure 11. Vertical displacement of surrounding rock for tunnel under different freeze-thaw cycles
图 13 不同冻融次数下隧道衬砌最大主应力
Figure 13. Maximum principal stress in tunnel lining under different freeze-thaw cycles
表 1 低孔隙率硬质花岗岩力学参数
Table 1. Mechanical parameters of hard granite with low porosity
N/次 围压/
MPa峰值应力/
MPa峰值点
应变残余强度/
MPa弹性模量/
GPa0 10 185.14 0.0052 40 36.71 50 160.02 0.0068 140 26.67 100 123.29 0.0058 90 22.99 150 106.24 0.0064 40 21.12 
下载: 导出CSV
表 2 高孔隙率软质红砂岩力学参数
Table 2. Mechanical parameters of soft red sandstone with high porosity
N/次 围压/MPa 峰值应力/MPa 峰值点应变 残余强度/MPa 弹性模量/GPa 泊松比 0 6 24.866 0.016 22.408 1.649 0.254 5 4 19.132 0.013 14.822 1.452 0.257 10 2 12.701 0.011 7.020 1.156 0.262
下载: 导出CSV
表 3 不同冻融次数下隧道围岩力学参数
Table 3. Mechanical parameters of surrounding rock for tunnel under different freeze-thaw cycles
N/次 弹性模量/GPa 黏聚力/MPa 泊松比 容重/（kN·m−3） 内摩擦角/（°） 0 1.387 1.91 0.258 22 30 5 1.295 1.62 0.259 10 1.156 1.41 0.262
下载: 导出CSV
-
[1] 张德,刘恩龙,刘星炎,等. 基于修正Mohr-Coulomb屈服准则的冻结砂土损伤本构模型[J]. 岩石力学与工程学报,2018,37(4): 978-986. doi: 10.13722/j.cnki.jrme.2017.1318ZHANG De, LIU Enlong, LIU Xingyan, et al. A damage constitutive model for frozen sandy soils based on modified Mohr-Coulomb yield criterion[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(4): 978-986. doi: 10.13722/j.cnki.jrme.2017.1318 [2] 李玉成,陈有亮,孙浩程,等. 寒区花岗岩冻融损伤破坏的试验研究[J]. 工业建筑,2019,49(7): 83-88,107.LI Yucheng, CHEN Youliang, SUN Haocheng, et al. An experimental study of granite’s damage caused by freeze-thaw cycles in cold region[J]. Industrial Construction, 2019, 49(7): 83-88,107. [3] 罗学东,黄成林,肜增湘,等. 冻融循环作用下蒙库铁矿边坡岩体物理力学特性研究[J]. 岩土力学,2011,32(增1): 155-159.LUO Xuedong, HUANG Chenglin, RONG Zengxiang, et al. Study of physico-mechanical characteristics of rocks in slope of Mengku iron mine under freezing-thawing cyclic effect[J]. Rock and Soil Mechanics, 2011, 32(S1): 155-159. [4] 吕文韬,杨龙,魏云杰,等. 新疆塔县地区片麻岩冻融劣化机理与规律试验研究[J]. 水文地质工程地质,2019,46(3): 95-100.LYU Wentao, YANG Long, WEI Yunjie, et al. Research on mechanism of freezing-thawing deterioration of gneisses in the Taxian area of Xinjiang[J]. Hydrogeology & Engineering Geology, 2019, 46(3): 95-100. [5] 方云,乔梁,陈星,等. 云冈石窟砂岩循环冻融试验研究[J]. 岩土力学,2014,35(9): 2433-2442.FANG Yun, QIAO Liang, CHEN Xing, et al. Experimental study of freezing-thawing cycles on sandstone in Yungang grottos[J]. Rock and Soil Mechanics, 2014, 35(9): 2433-2442. [6] 徐拴海,李宁,王晓东,等. 露天煤矿冻岩边坡饱和砂岩冻融损伤试验与劣化模型研究[J]. 岩石力学与工程学报,2016,35(12): 2561-2571.XU Shuanhai, LI Ning, WANG Xiaodong, et al. Damage test and degradation model of saturated sandstone due to cyclic freezing and thawing of rock slopes of open-pit coal mine[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(12): 2561-2571. [7] HUANG S B, LIU Q S, LIU Y Z, et al. Frost heaving and frost cracking of elliptical cavities (fractures) in low-permeability rock[J]. Engineering Geology, 2018, 234: 1-10. doi: 10.1016/j.enggeo.2017.12.024 [8] PARK J, HYUN C U, PARK H D. Changes in microstructure and physical properties of rocks caused by artificial freeze-thaw action[J]. Bulletin of Engineering Geology and the Environment, 2015, 74(2): 555-565. doi: 10.1007/s10064-014-0630-8 [9] 李新平,路亚妮,王仰君. 冻融荷载耦合作用下单裂隙岩体损伤模型研究[J]. 岩石力学与工程学报,2013,32(11): 2307-2315.LI Xinping, LU Yani, WENG Yangjun. Research on damage model of single jointed rock masses under coupling action of freeze-thaw and loading[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(11): 2307-2315. [10] 张慧梅,雷利娜,杨更社. 温度与荷载作用下岩石损伤模型[J]. 岩石力学与工程学报,2014,33(增2): 3391-3396.ZHANG Huimei, LEI Lina, YANG Gengshe. Damage model of rock under temperature and load[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(S2): 3391-3396. [11] HUANG S B, LIU Q S, CHENG A P, et al. A statistical damage constitutive model under freeze-thaw and loading for rock and its engineering application[J]. Cold Regions Science and Technology, 2018, 145: 142-150. doi: 10.1016/j.coldregions.2017.10.015 [12] JIA H L, XIANG W, KRAUTBLATTER M. Quantifying rock fatigue and decreasing compressive and tensile strength after repeated freeze-thaw cycles[J]. Permafrost and Periglacial Processes, 2015, 26(4): 368-377. doi: 10.1002/ppp.1857 [13] 袁超,张慧梅,孟祥振,等. 内摩擦角对冻融岩石损伤本构模型的影响探讨[J]. 煤田地质与勘探,2019,47(4): 138-143. doi: 10.3969/j.issn.1001-1986.2019.04.021YUAN Chao, ZHANG Huimei, MENG Xiangzhen, et al. Influence of internal friction angle on damage constitutive model of freeze-thaw rock[J]. Coal Geology & Exploration, 2019, 47(4): 138-143. doi: 10.3969/j.issn.1001-1986.2019.04.021 [14] 张峰瑞,姜谙男,杨秀荣,等. 冻融循环下花岗岩剪切蠕变试验与模型研究[J]. 岩土力学,2020,41(2): 509-519.ZHANG Fengrui, JIANG Annan, YANG Xiurong, et al. Experimental and model research on shear creep of granite under freeze-thaw cycles[J]. Rock and Soil Mechanics, 2020, 41(2): 509-519. [15] JEAN L. A continuous damage mechanics model for ductile fracture[J]. Journal of Engineering Materials and Technology, 1985, 107(1): 83-89. doi: 10.1115/1.3225775 [16] 张慧梅,孟祥振,彭川,等. 岩石变形全过程冻融损伤模型及其参数[J]. 西安科技大学学报,2018,38(2): 260-265.ZHANG Huimei, MENG Xiangzhen, PENG Chuan, et al. Freeze-thaw damage model and parameters of rock deformation in whole process[J]. Journal of Xi’an University of Science and Technology, 2018, 38(2): 260-265. [17] TAN X J, CHEN W Z, YANG J P, et al. Laboratory investigations on the mechanical properties degradation of granite under freeze-thaw cycles[J]. Cold Regions Science and Technology, 2011, 68(3): 130-138. doi: 10.1016/j.coldregions.2011.05.007 [18] 张慧梅,杨更社. 冻融岩石损伤劣化及力学特性试验研究[J]. 煤炭学报,2013,38(10): 1756-1762.ZHANG Huimei, YANG Gengshe. Experimental study of damage deterioration and mechanical properties for freezing-thawing rock[J]. Journal of China Coal Society, 2013, 38(10): 1756-1762. [19] 张慧梅,孟祥振,彭川,等. 冻融-荷载作用下基于残余强度特征的岩石损伤模型[J]. 煤炭学报,2019,44(11): 3404-3411.ZHANG Huimei, MENG Xiangzhen, PENG Chuan, et al. Rock damage constitutive model based on residual intensity characteristics under freeze-thaw and load[J]. Journal of China Coal Society, 2019, 44(11): 3404-3411. -
下载:
下载:
点击查看大图
计量
- 文章访问数: 495
- HTML全文浏览量: 170
- PDF下载量: 42
- 被引次数: 0