You do not have permission to edit this page, for the following reason:

You are not allowed to execute the action you have requested.


You can view and copy the source of this page.

x
 
1
<!-- metadata commented in wiki content
2
3
4
'''A BP neural network-based micro particle parameters calibration and an energy criterion for MatDEM coupled with strength reduction method to evaluate 3D slope stability'''
5
6
Wei Jiang*<sup>1, 2</sup>, Yihong Tan<sup>1</sup>, Jinzhou Yan<sup>3</sup>, Ye Ouyang<sup>1</sup>, Zhaoyu Fu<sup>1</sup>, Qiang Feng<sup>1</sup>
7
8
1.Hubei Key Laboratory of Disaster Prevention and Mitigation, China Three Gorges University, Yichang, Hubei, 443002, China
9
10
2. Key Laboratory of Geological Hazards on Three Gorges Reservoir Area, Ministry of Education, Yichang, Hubei, 443002, China
11
12
3. Yichang Construction Investment & Development Co., Ltd., Yichang, Hubei, 443002, China
13
14
Correspondence: Wei Jiang ([mailto:jiangweilion@ctgu.edu.cn jiangweilion@ctgu.edu.cn])
15
-->
16
17
'''Abstract: '''To enhance the applicability of discrete element method in 3D slope stability analysis, a BP neural network-based micro parameter calibration method and an energy criterion are proposed by taking MatDEM as an example. Firstly, the relationship between the micro particle parameters and the shear strengths of particle aggregate are represented by using the BP neural network. And then the micro particle parameters are obtained for the given shear strengths by using a correction calibration. Next, the energy conversions are investigated for the stable and instable slope models in MatDEM. From a view of practical application, the abrupt in variation tendency and magnitude of the kinetic energy is selected for indicating the emergence of the limit equilibrium state of a slope. Finally, the effectiveness of the proposed improvements is testified by taking Baijiabao landslide as an example. Results verify that the calibration method established in this study is applicable to provide the micro particle parameters when the shear strength is constantly reduced, and the factor of safety determined by the kinetic energy criterion reflects the landslide stability at the global level.
18
19
'''Key words:''' 3D slope stability analysis, MatDEM, strength reduction method, micro parameters, failure criterion
20
21
=1 Introduction=
22
23
Slope stability analysis is a classic problem in geotechnical engineering. Three dimensional (3D) analyses has some significant advantages over two dimensional (2D) analyses in the slope stability analysis that the geometry structure, boundary conditions and spatial variations of a slope is fully took into account. Therefore 3D slope analysis has become a popular research topic in recent years [1]. The limit equilibrium method (LEM) mainly developed in the last century is still the most popular method among engineers. While, its common demerits are that it involves many assumptions of force and it is extremely difficult to locate a critical slip surface [2]. Additionally, it is hard to generalize the procedure of LEM from 2D to 3D [3].
24
25
With the advances of computers, numerical analysis methods have become affordable to engineers, such as the finite element method (FEM) [4-5], boundary element method (BEM) [6], virtual element method (VEM) [7], discrete element method (DEM) [8-9], discontinuous deformation analysis (DDA) [10], numerical manifold method (NMM) [11], material point method (MPM) [12], smoothed particle hydrodynamics (SPH) [13], and some combined methods [14-15]. These numerical methods fall into two categories, continuum-based and discontinuum-based numerical methods. The latter kind of methods, e.g. DEM, overcome the limitation of the assumption of macro-continuity in the former kind of methods, and are capable of simulating the generation and development process of the slip surface. Due to these advantages, DEM is becoming increasing popular for slope stability analysis in situation where the continuum-based numerical model is inapplicable [16-18].
26
27
When analyzing the slope stability using numerical methods, the strength reduction method (SRM) [19] is a commonly used strategy for determining the factor of safety (FOS) of a slope. The SRM was initially proposed for FEM analysis of slope stability [20]. SRM continuously reduces the shear strength of soil and rock until a slope reaches a limit equilibrium, and then the reduction factor is considered as the FOS of a slope. Theoretically, SRM can be combined with DEM similarity to the coupling of FEM and SRM. But, two critical issues are expected to be addressed when analyzing the slope stability by integrating SRM and DEM.
28
29
The first critical issue is how to reduce the shear strength of soil and rock in the DEM modeling. DEM simulates the mechanical properties of soil and rock from the micro point of view, and the relationship between the macro parameters of soil and rock and the micro parameters of particles is complicated. Theoretical and numerical approaches were used to reveal this relationship for soil and rock simulated by different models, e.g. the bonded particle model [21-22], the flat jointed model [23], the close-packed lattice model [24], and the coarse grained particle model [25]. But, a rigid theory is still not available to ensure the magnitude of the relationship between macro parameters and micro parameters [16]. The trial-and-error calibration remains the usual approach to obtain the micro parameters of particles corresponding to the macro parameters [26], even if it exhibits high computing costs. In slope stability analysis by integrating SRM and DEM, the shear strength is constantly reduced, which means that the micro parameters are necessary to be adjusted correspondingly. To avoid the cumbersome calibration on the micro parameters, the shear strength reduction is usually simplified as the reduction of several micro parameters [27-28] based on effect analysis of micro parameter on macro property. The simplification may results in the FOS deviating from its original definition, in view of the complicated relationship between the shear strength and the micro parameters. Therefore, a reliable approach is expected to reduce the shear strength of soil and rock in the DEM modeling.
30
31
The second critical issue is how to determine the FOS when integrating SRM and DEM to analyze the slope stability. The target is to establish an appropriate criterion to indicate that a slope is in the limit equilibrium state. When analyzing the slope stability by using FEM with SRM, whether the slope is in the limit equilibrium state is usually evaluated by using three criteria: (1) penetration of the plastic zone in the slope, (2) convergence of the unbalance force, and (3) an abrupt change in the displacements of selected characteristic points. These criteria cannot be simply duplicated into the slope stability analysis by using DEM with SRM. In many DEM modelings, the particles are treated as rigid bodies, and thus there will be no plastic zone, which induces that it is impossible to employ the first criterion for the DEM analysis of the slope stability. In case some particles fall or a small region of the slope is deformed, the unbalance forces in the DEM modeling will not convergence, while actually the slope is still stable overall. Thus, it is unreasonable to judge whether the slope is unstable based on the second criterion [29]. Similarly, in views of the fact that the displacement characteristic of some selected points cannot reflect the overall behavior of the entire slope, it is improper to determine whether the entire slope is unstable based on the third criterion [30]. To avoid the limitations of the traditional criteria, the abrupt change in the vector graphic of global displacements [31], the variation coefficient of global displacements [29], and the slope failure extent [16] have been used as the failure criterion in the slope stability analysis by using DEM with SRM. These criteria are capable to represent the slope stability at the global level, but they are inapplicable in some cases. For example, the vector graphic of global displacements may be quite disordered for a complex slope, and it is difficult to evaluate significant changes. Thus, a reasonable criterion is expected to determine the FOS when analyzing the slope stability by using DEM and SRM.
32
33
The fast GPU matrix computing discrete element method (MatDEM) [32] is a newly developed code based on the DEM theory and MATLAB, which has been verified as an effective numerical simulations by many examples [33-36]. Since matrix operations and high-performance GPU calculations are used in the code, its computational efficiency is significantly improved. By taking MatDEM as an example, two improvements are proposed in this study to address the aforementioned two issues. Firstly, to precisely represent the reduction of the shear strengths of the geological material in the MatDEM modeling, the relationship between the macro shear strengths and the micro particle parameters is represented by an overdetermined BP neural network, and a correction calibration is developed to obtain the micro particle parameters corresponding to the original and reduced shear strengths. Then, to reasonably determine the FOS of a slope at the global level, the energy conversion and heat generation during the failure of a slope are fully investigated, and thus an energy criterion is developed for MatDEM coupled with SRM to evaluate the slope stability. Finally, the validity of the proposed strategies is testified by taking Baijiabao landslide situated in Zigui County, Hubei Province of China as a case study.
34
35
=2 Fundamental principles of the MatDEM=
36
37
==2.1 Contact model of particles==
38
39
The basic model of the MatDEM code is a series of 3D close-packed elastic particles as shown in Fig. 1a, and these particles interact through the spring forces (Fig. 1b).
40
41
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
42
 [[Image:Review_852894090255-image1.png|252px]] </div>
43
44
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
45
<span style="text-align: center; font-size: 75%;">(a) 3D close-packed elastic particles</span></div>
46
47
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
48
 [[Image:Review_852894090255-image2.png|474px]] </div>
