A unified formulation of the piecewise exact method for inelastic seismic demand analysis including the P‐delta effect

The non‐linear analysis of single‐degree‐of‐freedom (SDOF) systems provides the essential background information for both strength‐based design and displacement‐based evaluation/design methodologies through the development of the inelastic response spectra. The recursive solution procedure called the piecewise exact method, which is efficiently used for the response analysis of linear SDOF systems, is re‐formulated in this paper in a unified format to analyse the non‐linear SDOF systems with multi‐linear hysteresis models. The unified formulation is also capable of handling the P‐delta effect, which generally involves the negative post‐yield stiffness of the hysteresis loops. The attractiveness of the method lies in the fact that it provides the exact solution when the loading time history is composed of piecewise linear segments, a condition that is perfectly satisfied for the earthquake excitation. Based on simple recursive relationships given for positive, negative and zero effective stiffnesses, the unified form of the piecewise exact method proves to be an extremely powerful and probably the best tool for the SDOF inelastic time‐history and response spectrum analysis including the P‐delta effect. A number of examples are presented to demonstrate the implementation of the method. Copyright © 2003 John Wiley & Sons, Ltd.     

Characterisation of the 1997 Vulcanian explosions of Soufrière Hills Volcano,Montserrat, by video analysis

The activity of Convention at Montserrat Soufrière Hills Volcano, Montserrat, during the period 1995–1999 included numerous violent explosions. Two major cycles of Vulcanian explosions occurred in 1997: a first of 13 explosions between 4 and 12 August and a second of 75 between 22 September and 21 October. The explosions were short-lived events lasting a few tens of seconds during which partial fountain collapse generated pyroclastic surges and pyroclastic flows, and buoyant plumes ascended 3–15 km into the atmosphere. Each explosion discharged on average 3×10⁵ m³ (dense-rock equivalent, DRE) of magma, draining the conduit to depths of 1–2 km. The paper focuses on the first few seconds of three explosions of the 75 that occurred in September/October 1997: 6 October 1997 at 17:50, 7 October 1997 at 16:02 and 9 October 1997 at 12:32. Physical parameters such as exit velocities, magmatic water contents and magma pressures at fragmentation are estimated by following and modelling the ascent of individual momentum-dominated finger jets visible on videos during the initial stages of each explosion. The model treats each finger jet as an incompressible flow sustained by a steady flux of gas and particles during the few seconds of ascent, and produces results that compare favourably with those using a multiphase compressible code run using similar eruptive parameters. Each explosion reveals a progressive increase in eruptive intensity with time, jet exit velocities increasing from 40 m s^–1 at the beginning of the explosion up to 140 m s^–1 after a few seconds. Modelling suggests that the first magma to exit was largely degassed, whereas that discharged after a few seconds contained up to 2 wt% water. Magma overpressures up to ~10 MPa are estimated to have existed in the conduit immediately prior to each explosion. Progressive increases in jet exit velocity with time over the first few seconds of each explosion provide direct evidence for strong pre-eruptive gradients in water content and magma pressure in the upper reaches (probably 100–500 m) of the conduit. Fountain collapse occurred during the first 10–20 s of each explosion because the discharging jets had bulk densities up to 100 times that of the atmosphere and were unable to entrain enough air to become buoyant. Such high eruptive densities were due to the presence of partially degassed magma in the conduit.Editorial responsibility: A. Woods     

Landslide susceptibility analysis using GIS and artificial neural network

The purpose of this study is to develop landslide susceptibility analysis techniques using an arti?cial neural network and to apply the newly developed techniques to the study area of Yongin in Korea. Landslide locations were identi?ed in the study area from interpretation of aerial photographs, ?eld survey data, and a spatial database of the topography, soil type and timber cover. The landslide‐related factors (slope, curvature, soil texture, soil drainage, soil effective thickness, timber age, and timber diameter) were extracted from the spatial database. Using those factors, landslide susceptibility was analysed by arti?cial neural network methods. The landslide susceptibility index was calculated by the back‐propagation method, which is a type of arti?cial neural network method, and the susceptibility map was made with a geographic information system (GIS) program. The results of the landslide susceptibility analysis were veri?ed using landslide location data. The validation results showed satisfactory agreement between the susceptibility map and the existing data on landslide location. A GIS was used to ef?ciently analyse the vast amount of data, and an arti?cial neural network to be an effective tool to maintain precision and accuracy. The results can be used to reduce hazards associated with landslides and to plan land use and construction. Copyright © 2003 John Wiley & Sons, Ltd.     

Computer Simulation of Shaped Charge Detonation Using Meshless Particle Method

