# Study on the Nonlinear Damage Creep Model of the Weak Interlayer

Authors: Erjian Wei; Bin Hu; Kunyun Tian; Peishan Cen; Zhen Zhang; Zeqi Wang; Shuxiang Chang
Journal: Advances in Civil Engineering
 (2022)

---

## Abstract

The weak interlayer has become a weak link in slope engineering due to its rheological effect. It is of great significance to study the nonlinear creep model of weak interlayer for long-term stability of the slope. In this paper, based on the creep curve characteristics of weak interlayer and considering the influence of aging damage, the nonlinear improvement of a classical viscoplastic body under stress and time-double threshold conditions is carried out, so that it can more accurately reflect the accelerated creep characteristics of the weak interlayer. By analyzing the relationship between failure load and time, the accelerated creep time threshold of the weak interlayer is obtained. On this basis, a nonlinear damage creep constitutive model of the weak interlayer is constructed and its creep equation is derived. By using the self-defined function fitting tool of Origin software and the Levenberg–Marquardt optimization algorithm, the creep test data of weak interlayer are fitted and compared. The fitting curve is in good agreement with the test data, which shows the rationality and applicability of the nonlinear creep model. The results show that the nonlinear damage creep model constructed in this paper can well describe the creep characteristics of the weak interlayer and the model has important theoretical reference significance for the study of long-term stability of slope with the weak interlayer.


---

## Body

## 1. Introduction

Rock rheological effect is a common phenomenon in geotechnical engineering. A large amount of slope engineering and tunnel engineering damage is caused by rock rheological effect [1–5]. As a special structural plane, the weak interlayer has low mechanical strength and obvious rheological effect, which often constitutes the weak link in the slope, so it poses a serious threat to the slope stability [6–10]. Therefore, it is necessary to study the rheological mechanical properties of the weak interlayer, and the study of the creep constitutive model of the weak interlayer is the core content [11–15]. Therefore, the study of the creep constitutive model of the weak interlayer has important theoretical significance and practical value for ensuring the long-term stability of slope engineering [16–18].Generally speaking, there are two methods to establish the rheological constitutive model: the first one is directly fitting the rock rheological test curve with the empirical equation through the rheological test of rock. This method has a good fitting effect, but the physical meaning of the model is not clear. The second is based on the rheological test results, which is composed of series and parallel combinations of traditional model components, and then, the unknown rheological model component parameters are determined by identifying the component model and parameter inversion method [6]. Xia et al. [19] established a unified rheological mechanical model including 15 rheological mechanical properties. Nevertheless, since the traditional rheological model is composed of linear components, no matter how many components are in the model, the model is more complex. The final model can only reflect the characteristics of linear viscoelastoplasticity and cannot describe the accelerated rheological stage [20]. Therefore, more and more nonlinear rheological models are proposed. Yang et al. [21] proposed a new nonlinear rheological element NRC model by assuming that the nonlinear shear rheological model of rock is a Weibull distribution function of time and combined it with the time function to describe the accelerated rheological stage. Zhao et al. [22] established a new viscoelastic-plastic creep damage model by combining the Burgers model and nonlinear Mohr–Coulomb plastic element in series. Xu et al. [23] defined the piecewise function of the greenschist creep damage variable changing with time by analyzing the whole creep curve of greenschist, proposed the conjecture that creep damage only appeared in the accelerated rheological stage, and established the generalized Bingham model considering the damage variable. Zhang et al. [24] improved the NRC model proposed by Xu et al. combined with the traditional Kelvin model to form a four-element nonlinear viscoelastic-plastic rheological model. Zhu et al. [25] and Huang et al. [26] assumed the creep damage variable of rock as a negative exponential function and used the damage Burgers model to describe the creep characteristics of the rock.A large number of research results have been achieved in the nonlinear creep model of the rock. However, the research on the creep model of the weak interlayer is relatively rare, and there are few reports on the construction of a nonlinear creep model based on the double thresholds of stress and time. In view of this, based on the series-parallel connection of classical components, this paper introduces the damage variable in the accelerated rheological stage and considers the influence of the time threshold to establish the damage constitutive model that describes the nonlinear rheological properties of the weak interlayer. On this basis, the creep test results of the weak interlayer are nonlinearly fitted to verify the rationality and applicability of the constructed model.


## 2. Establishment of the Nonlinear Creep Model of the Weak Interlayer

