A finite element solution scheme for an elastic–viscoplastic soil model

A Newton–Raphson solution scheme with a stress point algorithm is presented for the implementation of an elastic–viscoplastic soilmodel in a finite element program. Viscoplastic strain rates are calculated using the stress and volumetric states of the soil. Sub-incrementsof time are defined for each iterative calculation of elastic–viscoplastic stress changes so that their sum adds up to the time incrementfor the load step. This carefully defined ‘iterative time’ ensures that the correct amount of viscoplastic straining is accumulated overthe applied load step. The algorithms and assumptions required to implement the solution scheme are provided. Verification of the solutionscheme is achieved by using it to analyze typical boundary value problems.

Finite element analysis of an embankment on a soft estuarine deposit using an elastic–viscoplastic soil model

A new elastic–viscoplastic (EVP) soil model has been used to simulate the measured deformation response of a soft estuarine soil loaded by a stage-constructed embankment. The simulation incorporates prefabricated vertical drains installed in the foundation soils and reinforcement installed at the base of the embankment. The numerical simulations closely matched the temporal changes in surface settlement beneath the centerline and shoulder of the embankment. More importantly, the elastic–viscoplastic model simulated the pattern and magnitudes of the lateral deformations beneath the toe of the embankment — a notoriously difficult aspect of modelling the deformation response of soft soils. Simulation of the excess pore-water pressure proved more difficult because of the heterogeneous nature of the estuarine deposit. Excess pore-water pressures were, however, mapped reasonably well at three of the six monitoring locations. The simulations were achieved using a small set of material constants that can easily be obtained from standard laboratory tests. This study validates the use of the EVP model for problems involving soft soil deposits beneath loading from a geotechnical structure.

An improved elastic-viscoplastic soil model

This paper develops an improved and accessible framework for modelling time-dependent behaviour of soils using the concepts of elasticity and viscoplasticity. The mathematical description of viscoplastic straining is formulated based on a purely viscoplastic and measurable phenomenon, namely creep. The resulting expression for the viscoplastic strain rates includes a measure of both effective stress and the corresponding volumetric packing of the soil particles. In this way, the model differs from some earlier viscoplastic models and arguably provides a better conceptual description of time-dependent behaviour. Analytical solutions are developed for the simulation of drained and undrained strain-controlled triaxial compression tests. The model is then used to back-analyze the measured response of normally consolidated to moderately overconsolidated specimens of a soft estuarine soil in undrained triaxial compression. The model captures aspects of soil behaviour that cannot be simulated using time-independent elastic–plastic models. Specifically, it can capture the dependence of stress–strain relationships and undrained shear strength on strain rate, the development of irrecoverable plastic strains at constant stress (creep), and the relaxation of stresses at constant strain

An application of a physically based SOIL model to the Bowen ratio data

Yhteenveto: Haihdunnan mittaamisesta ja mallintamisesta peltoalueella

Numerical study of soil heterogeneity effects on contaminant transport in unsaturated soil (model development and validation)

The movement of chemicals through soil to groundwater is a major cause of degradation of water resources. In many cases, serious human and stock health implications are associated with this form of pollution. The study of the effects of different factors involved in transport phenomena can provide valuable information to find the best remediation approaches. Numerical models are increasingly being used for predicting or analyzing solute transport processes in soils and groundwater. This article presents the development of a stochastic finite element model for the simulation of contaminant transport through soils with the main focus being on the incorporation of the effects of soil heterogeneity in the model. The governing equations of contaminant transport are presented. The mathematical framework and the numerical implementation of the model are described. The comparison of the results obtained from the developed stochastic model with those obtained from a deterministic method and some experimental results shows that the stochastic model is capable of predicting the transport of solutes in unsaturated soil with higher accuracy than deterministic one. The importance of the consideration of the effects of soil heterogeneity on contaminant fate is highlighted through a sensitivity analysis regarding the variance of saturated hydraulic conductivity as an index of soil heterogeneity. © 2011 John Wiley & Sons, Ltd.

Wave propagation in nonlinear one-dimensional soil model

