A mathematical hysteretic model for elastomeric isolation bearings

An analytical model for high damping elastomeric isolation bearings is presented in this paper. The model is used to describe mathematically the damping force and restoring force of the rubber material and bearing. Ten parameters to be identified from cyclic loading tests are included in the model. The sensitivity of the ten parameters in affecting the model is examined. These ten parameters are functions of a number of influence factors on the elastomer such as the rubber compound, Mullins effect, scragging effect, frequency, temperature and axial load. In this study, however, only the Mullins effect, scragging effect, frequency and temperature are investigated. Both material tests and shaking table tests were performed to validate the proposed model. Based on the comparison between the experimental and the analytical results, it is found that the proposed analytical model is capable of predicting the shear force–displacement hysteresis very accurately for both rubber material and bearing under cyclic loading reversals. The seismic response time histories of the bearing can also be captured, using the proposed analytical model, with a practically acceptable precision. Copyright © 2002 John Wiley & Sons, Ltd.     

Modelling evaporation using an artificial neural network algorithm

This paper investigates the prediction of Class A pan evaporation using the artificial neural network (ANN) technique. The ANN back propagation algorithm has been evaluated for its applicability for predicting evaporation from minimum climatic data. Four combinations of input data were considered and the resulting values of evaporation were analysed and compared with those of existing models. The results from this study suggest that the neural computing technique could be employed successfully in modelling the evaporation process from the available climatic data set. However, an analysis of the residuals from the ANN models developed revealed that the models showed significant error in predictions during the validation, implying loss of generalization properties of ANN models unless trained carefully. The study indicated that evaporation values could be reasonably estimated using temperature data only through the ANN technique. This would be of much use in instances where data availability is limited. Copyright © 2002 John Wiley & Sons, Ltd.     

Analogue study of particle segregation in pyroclastic density currents, with implications for the emplacement mechanisms of large ignimbrites

Abstract Analogue flume experiments were conducted to investigate the transport and sedimentation behaviour of turbulent pyroclastic density currents. The experimental currents were scaled approximately to the natural environment in three ways: (1) they were fully turbulent; (2) they had a very wide range of particle sizes and associated Rouse numbers (the ratio of particle settling velocity to effective turbulent eddy velocity in the current); and (3) they contained particles of two different densities. Two sets of surge‐type experiments were conducted in a 5 m long, water‐filled lock‐exchange flume at five different volumetric particle concentrations from 0·6% to 23%. In one set (one‐component experiments), the currents contained just dense particles; in the other set (two‐component experiments), they contained both light and dense particles in equal volume proportions. In both sets of experiments, the population of each component had a log‐normal size distribution. In the two‐component experiments, the size range of the light particle population was selected in order to be in hydrodynamic equivalence with that of the dense particles. Dense particles were normally graded, both vertically and downstream, in the deposits from both sets of experiments. The mass loading (normalized to the initial mass of the suspension) and grain size of the dense component in the deposits decreased with distance from the reservoir and were insensitive to initial total particle concentration in the currents. On the other hand, in the two‐component experiments, the light particles were extremely sensitive to concentration. They were deposited in hydrodynamic equivalence with the dense particles from dilute currents, but were segregated efficiently at concentrations higher than a few per cent. With increasing particle concentration, the large, light particles were carried progressively further down the flume because of buoyancy effects. Deposits from the high‐concentration currents exhibited reverse vertical grading of the large, light particles. Efficient segregation of the light component was observed even if the bulk density of the current was less than that of the light particles. In both sets of experiments, marked inflexions in the rate of downstream decline in mass loading and maximum grain size of the dense component can be attributed to the presence of two different particle settling regimes in the flow: (1) particles with Rouse numbers >2·5, which did not respond to the turbulence and settled rapidly; and (2) particles with Rouse numbers <2·5, which followed the turbulent eddies and settled slowly. The results are applied to the transport and sedimentation dynamics of pyroclastic density currents that generate large, widespread ignimbrites. Field data fail to reveal significant departures from aerodynamic equivalence between pumice and lithic clasts in three such ignimbrites: the particulate loads of some large ignimbrites are transported principally in turbulent suspensions of low concentration. In some ignimbrites, the well‐developed inflexions in curves of maximum lithic (ML) size vs. distance can be attributed to the existence of distinct high and low Rouse number particle settling regimes that mark the transition from an overcharged state to one in which the residual particulate load is transported more effectively by turbulence.     

High-T, low-P metamorphism in the Palaeoproterozoic Halls Creek Orogen, northern Australia: the middle crustal response to a mantle-related transient thermal pulse