### 2.1. Creep Curve Characteristics of the Weak Interlayer

The creep of the weak interlayer is a complex process in which multiple deformations such as elasticity, viscosity, plasticity, viscoelasticity, and viscoplasticity coexist. When the stress of the weak interlayer is less than its long-term strength, the creep curve of the weak interlayer is as shown in Figure1. The weak interlayer generates elastic strain ε0 at the moment of loading and then enters the attenuation creep stage. The creep deformation increases continuously, while the creep rate decreases continuously. At time ta, the creep rate attenuates to zero, and the strain is stable at εa. In this case, the classical element combination model can be used to describe the creep characteristics of the rock. When the stress of the weak interlayer exceeds its long-term strength, the creep curve of the weak interlayer is as shown in Figure 2. At the initial stage of loading, the creep characteristics of the weak interlayer are the same as those mentioned above. They all go through the instantaneous elasticity first and then enter the attenuation creep stage. However, starting from the ta moment, the weak interlayer enters the constant creep stage, the strain of the weak interlayer is still increasing at this stage and the creep rate is a constant value. When the time reaches ts, the weak interlayer enters the accelerated creep stage and εs is the critical strain value of accelerated creep initiation. Since the traditional creep components are linear components, the classical component combination model cannot describe the nonlinear characteristics of accelerated creep. In this paper, the classical components are improved to construct a creep model that can reflect the nonlinear characteristics of the weak interlayer.Figure 1
Steady-state creep curve.Figure 2
Creep failure curve.


### 2.2. Nonlinear Viscoplastic Body Based on Double Threshold Conditions

The viscoplastic body composed of classical elements is shown in Figure3, which is composed of a plastic switch element and a viscous element in parallel. In the figure, τ is the shear stress, η is the viscosity coefficient of the viscous element, and τs is the long-term shear strength of the rock. When the stress of the plastic element is not more than τs, the plastic element is not opened and the viscoplastic body does not produce strain. When the stress of the plastic element is more than τs, the plastic element is opened and the viscoplastic body produces strain. Since the viscous element of the viscoplastic body is a linear element, it cannot describe the nonlinear accelerated creep stage of rock, so it needs to be improved. The idea of improvement is to conduct nonlinear treatment of viscous components. A large number of tests show that when the stress of rock is less than its long-term strength, the rock is in the first two stages of creep and no accelerated creep failure occurs. When the stress of the rock exceeds its long-term strength, it is not immediately leading to accelerated creep failure, but at a certain time ts, the rock enters the accelerated creep stage. According to the current relevant research [6], the aging damage of rock will occur when it enters accelerated creep. Therefore, based on the double threshold conditions of stress and time and considering the influence of aging damage, this paper improves the classical viscoplastic body, and the improved nonlinear viscoplastic body is shown in Figure 4. In the figure, D is the damage variable, ts is the start time of accelerated creep, and the meaning of other parameters is the same as that of the corresponding parameters in Figure 3.Figure 3
Classical viscoplastic body.Figure 4
Nonlinear viscoplastic body.The damage variableD can be expressed as follows:(1)D=0t≤ts1−1aexp−at−tst>ts..Here,a is the material parameter, which can be determined by fitting test data; t is the creep time, and ts is the start time of accelerated creep. The damage variable D value is 0∼1. If t ≤ ts, then D = 0; that is, the weak interlayer is not damaged, When t approaches infinity, D = 1, indicating that the weak interlayer has been destroyed.Based on the above analysis, we can construct the creep equation of a nonlinear viscoplastic body.Whenτ ≤ τs, the plastic switch of the nonlinear viscoplastic body is closed and there is no strain in the viscoplastic body; that is, ε = 0. When τ ＞ τs and t ≤ ts, the plastic switch of the nonlinear viscoplastic body will be opened, but it has not entered the acceleration stage and the aging damage has not occurred. We regard it as a classical viscoplastic body, so its creep equation can be expressed as follows: (2)ε=τ−τsηt.Whenτ＞τs and t＞ts, the plastic switch has been opened and the aging damage has also occurred. The nonlinear viscoplastic body enters the accelerated creep stage. The constitutive equation of the nonlinear viscoplastic body can be expressed as follows:(3)τ=τs+η1−Dε˙=τs+η1aexp−at−tsε˙.From formula (3),(4)ε˙=τ−τsηaexpat−ts.The creep equation of the nonlinear viscoplastic body is obtained by integrating both sides of (4):(5)ε=τ−τsηexpat−ts.