The objective of the research conducted by the authors is to explore the feasibility of determining reliable in situ values of shear modulus as a function of strain. In this paper the meaning of the material stiffness obtained from impact and harmonic excitation tests on a surface slab is discussed. A one-dimensional discrete model with the nonlinear material stiffness is used for this purpose. When a static load is applied followed by an impact excitation, if the amplitude of the impact is very small, the measured wave velocity using the cross-correlation indicates the wave velocity calculated from the tangent modulus corresponding to the state of stress caused by the applied static load. The duration of the impact affects the magnitude of the displacement and the particle velocity but has very little effect on the estimation of the wave velocity for the magnitudes considered herein. When a harmonic excitation is applied, the cross-correlation of the time histories at different depths estimates a wave velocity close to the one calculated from the secant modulus in the stress-strain loop under steady-state condition. Copyright © 2008 John Wiley & Sons, Ltd.

A simple model for predicting soil temperature in snow-covered and seasonally frozen soil: model description and testing

Microbial processes in soil are moisture, nutrient and temperature dependent and, consequently, accurate calculation of soil temperature is important for modelling nitrogen processes. Microbial activity in soil occurs even at sub-zero temperatures so that, in northern latitudes, a method to calculate soil temperature under snow cover and in frozen soils is required. This paper describes a new and simple model to calculate daily values for soil temperature at various depths in both frozen and unfrozen soils. The model requires four parameters average soil thermal conductivity, specific beat capacity of soil, specific heat capacity due to freezing and thawing and an empirical snow parameter. Precipitation, air temperature and snow depth (measured or calculated) are needed as input variables. The proposed model was applied to five sites in different parts of Finland representing different climates and soil types. Observed soil temperatures at depths of 20 and 50 cm (September 1981-August 1990) were used for model calibration. The calibrated model was then tested using observed soil temperatures from September 1990 to August 2001. R-2-values of the calibration period varied between 0.87 and 0.96 at a depth of 20 cm and between 0.78 and 0.97 at 50 cm. R-2 -values of the testing period were between 0.87 and 0.94 at a depth of 20cm. and between 0.80 and 0.98 at 50cm. Thus, despite the simplifications made, the model was able to simulate soil temperature at these study sites. This simple model simulates soil temperature well in the uppermost soil layers where most of the nitrogen processes occur. The small number of parameters required means, that the model is suitable for addition to catchment scale models.

Pull-out behaviour of steel grid soil reinforcement embedded in silty sand

This paper investigates the pull-out behaviour (particularly the bearing resistance) of a steel grid reinforcement embedded in silty sand using laboratory tests and numerical analyses. It is demonstrated that the various common analytical equations for calculating the bearing component of pull-out resistance give a wide range of calculated values, up to about 200% disparity. The disparity will increase further if the issue of whether to use the peak or critical state friction angle is brought in. Furthermore, these equations suggest that the bearing resistance factor, N, is only a function of soil friction angle which is not consistent with some design guidelines. In this investigation, a series of large scale laboratory pull-out tests under different test pressures were conducted. The test results unambiguously confirmed that the N factor is a function of test pressure. A modified equation for calculating N is also proposed. To have more in-depth understanding of the pull-out behaviour, the tests were modelled numerically. The input parameters for the numerical analysis were obtained from laboratory triaxial tests. The analysis results were compared with the experimental results. Good agreement between experimental and numerical results was achieved if the strain-softening behaviour from peak strength to critical state condition was captured by the soil model used. © 2013 Elsevier Ltd.

Assessing the environmental fate of selected polybrominated diphenyl ethers in the region surrounding the Zhuoshui River of Taiwan based on an Equilibrium Constant fugacity model

Polybrominated diphenyl ethers (PBDEs) are a group of flame retardants that have been in use since the 1970s. They are included in the list of hazardous substances known as persistent organic pollutants (POPs) because they are extremely hazardous to the environment and human health. PBDEs have been extensively used in industry and manufacturing in Taiwan, thus its citizens are at high risk of exposure to these chemicals.
An assessment of the environmental fate of these compounds in the Zhuoshui river and Changhua County regions of western Taiwan, and also including the adjacent area of the Taiwan Strait, was conducted for three high risk congeners, BDE-47, -99 and -209, to obtain information regarding the partitioning, advection, transfer and long range transport potential of the PBDEs in order to identify the level of risk posed by the pollutants in this region.
The results indicate that large amounts of PBDEs presently reside in all model compartments – air, soil, water, and sediment – with particularly high levels found in air and especially in sediment. The high levels found in sediment, particularly for BDE-209, are significant, since there is the threat of these pollutants entering the food chain, either directly through benthic feeding, or through resuspension and subsequent feeding in the pelagic region of the water column which is a distinct possibility in the strong currents found within the Taiwan Strait. Another important result is that a substantial portion of emissions leave the model domain directly through advection, particularly for BDE-47 (58%) and BDE-209 (75%), thus posing a risk to adjacent communities.
Model results were generally in reasonable agreement with available measured concentrations. In air, model concentrations are in reasonably good agreement with available measured values. For both BDE-47 and -99, model concentrations are a factor of 2-3 higher and BDE-209 within the range of measured values. In soil, model results are somewhat less than measured values. In sediment, model results are at the high end of measured values.

