Transient analysis of wave propagation in non‐homogeneous elastic unbounded domains by using the scaled boundary finite‐element method

The scaled boundary finite‐element method has been developed for the dynamic analysis of unbounded domains. In this method only the boundary is discretized resulting in a reduction of the spatial dimension by one. Like the finite‐element method no fundamental solution is required. This paper extends the scaled boundary finite‐element method to simulate the transient response of non‐homogeneous unbounded domains with the elasticity modulus and mass density varying as power functions of spatial coordinates. To reduce the number of degrees of freedom and the computational cost, the technique of reduced set of base functions is applied. The scaled boundary finite‐element equation for an unbounded domain is reformulated in generalized coordinates. The resulting acceleration unit‐impulse response matrix is obtained and assembled with the equation of motion of standard finite elements. Numerical examples of non‐homogeneous isotropic and transversely isotropic unbounded domains demonstrate the accuracy of the scaled boundary finite‐element method. Copyright © 2006 John Wiley & Sons, Ltd.     

Ground viscosity and stiffness measurements for near surface seismic velocity

In land surveys, the weathering layer can often distort the seismic signal due to it passing through rapid velocity and density changes, dispersion, scattering and inelastic absorption. In a simple spring‐dashpot model for the earth response, an equivalent medium groups these complex phenomena into two parameters only; these are called ground viscosity and ground stiffness. The most recent controllers for vibrators can estimate both parameters. To validate these measurements, Saudi Aramco conducted an experiment measuring ground viscosity and stiffness from two different vibrator control systems over an area of varying terrain conditions, including unconsolidated sand and limestone outcrop. The two systems measured different values, but detected similar trends that correlated well with weathering conditions and surface geology, e.g. lower viscosity values on the outcrop than on the sand. The ratio of ground viscosity to ground stiffness can approximate the shallow S‐wave velocity, which we converted into P‐wave velocity through calibration with sparse uphole data. Static corrections incorporating this velocity information somewhat improved the focusing of seismic time sections. This new approach does not require additional acquisition efforts, and can model shallow complex formations in arid areas where classical static methods often fail.     

Early Cretaceous sedimentary provenance and structural evolution of the central Sanandaj–Sirjan Zone,Iran: implications for palaeogeographic reconstructions of the northern Neo-Tethyan margin

ABSTRACT

The Early Cretaceous was an important epoch in the evolution of the Earth system in which major tectonic episodes occurred, especially along the Alpine–Himalayan belt. The paucity of reliable palaeogeographic data from the central segment of this geological puzzle, however, hampers the reconstruction of a panoramic view of its Early Cretaceous palaeogeography and geodynamic setting. Here we present multidisciplinary provenance data from Lower Cretaceous strata of the overriding plate of the Neo-Tethyan subduction zone (the Sanandaj–Sirjan Zone; SSZ, of central Iran), including structural, basin-fill evolution, petrographic and geochemical analyses. Sandstone provenance analysis of Lower Cretaceous red beds suggests the occurrence of sub-mature litho-quartzose sandstones attributed to an active continental arc margin in convergent setting predominantly derived from plutonic, quartzose sedimentary and metamorphic rocks exposed in the central SSZ. Weathering indices indicate moderate chemical weathering in the source area which may be related to close source-to-sink relationships or arid climate. Our palaeogeographic reconstructions and original geological mapping indicate that the erosion of uplifted basement rocks exposed in horst blocks provided the sediment sources for the syn-extensional deposition of uppermost Jurassic–lowermost Cretaceous conglomerates and Lower Cretaceous siliciclastic red beds within a continental retro-arc basin during initiation of the ‘Neo-Tethys 2?. The polyphase tectonic reactivation along the principal fault of the study area controlled the syn- and post-extensional tectonostratigraphic evolution that reflect the corresponding mechanical decoupling/coupling along the northern Neo-Tethyan plate margin.     

Discussion of “Consequences of dike breaches and dike overflow in a bifurcating river system” by Anouk Bomers,Ralph M. J. Schielen and Suzanne J. M. H. Hulscher

In this discussion, the authors will point out that even if Bomers et al. (Nat Hazards 97:309–334, 2019) tackle an important problem, ignoring the uncertainties related to the roughness coefficients, Manning coefficients, the downstream boundary and most importantly the errors of the chosen software, HEC-RAS, are serious shortcomings of their study.

     

Three‐dimensional nonlinear finite element analysis of pile groups in saturated porous media using a new transmitting boundary

The dynamic behaviour of pile groups subjected to an earthquake base shaking is analysed. An analysis is formulated in the time domain and the effects of material nonlinearity of soil, pile–soil–pile kinematic interaction and the superstructure–foundation inertial interaction on seismic response are investigated. Prediction of response of pile group–soil system during a large earthquake requires consideration of various aspects such as the nonlinear and elasto‐plastic behaviour of soil, pore water pressure generation in soil, radiation of energy away from the pile, etc. A fully explicit dynamic finite element scheme is developed for saturated porous media, based on the extension of the original formulation by Biot having solid displacement (u) and relative fluid displacement (w) as primary variables (u–w formulation). All linear relative fluid acceleration terms are included in this formulation. A new three‐dimensional transmitting boundary that was developed in cartesian co‐ordinate system for dynamic response analysis of fluid‐saturated porous media is implemented to avoid wave reflections towards the structure. In contrast to traditional methods, this boundary is able to absorb surface waves as well as body waves. The pile–soil interaction problem is analysed and it is shown that the results from the fully coupled procedure, using the advanced transmitting boundary, compare reasonably well with centrifuge data. Copyright © 2007 John Wiley & Sons, Ltd.     

