An examination of strain space versus stress space for the formulation of elastoplastic constitutive models

Several elastoplastic soil models have been proposed over the years that are formulated in strain space rather than stress space due to certain analytical and computational advantages. One such model, BRICK (Simpson 1992), has been continuously utilized and developed for industrial applications within Arup Geotechnics for more than two decades. This paper aims to describe the advantages and difficulties associated with strain space modeling. In addition, it will show how recent advances in modeling the effects of stress history, stiffness anisotropy, strength anisotropy and time-dependence in conventional stress space models can be transferred to the BRICK model. © 2010 Taylor & Francis Group, London.

The high strain rate response of Ultra High Molecular-weight Polyethylene: From fibre to laminate

The effect of strain rate upon the uniaxial response of Ultra High Molecular-weight Polyethylene (UHMWPE) fibres, yarns and laminates of lay-up [0/90]48 has been measured in both the 0/90 and ±45 configurations. The tensile strength of the matrix-dominated ±45 laminate is two orders of magnitude less than that of the fibre-dominated 0/90 laminate, and is more sensitive to strain rate. A piezoelectric force sensor device was developed to obtain the high strain rate data, and this achieved a rise time of less than 1 μs. It is found that the failure strength (and failure strain) of the yarn is almost insensitive to strain rate within the range (10 -1-103 s-1). At low strain rates (below 10 -1 s-1), creep of the yarn dominates and the failure strain increases with diminishing strain rate. The tensile strength of the dry yarn exceeds that of the laminate by about 20%. Tests on single fibres exceed the strength of the yarn by 20%. © 2013 Elsevier Ltd. All rights reserved.

Mechanical and hydraulic properties of Nankai accretionary prism sediments: Effect of stress path

We have conducted triaxial deformation experiments along different loading paths on prism sediments from the Nankai Trough. Different load paths of isotropic loading, uniaxial strain loading, triaxial compression (at constant confining pressure, Pc), undrained Pc reduction, drained Pc reduction, and triaxial unloading at constant Pc, were used to understand the evolution of mechanical and hydraulic properties under complicated stress states and loading histories in accretionary subduction zones. Five deformation experiments were conducted on three sediment core samples for the Nankai prism, specifically from older accreted sediments at the forearc basin, underthrust slope sediments beneath the megasplay fault, and overthrust Upper Shikoku Basin sediments along the frontal thrust. Yield envelopes for each sample were constructed based on the stress paths of Pc-reduction using the modified Cam-clay model, and in situ stress states of the prism were constrained using the results from the other load paths and accounting for horizontal stress. Results suggest that the sediments in the vicinity of the megasplay fault and frontal thrust are highly overconsolidated, and thus likely to deform brittle rather than ductile. The porosity of sediments decreases as the yield envelope expands, while the reduction in permeability mainly depends on the effective mean stress before yield, and the differential stress after yield. An improved understanding of sediment yield strength and hydromechanical properties along different load paths is necessary to treat accurately the coupling of deformation and fluid flow in accretionary subduction zones. © 2012 American Geophysical Union All Rights Reserved.

Modeling of strain rate effects on clay in simple shear

The response of clay is highly dependent on straining and loading rate. To obtain a realistic prediction of the response for time dependent problems, it is essential to use a model that accounts for rate effects in the stress-strain-strength properties of soils. The proposed model has been expanded from the existing SIMPLE DSS framework to account for the strain rate effects on clays in simple shear conditions. In accordance with the findings in the existing literature, soil response displays a unique relationship between shear strength and strain rate. The predicting model is illustrated with a limited test data. Copyright ASCE 2006.

Plane-strain discrete dislocation plasticity with climb-assisted glide motion of dislocations

A small-strain two-dimensional discrete dislocation plasticity (DDP) framework is developed wherein dislocation motion is caused by climb-assisted glide. The climb motion of the dislocations is assumed to be governed by a drag-type relation similar to the glide-only motion of dislocations: such a relation is valid when vacancy kinetics is either diffusion limited or sink limited. The DDP framework is employed to predict the effect of dislocation climb on the uniaxial tensile and pure bending response of single crystals. Under uniaxial tensile loading conditions, the ability of dislocations to bypass obstacles by climb results in a reduced dislocation density over a wide range of specimen sizes in the climb-assisted glide case compared to when dislocation motion is only by glide. A consequence is that, at least in a single slip situation, size effects due to dislocation starvation are reduced. By contrast, under pure bending loading conditions, the dislocation density is unaffected by dislocation climb as geometrically necessary dislocations (GNDs) dominate. However, climb enables the dislocations to arrange themselves into lower energy configurations which significantly reduces the predicted bending size effect as well as the amount of reverse plasticity observed during unloading. The results indicate that the intrinsic plasticity material length scale associated with GNDs is strongly affected by thermally activated processes and will be a function of temperature. © 2013 IOP Publishing Ltd.

