This thesis presents the implementation of the upper bound limit analysis in combination with finite elements and linear optimization for solving different stability problems in geomechanics under plane strain conditions. Although the nonlinear optimization techniques are becoming quite popular, the linear optimization has been adopted due to its simplicity in implementation and ease in attaining the convergence while performing the analysis. The objectives of the present research work are (i) to reduce the computational effort while using an upper bound finite element limit analysis with linear programming in dealing with geotechnical stability problems, and (ii) to obtain solutions for a few important geotechnical stability problems associated with reinforced earth, unsupported tunnels and a group of anchors. It is also intended to examine the developments of the failure patterns in all the cases. For carrying out the analysis for different stability problems, three noded triangular elements have been used throughout the thesis. The nodal velocities are treated as basic unknown variables and the velocity discontinuities are employed along the interfaces of all the elements. The soil mass is assumed to obey the Mohr-Coulomb’s failure criterion and an associated flow rule. The Mohr-Coulomb yield surface is linearized by means of an exterior regular polygon circumscribing the actual yield circle so that the finite element formulation leads to a linear programming problem. A simple technique has been proposed for reducing the computational effort while solving any geotechnical stability problem by using the upper bound finite element limit analysis and linear optimization. In the proposed method, the problem domain has been discretized into a number of different regions in which a particular order (number of sides) of the polygon has been specified to linearize the Mohr-Coulomb yield criterion. A greater order of the polygon needs to be chosen only in that part of the domain wherein the rate of the plastic strains becomes higher. The computational effort required to solve the problem with this implementation reduces considerably. By using the proposed method, the bearing capacity has been computed for smooth as well as rough strip footings and the results obtained are found to be quite satisfactory. The ultimate bearing capacity of a rigid strip footing placed over granular, cohesive-frictional and purely cohesive soils, reinforced with single and a group of two horizontal layers of reinforcements has been determined. The necessary formulation has been introduced to incorporate the inclusion of reinforcement in the analysis. The efficiency factors, and , to be multiplied with Nc and Nγ for finding the bearing capacity of reinforced foundations, have been established. The results have been obtained (i) for different values of soil friction angles in case of granular and cohesive-frictional soils, and (ii) for different rates at which the cohesion increases with depth for purely cohesive soil under undrained condition. The optimum positions of the reinforcements' layers corresponding to which and becomes maximum, have been established. The effect of the length of the reinforcements on the results has also been analyzed. As compared to cohesive soil, the granular soils, especially with greater values of frictional angle, cause much more predominant increase in the bearing capacity. The stability of a long open vertical trench laid in a fully cohesive and cohesive-frictional soil has been determined with an inclusion of single and a group of two layers of horizontal reinforcements. For different positions of the reinforcement layers, the efficiency factor (ηs), has been determined for several combinations of H/B, m and where H and B refer to height and width of the trench, respectively, and m accounts for the rate at which the cohesion increases linearly with depth for a fully cohesive soil with = 0. The effect of height to width of the long vertical trench on the stability number has been examined for both unreinforced and reinforced soils. The optimal positions of the reinforcements layers, corresponding to which becomes maximum, have been established. The required length of reinforcements to achieve maximum efficiency factor corresponding to optimum depth of reinforcement has also been determined. The magnitude of the maximum efficiency factor increases continuously with an increase in both m and . The effect of pseudo-static horizontal earthquake body forces on the stability of a long unsupported circular tunnel (opening) formed in a cohesive frictional soil has been determined. The stability numbers have been obtained for various values of H/D (H = tunnel cover, D = diameter of the tunnel), internal friction angle of soil, and the horizontal earthquake acceleration coefficient The computations revealed that the values of the stability numbers (i) decreases quite significantly with an increase in , and (ii) become continuously higher for greater values of H/D and . The failure patterns have also been drawn for different combinations of H/D, and . The geometry of the failure zone around the periphery of the tunnel becomes always asymmetrical with an inclusion of horizontal seismic body forces. The interference effect on the stability of two closely spaced parallel (twin) long unsupported circular tunnels formed in fully cohesive and cohesive-frictional soils has been evaluated. The variation of the stability number with S/D has been established for different combinations of H/D, m and ; where D refers to the diameter of each tunnel, S is the clear spacing between the tunnels, and is the internal friction angle of soil and m accounts for the rate at which the cohesion increases linearly with depth for a soil with = 0. On account of the interference of two tunnels, the stability number reduces continuously with a decrease in the spacing between the tunnels. The minimum spacing between the two tunnels required to eliminate the interference effect increases with (i) an increase in H/D and (ii) a decrease in the values of both m and . The failure patterns have also been generated for a few cases with different values of S/D. The size of the failure zone is found to become smaller for greater values of m and . The horizontal pullout capacity of a group of two vertical strip anchors embedded, along the same vertical plane in sand, at shallow depths has been determined. At collapse, it is assumed that the anchor plates are subjected to the same uniform horizontal velocity without any bending or tilt. The pullout resistance increases invariably with increases in the values of embedment ratio, friction angle of the sand mass and anchor-soil interface friction angle. The effect of spacing (S) between the anchors on their group collapse load is examined in detail. For a given embedment ratio, the total group failure load becomes maximum corresponding to a certain optimal spacing (Sopt). The values of Sopt increases with an increase in the value of , but the changes in the value of H/B and do not have any significant effect on Sopt. The vertical uplift capacity of a group of two horizontal strip plate anchors with the common vertical axis buried in purely cohesive as well as in cohesive frictional soil has been computed. The variation of the uplift factors Fc, Fq and F , due to the contributions of soil cohesion, surcharge pressure and unit weight, respectively, has been evaluated for different combinations of S/B and H/B. As compared to a single isolated anchor, the group of two anchors generates significantly greater magnitude of Fc. On the other hand, the factors Fq and F , for a group of two anchors are found to become almost equal to that of a single isolated anchor as long as the levels of the lower plate in the group and the single isolated anchor are kept the same. For the group of two horizontal strip plate anchors in purely cohesive soil, an increase of cohesion of soil mass with depth and the effect of self weight of the soil have been incorporated. The uplift factor Fcy both due to cohesion and unit weight of the soil has also been computed for the anchors embedded in clay under undrained condition. For given embedment ratios, the factor Fcy increases linearly with an increase in the normalized unit weight of soil mass upto a certain value before attaining a certain maximum magnitude. The computational results obtained for different research problems would be useful for design.