Numerical simulation of fully saturated porous materials

Fully coupled, porous solid–fluid formulation, implementation and related modeling and simulation issues are presented in this work. To this end, coupled dynamic field equations with u?p?U formulation are used to simulate pore fluid and soil skeleton (elastic–plastic porous solid) responses. Present formulation allows, among other features, for water accelerations to be taken into account. This proves to be useful in modeling dynamic interaction of media of different stiffnesses (as in soil–foundation–structure interaction). Fluid compressibility is also explicitly taken into account, thus allowing excursions into modeling of limited cases of non‐saturated porous media. In addition to these features, present formulation and implementation models in a realistic way the physical damping, which dissipates energy. In particular, the velocity proportional damping is appropriately modeled and simulated by taking into account the interaction of pore fluid and solid skeleton. Similarly, the displacement proportional damping is physically modeled through elastic–plastic processes in soil skeleton. An advanced material model for sand is used in present work and is discussed at some length. Also explored in this paper are the verification and validation issues related to fully coupled modeling and simulations of porous media. Illustrative examples describing the dynamical behavior of porous media (saturated soils) are presented. The verified and validated methods and material models are used to predict the behavior of level and sloping grounds subjected to seismic shaking. Copyright © 2008 John Wiley & Sons, Ltd.     

Limiting equilibrium method for slope/drilled shaft system

The use of drilled shafts to stabilize an unstable slope has been a widely accepted practice. There are two basic design and analysis issues involved: one is to determine the global factor of safety of the drilled shafts stabilized slope and the other one is to determine the design earth thrust on the drilled shafts for structural design of the shafts. In this paper, a limiting equilibrium method of slices based solution for calculating global factor of safety (FS) of a slope with the presence of a row of drilled shafts is developed. The arching mechanisms due to the presence of the drilled shafts on slope were taken into account by a load transfer factor. The method for calculating the net force applied to the drilled shaft from the soil mass was also developed. The interrelationships among the drilled shaft location on the slope, the load transfer factor, and the global FS of the slope/shaft system were derived utilizing the developed numerical closed‐form solution. An illustrative example is presented to elucidate the use of the solution in optimizing the location of the drilled shafts on slope to achieve the desired global factor of safety of the slope/shaft system. Copyright © 2009 John Wiley & Sons, Ltd.     

Mineral chemical composition and geodynamic significance of peridotites from Nain ophiolite,central Iran

The peridotites from north of the town of Nain in central Iran consist of clinopyroxene-bearing harzburgite and lherzolite with small lenses of dunite and chromitite pods. The lherzolite contains aluminous spinel with a Cr number (Cr^# = Cr/[Cr + Al]) of 0.17. The Cr number of spinels in harzburgite and chromitite is 0.38–0.42 and 0.62, respectively. This shows that the lherzolite and harzburgite resulted from <18% of partial melting of the source materials. The estimated temperature is 1100 ± 200 °C for peridotites, the estimated pressure is <15 ± 2.3 kbar for harzburgites and >16 ± 2.3 kbar for lherzolites and estimated fo₂ is 10^?1±0.5 for peridotites. Discriminant geochemical diagrams based on mineral chemistry of harzburgites indicate a supra-subduction zone (SSZ) to mid-oceanic ridge (MOR) setting for these rocks. On the basis of their Cr^#, the harzburgite and lherzolite spinels are analogous to those from abyssal peridotites and oceanic ophiolites, whereas the chromites in the chromitite (on the basis of Cr^# and boninitic nature of parental melts) resemble those from SSZ ophiolitic sequences. Therefore, the Nain ophiolite complex most likely originated in an oceanic crust related to supra-subduction zone, i.e. back arc basin. Field observations and mineral chemistry of the Nain peridotites, indicating the suture between the central Iran micro-continent (CIM) block and the Sanandaj–Sirjan zone, show that these peridotites mark the site of the Nain–Baft seaway, which opened with a slow rate of ocean-floor spreading behind the Mesozoic arc of the Sanandaj–Sirjan zone as a result of change of Neo Tethyan subduction régime during middle Cretaceous.     

Development in modeling cyclic loading of sands based on kinematic hardening

In this paper, there is presented an elastoplastic constitutive model to predict sandy soils behavior under monotonic and cyclic loadings. This model is based on an existing model (Cambou‐Jafari‐Sidoroff) that takes into account deviatoric and isotropic mechanisms of plasticity. The flow rule used in the deviatoric mechanism is non‐associated and a mixed hardening law controls the evolution of the yield surface. In this research the critical state surface and history surface, which separates the virgin and cyclic states in stress space, are defined. Kinematic hardening modulus and stress–dilatancy law for monotonic and cyclic loadings are effectively modified. With taking hardening modulus as a function of deviatoric and volumetric plastic strain and with defining the history surface and stress reversal, the model has the ability to predict the sandy soils' behavior. All of the model parameters have clear physical meanings and can be determined from usual laboratory tests. In order to validate the model, the results of homogeneous tests on Hostun and Toyoura sands are used. The results of validation show a good capability of the proposed model. Copyright © 2009 John Wiley & Sons, Ltd.     

List of forthcoming papers

List of forthcoming papers