Effect of microstructure on stress–strain and pore-pressure response of sabarmati sand under the influence of mica

ABSTRACT

Micaceous soil is believed to be detrimental for civil engineering constructions due to the effect of high compressibility, low compacted density and low shear strength. Individual mica particle has numerous intact mica flakes foliated over each other making it flexible upon loading and rebound upon unloading due to its low hardness and resilient nature. Hence, micaceous soils with mica content more than 10% are considered undesirable for highway pavements, embankments and railway track constructions. When platy mica particles are sufficiently numerous to interact with spherical sand particles, bridging and ordering phenomena are augmented within the soil mass creating unique sand-mica particle orientation (MS microstructure) unlike sand-sand particle orientation (PS microstructure). The current experimental research was conducted to evaluate the variation in stress–strain, pore pressure and effective stress path response of Sabarmati sand under the influence of mica (sand with 30% mica and pure sand) with MS and PS microstructure respectively. Effect of particle crushing on stress–strain and pore pressure response was also studied on Sabarmati sand with MS and PS microstructure. Distinctive macroscopic response was observed in Sabarmati sand with MS microstructure under the influence of mica as well as mica particle crushing.     

Runoff quality prediction from small urban catchments using SWMM

The RUNOFF block of EPA's storm water management model (SWMM) was used to simulate the quantity and quality of urban storm water runoff from four relatively small sites (i.e. 5·97–23·56 ha) in South Florida, each with a specific predominant land use (i.e. low density residential, high density residential, highway and commercial). The objectives of the study were to test the applicability of this model in small subtropical urban catchments and provide modellers with a way to select appropriate input parameters to be used in planning studies. A total of 58 storm events, measured by the US Geological Survey (USGS), provided hyetographs, hydrographs and pollutant loadings for biological oxygen demand (BOD₅), total suspended solids (TSS), total Kjeldahl nitrogen (TKN) and lead (Pb), and were used for calibration of the model. Several other catchment characteristics, also measured or estimated by USGS, were used in model input preparation. Application of the model was done using the Green–Ampt equation for infiltration loss computation, a pollutant accumulation equation using a power build-up equation dependent on the number of dry days, and a power wash-off equation dependent on the predicted runoff rate. Calibrated quantity input parameters are presented and compared with suggested values in the literature. The impervious depression storage was generally found to be the most sensitive calibration parameter, followed by the Manning's roughness coefficients of conduit and overland flow, the Green–Ampt infiltration parameters and, finally, the pervious depression storage. Calibrated quality input parameters are presented in the form of regression equations, as a function of rainfall depth and the number of antecedent dry days. A total of 16 independent rainfall events were used for verification of the model, which showed a good comparison with observed data for both hydrographs and pollutant loadings. Average model predictions for the four constituent concentrations from the verification runs also showed good agreement with NURP published values in Florida and US sites. © 1998 John Wiley & Sons, Ltd.     

Development and application of a distributed catchment-scale hydrological model for the River Ythan,NE Scotland

Mathematical models are being used to develop a decision support system for integrated management of the Ythan catchment in NE Scotland. One component of this has involved the development of a distributed catchment-scale hydrological model. The model is based on subsurface flow routing and calculates the contribution to stream flow from each 50 m×50 m cell in the 548 km² catchment. It uses two topographic parameters, slope and distance to stream following the main line of flow, and five physical parameters. The topographic analysis and distributed flow accumulation are performed by linking the single cell model with a geographic information system. Preliminary results from a three-year simulation of daily flows indicate that the model successfully predicts the main characteristics of the catchment flow. © 1998 John Wiley & Sons, Ltd.     

Interaction between lateral sorting in river bends and vertical sorting in dunes

Sediment is sorted in river bends under the influence of gravity that pulls the heavier grains downslope and secondary flow that drags the finer grains upslope. Furthermore, when dunes are present, sediment is also sorted vertically at the dune lee side. However, sorting functions are poorly defined, since the relation to transverse bed slope and the interaction between lateral and vertical sorting is not yet understood for lack of data under controlled conditions. The objective of this study is to describe lateral sorting as a function of transverse bed slope and to gain an understanding of the interaction between lateral and vertical sorting in river bends. To this end, experiments were conducted with a poorly sorted sediment mixture in a rotating annular flume in which secondary flow intensity can be controlled separately from the main flow velocity, and therefore transverse bed slope towards the inner bend and dune dimensions can be systematically varied. Sediment samples were taken along cross-sections at the surface of dune troughs and dune crests, and over the entire depth at the location of dune crests (bulk samples), which enabled comparison of the relative contribution of vertical sorting by dunes to lateral sorting by the transverse bed slope. The data show that lateral sorting is always the dominant sorting mechanism in bends, and bulk samples showed minor effects of vertical sorting by dunes as long as all grain-size fractions are mobile. An empirical bend sorting model was fitted that redistributes the available sediment fractions over the cross-section as a function of transverse bed slope. Comparison with field data showed that the model accurately reproduces spatially-averaged trends in sorting at the bend apex in single-thread channels. The bend sorting model therefore provides a better definition of bend sorting with conservation of mass by size fraction and adds to current understanding of bend sorting. The implication for numerical modelling is that bend sorting mechanisms can be modelled independently of dunes, allowing the application of the active layer concept.     

