| 摘要: |
| 自然界中岩体破坏的主要形式是沿着滑面的剪切破坏，且主要表现为剪切流变损伤特征。根据软弱夹层不同剪切应力水平下的试验数据分析，引入可以表征其流变过程中参数损伤的变量D，提出了一个基于D的可以反映软弱夹层加速流变特性的黏弹塑性非线性流变模型，与伯格斯模型串联起来，构成了一个能全面反映三个流变阶段的新的软弱夹层剪切流变损伤模型，进而写出了不同应力水平下该模型的流变方程，并推导出了模型三维情况下的有限差分格式，基于FLAC3D使用C++语言对该流变损伤模型进行了二次开发，与室内流变试验进行了对比验证。结果表明：（1）每级剪应力下应变的增量值，其瞬时应变都较为接近，且经过相同时间的应力加载后，无论是稳定后的应变值还是加速阶段的应变值都较为接近。（2）数值模拟试验在施加的剪切荷载未到达屈服应力，其减速阶段历时较室内试验减速阶段历时短，即更快地达到稳定流变阶段，而在施加的剪切荷载超过屈服应力阈值时，其加速流变阶段位移变化速率较室内试验位移变化速率大，且最终的位移比室内试验位移稍大。说明了开发的流变损伤模型是可靠的，研究结果可为其他本构模型的二次开发提供一定参考。 |
| 关键词: 剪切流变 流变损伤 非线性流变 数值模拟 FLAC3D二次开发 |
| DOI： |
| 分类号:TU485 |
| 基金项目:国家自然科学基金资助项目(41672317,U1802243)；湖北省技术创新专项(重大项目)(2017ACA184)；武科大重大科技项目培育类创新团队(A类)(2018TDX01)。 |
|
| Study on shear rheology constitutive model of weak interlayer considering damage and its program realization |
|
ZHU Xin, HU Bin, LI Jing, CUI Kai, WEI Erjian, LIU Yang
|
|
School of Resources and Environmental Engineering,Wuhan University of Science and Technology,Wuhan Hubei
|
| Abstract: |
| The destruction of the rock mass is the shear damage that mainly 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, and 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 equations of the rheological damage model under different stress levels were written, and the finite difference scheme for the three-dimensional case of the model was derived, and the rheological damage model was redeveloped with C + + language based on FLAC3D, and compared with the indoor rheological test. The results showed that: (1) the instantaneous strain of the increment value of strain under each shear stress was 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 did not reach the yield stress, the deceleration period of the numerical simulation test was shorter than that of the indoor test, that was, it reached the stable rheological stage faster. When the applied shear load exceeded the yield stress threshold, the displacement change rate of the accelerated rheological stage was larger than that of the indoor test, and the final displacement ws slightly larger than the indoor test displacement. It showed that the developed rheological damage model was reliable, and the research results can provide some reference for the secondary development of other constitutive models. |
| Key words: shear rheology rheology damage nonlinear rheology numerical simulation secondary development of FLAC3D |
