Site characterization is an important task in Geotechnical Engineering. In situ tests based on standard penetration test (SPT), cone penetration test (CPT) and shear wave velocity survey are popular among geotechnical engineers. Site characterization using any of these properties based on finite number of in-situ test data is an imperative task in probabilistic site characterization. These methods have been used to design future soil sampling programs for the site and to specify the soil stratification. It is never possible to know the geotechnical properties at every location beneath an actual site because, in order to do so, one would need to sample and/or test the entire subsurface profile. Therefore, the main objective of site characterization models is to predict the subsurface soil properties with minimum in-situ test data. The prediction of soil property is a difficult task due to the uncertainities. Spatial variability, measurement ‘noise’, measurement and model bias, and statistical error due to limited measurements are the sources of uncertainities. Liquefaction in soil is one of the other major problems in geotechnical earthquake engineering. It is defined as the transformation of a granular material from a solid to a liquefied state as a consequence of increased pore-water pressure and reduced effective stress. The generation of excess pore pressure under undrained loading conditions is a hallmark of all liquefaction phenomena. This phenomena was brought to the attention of engineers more so after Niigata(1964) and Alaska(1964) earthquakes. Liquefaction will cause building settlement or tipping, sand boils, ground cracks, landslides, dam instability, highway embankment failures, or other hazards. Such damages are generally of great concern to public safety and are of economic significance. Site-spefific evaluation of liquefaction susceptibility of sandy and silty soils is a first step in liquefaction hazard assessment. Many methods (intelligent models and simple methods as suggested by Seed and Idriss, 1971) have been suggested to evaluate liquefaction susceptibility based on the large data from the sites where soil has been liquefied / not liquefied. The rapid advance in information processing systems in recent decades directed engineering research towards the development of intelligent models that can model natural phenomena automatically. In intelligent model, a process of training is used to build up a model of the particular system, from which it is hoped to deduce responses of the system for situations that have yet to be observed. Intelligent models learn the input output relationship from the data itself. The quantity and quality of the data govern the performance of intelligent model. The objective of this study is to develop intelligent models [geostatistic, artificial neural network(ANN) and support vector machine(SVM)] to estimate corrected standard penetration test (SPT) value, Nc, in the three dimensional (3D) subsurface of Bangalore. The database consists of 766 boreholes spread over a 220 sq km area, with several SPT N values (uncorrected blow counts) in each of them. There are total 3015 N values in the 3D subsurface of Bangalore. To get the corrected blow counts, Nc, various corrections such as for overburden stress, size of borehole, type of sampler, hammer energy and length of connecting rod have been applied on the raw N values. Using a large database of Nc values in the 3D subsurface of Bangalore, three geostatistical models (simple kriging, ordinary kriging and disjunctive kriging) have been developed. Simple and ordinary kriging produces linear estimator whereas, disjunctive kriging produces nonlinear estimator. The knowledge of the semivariogram of the Nc data is used in the kriging theory to estimate the values at points in the subsurface of Bangalore where field measurements are not available. The capability of disjunctive kriging to be a nonlinear estimator and an estimator of the conditional probability is explored. A cross validation (Q1 and Q2) analysis is also done for the developed simple, ordinary and disjunctive kriging model. The result indicates that the performance of the disjunctive kriging model is better than simple as well as ordinary kriging model. This study also describes two ANN modelling techniques applied to predict Nc data at any point in the 3D subsurface of Bangalore. The first technique uses four layered feed-forward backpropagation (BP) model to approximate the function, Nc=f(x, y, z) where x, y, z are the coordinates of the 3D subsurface of Bangalore. The second technique uses generalized regression neural network (GRNN) that is trained with suitable spread(s) to approximate the function, Nc=f(x, y, z). In this BP model, the transfer function used in first and second hidden layer is tansig and logsig respectively. The logsig transfer function is used in the output layer. The maximum epoch has been set to 30000. A Levenberg-Marquardt algorithm has been used for BP model. The performance of the models obtained using both techniques is assessed in terms of prediction accuracy. BP ANN model outperforms GRNN model and all kriging models. SVM model, which is firmly based on the theory of statistical learning theory, uses regression technique by introducing -insensitive loss function has been also adopted to predict Nc data at any point in 3D subsurface of Bangalore. The SVM implements the structural risk minimization principle (SRMP), which has been shown to be superior to the more traditional empirical risk minimization principle (ERMP) employed by many of the other modelling techniques. The present study also highlights the capability of SVM over the developed geostatistic models (simple kriging, ordinary kriging and disjunctive kriging) and ANN models. Further in this thesis, Liquefaction susceptibility is evaluated from SPT, CPT and Vs data using BP-ANN and SVM. Intelligent models (based on ANN and SVM) are developed for prediction of liquefaction susceptibility using SPT data from the 1999 Chi-Chi earthquake, Taiwan. Two models (MODEL I and MODEL II) are developed. The SPT data from the work of Hwang and Yang (2001) has been used for this purpose. In MODEL I, cyclic stress ratio (CSR) and corrected SPT values (N1)60 have been used for prediction of liquefaction susceptibility. In MODEL II, only peak ground acceleration (PGA) and (N1)60 have been used for prediction of liquefaction susceptibility. Further, the generalization capability of the MODEL II has been examined using different case histories available globally (global SPT data) from the work of Goh (1994). This study also examines the capabilities of ANN and SVM to predict the liquefaction susceptibility of soils from CPT data obtained from the 1999 Chi-Chi earthquake, Taiwan. For determination of liquefaction susceptibility, both ANN and SVM use the classification technique. The CPT data has been taken from the work of Ku et al.(2004). In MODEL I, cone tip resistance (qc) and CSR values have been used for prediction of liquefaction susceptibility (using both ANN and SVM). In MODEL II, only PGA and qc have been used for prediction of liquefaction susceptibility. Further, developed MODEL II has been also applied to different case histories available globally (global CPT data) from the work of Goh (1996). Intelligent models (ANN and SVM) have been also adopted for liquefaction susceptibility prediction based on shear wave velocity (Vs). The Vs data has been collected from the work of Andrus and Stokoe (1997). The same procedures (as in SPT and CPT) have been applied for Vs also. SVM outperforms ANN model for all three models based on SPT, CPT and Vs data. CPT method gives better result than SPT and Vs for both ANN and SVM models. For CPT and SPT, two input parameters {PGA and qc or (N1)60} are sufficient input parameters to determine the liquefaction susceptibility using SVM model. In this study, an attempt has also been made to evaluate geotechnical site characterization by carrying out in situ tests using different in situ techniques such as CPT, SPT and multi channel analysis of surface wave (MASW) techniques. For this purpose a typical site was selected wherein a man made homogeneous embankment and as well natural ground has been met. For this typical site, in situ tests (SPT, CPT and MASW) have been carried out in different ground conditions and the obtained test results are compared. Three CPT continuous test profiles, fifty-four SPT tests and nine MASW test profiles with depth have been carried out for the selected site covering both homogeneous embankment and natural ground. Relationships have been developed between Vs, (N1)60 and qc values for this specific site. From the limited test results, it was found that there is a good correlation between qc and Vs. Liquefaction susceptibility is evaluated using the in situ test data from (N1)60, qc and Vs using ANN and SVM models. It has been shown to compare well with “Idriss and Boulanger, 2004” approach based on SPT test data. SVM model has been also adopted to determine over consolidation ratio (OCR) based on piezocone data. Sensitivity analysis has been performed to investigate the relative importance of each of the input parameters. SVM model outperforms all the available methods for OCR prediction.

