A micromechanical finite element model for linear and damage‐coupled viscoelastic behaviour of asphalt mixture

This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd.     

Wavelet analysis of inter‐annual variability in the runoff regimes of glacial and nival stream catchments,Bow Lake,Alberta

Continuous wavelet analyses of hourly time series of air temperature, stream discharge, and precipitation are used to compare the seasonal and inter‐annual variability in hydrological regimes of the two principal streams feeding Bow Lake, Banff National Park, Alberta: the glacial stream draining the Wapta Icefields, and the snowmelt‐fed Bow River. The goal is to understand how water sources and flow routing differ between the two catchments. Wavelet spectra and cross‐wavelet spectra were determined for air temperature and discharge from the two streams for summers (June–September) 1997–2000, and for rainfall and discharge for the summers of 1999 and 2000. The diurnal signal of the glacial runoff was orders of magnitude higher in 1998 than in other years, indicating that significant ice exposure and the development of channelized glacial drainage occurred as a result of the 1997–98 El Niño conditions. Early retreat of the snowpack in 1997 and 1998 led to a significant summer‐long input of melt runoff from a small area of ice cover in the Bow River catchment; but such inputs were not apparent in 1999 and 2000, when snow cover was more extensive. Rainfall had a stronger influence on runoff and followed quicker flow paths in the Bow River catchment than in the glacial catchment. Snowpack thickness and catchment size were the primary controls on the phase relationship between temperature and discharge at diurnal time scales. Wavelet analysis is a fast and effective means to characterize runoff, temperature, and precipitation regimes and their interrelationships and inter‐annual variability. The technique is effective at identifying inter‐annual and seasonal changes in the relative contributions of different water sources to runoff, and changes in the time required for routing of diurnal meltwater pulses through a catchment. However, it is less effective at identifying changes/differences in the type of the flow routing (e.g. overland flow versus through flow) between or within catchments. Copyright © 2003 John Wiley & Sons, Ltd.     

Delamerian‐aged metamorphism in the southern Curnamona Province,Australia: implications for the evolution of the Mesoproterozoic Olarian Orogeny

Monazite electron microprobe U–Th–Pb and garnet Sm–Nd isotopic data from metapelitic assemblages in the Willyama Supergroup in the southern Curnamona Province, south‐central Australia, indicate that the terrain underwent regional greenschist to amphibolite‐grade metamorphism during the c. 500 Ma Delamerian Orogeny. The Delamerian‐aged mineral assemblages include prograde garnet–staurolite and kyanite‐bearing associations that overprint andalusite‐ and sillimanite‐bearing assemblages that are interpreted to have developed during the c. 1600 Ma Olarian Orogeny. Importantly, the development of secondary kyanite‐bearing assemblages in the southern Curnamona Province has been used previously to suggest that the Olarian Orogeny followed an anticlockwise P–T evolution. If such assemblages are the product of c. 500 Ma metamorphism, then the anticlockwise P–T path is an apparent path, due to the overprint of a distinct metamorphic cycle c. 1100 Ma later. Making such distinctions is therefore extremely important when using the textural and metamorphic evolution of polycyclic terrains to model the thermal behaviour of the crust during orogeny. This study highlights the utility of in situ geochronology, linking age data to petrologically important phases and assemblages.     

SANICLAY: simple anisotropic clay plasticity model

SANICLAY is a new simple anisotropic clay plasticity model that builds on a modification of an earlier model with an associated flow rule, in order to include simulations of softening response under undrained compression following K_o consolidation. Non‐associativity is introduced by adopting a yield surface different than the plastic potential surface. Besides, the isotropic hardening of the yield surface both surfaces evolve according to a combined distortional and rotational hardening rule, simulating the evolving anisotropy. Although built on the general premises of critical state soil mechanics, the model induces a critical state line in the void ratio–mean effective stress space, which is a function of anisotropy. To ease interpretation, the model formulation is presented firstly in the triaxial stress space and subsequently, its multiaxial generalization is developed systematically, in a form appropriate for implementation in numerical codes. The SANICLAY is shown to provide successful simulation of both undrained and drained rate‐independent behaviour of normally consolidated sensitive clays, and to a satisfactory degree of accuracy of overconsolidated clays. The new model requires merely three constants more than those of the modified Cam clay model, all of which are easily calibrated from well‐established laboratory tests following a meticulously presented procedure. Copyright © 2006 John Wiley & Sons, Ltd.     

WEKU – a GIS-Supported stochastic model of groundwater residence times in upper aquifers for the supraregional groundwater management

 The supraregional GIS-supported stochastical model, WEKU, for the determination of groundwater residence times in the upper aquifers of large groundwater provinces is presented. Using a two-dimensional analytical model of groundwater flow, groundwater residence times are determined within two extreme cases. In the first case, maximal groundwater residence times are calculated, representing the part of groundwater, that is drained by the main surface water of a groundwater catchment area. In the second case, minimal groundwater residence times for drainage into the nearest surface water are determined. Using explicit distribution functions of the input parameters, mean values as well as potential ranges of variations of the groundwater residence times are derived. The WEKU model has been used for the determination of groundwater residence times throughout Germany. The model results – mean values and deviations of the groundwater velocity and the maximal and minimal groundwater residence times in the upper aquifers – are presented by general maps and discussed in detail. It is shown that the groundwater residence times in the upper aquifer vary regionally, differentiated between less than 1 year and more than 2000 years. Using this information, the time scales can be specified, until measures to remediate polluted groundwater resources may lead to a substantial groundwater quality improvement in the different groundwater provinces of Germany. With respect to its supraregional scale of application, the WEKU model may serve as a useful tool for the supraregional groundwater management on a state, federal or international level. Received: 15 August 1995 · Accepted: 15 October 1995     

Radon variations in a Hungarian village

 A steady radon exhalation is assumed in most publications. In a village of North-East Hungary, however, high radon concentrations have been measured, differing strongly in neighbouring houses and varying in time, due to the interplay of geochemical phenomena. Received: 20 November 1995 · Accepted: 18 June 1996     

Reliability assessment of serviceability performance of braced retaining walls using a neural network approach

In urban environments, one major concern with deep excavations in soft clay is the potentially large ground deformations in and around the excavation. Excessive movements can damage adjacent buildings and utilities. There are many uncertainties associated with the calculation of the ultimate or serviceability performance of a braced excavation system. These include the variabilities of the loadings, geotechnical soil properties, and engineering and geometrical properties of the wall. A risk‐based approach to serviceability performance failure is necessary to incorporate systematically the uncertainties associated with the various design parameters. This paper demonstrates the use of an integrated neural network–reliability method to assess the risk of serviceability failure through the calculation of the reliability index. By first performing a series of parametric studies using the finite element method and then approximating the non‐linear limit state surface (the boundary separating the safe and ‘failure’ domains) through a neural network model, the reliability index can be determined with the aid of a spreadsheet. Two illustrative examples are presented to show how the serviceability performance for braced excavation problems can be assessed using the reliability index. Copyright © 2005 John Wiley & Sons, Ltd.