PDF(755 KB)
PDF(755 KB)
PDF(755 KB)
圆形刚性板岩体蠕变试验五参量蠕变模型推导及应用
Derivation and Application of a Generalized Five-Parameter Creep Model for Creep Test of Circular Rigid Plate Rock Mass
岩体蠕变时间一般长达几年,而岩块蠕变达到稳定的时间仅需几小时,两者相差显著。现场岩体蠕变试验蠕变持续时间长,能较好地反映岩体的蠕变时效性。详细推导了圆形刚性承压板岩体蠕变试验广义Kelvin五参量模型变形表达式,该表达式在特定条件下可以退化为广义Kelvin三参量模型或Burgers模型,对岩体蠕变试验的参数辨识具有更好地适应性。对钙质板岩的现场蠕变试验按五参量模型进行了参数辨识,与三参量模型拟合参数对比,五参量模型拟合精度更高,得到的长期剪切模量更小,且能间接反映钙质板岩的加载特征,因此五参量模型能更好地反映岩体的蠕变特性。
[Objective] The creep time of rock mass generally lasts for several years, while the time for rock block to reach creep stability only takes a few hours. There is a significant difference in the creep timeliness between rock mass and rock block. On-site rock mass creep tests generally last not less than 10 days, and in some cases can reach 1 to 2 months, which can better reflect the creep timeliness of the rock mass. Therefore, on-site rock mass creep tests more accurately represent the creep timeliness of rock mass. This study aims to simulate the creep test process of on-site circular rigid plate rock mass creep test through the generalized five-parameter Kelvin creep model and to obtain the deformation expression of the on-site rigid plate creep test:$w(t)=\frac{\pi p R}{8}\left\{\frac{1}{G_{\mathrm{H}}}+\frac{1-\mathrm{e}^{\left(-G_{1} t / \eta_{1}\right)}}{G_{1}}+\frac{1-\mathrm{e}^{\left(-G_{2} t / \eta_{2}\right)}}{G_{2}}+\frac{3}{B}+\frac{3 G_{\mathrm{H}}^{2}}{B b_{2}\left(a_{1}-a_{2}\right)}\right.$$\left.\left[\frac{G_{1}+G_{2}}{a_{1} a_{2}}\left(a_{1}-a_{2}+a_{2} \mathrm{e}^{\left(-a_{1} t\right)}-a_{1} \mathrm{e}^{\left(-a_{2} t\right)}\right)-\left(\eta_{1}+\eta_{2}\right)\left(\mathrm{e}^{\left(-a_{1} t\right)}-\mathrm{e}^{\left(-a_{2} t\right)}\right)\right]\right\}$. [Methods] The generalized five-parameter Kelvin creep model consists of one spring element in series with two Kelvin elements. According to the deformation expression of the generalized five-parameter Kelvin creep model, when G1 or G2 equals 0, the generalized five-parameter Kelvin creep model degenerates into the Burgers model; when G1 or G2 approaches infinity, the generalized five-parameter Kelvin creep model degenerates into the generalized three-parameter Kelvin creep model. The generalized five-parameter Kelvin creep model is a stable model, with a long-term shear modulus of G∞=1/(1/GH+1/G1+1/G2).Therefore, the generalized five-parameter Kelvin creep model has higher fitting accuracy and wider applicability compared to the generalized three-parameter Kelvin creep model. For the generalized Kelvin five-parameter creep model, since the four variables of G∞,GH,G1,G2 and E∞,EH,E1,E2 don’t correspond one-to-one, it means that in three-dimensional rock engineering problems, the viscosity coefficient derived from the shear modulus isn’t the same as the one derived from the elastic modulus. [Results] On-site rigid plate creep tests were conducted on calcareous slate, and parameter fitting was carried out using both the generalized five-parameter Kelvin creep model and the generalized three-parameter Kelvin creep model. Key indicators during the fitting of the creep expression, such as Reduced Chi-Sqr, Adj. R2, and RMSE, show that the generalized five-parameter Kelvin creep model outperforms the generalized three-parameter Kelvin creep model. Under different stress conditions, the shear modulus G1 of the generalized five-parameter Kelvin creep model increases with increasing stress, while G2 first decreases and then increases with increasing stress. The parameter G1 of the generalized three-parameter Kelvin creep model also slightly increases with increasing stress, but the change is minor. The parameters of the generalized five-parameter Kelvin creep model indirectly reflect the variations during the rock mass fracture compaction stage and fracture closure stage. The long-term shear modulus of the generalized five-parameter Kelvin creep model is lower than that of the generalized three-parameter creep model. Under low stress, the long-term deformation of the generalized five-parameter Kelvin creep model increases by about 50.8% from that of the generalized three-parameter creep model; under high stress, it only increases by about 2.71%. Therefore, the generalized five-parameter Kelvin creep model is capable of expressing greater deformation, which has benefit for engineering safety. [Conclusions] The generalized five-parameter Kelvin creep model has higher fitting accuracy and can be compatible with both the generalized three-parameter Kelvin creep model and the Burgers creep model. According to the generalized five-parameter Kelvin creep model, the long-term shear modulus is lower, and long-term deformation is greater, providing a conservative evaluation of long-term rock mass deformation, which can reduce potential safety risks. In summary, the generalized five-parameter Kelvin creep model has good fitting accuracy and can provide rich information, allowing it better at reflecting the creep characteristics of rock mass.
| [1] |
|
| [2] |
Failure to properly characterize time-dependent behavior (i.e. creep) of rocks and the use of inappropriate behavior models when analyzing stability of underground spaces such as mining drifts and stopes, road tunnels, different types of caverns and reservoirs for storing natural gas, petroleum fluids, and compressed air energy as well as nuclear waste disposal caverns will end up with costly losses. The most important factors affecting creep behavior of rock masses include composition of the minerals composing the rock, size of the grains or crystals of which the rock is composed, humidity, temperature, time, loading scheme, loading rate, strain rate and loading frequency. On this basis, in order to attain a proper understanding of creep behavior of rocks, it is necessary to determine creep coefficients of rocks via in-lab experiments. In the present research, rock salt samples were prepared in the form of cylinders at length-to-diameter ratios of greater than 2. Performed experiments in this study were long-term creep (LTC) tests at 6 stress levels of 5.5, 7.5, 10, 12, 14.02, and 18MPa and stepwise short-term creep (STC) tests (each at 3 stress levels, namely 4.4, 10.1, and 11.9MPa, and 7.5, 12, and 17MPa, respectively). The tests were performed at ambient conditions at a constant temperature of 22 degrees C and relative humidity of 23%. Then, according to Burger's rheological model and the Lubby 2 constitutive model and analyzing the obtained information from the tests, Kelvin's viscoelastic coefficient and shear modulus as well as Maxwell's elastic and viscoplastic coefficients were calculated. Also, obtained creep coefficients were improved using linear and nonlinear regression analysis of experiment dates.
|
| [3] |
项明, 张喆, 李梦博, 等. 盐岩岩芯制备方法及蠕变特性[J]. 科学技术与工程, 2025, 25(3): 1018-1027.
(
|
| [4] |
郭燕忠. 基于参数修正模型的JM矿岩石流变特性模拟[J]. 矿业研究与开发, 2024, 44(7):182-190.
(
|
| [5] |
|
| [6] |
Creep behaviour in rocks is a typical mechanical property that is directly linked to the stability of underground engineering. The deformations and rate of rocks creep are not only influenced by time but also by the loading and unloading history. To more accurately predict creep mechanical behaviour of salt rocks, the rocks hardening is described by introducing a state variable. A new three-dimensional creep constitutive model of salt rocks was established to describe the loading and unloading history effect of the rheological properties. In this paper, salt rocks creep tests under various loading and unloading histories were conducted to investigate how different loading routes affect the creep behaviour of salt rocks. The effects of the model state variables were analysed through different indicators. An example verification was carried out with the results of plastic deformation tests performed at different loading paths. The findings indicated the creep rate of stepped loading and the stepped unloading under the same stress level were significantly affected by the loading history. The proposed constitutive model can accurately fit the creep test curves of different loading paths, indicate that it can provided a prediction of the historical effect of the creep behaviour of salt rocks. Different parameters affect the different phases of the creep curve. The parameter k primarily affects the overall shape of the creep curves. Parameters m and c primarily influence the steady-state creep length and creep rates, excluding the initial cycle.
|
| [7] |
陈卫忠, 谭贤君, 吕森鹏, 等. 深部软岩大型三轴压缩流变试验及本构模型研究[J]. 岩石力学与工程学报, 2009, 28(9): 1735-1744.
(
|
| [8] |
冯杰, 钟志彬, 吕蕾, 等. 侧向约束下原状软碎岩体变形特性试验研究[J]. 地下空间与工程学报, 2025, 21(3): 950-958.
(
|
| [9] |
夏维学, 张巨会. 原状软碎岩体室内被动三轴变形特性研究[J]. 四川水力发电, 2021, 40(1): 1-6.
(
|
| [10] |
白高峰, 赵丰博, 霍雨佳, 等. 某工程地基岩体流变特性现场试验研究[J]. 地基处理, 2021, 3(5): 388-396.
(
|
| [11] |
张强勇, 张建国, 杨文东, 等. 软弱岩体蠕变模型辨识与参数反演[J]. 水利学报, 2008, 39(1):66-72.
(
|
| [12] |
熊诗湖, 周火明, 钟作武. 岩体载荷蠕变试验方法研究[J]. 岩石力学与工程学报, 2009, 28(10):2121-2127.
(
|
| [13] |
郭向龙. 粘弹性地基计算模型综述[J]. 山西建筑, 2010, 36(32):117-119.
(
|
| [14] |
王芝银, 李云鹏. 岩体流变理论及其数值模拟[M]. 北京: 科学出版社, 2008.
(
|
| [15] |
范广勤. 岩土工程流变力学[M]. 北京: 煤炭工业出版社, 1993.
(
|
| [16] |
于学馥, 郑颖人, 刘怀恒, 等. 地下工程围岩稳定分析[M]. 北京: 煤炭工业出版社, 1983.
(
|
| [17] |
丁秀丽. 岩体流变特性的试验研究及模型参数辨识[D]. 武汉: 中国科学院武汉岩土力学研究所, 2005.
(
|
| [18] |
徐芝纶. 弹性力学[M]. 2版. 北京: 人民教育出版社,1982:228.
(
|
| [19] |
刘允芳, 周密. 承压垫板法试验的研究[J]. 岩石力学与工程学报, 2005, 24(2): 257-262.
(
|
/
| 〈 |
|
〉 |