A novel predictor–corrector algorithm for sub‐structure pseudo‐dynamic testing

A new predictor–corrector (P–C) method for multi‐site sub‐structure pseudo‐dynamic (PSD) test is proposed. This method is a mixed time integration method in which computational components separable from experimental components are solved by implicit time integration method (Newmark β method). The experiments are performed quasi‐statically based on explicit prediction of displacement. The proposed P–C method has an important advantage as it does not require the determination of the initial stiffness values of experimental components and is thus suitable for representing elastic and inelastic systems. A parameter relating to quality of displacement prediction at boundaries nodes is introduced. This parameter is determined such that P–C method can be applicable to many practical problems. Error‐propagation characteristics of P–C method are also presented. A series of examples including linear and non‐linear soil–foundation–structure interaction problem demonstrate the performance of the proposed method. Copyright © 2005 John Wiley & Sons, Ltd.     

Pseudo‐dynamic testing with substructuring at the ELSA Laboratory

The paper describes the distinctive features of the pseudo‐dynamic test method as implemented at the ELSA reaction‐wall facility. Both hardware and software aspects are considered. Particular attention is devoted to the digital control system and to a coupled numerical–experimental substructuring technique allowing realistic earthquake testing of very large structures. Mathematical and implementation details corresponding to this testing technique are given for both synchronous and asynchronous input motion. Selected test results illustrate the advantages of the presented features. Copyright © 2000 John Wiley & Sons, Ltd.     

A temporal stable isotopic (δ18O, δD,d‐excess) comparison in glacier meltwater streams,Taylor Valley,Antarctica

In this paper, we describe the importance of hyporheic dynamics within Andersen Creek and Von Guerard Stream, Taylor Valley, Antarctica, from the 2010–2011 melt season using natural tracers. Water collection started at flow onset and continued, with weekly hyporheic‐zone sampling. The water δ¹⁸O and δD values were isotopically lighter in the beginning and heavier later in the season. D‐excess measurements were used as an indicator of mixing because an evaporative signature was evident and distinguishable between 2 primary end‐members (glacier meltwater and hyporheic zone). Hyporheic‐zone influence on the channel water was variable with a strong control on streamwater chemistry, except at highest discharges. This work supports previous research indicating that Von Guerard Stream has a large, widespread hyporheic zone that varies in size with time and discharge. Andersen Creek, with a smaller hyporheic zone, displayed hyporheic‐zone solute interaction through the influence from subsurface hypersaline flow. Overall, the evolution of Taylor Valley hyporheic‐zone hydrology is described seasonally. In mid‐December, the hyporheic zone is a dynamic system exchanging with the glacier meltwater in the channel, and with diminishing flow in January, the hyporheic zone drains back into the channel flow also impacting stream chemistry. This work adds new information on the role of hyporheic zone–stream interaction in these glacier meltwater streams.     

Direct finite element method for nonlinear earthquake analysis of 3‐dimensional semi‐unbounded dam–water–foundation rock systems

A direct finite element method for nonlinear earthquake analysis of 2‐dimensional dam–water–foundation rock systems has recently been presented. The analysis procedure uses standard viscous‐damper absorbing boundaries to model the semi‐unbounded foundation‐rock and fluid domains and specifies the seismic input as effective earthquake forces at these boundaries. Presented in this paper is a generalization of the direct finite element method with viscous‐damper boundaries to 3‐dimensional dam–water–foundation rock systems. Step‐by‐step procedures for determining the effective earthquake forces starting from a ground motion specified at a control point on the foundation‐rock surface is developed, and several numerical examples are computed and compared with independent benchmark solutions to demonstrate the effectiveness of the analysis procedure for modeling 3‐dimensional systems.     

A fuzzy–random analysis model for seismic performance of framed structures incorporating structural and non‐structural damage