Study on small-strain behaviours of methane hydrate sandy sediments using discrete element method

Methane hydrate bearing soil has attracted increasing interest as a potential energy resource where methane gas can be extracted from dissociating hydrate-bearing sediments. Seismic testing techniques have been applied extensively and in various ways, to detect the presence of hydrates, due to the fact that hydrates increase the stiffness of hydrate-bearing sediments. With the recognition of the limitations of laboratory and field tests, wave propagation modelling using Discrete Element Method (DEM) was conducted in this study in order to provide some particle-scale insights on the hydrate-bearing sandy sediment models with pore-filling and cementation hydrate distributions. The relationship between shear wave velocity and hydrate saturation was established by both DEM simulations and analytical solutions. Obvious differences were observed in the dependence of wave velocity on hydrate saturation for these two cases. From the shear wave velocity measurement and particle-scale analysis, it was found that the small-strain mechanical properties of hydrate-bearing sandy sediments are governed by both the hydrate distribution patterns and hydrate saturation. © 2013 AIP Publishing LLC.

Stiffness of clays and silts: Normalizing shear modulus and shear strain

An analysis is presented of a database of 67 tests on 21 clays and silts of undrained shear stress-strain data of fine-grained soils. Normalizations of secant G in terms of initial mean effective stress p9 (i.e., G=p9 versus log g) or undrained shear strength cu (i.e., G=cu versus log g) are shown to be much less successful in reducing the scatter between different clays than the approach that uses the maximum shear modulus,Gmax, a technique still not universally adopted by geotechnical researchers and constitutive modelers. Analysis of semiempirical expressions forGmax is presented and a simple expression that uses only a void-ratio function and a confining-stress function is proposed. This is shown to be superior to a Hardin-style equation, and the void ratio function is demonstrated as an alternative to an overconsolidation ratio (OCR) function. To derive correlations that offer reliable estimates of secant stiffness at any required magnitude of working strain, secant shear modulus G is normalized with respect to its small-strain value Gmax, and shear strain g is normalized with respect to a reference strain gref at which this stiffness has halved. The data are corrected to two standard strain rates to reduce the discrepancy between data obtained from static and cyclic testing. The reference strain gref is approximated as a function of the plasticity index.Aunique normalized shear modulus reduction curve in the shape of a modified hyperbola is fitted to all the available data up to shear strains of the order of 1%. As a result, good estimates can be made of the modulus reduction G/Gmax ±30% across all strain levels in approximately 90% of the cases studied. New design charts are proposed to update the commonly used design curves. © 2013 American Society of Civil Engineers.

Grain boundary interface mechanics in strain gradient crystal plasticity

Interactions between dislocations and grain boundaries play an important role in the plastic deformation of polycrystalline metals. Capturing accurately the behaviour of these internal interfaces is particularly important for applications where the relative grain boundary fraction is significant, such as ultra fine-grained metals, thin films and microdevices. Incorporating these micro-scale interactions (which are sensitive to a number of dislocation, interface and crystallographic parameters) within a macro-scale crystal plasticity model poses a challenge. The innovative features in the present paper include (i) the formulation of a thermodynamically consistent grain boundary interface model within a microstructurally motivated strain gradient crystal plasticity framework, (ii) the presence of intra-grain slip system coupling through a microstructurally derived internal stress, (iii) the incorporation of inter-grain slip system coupling via an interface energy accounting for both the magnitude and direction of contributions to the residual defect from all slip systems in the two neighbouring grains, and (iv) the numerical implementation of the grain boundary model to directly investigate the influence of the interface constitutive parameters on plastic deformation. The model problem of a bicrystal deforming in plane strain is analysed. The influence of dissipative and energetic interface hardening, grain misorientation, asymmetry in the grain orientations and the grain size are systematically investigated. In each case, the crystal response is compared with reference calculations with grain boundaries that are either 'microhard' (impenetrable to dislocations) or 'microfree' (an infinite dislocation sink). © 2013 Elsevier Ltd. All rights reserved.