### 2.3. Determination of Accelerated Creep Time Threshold

From the previous section, we can see that the accelerated creep time thresholdts reflects the opening time of the accelerated creep of the weak interlayer. It is an important parameter of the creep characteristics of the weak interlayer, which can be determined by the following method.According to the research of some scholars [27], the failure load of the weak interlayer decreases with the increase in failure time, as shown in Figure 5. τ0 is the instantaneous strength of the weak interlayer, τ∞ is the long-term strength of the weak interlayer, and the stress corresponding to the accelerated creep start time ts is τ. When the stress of the weak interlayer is higher than its long-term strength, the relationship between accelerated creep time threshold and stress can be established from Figure 5, as shown in the following equation:Figure 5
Relationship between failure load and time.(6)τ=exp−αts+β.From formula (6),(7)ts=β−lnτα.Here,α and β are undetermined parameters, which can be determined by fitting test data. The accelerated creep time threshold of the weak interlayer can be determined by formula (7).


### 2.4. Establishment of the Nonlinear Creep Model

According to the creep characteristics of the weak interlayer, when the stress is less than the stress threshold of accelerated creep, the weak interlayer generates instantaneous elastic strain, decay creep, and steady creep. At this time, the Bergs model can be used to describe the creep characteristics of the weak interlayer. When the stress exceeds the stress threshold of accelerated creep, the rock will eventually enter the nonlinear accelerated creep stage at a certain time point after experiencing the creep deformation of the first two stages. A nonlinear viscoelastic-plastic damage creep model which can describe the whole creep process of the weak interlayer can be constructed by using the improved nonlinear viscoplastic body and Bergs body in series. The model is shown in Figure6. In the figure, I describes the instantaneous elastic strain of the weak interlayer, II reflects the viscoelasticity of the weak interlayer, III reflects the viscosity of the weak interlayer, II and III describe the attenuation creep stage and steady creep stage of the weak interlayer, and IV reflects the nonlinear viscoplasticity of the weak interlayer, which describes the nonlinear accelerated creep stage of the weak interlayer.Figure 6
Nonlinear damage creep model.It can be seen from Figure6 that when τ ≤ τs, I, II, and III are all involved in creep deformation. The state equations of the creep model are as follows:(8)τ=τ1=τ2=τ3ε=ε1+ε2+ε3τ1=G1ε1τ2=τ21+τ22=G2ε2+η1ε˙2τ3=η2ε˙3..From formula (8),(9)ε=τG1+τG21−exp−G2η1t+τη2t.Whenτ＞τs, and t ≤ ts, I, II, III, and IV all participate in creep deformation; however, the viscoplastic body has not entered the accelerated creep stage, and it has not been damaged. At this time, the state equation of the creep model is as follows:(10)τ=τ1=τ2=τ3=τ4ε=ε1+ε2+ε3+ε4τ1=G1ε1τ2=τ21+τ22=G2ε2+η1ε˙2τ3=η2ε˙3τ4=τs+η3ε˙4..From formula (10),(11)ε=τG1+τG21−exp−G2η1t+τη2t+τ−τsη3t.Whenτ＞τs, and τ＞τs, I, II, III, and IV all participate in creep deformation, the viscoplastic body has been damaged, and it has entered the accelerated creep stage. At this time, the state equation of the creep model is as follows:(12)τ=τ1=τ2=τ3=τ4ε=ε1+ε2+ε3+ε4τ1=G1ε1τ2=τ21+τ22=G2ε2+η1ε˙2τ3=η2ε˙3τ4=τs+η31−Dε˙4=τs+η31aexp−at−tsε˙4..From formula (12),(13)ε=τG1+τG21−exp−G2η1t+τη2t+τ−τsη3expat−ts.Therefore, the nonlinear damage creep equation of the weak interlayer is as follows:(14)ε=τG1+τG21−exp−G2η1t+τη2tτ≤τs,ε=τG1+τG21−exp−G2η1t+τη2t+τ−τsη3tτ<τs,t≤tsε=τG1+τG21−exp−G2η1t+τη2t+τ−τsη3expat−tsτ<τs,t≤ts..


## 2.1. Creep Curve Characteristics of the Weak Interlayer

