Generalized strain probing of constitutive models

Advanced material constitutive models are used to describe complex soil behaviour. These models are often used in the solution of boundary value problems under general loading conditions. Users and developers of constitutive models need to methodically investigate the represented soil response under a wide range of loading conditions. This paper presents a systematic procedure for probing constitutive models. A general incremental strain probe, 6D hyperspherical strain probe (HSP), is introduced to examine rate‐independent model response under all possible strain loading conditions. Two special cases of HSP, the true triaxial strain probe (TTSP) and the plane‐strain strain probe (PSSP), are used to generate 3‐D objects that represent model stress response to probing. The TTSP, PSSP and general HSP procedures are demonstrated using elasto‐plastic models. The objects resulting from the probing procedure readily highlight important model characteristics including anisotropy, yielding, hardening, softening and failure. The PSSP procedure is applied to a Neural Network (NN) based constitutive model. It shows that this probing is especially useful in understanding NN constitutive models, which do not contain explicit functions for yield surface, hardening, or anisotropy. Copyright © 2004 John Wiley & Sons, Ltd.     

A one-dimensional model for the interaction between continental-scale ice sheets and atmospheric stationary waves

 The great continental ice sheets of the Pleistocene represented a significant topographic obstacle to the westerly winds in northern midlatitudes. This work explores how consequent changes in the atmospheric stationary wave pattern might have affected the shape and growth of the ice sheets themselves. A one dimensional (1-D) model is developed which permits an examination of the types and magnitudes of the feedbacks that might be expected. When plausible temperature perturbations are introduced at the ice-sheet margin which are proportional to the stationary wave amplitude, the equilibrium shape of the ice sheet is significantly altered, and depends on the sign of the perturbation. The proposed feedback also affects the response of the ice sheet to time-varying climate forcing. The results suggest that the evolution of a continental-scale ice sheet with a land-based margin may be significantly determined by the changes it induces in the atmospheric circulation. Received: 1 October 1999 / Accepted: 17 July 2000     

Theoretical accuracy of synthetic aperture sonar micronavigation using a displaced phase-center antenna

The Cramer-Rao lower bounds on the cross-track translation and rotation of a displaced phase-center antenna (DPCA) in the slant range plane between two successive pings (known as DPCA sway and yaw in what follows) are computed, assuming statistically homogeneous backscatter. These bounds are validated using experimental data from a 118-182-kHz sonar, showing an accuracy of the order of 20 microns on the ping-to-ping cross-track displacements. Next, the accuracy required on the DPCA sway and yaw in order to achieve a given synthetic aperture sonar (SAS) beampattern specification, specified by the expected SAS array gain, is computed as a function of the number P of pings in the SAS. Higher accuracy is required when P increases to counter the accumulation of errors during the integration of the elementary ping-to-ping estimates: the standard deviation must decrease as P/sup -1/2/ for the DPCA sway and P/sup -3/2/ for the yaw. Finally, by combining the above results, the lower bounds on DPCA micronavigation accuracy are established. These bounds set an upper limit to the SAS length achievable in practice. The maximum gain Q in cross-range resolution achievable by a DPCA micronavigated SAS is computed as a function of the key SAS parameters. These theoretical predictions are compared with simulations and experimental results.     

Dissolution and depletion of ferromagnesian minerals from Holocene tephra layers in an acid bog,New Zealand,and implications for tephra correlation

This study examines the depletion of ferromagnesian silicate minerals from a sequence of thin, distal, mainly rhyolitic tephra layers of Holocene age preserved in an acid peat bog (Kopouatai), North Island, New Zealand. The rate of such depletion has been fast, as indicated by the complete loss of biotite from one tephra layer (Kaharoa Tephra), in which it is normally dominant, in only ca. 770 yr. Chemical dissolution is advocated as the likely cause for the depletion, with amphiboles and other mineral grains commonly showing etch pits, microcaves, and other characteristic surface solution features. Theoretical thermodynamic and kinetic models show a marked increase in the rate of dissolution of all ferromagnesian minerals under conditions of low pH (< 4), but that where silica concentrations in solution are high the relative proportions of minerals remaining are unaffected. However, where concentrations of dissolved silica are low, as in most bog environments, the relative proportions of ferromagnesian minerals are affected as well as absolute amounts being decreased. Amphiboles are depleted relative to pyroxenes, consistent with kinetic studies. The results show that the identification and correlation of tephras on the basis of relative abundances of ferromagnesian minerals alone may be unreliable, and emphasise the need to use multiple criteria in such studies.     

Response of circular footings and anchor plates in non-homogeneous elastic soils

