Implementation of a Critical State Two-Surface Model to Evaluate the Response of Geosynthetic Reinforced Pavements

A finite-element model was developed using ABAQUS software package to investigate the effect of placing geosynthetic reinforcement within the base course layer on the response of a flexible pavement structure. A critical state two-surface constitutive model was first modified to represent the behavior of base course materials under the unsaturated field conditions. The modified model was then implemented into ABAQUS through a user defined subroutine, UMAT. The implemented model was validated using the results of laboratory triaxial tests. Finite-element analyses were then conducted on different unreinforced and geosynthetic reinforced flexible pavement sections. The results of this study demonstrated the ability of the modified critical state two-surface constitutive model to predict, with good accuracy, the response of the considered base course material at its optimum field conditions when subjected to cyclic as well as static loads. The results of the finite-element analyses showed that the geosynthetic reinforcement reduced the lateral strains within the base course and subgrade layers. Furthermore, the inclusion of the geosynthetic layer resulted in a significant reduction in the vertical and shear strains at the top of the subgrade layer. The improvement of the geosynthetic layer was found to be more pronounced in the development of the plastic strains rather than the resilient strains. The reinforcement benefits were enhanced as its elastic modulus increased.     

Shakedown Approaches to Rut Depth Prediction in Low-Volume Roads

Rutting, due to permanent deformations of unbound materials, is one of the principal damage modes of low traffic pavements. Flexible pavement design methods remain empirical; they do not take into account the inelastic behavior of pavement materials and do not predict the rutting under cyclic loading. A finite-element program, based on the concept of the shakedown theory developed by Zarka for metallic structures under cyclic loadings, has been used to estimate the permanent deformations of unbound granular materials subjected to traffic loading. Based on repeated load triaxial tests, a general procedure has been developed for the determination of the material parameters of the constitutive model. Finally, the results of a finite-element modeling of the long-term behavior of a flexible pavement with the simplified method are presented and compared to the results of a full-scale flexible pavement experiment performed by Laboratoire Central des Ponts et Chaussées. Finally, the calculation of the rut depth evolution with time is carried out.     

Expansive Soil Test Site Near Newcastle

A field site was established in 1993 near Newcastle, Australia, as part of a long-term study into expansive soil behavior. The primary objectives in establishing the site were to collect high quality data with which to check current design methods for lightly loaded building foundations and to develop improved understanding of the physical processes that drive unsaturated expansive soil behavior. The site was instrumented to allow measurement of soil water content, soil moisture suction, and ground movement to depths of 3 m. The site was provided with two ground covers to simulate moisture boundary conditions due to the presence of typical structures. This paper presents some of the more important findings from the 7 years of data acquired so far. These include a qualitative assessment of the overall site behavior, and quantitative observations of the range of total suction and water content changes with depth, the depth to which moisture changes occur, and the contributions to surface movement from ground movement at various depths. The shape of a mound developed beneath a flexible cover on an initially dry site is examined, and the effects of a large tree on moisture changes are reported.     

Analysis and Implementation of Resilient Modulus Models for Granular Solids

Constitutive equations based upon stress dependent moduli, like K-θ and Uzan-Witczak, are widely used to characterize the resilient response of granular materials for the analysis and design of pavement systems. These constitutive models are motivated by the observation that the granular layers used in pavement structures shake down to (nonlinear) elastic response under construction loads and will, therefore, respond elastically under service loads typically felt by these systems. Due to their simplicity, their great success in organizing the response data from cyclic triaxial tests, and their success relative to competing material models in predicting the behavior observed in the field, these resilient modulus constitutive models have been implemented in many computer programs used by researchers and design engineers. This paper provides an analysis of the nonlinear solution algorithms that have been used in implementing these models in a conventional nonlinear 3D finite-element framework. The analysis shows that these conventional algorithms are destined to fail at higher load levels. The paper offers two competitive methods for global analysis with these models. A comparative study of eight possible implementations of the algorithms described in the paper is made through two simulation examples.     