In recent years, the increasing impacts of landslide hazards on human lives and lifeline facilities worldwide has advanced the necessity to find out both economically acceptable and useful techniques to predict the occurrence and destructive power of landslides. Though many projects exist to attain this goal, the current investigation set out to establish an understanding of the initiation and propagation mechanisms of landslides via numerical simulations, so that mitigation strategies to reduce the long-term losses from landslide hazards can be made. In this research, the Discrete Element Method (DEM) and Computational Fluid Dynamics (CFD) have been used to investigate the mechanical and hydraulic behaviour of granular materials involved in landslides. The main challenge is to provide rational analyses of large scale landslides via small scale numerical simulations. To solve this problem, dimensional analyses have been performed on a simple granular column collapse model. The influence of governing dimensionless groups on the debris runout distance and deposit height has been studied for the terrestrial and submerged granular flows. 3D DEM investigations of granular flows in plane strain conditions have been performed in this research. The input parameters of the DEM model have been calibrated by the numerical triaxial tests, based on which, the relationships between the microscopic variables and the macroscopic soil strength properties are analysed. Using the simple granular column collapse model, the influences of column aspect ratio, characteristic strain, model size ratio and material internal friction angle on the runout distance and deposit height of granular materials have been examined. Additionally, the deformation and energy evolution of dry granular materials are also discussed. The DEM-CFD coupling model has been employed to study the mechanical and hydraulic behaviour of highly mobilized terrestrial / submarine landslides. This model has been validated via numerical simulations of fluid flow through a porous soil sample and grain batch sedimentations. The simulations of granular flows in the submerged environment have led to some meaningful insights into the flow mechanisms, such as the mobilization of sediments, the generation and dissipation of excess pore water pressures and the evolution of effective stresses. Overall, this study shows that the proposed numerical tools are capable of modelling the mechanical and hydraulic behaviour of terrestrial and submarine landslides.