The creep of the weak interlayer is a complex process in which multiple deformations such as elasticity, viscosity, plasticity, viscoelasticity, and viscoplasticity coexist. When the stress of the weak interlayer is less than its long-term strength, the creep curve of the weak interlayer is as shown in Figure1. The weak interlayer generates elastic strain ε0 at the moment of loading and then enters the attenuation creep stage. The creep deformation increases continuously, while the creep rate decreases continuously. At time ta, the creep rate attenuates to zero, and the strain is stable at εa. In this case, the classical element combination model can be used to describe the creep characteristics of the rock. When the stress of the weak interlayer exceeds its long-term strength, the creep curve of the weak interlayer is as shown in Figure 2. At the initial stage of loading, the creep characteristics of the weak interlayer are the same as those mentioned above. They all go through the instantaneous elasticity first and then enter the attenuation creep stage. However, starting from the ta moment, the weak interlayer enters the constant creep stage, the strain of the weak interlayer is still increasing at this stage and the creep rate is a constant value. When the time reaches ts, the weak interlayer enters the accelerated creep stage and εs is the critical strain value of accelerated creep initiation. Since the traditional creep components are linear components, the classical component combination model cannot describe the nonlinear characteristics of accelerated creep. In this paper, the classical components are improved to construct a creep model that can reflect the nonlinear characteristics of the weak interlayer.Figure 1
Steady-state creep curve.Figure 2
Creep failure curve.


## 2.2. Nonlinear Viscoplastic Body Based on Double Threshold Conditions

The viscoplastic body composed of classical elements is shown in Figure3, which is composed of a plastic switch element and a viscous element in parallel. In the figure, τ is the shear stress, η is the viscosity coefficient of the viscous element, and τs is the long-term shear strength of the rock. When the stress of the plastic element is not more than τs, the plastic element is not opened and the viscoplastic body does not produce strain. When the stress of the plastic element is more than τs, the plastic element is opened and the viscoplastic body produces strain. Since the viscous element of the viscoplastic body is a linear element, it cannot describe the nonlinear accelerated creep stage of rock, so it needs to be improved. The idea of improvement is to conduct nonlinear treatment of viscous components. A large number of tests show that when the stress of rock is less than its long-term strength, the rock is in the first two stages of creep and no accelerated creep failure occurs. When the stress of the rock exceeds its long-term strength, it is not immediately leading to accelerated creep failure, but at a certain time ts, the rock enters the accelerated creep stage. According to the current relevant research [6], the aging damage of rock will occur when it enters accelerated creep. Therefore, based on the double threshold conditions of stress and time and considering the influence of aging damage, this paper improves the classical viscoplastic body, and the improved nonlinear viscoplastic body is shown in Figure 4. In the figure, D is the damage variable, ts is the start time of accelerated creep, and the meaning of other parameters is the same as that of the corresponding parameters in Figure 3.Figure 3
Classical viscoplastic body.Figure 4
Nonlinear viscoplastic body.The damage variableD can be expressed as follows:(1)D=0t≤ts1−1aexp−at−tst>ts..Here,a is the material parameter, which can be determined by fitting test data; t is the creep time, and ts is the start time of accelerated creep. The damage variable D value is 0∼1. If t ≤ ts, then D = 0; that is, the weak interlayer is not damaged, When t approaches infinity, D = 1, indicating that the weak interlayer has been destroyed.Based on the above analysis, we can construct the creep equation of a nonlinear viscoplastic body.Whenτ ≤ τs, the plastic switch of the nonlinear viscoplastic body is closed and there is no strain in the viscoplastic body; that is, ε = 0. When τ ＞ τs and t ≤ ts, the plastic switch of the nonlinear viscoplastic body will be opened, but it has not entered the acceleration stage and the aging damage has not occurred. We regard it as a classical viscoplastic body, so its creep equation can be expressed as follows: (2)ε=τ−τsηt.Whenτ＞τs and t＞ts, the plastic switch has been opened and the aging damage has also occurred. The nonlinear viscoplastic body enters the accelerated creep stage. The constitutive equation of the nonlinear viscoplastic body can be expressed as follows:(3)τ=τs+η1−Dε˙=τs+η1aexp−at−tsε˙.From formula (3),(4)ε˙=τ−τsηaexpat−ts.The creep equation of the nonlinear viscoplastic body is obtained by integrating both sides of (4):(5)ε=τ−τsηexpat−ts.


## 2.3. Determination of Accelerated Creep Time Threshold

