Study on the Seismic Active Earth Pressure by Variational Limit Equilibrium Method

Abstract

In the framework of limit equilibrium theory, the isoperimetric model of functional extremum regarding the seismic active earth pressure is deduced according to the variational method. On this basis, Lagrange multipliers are introduced to convert the problem of seismic active earth pressure into the problem on the functional extremum of two undetermined function arguments. Based on the necessary conditions required for the existence of functional extremum, the function of the slip surface and the normal stress distribution on the slip surface is obtained, and the functional extremum problem is further converted into a function optimization problem with two undetermined Lagrange multipliers. The calculated results show that the slip surface is a plane and the seismic active earth pressure is minimal when the action point is at the lower limit position. As the action point moves upward, the slip surface becomes a logarithmic spiral and the corresponding value of seismic active earth pressure increases in a nonlinear manner. And the seismic active earth pressure is maximal at the upper limit position. The interval estimation constructed by the minimum and maximum values of seismic active earth pressure can provide a reference for the aseismic design of gravity retaining walls.

[ProQuest: [...] denotes non US-ASCII text; see PDF]

Jiangong Chen 1,2 and Zejun Yang 1,2 and Richeng Hu 1,2 and Haiquan Zhang 1,2

Academic Editor:Sergio De Rosa

1, College of Civil Engineering, Chongqing University, Chongqing 400045, China

2, Key Laboratory of New Technology for Construction of Cities in Mountain Area, Chongqing University, Ministry of Education, Chongqing 400045, China

Received 8 January 2016; Revised 19 April 2016; Accepted 8 May 2016

This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1. Introduction

The magnitude and distribution of active earth pressure on the retaining wall under the seismic loading are the theoretical premises of the aseismic design for the retaining wall and play a vital role in evaluating the stability of the retaining wall in the seismic area. The common calculation method of seismic active earth pressure is the Mononobe-Okabe (M-O) theory which is based on the Coulomb theory and believes that the sliding soil wedge is a rigid body under the seismic loading. The pseudo-static approach is adopted to simplify the seismic force into an inertia force acting on the sliding soil wedge and transcribe the dynamic problem as a static problem [1]. The M-O theory assumes that the backfill behind the wall is cohesion-less soil and the slip surface is a plane, and the theory cannot obtain the real action point of resultant force without considering the equation of moment equilibrium; thus, restrictions are formed on this theory. The M-O theory has been improved by a number of scholars to expand its application scope [2-6].

Both the M-O theory and the improvement method based thereon are on the basis of the assumption that the slip surface is a plane, which does not conform to the practical situation. One rigorous method in math is the variational limit equilibrium method, by which the seismic active earth pressure on the retaining wall is attributed as the functional extremum problem of two undetermined functions. One undetermined function is expressed by the shape of the slip surface, while the other is the function of the normal stress distribution on the slip surface. They are numerically solved by the variational method.

The variational limit equilibrium method was first proposed by Kopascy [7-9]. Then it was introduced to the stability analysis of slope and foundation [10-22]. And some scholars utilized the variational limit equilibrium method to study the lateral earth pressure on the retaining wall [23-26]. According to the variational method, Shaojun studied the shape of the slip surface of the sliding soil wedge behind the retaining wall and obtained the analytical solution of the shape of the slip surface and the magnitude of seismic active earth pressure [27]. However, only the cases where the retaining wall is vertical and the backfill surface is horizontal without surcharge are considered in the calculation models, while the influences of wall-movement modes of the retaining wall on the magnitude and the action point position of seismic active earth pressure fail to be taken into account. In fact, the distribution of earth pressure on the retaining wall is nonlinear. The magnitude and the action point position of seismic active earth pressure depend on the coordinated deformation of soil-wall contact surface and vary with the change of wall-movement modes of the retaining wall [26]. However, in the design process, it is often hard to accurately estimate the wall-movement modes of the retaining wall. For the static and dynamic ultimate load acting on the retaining wall, a reasonable approach is to contain the seismic active earth pressure under different wall-movement modes in a certain range as possible for the engineering designers to select and use. In this paper, the variational limit equilibrium method is used to study the seismic active earth pressure on the gravity retaining wall under general conditions (the retaining wall is inclined and coarse; the backfill is cohesive soil; the backfill surface is a curved surface with nonuniform surcharge). The interval of the seismic active earth pressure under different wall-movement modes can be effectively estimated by the proposed approach.