The Canterbury Earthquakes of 2010-2011, in particular the 4th September 2010 Darfield earthquake and the 22nd February 2011 Christchurch earthquake, produced severe and widespread liquefaction in Christchurch and surrounding areas. The scale of the liquefaction was unprecedented, and caused extensive damage to a variety of man-made structures, including residential houses. Around 20,000 residential houses suffered serious damage as a direct result of the effects of liquefaction, and this resulted in approximately 7000 houses in the worst-hit areas being abandoned. Despite the good performance of light timber-framed houses under the inertial loads of the earthquake, these structures could not withstand the large loads and deformations associated with liquefaction, resulting in significant damage. The key structural component of houses subjected to liquefaction effects was found to be their foundations, as these are in direct contact with the ground. The performance of house foundations directly influenced the performance of the structure as a whole. Because of this, and due to the lack of research in this area, it was decided to investigate the performance of houses and in particular their foundations when subjected to the effects of liquefaction. The data from the inspections of approximately 500 houses conducted by a University of Canterbury summer research team following the 4th September 2010 earthquake in the worst-hit areas of Christchurch were analysed to determine the general performance of residential houses when subjected to high liquefaction loads. This was followed by the detailed inspection of around 170 houses with four different foundation types common to Christchurch and New Zealand: Concrete perimeter with short piers constructed to NZS3604, concrete slab-on-grade also to NZS3604, RibRaft slabs designed by Firth Industries and driven pile foundations. With a focus on foundations, floor levels and slopes were measured, and the damage to all areas of the house and property were recorded. Seven invasive inspections were also conducted on houses being demolished, to examine in more detail the deformation modes and the causes of damage in severely affected houses. The simplified modelling of concrete perimeter sections subjected to a variety of liquefaction-related scenarios was also performed, to examine the comparative performance of foundations built in different periods, and the loads generated under various bearing loss and lateral spreading cases. It was found that the level of foundation damage is directly related to the level of liquefaction experienced, and that foundation damage and liquefaction severity in turn influence the performance of the superstructure. Concrete perimeter foundations were found to have performed most poorly, suffering high local floor slopes and being likely to require foundation repairs even when liquefaction was low enough that no surface ejecta was seen. This was due to their weak, flexible foundation structure, which cannot withstand liquefaction loads without deforming. The vulnerability of concrete perimeter foundations was confirmed through modelling. Slab-on-grade foundations performed better, and were unlikely to require repairs at low levels of liquefaction. Ribraft and piled foundations performed the best, with repairs unlikely up to moderate levels of liquefaction. However, all foundation types were susceptible to significant damage at higher levels of liquefaction, with maximum differential settlements of 474mm, 202mm, 182mm and 250mm found for concrete perimeter, slab-on-grade, ribraft and piled foundations respectively when subjected to significant lateral spreading, the most severe loading scenario caused by liquefaction. It was found through the analysis of the data that the type of exterior wall cladding, either heavy or light, and the number of storeys, did not affect the performance of foundations. This was also shown through modelling for concrete perimeter foundations, and is due to the increased foundation strengths provided for heavily cladded and two-storey houses. Heavy roof claddings were found to increase the demands on foundations, worsening their performance. Pre-1930 concrete perimeter foundations were also found to be very vulnerable to damage under liquefaction loads, due to their weak and brittle construction.