Shaped charges are widely applied in the defense industry and civil engineering for penetrating hardened targets. Simulations of shaped charges are generally based on traditional grid-based numerical methods. Different from the existing literature, the SPH method is applied in this paper to simulate the detonation process of shaped charge explosive. The combination of the meshless, Lagrangian and particle nature inherent in the SPH method is very attractive in dealing with large deformations and large inhomogeneities in extremely transient situations such as detonation process. The SPH method and some related numerical aspects are briefly discussed first, followed by applications in simulating shaped charge detonations. The effects of different detonation cavity angles and different charge head lengths are investigated. The results demonstrate that major shaped charge detonation phenomena are captured with revealing observations.     

Spatial probabilistic modeling of slope failure using an integrated GIS Monte Carlo simulation approach

Spatial probabilistic modeling of slope failure using a combined Geographic Information System (GIS), infinite-slope stability model and Monte Carlo simulation approach is proposed and applied in the landslide-prone area of Sasebo city, southern Japan. A digital elevation model (DEM) for the study area has been created at a scale of 1/2500. Calculated results of slope angle and slope aspect derived from the DEM are discussed. Through the spatial interpolation of the identified stream network, the thickness distribution of the colluvium above Tertiary strata is determined with precision. Finally, by integrating an infinite-slope stability model and Monte Carlo simulation with GIS, and applying spatial processing, a slope failure probability distribution map is obtained for the case of both low and high water levels.     

Famennian microbial reef facies, Napier and Oscar Ranges, Canning Basin, western Australia

Following the Frasnian–Famennian mass extinction, which eliminated most skeletal reef-building fauna, the early Famennian reefs of the Canning Basin were constructed primarily by reef-framework microbial communities. In the Napier and Oscar Ranges, the Famennian reef complexes had high-energy, reef-flat depositional environments on a reef-rimmed platform that transitioned into low-energy, deep-water reefs growing in excess of 50 m below sea level. High-energy, reef-flat depositional environments contain doming fenestral stromatolites that grade into porous thrombolites and are associated with coarse-grained sandstones and grainstones. The reef-margin subfacies contains mounds of microdigitate thrombolites, which are more delicate than the reef-flat thrombolites and locally contain abundant red algae, Girvanella, renalcids and sediment-filled tubes. Within the thrombolites, the red algae are in upright growth positions, suggesting that the thrombolites are largely composed of carbonate that precipitated in situ. Reefal-slope environments are dominated by Wetheredella and Rothpletzella with locally abundant Girvanella, renalcids and Uralinella. In reefal-slope strata, delicate fans and microdigitate stromatolites of Wetheredella and Rothpletzella are often oriented horizontal or diagonal to bedding and are interpreted as syndepositionally toppled over. Most mesoscale microbial community structures contain several species of microbial fossils, and no single microbial species appears to have controlled the morphology of the community structure. Therefore, the depositional environment must have determined the distribution and morphology of the stromatolites, thrombolites and other microbial community structures. The adaptability of microbial communities to various reef environments allowed them to fill ecological niches opportunistically after the Frasnian–Famennian mass extinction.     

Hydrochemical and isotopic characteristics of groundwater in the Souss Upstream Basin, southwestern Morocco

Heterogeneous shallow Plio-Quaternary formations of the Souss Plain represent the most important aquifer in southern High Atlas Mountains in Morocco. The present work was conducted in the Souss Upstream Basin to identify the chemical characteristics and the origin of groundwater in an aquifer under semi-arid climate. Isotopic and hydrochemical compositions combined with geological and hydrogeological data were used for this purpose. The total dissolved solids vary from 239 to 997 mg l⁻¹, and the following groundwater types are recognized: Ca²⁺–Mg²⁺–HCO₃⁻, Ca²⁺–Mg²⁺–SO₄²⁻ and Ca²⁺–Mg²⁺–Cl⁻. The groundwater is saturated and slightly supersaturated with respect to carbonate minerals and undersaturated with respect to evaporite minerals, which means that the groundwater composition is largely controlled by the dissolution of carbonate rocks known in the basin. The isotopic contents of groundwaters ranged from −8‰ to −5.2‰ for δ¹⁸O, from −52‰ to −34‰ for δD, and from 0 to 5.5 TU for tritium. The hydrogen (δD) and oxygen (δ¹⁸O) isotope signatures reveal a significant infiltration before evaporation takes place, indicating a major recharge directly from fractures in the crystalline and limestone formations of Atlas Mountains (above 800 m a.s.l.) and infiltration of surface water in the alluvial cones at the border of the Atlas basins. The very low tritium values suggest that the groundwater recharge follows a long flow path and a mixing between old and modern water is shown. However, a slight evaporation effect is noted in the southern part of the basin close to the Anti-Atlas Mountains.