The axisymmetric elastic response of circular footings and anchor plates in a linearly non-homogeneous elastic soil is analysed. It is assumed that footings/anchors are flexible and subjected to axisymmetric vertical loads. The response of the footing/anchor is modelled by using the classical Poisson–Kirchhoff thin plate theory. A variational technique is used to analyse the interaction problem. A representation for the contact stress is established by using a fundamental solution corresponding to a unit vertical pressure acting over an annular region in the interior of the non-homogeneous soil. The fundamental solution can be derived by using rigorous analytical procedures. The influence of the footing flexibility and the degree of soil non-homogeneity on the displacements, bending moments and contact stresses of a surface footing is examined over a wide range of governing parameters. In the case of anchor plates the influence of depth of embedment, degree of soil non-homogeneity and anchor flexibility on the anchor displacement is investigated.     

Sources and sinks of sea salt in a newfoundland blanket bog

Over an oceanic peatland, the concentration of Na in fog averaged 38.1 mgl^?1 compared with 1.8 mgl^?1 in rain, resulting in a significant flux of mineral elements to the surface. Between 16 May and 20 June 1990 the average mass flux of Na to the bog surface by fog, rain, and dry deposition was 21.9, 10.4 and 7.0 mg m^?2 d^?1. There was little long-term storage of Na within the peatland system, where Na losses measured in stream runoff averaged 34.8 mg m² d^?1, and deep groundwater losses 4 mg m^?2 d^?1. Calcium and Mg were preferentially retained in the organic soil, whereas K was relatively mobile. Potassium tended to become concentrated in the unsaturated zone. Stream runoff had a consistently higher pH than groundwater, corresponding to higher Ca and Mg concentrations, which may have been from mineral sources in the headwater ponds. Otherwise, the stream water chemistry was closely related to groundwater in the upper layers of the peat deposit.     

Comparison of results from several AOGCMs for global and regional sea-level change 1900–2100

 Sea-level rise is an important aspect of climate change because of its impact on society and ecosystems. Here we present an intercomparison of results from ten coupled atmosphere-ocean general circulation models (AOGCMs) for sea-level changes simulated for the twentieth century and projected to occur during the twenty first century in experiments following scenario IS92a for greenhouse gases and sulphate aerosols. The model results suggest that the rate of sea-level rise due to thermal expansion of sea water has increased during the twentieth century, but the small set of tide gauges with long records might not be adequate to detect this acceleration. The rate of sea-level rise due to thermal expansion continues to increase throughout the twenty first century, and the projected total is consequently larger than in the twentieth century; for 1990–2090 it amounts to 0.20–0.37 m. This wide range results from systematic uncertainty in modelling of climate change and of heat uptake by the ocean. The AOGCMs agree that sea-level rise is expected to be geographically non-uniform, with some regions experiencing as much as twice the global average, and others practically zero, but they do not agree about the geographical pattern. The lack of agreement indicates that we cannot currently have confidence in projections of local sea-level changes, and reveals a need for detailed analysis and intercomparison in order to understand and reduce the disagreements. Received: 1 September 2000 / Accepted: 20 April 2001     

Grain boundary diffusion in enstatite

We have investigated grain boundary diffusion rates in enstatite by heating single crystals of quartz packed in powdered San Carlos olivine (Mg_0.90Fe_0.10)₂SiO₄ at controlled oxygen fugacities in the range 10^?5.7 to 10^?8.7?atm and temperatures from 1350° to 1450?°C for times from 5 to 100?h at 1?atm total pressure. Following the experiments, the thickness of the coherent polycrystalline reaction rim of pyroxene that had formed between the quartz and olivine was measured using backscatter scanning imaging in the electron microprobe. Quantitative microprobe analysis indicated that the composition of this reaction phase is (Mg_0.92Fe_0.08)₂Si₂O₆. The rate of growth of the pyroxene increases with increasing temperature, is independent of the oxygen fugacity, and is consistent with a parabolic rate law, indicating that the growth rate is controlled by ionic diffusion through the pyroxene rim. Microstructural observations and platinum marker experiments suggest that the reaction phase is formed at the olivine-pyroxene interface, and is therefore controlled by the diffusion of silicon and oxygen. The parabolic rate constants determined from the experiments were analyzed in terms of the oxide activity gradient across the rim to yield mean effective diffusivities for the rate-limiting ionic species, assuming bulk transport through the pyroxene layer. These effective diffusivities are faster than the lattice diffusivities for the slowest species (silicon) calculated from creep experiments, but slower than measured lattice diffusivities for oxygen in enstatite. Thus, silicon grain boundary diffusion is most likely to be the rate-limiting process in the growth of the pyroxene rims. Also, as oxygen transport through the pyroxene rims must be faster than silicon transport, diffusion of oxygen along the grain boundaries must be faster than through the lattice. The grain boundary diffusivity for silicon in orthopyroxenite is then given by D¯gbSiδ=(3.3±3.0)×10^?9f0.0O₂e^?400±65/RT?m³s^?1, where the activation energy for diffusion is in kJ/mol, and δ is the grain boundary width in m. Calculated growth rates for enstatite under these conditions are significantly slower than predicted by an extrapolation from similar experiments performed at 1000?°C under high pressure (hydrous) conditions by Yund and Tullis (1992), perhaps due to water-enhancement of diffusion in their experiments.