The distinct element method (DEM) is a useful tool in rock engineering to model jointed rock masses. To simulate a jointed rock mass realistically, the main challenge is to be able to capture its complex geometry which consists of blocks with various shapes and sizes, and to model the interactions between these blocks. The main contribution of this thesis is the development of novel algorithms in the DEM to model jointed rock masses, namely rock slicing procedures for block generation, and algorithms for contact detection between polygonal blocks in 2-D or polyhedral blocks in 3-D. These algorithms make use of convex optimisation techniques, for which there exist efficient solution procedures. They do not rely on conventional vertex-edge-face hierarchical data structures and tedious housekeeping algorithms. The algorithms have been verified against analytical and numerical solutions, as well as validated against experimental results published in the literature. Among those, the results of DEM simulations were compared against the experimental model tests and numerical simulations of jointed beams carried out by Talesnick et al. (2007) and Tsesarsky & Talesnick (2007) respectively. Emphasis was placed on modelling the stiffness of the block interfaces accurately, and this was accomplished by reinterpreting the laboratory data published by the investigators. The capabilities of the numerical tools are also examined and demonstrated in areas for which the DEM has found practical application. A substantial fraction of this thesis is devoted to illustrating how these tools can assist the engineer in designing support systems; for example, designing the length and spacing of rock bolts and the lining thickness for a tunnel. Algorithms to model rock bolt and lining support were implemented for this purpose. Interesting comparisons with elastic solutions for supported openings were obtained. Further, it is shown that the relative benefit of introducing more rock bolts or thicker lining can be evaluated using the numerical tools with the aid of an interaction diagram. In the final part of this thesis, the case history of the 1963 Vaiont rock slide in Italy is studied. The 2-D analyses led to useful insights concerning the influence of the reservoir water level, the rock mass strength and deformability, and the slide surface shear stiffness. 3-D analyses were undertaken to investigate the influence of the eastern boundary of the slope, and interesting insights were obtained concerning the slope kinematics. Overall, the case study shows that the tools are capable of modelling problems with specific physical and geometrical detail in both 2-D and 3-D.

Vibratory driven sheet piles are a cost-effective retaining wall structure, and in coming decades the continued use of this method will be crucial for minimising costs within the construction sector. However, vibratory driven sheet piles are a source of ground vibrations, which may harm structures or induce disturbance. Most urban construction projects face strict limits on permissible vibration level. Being able to reliably predict the expected vibration level prior to construction is therefore highly important. Reliable prediction demands a profound knowledge of the vibration transfer process, from source to point of interest. This thesis focuses on clarifying the vibration transfer process and will serve as a platform for the future development of a reliable prediction model. The vibration transfer process is divided into two main parts: vibration source and vibrations in soil. The different parts in the vibration transfer process are studied and investigated with the help of a literature review, field tests and numerical modelling. Within the scope of this thesis, three field tests have been conducted and a new instrumentation system has been developed. The new instrumentation system enables recording of both sheet pile vibrations and ground vibrations at depth during the entire driving. The field tests aimed to study the vibration transfer from sheet pile to soil and the vibration transfer within a sheet pile wall, as well as the wave pattern in soil. To study sheet pile behaviour during driving a numerical model was developed, which is also meant to serve as a basis for further studies. The main scientific contribution of this thesis is the identification of the sheet pile behaviour during driving. For practical application, the main contribution is the development of an increased knowledge of the vibration transfer process from source to soil, together with the new instrumentation system and the development of the numerical model. / Vibrodriven spont är en kostnadseffektiv stödkonstruktion och i framtiden kommer den fortsatta användningen av denna metod att vara nödvändig för att minimera kostnader för byggprojekt. Vibrodriven spont är dock en källa till markvibrationer, som kan skada byggnader eller orsaka störningar. De flesta byggprojekt måste förhålla sig till strikta krav gällande vibrationsnivåer. Möjligheten att på ett tillförlitligt sätt förutsäga vibrationsnivåerna innan bygget startar är därför av största vikt. Tillförlitlig prognos av vibrationsnivåer i samband med vibrodrivning av spont kräver god kännedom om vibrationsöverföringsprocessen, från källan till det potentiella skadeobjektet. Denna avhandling fokuserar på att förtydliga vibrationsöverföringsprocessen och fungera som en plattform för framtida utveckling av en tillförlitlig prognosmodell. Vibrationsöverföringsprocessen delas in i två huvuddelar; vibrationskällan och vibrationer i jord. De olika delarna av vibrationsöverföringsprocessen studeras och undersöks med hjälp av litteraturstudie, fältförsök och numerisk modellering. Inom ramarna för denna avhandling har tre fältförsök utförts och ett nytt instrumenteringssystem har utvecklats. Det nya instrumenteringssystemet möjliggör mätning av både spontvibrationer och vibrationer på djup i jorden, under hela neddrivningsfasen. Fältförsöken syftade till att studera vibrationsöverföringen mellan spont och jord, vibrationsöverföringen inom en spontvägg samt vågmönstret i jorden under drivning. För att studera spontens beteende under neddrivning utvecklades en numerisk modell, som också kan fungera som en bas för framtida studier. Avhandlingens huvudsakliga vetenskapliga bidrag är identifieringen av spontens beteende under neddrivning. För praktisk tillämpning är det huvudsakliga bidraget förklaringen av vibrationsöverföringsprocessen från källa till jord, det nya instrumenteringssystemet samt utvecklingen av den numeriska modellen.