2. Variational Analysis of Seismic Active Earth Pressure

2.1. Basic Assumption

(1) The research problem is a plane strain problem; (2) the soil behind the wall is Coulomb material, which can be represented by the intensity parameter cohesion c and the internal friction angle [straight phi] ; (3) when the backfill soil is in the critical active state, a sliding soil wedge is formed and its slip surface passes through the wall heel; (4) the retaining wall is rigid and its motion forms are unconstrained, wherein the motion displacement can be ignored compared with the wall height; (5) the seismic action is simplified to static load acting on the sliding wedge, with horizontal seismic coefficient _{k H} and vertical seismic coefficient _{k V} .

2.2. Limit Equilibrium Equation of Sliding Soil Wedge

The calculation model of active earth pressure under the seismic loading is shown in Figure 1, wherein the height of the retaining wall is H ; the retaining wall is inclined and coarse; the slope angle of the wall to vertical is α ; the friction angle between soil and wall is δ ; the unit weight of soil is γ ; the cohesion is c ; the internal fraction angle is [straight phi] , the expression of the backfill surface is g ( x ) , and the expression of the slip surface is s ( x ) ; the expression of the vertical surcharge distribution on the backfill surface is q ( x ) , the expression of the tangential stress distribution on the slip surface is τ ( x ) , and the expression of the normal stress distribution on the slip surface is σ ( x ) ; the seismic soil pressure on the wall is _{E a} .

Figure 1: Calculation model of seismic active earth pressure.

[figure omitted; refer to PDF]

The equilibrium equation of the sliding soil wedge OAB under the limit equilibrium state can be written as follows: [figure omitted; refer to PDF] wherein ^{s [variant prime]} = d s / d x , k = - cot [...] α , _{x 1} is X coordinate of B , _{x 2} is X coordinate of A , and _{x 2} = - H tan [...] α . ξ is the position coefficient of the action point of seismic active earth pressure, and its value is the ratio of the vertical distance (from the action point position to the wall heel) to the height of the retaining wall.

It is assumed that the tangential stress τ and the normal stress σ on the slip surface are subject to the Mohr-Coulomb failure criterion: [figure omitted; refer to PDF] wherein _{n 1} = tan [...] [straight phi] .

2.3. Variational Limit Equilibrium Model of Seismic Active Earth Pressure

The problem to solve the seismic active earth pressure _{E a e} can now be stated as follows. The seismic active earth pressure _{E a e} is the maximum value of _{E a} , corresponding to the most dangerous (critical) slip surface. Realizing the maximum value of _{E a} is to find the function of the slip surface s ( x ) and the function of normal stress distribution σ ( x ) on the slip surface, subject to the three equations of LE ((1)~(3)). Among these three equations, (3) is changed as objective function and the other equations are the constraint conditions according to the variational method.

The following objective function J can be obtained after (4) is substituted into (3). Obviously, to solve the maximum of _{E a} means the same as to solve the maximum of objective function J : [figure omitted; refer to PDF] wherein [figure omitted; refer to PDF] Two constraint conditions can be obtained after (4) and (5) are substituted into (1) and (2): [figure omitted; refer to PDF] wherein [figure omitted; refer to PDF] Equations (5) and (7) show the isoperimetric model of the constrained variational extremum with an undetermined boundary. The starting point of the slip surface O is a fixed point on the coordinates of ( _{x 0} = 0 , _{y 0} = 0 ) and the end point B moves on the backfill surface g ( x ) with the undetermined coordinates of ( _{x 1} , g ( _{x 1} ) ) . The undetermined boundary is a variational boundary. The variational boundary refers to an unfixed boundary on which one or two ends, respectively, move on the given function.