49
50
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
51
<span style="text-align: center; font-size: 75%;"> (b) Normal and tangential springs between two particles</span></div>
52
53
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
54
<span style="text-align: center; font-size: 75%;">Fig.1 The contact model of particles in MatDEM</span></div>
55
56
The normal force ''F<sub>n</sub>'' and the normal relative displacement ''X<sub>n</sub>'' between two adjacent particles are simulated by a normal spring:
57
58
{| class="formulaSCP" style="width: 100%; text-align: center;" 
59
|-
60
| 
61
{| style="text-align: center; margin:auto;" 
62
|-
63
| <math>F_n=\lbrace \begin{array}{c}
64
K_nX_n\\
65
0\\
66
K_nX_n
67
\end{array}\mbox{ }\begin{array}{c}
68
X_n\leq X_b\mbox{, intact }\\
69
X_n\geq 0,\mbox{ broken}\\
70
X_n<0,\mbox{ broken}
71
\end{array}\mbox{ }\begin{array}{c}
72
\mbox{(a)}\\
73
\mbox{(b)}\\
74
\mbox{(c)}
75
\end{array}</math>
76
|}
77
| style="width: 5px;text-align: right;white-space: nowrap;" | (1)
78
|}
79
80
81
in which ''K<sub>n</sub>'' is the stiffness of the normal spring, and ''X''<sub>b</sub> is the breaking displacement. Initially, the particles are interconnected with their adjacent particles and subjected to tensile or compressive spring forces as Eq. 1a. In case of ''X<sub>n</sub>'' between the particle pair exceeding ''X''<sub>b</sub>, the spring breaks and the tensile force ceases to exist between them (Eq. 1b). However, the compressive force may act between them when they return to a compressive status (Eq. 1c). It is noticed that the tensile force is positive, and the compressive force is negative in Eq. 1.
82
83
The shear force ''F''<sub>s</sub>'' ''and the tangential relative displacement ''X''<sub>s</sub> between two particles are simulated by a tangential spring:
84
85
{| class="formulaSCP" style="width: 100%; text-align: center;" 
86
|-
87
| 
88
{| style="text-align: center; margin:auto;" 
89
|-
90
| <math>F_\mbox{s}=K_sX_s</math>
91
|}
92
| style="width: 5px;text-align: right;white-space: nowrap;" | (2)
93
|}
94
95
96
where ''K''<sub>s</sub> is the stiffness of the tangential spring. The tangential spring also has a fracture criterion based on the Mohr-Coulomb yield criterion:
97
98
{| class="formulaSCP" style="width: 100%; text-align: center;" 
99
|-
100
| 
101
{| style="text-align: center; margin:auto;" 
102
|-
103
| <math>F_n=\lbrace \begin{array}{c}
104
F_s0X_n\\
105
0\\
106
K_nX_n
107
\end{array}\mbox{ }\begin{array}{c}
108
X_n\leq X_b\mbox{, intact }\\
109
X_n\geq 0,\mbox{ broken}\\
110
X_n<0,\mbox{ broken}
111
\end{array}\mbox{ }\begin{array}{c}
112
\mbox{(a)}\\
113
\mbox{(b)}\\
114
\mbox{(c)}
115
\end{array}</math>
116
|}
117
| style="width: 5px;text-align: right;white-space: nowrap;" | (3)
118
|}
119
120
121
where ''F<sub>s</sub>''<sub>max</sub> is the inter-particle shear resistance, ''F<sub>s0</sub>'' is the initial inter-particle shear resistance when no normal force exists, and ''&#x03bc;<sub>p</sub>'' is the inter-particle coefficient of friction. When the tangential force ''F<sub>s</sub>'' exceeds ''F<sub>s</sub>''<sub>max</sub> in Eq. 3a, the tangential spring will break, whereupon the shear force ''F''<sub>s</sub> is limited to be less than or equal to the inter-particle shear resistance ''F<sub>s</sub>''<sub>max</sub> of the broken spring. If the tangential spring is broken and the magnitude of external shear force exceeds the limit ''F<sub>s</sub>''<sub>max</sub> in Eq. 3b, two particles begin slipping, and the slipping friction between them is ''F<sub>s</sub>''<sub>max</sub> in Eq. 3b.
122
123
==2.2 Energy system of the MatDEM modeling==
124
125
The MatDEM modeling follows the law of energy conversation and Newton’s law of motion. The energies of a model fall into two categories, mechanical energy and heat. A brief introduction of their compositions is given as follows.
126
127
(1) The mechanical energy in the system consists of elastic potential energy, gravitational potential energy and kinetic energy. The elastic potential energy ''E''<sub>e</sub> is the sum of the strain energy of normal and tangential springs between particles:
128
129
{| class="formulaSCP" style="width: 100%; text-align: center;" 
130
|-
131
| 
132
{| style="text-align: center; margin:auto;" 
133
|-
134
| <math>E_e=\frac{1}{2}K_nX_n^2+\frac{1}{2}K_sX_s^2</math>
135
|}
136
| style="width: 5px;text-align: right;white-space: nowrap;" | (4)
137
|}
138
139
140
The gravitational potential energy ''E''<sub>g</sub> of a particle is:
141
142
{| class="formulaSCP" style="width: 100%; text-align: center;" 
143
|-
144
| 
145
{| style="text-align: center; margin:auto;" 
146
|-
147
| <math>E_g=mgh</math>
148
|}
149
| style="width: 5px;text-align: right;white-space: nowrap;" | (5)
150
|}
151
152
153
where ''m'' is the particle mass, ''g'' denotes the gravity acceleration, and ''h'' represents the height above the reference level. The kinetic energy ''E''<sub>k</sub> of a particle is computed as:
154
155
{| class="formulaSCP" style="width: 100%; text-align: center;" 
156
|-
157
| 
158
{| style="text-align: center; margin:auto;" 
159
|-
160
| <math>E_k=\frac{1}{2}mv^2</math>
161
|}
162
| style="width: 5px;text-align: right;white-space: nowrap;" | (6)
163
|}
164
165
166
in which ''v'' is the scalar of particle velocity.
167
168
(2) The heat in the system comes from the viscous damping, the spring breaking and the friction. Damping is used in the MatDEM to weaken the elastic wave energy of the model and dissipate the kinetic energy in the system. The damping force '''F'''<sub>d</sub> is given by
169
170
{| class="formulaSCP" style="width: 100%; text-align: center;" 
171
|-
172
| 
173
{| style="text-align: center; margin:auto;" 
174
|-
175
| <math>\boldsymbol{F}_d=-\eta \boldsymbol{v}</math>
176
|}
177
| style="width: 5px;text-align: right;white-space: nowrap;" | (7)
178
|}
179
180
181
where ''η'' is the damping coefficient and '''v''' is the vector of particle velocity. Because the time step of the simulation is very small, the particle velocity is assumed to be constant in a step. Viscous heat ''Q''<sub>d</sub> generated during viscous damping is calculated by
182
183
{| class="formulaSCP" style="width: 100%; text-align: center;" 
184
|-
185
| 
186
{| style="text-align: center; margin:auto;" 
187
|-
188
| <math>Q_d=-\boldsymbol{F}_d\cdot d\boldsymbol{x}</math>
189
|}
190
| style="width: 5px;text-align: right;white-space: nowrap;" | (8)
191
|}
192
193
194
in which ''d'' '''x''' denotes the particle displacement in the current time step.
195
196
When an intact spring breaks, the spring force reduces, and the elastic potential energy of the inter-particle normal spring and/or tangential spring will dissipate into heat. Thus, breaking heat is equal to the reduction of the elastic potential energy ''E''<sub>e</sub>. The calculation of breaking heat is based on the status of the inter-particle normal force. If the inter-particle normal force is tensile, both the normal and tangential spring forces will reduce to zero, and thus the breaking heat ''Q''<sub>b</sub> is the sum of elastic potential energy of both normal and tangential springs:
197
198
{| class="formulaSCP" style="width: 100%; text-align: center;" 
199
|-
200
| 
201
{| style="text-align: center; margin:auto;" 
202
|-
203
| <math>Q_b=\frac{1}{2}K_nX_n^2+\frac{1}{2}K_sX_s^2</math>
204
|}
205
| style="width: 5px;text-align: right;white-space: nowrap;" | (9)
206
|}
207
208
209
If the normal force is compressive, the normal spring force will not be changed, and the tangential spring force reduced from ''F<sub>s</sub>''<sub>max</sub> in Eq. 3a to ''F<sub>s</sub>''<sub>max</sub> in Eq. 3b. The breaking heat ''Q''<sub>b</sub> is the reduction of the elastic potential energy of the tangential spring, which can be calculated by
210
211
{| class="formulaSCP" style="width: 100%; text-align: center;" 
212
|-
213
| 
214
{| style="text-align: center; margin:auto;" 
215
|-
216
| <math>Q_b=\frac{\left({\left(F_{s0}-{\mu }_pF_n\right)}^2-{\left({\mu }_pF_n\right)}^2\right)}{2K_s}=</math><math>\frac{F_{s0}^2-2{\mu }_pF_nF_{s0}}{2K_s}</math>
217
|}
218
| style="width: 5px;text-align: right;white-space: nowrap;" | (10)
219
|}
220
221
222
Friction heat ''Q''<sub>f</sub> generated during the sliding process is defined as the product of the average sliding friction and the effective sliding distance:
223
224
{| class="formulaSCP" style="width: 100%; text-align: center;" 
225
|-
226
| 
227
{| style="text-align: center; margin:auto;" 
228
|-
229
| <math>Q_f=\vert 0.5\left(\boldsymbol{F}_{s1}+\boldsymbol{F}_{s1}\right)\cdot d\boldsymbol{S}\vert </math>
230
|}
231
| style="width: 5px;text-align: right;white-space: nowrap;" | (11)
232
|}
233
234
235
in which '''F'''<sub>s1</sub> and '''F'''<sub>s2</sub> are the sliding friction forces at respectively the beginning and the end of the current time step, and d'''S''' denotes the effective sliding distance. More details about Eq. 11 can be found in [32].
236
237
The total energy of a MatDEM model is the sum of all mechanical energy and heat:
238
239
{| class="formulaSCP" style="width: 100%; text-align: center;" 
240
|-
241
| 
242
{| style="text-align: center; margin:auto;" 
243
|-
244
| <math>E_{total}=E_e+E_g+E_k+Q</math>
245
|}
246
| style="width: 5px;text-align: right;white-space: nowrap;" | (12)
247
|}
248
249
250
where ''Q'' represents the sum of viscous heat ''Q''<sub>d</sub>, breaking heat ''Q''<sub>b</sub> and friction heat ''Q''<sub>f</sub>:
251
252
{| class="formulaSCP" style="width: 100%; text-align: center;" 
253
|-
254
| 
255
{| style="text-align: center; margin:auto;" 
256
|-
257
| <math>Q=Q_d+Q_b+Q_f</math>
258
|}
259
| style="width: 5px;text-align: right;white-space: nowrap;" | (13)
260
|}
261
262
263
Based on the law of energy conversation, the total energy of an isolated system is constant. If the system deforms under external force, the increment of the total energy must be equal to the work done by the external force.
264
265
=3 A BP neural network-based micro particle parameters calibration=
266
267
==3.1 The structure of the BP neural network to predict the shear strength==
268
269
With the rapid development of neural network technology, recently many researchers attempted to establish an intelligent method to represent the relationship between macro and micro parameters for DEM simulation. Albeit some sophisticated algorithms have been employed for accomplishing this object, e.g. the improved simulated annealing algorithm [37], a sequential quasi-Monte Carlo filter [38] and the non-dominated sorting genetic algorithm [39], the backpropagation (BP) neural network is still the most commonly used algorithm in the previous attempts [40-43] due to its simplicity and widely available software. In this study, the BP neural network is used to establish the relationship between the shear strength and micro parameters for MatDEM modeling.
270
271
BP neural network is a multilayer feed-forward neural network which consists of an input layer, implicit layer and output layer as shown in Fig.2.
272
273
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
274
 [[Image:Review_852894090255-image16.png|354px]] </div>
275
276
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
277
<span style="text-align: center; font-size: 75%;">Fig.2 Typical structure of BP neural network</span></div>
278
279
When the value of neurons in the input layer is given, neurons in the implicit layer and the output layer are calculated according to the following forward propagation:
280
281
{| class="formulaSCP" style="width: 100%; text-align: center;" 
282
|-
283
| 
284
{| style="text-align: center; margin:auto;" 
285
|-
286
| <math>u_i^l={\sigma }_i^l\left(\sum_{j=1}^{k_{l-1}}w_{ji}^lu_j^{l-1}+\right. </math><math>\left. b_i^l\right)</math>
287
|}
288
| style="width: 5px;text-align: right;white-space: nowrap;" | (14)
289
|}
290
291
292
in which ''l'' is the number of layers in the implicit layer, ''ul i'' denotes the ''i''-th neuron in the ''l''-th layer, ''σl i'' is the activation function, ''k<sub>l</sub>'' is the number of neurons in the ''l''-th layer, ''wl ji'' is the weight and ''bl i''  is the deviation. The training of BP neural network indicates constantly adjusting the weights and deviations to minimize the differences between the output value and the true value.
293
294
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
295
 [[Image:Review_852894090255-image18.png|396px]] </div>