Past earthquake experiences indicate that most buildings designed in accordance with modern seismic design codes could survive moderate‐to‐strong earthquakes; however, the financial loss due to repairing cost and the subsequent business interruption can be unacceptable. Designing building structures to meet desired performance targets has become a clear direction in future seismic design practice. As a matter of fact, the performance of buildings is affected by structural as well as non‐structural components, and involves numerous uncertainties. Therefore, appropriate probabilistic approach taking into account structural and non‐structural damages is required. This paper presents a fuzzy–random model for the performance reliability analysis of RC framed structures considering both structural and non‐structural damages. The limit state for each performance level is defined as an interval of inter‐storey drift ratios concerning, respectively, the non‐structural and structural damage with a membership function, while the relative importance of the two aspects is reflected through the use of an appropriate cost function. To illustrate the methodology, herein the non‐structural damage is represented by infill masonry walls. The probabilistic drift limits for RC components and masonry walls from the associated studies are employed to facilitate the demonstration of the proposed model in an example case study. The results are compared with those obtained using classical reliability model based on single‐threshold performance definition. The proposed model provides a good basis for incorporating different aspects into the performance assessment of a building system. Copyright © 2005 John Wiley & Sons, Ltd.     

Estimation of uncertainty propagation in flood inundation mapping using a 1‐D hydraulic model

A need for more accurate flood inundation maps has recently arisen because of the increasing frequency and extremity of flood events. The accuracy of flood inundation maps is determined by the uncertainty propagated from all of the variables involved in the overall process of flood inundation modelling. Despite our advanced understanding of flood progression, it is impossible to eliminate the uncertainty because of the constraints involving cost, time, knowledge, and technology. Nevertheless, uncertainty analysis in flood inundation mapping can provide useful information for flood risk management. The twin objectives of this study were firstly to estimate the propagated uncertainty rates of key variables in flood inundation mapping by using the first‐order approximation method and secondly to evaluate the relative sensitivities of the model variables by using the Hornberger–Spear–Young (HSY) method. Monte Carlo simulations using the Hydrologic Engineering Center's River Analysis System and triangle‐based interpolation were performed to investigate the uncertainty arising from discharge, topography, and Manning's n in the East Fork of the White River near Seymour, Indiana, and in Strouds Creek in Orange County, North Carolina. We found that the uncertainty of a single variable is propagated differently to the flood inundation area depending on the effects of other variables in the overall process. The uncertainty was linearly/nonlinearly propagated corresponding to valley shapes of the reaches. In addition, the HSY sensitivity analysis revealed the topography of Seymour reach and the discharge of Strouds Creek to be major contributors to the change of flood inundation area. Copyright © 2014 John Wiley & Sons, Ltd.     

Study on effects of damping in laminated rubber bearings on seismic responses for a ⅛ scale isolated test structure

The effects of damping in various laminated rubber bearings (LRB) on the seismic response of a ?‐scale isolated test structure are investigated by shaking table tests and seismic response analyses. A series of shaking table tests of the structure were performed for a fixed base design and for a base isolation design. Two different types of LRB were used: natural rubber bearings (NRB) and lead rubber bearings (LLRB). Three different designs for the LLRB were tested; each design had a different diameter of lead plug, and thus, different damping values. Artificial time histories of peak ground acceleration 0.4g were used in both the tests and the analyses. In both shaking table tests and analyses, as expected, the acceleration responses of the seismically isolated test structure were considerably reduced. However, the shear displacement at the isolators was increased. To reduce the shear displacement in the isolators, the diameter of the lead plug in the LLRB had to be enlarged to increase isolator damping by more than 24%. This caused the isolator stiffness to increase, and resulted in amplifying the floor acceleration response spectra of the isolated test structure in the higher frequency ranges with a monotonic reduction of isolator shear displacement. Copyright © 2002 John Wiley & Sons, Ltd.     

Hybrid formulation of operator splitting (OS) and Newmark‐β methods for collaborative structural analysis (CSA)