High‐T, low‐P metamorphic rocks of the Palaeoproterozoic central Halls Creek Orogen in northern Australia are characterised by low radiogenic heat production, high upper crustal thermal gradients (locally exceeding 40 °C km^?1) sustained for over 30 Myr, and a large number of layered mafic‐ultramafic intrusions with mantle‐related geochemical signatures. In order to account for this combination of geological and thermal characteristics, we model the middle crustal response to a transient mantle‐related heat pulse resulting from a temporary reduction in the thickness of the mantle lithosphere. This mechanism has the potential to raise mid‐crustal temperatures by 150–400 °C within 10–20 Myr following initiation of the mantle temperature anomaly, via conductive dissipation through the crust. The magnitude and timing of maximum temperatures attained depend strongly on the proximity, duration and lateral extent of the thermal anomaly in the mantle lithosphere, and decrease sharply in response to anomalies that are seated deeper than 50–60 km, maintained for <5 Myr in duration and/or have half‐widths <100 km. Maximum temperatures are also intimately linked to the thermal properties of the model crust, primarily due to their influence on the steady‐state (background) thermal gradient. The amplitudes of temperature increases in the crust are principally a function of depth, and are broadly independent of crustal thermal parameters. Mid‐crustal felsic and mafic plutonism is a predictable consequence of perturbed thermal regimes in the mantle and the lowermost crust, and the advection of voluminous magmas has the potential to raise temperatures in the middle crust very quickly. Although pluton‐related thermal signatures significantly dissipate within <10 Myr (even for very large, high‐temperature intrusive bodies), the interaction of pluton‐ and mantle‐related thermal effects has the potential to maintain host rock temperatures in excess of 400–450 °C for up to 30 Myr in some parts of the mid‐crust. The numerical models presented here support the notion that transient mantle‐related heat sources have the capacity to contribute significantly to the thermal budget of metamorphism in high‐T, low‐P metamorphic belts, especially in those characterised by low surface heat flow, very high peak metamorphic geothermal gradients and abundant mafic intrusions.     

The geological and geodynamic evolution of the eastern Black Sea basin: insights from 2-D and 3-D tectonic modelling

Subsidence mechanisms that may have controlled the evolution of the eastern Black Sea have been studied and simulated using a numerical model that integrates structural, thermal, isostatic and surface processes in both two- (2-D) and three-dimensions (3-D). The model enables the forward modelling of extensional basin evolution followed by deformation due to subsequent extensional and compressional events. Seismic data show that the eastern Black Sea has evolved via a sequence of interrelated tectonic events that began with early Tertiary rifting followed by several phases of compression, mainly confined to the edges of the basin. A large magnitude (approximately 12 km) of regional subsidence also occurred in the central basin throughout the Tertiary. Models that simulate the magnitude of observed fault controlled extension (β=1.13) do not reproduce the total depth of the basin. Similarly, the modelling of compressional deformation around the edges of the basin does little to enhance subsidence in the central basin. A modelling approach that quantifies lithosphere extension according to the amount of observed crustal thinning and thickening across the basin provides the closest match to overall subsidence. The modelling also shows that deep crustal and mantle–lithosphere processes can significantly influence the rate and magnitude of syn- to post-rift subsidence and shows that such mechanisms may have played an important role in forming the anomalously thin syn-rift and thick Miocene–Quaternary sequences observed in the basin. It is also suggested that extension of a 40–45 km thick pre-rift crust is required to generate the observed magnitude of total subsidence when considering a realistic bathymetry.     

Strike-slip fault bridge fluid pumping mechanism: insights from field-based palaeostress analysis and numerical modelling

We present a finite-element study of stress perturbation in evolving compressive and extensional strike-slip fault bridges. The results are compared with a fracture study of a compressive bridge at St Donats, UK. Horizontally interbedded calcareous mudstone and bioclastic calcilutite at St Donats have a distinct vertical permeability anisotropy. This sedimentary sequence behaves as a set of horizontal aquifers. The fluid flow in these aquifers is sensitive to mean stress gradients. Paleostress analysis of field fracture data, verified by finite-element modelling, indicates a rotation of σ₁ towards parallelism with boundary faults inside the growing compressive bridge. Boundary faults and bridge faults recorded numerous fluid flow events. The modelled mean stress pattern shows a regional maximum within the bridge and local maxima/minima pairs at boundary fault tips.Finite-element modelling of an extensional bridge indicates that σ₃ rotates towards parallelism with boundary faults. The mean stress pattern is similar to the pattern in compressive bridge but with maxima and minima locations interchanged. The stress patterns are reestablished by each stress build-up preceding the rupturation of the boundary faults throughout the development stages of strike-slip fault bridges. Mean stress gradients developed pre-failure control the fluid flow in fractures of the strike-slip fault system at and after the end of each stress build-up and the fluid flow in boundary faults post-failure. Fracture reactivation and new fracture generation within an evolving bridge is a process consisting of multiple successive events that retain the storage capacity of the bridge. Rupture and sealing of the main bounding-faults is a step-wise process that opens and closes fluid conduits between areas with different pressures.     