Unsaturated Constitutive Surfaces from Pressuremeter Tests

A methodology to identify the collapse potential of unsaturated soils is proposed in this paper on the basis of pressuremeter test results associated with independent measurements of the in situ matric suction. A solution combining the expansion of a cylindrical cavity to a modified Cam clay critical state model has been introduced and accommodated to the framework of unsaturated soil behavior. This accounts for changes in soil properties induced by suction changes. Interpretation of pressuremeter tests performed under unsaturated and soaked conditions links the amount of collapse to strength and stiffness changes and provides assessment to the constitutive soil parameters that are necessary to define the yield envelopes of the soil. A comprehensive site investigation program comprising field and laboratory tests carried out in two residual soil sites is discussed in order to validate the proposed methodology. Values of shear strength, in situ stress, and yield pressure derived from both field and laboratory data are used as input parameters of a constitutive model adopted for describing the yield envelopes of these unsaturated residual soil sites.     

Behavior of a Loosely Compacted Unsaturated Volcanic Soil

In this paper, the stress-strain relationship and volumetric behavior of a loosely compacted unsaturated decomposed volcanic soil (fill) were studied by conducting three series of triaxial stress path tests: (1) consolidated undrained on the saturated fill; (2) constant water content; and (3) a reducing suction under constant deviator stress on the unsaturated fill. The last two series of tests were designed to simulate the effects of undrained response and rainfall infiltration in initially unsaturated slopes, respectively. It was found that the saturated loose volcanic soil behaves like clay under isotropic compression but it resembles sand behavior when it was subjected to undrained shear. For isotropically consolidated unsaturated specimens sheared under a constant water content, a hardening stress-strain and a nonlinear shear strength-suction relationship are observed. At relatively high suctions, both angle of friction and apparent cohesion appear to be independent of suction. Volumetric contraction during shear is observed in this series of tests. On the other hand, anisotropically consolidated loose unsaturated specimens subjected to a reducing suction change from contractive to dilative behavior as the net mean stress increases. This observed volumetric behavior, unlike the shear strength, is stress path-dependent and cannot be explained by using the existing elastoplastic critical state theoretical framework extended for unsaturated soils.     

Volume Change Behavior of Collapsible Compacted Gneiss Soil

The mechanical behavior of a compacted metastable-structured residual gneiss soil was experimentally investigated. Volume changes were investigated using both a conventional oedometer cell and a triaxial permeameter system where the stress state variables were independently controlled. Wetting stress paths were utilized to reflect field conditions associated with collapsing earth structures. The compacted specimens were consolidated under both isotropic and k0-oedometer conditions. Measurements of total volume change and water content change were made at specified matric suction values following a wetting stress path under a given loading. The experimental data were analyzed to define volume change constitutive relationships for the metastable-structured soil. Best-fit modeling was used for predicting the volume change behavior of the compacted metastable-structured residual soil during wetting-induced collapse. The measured data are discussed from both phenomenological and microscopic viewpoints. Predictions of the best-fit models are compared to experimental results.     

Impact of Soil Type and Compaction Conditions on Soil Water Characteristic

Tests were conducted to determine the variation of water content and pore water suction for compacted clayey soils. The soils had varying amounts of clay fraction with plasticities ranging from low to high plasticity. The unsaturated soil behavior was investigated for six conditions, covering a range of compactive efforts and water contents. The experimental data were fit to four commonly used models for the water content-pore water suction relationship. Each model provided a satisfactory fit to the experimental data. However, the individual parameters obtained from the curve fits varied significantly between models. The soil water characteristic curves (SWCCs) were more sensitive to changes in compaction effort than changes in compaction water content. At similar water contents, the pore water suction increased with increasing compaction effort for each compaction condition and soil type. For all compaction conditions, the lowest plasticity soils retained the smallest water content and the highest plasticity soils retained the highest water content at a specified suction. In addition, SWCCs for soils compacted in the laboratory and in the field were similar.     