A collaborative structural analysis (CSA) system is developed, which is capable of performing highly sophisticated structural analyses utilizing beneficial features of existing individual structural analysis programs. In the system, the global equations of motion for the overall structural system are formulated in the host program. Some substructures, whose behaviors are relatively simple, are directly solved in the host program, whereas those having complex behavior are analyzed by the station programs. A time‐consuming static condensation procedure is needed for the substructures analyzed by the station programs if adopting an implicit integration scheme. The operator splitting (OS) method, which does not require tangential stiffness, can be used to improve the system efficiency. To this end, a hybrid formulation of the Newmark‐β and OS methods is proposed, and a CSA scheme based on the hybrid formulation is developed. In the CSA system adopting the hybrid formulation, the degrees of freedom whose tangential stiffness are unavailable are formulated by the OS method, whereas the rest are still formulated by the commonly used Newmark‐β method. Using the system, analyses of a three‐story‐braced steel moment‐resisting frame are conducted. In the analyses, the column bases are analyzed using the commercial finite element method software ABAQUS, and the remaining structural elements are analyzed using a frame analysis program called NETLYS. Results suggest that the hybrid formulation is very effective for the CSA system. Copyright © 2008 John Wiley & Sons, Ltd.     

Physical constraints on c13 and δ for transversely isotropic hydrocarbon source rocks

Based on the theory of anisotropic elasticity and observation of static mechanic measurement of transversely isotropic hydrocarbon source rocks or rock‐like materials, we reasoned that one of the three principal Poisson's ratios of transversely isotropic hydrocarbon source rocks should always be greater than the other two and they should be generally positive. From these relations, we derived tight physical constraints on c₁₃, Thomsen parameter δ, and anellipticity parameter η. Some of the published data from laboratory velocity anisotropy measurement are lying outside of the constraints. We analysed that they are primarily caused by substantial uncertainty associated with the oblique velocity measurement. These physical constraints will be useful for our understanding of Thomsen parameter δ, data quality checking, and predicting δ from measurements perpendicular and parallel to the symmetrical axis of transversely isotropic medium. The physical constraints should also have potential application in anisotropic seismic data processing.     

Characteristics of δ13CDIC in lakes on the Tibetan Plateau and its implications for the carbon cycle

Dissolved inorganic carbon isotope (δ¹³C_DIC) is an important tool to reveal the carbon cycle in lake systems. However, there are only few studies focusing on the spatial variation of δ¹³C_DIC of closed lakes. Here we analyze the characteristics of δ¹³C_DIC of 24 sampled lakes (mainly closed lakes) across the Qiangtang Plateau (QTP) and identify the driving factors for its spatial variation. The δ¹³C_DIC value of these observed lakes varies in the range of ? 15·0 to 3·2‰, with an average value of ? 1·2‰. The δ¹³C_DIC value of closed lakes is close to the atmospheric isotopic equilibrium value, much higher than that in rivers and freshwater lakes reported before. The high δ¹³C_DIC value of closed lakes is mainly attributed to the significant contribution of carbonate weathering in the catchment and the evasion of dissolved CO₂ induced by the strong evaporation of lake water. The δ¹³C_DIC value of closed lakes has a logarithmic correlation with water chemistry (TDS, DIC and pCO₂), also suggesting that the evapo‐concentration of lake water can influence the δ¹³C_DIC value. The δ¹³C_DIC value shows two opposite logarithmic correlations with lake size depending on the δ¹³C_DIC range. This study suggests that the δ¹³C in carbonates in lacustrine sediments can be taken as an indicator of lake volume variation in closed lakes on QTP. Copyright © 2011 John Wiley & Sons, Ltd.     

Characterization of surface and ground water δ18O seasonal variation and its use for estimating groundwater residence times

