The influence of particle shape on the stress-strain and creep response of a fine silica sand

The influence of particle shape on the stress-strain response of fine silica sand is investigated experimentally. Two sands from the same source and with the same particle size distribution were examined using Fourier descriptor analysis for particle shape. Their grains were, on average, found to have similar angularity but different elongation. During triaxial stress path testing, the stress-strain behavior of the sands for both loading and creep stages were found to be influenced by particle elongation. In particular, the behavior of the sand with less elongated grains was more like that of rounded glass beads during creep. The results highlight the role of particle shape in stress transmission in granular packings and suggest that shape should be taken more rigorously into consideration in characterizing geomaterials. © 2005 Taylor & Francis Group.

Stress-strain response of hydrate-bearing sands: Numerical study using discrete element method simulations

Gas hydrate is a crystalline solid found within marine and subpermafrost sediments. While the presence of hydrates can have a profound effect on sediment properties, the stress-strain behavior of hydrate-bearing sediments is poorly understood due to inherent limitations in laboratory testing. In this study, we use numerical simulations to improve our understanding of the mechanical behavior of hydrate-bearing sands. The hydrate mass is simulated as either small randomly distributed bonded grains or as "ripened hydrate" forming patchy saturation, whereby sediment clusters with 100% pore-filled hydrate saturation are distributed within a hydrate-free sediment. Simulation results reveal that reduced sand porosity and higher hydrate saturation cause an increase in stiffness, strength, and dilative tendency, and the critical state line shifts toward higher void ratio and higher shear strength. In particular, the critical state friction angle increases in sands with patchy saturation, while the apparent cohesion is affected the most when the hydrate mass is distributed in pores. Sediments with patchy hydrate distribution exhibit a slightly lower strength than sediments with randomly distributed hydrate. Finally, hydrate dissociation under drained conditions leads to volume contraction and/or stress relaxation, and pronounced shear strains can develop if the hydrate-bearing sand is subjected to deviatoric loading during dissociation.

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.

Study of the turbulence in the air-side and water-side boundary layers in experimental laboratory wind induced surface waves

This study detailed the structure of turbulence in the air-side and water-side boundary layers in wind-induced surface waves. Inside the air boundary layer, the kurtosis is always greater than 3 (the value for normal distribution) for both horizontal and vertical velocity fluctuations. The skewness for the horizontal velocity is negative, but the skewness for the vertical velocity is always positive. On the water side, the kurtosis is always greater than 3, and the skewness is slightly negative for the horizontal velocity and slightly positive for the vertical velocity. The statistics of the angle between the instantaneous vertical fluctuation and the instantaneous horizontal velocity in the air is similar to those obtained over solid walls. Measurements in water show a large variance, and the peak is biased towards negative angles. In the quadrant analysis, the contribution of quadrants Q2 and Q4 is dominant on both the air side and the water side. The non-dimensional relative contributions and the concentration match fairly well near the interface. Sweeps in the air side (belonging to quadrant Q4) act directly on the interface and exert pressure fluctuations, which, in addition to the tangential stress and form drag, lead to the growth of the waves. The water drops detached from the crest and accelerated by the wind can play a major role in transferring momentum and in enhancing the turbulence level in the water side.On the air side, the Reynolds stress tensor's principal axes are not collinear with the strain rate tensor, and show an angle α σ≈=-20°to-25°. On the water side, the angle is α σ≈=-40°to-45°. The ratio between the maximum and the minimum principal stresses is σ a/σ b=3to4 on the air side, and σ a/σ b=1.5to3 on the water side. In this respect, the air-side flow behaves like a classical boundary layer on a solid wall, while the water-side flow resembles a wake. The frequency of bursting on the water side increases significantly along the flow, which can be attributed to micro-breaking effects - expected to be more frequent at larger fetches. © 2012 Elsevier B.V.

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.

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.