Degradation of a Granular Base under a Flexible Pavement: DEM Simulation

Flexible pavements are composed by an asphalt concrete layer, granular base and subbase layers, and a natural subgrade. The granular materials forming part of the granular layers are subjected to static and dynamic loads during their engineering life. As a result of these loads particle crushing may occur depending on the strength of the particles forming the granular layers. Particle crushing is important since it is associated with several detrimental effects such as settlements and a reduction in hydraulic conductivity. A computer simulation using the discrete element method (DEM) is presented in order to understand and visualize how crushing initiates and develops inside a simulated pavement structure.     

Impact of Bituminous Subballast on Railroad Track Deformation Considering Atmospheric Actions

High-speed ballasted track design standards usually require the use of sand and gravel layers as subballast to fulfill an accurate protection of the formation layers, not only against traffic loads, but also against the effects of weather. Seasonal changes in soil water content, or suction changes, are responsible for cyclic volumetric strains on railroad trackbed layers, thus on the infrastructure. Being almost completely water-resistant when compared with granular-only layers, bituminous subballast offers a higher protection of the subgrade, consequently improving its behavior along the infrastructure life cycle. This question is investigated through the comparison of the performance of the track formation against atmospheric actions, taking into consideration the unsaturated state of the geomaterials. The method adopted consists of modeling the vertical displacements of both bituminous and granular subballast designs through a finite-element coupled thermo-hydro-mechanical (THM) analysis. The comparison of the two design solutions confirms that the adoption of a bituminous subballast layer might allow important reductions in seasonal vertical displacements.     

Critical Review of the Methodologies Employed for Soil Suction Measurement

Modeling the behavior of unsaturated soils necessitates the measurement of soil suction and the establishment of its variation with the water content, which is commonly known as the soil-water characteristic curve (SWCC). Several methodologies have been developed for measuring either total suction ψ (sum of matric suction ψm and osmotic suction ψo) or ψm. While employing different methodologies for suction measurement, there is a possibility that various factors (viz., type of the soil, measurement methodology, range of the suction measurement, equilibration time, and presence of salts or contaminants in the soil) may influence the results and hence the SWCC. Therefore, it is essential to investigate the uniqueness of SWCC, determined by using some commonly adopted suction measurement methodologies. This study indicates that the SWCC established by adopting different methodologies may not be unique and is primarily influenced by the range of suction measurement. As such, it is essential to highlight the range of suction values involved for establishing the SWCC, to facilitate unambiguous modeling and to precisely understand the behavior of unsaturated soil.     

Triaxial Apparatus for Applying Liquid Infiltration with Controlled Boundary Conditions and Internal Suction Measurement

This note describes a triaxial apparatus that applies liquid infiltration under automatically controlled boundary conditions while using the new Xeritron sensor to measure suction. Three different infiltration boundary conditions of interest including constant mean stress (CMS), constant volume (CV), and constant stiffness (CS) are applied to examine the influence of boundary conditions on the mechanical and hydraulic behavior of unsaturated clay materials. CMS and CV tests provide limits of infiltration boundary conditions while CS tests represent a flexible spring-type boundary condition. Spatial distribution of gravimetric water content and bulk density are measured to calculate dry density and degree of saturation following testing to determine internal changes to the specimens. The results from initial testing show that the apparatus is providing experimental evidence to evaluate unsaturated flow models and elastic-plastic constitutive models to examine the behavior of swelling clay soils under varying confinement conditions.     

Application of Graded Finite Elements for Asphalt Pavements

Asphalt paving layers, particularly the surface course, exhibit vertically graded material properties. This grading is caused primarily by temperature gradients and aging related stiffness gradients. Most conventional existing analysis models do not directly account for the continuous grading of properties in flexible pavement layers. As a result, conventional analysis methods may lead to inaccurate prediction of pavement responses and distress under traffic and environmental loading. In this paper, a theoretical formulation for the graded finite element method is provided followed by an implementation using the user material subroutine (UMAT) capability of the finite element software ABAQUS. Numerical examples using the UMAT are provided to illustrate the benefits of using graded elements in pavement analysis.     