A REVIEW OF SOIL CONSTITUTIVE MODELS FOR SOIL STRUCTURE INTEREACTION ANALYSIS

Soils are multiphase materials comprised of mineral grains, air voids and water. Soils are not linearly elastic or perfectly plastic for external loading. Various constitutive models are available to describe the various aspects of soil behaviour. But no single soil model can completely describe the behaviour of real soil under all conditions. This paper attempts to compare various soil models and suggest a suitable model for the Soil Structure Interaction analysis especially for Kochi marine clay.

Community Structure in Soil Porous System

Diffusion controls the gaseous transport process in soils when advective transport is almost null. Knowledge of the soil structure and pore connectivity are critical issues to understand and modelling soil aeration, sequestration or emission of greenhouse gasses, volatilization of volatile organic chemicals among other phenomena. In the last decades these issues increased our attention as scientist have realize that soil is one of the most complex materials on the earth, within which many biological, physical and chemical processes that support life and affect climate change take place. A quantitative and explicit characterization of soil structure is difficult because of the complexity of the pore space. This is the main reason why most theoretical approaches to soil porosity are idealizations to simplify this system. In this work, we proposed a more realistic attempt to capture the complexity of the system developing a model that considers the size and location of pores in order to relate them into a network. In the model we interpret porous soils as heterogeneous networks where pores are represented by nodes, characterized by their size and spatial location, and the links representing flows between them. In this work we perform an analysis of the community structure of porous media of soils represented as networks. For different real soils samples, modelled as heterogeneous complex networks, spatial communities of pores have been detected depending on the values of the parameters of the porous soil model used. These types of models are named as Heterogeneous Preferential Attachment (HPA). Developing an exhaustive analysis of the model, analytical solutions are obtained for the degree densities and degree distribution of the pore networks generated by the model in the thermodynamic limit and shown that the networks exhibit similar properties to those observed in other complex networks. With the aim to study in more detail topological properties of these networks, the presence of soil pore community structures is studied. The detection of communities of pores, as groups densely connected with only sparser connections between groups, could contribute to understand the mechanisms of the diffusion phenomena in soils.

Community Structurein Soil Porous System.

Soil is well recognized as a highly complex system. The interaction and coupled physical, chemical, and biological processes and phenomena occurring in the soil environment at different spatial and temporal scales are the main reasons for such complexity. There is a need for appropriate methodologies to characterize soil porous systems with an interdisciplinary character. Four different real soil samples, presenting different textures, have been modeled as heterogeneous complex networks, applying a model known as the heterogeneous preferential attachment. An analytical study of the degree distributions in the soil model shows a multiscaling behavior in the connectivity degrees, leaving an empirically testable signature of heterogeneity in the topology of soil pore networks. We also show that the power-law scaling in the degree distribution is a robust trait of the soil model. Last, the detection of spatial pore communities, as densely connected groups with only sparser connections between them, has been studied for the first time in these soil networks. Our results show that the presence of these communities depends on the parameter values used to construct the network. These findings could contribute to understanding the mechanisms of the diffusion phenomena in soils, such as gas and water diffusion, development and dynamics of microorganisms, among others.

An agro-ecological simulation model system

In this paper five different models, as five modules of a complex agro-ecosystem are investigated. The water and nutrient flow in soil is simulated by the nutrient-in-soil model while the biomass change according to the seasonal weather aspects, the nutrient content of soil and the biotic interactions amongst the other terms of the food web are simulated by the food web population dynamical model that is constructed for a piece of homogeneous field. The food web model is based on the nutrient-in-soil model and on the activity function evaluator model that expresses the effect of temperature. The numbers of individuals in all phenological phases of the different populations are given by the phenology model. The food web model is extended to an inhomogeneous piece of field by the spatial extension model. Finally, as an additional module, an application of the above models for multivariate state-planes, is given. The modules built into the system are closely connected to each other as they utilize each other’s outputs, nevertheless, they work separately, too. Some case studies are analysed and a summarized outlook is given.