Nonlinear modeling of cyclic response of RC beam–column joints reinforced by plain bars

Experience of previous earthquakes shows that a considerable portion of buildings reinforced with plain bars sustain relatively large damages especially at the beam–column joints where the damages are mostly caused by either diagonal shear cracks or intersectional cracks caused by bar slippage. While previous works mainly focus on shear failure mode, in this study, the emphasis is placed on slip based cracks as the dominant failure mode. A systematic procedure is introduced to predict the dominant failure mode at the joint which is based on the dimensional properties, reinforcement details, and axial and shear load at the joint. In addition, a relatively simple and efficient nonlinear model is proposed to simulate pre- and post-elastic behavior of the joints which fail under bar slippage mode. In this model, beam and column components are represented by linear elastic elements, dimensions of the joint panel are defined by rigid elements, and effect of slip is taken into account by a nonlinear rotational spring at the end of the beam. The proposed method is validated by experimental results for both internal and external joints .     

Winkler model for dynamic response of composite caisson–piles foundations: Seismic response

A simplified method with a dynamic Winkler model to study the seismic response of composite caisson–piles foundations (CCPF¹) is developed. Firstly, with the dynamic Winkler model, the kinematic response of the CCPF subjected to vertically propagating seismic S-wave is analyzed by coupling the responses of caisson part and pile part. Secondly, a simplified model for the foundation–structure system is created with the structure simplified as a lumped mass connected to the foundation with an elastic column, and through the Fast Fourier Transformation (FFT) this model is enabled to solve transient seismic problems. Thirdly, the proposed method for the seismic response of CCPF-structure systems is verified by comparison against 3D dynamic finite element simulation, in which the Domain Reduction Method (DRM²) is utilized. Lastly, the mechanism and significance of adding piles in improving the earthquake resistance of the foundation and structure is analyzed through an example with different soil conditions. Discovered in this study is that adding piles under the caisson is an efficient way to increase seismic resistant capability of the soil–foundation–structure system, and the main mechanism of that is the elimination of the pseudo-resonance.     

Nonlinear extensions of a fractal–multifractal approach for environmental modeling

We present the extension of a deterministic fractal geometric procedure aimed at representing the complexity of patterns encountered in environmental applications. The procedure, which is based on transformations of multifractal distributions via fractal functions, is extended through the introduction of nonlinear perturbations in the generating iterated linear maps. We demonstrate, by means of various simulations based on changes in parameters, that the nonlinear perturbations generate yet a richer collection of interesting patterns, as reflected by their overall shapes and their statistical and multifractal properties. It is shown that the nonlinear extensions yield structures that closely resemble complex hydrologic spatio-temporal datasets, such as rainfall and runoff time series, and width-functions of river networks. The implications of this nonlinear approach for environmental modeling and prediction are discussed.     

An efficient fiber beam-column element considering flexure–shear interaction and anchorage bond-slip effect for cyclic analysis of RC structures

In this paper, a fiber beam-column element considering flexure–shear interaction and bond-slip effect is developed for cyclic analysis of reinforced concrete (RC) structures. The element is based on conventional displacement-based Timoshenko beam theory, where the transverse shear deformation is included, and adopts the fiber model to describe the section force–deformation behavior. In the fiber model, shear deformation is assumed to be uniformly distributed along the section and is only resisted by concrete, thus the multi-dimensional concrete damage model is used for concrete fibers and therefore flexure–shear interaction is reflected naturally at the material level. Meanwhile, to account for the significant bond-slip effect at critical regions, the anchorage slip of bars at these regions is analytically derived. Then it is used to modify the uniaxial stress–strain model for steel fibers by assuming that the total strain can be treated as the sum of the bar deformation and anchorage slip, therefore the bond-slip effect is implicitly but simply represented. To validate the proposed element, a series of RC member and structure tests under cyclic loading are simulated. The results indicate that the proposed element can predict cyclic responses of RC structures, and can be used as a reliable tool for analysis of RC structures.     

Simplified seismic analysis of soil–well–pier system for bridges

Seismic analysis of soil–well–pier system was carried out using three different approaches to evaluate their comparative performance and associated complexities. These approaches were (a) two-dimensional nonlinear (2D-NL), (b) two-dimensional equivalent-linear (2D-EqL), and (c) one-dimensional spring–dashpot (1D). Soil was modeled as 2D plane-strain elements in the 2D-NL and 2D-EqL approaches, and as springs and dashpots in the 1D approach. Nonlinear behavior of soil was captured rigorously in the 2D-NL approach and approximately in the remaining two approaches. Results of the two approximate analyses (i.e., 2D-EqL and 1D) were compared with those of the 2D-NL analysis with the objective to assess suitability of approximate analysis for practical purposes. In the 1D approach, several combinations of Novak's and Veletsos' springs were used to come up with a simplified 1D model using three types of spring–dashpots. The proposed model estimates the displacement and force resultants relatively better than the other 1D models available in literature.     