296
297
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
298
<span style="text-align: center; font-size: 75%;">Fig.3 BP neural network prediction model for the shear strength</span></div>
299
300
The structure of the BP neural network prediction model for the shear strength is plotted in Fig. 3. The micro parameters are placed in the input layer and the shear strengths are placed in the output layer. It is noticed that only the micro parameters in terms of the particle contact model is involved in the BP neural network. Although some other micro parameters, e.g. the particle radius and the particle density, have been verified to influence the macro mechanical properties of geotechnical materials in the numerical model [39], their effect is ignored in this study. The reason is that the particle radius and the particle density will be constant when the numerical model of the slope has been established in MatDEM.
301
302
Considering that an ultimate aim of this study is to acquire the micro parameters corresponding to the original and reduced shear strength, the established BP neural network seems unfavorable. The reasons are explained as follows. If placing the micro parameters in the output layer and the shear strengths in the input layer, the output layer have more neurons than the input layer, which forms an underdetermined mathematical problem. In this situation, the BP neural network is possible to result in a result of low precision. For example, to calibrate micro particle parameters for Particle Flow Code (PFC) modeling, a BP neural network was built by placing three macro geotechnical mechanical parameters in the input layer and four micro parameters in the output layer [43]. Although the BP neural network database consists of four hundred samples, results of some testing samples were still unsatisfactory. Thus, the micro parameters and the shear strengths are placed in the input layer and the output layer respectively, to promise the validity of the BP neural network. Base on the BP neural network predict model of the shear strength, a calibration of the micro particle parameters will be developed similar to the work in [42].
303
304
==3.2 Virtual 3D direct shear tests to measure the shear strength==
305
306
The direct shear test and the triaxial compression test are two commonly used tests to measure the shear strength of geotechnical materials in laboratory, and their virtual counterparts have been established in DEM simulation by researchers [38, 44]. Under the same micro parameter conditions, the two virtual tests are possible to result in different values of the shear strength [44]. Considering that the shear strength of geotechnical materials in the latter case study is obtained by using direct shear test, virtual 3D direct shear test is used to measure the shear strength in this study.
307
308
As Fig.4 shows, the virtual 3D direct shear test model consists of two parts, the specimen and the shear box. The specimen is composed of 7636 spherical particles, whose radiuses vary from 0.83 mm to 1.20 mm. The radius of the specimen is 30 mm and the height is 20 mm. The shear box consists of 7374 walls. To simulate the shearing process, the imposed boundary conditions are the horizontal velocity of the lower box and the normal pressure on the top wall. By prescribing the horizontal displacement of the lower box in each time step to be 0.05 mm, the shear stress was recorded until the specimen was finally damaged.
309
310
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
311
 [[Image:Review_852894090255-image19.png|360px]] </div>