Thermo-hydro-mechanical modeling of unsaturated soils using isogeometric analysis: Model development and application to strain localization simulation

This study presents a thermo-hydro-mechanical (THM) model of unsaturated soils using isogeometric analysis (IGA). The framework employs Bézier extraction to connect IGA to the conventional finite element analysis (FEA), featuring the current study as one of the first attempts to develop an IGA-FEA framework for solving THM problems in unsaturated soils. IGA offers higher levels of interelement continuity making it an attractive method for solving highly nonlinear problems. The governing equations of linear momentum, mass, and energy balance are coupled based on the averaging procedure within the hybrid mixture theory. The Drucker-Prager yield surface is used to limit the modified effective stress where the model follows small strain, quasi-static loading conditions. Temperature dependency of the surface tension is implemented in the soil-water retention curve. Nonuniform rational B-splines (NURBS) basis functions are used in the standard Galerkin method and weak formulations of the balance equations. Displacement, capillary pressure, gas pressure, and temperature are four independent quantities that are approximated by NURBS in spatial discretization. The framework is used to simulate strain localization in an undrained dense sand subjected to plane strain biaxial compression under different temperatures and displacement velocities. Results show that an increase in the displacement rate leads to reduction in the equivalent plastic strain while an increase in the temperature leads to an increase in the equivalent plastic strain. The findings suggest that the proposed IGA-based framework offers a viable alternative for solving THM problems in unsaturated soils.     

An examination of the validation of a model of the hydro/thermo/mechanical behaviour of engineered clay barriers

This paper focuses attention on the development of a numerical model of the hydro/thermo/mechanical behaviour of unsaturated clay and its consequent verification and validation. The work presented describes on-going collaboration between the Cardiff School of Engineering and Atomic Energy of Canada. The model development, which was carried out at Cardiff, can be described as being based on a mechanistic approach to coupled heat, moisture and air flow. This is then linked to a deformation analysis of the material within a ‘consolidation’ type of model. The whole is solved via the finite element method to yield a computer software code named COMPASS (COde for Modelling PArtly Saturated Soil). Some aspects of verification and validation of the model have been addressed in-house. However, the purpose of current AECL work is to provide an independent, rigorous, structured programme of validation and the paper will also explore the further validation of COMPASS within this context. © 1998 by John Wiley & Sons, Ltd.     

CASM: a unified state parameter model for clay and sand

The purpose of this paper is to present a simple, unified critical state constitutive model for both clay and sand. The model, called CASM (Clay And Sand Model), is formulated in terms of the state parameter that is defined as the vertical distance between current state (v, p′) and the critical state line in v–ln p′ space. The paper first shows that the standard Cam-clay models (i.e. the original and modified Cam-clay models) can be reformulated in terms of the state parameter. Although the standard Cam-clay models prove to be successful in modelling normally consolidated clays, it is well known that they cannot predict many important features of the behavior of sands and overconsolidated clays. By adopting a general stress ratio-state parameter relation to describe the state boundary surface of soils, it is shown that a simple, unified constitutive model (CASM) can be developed for both clay and sand. It is also demonstrated that the standard Cam-clay yield surfaces can be either recovered or approximated as special cases of the yield locus assumed in CASM. The main feature of the proposed model is that a single set of yield and plastic potential functions has been used to model the behaviour of clay and sand under both drained and undrained loading conditions. In addition, it is shown that the behaviour of overconsolidated clays can also be satisfactorily modelled. Simplicity is a major advantage of the present state parameter model, as only two new material constants need to be introduced when compared with the standard Cam-clay models. © 1998 John Wiley & Sons, Ltd.     

GEOMORPHOLOGICAL TRANSFER FUNCTION CALCULATED FROM DIGITAL ELEVATION MODELS FOR DISTRIBUTED HYDROLOGICAL MODELLING

Recently, several attempts have been made to relate the hydrological response of a catchment to its morphological and topographical features using different hypotheses to model the effect of the drainage network. Several transfer functions were developed and some of these are based on the theory of a linear model, the geomorphological unit hydrograph. The aim of this paper is to present a methodology to automatically identify the transfer function, using digital elevation models for applications in distributed hydrological modelling. The transfer function proposed herein is based on the Hayami approximation solution of the diffusive wave equation especially adapted for the routing hydrograph through a channel network. The Gardon d’Anduze basin, southern France, was retained for applications. Digital elevation models were used to extract the channel network and divide the basin into subcatchments. Each subcatchment produces, at its own outlet, an impulse response which is routed to the outlet of the whole catchment using the diffusive wave model described by two parameters: celerity and diffusivity functions of geometrical characteristics of the channel network. Firstly, a geomorphological unit hydrograph obtained by routing a homogeneous effective rainfall was compared with the unit hydrograph identified by a lumped model scheme, then the distributed model was applied to take into account the spatial variability of effective rainfall in the catchment. Results show that this new method seems to be adapted for distributed hydrological modelling; it enables identification of a transfer function response for each hydrological unit, here subcatchments, and then simulation of the contribution of each unit to the hydrograph at the outlet. © 1997 by John Wiley & Sons, Ltd.