From the previous section, we can see that the accelerated creep time thresholdts reflects the opening time of the accelerated creep of the weak interlayer. It is an important parameter of the creep characteristics of the weak interlayer, which can be determined by the following method.According to the research of some scholars [27], the failure load of the weak interlayer decreases with the increase in failure time, as shown in Figure 5. τ0 is the instantaneous strength of the weak interlayer, τ∞ is the long-term strength of the weak interlayer, and the stress corresponding to the accelerated creep start time ts is τ. When the stress of the weak interlayer is higher than its long-term strength, the relationship between accelerated creep time threshold and stress can be established from Figure 5, as shown in the following equation:Figure 5
Relationship between failure load and time.(6)τ=exp−αts+β.From formula (6),(7)ts=β−lnτα.Here,α and β are undetermined parameters, which can be determined by fitting test data. The accelerated creep time threshold of the weak interlayer can be determined by formula (7).


## 2.4. Establishment of the Nonlinear Creep Model

According to the creep characteristics of the weak interlayer, when the stress is less than the stress threshold of accelerated creep, the weak interlayer generates instantaneous elastic strain, decay creep, and steady creep. At this time, the Bergs model can be used to describe the creep characteristics of the weak interlayer. When the stress exceeds the stress threshold of accelerated creep, the rock will eventually enter the nonlinear accelerated creep stage at a certain time point after experiencing the creep deformation of the first two stages. A nonlinear viscoelastic-plastic damage creep model which can describe the whole creep process of the weak interlayer can be constructed by using the improved nonlinear viscoplastic body and Bergs body in series. The model is shown in Figure6. In the figure, I describes the instantaneous elastic strain of the weak interlayer, II reflects the viscoelasticity of the weak interlayer, III reflects the viscosity of the weak interlayer, II and III describe the attenuation creep stage and steady creep stage of the weak interlayer, and IV reflects the nonlinear viscoplasticity of the weak interlayer, which describes the nonlinear accelerated creep stage of the weak interlayer.Figure 6
Nonlinear damage creep model.It can be seen from Figure6 that when τ ≤ τs, I, II, and III are all involved in creep deformation. The state equations of the creep model are as follows:(8)τ=τ1=τ2=τ3ε=ε1+ε2+ε3τ1=G1ε1τ2=τ21+τ22=G2ε2+η1ε˙2τ3=η2ε˙3..From formula (8),(9)ε=τG1+τG21−exp−G2η1t+τη2t.Whenτ＞τs, and t ≤ ts, I, II, III, and IV all participate in creep deformation; however, the viscoplastic body has not entered the accelerated creep stage, and it has not been damaged. At this time, the state equation of the creep model is as follows:(10)τ=τ1=τ2=τ3=τ4ε=ε1+ε2+ε3+ε4τ1=G1ε1τ2=τ21+τ22=G2ε2+η1ε˙2τ3=η2ε˙3τ4=τs+η3ε˙4..From formula (10),(11)ε=τG1+τG21−exp−G2η1t+τη2t+τ−τsη3t.Whenτ＞τs, and τ＞τs, I, II, III, and IV all participate in creep deformation, the viscoplastic body has been damaged, and it has entered the accelerated creep stage. At this time, the state equation of the creep model is as follows:(12)τ=τ1=τ2=τ3=τ4ε=ε1+ε2+ε3+ε4τ1=G1ε1τ2=τ21+τ22=G2ε2+η1ε˙2τ3=η2ε˙3τ4=τs+η31−Dε˙4=τs+η31aexp−at−tsε˙4..From formula (12),(13)ε=τG1+τG21−exp−G2η1t+τη2t+τ−τsη3expat−ts.Therefore, the nonlinear damage creep equation of the weak interlayer is as follows:(14)ε=τG1+τG21−exp−G2η1t+τη2tτ≤τs,ε=τG1+τG21−exp−G2η1t+τη2t+τ−τsη3tτ<τs,t≤tsε=τG1+τG21−exp−G2η1t+τη2t+τ−τsη3expat−tsτ<τs,t≤ts..


## 3. Verification of the Creep Model