According to the variational method of the functional with constraints, the following auxiliary functional ^{J [low *]} is constructed by Lagrange multipliers to convert the functional extremum problem under constraint conditions into a functional extremum problem without constraint conditions: [figure omitted; refer to PDF] wherein _{λ 1} and _{λ 2} are Lagrange multipliers.

Now, the solution of the maximum of _{E a} is converted into the solution of the maximum of auxiliary function ^{J [low *]} . According to the necessary conditions for the existence of the extremum of auxiliary functional ^{J [low *]} , the function of the slip surface s ( x ) and the function of the normal stress distribution σ ( x ) on the slip surface must meet the following requirements.

(1) The Euler differential equation of the auxiliary function F is [figure omitted; refer to PDF]

(2) The integral constraint equations are (7).

(3) Boundary conditions are [figure omitted; refer to PDF]

(4) Transversality condition at the variational boundary is [figure omitted; refer to PDF]

3. Variational Solution of Seismic Active Earth Pressure

3.1. Shape Function of the Slip Surface

The following equation can be obtained after (5), (7), and (9) are substituted into (10): [figure omitted; refer to PDF] Polar coordinate transformations r cos [...] θ = x + _{λ 2} and r sin [...] θ = s - _{λ 1} are brought in. For the calculation model in polar coordinates, see Figure 2. Equation (14) is transcribed as a differential equation: [figure omitted; refer to PDF] The general solution of the differential equation is [figure omitted; refer to PDF] wherein _{z 5} is an arbitrary integration constant.

Figure 2: Calculation model of seismic active earth pressure in polar coordinates.

[figure omitted; refer to PDF]

The fixed boundary condition s ( 0 ) = 0 shows that the pole of the polar coordinate is ( _{r 0} , _{θ 0} ) . After it is substituted into (16), the expression of the logarithmic spiral of the slip surface can be obtained: [figure omitted; refer to PDF] wherein [figure omitted; refer to PDF]

3.2. Normal Stress Distribution on the Slip Surface

After (5), (7), and (9) are substituted into (11), the following equation can be obtained: [figure omitted; refer to PDF] Introducing the Polar coordinate transformations r cos [...] θ = x + _{λ 2} and r sin [...] θ = s - _{λ 1} , (19) is transcribed as a differential equation: [figure omitted; refer to PDF] The general solution of the differential equation is (when _{n 1} ≠ 0 ) [figure omitted; refer to PDF] wherein _{z 6} is an integration constant; _{θ 1} is an arbitrary angle; _{θ 1} = 0 is available [figure omitted; refer to PDF] According to the transversality condition at the variational boundary, the normal stress distribution at B of the slip surface can be obtained: [figure omitted; refer to PDF] After it is substituted into (22), the following equation is obtained: [figure omitted; refer to PDF]

3.3. Optimal Solution of Seismic Active Earth Pressure

The analysis carried out so far indicates that the maximum value of _{E a} can be obtained according to a pair of functions r ( θ ) and σ ( θ ) determined by two undetermined Lagrange multipliers _{λ 1} and _{λ 2} . The problem is converted into the search of a pair of constants _{λ 1} and _{λ 2} through the two constraint equations. The two constraint equations-(7)-must be satisfied: [figure omitted; refer to PDF] The following equation is always correct: [figure omitted; refer to PDF] Therefore, two integral equations-(25)-are equivalent to the following equation: [figure omitted; refer to PDF] For Φ ( _{λ 1} , _{λ 2} ) ≥ 0 , so [figure omitted; refer to PDF] The above analysis indicates that the solution of the two constraint equations-(7)-can be obtained through the solution of the minimum value of the function Φ (the minimum value is 0). When the minimum value of Φ is not 0, it shows that the sliding soil wedge cannot keep balance and is improper.