¹⁸O is an ideal tracer for characterizing hydrological processes because it can be reliably measured in several watershed hydrological compartments. Here, we present multiyear isotopic data, i.e. ¹⁸O variations (δ¹⁸O), for precipitation inputs, surface water and groundwater in the Shingobee River Headwaters Area (SRHA), a well‐instrumented research catchment in north‐central Minnesota. SRHA surface waters exhibit δ¹⁸O seasonal variations similar to those of groundwaters, and seasonal δ¹⁸O variations plotted versus time fit seasonal sine functions. These seasonal δ¹⁸O variations were interpreted to estimate surface water and groundwater mean residence times (MRTs) at sampling locations near topographically closed‐basin lakes. MRT variations of about 1 to 16 years have been estimated over an area covering about 9 km² from the basin boundary to the most downgradient well. Estimated MRT error (±0·3 to ±0·7 years) is small for short MRTs and is much larger (±10 years) for a well with an MRT (16 years) near the limit of the method. Groundwater transit time estimates based on Darcy's law, tritium content, and the seasonal δ¹⁸O amplitude approach appear to be consistent within the limits of each method. The results from this study suggest that use of the δ¹⁸O seasonal variation method to determine MRTs can help assess groundwater recharge areas in small headwaters catchments. Copyright © 2005 John Wiley & Sons, Ltd.     

Experimental multi‐level seismic performance assessment of 3D RC frame designed for damage avoidance

This paper experimentally investigates the application of damage avoidance design (DAD) philosophy to moment‐resisting frames with particular emphasis on detailing of rocking interfaces. An 80% scale three‐dimensional rocking beam–column joint sub‐assembly designed and detailed based on damage avoidance principles is constructed and tested. Incremental dynamic analysis is used for selecting ground motion records to be applied to the sub‐assembly for conducting a multi‐level seismic performance assessment (MSPA). Analyses are conducted to obtain displacement demands due to the selected near‐ and medium‐field ground motions that represent different levels of seismic hazard. Thus, predicted displacement time histories are applied to the sub‐assembly for conducting quasi‐earthquake displacement tests. The sub‐assembly performed well reaching drifts up to 4.7% with only minor spalling occurring at rocking beam interfaces and minor flexural cracks in beams. Yielding of post‐tensioning threaded bars occurred, but the sub‐assembly did not collapse. The externally attached energy dissipators provided large hysteretic dissipation during large drift cycles. The sub‐assembly satisfied all three seismic performance requirements, thereby verifying the superior performance of the DAD philosophy. Copyright © 2007 John Wiley & Sons, Ltd.     

Effects of 17β-estradiol on the reproduction of Java-medaka (Oryzias javanicus), a new test fish species

17β-estradiol (E2) is a female hormone which is known to be one of the strongest estrogenic chemicals in the environment. The present study investigated the effects of E2 on the reproduction of the estuarine fish, Java-medaka (Oryzias javanicus). Starting from the embryonic stages, Java-medaka was exposed for 6 months to 9.5, 16, 68, 159 and 243 ng/l of E2. The fecundity of Java-medaka exposed to E2 levels >16 ng/l was significantly lower than that of the control. The appearance of secondary sexual characteristics seemed to be inhibited by exposure to 159 and 243 ng/l E2. Vitellogenin concentrations in the liver of male fish exposed to 68, 159 and 243 ng/l were significantly higher than that of the control. Thirty-three and sixty percent of male fishes exposed to 159 and 243 ng/l, respectively, had testis-ova. These results suggest that reproduction in the Java-medaka, especially the male fish, could be affected by exposure to E2 concentrations greater than 16 ng/l.     

A multi‐geochronological study of the Hakusan volcano,central Japan

The Hakusan volcano, central Japan, is located in a region where two subducting plates (the Pacific Plate and the Philippine Sea Plate) overlap near the junction of four plates adjacent to the Japanese Islands (the Pacific Plate, the Philippine Sea Plate, the Eurasia Plate, and the North American Plate). The Hakusan volcano consists of products from four major volcanic episodes: Kagamuro, Ko‐hakusan, and Shin‐Hakusan I and II. To date the eruption events of the Hakusan volcano we applied thermoluminescence and fission track methods. ²³⁸U(²³⁴U)–²³⁰Th disequilibrium and ²⁰⁶Pb/²³⁸U methods were applied to date the zircon crystallization ages for estimating the magma residence time before the eruptions. The eruption ages we obtained are ca 250 ka for Kagamuro, ca 100 ka and ca 60 ka for Ko‐Hakusan, ca 50 ka for Shin‐Hakusan I, and <10 ka for Shin‐Hakusan II. They are concordant with previous reports based on K–Ar dating. Some of the pyroclastic rocks, possibly originating from Shin‐Hakusan II activities, are dated to be ca 36 ka or 50 ka, and belong to the Shin‐Hakusan I activity. The zircon crystallization ages show several clusters prior to eruption. The magma residence time was estimated for each volcanic activity by comparing the major crystallization events and eruption ages, and we found a gradual decrease from ca. 500 ky for the Kagamuro activity to ca. 5 ky for the Shin‐Hakusan II activity. This decrease in residence time may be responsible for the decrease in volume of erupted material estimated from the current topography of the region. The scale of volcanic activity, which was deduced from the number of crystallized zircons, is more or less constant throughout the Hakusan volcanic activity. Therefore, the decrease in magma residence time is most likely the result of stress field change.     

