祝鑫,胡斌,李京,等. 考虑损伤的软弱夹层剪切流变模型及程序实现[J]. 水利水运工程学报,2021(6):124-132. doi: 10.12170/20201228002
引用本文: 祝鑫,胡斌,李京,等. 考虑损伤的软弱夹层剪切流变模型及程序实现[J]. 水利水运工程学报,2021(6):124-132. doi: 10.12170/20201228002
(ZHU Xin, HU Bin, LI Jing, et al. Study on shear rheology constitutive model of weak interlayer considering damage and its program realization[J]. Hydro-Science and Engineering, 2021(6): 124-132. (in Chinese)). doi: 10.12170/20201228002
Citation: (ZHU Xin, HU Bin, LI Jing, et al. Study on shear rheology constitutive model of weak interlayer considering damage and its program realization[J]. Hydro-Science and Engineering, 2021(6): 124-132. (in Chinese)). doi: 10.12170/20201228002

考虑损伤的软弱夹层剪切流变模型及程序实现

Study on shear rheology constitutive model of weak interlayer considering damage and its program realization

  • 摘要: 自然界中岩体破坏的主要形式是沿着滑面的剪切破坏,且主要表现为剪切流变损伤特征。根据软弱夹层不同剪切应力水平下的试验数据分析,引入可以表征其流变过程中参数损伤的变量D,提出了一个基于D的可以反映软弱夹层加速流变特性的黏弹塑性非线性流变模型,与伯格斯模型串联构成了一个能全面反映3个流变阶段的新的软弱夹层剪切流变损伤模型。基于FLAC3D使用C++语言对该流变损伤模型进行了二次开发,与室内流变试验进行了对比验证。结果表明:(1)每级剪应力下,其瞬时应变均较为接近,且经过相同时间的应力加载后,无论是稳定后的应变值还是加速阶段的应变值也均较为接近。(2)数值模拟试验中,当施加的剪切荷载未达到屈服应力时,其减速阶段历时较室内试验的历时短,即更快地达到稳定流变阶段;而当施加的剪切荷载超过屈服应力阈值时,其加速流变阶段位移变化速率大于室内试验的,且最终的位移比室内试验位移稍大。

     

    Abstract: The destruction of the rock mass is the shear damage that mainly occurs along the failure surface in nature, and it is mainly characterized by shear rheological damage. Based on the analysis of experimental data of weak intercalation under different shear stress levels, the variable D which can characterize the parameter damage in the rheological process was introduced. A viscoelastic plastic nonlinear rheological model based on D, which can reflect the accelerated rheological characteristics of weak intercalation, was proposed. And a new shear rheological damage model of weak intercalation was constructed by connecting it with the Burgers model, which can fully reflect the three rheological stages. Then the rheological damage model was redeveloped with C++ language based on FLAC3D, and compared with the indoor rheological test. The results show that: (1) The instantaneous strain of the increment value of strain under each shear stress is quite close, and both the stable strain value and the accelerated strain value are very close after the same time of stress loading. (2) When the applied shear load does not reach the yield stress, the deceleration period of the numerical simulation test is shorter than that of the indoor test, that is, it reaches the stable rheological stage faster. When the applied shear load exceeds the yield stress threshold, the displacement change rate of the accelerated rheological stage is larger than that of the indoor test, and the final displacement is slightly larger than the indoor test displacement.

     

/

返回文章
返回