A simplified approximation for seismic analysis of silo–bulk material system

In order to estimate the distribution, as well as the magnitude, of dynamic material pressures on ground-supported silos a simplified seismic analysis procedure was utilized. The seismic analysis of silos can be complex, as the evaluation of several parameters must be taken into consideration, including the properties of bulk materials used and how the bulk materials and silo wall are joined together. It is therefore useful to develop an analytical approximation in order to better assess results. In addition to a simplified model for the seismic analysis of a silo–bulk material system being utilized, a three-dimensional finite element model was also incorporated. Using the finite element method, a more realistic representation of the structure is possible. Moreover, the finite element method also takes into consideration contact problems between the bulk material and the silo wall, which results in easier analyses. Both a squat and a slender silo were selected for this study. The results obtained in the study of selected examples were compared with those findings obtained via EN1998-4. Modified Veletsos and Younan approximations, which are commonly used for the analysis of grain silos, were also used. Results and analysis concluded that the proposed analytical model provided, overall, a good outcome, especially in regards to the analysis of dynamic material pressure. It should be noted that using the analytical method as proposed in Eurocode, the dynamic material pressure for squat silos can be underestimated, but the results for slender silos are stronger.     

Bending–buckling interaction as a failure mechanism of piles in liquefiable soils

This paper investigates the importance of bending–buckling interaction in seismic design of piles in liquefiable soils using numerical techniques. A pseudo-static analysis has been performed using a well documented case history, where the pile–soil interaction is modelled as a beam on nonlinear winkler foundation (BNWF). Six possible analytical methods, three force based and three displacement based, are performed in which the pile is subjected to both lateral and axial load. Three out of six analysis cases did not predict the failure of the piles when analysed only for bending (i.e., lateral loads only). The buckling analysis showed that the pile was also safe against pure buckling during full liquefaction. Further, two out of those three cases which did not predict failure in bending were reanalysed for bending–buckling interaction (i.e., lateral and axial loads acting simultaneously). These combined analyses showed a more realistic behaviour of pile response and did predict the pile failure. Hence, it can be concluded that if a pile is designed for bending and buckling criteria separately and safe for these individual design criteria, it may fail due to their combined effect.     

Using GPS–GLONASS–GALILEO data and IRI modeling for ionospheric calibration of radio telescopes and radio interferometers

VHF and HF radio signals are widely used to observe the Sun and pulsars. Nowadays, large low-frequency radio astronomical arrays (LOFAR, 30–240 MHz; MIRA, 80–300 MHz) are being constructed to record radiation of pulsars at the maximum distance. registration of the solar radio emission intensity at fixed frequencies and in the spectral VHF band is very important along with other methods of monitoring of coronal mass ejections. Interpreting radio astronomical data is known to be necessary to take into account possible distortions of these signals in the Earth ionosphere. However, in contrast to modern navigation systems (Global Position System (GPS), GLObal NAvigation Satellite System (GLONASS), GALILEO), in which a very accurate reconstruction of ionosphere parameters is a built-in function, in present-day radio astronomy a retrieve of ionosphere transfer characteristics has not been appropriately worked out. This collides with increasing requirements to accuracy of the analysis of a radio emission amplitude profile and to the angular and polarizing resolution of radio telescopes of new generation (LOFAR, SKA, etc.). We have developed a method and software to calculate the ionosphere rotation measure (RM) and dispersion measure (DM). We used the ionosphere model IRI-2001, magnetic field model IGRF-10, and the ionosphere total electron content values obtained from GPS measurements. The obtained values of DM and RM were recalculated into characteristics of the phase delay, Faraday amplitude modulation, and polarization changes. We calculated ones for different levels of geomagnetic activity as well as different angular positions of radio sources. Our main idea is to use a signal of navigation satellites (GPS, GLONASS, GALILEO) as a testing signal from a “reference” source located at minimal angle distance from a source studied. Our project allows development of methods and systems of ADAPTIVE RADIO ASTRONOMY, adaptive to the non-uniform and non-stationary ionosphere, by analogy with known systems of adaptive optics intended to adapt optical telescopes to varying conditions of the optically non-uniform and non-stationary troposphere.     

Matrix–fracture transfer shape factor for modeling flow of a compressible fluid in dual-porosity media