In this paper, the fminsearch function provided by MATLAB is used to find the optimal solution. For this function, it only requires establishing a function module of Φ and giving an initial value to search the minimum value of the function. Fminsearch applies the derivative-free method to find the minimum value of the unconstrained multivariable function, which is generally called unconstrained nonlinear optimization. Fminsearch finds the minimum of a scalar function of several variables starting at an initial estimate. Therefore, one major drawback of this optimization tool is that the search would be very slow if the initial estimate is improper. In order to overcome this drawback and in the consideration of the condition that the center of the logarithmic spiral is usually at the top left corner of the retaining wall, the initial value is _{λ 1} = _{λ 2} = H in the optimization analysis of this paper and the search is quite rapid.

4. Calculated Result and Parameter Analysis

4.1. Interval Estimation of Seismic Active Earth Pressure

Under the general conditions that the retaining wall is inclined and coarse, the backfill is cohesive soil, and the backfill surface is inclined with surcharge, the interval of the position coefficient ξ of the action point under seismic loading is numerically solved. Besides, the influence of the action point position on the shape of the slip surface and the magnitude of seismic active earth pressure is studied.

The retaining wall has a height of 10 m, with the slope angle of the wall to vertical α = 1 ^{0 [composite function]} , the friction angle between soil and wall δ = 1 ^{5 [composite function]} , the unit weight of soil γ = 18 k N / ^{m 3} , cohesion c = 10 k P a , internal friction angle [straight phi] = ^{20 °} , slope angle of the backfill soil β = ^{5 °} , surcharge q = 50 k N / m , and seismic coefficients _{k H} = 0.05 and _{k V} = 0.0 .

Figure 3 shows the curve of the minimum value of the function Φ varying with the position coefficient ξ of the action point. As shown in the figure, the position coefficient ξ of the action point of seismic active earth pressure is in an interval (lower limit value _{ξ d} = 0.449 ; upper limit value _{ξ u} = 0.612 ). In the interval, the minimum value of the function Φ is 0. In other words, when the action point of seismic active earth pressure falls within the range of 4.49~6.12 m from the wall bottom, the counterforce provided by the retaining wall can balance the soil mass behind the wall; otherwise, the soil mass will inevitably lose its balance no matter how big the counterforce provided by the retaining wall is.

Figure 3: The minimum value of the function Φ .

[figure omitted; refer to PDF]

There is a one-to-one corresponding relation between the magnitude of seismic active earth pressure and the shape of slip surface (Figure 5). For the curves varying with the position coefficient ξ of the action point of seismic active earth pressure, see Figures 3 and 4. At the lower limit of ξ , the slip surface is a plane, with the minimal value of seismic active earth pressure. As the action point moves upward, the slip surface becomes a logarithmic spiral surface and its curvature is gradually increased; that is, in the logarithmic spiral equation, _{r 0} , is gradually decreased, and the corresponding value of seismic active earth pressure is increased in a nonlinear manner. At the upper limit of ξ , the seismic active earth pressure is maximal. The minimum and maximum values can be as an interval estimation of the seismic active earth pressure under different wall-movement modes of the retaining wall.

Figure 4: The shape of the slip surface.

[figure omitted; refer to PDF]

Figure 5: The magnitude of seismic active earth pressure.

[figure omitted; refer to PDF]

4.2. Comparison with the M-O Theory Solution

To verify the correctness of the calculation method proposed in this paper, a comparison is made between the calculated result of the proposed method and the solution of M-O theory. According to the assumption of M-O theory, the values of cohesion c and the surcharge q are both 0 in the contrastive analysis. For other parameters and specific comparison results, see Table 1.

Table 1: Comparison of calculated results of active earth pressure with that of M-O.

α /°	β /°	k H	k V	The proposed method	M-O method
E a e /kN	ξ	E a /kN	ξ
0	0	0	0	391.0~412.1	0.265~0.408	391.0	1/3
0	5	0.05	0.05	424.6~450.2	0.286~0.449	424.5	1/3
0	10	0.1	0.1	470.8~495.0	0.306~0.490	470.7	1/3
10	0	0.05	0.1	460.5~475.2	0.327~0.408	460.1	1/3
10	5	0.1	0	504.3~525.3	0.327~0.449	504.3	1/3
10	10	0	0.05	559.9~572.8	0.367~0.469	559.6	1/3
20	0	0.1	0.05	550.4~562.8	0.347~0.408	550.3	1/3
20	5	0	0.1	604.9~616.1	0.367~0.429	604.7	1/3
20	10	0.05	0	682.7~696.3	0.388~0.469	682.5	1/3