Correct modelling of embankment breach formation is essential for an accurate assessment of the associated flood risk. Modelling breach formation due to overflow requires a thorough understanding of the geotechnical processes in unsaturated soils as well as erosion processes under supercritical flow conditions. This thesis describes 1D slope stability analysis performed for unsaturated soils whose moisture content changes with time. The analysis performed shows that sediment-laden gravity flows play an important role in the erosion behaviour of embankments. The thesis also describes a practical, fast breach model based on a simplified description of the physical processes that can be used in modelling and decision support frameworks for flooding. To predict the breach hydrograph, the rapid model distinguishes between breach formation due to headcut erosion and surface erosion in the case of failure due to overflow. The model also predicts the breach hydrograph in the case of failure due to piping. The assumptions with respect to breach flow modelling are reviewed, and result in a new set of breadth-integrated Navier-Stokes equations, that account for wall shear stresses and a variable breadth geometry. The vertical 2D flow field described by the equations can be used to calculate accurately the stresses on the embankment during the early stages of breach formation. Pressure-correction methods are given for solving the 2D Navier-Stokes equations for a variable breadth, and good agreement is found when validating the flow model against analytical solutions.

Microbial Induced Calcite Precipitation (MICP) is an attractive alternative for a variety geotechnical ground improvement practices commonly used today and has a variety of potential applications. This research focuses primarily on its use as a soil stabilization technique using the bacteria Sporosarcina Pasteurii and a single injection point percolation method adapted from previous research in granular soils. This method, and most published data, show an inherent variability in both physical and engineering properties due to the distribution of precipitated calcite within the specimen. The focus of this research is on the quantification of the variability in shear strength parameters induced by MICP treatment in sand. Also, on the initial development of a new treatment method which aims to reduce this inherent variability and offer a more feasible option for field applications. The MICP treated soil columns were sampled at constant intervals from the injection point and then subject to direct shear testing (DST) and calcite distribution analysis. This analysis reiterates previously documented reduction in cementation as distance from injection point increases. The reduction in cementation results in reduced shear strength parameter improvements. This research also concluded a minimum of two percent mass of calcite per total mass of treated soil for significant strength improvements.

The United States Navy (USN) currently utilizes a Rapid Penetration Test (RPT) on both land and in water as the means to determine whether sufficient soil bearing capacity exists for piles in axial compression, prior to construction of the Elevated Causeway System (Modular) [ELCAS(M)] pile-supported pier system. The USN desires a replacement for the RPT because of issues with the method incorrectly classifying soils as well as the need to have a less labor-and-equipment-intensive method for geotechnical investigation. The Multichannel Analysis of Surface Waves (MASW) method is selected herein as the potential replacement for the RPT. The MASW method is an existing, geophysical method for determining soil properties based upon the acquisition and analysis of seismic surface waves used to develop shear wave velocity profiles for the soils at specific sites. Correlations between shear wave velocity and Cone Penetration Testing are utilized to classify soils, develop pile blow count estimates, and calculate soil bearing capacity. This researcher found that the MASW method was feasible and reliable in predicting the required properties for terrestrial sites. However, it was not successful in predicting those properties for underwater marine sites due to issues with equipment and field setup. Future areas of improvement are recommended to address these issues and, due to the success of the method on land, it is expected that once the issues are addressed the MASW method will be a reliable replacement for the RPT method across the entire subaerial and subaqueous profile.

