Micromechanical Modelling of Granular Materials: Effect of Particle Size and Gradation

In this paper the effects of maximum particle size, particle gradation/sorting and fabric on bulk mechanical behaviour of granular materials such as coarse grained soils and rockfills are investigated" from micromechanical considerations starting from the grain scale level, using numerical" simulations based on Discrete Element Modelling (DEM). Hydrostatic compaction and biaxial tests on 2-dimensional assemblies of discs with varying particle sizes and gradations were modelled using DEM. An examination of the constitutive behaviour of granular media considering" the particulate nature of the medium has been attempted to explain the effect of particle size and gradation. Simulation results on perfectly parallel graded assemblies indicate that with increase in the size of the particles, a marginal increase (or no increase) in the angle of internal friction is observed during biaxial loading conditions. A change to a wider gradation (keeping the minimum grain size the same) results in a decrease in the angle of internal friction and an increase in volumetric strain to a considerable extent. Based on micromechanical force and fabric parameters, the basis for the physical behaviour was established. This helps in understanding the physics of parallel gradation techniques.     

Evaluation of spatial variation of peak horizontal acceleration and spectral acceleration for south India: a probabilistic approach

In this work, an attempt has been made to evaluate the spatial variation of peak horizontal acceleration (PHA) and spectral acceleration (SA) values at rock level for south India based on the probabilistic seismic hazard analysis (PSHA). These values were estimated by considering the uncertainties involved in magnitude, hypocentral distance and attenuation of seismic waves. Different models were used for the hazard evaluation, and they were combined together using a logic tree approach. For evaluating the seismic hazard, the study area was divided into small grids of size 0.1° × 0.1°, and the hazard parameters were calculated at the centre of each of these grid cells by considering all the seismic sources within a radius of 300 km. Rock level PHA values and SA at 1 s corresponding to 10% probability of exceedance in 50 years were evaluated for all the grid points. Maps showing the spatial variation of rock level PHA values and SA at 1 s for the entire south India are presented in this paper. To compare the seismic hazard for some of the important cities, the seismic hazard curves and the uniform hazard response spectrum (UHRS) at rock level with 10% probability of exceedance in 50 years are also presented in this work.     

Non-linear analysis of geomechanical problems using coupled finite and infinite elements

In modeling of many geomechanics problems such as underground openings, soil-foundation structure interaction problems, and in wave propagation problems through semi-infinite soil medium the soil is represented as a region of either infinite or semi-infinite extent. Numerical modeling of such problems using conventional finite elements involves a truncation of the far field in which the infinite boundary is terminated at a finite distance. In these problems, appropriate boundary conditions are introduced to approximate the solution of the infinite or semi-infinite boundaries as closely as possible. However, the task of positioning the finite boundary in conventional finite element discretization and the definition of the boundary and its conditions is very delicate and depends on the modeller's skill and intuition. Moreover, such a choice is influenced by the size of the domain to be discretized. Consequently, the dimensions of the global matrices and the time required for solution of the problem will increase considerably and also selection of the arbitrary location of truncated boundary may lead to erroneous result. In order to over come these problems, mapped infinite elements have been developed by earlier researchers (Simoni and Schrefier, 1987). In the present work the applicability of infinite element technique is examined for different geomechanics problems. A computer program INFEMEP is developed based on the conventional finite element and mapped infinite element technique. It is then validated using selected problems such as strip footing and circular footing. CPU time taken to obtain solutions using finite element approach and infinite element approach was estimated and presented to show the capability of coupled modeling in improving the computational efficiency. Mesh configurations of different sizes were used to explore the enhancement of both computational economy and solution accuracy achieved by incorporation of infinite elements to solve elastic and elasto-plastic problems in semi-infinite/finite domain as applied to geotechnical engineering. © Rapid Science Ltd. 1998     

Seismic microzonation of a nuclear power plant site with detailed geotechnical,geophysical and site effect studies