The matrix–fracture transfer shape factor is one of the important parameters in the modeling of fluid flow in fractured porous media using a dual-porosity concept. Warren and Root [36] introduced the dual-porosity concept and suggested a relation for the shape factor. There is no general relationship for determining the shape factor for a single-phase flow of slightly compressible fluids. Therefore, different studies reported different values for this parameter, as an input into the flow models. Several investigations have been reported on the shape factor for slightly compressible fluids. However, the case of compressible fluids has not been investigated in the past. The focus of this study is, therefore, to find the shape factor for the single-phase flow of compressible fluids (gases) in fractured porous media. In this study, a model for the determination of the shape factor for compressible fluids is presented; and, the solution of nonlinear gas diffusivity equation is used to derive the shape factor. The integral method and the method of moments are used to solve the nonlinear governing equation by considering the pressure dependency of the viscosity and isothermal compressibility of the fluid. The approximate semi-analytical model for the shape factor presented in this study is verified using single-porosity, fine-grid, numerical simulations. The dependency of the shape factor on the gas specific gravity, pressure and temperature are also investigated. The theoretical analysis presented improves our understanding of fluid flow in fractured porous media. In addition, the developed matrix–fracture transfer shape factor can be used as an input for modeling flow of compressible fluids in dual-porosity systems, such as naturally fractured gas reservoirs, coalbed methane reservoirs and fractured tight gas reservoirs.     

An efficient approach for frequency-domain and time-domain hydrodynamic analysis of dam–reservoir systems

An efficient procedure is developed for the hydrodynamic analysis of dam–reservoir systems. The governing equations of hydrodynamic pressure in the frequency as well as time domain are derived in the framework of the scaled boundary finite element method. The water compressibility and absorption of reservoir sediments can be conveniently taken into consideration. By extending the reservoir to infinity with uniform cross-section, only the dam–reservoir interface needs to be discretized to model the fluid domain, and the hydrodynamic pressure in the stream direction is solved analytically. Several numerical examples including a gravity dam with an inclined upstream face and an arch dam with a reservoir of arbitrary cross-section are provided to demonstrate the computational efficiency and accuracy of the proposed method. Copyright © 2012 John Wiley & Sons, Ltd.     

Soil–structure interaction analysis of cable-stayed bridges for spatially varying ground motion components

In this study, it is intended to determine the effects of soil–structure interaction (SSI) and spatially varying ground motion on the dynamic characteristics of cable-stayed bridges. For this purpose, ground motion time histories are simulated for spatially varying ground motions, depending on its components of incoherence, wave-passage and site-response effects. The substructure method, which partitions the total soil–structure system into the structural system and the soil system, is used to treat the soil–structure interaction problem. To emphasize the relative importance of the spatial variability effects of earthquake ground motion, bridge responses are determined for the fixed base bridge model, which neglects the soil–structure interaction (no SSI) and for the bridge model including the soil–structure interaction (SSI). This parametric study concerning the relative importance of the soil–structure interaction and spatially varying ground motion shows that these effects should be considered in the dynamic analyses of cable-stayed bridges.     

Following the curve? Reviewing the physical basis of the SCS curve number method for estimating storm runoff

Less attention has been paid to runoff generation from semi-arid than from humid-temperate catchments. The SCS curve number approach is simple to apply and widely used, but lacks physical underpinning. Here output from a runoff generation models is compared with data from field measurements, making use of 11 years data from rainfall and runoff events at the Sierra de Enguera Soil Erosion Experimental Station in Eastern Spain. Runoff from natural rainfall events was monitored for 10 years on bare plots of 1–16 m length. The largest storm event was of 142 mm, generating runoff of up to 115 mm on the smallest plots. The model presented simulates overland storm flow on a sloping rough and unvegetated surface, representing an area of 320 × 320 m. Green-Ampt infiltration constants are randomly assigned to each cell in a 128 × 128 grid, and rectangular storms applied at a range of total amounts and intensities to simulate runoff at each transect across the area. A simple algebraic expression is developed to estimate total runoff and storage in terms of storm size and duration, and plot length, with parameters that reflect infiltration behaviour, and this expression is compared with the SCS curve number approach. For the very largest storms, both expressions converge asymptotically towards 100% runoff, but the revised expression greatly improves estimates of runoff from smaller events. Output of these simulations is compared with measured storm runoff data on bare runoff plots at the Sierra de Enguera experimental Station in SE Spain and gives further support to the proposed expression for storm runoff.     

High-precision modeling of tide-induced 3-D magnetic field and analysis of geomagnetic satellite orbit requirements

Ground-based magnetic observatories and geomagnetic satellites can observe the induced magnetic field generated by the motion of seawater containing sodium and chlorine ions.Calculating the three-dimensional (3-D) spatial distribution of tide-induced magnetic fields (TIMF) is crucial for inverting the electrical conductivity structure of the oceanic lithosphere.It also serves as an essential basis for designing optimal geomagnetic observatories and satellite orbits.However,existing methods for s...     

Convergence analyses of different modeling schemes for generalized Lippmann–Schwinger integral equation in piecewise heterogeneous media