Microbially induced calcite precipitation (MICP) has been used for a number of years as a technique for the improvement of various geological materials. MICP has been used in a limited capacity in organic rich soils with varying degrees of success. Investigators hypothesized that microbially-induced cementation could be improved in organic soils by using a surfactant. Varying amounts of Sodium Dodecyl Sulfate (SDS) were added to soils of varying organic content and a mixing procedure was used to treat these soils via MICP. Treated specimens were tested for unconfined compressive strength (UCS). Results appeared to show direct relationships between SDS content and treated specimen strength although significant variability was present in the data. In addition, results also indicated that while addition of SDS during MICP treatment strengthens soil, the strengthening is likely from the formation of a calcium dodecyl sulfate (CDS) complex in which the CDS surrounds the soil in a matrix, and formation of MICP-induced calcite has very little to do with overall soil performance. As such, a new method for stabilizing loose soils dubbed ‘Surfactant-induced soil stabilization’ (SISS) was further explored by treating additional soil specimens. Samples treated using this technique showed increases in strength when compared to untreated specimens. In addition, preliminary data indicated that SISS treated specimens were insoluble. The SISS technique presents a number of advantages when compared to traditional soil stabilization techniques. In particular it should be relatively low-cost and simple to administer since its only components are SDS and calcium chloride. Additionally, these constituents are relatively more sustainable than chemicals associated with more-traditional loose soil stabilization techniques.

Use of soil retaining walls for roads, embankments and bridges is increasing with time and reinforced soil retaining walls are found to be very efficient even under critical conditions compared to unreinforced walls. They offer competitive solutions to earth retaining problems associated with less space and more loads posed by tremendous growth in infrastructure, in addition to the advantages in ease and cost of construction compared to conventional retaining wall systems. The study of seismic performance of reinforced soil retaining walls is receiving much attention in the light of lessons learned from past failures of conventional retaining walls. Laboratory model studies on these walls under controlled seismic loading conditions help to understand better how these walls actually behave during earthquakes. The objective of the present study is to investigate the seismic response of geosynthetic reinforced soil wall models through shaking table tests. To achieve this, wrap faced and rigid faced reinforced soil retaining walls of size 750 × 500 mm in plan and 600 mm height are built in rigid and flexible containers and tested under controlled dynamic conditions using a uni-axial shaking table. The effects of frequency and acceleration of the base motion, surcharge pressure on the crest, number of reinforcing layers, container boundary, wall structure and reinforcement layout on the seismic performance of the retaining walls are studied through systematic series of shaking table tests. Results are analyzed to understand the effect of each of the considered parameters on the face displacements, acceleration amplifications and soil pressures on facing at different elevations of the walls. A numerical model is developed to simulate the shaking table tests on wrap faced reinforced soil walls using a computer program FLAC (Fast Lagrangian Analysis of Continua). The experimental data are used to validate the numerical model and parametric studies are carried out on 6 m height full-scale wall using this model. Thus, the study deals with the shaking table tests, dynamic response of reinforced walls and their numerical simulation. The thesis presents detailed description of various features and various parts of the shaking table facility along with the instrumentation and model containers. Methodology adopted for the construction of reinforced soil model walls and testing procedures are briefly described. Scaling and stability issues related to the model wall size and reinforcement strength are also discussed. From the study, it is observed that the displacements are decreasing with the increase in relative density of backfill, increase in surcharge pressure and increase in number of reinforcing layers; In general, accelerations are amplified to the most at the top of the wall; Behaviour of model walls is sensitive to model container boundary. The frequency content is very important parameter affecting the model response. Further, it is noticed that the face displacements are significantly affected by all of the above parameters, while the accelerations are less sensitive to reinforcement parameters. Even very low strength geonet and geotextile are able to reduce the displacements by 75% compared to unreinforced wall. The strain levels in the reinforcing elements are observed to be very low, in the order of ±150 micro strains. A random dynamic event is also used in one of the model tests and the resulted accelerations and displacements are presented. Numerical parametric studies provided important insight into the behaviour of wrap faced walls under various seismic loading conditions and variation in physical parameters.