A new family of generalized‐α time integration algorithms without overshoot for structural dynamics

A new family of generalized‐α (G‐α) algorithm without overshoot is presented by introducing seven free parameters into the single‐step three‐stage formulation for solution of structural dynamic problems. It is proved through finite difference analysis that these algorithms are unconditionally stable, second‐order accurate and numerical dissipation controllable. The comparison of the new G‐α algorithms with the commonly used G‐α algorithms shows that the newly developed algorithms have the advantage of eliminating the overshooting characteristics exhibited by the commonly used algorithms while their excellent property of dissipation is preserved. The numerical simulation results obtained using a single‐degree‐of‐freedom system and a two‐degree‐of‐freedom system to represent the character of typical large systems coincide well with the results of theoretical analyses. Copyright © 2008 John Wiley & Sons, Ltd.     

On‐line hybrid test method for evaluating the performance of structural details to failure

A test environment to evaluate the seismic performance of gusset plate connections intended for steel braced frames is proposed. The developed test method combines the substructuring techniques with finite element analysis methods in an on‐line hybrid scheme. Numerical substructure analysis is conducted on bracing members, while bracing connections are treated as experimental substructures. A force‐displacement combined control imposed with the aid of 2 jacks ensures physical continuity between the analysis and test. The rotational behavior of gusset plate connections subjected to large inelasticity and varying axial loading until fracture is investigated. Two gusset plate details were designed and tested to verify the efficiency of the proposed method. The test method is rational, and smooth operation is achieved. The test results revealed the advantage of the developed on‐line hybrid test method in exploring the ultimate capacity of bracing connections.     

Combining isotopic tracers (222Rn and δ13C) for improved modelling of groundwater discharge to small rivers

In regions where aquifers sustain rivers, the location and quantification of groundwater discharge to surface water are important to prevent pollution hazards, to quantify and predict low flows and to manage water supplies. ²²²Rn is commonly used to determine groundwater discharge to rivers. However, using this isotopic tracer is challenging because of the high diffusion capacity of ²²²Rn in open water. This study illustrates how a combination of isotopic tracers can contribute to an enhanced understanding of groundwater discharge patterns in small rivers. The aim of this paper is to combine ²²²Rn and δ¹³C_DIC to better constrain the physical parameters related to the degassing process of these tracers in rivers. The Hallue River (northern France) was targeted for this study because it is sustained almost exclusively by a fractured chalk aquifer. The isotopes ²²²Rn, δ¹³C_DIC, δ²H and δ¹⁸O were analysed along with other natural geochemical tracers. A mass balance model was used to simulate ²²²Rn and δ¹³C_DIC. The results of δ²H and δ¹⁸O analyses prove that evaporation did not occur in the river. The calibration of a numerical model to reproduce ²²²Rn and δ¹³C_DIC provides a best‐fit diffusive layer thickness of 3.21 × 10^?5 m. This approach is particularly useful for small rivers flowing over carbonate aquifers with high groundwater DIC where the evolution of river DIC reflects the competing processes of groundwater inflow and CO₂ degassing. This approach provides a means to evaluate groundwater discharge in small ungauged rivers. Copyright © 2014 John Wiley & Sons, Ltd.