This paper highlights the seismic microzonation carried out for a nuclear power plant site. Nuclear power plants are considered to be one of the most important and critical structures designed to withstand all natural disasters. Seismic microzonation is a process of demarcating a region into individual areas having different levels of various seismic hazards. This will help in identifying regions having high seismic hazard which is vital for engineering design and land-use planning. The main objective of this paper is to carry out the seismic microzonation of a nuclear power plant site situated in the east coast of South India, based on the spatial distribution of the hazard index value. The hazard index represents the consolidated effect of all major earthquake hazards and hazard influencing parameters. The present work will provide new directions for assessing the seismic hazards of new power plant sites in the country. Major seismic hazards considered for the evaluation of the hazard index are (1) intensity of ground shaking at bedrock, (2) site amplification, (3) liquefaction potential and (4) the predominant frequency of the earthquake motion at the surface. The intensity of ground shaking in terms of peak horizontal acceleration (PHA) was estimated for the study area using both deterministic and probabilistic approaches with logic tree methodology. The site characterization of the study area has been carried out using the multichannel analysis of surface waves test and available borehole data. One-dimensional ground response analysis was carried out at major locations within the study area for evaluating PHA and spectral accelerations at the ground surface. Based on the standard penetration test data, deterministic as well as probabilistic liquefaction hazard analysis has been carried out for the entire study area. Finally, all the major earthquake hazards estimated above, and other significant parameters representing local geology were integrated using the analytic hierarchy process and hazard index map for the study area was prepared. Maps showing the spatial variation of seismic hazards (intensity of ground shaking, liquefaction potential and predominant frequency) and hazard index are presented in this work.     

Assessment of Seismically Induced Landslide Hazard for the State of Karnataka Using GIS Technique

Landslide hazards are a major natural disaster that affects most of the hilly regions around the world. In India, significant damages due to earthquake induced landslides have been reported in the Himalayan region and also in the Western Ghat region. Thus there is a requirement of a quantitative macro-level landslide hazard assessment within the Indian subcontinent in order to identify the regions with high hazard. In the present study, the seismic landslide hazard for the entire state of Karnataka, India was assessed using topographic slope map, derived from the Digital Elevation Model (DEM) data. The available ASTER DEM data, resampled to 50 m resolution, was used for deriving the slope map of the entire state. Considering linear source model, deterministic seismic hazard analysis was carried out to estimate peak horizontal acceleration (PHA) at bedrock, for each of the grid points having terrain angle 10° and above. The surface level PHA was estimated using nonlinear site amplification technique, considering B-type NEHRP site class. Based on the surface level PHA and slope angle, the seismic landslide hazard for each grid point was estimated in terms of the static factor of safety required to resist landslide, using Newmark’s analysis. The analysis was carried out at the district level and the landslide hazard map for all the districts in the Karnataka state was developed first. These were then merged together to obtain a quantitative seismic landslide hazard map of the entire state of Karnataka. Spatial variations in the landslide hazard for all districts as well as for the entire state Karnataka is presented in this paper. The present study shows that the Western Ghat region of the Karnataka state is found to have high landslide hazard where the static factor of safety required to resist landslide is very high.     

Undrained Cyclic Pore Pressure Response of Sand–Silt Mixtures: Effect of Nonplastic Fines and Other Parameters

Literature regarding the pore pressure generation characteristics and in turn the cyclic resistance behaviour of silty sand deposits is confusing. In an attempt to clarify the effect of nonplastic fines on undrained cyclic pore pressure response of sand–silt mixtures, an experimental programme utilising around 289 stress-controlled cyclic triaxial tests on specimens of size 50 mm diameter and 100 mm height was carried out at a frequency of 0.1 Hz. Specimens were prepared to various measures of density through constant gross void ratio approach, constant relative density approach, constant sand skeleton void ratio approach, and constant interfine void ratio approach to study the effect of nonplastic fines on pore pressure response of sand–silt mixtures. The effect of relative density, confining pressure as well as the frequency and magnitude of cyclic loading was also studied. It was observed that the pore pressure response is greatly influenced by the limiting silt content and the relative density of a specimen corresponding to any approach. The influence of other parameters such as relative density, confining pressure and magnitude of cyclic loading was as usual but an increase in frequency of cyclic loading was seen to generate excess pore pressure at a higher rate indicating an impact load type of behaviour at higher frequency. Utilising the entire test results over a wide range of parameters a new pore pressure band for sand–silt mixtures in line with Lee and Albaisa (1974) has been proposed. Similarly another pore pressure band corresponding to 10th cycle of loading as suggested by Dobry (1985) and up to a shear strain of around 25% has been proposed. These two bands can readily be used by researchers and field engineers to readily assess the pore pressure response of sand–silt mixtures.     