Modulus-Suction-Moisture Relationship for Compacted Soils in Postcompaction State

Despite clear evidence, changes in mechanical properties (i.e., stiffness or modulus) of compacted subgrades in response to subgrade moisture regime changes after construction have rarely been investigated in the geotechnical profession. In particular, when in-service assessment of pavement subgrade is made, the modulus-moisture variation should be addressed on the basis of unsaturated soil mechanics. This study presents the unsaturated small-strain modulus behavior of five predominately fine-grained compacted subgrade soils. The small-strain shear modulus (Go) of saturated compacted specimens subjected to a desorption soil-water characteristic curve (SWCC) was evaluated using bender elements. A test apparatus was designed to apply two stress state variables, the net confining pressure and matric suction, during the Go measurements. The relationship between Go and the SWCC under a constant mean net stress was developed. Additionally, the effect of compaction moisture content, compaction energy, and soil type on the Go-SWCC relationship was investigated. Finally, a relationship describing the small-strain modulus behavior of unsaturated compacted soils is proposed.     

Generalized Model for Geosynthetic-Reinforced Granular Fill-Soft Soil with Stone Columns

This paper pertains to the development of a mechanical model to predict the behavior of a geosynthetic-reinforced granular fill over soft soil improved with stone columns. The saturated soft soil has been idealized by Kelvin–Voight model to represent its consolidation behavior. The stone columns are idealized by stiffer springs. Pasternak shear layer and rough elastic membrane represent the granular fill and geosynthetic reinforcement layer, respectively. The nonlinear behavior of the granular fill and the soft soil is considered. Effect of consolidation of the soft soil due to inclusion of the stone columns has also been included in the model. Plane strain conditions are considered for the loading and reinforced foundation soil system. An iterative finite difference scheme is applied for obtaining the solution, and results are presented in nondimensional form. Comparison between the results from the present study and the analytical solution using theory of elasticity shows reasonable agreement. The advantage of using geosynthetic reinforcement is highlighted. Results indicate that inclusion of the geosynthetic layer effectively reduces the settlement. Nonlinearity in the behavior of the soft soil and the granular fill is reduced due to the use of geosynthetic reinforcement layer.     

Hysteresis of Capillary Stress in Unsaturated Granular Soil

Constitutive relationships among water content, matric suction, and capillary stress in unsaturated granular soils are modeled using a theoretical approach based on the changing geometry of interparticle pore water menisci. A series of equations is developed to describe the net force among particles attributable to the combined effects of negative pore water pressure and surface tension for spherical grains arranged in simple-cubic or tetrahedral packing order. The contact angle at the liquid–solid interface is considered as a variable to evaluate hysteretic behavior in the soil–water characteristic curve, the effective stress parameter χ, and capillary stress. Varying the contact angle from 0 to 40° to simulate drying and wetting processes, respectively, is shown to have an appreciable impact on hysteresis in the constitutive behavior of the modeled soils. A boundary between regimes of positive and negative pore water pressure is identified as a function of water content and contact angle. Results from the analysis are of practical importance in understanding the behavior of unsaturated soils undergoing natural wetting and drying processes, such as infiltration, drainage, and evaporation.     

Role of Matric Suction in Collapse of Compacted Clay Soil

This paper examines the influence of variations in matric suction on the collapse behavior of compacted Bangalore clay soil. The ASTM filter paper method measured the matric suctions of the compacted soil specimens. The matric suction of the compacted clay soil specimens ranged between 50 and 8,000 kPa at the as-compacted Sr values of 90 and 35%. Comparison of the matric suction-gravimetric water content relations of various compacted soils showed that the soil with a higher liquid limit has a higher matric suction at a given gravimetric water content. Variations in as-compacted degrees of saturation at a constant relative compaction or variations in relative compaction at a constant as-compacted degree of saturation notably affected the matric suction of the Bangalore clay soil. Experimental results also showed that the influence of matric suction on the collapse behavior of this compacted clay soil greatly depended on the relative compaction of the specimens and the pressure at which they were inundated.     

