Finite Element Limit Analysis for Solving Different Axisymmetric Stability Problems in Geomechanics : Formulations and Solutions

Abstract

Limit analysis is a very powerful tool to find accurate solutions of several geotechnical stability problems. This analysis is based on the theory of the plasticity and it provides two limiting solutions within lower and upper bounds. With the advancement of the finite elements and different robust optimization techniques, the numerical limit analysis approach in association with finite elements is becoming quite popular to assess the stability of various complicated structures. The present thesis deals with the formulations and the implementation of the finite element limit analysis to obtain the solutions of different geotechnical axisymmetric stability problems. The objectives of the present thesis are twofold: (a) developing limit analysis formulations in conjunction with linear and nonlinear optimizations for solving axisymmetric stability problems related with soil and rock mechanics, and then (b) implementing these axisymmetric formulations for solving various important axisymmetric stability problems in geomechanics. Three noded linear triangular elements have been used throughout the thesis. In order to solve the different problems, the associated computer programs have been written in MATLAB. With reference to the first objective of the thesis, the existing finite element lower bound axisymmetric formulation with linear programming has been presented. A new technique has also been proposed for solving an axisymmetric geomechanics stability problem by employing an upper bound limit analysis in combination with finite elements and linear programming. The method is based on the application of the von-Karman hypothesis to fix the constraints associated with the magnitude of the circumferential stress (), and finally the method involves only the nodal velocities as the basic unknown variables. The required computational effort becomes only marginally greater than that needed for an equivalent plane strain problem. The proposed methodology has been found to be computationally quite efficient. A new lower bound axisymmetric limit analysis formulation, by using two dimensional finite elements, the three dimensional Mohr-Coulomb (MC) yield criterion, and nonlinear optimization has also been presented for solving different axisymmetric stability problems in geomechanics. The nonlinear optimization was carried out by employing an interior point method based on the logarithmic barrier function. The yield surface was smoothened (i) by removing the tip singularity at the apex of the pyramid in the meridian plane, and (ii) by eliminating the stress discontinuities at the corners of the yield hexagon in the plane. No inherent assumption concerning with the hoop stress needs to be made in this formulation. The Drucker-Prager (DP) yield criterion was also used for computing the lower bound axisymmetric collapse load. The advantage of using the DP yield criterion is that it does not exhibit any singularity in the plane. A new proposal has also been given to simulate the DP yield cone with the MC hexagonal yield pyramid. The generalized Hoek-Brown (HB) yield criterion has also been used. This criterion has been smoothened both in the meridian and  planes and a new formulation is prescribed for obtaining the lower bound axisymmetric problems in rock media in combination with finite elements and nonlinear optimization. With reference to the second objective, a few important axisymmetric stability problems in soil mechanics associated with footings and excavations have been solved in the present thesis. In all these problems, except that of a flat circular footing lying over either homogeneous soil or rock media, it is assumed that the medium is governed by the MC failure criterion and it follows an associated flow rule. For determining the collapse loads for a circular footing over homogenous soil and rock media, the problem has been solved with the usage of Drucker-Prager, Mohr-Coulomb and Hoek-Brown criteria. The bearing capacity of a circular footing lying over fully cohesive strata, with an inclusion of a sand layer is evaluated. The effects of the thickness and internal friction angle of the sand layer () on the bearing capacity have been examined for different combinations of cu/(b) and q; where (i) cu defines the undrained shear strength, (ii)  is the unit weight of sand, (iii) b corresponds to the footing radius, and (iv) q is the surcharge pressure. The results have been presented in the form of a ratio () of the bearing capacities with an insertion of the sand layer to that for a footing lying directly over clayey strata. It is noted that an introduction of a layer of medium dense to dense sand over soft clay improves considerably the bearing capacity of the foundation. The improvement in the bearing capacity increases continuously (i) with decreases in cu/(b), and (ii) increases in  and q/(b). The bearing capacity factors, Nc, Nq and N, for a conical footing are obtained in a bound form for a wide range of the values of cone apex angle () and with  = 0, 0.5 and . The bearing capacity factors for a perfectly rough ( = conical footing generally increase with a decrease in . On contrary for  = 0, the factors Nc and Nq reduce gradually with a decrease in . For  = 0, the factor N for  ≥ 35o becomes minimum for  approximately equal to 90o. For  = 0, the factor N for  ≤ 30o, like in the case of  = , generally reduces with an increase in . It has also been intended to compute the bearing capacity factors Nc, Nq and N, for smooth and rough ring footing for different combinations of ri/ro and ; where ri and ro refer to inner and outer radii of the ring, respectively. It is observed that for a smooth footing, with a given value of ro, the magnitude of the collapse load decreases continuously with an increase in ri. On the other hand, for a rough base, for a given value of ro, hardly any reduction occurs in the magnitude of collapse load up to ri/ro ≈ 0.2, whereas beyond this ri/ro, the magnitude of the collapse load, similar to that of a smooth footing, decreases continuously with an increase in ri/ro. An attempt has also been made to determine the ultimate bearing capacity of a circular footing, placed over a soil mass which is reinforced with horizontal layers of circular reinforcement sheets. For performing the analysis, three different soil media have been separately considered, namely, (i) fully granular, (ii) cohesive frictional, and (iii) fully cohesive with an additional provision to account for an increase of cohesion with depth. The reinforcement sheets are assumed to be structurally strong to resist axial tension but without having any resistance to bending; such an approximation usually holds good for geogrid sheets. The shear failure between the reinforcement sheet and adjoining soil mass has been considered. The increase in the magnitudes of the bearing capacity factors (Nc and N) with an inclusion of the reinforcement has been computed in terms of the efficiency factors c and . The critical positions and corresponding optimum diameter of the reinforcement sheets, for achieving the maximum bearing capacity, have also been established. The increase in the bearing capacity with an employment of the reinforcement increases continuously with an increase in . The improvement in the bearing capacity becomes quite extensive for two layers of the reinforcements as compared to the single layer of the reinforcement. The stability of an unsupported vertical cylindrical excavation has been assessed. For the purpose of design, stability numbers (Sn) have been generated for both (i) cohesive frictional soils, and (ii) pure cohesive soils with an additional provision to account for linearly increasing cohesion with depth by using a non-dimensional factor m. The variation of Sn with H/b has been established for different values of m and ; where H and b refer to height and radius of the cylindrical excavation. A number of useful observations have been drawn about the variation of the stability number and nodal velocity patterns with changes in H/b,  and m. In the last, by using the smoothened generalized HB yield criterion, the ultimate bearing capacity of a circular footing placed over a rock mass is evaluated in a non-dimensional form for different values of GSI, mi, ci/(b) and q/ci. For validating the results, computations were exclusively performed for a strip footing as well. For the various problems selected in the present thesis, the failure and nodal velocity patterns have been examined. The results obtained from the analysis have been thoroughly compared with that reported from literature. It is expected that the various design charts presented here will be useful for the practicing engineers. The formulations given in the thesis can also be further used for solving various axisymmetric stability problems in geomechanics.