When the slip surface is a plane, the seismic active earth pressure is the lower limit value of the interval in the table. This value is equivalent to the solution of seismic active earth pressure calculated by the M-O theory. The solution obtained in the proposed method is degraded into M-O theory solution, thereby proving the correctness of the proposed method.

When the slip surface is a logarithmic spiral surface, the seismic active earth pressure calculated by the M-O theory is relatively small, and the assumed action point position is lower than the actual action point position under most possible wall-movement modes. The underestimation of the magnitude and the action point position of seismic active earth pressure cause the potential safety hazard of overturning of the retaining wall designed by M-O theory.

4.3. Parameter Analysis

In this section, discussions are made regarding the influence of relevant parameters on the magnitude and the action point position of seismic active earth pressure, as well as the overturning moment on the retaining wall. Such parameters include the slope angle α of the retaining wall, the slope angle β of backfill surface, the surcharge q on the backfill surface, the horizontal seismic coefficient _{k H} , and the vertical seismic coefficient _{k V} .

4.3.1. Influence of Parameters on the Action Point Position

The influence of relevant parameters on the position coefficient ξ of action point of seismic active earth pressure is shown in Figures 6(a)-6(c). With the increase of the slope angle β of backfill surface, the action point position moves up; with the increase of α , the action point position slightly moves down; with the increase of the horizontal seismic coefficient _{k H} and the vertical seismic coefficient _{k V} , the action point position moves up, wherein the horizontal seismic coefficient _{k H} has a significant influence on it; with the increase of the surcharge q , the action point position moves up in a nonlinear way.

Figure 6: Effects of relevant parameters on ξ .

(a) α and β

[figure omitted; refer to PDF]

(b) _{k H} and _{k V}

[figure omitted; refer to PDF]

(c) q

[figure omitted; refer to PDF]

4.3.2. Influence of Parameters on the Seismic Active Earth Pressure _{E a e}

Figures 7(a)-7(c) show the influence of relevant parameters on the dimensionless seismic active earth pressure _{E a e} / ( γ ^{H 2} ) . The seismic active earth pressure _{E a e} is increased with the increase of the slope angle α of the retaining wall and the slope angle β of backfill surface; _{E a e} is increased with the increase of the horizontal seismic coefficient _{k H} and decreased with the increase of the vertical seismic coefficient _{k V} ; _{E a e} is increased linearly with the increase of the dimensionless surcharge q / ( γ H ) of slope surface.

Figure 7: Effects of relevant parameters on _{E a e} .

(a) α and β

[figure omitted; refer to PDF]

(b) _{k H} and _{k V}

[figure omitted; refer to PDF]

(c) q

[figure omitted; refer to PDF]

4.3.3. Influence of Parameters on the Overturning Moment _{M a}

Figures 8(a)-8(c) show the influence of relevant parameters on the dimensionless overturning moment _{M a} / γ ^{H 3} . The overturning moment _{M a} is increased with the increase of the slope angle α of the retaining wall and the slope angle β of backfill surface; _{M a} is increased with the increase of the horizontal seismic coefficient _{k H} and decreased with the increase of the vertical seismic coefficient _{k V} ; _{M a} is increased linearly with the increase of the dimensionless surcharge load q / ( γ H ) of slope surface.

Figure 8: Effects of relevant parameters on _{M a} .

(a) α and β

[figure omitted; refer to PDF]

(b) _{k H} and _{k V}

[figure omitted; refer to PDF]

(c) q

[figure omitted; refer to PDF]

5. Conclusion

The seismic active earth pressure is studied in the paper based on the variational limit equilibrium method, and the following conclusions are obtained.

(1) In the framework of limit equilibrium theory, the isoperimetric model of functional extremum regarding the seismic active earth pressure under general conditions (the retaining wall is inclined and coarse; the backfill is cohesive soil; the backfill surface is a curved surface with nonuniform surcharge) is deduced.