By fitting the creep equation derived in the above section to the creep test data of the weak interlayer, the rationality and applicability of the model constructed in this paper can be verified. In this paper, data of Zhu et al. [28] were used to conduct shear creep test of the weak interlayer. In the first group of tests, creep tests were carried out on the rock samples with weak interlayers under the normal stress of 0.5 MPa by applying the shear stress step by step. The shear stresses applied at all levels were 0.10 MPa, 0.19 MPa, 0.29 MPa, 0.39 MPa, and 0.59 MPa, respectively. According to the test results, the long-term shear strength of the weak interlayer under this state was 0.423 MPa. In the second group, the creep test was carried out on the rock samples of weak interlayers under a normal stress of 0.7 MPa by applying shear stress in different grades. The applied shear stresses at different levels were 0.216 MPa, 0.432 MPa, 0.648 MPa, 0.864 MPa, 1.080 MPa, and 1.296 MPa, respectively. According to the test results, the long-term shear strength of the weak interlayer under this condition was 1.080 MPa. Based on the test results, using the custom function fitting tool of Origin software and Levenberg–Marquardt optimization algorithm, the creep test data of weak interlayers under different shear stresses are fitted and analyzed.It should be pointed out that for the three stages of complete creep, because the creep equation is more complex, the creep parameters are more and the initial value of the parameters to be optimized is difficult to determine; if the whole method is used for fitting, the results will be as shown in Figure7, indicating that the fitting has failed. Therefore, the fitting method needs to be improved. Firstly, the complete creep curve is divided into two parts, which are the first two stages of the creep curve and the accelerated creep curve. They are fitted, respectively, and the initial values of each creep parameter are obtained, as shown in Figure 8 and Figure 9. Secondly, by using these initial creep parameters and using the integral method to fit the complete creep curve, the ideal fitting effect can be obtained, as shown in Figure 10.Figure 7
Fitting of the integral method.Figure 8
Fitting of the first two creep stages.Figure 9
Fitting of the accelerated creep stage.Figure 10
Fitting of the improved integral method.This method is used to fit the creep test data of these two groups of weak interlayers. The fitting curves are shown in Figures11 and 12. The creep model parameters are obtained by fitting, as shown in Tables 1 and 2.Figure 11
Experimental data and the fitting curves (σ = 0.5 MPa; τ = 0.10∼0.59 MPa).Figure 12
Experimental data and the fitting curves (σ = 0.7 MPa; τ = 0.216∼1.296 MPa).Table 1
Creep parameter identification results (σ = 0.5 MPa; τ = 0.10∼0.59 MPa).
τ/MPaG1/MPaG2/MPaη1/MPa·mm-1 hη2/MPa·mm-1 hη3/MPa·mm-1 haR20.13.012.981.53127.340.99370.191.992.494.77239.540.98510.291.553.381.64446.510.98930.391.244.211.66505.750.97580.591.341.641.6579.415.610.510.9963Table 2
Creep parameter identification results (σ = 0.7 MPa; τ = 0.216∼1.296 MPa).
τ/MPaG1/MPaG2/MPaη1/MPa·mm-1 hη2/MPa·mm-1 hη3/MPa·mm-1 haR20.21649.1110.3453.49402.270.99420.43210.0113.9161.043256.510.99540.6488.1511.6128.679076.810.99680.8645.7911.3434.4421063.370.99891.0804.6921.2145.6114321.180.99831.2964.4725.0558.413.938.780.380.9931From Figures11 and 12, it can be seen that the creep test data of weak interlayers under various loads are in good agreement with their fitting curves and the correlation coefficients in Tables 1 and 2 are basically above 0.95, which indicates that the nonlinear damage creep model constructed in this paper can well describe the instantaneous deformation, attenuation creep stage, steady creep stage, and accelerated creep stage of the weak interlayer, which further illustrates the rationality and applicability of the model.


## 4. Conclusions

(1)
Based on the double threshold conditions of stress and time and considering the influence of aging damage, the classical viscoplastic body is improved in this paper. The improved nonlinear viscoplastic body can more accurately reflect the characteristics of the accelerated creep stage.(2)
In the fitting analysis of creep curves containing complete three stages, the effect of the complete method is usually poor. In this paper, a piecewise fitting method with a good fitting effect is innovatively proposed.(3)
The creep test data of weak interlayers are fitted and analyzed by the creep model constructed in this paper. The results show that the fitting curve is in good agreement with the experimental data, indicating that the nonlinear damage creep model constructed in this paper can well describe the creep characteristics of the weak interlayer. This model can provide important theoretical support for the study of long-term stability of slopes with weak interlayers.



---

*Source: 3566521-2022-04-08.xml*