For wave propagation simulation in piecewise heterogeneous media, Gaussian-elimination-based full-waveform solutions to the generalized Lippmann–Schwinger integral equation (GLSIE) are highly accurate, but involved with extremely time-consuming computations because of the very large size of the resulting boundary–volume integral equation matrix to be inverted. Several flexible approximations to the GLSIE are scaled in an iterative way to adapt numerical solutions to the smoothness of heterogeneous media in terms of incident wavelengths, with a great saving of computing time and memory. Among various typical iterative schemes to the GLSIE matrix, the generalized minimal residual method (GMRES) is an efficient approach to reduce the computational intensity to some degree. The most efficient approximation can be obtained using a Born series, as an alternative iterative solution, to both the boundary-scattering and volume-scattering waves, leading to the Born-series approximation (BSA) scheme and the improved Born-series approximation (IBSA) scheme. These iteration schemes are validated by dimensionless frequency responses to a heterogeneous semicircular alluvial valley, and then applied to a heterogeneous multilayered model by calculating synthetic seismograms to evaluate approximation accuracies. Numerical experiments, compared with the full-waveform numerical solution, indicate that the convergence rates of these methods decrease gradually with increasing velocity perturbations. The comparison also shows that the BSA scheme has a faster convergence than the GMRES method for velocity perturbations less than 10 percent, but converges slowly and even hardly achieves convergence for velocity perturbations greater than 15 percent. The IBSA scheme gives a superior performance over the other methods, with the least iterations to achieve the necessary convergence.     

Time–frequency analysis of VLF for seismic-ionospheric precursor detection: Evaluation of Zhao-Atlas-Marks and Hilbert-Huang Transforms

This work investigates the application of two post-processing methods of extracting spectra from VLF signals in order to detect disturbances that could be attributed to seismic-ionospheric precursory phenomena. Although precursory phenomena have been investigated in detail in past studies, a different application of time–frequency analysis methods may produce distinct patterns, which reveal disturbances in the VLF spectra received from stations that are in the propagation path over preparation zones, and also pinpoint disturbances that could be attributed to seismic-ionospheric precursors.To this purpose, three different methods of post processing are compared. These are the Wavelet Transform as a benchmark method in the form of the Continuous Wavelet Transform, a noise-assisted variant of the Hilbert-Huang Transform and the Zhao-Atlas-Marks Distribution. Comparative diagrams are presented and the advantages and weaknesses of each method are presented.     

Temperature–salinity modeling for Ruwais coastal area in United Arab Emirates

A three dimensional rectangular grid model is applied to resolve the temperature–salinity dynamics of Ruwais, a segment of the UAE coast which is well known as dense water formation zone. The model employs a heat flux module and a turbulence closure scheme that facilitate realistic calculation of temperature–salinity dynamics. A field survey campaign is carried out to support the modeling study, involving measurements of tide, currents, temperature, and salinity. Investigation is done for two meteorologically extreme conditions, i.e. summer and winter. The model study showed that the western flux develops an anticlockwise circulation in the study area. The water industrial discharges elevated the temperature and salinity of the water near the southeastern shoreline. This water mass propagated towards north under the influence of gravity.     

Sampling soil water along the pF curve for δ2H and δ18O analysis

Soil water stable isotopes are widely used across disciplines (e.g., hydrology, ecology, soil science, and biogeochemistry). However, the full potential of stables isotopes as a tool for characterizing the origin, flow path, transport processes and residence times of water in different eco-, hydro-, and geological compartments has not yet been exploited. This is mainly due to the large variety of different methods for pore water extraction. While recent work has shown that matric potential affects the equilibrium fractionation, little work has examined how different water retention characteristics might affect the sampled water isotopic composition. Here, we present a simple laboratory experiment with two well-studied standard soils differing in their physico-chemical properties (e.g., clayey loam and silty sand). Samples were sieved, oven-dried and spiked with water of known isotopic composition to full saturation. For investigating the effect of water retention characteristics on the extracted water isotopic composition, we used pressure extractors to sample isotopically labelled soil water along the pF curve. After pressure extraction, we further extracted the soil samples via cryogenic vacuum extraction. The null hypothesis guiding our work was that water held at different tensions shows the same isotopic composition. Our results showed that the sampled soil water differed isotopically from the introduced isotopic label over time and sequentially along the pF curve. Our and previous studies suggest caution in interpreting isotope results of extracted soil water and a need to better characterize processes that govern isotope fractionation with respect to soil water retention characteristics. In the future, knowledge about soil water retention characteristics with respect to soil water isotopic composition could be applied to predict soil water fractionation effects under natural and non-stationary conditions. In this regard, isotope retention characteristics as an analog to water retention characteristics have been proposed as a way forward since matric potential affects the equilibrium fractionation between the bound water and the water vapour.