(2) With the introduction of Lagrange multipliers and in combination with the necessary conditions required for the existence of functional extremum, the solution of seismic active earth pressure is transcribed as an optimization problem of two undetermined Lagrange multipliers; with the fminsearch function provided by MATLAB, an optimization solution of two Lagrange multipliers is founded to obtain the magnitude and the action point position of seismic active earth pressure, as well as the shape of slip surface.

(3) With the continuous changes of wall-movement modes of the retaining wall, the position coefficient ξ of the action point of seismic active earth pressure is an interval with lower and upper limit values. When ξ is the lower limit value, the slip surface is a plane and the seismic active earth pressure is minimal. As the action point moves upward, the slip surface becomes a logarithmic spiral surface and its curvature is gradually increased, with the corresponding value of seismic active earth pressure increased in a nonlinear manner. At the upper limit of ξ , the seismic active earth pressure is maximal; the minimum and maximum values can be estimated as an interval of the seismic active earth pressure under different wall-movement modes of the retaining wall.

(4) The contrastive analysis shows that when the slip surface is a plane, the solution of seismic active earth pressure obtained by the proposed method is in line with the M-O theory solution, verifying the correctness of the proposed method; when the slip surface is a logarithmic spiral surface, the underestimation of the magnitude and the action point position of seismic active earth pressure cause the potential safety hazard of overturning of the retaining wall designed by M-O theory.

(5) The parameter analysis shows that the action point position moves up with the increase of the slope angle β of backfill surface, the horizontal seismic coefficient _{k H} , and the surcharge q on the backfill surface. However, the vertical seismic coefficient _{k V} and the slope angle α of the retaining wall have little influence on it, so they can be neglected. The seismic active earth pressure _{E a e} and the overturning moment _{M a} are increased with the increase of the slope angle α of the retaining wall, the slope angle β of backfill surface, the horizontal seismic coefficient _{k H} , and the surcharge q on the backfill surface and decreased with the increase of the vertical seismic coefficient _{k V} .

References

[1] S. Caltabiano, E. Cascone, M. Maugeri, "Static and seismic limit equilibrium analysis of sliding retaining walls under different surcharge conditions," Soil Dynamics and Earthquake Engineering , vol. 37, pp. 38-55, 2012.

[2] S. Saran, R. P. Gupta, "Seismic earth pressures behind retaining walls," Indian Geotechnical Journal , vol. 33, no. 3, pp. 195-213, 2003.

[3] G. Mylonakis, P. Kloukinas, C. Papantonopoulos, "An alternative to the Mononobe-Okabe equations for seismic earth pressures," Soil Dynamics and Earthquake Engineering , vol. 27, no. 10, pp. 957-969, 2007.

[4] S. Ghosh, S. Sengupta, "Extension of Mononobe-Okabe theory to evaluate seismic active earth pressure supporting c-[straight phi] backfill," Electronic Journal of Geotechnical Engineering , vol. 17, pp. 495-504, 2012.

[5] M. Yazdani, A. Azad, A. H. Farshi, S. Talatahari, "Extended 'mononobe-okabe' method for seismic design of retaining walls," Journal of Applied Mathematics , vol. 2013, 2013.

[6] Y.-L. Lin, W.-M. Leng, G.-L. Yang, L.-H. Zhao, L. Li, J.-S. Yang, "Seismic active earth pressure of cohesive-frictional soil on retaining wall based on a slice analysis method," Soil Dynamics and Earthquake Engineering , vol. 70, pp. 133-147, 2015.

[7] J. Kopascy, K. Szechy, "Über die Bruchflachen und Bruchspannungen in den Erdbauten," Gedenkbuch Für Prof. Dr. J. Jaky , pp. 81-99, Akademiaikiado, Budapest, Hungary, 1955.

[8] J. Kopascy, "Three-dimensional stress distribution and slip surfaces in earth works at rupture," in Proceedings of the 4th International Conference on Soil Mechanics and Foundations Engineering, vol. 1, pp. 339-342, London, UK, 1957.