Model for Granular Materials with Surface Energy Forces

In light of environmental differences (such as gravitational fields, surface temperatures, atmospheric pressures, etc.), the mechanical behavior of the subsurface soil on the Moon is expected to be different from that on the Earth. Before any construction on the Moon can be envisaged, a proper understanding of soil properties and its mechanical behavior in these different environmental conditions is essential. This paper investigates the possible effect of surface-energy forces on the shear strength of lunar soil. All materials, with or without a net surface charge, exhibit surface-energy forces, which act at a very short range. Although, these forces are negligible for usual sand or silty sand on Earth, they may be important for surface activated particles under extremely low lunar atmospheric pressure. This paper describes a constitutive modeling method for granular material considering particle level interactions. Comparisons of numerical simulations and experimental results on Hostun sand show that the model can accurately reproduce the overall mechanical behavior of soils under terrestrial conditions. The model is then extended to include surface-energy forces between particles in order to describe the possible behavior of lunar soil under extremely low atmospheric pressure conditions. Under these conditions, the model shows that soil has an increase of shear strength due to the effect of surface-energy forces. The magnitude of increased shear strength is in reasonable agreement with the observations of lunar soil made on the Moon’s surface.     

Behavior of a Compacted Completely Decomposed Granite Soil from Suction Controlled Direct Shear Tests

A series of single-staged consolidated drained direct shear tests are carried out on recompacted completely decomposed granite (CDG) soil—a typical residual soil in Hong Kong, under different matric suctions and net normal stresses. Matric suction is controlled by applying air pressure in the pressure chamber and water pressure at the bottom of the high air-entry ceramic disk. The experimental results show that the contribution of suction to shear strength is significant. Shear strength of CDG soil increases with the increase of matric suction. Net normal stress has a remarkable influence on the shear strength of unsaturated CDG soil. The increase in shear strength due to an increase in matric suction (suction envelope) is observed as nonlinear i.e., ?b value varies with matric suction. No soil dilatancy is observed for zero matric suction (saturated case) but as the suction value is increased, higher soil dilatancy is obvious in lower net normal stresses. The rate of increase of soil dilatancy is greater in lower suction range than in higher suction range. The experimental shear strength data match closely with the shear strength predicted by existing shear strength model considering the soil-dilation effect.     

Unified DSC Constitutive Model for Pavement Materials with Numerical Implementation

Need for unified and mechanistic constitutive models for pavement materials for evaluation of various distresses has been recognized; however, such models are not yet available. There have been efforts to develop unified models; however, they have been based usually on ad hoc combinations of models for special properties such as elastic, plastic, creep and fracture, often without appropriate connections to various coupled responses of bound and unbound materials, they may result and in a large number of parameters, often without physical meanings. The disturbed state concept (DSC) provides a modeling approach that includes various responses such as elastic, plastic, creep, microcracking and fracture, softening and healing under mechanical and environmental (thermal, moisture, etc.) within a single unified and coupled framework. A brief review is presented to identify the advantages of the DSC compared to other available models. The DSC has been validated and applied to a wide range of materials: geologic, asphalt, concrete, ceramic, metal alloys, and silicon. It allows for evaluation of various distresses such as permanent deformations (rutting), microcracking and fracture, reflection cracking, thermal cracking, and healing. The DSC is implemented in two- and three-dimensional finite-element (FE) procedures, which allow static, repetitive, and dynamic loads including elastic, plastic, creep, microcracking leading to fracture and failure. A number of examples are solved for various distresses considering flexible (asphalt) pavements; however, the DSC model is applicable to rigid (concrete) pavements also. It is felt that the DSC and the FE computer programs provide unique and novel approaches for pavement engineering. It is desirable to perform further research and applications including validation with respect to simulated and field behavior of pavements.