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.     

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.     

A sapphirine–phlogopite–cordierite paragenesis in a low-P amphibolite facies terrain, Arunta Inlier, Australia

Summary Silica-undersaturated phlogopite schists from the Cackleberry Metamorphics, Arunta Inlier, central Australia, preserve relatively low-temperature sapphirine-bearing parageneses that developed during low-pressure upper amphibolite facies metamorphism. Peak metamorphic phlogopite–cordierite–sapphirine assemblages are interpreted to have formed during the same event recorded in nearby metapelites, at c.3 kbar and 650–700 °C. Initial cooling of the terrain resulted in the breakdown of sapphirine to corundum–chlorite–phlogopite and corundum–spinel–chlorite assemblages. Further retrogression at greenschist facies conditions resulted in the replacement of sapphirine by diaspore–chlorite intergrowths. The reaction textures are consistent with a near-isobaric heating-cooling path at low-pressure, and provide evidence for the stability of sapphirine at c.700 °C at low pressures in rocks of an appropriate Mg- and Fe³⁺-rich bulk composition. Received August 15, 2001 accepted December 27, 2001     

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.     

Solid bitumens: an assessment of their characteristics, genesis, and role in geological processes

Familiar since antiquity, and subject in contemporary times to various characterization schemes, the exact nature of solid bitumen is not yet fully known. Bitumens have ‘random polymer-like’ molecular structures, are mobile as highly viscous fluids or were once fluids but have since turned into solids. Solid bitumens consist mainly of large moieties, of polyclyclic aromatic hydrocarbons, occasionally with finely admixed, fine-grained cryptocrystalline graphite. Solid bitumens are distinguished from kerogen, which is the syngenetic and generally finely dispersed particulate organic matter in sedimentary rock that virtually does not migrate following its deposition. Occurrences of solid bitumens are relevant to petroleum exploration as well as the search for, and evaluation of, a variety of metallic mineral deposits. Genesis of bitumen is in many cases linked to the thermal and hydrothermal history of organic matter in sedimentary rock. Apparently bitumen, or more specifically organic acids generated along with bitumen during diagenesis, may alter porosity of reservoir rocks or otherwise prepare the ground for ore deposition. Bitumen is also relatively sensitive to alteration processes, some of which, such as oxidative weathering, water leaching, biodegradation (contact) metamorphism and ionizing radiation may likewise affect its nature. Elemental composition of bitumen commonly reflects the nature of mineral deposits. Is is possible that in petroleum exploration, trace metal abundances of bitumen may eventually allow prediction of crude oil types and volumes anticipated from a given source rock? Beside transition elements, notably Ni and V, highly anomalous concentrations of U, Pt and Au occur in some solid bitumens. During the generation of petroleum from kerogen, the trend in δ¹³C is toward lighter values. The opposite seems to occur when liquid petroleum is subjected to thermal cracking (and /or related processes) yielding solid bitumen enriched in ¹³C, and isotopically light methane. In fact, except for deasphalting and possibly some irradiation processes, the result of thermal cracking, oxidation, water leaching, inspissation (drying) and bacterial degradation of crude oil is that lower molecular weight hydrocarbons are removed leaving bitumen residues enriched in aromatic hydrocarbons, heteroatomic compounds (NSO) and ¹³C. Such phenomena are relevant to bitumen paragenesis in petroleum reservoir rocks, to certain Phanerozoic occurrences of multiple generations of bitumens, and to bitumens in mineral deposits.     

Geophysical investigations over a segment of the Central Indian Ridge, Indian Ocean

 Swath bathymetric, gravity, and magnetic studies were carried out over a 55 km long segment of the Central Indian Ridge. The ridge is characterized by 12 to 15 km wide rift valley bounded by steep walls and prominent volcanic constructional ridges on either side of the central rift valley. A transform fault at 7°45′S displaces the ridge axis. A mantle Bouguer anomaly low of −14 mGals and shallowing of rift valley over the middle of the ridge segment indicate along axis crustal thickness variations. A poorly developed neovolcanic zone on the inner rift valley floor indicate dominance of tectonic extension. The off-axis volcanic ridgs suggest enhanced magmatic activity during the recent past. Received: 24 May 1996 / Rivision received: 13 January 1997