Spatial Variability of the Depth of Weathered and Engineering Bedrock using Multichannel Analysis of Surface Wave Method

In this paper an attempt has been made to evaluate the spatial variability of the depth of weathered and engineering bedrock in Bangalore, south India using Multichannel Analysis of Surface Wave (MASW) survey. One-dimensional MASW survey has been carried out at 58 locations and shear-wave velocities are measured. Using velocity profiles, the depth of weathered rock and engineering rock surface levels has been determined. Based on the literature, shear-wave velocity of 330 ± 30 m/s for weathered rock or soft rock and 760 ± 60 m/s for engineering rock or hard rock has been considered. Depths corresponding to these velocity ranges are evaluated with respect to ground contour levels and top surface levels have been mapped with an interpolation technique using natural neighborhood. The depth of weathered rock varies from 1 m to about 21 m. In 58 testing locations, only 42 locations reached the depths which have a shear-wave velocity of more than 760 ± 60 m/s. The depth of engineering rock is evaluated from these data and it varies from 1 m to about 50 m. Further, these rock depths have been compared with a subsurface profile obtained from a two-dimensional (2-D) MASW survey at 20 locations and a few selected available bore logs from the deep geotechnical boreholes.     

Effect of infill materials on the performance of geocell reinforced soft clay beds

The present study deals with model plate load tests conducted on geocell reinforced soft clay beds to evaluate the effect of infill materials on the performance of the geocell. Commercially available Neoweb geocells are used in the study. Three different infill materials namely aggregate, sand and local red soil were used in the study. The load carrying capacity of the geocell reinforced bed (as compared to an unreinforced bed) was found to be increased by 13 times for the aggregate infill, 11 times for the sand infill and 10 times for the red soil infill. Similarly the reduction in the settlement was in the order of 78%, 73% and 70% aggregate, sand and the red soil infill materials respectively. Results suggest that the performance of the geocell was not heavily influenced by the infill materials. Further, numerical simulations were carried out using FLAC^2D to validate the experimental findings. The results from numerical studies are in reasonably good agreement with the experimental findings. The outcome of this work is successfully implemented in the construction of the geocell foundation to support a 3 m high embankment in the settled red mud in Lanjighar (Orissa) in India.     

Analysis of strength and moduli of jointed rocks

This paper deals with two aspects of jointed rock mass behavior, first the finite element modeling of a jointed rock mass as an equivalent continuum, second the comparison of empirical strength criteria of a jointed rock mass. In finite element modeling the jointed rock properties are represented by a set of empirical relationships, which express the properties of the jointed medium as a function of joint factor and the properties of the intact rock. These relationships have been derived from a large set of experimental data of tangent elastic modulus. It is concluded that equivalent continuum analysis gives the best results for both single and multiple jointed rock. The reliability of the analysis depends on the estimation of joint factor, which is a function of the joint orientation, joint frequency and joint strength.Empirical strength criteria for jointed rocks, namely Hoek and Brown, Yudhbir et al., Ramamurthy and Arora, Mohr–Coulomb have been incorporated in a nonlinear finite element analysis of jointed rock using the equivalent continuum approach, to determine the failure stress. The major principal stress at failure, obtained using Ramamurthy's criteria, compares very well with experimental results.