312
313
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
314
<span style="text-align: center; font-size: 75%;">Fig.4 Virtual 3D direct shear test model for measuring the shear strength of particle aggregate</span></div>
315
316
Given the micro particle parameters as ''K<sub>n</sub>''=5.0&#x00d7;10<sup>6</sup> N/mm, ''K<sub>s</sub>''=1.0&#x00d7;10<sup>6</sup> N/mm, ''X<sub>b</sub>''=1.0&#x00d7;10<sup>-5</sup> m, ''F<sub>s0</sub>''=1.0&#x00d7;10<sup>6</sup> Pa and'' u<sub>p</sub>''=0.20, four normal pressures, 100 kPa, 200kPa, 300kPa and 400 kPa, are applied on the top wall to perform the virtual 3D direct shear test. The shear stress−shear displacement curves under the four normal pressures are plotted in Fig. 5a. The trend of the four curves is consistent with each other. The increasing of the imposed normal pressures induces a higher peak value of the shear stress, and the shear displacement where the shear stress reaches its peak value increases slightly, which is similar to the result of the previous virtual direct shear test [44]. The Coulomb formula
317
318
{| class="formulaSCP" style="width: 100%; text-align: center;" 
319
|-
320
| 
321
{| style="text-align: center; margin:auto;" 
322
|-
323
| <math>\tau \mbox{=}C\mbox{+}\sigma tan\varphi </math>
324
|}
325
| style="width: 5px;text-align: right;white-space: nowrap;" | (15)
326
|}
327
328
329
is employed to determine the shear strength of the particle aggregate as shown in Fig. 5b. Finally, the friction angle ''φ'' is 32.4°, and the cohesion ''C'' is 76 kPa.<span style="text-align: center; font-size: 75%;">'''  '''</span>
330
331
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
332
<span style="text-align: center; font-size: 75%;">'''  [[Image:Review_852894090255-image21.png|270px]] '''</span></div>
333
334
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
335
<span style="text-align: center; font-size: 75%;">(a)The shear stress−shear displacement curves</span></div>
336
337
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
338
<span style="text-align: center; font-size: 75%;">''' [[Image:Review_852894090255-image22.png|246px]] '''</span></div>
339
340
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
341
<span style="text-align: center; font-size: 75%;">          (b) Determination of the shear strengths</span></div>
342
343
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
344
<span style="text-align: center; font-size: 75%;">Fig.5 Virtual 3D direct shear test model for the shear strength measurement</span></div>
345
346
==3.3 A correction calibration of micro particle parameters based on the BP neural network ==
347
348
Based on the previous studies on the influence of the micro particle parameters in MatDEM on the mechanical properties of the particle aggregate [32], four levels are considered for ''F<sub>s0</sub>'' and ''&#x03bc;<sub>p</sub>'', and three levels are considered for ''K<sub>n</sub>'','' K<sub>s</sub>'' and ''X''<sub>b</sub> when establishing the BP neural network database. Table 1 lists the level values of micro particle parameters. By using the full orthogonal combination, 432 combinations of micro particle parameters are generated. And then, the virtual 3D direct shear test is executed to obtain their corresponding shear strengths. Finally, the BP neural network database is composed by 432 groups of micro particle parameters and the corresponding shear strengths.
349
350
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
351
<span style="text-align: center; font-size: 75%;">'''Table 1 The level values of micro particle parameters to establish the neural network database'''</span></div>
352
353
{| style="width: 100%;border-collapse: collapse;" 
354
|-
355
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''K''<sub>s</sub></span><span style="text-align: center; font-size: 75%;">(N/mm)</span>
356
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''K''<sub>n</sub></span><span style="text-align: center; font-size: 75%;">(N/mm)</span>
357
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''F<sub>s0</sub>''</span><span style="text-align: center; font-size: 75%;">(P<sub>a</sub></span><span style="text-align: center; font-size: 75%;"></span>
358
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''u<sub>p</sub>''</span>
359
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''X''<sub>b</sub></span><span style="text-align: center; font-size: 75%;">(m)</span>
360
|-
361
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.0&#x00d7;10<sup>5</sup></span>
362
363
<span style="text-align: center; font-size: 75%;">1.0&#x00d7;10<sup>6</sup></span>
364
365
<span style="text-align: center; font-size: 75%;">1.0&#x00d7;10<sup>7</sup></span>
366
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.0&#x00d7;10<sup>5</sup></span>
367
368
<span style="text-align: center; font-size: 75%;">5.0&#x00d7;10<sup>6</sup></span>
369
370
<span style="text-align: center; font-size: 75%;">5.0&#x00d7;10<sup>7</sup></span>
371
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.0&#x00d7;10<sup>4</sup></span>
372
373
<span style="text-align: center; font-size: 75%;">1.0&#x00d7;10<sup>5</sup></span>
374
375
<span style="text-align: center; font-size: 75%;">5.0&#x00d7;10<sup>5</sup></span>
376
377
<span style="text-align: center; font-size: 75%;">1.0&#x00d7;10<sup>6</sup></span>
378
|  style="border: 1pt solid black;text-align: center;vertical-align: top;"|<span style="text-align: center; font-size: 75%;">0.10</span>
379
380
<span style="text-align: center; font-size: 75%;">0.20</span>
381
382
<span style="text-align: center; font-size: 75%;">0.30</span>
383
384
<span style="text-align: center; font-size: 75%;">0.40</span>
385
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.0&#x00d7;10<sup>-6</sup></span>
386
387
<span style="text-align: center; font-size: 75%;">1.0&#x00d7;10<sup>-5</sup></span>
388
389
<span style="text-align: center; font-size: 75%;">1.0&#x00d7;10<sup>-4</sup></span>
390
|}
391
392
393
The training of the BP neural network actually successfully establishes the mapping relationship between the micro particle parameters in the input layer and the shear strengths in the output layer, which can be denoted as the following formula.
394
395
{| class="formulaSCP" style="width: 100%; text-align: center;" 
396
|-
397
| 
398
{| style="text-align: center; margin:auto;" 
399
|-
400
| <math>\left(C,tan\varphi \right)=f\left(K_n,K_s,X_b,F_{s0},{\mu }_p\right)</math>
401
|}
402
| style="width: 5px;text-align: right;white-space: nowrap;" | (16)
403
|}
404
405
406
Assuming that the shear strengths (outputs) are continuous in the micro particle parameter space, the micro particle parameters (inputs) corresponding to the target shear strengths could be resolved similarly to the resolution of the nonlinear equations, because the gradient at any given input could be computed in the BP neural network. In the first place, an objective function is defined as
407
408
{| class="formulaSCP" style="width: 100%; text-align: center;" 
409
|-
410
| 
411
{| style="text-align: center; margin:auto;" 
412
|-
413
| <math>Error={\left(\frac{C-C^\ast }{C^\ast }\right)}^2+{\left(\frac{tan\varphi -tan{\varphi }^\ast }{tan{\varphi }^\ast }\right)}^2</math>
414
|}
415
| style="width: 5px;text-align: right;white-space: nowrap;" | (17)
416
|}
417
418
419
where ''C''<sup>*</sup> and ''&#x03c6;''<sup>*</sup>'' ''are the target cohesion and the target friction angle, respectively. And then, the gradient descent method is used to adjust the inputs until the objective function have a value smaller than the prescribed threshold value. For a simplicity expression, the input (''K<sub>n</sub>'','' K<sub>s</sub>'', ''X''<sub>b</sub>, ''F<sub>s0</sub>'','' &#x03bc;<sub>p</sub>'') is denoted by '''X''', the output (''C'','' ''tan''&#x03c6;'') is denoted by '''Y''', and thus Eq. 16 can be written as '''Y '''= ''f''('''X'''). The resolution process is explained as follows:
420
421
(1) Determine the initial exploration input '''X'''<sub>0</sub> by using the K-means clustering analysis [45]. K-means clustering is a vector quantization method originally from signal processing and can partition ''n'' observations into ''k'' clusters where each observation belongs to the cluster with the nearest mean. The samples in the neural network database are partitioned into ''k'' clusters by using K-means clustering analysis. Distances between the cluster centers and the target are computed by Eq. 17, and the nearest cluster center is selected as '''X'''<sub>0</sub>.
422
423
(2) Calculate the response outputs '''Y''' ''<sub>i</sub>'' for the exploration point of '''X''' ''<sub>i</sub>'' at step ''i'', and then invoke the objective function. When the objective function has a value smaller than the threshold value, the resolution is terminated and '''X''' takes the value of '''X''' ''<sub>i</sub>''. If not, continue to perform step (3).
424
425
(3) Calculate the gradient of the objective function at '''X''' ''<sub>i</sub>'', and update '''X''' along the descent direction of the objective function. The formulation to update '''X''' is
426
427
{| class="formulaSCP" style="width: 100%; text-align: center;" 
428
|-
429
| 
430
{| style="text-align: center; margin:auto;" 
431
|-
432
| <math>\Delta \boldsymbol{X}=-\eta \frac{\partial \left(Error\right)}{\partial \boldsymbol{X}}\vert _{\boldsymbol{X}=\boldsymbol{X}_i},\boldsymbol{X}_{i+1}=</math><math>\boldsymbol{X}_i+\Delta \boldsymbol{X}</math>
433
|}
434
| style="width: 5px;text-align: right;white-space: nowrap;" | (18)
435
|}
436
437
438
in which Δ'''X''' is the increment used to update '''X''' ''<sub>i</sub>'' and ''η'' is the step length. Δ'''X''' is obtained by multiplying the negative gradient of the objective function by ''η''. On account of the nonlinearity of forward propagation rule, ''η ''is introduced for performing a step search approximation along the error gradient descent direction.
439
440
(4) Let ''i'' =''i''+1 and return to step (2).
441
442
Forward propagation of the established BP neural network can be regarded as a substitution of the virtual 3D direct shear test. The samples in the neural network database are always finite, which indicates that the substitution is certain to be approximate. Therefore, the inputs obtained from the mathematical resolution may have a low accuracy. An additional step called as “verifying test” is used to determine whether the inputs are acceptable. The verifying test performs virtual 3D direct shear test by using the inputs obtained from the mathematical resolution, and calculates the error between the resulted shear strength and the target by using Eq. 17. If the objective function has a value smaller than the threshold value, the inputs are acceptable. Otherwise, a correction should be introduced for obtain a more reasonable input.
443
444
The failure of the verifying test indicates that the current BP neural network is still not accurate enough to obtain a fine mathematical resolution. A feasible strategy to improve the precision of BP neural network is adding some new samples to the database. Although the current BP neural network is not accurate enough, it provides a rough search direction for the resolution and the inputs resulted by the last round of the mathematical resolution must be a bit closer to the optimal solution than the initial inputs. Considering that, a new sample consisting of the resulted inputs in the last round and the corresponding shear strength is added to the database, and then the neural network is retrained. After the retraining of the BP neural network, the mapping relationship in Eq. 16 has been replaced by a new relationship. Then, the mathematical resolution is executed again by selecting the resulted inputs in the last round as the initial exploration input, and the verifying test is performed again. The “mathematical resolution−verifying test−neural network correction” flow is repeated till the objective function has a value smaller than the threshold value in verifying test. The entire flow of the correction method to calibrate the micro particle parameters is plotted in Fig. 6.
445
446
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
447
 [[Image:Review_852894090255-image26.png|600px]] </div>
448
449
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
450
<span style="text-align: center; font-size: 75%;">Fig.6 The entire flow of the correction method to calibrate the micro particle parameters</span></div>
451
452
==3.5 Validity of the proposed correction calibration ==
453
454
To verify the validity of the proposed calibration method, the target shear strengths are given as ''C'' =30.0 kPa and ''&#x03c6; ''=19.6<sup>o</sup>, and the optimal micro particle parameters are searched by prescribing the threshold value as 2.0×10<sup>−4</sup> for Eq. 17. After four rounds of mathematical resolution, the resulted micro particle parameters passed the verifying test. Table 2 lists the micro parameters provided by the mathematical resolution and the shear strengths provided by the verifying test in each round.
455
456
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
457
<span style="text-align: center; font-size: 75%;">'''Table 2 Micro particle parameters and shear strengths in four rounds'''</span></div>
458
459
{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
460
|-
461
|  colspan='2'  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Parameter</span>
462
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''K<sub>n</sub>''</span><span style="text-align: center; font-size: 75%;">(N/m)</span>
463
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''K<sub>s</sub>''</span><span style="text-align: center; font-size: 75%;">(N/m)</span>
464
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''X<sub>b</sub>''</span><span style="text-align: center; font-size: 75%;">(m)</span>
465
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''F<sub>s0</sub>''</span><span style="text-align: center; font-size: 75%;">(Pa)</span>
466
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''&#x03bc;<sub>p</sub>''</span>
467
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''C''(kPa)</span>
468
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''&#x03c6;''(&#x00b0;)</span>
469
|-
470
|  rowspan='4' style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Round</span>
471
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1</span>
472
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.27&#x00d7;10<sup>7</sup></span>
473
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.27&#x00d7;10<sup>6</sup></span>
474
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.95&#x00d7;10<sup>-5</sup></span>
475
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.25&#x00d7;10<sup>5</sup></span>
476
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.320</span>
477
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">21.7</span>
478
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">30.2</span>
479
|-
480
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2</span>
481
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.47&#x00d7;10<sup>7</sup></span>
482
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.84&#x00d7;10<sup>6</sup></span>
483
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.72&#x00d7;10<sup>-5</sup></span>
484
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.67&#x00d7;10<sup>5</sup></span>
485
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.220</span>
486
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">28.4</span>
487
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">22.8</span>
488
|-
489
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3</span>
490
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.03&#x00d7;10<sup>7</sup></span>
491
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.32&#x00d7;10<sup>6</sup></span>
492
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.84&#x00d7;10<sup>-5</sup></span>
493
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">6.42&#x00d7;10<sup>5</sup></span>
494
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.143</span>
495
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">29.5</span>
496
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">18.9</span>
497
|-
498
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">4</span>
499
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.27&#x00d7;10<sup>7</sup></span>
500
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.09&#x00d7;10<sup>6</sup></span>
501
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.04&#x00d7;10<sup>-5</sup></span>
502
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.95&#x00d7;10<sup>5</sup></span>
503
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.141</span>
504
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">30.2</span>
505
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">19.6</span>
506
|}
507
508
509
Taking the target shear strengths as the reference, the relative errors of the resulted shear strengths are 27.7% and 63.5% for the cohesion and the friction coefficient in the first round, respectively. In the second round, they decrease to 5.3% and 18.1% for the cohesion and the friction coefficient, respectively. Ultimately, they decrease to 0.7% and 0.2% in the fourth round for the cohesion and the friction coefficient, respectively. Thus, the proposed method has the ability to obtain the micro particle parameters that well reflects the shear strengths of particle assembly.
510
511
When performing the slope stability analysis by MatDEM coupled with SRM, the shear strengths of the geological material in the MatDEM modeling are constantly adjusted. Once a new value of the shear strengths is given, a new combination of micro particle parameters is to be searched. At the moment, the established neural network database will provide a new initial exploration input '''X'''<sub>0</sub> by using the K-means clustering analysis, and then the proposed method is performed to obtain the micro particle parameters.
512
513
=4 An energy criterion to evaluate the slope stability=
514
515
As mentioned in Section 2.2, the MatDEM modeling follows the law of energy conversation and Newton’s law of motion. When slope instability happens, the displacements of soil and rock will induce the conversion between various kinds of energies. Consequently, the energy conversion of an instable slope in MatDEM must be different with that of a stable slope. In this section, a simple 2D slope is taken as an example to investigate the energy conversion difference between the stable and instable slope. And then, an energy criterion is established to evaluate the slope stability.
516
517
==4.1 Investigation on the energy conversion difference between the stable and instable slope==
518
519
As Fig. 7 shows, the slope is 7 m height and the ratio of slope is 7: 3. And, the slope is assumed to be homogeneous and composed by one kind of geotechnical material for simplicity. The MatDEM model of the slope contains 6305 particles. The radius of particles varies from 0.064 m to 0.096 m. The density is taken as 2.04&#x00d7;10<sup>3</sup> kg/m<sup>3</sup>, and the gravity acceleration ''g'' = 9.8 m/s<sup>2</sup>. All particles at the left, right and bottom boundaries are fixed as the boundary conditions of the model.
520
521
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
522
 [[Image:Review_852894090255-image27.png|300px]] </div>
523
524
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
525
<span style="text-align: center; font-size: 75%;">Fig. 7 A simple homogeneous slope 2D model</span></div>
526
527
In the numerical simulation of the slope model, three groups of micro particle parameters listed in Table 3 are selected, which are denoted as Case 1, 2 and 3 respectively. Considering that the chief purpose here is to reveal the energy conversion difference between the stable and instable slope, the shear strengths are not measured for the three groups of micro particle parameters by using the virtual numerical test.
528
529
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
530
<span style="text-align: center; font-size: 75%;">'''Table 3 Micro particle parameters to be used in the investigation of energy conversion'''</span></div>
531
532
{| style="width: 83%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
533
|-
534
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Parameter</span>
535
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''K<sub>n</sub>''</span><span style="text-align: center; font-size: 75%;">(N/m)</span>
536
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''K<sub>s</sub>''</span><span style="text-align: center; font-size: 75%;">(N/m)</span>
537
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''X<sub>b</sub>''</span><span style="text-align: center; font-size: 75%;">(m)</span>
538
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''F<sub>s0</sub>''</span><span style="text-align: center; font-size: 75%;">(Pa)</span>
539
|  style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''&#x03bc;<sub>p</sub>''</span>
540
|-
541
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Case 1</span>
542
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">7.80&#x00d7;10<sup>6</sup></span>
543
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.10&#x00d7;10<sup>6</sup></span>
544
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.10&#x00d7;10<sup>-4</sup></span>
545
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.46&#x00d7;10<sup>5</sup></span>
546
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.020</span>
547
|-
548
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">Case 2</span>
549
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.80&#x00d7;10<sup>8</sup></span>
550
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.50&#x00d7;10<sup>8</sup></span>
551
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.60&#x00d7;10<sup>-</sup></span><span style="text-align: center; font-size: 75%;"><sup>5</sup></span>
552
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4.90&#x00d7;10<sup>5</sup></span>
553
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.120</span>
554
|-
555
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Case 3</span>
556
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">9.60&#x00d7;10<sup>7</sup></span>
557
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.80&#x00d7;10<sup>7</sup></span>
558
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">9.30&#x00d7;10<sup>-6</sup></span>
559
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.90&#x00d7;10<sup>5</sup></span>
560
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.309</span>
561
|}
562
563
564
After 10 time steps of the computation, the slope is verified to be instable in Case 1 but stable in Case 2 and 3. The displacement cloud charts are plotted in Fig. 8a and Fig. 9a for Case 1 and Case 2, respectively. The displacement cloud chart of Case 3 is not provided because the displacement difference between Case 2 and Case 3 is not noticeable. Fig. 8b and Fig. 9b illustrate the variations of the total energy, elastic potential energy, kinetic energy, gravitational potential energy and heat with the time steps for Case 1 and Case 2, respectively. Some interesting phenomena are observed in the energy-time steps curves, and a further investigation is valuable for the establishment of an energy criterion to evaluate the slope stability.
565
566
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
567
 [[Image:Review_852894090255-image28-c.png|324px]] </div>
568
569
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
570
<span style="text-align: center; font-size: 75%;">(a) Displacement cloud chart</span></div>
571
572
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
573
<span style="text-align: center; font-size: 75%;"> [[Image:Review_852894090255-image29.png|312px]] </span></div>
574
575
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
576
<span style="text-align: center; font-size: 75%;">(b) Energy variations with time step</span></div>
577
578
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
579
<span style="text-align: center; font-size: 75%;">Fig.8 Displacements and the energy conversion for the slope in Case 1</span></div>
580
581
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
582
 [[Image:Review_852894090255-image30-c.png|330px]] </div>
583
584
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
585
<span style="text-align: center; font-size: 75%;">(a) Displacement cloud chart                  </span></div>
586
587
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
588
<span style="text-align: center; font-size: 75%;"> [[Image:Review_852894090255-image31.png|306px]] </span></div>
589
590
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
591
<span style="text-align: center; font-size: 75%;">(b) Energy variations with time step</span></div>
592
593
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
594
<span style="text-align: center; font-size: 75%;">Fig.9 Displacements and the energy conversion for the slope in Case 2</span></div>
595
596
Firstly, the conclusion that MatDEM follows the law of energy conversation is confirmed again by the energy-time steps curves. For both two cases, the total energies are almost constant during the computation, no matter whether the slope is stable or not. However, the total energy of Case 2 has a larger value than that of Case 1, which can be explained as follows. In the MatDEM modeling, the gravitational potential energy is evaluated by the relative vertical displacement of particles. It is initialized to be zero for a particle, and becomes negative when the particle falls while positive when the particle rises. The kinetic energy is zero for a particle unless the motion of the particle happens. The heat induced by the viscous damping, the spring breaking and the friction is zero in the initial state. Therefore, the total energy should be equal to the elastic potential energy initially. Since the elastic potential energy is computed as Eq. 4, ''K''<sub>n</sub> and ''K''<sub>s</sub> dominate the value of the elastic potential energy. Because the values of ''K''<sub>n</sub> and ''K''<sub>s</sub> in Case 2 is greater than that in Case 1, the total energy of Case 2 has a larger value than that of Case 1.
597
598
Secondly, the variations of the kinetic energy and the gravitational potential energy in Fig. 8b and Fig. 9b are apparently different. With the time step increasing, the kinetic energy increases but the gravitational potential energy decreases in Fig. 8b. But, the gravitational potential energy and the kinetic energy are nearly constant in Fig. 9b. Fig. 10 plots the variations of the kinetic energy and the gravitational potential energy with time steps for all three cases.
599
600
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
601
 [[Image:Review_852894090255-image32.png|306px]] </div>
602
603
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
604
<span style="text-align: center; font-size: 75%;">(a) The kinetic energy                  </span></div>
605
606
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
607
<span style="text-align: center; font-size: 75%;"> [[Image:Review_852894090255-image33.png|306px]] </span></div>
608
609
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
610
<span style="text-align: center; font-size: 75%;">(b) The gravitational potential energy</span></div>
611
612
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
613
<span style="text-align: center; font-size: 75%;">Fig.10 Variations of the kinetic energy and the gravitational potential energy with time steps for all three cases</span></div>
614
615
For Case 3, the kinetic energy is observed, and the gravitational potential energy is positive. That may be caused by the slope deformation. The kinetic energy and the gravitational potential energy caused by the deformation can be easily distinguished from that caused by the failure. Above all, the magnitude of the energy induced by the deformation is tiny compared to that induced by the failure. For Case 3, the kinetic energy and the gravitational potential energy have peak values of 2800 J and 4700 J, respectively. For Case 2, they have peak values of 116 J and 945 J, respectively. Next, when the two energies are caused by deformation, their variations show the convergent trend. But, when the two energies are induced by the failure, their variations show the divergent trend. Since the clear difference in the variation of the kinetic energy and the gravitational potential energy can be observed between a stable slope and an instable slope, they have the potential to be used in the evaluation of the slope stability.
616
617
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
618
 [[Image:Review_852894090255-image34.png|312px]] </div>
619
620
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
621
<span style="text-align: center; font-size: 75%;">(a) Three cases                           </span></div>
622
623
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
624
<span style="text-align: center; font-size: 75%;"> [[Image:Review_852894090255-image35.png|300px]] </span></div>
625
626
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
627
<span style="text-align: center; font-size: 75%;">(b) Case 1</span></div>
628
629
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
630
<span style="text-align: center; font-size: 75%;">Fig.11 Variations of the heat with time steps</span></div>
631
632
Finally, the heat is investigated. The variations of heat are illegible in Fig. 8b and Fig. 9b because the magnitude of heat is too small. Fig. 11a plots the variations of the heat with time steps for all three cases. Obviously, the heats in Case 2 and Case 3 are much greater than that in Case 1. But it is incorrect to have a conclusion that the heat in a stable MatDEM slope model is more than that in an instable MatDEM slope model. The heat in the system is composed of the contributions of the viscous damping, the spring breaking and the friction force, in which the latter two sources are calculated based on the micro particle parameters. In Fig. 11a, the heat curve for Case 1 is situated in the lowest position, which maybe resulted by the micro particle parameters in Table 3. Thus, the magnitude of the heat is unreasonable to be chosen as a reference in the slope stability assessment. But, the shapes of the three heat curves are quite different. For a better demonstration, the heat curve is repainted for Case 1 in Fig. 11b. The heat curves of Case 2 and Case 3 are roughly convex, but the heat curve of Case 1 is concave. If the slope is stable, the position change of the particles induced by the deformation will stop rapidly, which promises the heat to be stable after several time steps. If the slope is sliding, the positions of the particles are still continuously changed, which results in the on-going increasing of the heat. Therefore, the shape of the heat curve may be effective in the evaluation of the slope stability.
633
634
==4.2 A criterion based on the kinetic energy of the MatDEM model==
635
636
As mentioned in the Section 4.1, the variation of the kinetic energy and the gravitational potential energy and the shape of the heat curve are verified to be different between a stable slope and an instable slope when performing the simulation by MatDEM. Their validities are compared in this section from a view of practical application.
637
638
When performing the slope stability analysis by using MatDEM with SRM, the micro parameters of the particles will be modified according to the reduction factor. Both the gravitational potential energy and the kinetic energy are independent of the micro particle parameters. However, the heat is calculated based on the micro particle parameters. When the energy variations resulted by using different micro particle parameters are collected together, it is more suitable to compare the results of the gravitational potential energy and the kinetic energy than that of the heat. Additionally, the number of the time steps has a great influence on the reliability of the conclusion if estimating the stable state of a slope by using the shape of the heat curves. As Fig. 11a shows, the overall trend of the heat curve is convergent for Case 2, while a slight turning happens at the sixth time step. A large number of time steps seem to be feasible to overcome the effect of the local abnormality, but it will cause another issue. In fact, the heat curve must be convergent ultimately no matter whether the slope is stable or not. In case of an instable slope, the position of the particles will be relocated finally when the sliding is terminated, and thus the heat will be also stable. Therefore, an incorrect conclusion may be drawn based on the shape of the heat curve when the number of time steps is inappropriate. In short, the gravitational potential energy and the kinetic energy are more suitable to be compared than the heat when estimating the stable state of a slope.
639
640
There is an essential association between the gravitational potential energy and the kinetic energy. The motion of the particles results in the change of their positions, so the kinetic energy always emerges before the gravitational potential energy in MatDEM. This deduction can be verified by the variations of the kinetic energy and the gravitational potential energy for Case 3. As shown in Fig. 10, for Case 3 it is the fourth time step that the kinetic energy begins to decrease, while it is the seventh time step that the gravitational potential energy begins to decrease. This phenomenon indicates that the convergence of the kinetic energy should be earlier than that of the gravitational potential energy for a stable slope. Moreover, the variation of the gravitational potential energy is noticed to be small when the time step ranges from four to seven for Case 1. The reason is that the motion of the particles in the model maybe temporarily dominated by the horizontal motion during the failure, but only the vertical motion of the particles induces the variation of the gravitational potential energy. Thus, the gravitational potential energy may have a more complicated curve than the kinetic energy. Based on the above analysis, the kinetic energy is taken as the main reference to estimate the stable state of a slope.
641
642
When performing the slope stability analysis by using MatDEM coupled with SRM, the variations of the kinetic energy are collected under different reduction factors. Then, the variation tendencies and the magnitude of the kinetic energies are compared. The abrupt in the variation tendencies and the magnitude is considered as a sign for the emergence of the critical point, and the safe factor of the slope is determined finally. At this point, a criterion has been proposed to evaluate the slope stability based on the kinetic energy.
643
644
=5 Case Study: Baijiabao landslide=
645
646
==5.1 Geological Background==
647
648
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
649
 [[Image:Review_852894090255-image36.png|600px]] </div>
650
651
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
652
<span style="text-align: center; font-size: 75%;">Fig.12 A photograph of Baijiabao landslide, taken from the opposite bank of the Xiangxi River</span></div>
653
654
Baijiabao landslide is situated in Guizhou village of Zigui County, Hubei Province of China on the right bank of the Xiangxi River (30<sup>o</sup>58′59.9′′N, 110<sup>o</sup>45′33.4′′E). A photograph of the landslide that was taken from the opposite bank of the Xiangxi River is shown in Fig. 12. The elevation of the landslide ranges from 125 m to 265 m above MSL, and the slope of the landslide ranges from 10<sup>o</sup> to 20<sup>o</sup>. The upper boundary of the landslide is defined by the interface between the bedrock and the soil. The left and right boundaries are defined by two natural gullies that contain many surface cracks. The toe of the landslide varies in elevation between 125 m to 135 m. The primary sliding direction of Baijiabao landslide varies between 75<sup>o</sup> and 85<sup>o</sup> (SW-NE). This landslide is 550 m long and 400 m wide, covers an area of 20 ha, and has an estimated volume of nearly 1 million m<sup>3</sup>. The landslide mass is primarily composed of loose Quaternary deposits. Its sliding zone is mainly silty clay. And the underlying bedrock contains quartz sandstone and argillaceous siltstone of the Early Jurassic Xiangxi Formation, which dips into the hill at an angle of 30<sup>o</sup> to 40<sup>o</sup> (Fig. 13).
655
656
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
657
<span style="text-align: center; font-size: 75%;"> [[Image:Review_852894090255-image37.png|414px]] </span></div>
658
659
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
660
<span style="text-align: center; font-size: 75%;">(a) Engineering geological map                           </span></div>
661
662
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
663
<span style="text-align: center; font-size: 75%;"> [[Image:Review_852894090255-image38.png|414px]] </span></div>
664
665
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
666
<span style="text-align: center; font-size: 75%;">(b) Geological cross-section</span></div>
667
668
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
669
<span style="text-align: center; font-size: 75%;">Fig.13 Geological background of Baijiabao landslide [46]</span></div>
670
671
Baijiabao landslide poses significant threats to public safety. The Zi-Xing road crosses the central part of the landslide mass as Fig. 13a shows. Before Three Gorges Reservoir (TGR) impoundment, 165 residents lived in the landslide area, while only 20 residents live there today. The deformation rate of the landslide was observed to increase significantly after TGR impoundment. Monitoring system installed in October 2006 consists of four GPS monitoring sites, which are noted as ZG323, ZG324, ZG325 and ZG326 in Fig. 13a. The detailed deformation history of Baijiabao landslide can be referred to the related literatures [46-48].
672
673
==5.2 Landslide model in the MatDEM modeling==
674
675
This subsection describes the establishment of the 3D MatDEM model for the stability analysis of Baijiabao landslide. Firstly, the geomorphy of the landslide is simulated by using ArcGIS software based on the contour line measured in 2006 (Fig. 14a). Then, a 3D model of Baijiabao landslide is built in MatDEM. The model is 564 m long, 497 m wide and 186 m height, covering the whole landslide mass (Fig. 14b). Next, geotechnical materials are allocated according to the engineering geologic investigation data. The sliding zone in the model (Fig. 14c) is positioned by using the interpolation of locations where the sliding zone is encountered in the drilling exploration. The landslide mass has a thickness varying between 20 m to 30 m in the front area, and has a thickness varying between 10 m to 40 m in the trailing area.
676
677
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
678
<span style="text-align: center; font-size: 75%;"> [[Image:Review_852894090255-image39.png|198px]] </span></div>
679
680
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
681
<span style="text-align: center; font-size: 75%;">(a)</span> <span style="text-align: center; font-size: 75%;">The geomorphy model        </span></div>
682
683
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
684
<span style="text-align: center; font-size: 75%;"> [[Image:Review_852894090255-image40.png|360px]] </span></div>
685
686
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
687
<span style="text-align: center; font-size: 75%;">(b) 3D MatDEM model</span></div>
688
689
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
690
<span style="text-align: center; font-size: 75%;"> [[Image:Review_852894090255-image41.png|318px]] </span></div>
691
692
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
693
<span style="text-align: center; font-size: 75%;">(c) An enlargement to exhibit the sliding zone</span></div>
694
695
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
696
<span style="text-align: center; font-size: 75%;">Fig.14 3D MatDEM model for the stability analysis of Baijiabao landslide</span></div>
697
698
Finally, the 3D MatDEM model of Baijiabao landslide is composed by 1,207,985 particles, in which the landslide mass and the sliding zone are represented by 733,104 removable particles, and the bedrock is represented by 474,881 fixed particles. The radius of particles varies between 0.55 m and 0.83 m. In Fig. 14c, the landslide mass, the sliding zone and the bedrock are presented as brown particles, green particles and blue particles, respectively.
699
700
==5.3 Mechanical parameters and micro particle parameters==
701
702
According to the report of engineering geological exploration, the unit weight and shear strengths are listed in Table 4 for the geotechnical materials in the landslide mass, the sliding zone and the bedrock, respectively. All parameters are measured in the native state, and taken their statistical mean value. It is noticed that the stability of Baijiabao landslide is actually affected by the rainfall and the reservoir water fluctuation. In view that the main purpose of this study is to improve the ability of MatDEM in 3D slope stability evaluation, the rainfall and the reservoir water fluctuation are ignored temporarily.
703
704
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
705
<span style="text-align: center; font-size: 75%;">'''Table 4 Unit weight and shear strengths of geotechincal materials in different regions'''</span></div>
706
707
{| style="width: 100%;border-collapse: collapse;" 
708
|-
709
|  style="border-top: 2pt solid black;border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Region</span>
710
|  style="border-top: 2pt solid black;border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Unit weight ''γ ''(kN/m<sup>3</sup></span><span style="text-align: center; font-size: 75%;">)</span>
711
|  style="border-top: 2pt solid black;border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''C'' (kPa)</span>
712
|  style="border-top: 2pt solid black;border-left: 1pt solid black;border-bottom: 1pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''&#x03c6;'' (<sup>o</sup></span><span style="text-align: center; font-size: 75%;">)</span>
713
|-
714
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Landslide mass</span>
715
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">18.0</span>
716
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">17.1</span>
717
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">20.2</span>
718
|-
719
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Sliding zone</span>
720
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">20.5</span>
721
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">24.6</span>
722
|  style="border: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">15.2</span>
723
|-
724
|  style="border-top: 1pt solid black;border-left: 1pt solid black;border-bottom: 2pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Bedrock</span>
725
|  style="border-top: 1pt solid black;border-left: 1pt solid black;border-bottom: 2pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">26.0</span>
726
|  style="border-top: 1pt solid black;border-left: 1pt solid black;border-bottom: 2pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1,200</span>
727
|  style="border-top: 1pt solid black;border-left: 1pt solid black;border-bottom: 2pt solid black;border-right: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">32.1</span>
728
|}
729
730
731
Using the BP neural network-based calibration method established in Section 3, micro parameters of particles representing various geotechnical materials are obtained and listed in Table 5. Because the particles have a much larger radius than the particles used in Section 3, the size of the virtual 3D direct shear test model is adjusted. Previous studies [21, 49] verified the effect of the particle size on the macro properties in DEM modeling, and declared that the ratio ''L''/''R'' is a critical variable to be carefully selected when performing virtual strength test. Here ''L'' is the size of the strength test model, and ''R'' denotes the average radius of particles. According to Su’s suggestion [49], ''L''/''R'' should be more than 200 to avoid the effect of the particle size on the macro properties of the particle assembly in DEM. Finally, the specimen to be used in virtual 3D direct shear test has a radius of 21 m and a height of 14 m.
732
733
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
734
<span style="text-align: center; font-size: 75%;">'''Table 5 The resulted micro particle parameters and corresponding shear strangths'''</span></div>
735
736
{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
737
|-
738
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Region</span>
739
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''K<sub>n</sub>''</span><span style="text-align: center; font-size: 75%;"> (N/m)</span>
740
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''K<sub>s</sub>''</span><span style="text-align: center; font-size: 75%;"> (N/m)</span>
741
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''X<sub>b</sub>''</span><span style="text-align: center; font-size: 75%;"> (m)</span>
742
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''F<sub>s0</sub>''</span><span style="text-align: center; font-size: 75%;"> (Pa)</span>
743
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''&#x03bc;<sub>p</sub>''</span>
744
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''C'' (kPa)</span>
745
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''&#x03c6;'' (&#x00b0;)</span>
746
|-
747
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Landslide mass</span>
748
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.27&#x00d7;10<sup>7</sup></span>
749
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.27&#x00d7;10<sup>6</sup></span>
750
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.95&#x00d7;10<sup>-5</sup></span>
751
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.25&#x00d7;10<sup>5</sup></span>
752
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.220</span>
753
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">17.0</span>
754
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">20.2</span>
755
|-
756
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">Sliding zone</span>
757
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.49&#x00d7;10<sup>7</sup></span>
758
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.84&#x00d7;10<sup>6</sup></span>
759
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.72&#x00d7;10<sup>-5</sup></span>
760
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.67&#x00d7;10<sup>5</sup></span>
761
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.120</span>
762
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">24.6</span>
763
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">15.2</span>
764
|-
765
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Bedrock</span>
766
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">4.27&#x00d7;10<sup>7</sup></span>
767
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.09&#x00d7;10<sup>6</sup></span>
768
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.04&#x00d7;10<sup>-5</sup></span>
769
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.95&#x00d7;10<sup>5</sup></span>
770
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.441</span>
771
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">37.5</span>
772
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">32.0</span>
773
|}
774
775
776
With respect to the bedrock, the micro particle parameters in Table 5 resulted in the shear strengths different to the actual shear strengths in Table 4. The reason is that it is hard to represent the three geotechnical materials by using micro parameters of the same order of magnitude. The cohesion of the bedrock is nearly 100 times that of the landslide mass and the sliding zone. Because all the particles representing the bedrock are fixed in the model, the difference between the cohesion of the bedrock in Table 4 and that in Table 5 will have little influence on the computation. Moreover, to reduce the workload, the micro particle parameters of the bedrock will not be changed when the stability of the landslide is evaluated by using SRM.
777
778
==5.4 Computational results==
779
780
As concluded in Section 4, the abrupt in the variation tendencies and the magnitude of the kinetic energy can be utilized to determine the FOS when estimating the slope stability by using MatDEM and SRM. Under different reduction factors, the variations of the kinetic energy with time steps are recorded for the landslide model and plotted in Fig. 15. The kinetic energy curve resulted by a greater reduction factor always has a higher situation in Fig. 15. But, if the reduction factor is less than 1.06, the difference in the magnitude of the kinetic energy between different reduction factors is relatively small at each time step. When the reduction factor is increasing from 1.06 to 1.07, the magnitude difference has a dramatically increasing. Meanwhile, the kinetic energy curves present a convergent trend when the reduction factor is less than 1.06, but present a divergent trend when the reduction factor arrives at 1.07. Thus, the FOS of Baijiabao landslide is considered as 1.06.
781
782
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
783
 [[Image:Review_852894090255-image42.png|318px]] </div>
784
785
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
786
<span style="text-align: center; font-size: 75%;">Fig.15 Variation of kinetic energy with time steps under different reduction factors</span></div>
787
788
<span id='_GoBack'></span>By using the proposed BP neural network-based calibration method, the micro particle parameters are constantly adjusted during the computation. Table 6 lists the micro particle parameters and the shear strengths under different reduction factors.
789
790
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
791
<span style="text-align: center; font-size: 75%;">'''Table 6 Micro particle parameters and corresponding shear strangths under different reduction factors'''</span></div>
792
793
{| style="width: 100%;margin: 1em auto 0.1em auto;border-collapse: collapse;" 
794
|-
795
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Region</span>
796
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">FOS</span>
797
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''K<sub>n</sub>''</span><span style="text-align: center; font-size: 75%;"> (N/m)</span>
798
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''K<sub>s</sub>''</span><span style="text-align: center; font-size: 75%;"> (N/m)</span>
799
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''X<sub>b</sub>''</span><span style="text-align: center; font-size: 75%;"> (m)</span>
800
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''F<sub>s0</sub>''</span><span style="text-align: center; font-size: 75%;"> (Pa)</span>
801
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''&#x03bc;<sub>p</sub>''</span>
802
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''C'' (kPa)</span>
803
|  style="border-top: 2pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">''&#x03c6;'' (&#x00b0;)</span>
804
|-
805
|  rowspan='7' style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Landslide mass</span>
806
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.00</span>
807
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.27&#x00d7;10<sup>7</sup></span>
808
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">5.27&#x00d7;10<sup>6</sup></span>
809
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.95&#x00d7;10<sup>-5</sup></span>
810
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.25&#x00d7;10<sup>5</sup></span>
811
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.220</span>
812
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">17.0</span>
813
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">20.2</span>
814
|-
815
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.10</span>
816
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.35&#x00d7;10<sup>7</sup></span>
817
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4.89&#x00d7;10<sup>6</sup></span>
818
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.61&#x00d7;10<sup>-5</sup></span>
819
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.03&#x00d7;10<sup>5</sup></span>
820
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.175</span>
821
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">15.5</span>
822
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">18.5</span>
823
|-
824
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.03</span>
825
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.29&#x00d7;10<sup>7</sup></span>
826
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">5.13&#x00d7;10<sup>6</sup></span>
827
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.86&#x00d7;10<sup>-5</sup></span>
828
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.18&#x00d7;10<sup>5</sup></span>
829
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.207</span>
830
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">16.5</span>
831
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">19.7</span>
832
|-
833
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.04</span>
834
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.18&#x00d7;10<sup>7</sup></span>
835
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">4.97&#x00d7;10<sup>6</sup></span>
836
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.82&#x00d7;10<sup>-5</sup></span>
837
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.16&#x00d7;10<sup>5</sup></span>
838
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.202</span>
839
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">16.3</span>
840
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">19.5</span>
841
|-
842
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.05</span>
843
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.22&#x00d7;10<sup>7</sup></span>
844
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">5.05&#x00d7;10<sup>6</sup></span>
845
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.77&#x00d7;10<sup>-5</sup></span>
846
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.14&#x00d7;10<sup>5</sup></span>
847
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.198</span>
848
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">16.2</span>
849
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">19.3</span>
850
|-
851
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.06</span>
852
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.53&#x00d7;10<sup>7</sup></span>
853
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">5.02&#x00d7;10<sup>6</sup></span>
854
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.75&#x00d7;10<sup>-5</sup></span>
855
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.12&#x00d7;10<sup>5</sup></span>
856
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.194</span>
857
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">16.0</span>
858
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">19.1</span>
859
|-
860
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.07</span>
861
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.43&#x00d7;10<sup>7</sup></span>
862
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">4.95&#x00d7;10<sup>6</sup></span>
863
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.71&#x00d7;10<sup>-5</sup></span>
864
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.09&#x00d7;10<sup>5</sup></span>
865
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.190</span>
866
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">15.9</span>
867
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">19.0</span>
868
|-
869
|  rowspan='7' style="border-top: 1pt solid black;border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">Sliding zone</span>
870
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.00</span>
871
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.49&#x00d7;10<sup>7</sup></span>
872
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.84&#x00d7;10<sup>6</sup></span>
873
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.72&#x00d7;10<sup>-5</sup></span>
874
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.67&#x00d7;10<sup>5</sup></span>
875
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.120</span>
876
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">24.6</span>
877
|  style="border-top: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">15.2</span>
878
|-
879
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.10</span>
880
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.45&#x00d7;10<sup>7</sup></span>
881
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.81&#x00d7;10<sup>6</sup></span>
882
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.36&#x00d7;10<sup>-5</sup></span>
883
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.05&#x00d7;10<sup>5</sup></span>
884
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.098</span>
885
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">22.4</span>
886
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">13.9</span>
887
|-
888
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.03</span>
889
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.69&#x00d7;10<sup>7</sup></span>
890
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.86&#x00d7;10<sup>6</sup></span>
891
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.61&#x00d7;10<sup>-5</sup></span>
892
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.48&#x00d7;10<sup>5</sup></span>
893
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.113</span>
894
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">23.9</span>
895
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">14.8</span>
896
|-
897
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.04</span>
898
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.38&#x00d7;10<sup>7</sup></span>
899
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.83&#x00d7;10<sup>6</sup></span>
900
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.57&#x00d7;10<sup>-5</sup></span>
901
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.42&#x00d7;10<sup>5</sup></span>
902
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.111</span>
903
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">23.7</span>
904
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">14.6</span>
905
|-
906
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.05</span>
907
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.33&#x00d7;10<sup>7</sup></span>
908
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.96&#x00d7;10<sup>6</sup></span>
909
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.54&#x00d7;10<sup>-5</sup></span>
910
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.36&#x00d7;10<sup>5</sup></span>
911
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.109</span>
912
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">23.4</span>
913
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">14.5</span>
914
|-
915
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.06</span>
916
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.41&#x00d7;10<sup>7</sup></span>
917
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">1.92&#x00d7;10<sup>6</sup></span>
918
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">2.50&#x00d7;10<sup>-5</sup></span>
919
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">3.29&#x00d7;10<sup>5</sup></span>
920
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">0.107</span>
921
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">23.2</span>
922
|  style="text-align: center;"|<span style="text-align: center; font-size: 75%;">14.4</span>
923
|-
924
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.07</span>
925
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.46&#x00d7;10<sup>7</sup></span>
926
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">1.87&#x00d7;10<sup>6</sup></span>
927
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">2.47&#x00d7;10<sup>-5</sup></span>
928
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">3.23&#x00d7;10<sup>5</sup></span>
929
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">0.104</span>
930
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">23.0</span>
931
|  style="border-bottom: 1pt solid black;text-align: center;"|<span style="text-align: center; font-size: 75%;">14.2</span>
932
|}
933
934
935
Fig. 16 illustrates the final displacement of Baijiabao landslide model when the reduction factor is taken as 1.06. Some particles in the front area of the model have the displacements of nearly 60 m, and most particles have the displacements less than 10 m. Thus, the FOS determined by the abrupt in the variation tendencies and the kinetic energy magnitude in fact reveals the landslide stability at the global level. The displacements when the landslide is close to its limit equilibrium state may be valuable in the landslide failure prediction. The field data [46] verified that at the monitoring period from January 2007 to December 2019, the ZG326 site (Fig. 13a) in the middle area of Baijiabao landslide recorded the max cumulative displacement of 1.8 m. Its displacements are still observed when the reservoir water level changes and the rainfall happen, while Baijiabao landslide hasn’t slipped entirely till now. Therefore, the surface movements before its failure are uncertain. In views that this study aims to establish a micro particle parameter calibration method and a criterion for MatDEM coupled with SRM to analysis 3D slope stability, the influence of the reservoir water and the rainfall on the stability of Baijiabao landslide is left out of consideration. Hence, the displacements in the limit equilibrium state obtained in this study are not suitable for the failure prediction for Baijiabao landslide.
936
937
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
938
 [[Image:Review_852894090255-image43-c.png|324px]] </div>
939
940
<div class="center" style="width: auto; margin-left: auto; margin-right: auto;">
941
<span style="text-align: center; font-size: 75%;">Fig.16 The displacement of Baijiabao landslide model with the reduction factor 1.06</span></div>
942
943
=6. Conclusions=
944
945
By taking MatDEM as an example, a calibration method based on BP neural network is proposed to acquire the micro particle parameters corresponding to the given shear strengths in this paper, and an energy criterion is proposed for determining the FOS of the slope when coupling DEM with SRM. And the proposed enhancements are applied in the 3D stability analysis of Baijiabao landslide. The following conclusions can be made:
946
947
(1) BP neural network can well represent the complicated relationship between the micro particle parameters and the shear strengths of the particle aggregate in MatDEM. And the proposed BP neural network-based micro particle parameters calibration method is applicable to provide the micro particle parameters when the shear strength is constantly reduced.
948
949
(2) The variations of the kinetic energy, the gravitational potential energy and the heat are noticed to be quite different between the stable and instable slope. However, the abrupt in the variation tendency and the magnitude of the kinetic energy is more suitable for determining the FOS of the slope.
950
951
(3) When determining the FOS by using the kinetic energy criterion, the displacements of particles in the slope model are permitted, so the FOS actually reveals the slope stability at the global level.
952
953
==Acknowledgments==
954
955
The work is supported by the Open Research Programme of the Hubei Key Laboratory of Disaster Prevention and Mitigation (China Three Gorges University) (no. 2022KJZ07), the Open Research Fund of Key Laboratory of Geological Hazards on Three Gorges Reservoir Area of Ministry of Education (2020KDZ10), National Natural Science Funds (Project No. 52079070), and the Open Research Fund of State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin (China Institute of Water Resources and Hydropower Research), Grant NO:IWHR-SKL-202020.
956
957
==References==
958
959
<span id='_ENREF_1'></span>[1]Ma, Y., Su, P. and Li, Y. Three-dimensional nonhomogeneous slope failure analysis by the strength reduction method and the local strength reduction method. Arabian Journal of Geosciences, 13(1): 21, 2020
960
961
<span id='_ENREF_1'></span>[2]Sun, G., Jiang, W., Cheng, S. and Zheng, H. Optimization Model for Determining Safety Factor and Thrust Line in Landslide Assessments. International Journal of Geomechanics, 17(4): 04016091, 2017.
962
963
<span id='_ENREF_1'></span>[3]Sloan, S.W. Geotechnical stability analysis. Géotechnique, 63(7): 531-571, 2013.
964
965
<span id='_ENREF_1'></span>[4]Lv, Q., Liu, Y. and Yang, Q. Stability analysis of earthquake-induced rock slope based on back analysis of shear strength parameters of rock mass. Engineering Geology, 228: 39-49, 2017.
966
967
<span id='_ENREF_1'></span>[5]Zheng, H., Sun, G. and Liu, D. A practical procedure for searching critical slip surfaces of slopes based on the strength reduction technique. Computers and Geotechnics, 36(1-2): 1-5, 2009.
968
969
<span id='_ENREF_1'></span>[6]Nie, Z., Zhang, Z., Zheng, H. and Lin, S. Stability analysis of landslides using BEM and variational inequality based contact model. Computers and Geotechnics, 123: 103575, 2020.
970
971
<span id='_ENREF_1'></span>[7]Sun, G., Lin, S., Zheng, H., Tan, Y. and Sui, T. The virtual element method strength reduction technique for the stability analysis of stony soil slopes. Computers and Geotechnics, 119: 103349, 2020.
972
973
<span id='_ENREF_1'></span>[8]Espada, M., Muralha, J., Lemos, J.V., Jiang, Q., Feng, X.-T., Fan, Q. and Fan, Y. Safety Analysis of the Left Bank Excavation Slopes of Baihetan Arch Dam Foundation Using a Discrete Element Model. Rock Mechanics and Rock Engineering, 51(8): 2597-2615, 2018.
974
975
<span id='_ENREF_1'></span>[9]Lu, Y., Tan, Y. and Li, X. Stability analyses on slopes of clay-rock mixtures using discrete element method. Engineering Geology, 244: 116-124, 2018.
976
977
<span id='_ENREF_1'></span>[10]Xu, D., Wu, A., Yang, Y., Lu, B., Liu, F. and Zheng, H. A new contact potential based three-dimensional discontinuous deformation analysis method. International Journal of Rock Mechanics and Mining Sciences, 127: 104206, 2020.
978
979
<span id='_ENREF_1'></span>[11]Yang, Y., Xia, Y., Zheng, H. and Liu, Z. Investigation of rock slope stability using a 3D nonlinear strength-reduction numerical manifold method. Engineering Geology, 292: 106285, 2021.
980
981
<span id='_ENREF_1'></span>[12]Wang, B., Vardon, P.J. and Hicks, M.A. Rainfall-induced slope collapse with coupled material point method. Engineering Geology, 239: 1-12, 2018.
982
983
<span id='_ENREF_1'></span>[13]Bui, H.H., Fukagawa, R., Sako, K. and Wells, J.C. Slope stability analysis and discontinuous slope failure simulation by elasto-plastic smoothed particle hydrodynamics (SPH). Géotechnique, 61(7): 565-574, 2011.
984
985
<span id='_ENREF_1'></span>[14]He, L., Tian, Q., Zhao, Z., Zhao, X., Zhang, Q. and Zhao, J. Rock Slope Stability and Stabilization Analysis Using the Coupled DDA and FEM Method: NDDA Approach. International Journal of Geomechanics, 18(6): 04018044, 2018.
986
987
<span id='_ENREF_1'></span>[15]Sun, L., Liu, Q., Abdelaziz, A., Tang, X. and Grasselli, G. Simulating the entire progressive failure process of rock slopes using the combined finite-discrete element method. Computers and Geotechnics, 141: 104557, 2022.
988
989
<span id='_ENREF_1'></span>[16]Bao, Y., Sun, X., Chen, J., Zhang, W., Han, X. and Zhan, J. Stability assessment and dynamic analysis of a large iron mine waste dump in Panzhihua, Sichuan, China. Environmental Earth Sciences, 78(2): 48, 2019.
990
991
<span id='_ENREF_1'></span>[17]Jiang, M., Niu, M., Zhang, F., Wang, H. and Liao, Z. Instability analysis of jointed rock slope subject to rainfall using DEM strength reduction technique. European Journal of Environmental and Civil Engineering, 26(10): 4664-4686, 2021.
992
993
<span id='_ENREF_1'></span>[18]Shi, C., Yang, B., Zhang, Y. and Yang, J. Application of Discrete-Element Numerical Simulation for Calculating the Stability of Dangerous Rock Mass: A Case Study. International Journal of Geomechanics, 20(12): 04020231, 2020.
994
995
<span id='_ENREF_1'></span>[19]Zienkiewicz, O.C., Humpheson, C. and Lewis, R.W. Associated and non-associated visco-plasticity and plasticity in soil mechanics. Geotechnique, 25(4): 671-689, 1975.
996
997
<span id='_ENREF_1'></span>[20]Griffiths, D.V. and Lane, P.A. Slope stability analysis by finite elements. Géotechnique, 49(3): 387-403, 1999.
998
999
<span id='_ENREF_1'></span>[21]Yang, B., Jiao, Y. and Lei, S. A study on the effects of microparameters on macroproperties for specimens created by bonded particles. Engineering Computations, 23(6): 607-631, 2006.
1000
1001
<span id='_ENREF_1'></span>[22]Zhou, Z.H., Wang, H.N. and Jiang, M.J. Macro- and micro-mechanical relationship of the anisotropic behaviour of a bonded ellipsoidal particle assembly in the elastic stage. Acta Geotechnica, 16(12): 3899-3921, 2021.
1002
1003
<span id='_ENREF_1'></span>[23]Cheng, Y. and Wong, L.N.Y. A study on mechanical properties and fracturing behavior of Carrara marble with the flat‐jointed model. International Journal for Numerical and Analytical Methods in Geomechanics, 44(6): 803-822, 2020.
1004
1005
<span id='_ENREF_1'></span>[24]Liu, C., Pollard, D.D. and Shi, B. Analytical solutions and numerical tests of elastic and failure behaviors of close-packed lattice for brittle rocks and crystals. Journal of Geophysical Research: Solid Earth, 118(1): 71-82, 2013.
1006
1007
<span id='_ENREF_1'></span>[25]Rackl, M. and Hanley, K.J. A methodical calibration procedure for discrete element models. Powder Technology, 307: 73-83, 2017.
1008
1009
<span id='_ENREF_1'></span>[26]Hou, J., Zhang, M., Chen, Q., Wang, D., Javadi, A. and Zhang, S.  Failure-mode analysis of loose deposit slope in Ya’an-Kangding Expressway under seismic loading using particle flow code. Granular Matter, 21(1): 8, 2018.
1010
1011
<span id='_ENREF_1'></span>[27]Su, H., Fu, Z., Gao, A. and Wen, Z. Numerical Simulation of Soil Levee Slope Instability Using Particle-Flow Code Method. Natural Hazards Review, 20(2): 04019001, 2019.
1012
1013
<span id='_ENREF_1'></span>[28]Xu, G.J., Zhong, K.Z., Fan, J.W., Zhu, Y.J. and Zhang, Y.Q. Stability analysis of cohesive soil embankment slope based on discrete element method. Journal of Central South University, 27(7): 1981-1991, 2020.
1014
1015
<span id='_ENREF_1'></span>[29]Wang, H., Zhang, B., Mei, G. and Xu, N. A statistics-based discrete element modeling method coupled with the strength reduction method for the stability analysis of jointed rock slopes. Engineering Geology, 264: 105247, 2020.
1016
1017
<span id='_ENREF_1'></span>[30]Deluzarche, R. and Cambou, B. Discrete numerical modelling of rockfill dams. International Journal for Numerical and Analytical Methods in Geomechanics, 30(11): 1075-1096, 2006.
1018
1019
<span id='_ENREF_1'></span>[31]Shen, H. and Abbas, S.M. Rock slope reliability analysis based on distinct element method and random set theory. International Journal of Rock Mechanics and Mining Sciences, 61: 15-22, 2013.
1020
1021
<span id='_ENREF_1'></span>[32]Liu, C., Xu, Q., Shi, B., Deng, S. and Zhu, H. Mechanical properties and energy conversion of 3D close-packed lattice model for brittle rocks. Computers & Geosciences, 103: 12-20,, 2017.
1022
1023
<span id='_ENREF_1'></span>[33]Liu, Y., Zhang, D., Wang, G.-y., Liu, C. and Zhang, Y. Discrete element method-based prediction of areas prone to buried hill-controlled earth fissures. Journal of Zhejiang University-SCIENCE A, 20(10): 794-803, 2019.
1024
1025
<span id='_ENREF_1'></span>[34]Luo, H., Xing, A., Jin, K., Xu, S. and Zhuang, Y. Discrete Element Modeling of the Nayong Rock Avalanche, Guizhou, China Constrained by Dynamic Parameters from Seismic Signal Inversion. Rock Mechanics and Rock Engineering, 54(4): 1629-1645, 2021.
1026
1027
<span id='_ENREF_1'></span>[35]Scaringi, G., Fan, X., Xu, Q., Liu, C., Ouyang, C., Domènech, G., Yang, F. and Dai, L. Some considerations on the use of numerical methods to simulate past landslides and possible new failures: the case of the recent Xinmo landslide (Sichuan, China). Landslides, 15(7): 1359-1375, 2018.
1028
1029
<span id='_ENREF_1'></span>[36]Zhan, Q., Wang, S., Wang, L., Guo, F., Zhao, D. and Yan, J. Analysis of Failure Models and Deformation Evolution Process of Geological Hazards in Ganzhou City, China. Frontiers in Earth Science, 9: 731447, 2021.
1030
1031
<span id='_ENREF_1'></span>[37]Wang, M. and Cao, P. Calibrating the Micromechanical Parameters of the PFC2D(3D) Models Using the Improved Simulated Annealing Algorithm. Mathematical Problems in Engineering, 2017: 6401835, 2017.
1032
1033
<span id='_ENREF_1'></span>[38]Cheng, H., Shuku, T., Thoeni, K. and Yamamoto, H. Probabilistic calibration of discrete element simulations using the sequential quasi-Monte Carlo filter. Granular Matter, 20(1): 11, 2018.
1034
1035
<span id='_ENREF_1'></span>[39]Do, H.Q., Aragón, A.M. and Schott, D.L. A calibration framework for discrete element model parameters using genetic algorithms. Advanced Powder Technology, 29(6): 1393-1403, 2018.
1036
1037
<span id='_ENREF_1'></span>[40]Benvenuti, L., Kloss, C. and Pirker, S. Identification of DEM simulation parameters by Artificial Neural Networks and bulk experiments. Powder Technology, 291: 456-465, 2016.
1038
1039
<span id='_ENREF_1'></span>[41]He, P., Fan, Y., Pan, B., Zhu, Y., Liu, J. and Zhu, D. Calibration and Verification of Dynamic Particle Flow Parameters by the Back-Propagation Neural Network Based on the Genetic Algorithm: Recycled Polyurethane Powder. Materials, 12(20): 3350, 2019.
1040
1041
<span id='_ENREF_1'></span>[42]Ye, F., Wheeler, C., Chen, B., Hu, J., Chen, K. and Chen, W. Calibration and verification of DEM parameters for dynamic particle flow conditions using a backpropagation neural network. Advanced Powder Technology, 30(2): 292-301, 2019.
1042
1043
<span id='_ENREF_1'></span>[43]Zhou, Y., Wu, S., Jiao, J. and Zhang, X. Research on mesomechanical parameters of rock and soil mass based on BP neural network. Rock Soil Mech, 32(12): 3821-3826, 2011.
1044
1045
<span id='_ENREF_1'></span>[44]Graziani, A., Rossini, C. and Rotonda, T. Characterization and DEM Modeling of Shear Zones at a Large Dam Foundation. International Journal of Geomechanics, 12(6): 648-664, 2012.
1046
1047
<span id='_ENREF_1'></span>[45]Peña, J.M., Lozano, J.A. and Larrañaga, P. An empirical comparison of four initialization methods for the K-Means algorithm. Pattern Recognition Letters, 20(10): 1027-1040, 1999.
1048
1049
<span id='_ENREF_1'></span>[46]Yao, W.M., Li, C.D., Guo, Y.C., Criss, R.E., Zuo, Q.J. and Zhan, H.B. Short-term deformation characteristics, displacement prediction, and kinematic mechanism of Baijiabao landslide based on updated monitoring data. Bulletin of Engineering Geology and the Environment, 81(9): 393, 2022.
1050
1051
<span id='_ENREF_1'></span>[47]Xu, S.L. and Niu, R.Q. Displacement prediction of Baijiabao landslide based on empirical mode decomposition and long short-term memory neural network in Three Gorges area, China. Computers & Geosciences, 111: 87-96, 2018.
1052
1053
<span id='_ENREF_1'></span>[48]Yao, W.M., Li, C.D., Zuo, Q.J., Zhan, H.B. and Criss, R.E. Spatiotemporal deformation characteristics and triggering factors of Baijiabao landslide in Three Gorges Reservoir region, China. Geomorphology, 343: 34-47, 2019.
1054
1055
<span id='_ENREF_1'></span>[49]Su, H., Yang, J., Hu, B., Gao, X. and Ma, H. Study of particle size effect of rock model based on particle discrete element method. Rock and Soil Mechanics, 39(12): 4642-4650, 2018.
1056

Return to Jiang et al 2022b.

Back to Top

Document information

Published on 19/01/23
Accepted on 13/01/23
Submitted on 03/12/22

Volume 39, Issue 1, 2023
DOI: 10.23967/j.rimni.2023.01.003
Licence: CC BY-NC-SA license

Document Score

0

Views 190
Recommendations 0

Share this document

claim authorship

Are you one of the authors of this document?