The present thesis is an attempt to investigate the interference effects on the magnitudes of the ultimate failure loads for a group of closely spaced strip footings and strip plate anchors. On account of an increase in the number of different civil engineering structures, footings and anchors are often need to be placed very close to each other. In such a situation, the ultimate bearing capacity/pullout capacity of an interfering footing/anchor becomes significantly different from that of a single isolated footing/anchor. The effect of interference on the magnitude of failure load is usually expressed in terms of an efficiency factor (%y); where £,y is defined as the ratio of the magnitude of the failure load for a strip footing/anchor of a given width in the presence of other footings/anchors to that of the magnitude of the failure load for an isolated single strip footing/anchor having exactly the same width. No rigorous analysis seems to have been carried out so far in literature to investigate the interference effect for a group of footings and anchors. In the present study, it is intended to use rigorous numerical upper bound limit analysis in combination with finite elements and linear programming in order to determine the collapse loads for the problems of both isolated and a group of footings and anchors. Three noded triangular elements are used throughout the thesis for carrying out the analysis for different problems. The velocity discontinuities are employed along the interfaces of all the elements. The plastic strains within the elements are incorporated by using an associated flow rule. The Mohr Coulomb yield surface is linearised by means of an exterior regular polygon circumscribing the actual failure surface so that the finite element formulation leads to a linear programming problem. In solving the different problems taken in this thesis, computer programs were developed using 'MATLAB' with the usage of 'LINPROG' - a library subprogram for doing the necessary optimization. The bearing capacity factor Ny for an isolated single rigid strip footing placed on a cohesionless ground surface has been computed and its variation with respect to the footing-soil roughness angle (8) has been examined in detail. It is clearly noted that an increase in 8 leads to a continuous increase in Ny. The solution is also obtained for a perfectly rough footing without considering any velocity discontinuity surface along the footing-soil interface. With 5 = <|), the magnitude of NY becomes almost the same as that for a perfectly rough footing. The size of the plastic zone increases with an increase in the values of 8 and <j). The obtained values of Ny for 5=0 and § compare quite favorably with the solutions reported earlier in literature. The ultimate bearing capacity for a group of two and an infinite number of multiple interfering rough strip footings placed on a cohesionless medium has been computed; all the footings are assumed to be perfectly rigid. It is specified that the footings are loaded simultaneously to failure exactly at the same magnitude of the failure load. For different clear spacing (S) between the adjacent footings, the magnitude of the efficiency factor (£,y) is determined. In the case of two footings, the value of E,y at S/B = 0 becomes exactly equal to 2.0, and the maximum ^occurs at a critical spacing (Scr). For S/B < Sor/B, the ultimate bearing pressure for a footing becomes equal to that of an isolated footing having the width (2B+S), and the ground mass encompassed between the two footings deforms mainly in the downward direction. In contrast, for S/B > Scr/B, ground heave is noticed along both the sides of the footing. As compared to the available theories in literature, the analysis presented in this thesis provides generally lower values of ^y for S/B > Scr/B. ' In the case of a group of multiple strip footings, the value of £y is found to increase continuously with a decrease in S/B. The effect of the variation of spacing on §y is found to be very extensive for small values of S/B; the magnitude of ^y approaches infinity at S/B = 0. For all the values of S/B ground heave is invariably observed on both the sides of the footings. The magnitudes of ^Y for given values of S/B and <|) for the two footings case are found to be smaller than the multiple footings case. The vertical uplift capacity of an isolated strip anchor embedded horizontally at shallow depths in sand has been examined; the anchor plate is assumed to be perfectly rigid and rough. The collapse load is expressed in terms of a non-dimensional uplift factor FY, the value of which needs to be known before calculating the failure load for an interfering anchor. The magnitude of Fr is found to increase continuously with increase in both embedment ratio (k) and the friction angle (<|>) of sand. Even though the analysis considers the development of plastic strain within all elements, however, at collapse, the soil mass just above the anchor is found to move as a single rigid block bounded by planar rupture surfaces; the rupture surfaces emerging from the anchor edges are seen to make approximately an angle <|> with the vertical. The vertical uplift capacity of a group of two and an infinite number of multiple interfering rigid rough strip anchors embedded horizontally in sand at shallow depths has been examined. At collapse, it is specified that all the anchors in the group are loaded to failure simultaneously exactly at the same magnitude of the failure load. For different clear spacing (S) between the anchors, the magnitude of the efficiency factor (£Y) is determined. On account of interference, the magnitude of 4y is found to reduce continuously with a decrease in the spacing between the anchors. For all values of X and §, the magnitude of ^y for the multiple anchors case is found to be always smaller than that for the two anchors case. In contrast to a group of footings under compression, the magnitude of ^v for a group of anchors is found to decrease invariably with an increase in $ for a given value of S/B. For S > 2c/tan<j) , the uplift resistance of anchors in the group becomes equal to that of an isolated anchor, and no interference is seen to exist; where d is the depth of anchor. By examining the nodal velocity patterns, it was noted that in the event of collapse, a wedge of soil mass just above the anchors and encompassed within linear rupture surfaces moves vertically upward almost as a single rigid unit with the velocity same as that of the anchor plate itself. On this basis, a closed form solution of the problem has been developed. The results from the closed form solution for the group of two anchors as well as for multiple anchors are found to provide an excellent comparison with the rigorous upper bound numerical solution especially for the value of § greater than or equal to about 35°. For all the problems taken in this study, it has been seen that an upper bound limit analysis in combination with finite elements and linear programming is a very useful numerical tool for determining the magnitudes of collapse loads.