[9] J. Kopascy, "Distribution des contraintesala rupture, forme de la surface de glissement et hauteur theorique des talus," in Proceedings of the 5th International Conference on Soil Mechanics and Foundations Engineering, vol. 2, pp. 641-650, Paris, France, 1961.

[10] W. F. Chen, N. Snitbhan, "On slip surface and slope stability analysis," Soils and Foundations , vol. 15, no. 3, pp. 41-49, 1975.

[11] M. Garber, R. Baker, "Bearing capacity by variational method," Journal of Geotechnical Engineering Division , vol. 103, no. 11, pp. 1209-1225, 1977.

[12] M. Garber, R. Baker, "Extreme-value problems of limiting equilibrium," Journal of the Geotechnical Engineering Division , vol. 105, no. 10, pp. 1155-1171, 1979.

[13] A. J. Spencer, T. C. O'Mahony, "An application of the calculus of variations to rectilinear flow of granular materials," International Journal for Numerical and Analytical Methods in Geomechanics , vol. 9, no. 3, pp. 225-235, 1985.

[14] D. Leshchinsky, "Slope stability analysis: generalized approach," Journal of Geotechnical Engineering , vol. 116, no. 5, pp. 851-867, 1990.

[15] D. Leshchinsky, C.-C. Huang, "Generalized three-dimensional slope-stability analysis," Journal of Geotechnical Engineering , vol. 118, no. 11, pp. 1748-1764, 1992.

[16] R. Baker, "Sufficient conditions for existence of physically significant solutions in limiting equilibrium slope stability analysis," International Journal of Solids and Structures , vol. 40, no. 13-14, pp. 3717-3735, 2003.

[17] R. Baker, "Stability chart for zero tensile strength Hoek-Brown materials-the variational solution and its engineering implications," Soils and Foundations , vol. 44, no. 3, pp. 125-132, 2004.

[18] R. Baker, "Variational slope stability analysis of materials with nonlinear failure criterion," Electronic Journal of Geotechnical Engineering , vol. 10, pp. 1-22, 2005.

[19] R. Baker, "A relation between safety factors with respect to strength and height of slopes," Computers and Geotechnics , vol. 33, no. 4-5, pp. 275-277, 2006.

[20] L.-Y. Wu, Y.-F. Tsai, "Variational stability analysis of cohesive slope by applying boundary integral equation method," Journal of Mechanics , vol. 21, no. 3, pp. 187-198, 2005.

[21] Y. M. Cheng, D. Z. Li, N. Li, Y. Y. Lee, S. K. Au, "Solution of some engineering partial differential equations governed by the minimal of a functional by global optimization method," Journal of Mechanics , vol. 29, no. 3, pp. 507-516, 2012.

[22] X. Li, "Bearing capacity factors for eccentrically loaded strip footings using variational analysis," Mathematical Problems in Engineering , vol. 2013, 2013.

[23] A. H. Soubra, R. Kastner, "Influence of seepage flow on the passive earth pressures," in Proceeding of the International Conference on Retaining Structures, pp. 67-76, ICE, Amsterdam, The Netherlands, 1992.

[24] A.-H. Soubra, R. Kastner, A. Benmansour, "Passive earth pressures in the presence of hydraulic gradients," Geotechnique , vol. 49, no. 3, pp. 319-330, 1999.

[25] O. Pula, W. Pula, A. Wolny, "On the variational solution of a limiting equilibrium problem involving an anchored wall," Computers and Geotechnics , vol. 32, no. 2, pp. 107-121, 2005.

[26] L. Xinggao, L. Weining, "Study on the action of the active earth pressure by variational limit equilibrium method," International Journal for Numerical and Analytical Methods in Geomechanics , vol. 34, no. 10, pp. 991-1008, 2010.

[27] M. Shaojun 2012 (Chinese) Study on Calculation of Earth Pressure behind Retaining Wall under Static and Dynamic Loads and Its Related Issues , Zhejiang University, Zhejiang, China