Cation distribution and structure modelling of spinel solid solutions

 This paper presents an improved generalisation of cation distribution determination based on an accurate fit of all crystal-chemical parameters. Cations are assigned to the tetrahedral and octahedral sites of the structure according to their scattering power and a set of bond distances optimised for spinel structure. A database of 295 spinels was prepared from the literature and unpublished data. Selected compositions include the following cations: Mg²⁺, Al³⁺, Si⁴⁺, Ti⁴⁺, V³⁺, Cr³⁺, Mn²⁺, Mn³⁺, Fe²⁺, Fe³⁺, Co²⁺, Ni²⁺, Zn²⁺ and vacancies. Bond distance optimisation reveals a definite lengthening in tetrahedral distance when large amounts of Fe³⁺ or Ni²⁺ are present in the octahedral site. This means that these cations modify the octahedral angle and hence the shared octahedral edge, causing an increase in the tetrahedral distance with respect to the size of the cations entering it. Some applications to published data are discussed, showing the capacity and limitations of the method for calculating cation distribution, and for identifying inconsistencies and inaccuracies in experimental data. Received: 19 February 2001 / Accepted: 1 June 2001     

Post-mid-Cretaceous eastern North Sea evolution inferred from 3D thermo-mechanical modelling

The Late Cretaceous–Cenozoic evolution of the eastern North Sea region is investigated by 3D thermo-mechanical modelling. The model quantifies the integrated effects on basin evolution of large-scale lithospheric processes, rheology, strength heterogeneities, tectonics, eustasy, sedimentation and erosion.

The evolution of the area is influenced by a number of factors: (1) thermal subsidence centred in the central North Sea providing accommodation space for thick sediment deposits; (2) 250-m eustatic fall from the Late Cretaceous to present, which causes exhumation of the North Sea Basin margins; (3) varying sediment supply; (4) isostatic adjustments following erosion and sedimentation; (5) Late Cretaceous–early Cenozoic Alpine compressional phases causing tectonic inversion of the Sorgenfrei–Tornquist Zone (STZ) and other weak zones.

The stress field and the lateral variations in lithospheric strength control lithospheric deformation under compression. The lithosphere is relatively weak in areas where Moho is deep and the upper mantle warm and weak. In these areas the lithosphere is thickened during compression producing surface uplift and erosion (e.g., at the Ringkøbing–Fyn High and in the southern part of Sweden). Observed late Cretaceous–early Cenozoic shallow water depths at the Ringkøbing–Fyn High as well as Cenozoic surface uplift in southern Sweden (the South Swedish Dome (SSD)) are explained by this mechanism.

The STZ is a prominent crustal structural weakness zone. Under compression, this zone is inverted and its surface uplifted and eroded. Contemporaneously, marginal depositional troughs develop. Post-compressional relaxation causes a regional uplift of this zone.

The model predicts sediment distributions and paleo-water depths in accordance with observations. Sediment truncation and exhumation at the North Sea Basin margins are explained by fall in global sea level, isostatic adjustments to exhumation, and uplift of the inverted STZ. This underlines the importance of the mechanisms dealt with in this paper for the evolution of intra-cratonic sedimentary basins.     

Lithosphere structure of European sedimentary basins from regional three-dimensional gravity modelling

A three-dimensional (3D) density model, approximated by two regional layers—the sedimentary cover and the crystalline crust (offshore, a sea-water layer was added), has been constructed in 1° averaging for the whole European continent. The crustal model is based on simplified velocity model represented by structure maps for main seismic horizons—the “seismic” basement and the Moho boundary. Laterally varying average density is assumed inside the model layers. Residual gravity anomalies, obtained by subtraction of the crustal gravity effect from the observed field, characterize the density heterogeneities in the upper mantle. Mantle anomalies are shown to correlate with the upper mantle velocity inhomogeneities revealed from seismic tomography data and geothermal data. Considering the type of mantle anomaly, specific features of the evolution and type of isostatic compensation, the sedimentary basins in Europe may be related into some groups: deep sedimentary basins located in the East European Platform and its northern and eastern margins (Peri-Caspian, Dnieper–Donets, Barents Sea Basins, Fore–Ural Trough) with no significant mantle anomalies; basins located on the activated thin crust of Variscan Western Europe and Mediterranean area with negative mantle anomalies of −150 to −200×10⁻⁵ ms⁻² amplitude and the basins associated with suture zones at the western and southern margins of the East European Platform (Polish Trough, South Caspian Basin) characterized by positive mantle anomalies of 50–150×10⁻⁵ ms⁻² magnitude. An analysis of the main features of the lithosphere structure of the basins in Europe and type of the compensation has been carried out.