A gradient‐ratio method for the semi‐automatic interpretation of gravity map data sets

The application of semi‐automatic interpretation techniques to potential field data can be of significant assistance to a geophysicist. This paper generalizes the magnetic vertical contact model tilt‐depth method to gravity data using a vertical cylinder and buried sphere models. The method computes the ratio of the vertical to the total horizontal derivative of data and then identifies circular contours within it. Given the radius of the contour and the contour value itself, the depth to the source can be determined. The method is applied both to synthetic and gravity data from South Africa. The Matlab source code can be obtained from the author upon request.     

Frequency‐domain acoustic‐elastic coupled waveform inversion using the Gauss‐Newton conjugate gradient method

We developed a frequency‐domain acoustic‐elastic coupled waveform inversion based on the Gauss‐Newton conjugate gradient method. Despite the use of a high‐performance computer system and a state‐of‐the‐art parallel computation algorithm, it remained computationally prohibitive to calculate the approximate Hessian explicitly for a large‐scale inverse problem. Therefore, we adopted the conjugate gradient least‐squares algorithm, which is frequently used for geophysical inverse problems, to implement the Gauss‐Newton method so that the approximate Hessian is calculated implicitly. Thus, there was no need to store the Hessian matrix. By simultaneously back‐propagating multi‐components consisting of the pressure and displacements, we could efficiently extract information on the subsurface structures. To verify our algorithm, we applied it to synthetic data sets generated from the Marmousi‐2 model and the modified SEG/EAGE salt model. We also extended our algorithm to the ocean‐bottom cable environment and verified it using ocean‐bottom cable data generated from the Marmousi‐2 model. With the assumption of a hard seafloor, we recovered both the P‐wave velocity of complicated subsurface structures as well as the S‐wave velocity. Although the inversion of the S‐wave velocity is not feasible for the high Poisson's ratios used to simulate a soft seafloor, several strategies exist to treat this problem. Our example using multi‐component data showed some promise in mitigating the soft seafloor effect. However, this issue still remains open.     

Dynamic response and optimal design of structures with large mass ratio TMD

This paper investigates the dynamic behavior and the seismic effectiveness of a non‐conventional Tuned Mass Damper (TMD) with large mass ratio. Compared with conventional TMD, the device mass is increased up to be comparable with the mass of the structure to be protected, aiming at a better control performance. In order to avoid the introduction of an excessive additional weight, masses already present on the structure are converted into tuned masses, retaining structural or architectural functions beyond the mere control function. A reduced order model is introduced for design purposes and the optimal design of a large mass ratio TMD for seismic applications is then formulated. The design method is specifically developed to implement High‐Damping Rubber Bearings (HDRB) to connect the device mass to the main structure, taking advantage of combining stiffness and noticeable damping characteristics. Ground acceleration is modeled as a Gaussian random process with white noise power spectral density. A numerical searching technique is used to obtain the optimal design parameter, the frequency ratio alpha, which minimizes the root‐mean‐square displacement response of the main structure. The study finally comprises shaking table tests on a 1:5 scale model under a wide selection of accelerograms, both artificial and natural, to assess the seismic effectiveness of the proposed large mass ratio TMD. Copyright © 2011 John Wiley & Sons, Ltd.     

Analysis of implicit HHT‐α integration algorithm for real‐time hybrid simulation

Real‐time hybrid simulation is a viable experiment technique to evaluate the performance of structures equipped with rate‐dependent seismic devices when subject to dynamic loading. The integration algorithm used to solve the equations of motion has to be stable and accurate to achieve a successful real‐time hybrid simulation. The implicit HHT α‐algorithm is a popular integration algorithm for conducting structural dynamic time history analysis because of its desirable properties of unconditional stability for linear elastic structures and controllable numerical damping for high frequencies. The implicit form of the algorithm, however, requires iterations for nonlinear structures, which is undesirable for real‐time hybrid simulation. Consequently, the HHT α‐algorithm has been implemented for real‐time hybrid simulation using a fixed number of substep iterations. The resulting HHT α‐algorithm with a fixed number of substep iterations is believed to be unconditionally stable for linear elastic structures, but research on its stability and accuracy for nonlinear structures is quite limited. In this paper, a discrete transfer function approach is utilized to analyze the HHT α‐algorithm with a fixed number of substep iterations. The algorithm is shown to be unconditionally stable for linear elastic structures, but only conditionally stable for nonlinear softening or hardening structures. The equivalent damping of the algorithm is shown to be almost the same as that of the original HHT α‐algorithm, while the period elongation varies depending on the structural nonlinearity and the size of the integration time‐step. A modified form of the algorithm is proposed to improve its stability for use in nonlinear structures. The stability of the modified algorithm is demonstrated to be enhanced and have an accuracy that is comparable to that of the existing HHT α‐algorithm with a fixed number of substep iterations. Both numerical and real‐time hybrid simulations are conducted to verify the modified algorithm. The experimental results demonstrate the effectiveness of the modified algorithm for real‐time testing. Copyright © 2011 John Wiley & Sons, Ltd.     

Enrichment of seismic ground motion databases using Karhunen–Loève expansion

Best estimate seismic analysis are generally based on time‐domain simulations of structural responses. The seismic load is then modeled by a stochastic process representing ground motion. For this purpose, the analyst can use recorded accelerograms or work with synthetically generated ones. The number of ground motion time‐histories available for a given scenario and site condition is limited and generally not sufficient for carrying out more advanced probabilistic structural response analysis. It is then necessary to have at our disposal methods that allow for generating synthetic accelerograms that realistically characterize earthquake ground motions. However, most of the methods proposed in literature for generating synthetic accelerograms do not accurately reproduce the natural variability of ground motion parameters (such as PGA, cumulative absolute velocity, and Arias intensity) observed for recorded time histories. In this paper, we introduce a new method for generating synthetic ground motion, based on Karhunen‐Loève decomposition and a non‐Gaussian stochastic model. The proposed method enables the structural analyst to simulate ground motion time histories featuring the properties mentioned above. To demonstrate its capability, we study the influence of the simulation method on different ground motion parameters and on soil response spectra. We finally compute fragility curves to illustrate the practical application of the proposed method. Copyright © 2012 John Wiley & Sons, Ltd.     

Transverse free vibrations of continuous bridges with abutment restraint

This paper analyzes the dynamic behaviour of continuous bridges with transverse abutment restraint. These are a particular class of multi‐span bridges whose transverse motion is restrained at the abutments. The transverse dynamic behaviour of these partially‐restrained bridges is described by a model consisting of a two‐dimensional simply‐supported beam with intermediate visco‐elastic restraints. A minimal set of characteristic problem parameters that completely describe the dynamic response is initially identified and an analytical solution is developed for particular configurations. The influence of these characteristic parameters on the properties of the dynamic systems are analyzed by considering the undamped and the damped free vibrations separately. Successively, in order to cover more general and complex configurations of realistic bridges, an extension of the definition of the characteristic parameters is proposed. The results obtained considering selected case studies are finally presented. Copyright © 2011 John Wiley & Sons, Ltd.     

Dynamic interaction of bridge–train system under non‐uniform seismic ground motion

The probability that an earthquake occurs when a train is running over a bridge in earthquake‐prone regions is much higher than before, for high‐speed railway lines are rapidly developed to connect major cities worldwide. This paper presents a finite element method‐based framework for dynamic analysis of coupled bridge–train systems under non‐uniform seismic ground motion, in which rail–wheel interactions and possible separations between wheels and rails are taken into consideration. The governing equations of motion of the coupled bridge–train system are established in an absolute coordinate system. Without considering the decomposition of seismic responses into pseudo‐static and inertia‐dynamic components, the equations of motion of the coupled system are formed in terms of displacement seismic ground motions. The mode superposition method is applied to the bridge structure to make the problem manageable while the Newmark‐β method with an iterative computation scheme is used to find the best solution for the problem concerned. Eight high‐speed trains running over a multi‐span steel truss‐arch bridge subject to earthquakes are taken as a case study. The results from the case study demonstrate that the spatial variation of seismic ground motion affects dynamic responses of the bridge–train system. The ignorance of pseudo‐static component when using acceleration seismic ground motions as input may underestimate seismic responses of the bridge–train system. The probability of separation between wheels and rails becomes higher with increasing train speed. Copyright © 2011 John Wiley & Sons, Ltd.     

Fragility analysis of low‐rise masonry in‐filled reinforced concrete buildings by a coefficient‐based spectral acceleration method

This study presents a seismic fragility analysis of low‐rise masonry in‐filled (MI) reinforced concrete (RC) buildings using a proposed coefficient‐based spectral acceleration method. The coefficient‐based method, without requiring any complicated finite element analysis, is a simplified procedure for assessing the spectral acceleration demand (or capacity) of buildings subjected to earthquakes. This paper begins with a calibration of the proposed coefficient‐based method for low‐rise MI RC buildings using published experimental results obtained from shaking table tests. Spectral acceleration‐based fragility curves for low‐rise MI RC buildings under various inter‐story drift limits are then constructed using the calibrated coefficient‐based method. A comparison of the experimental and estimated results indicates that the simplified coefficient‐based method can provide good approximations of the spectral accelerations at peak loads of low‐rise MI RC buildings, if a proper set of drift‐related factors and initial fundamental periods of structures are used. Moreover, the fragility curves constructed using the coefficient‐based method can provide a satisfactory vulnerability evaluation for low‐rise MI RC buildings under a given performance level. Copyright © 2011 John Wiley & Sons, Ltd.     

Inelastic displacement demands in steel structures and their relationship with earthquake frequency content parameters

This paper deals with the estimation of peak inelastic displacements of SDOF systems, representative of typical steel structures, under constant relative strength scenarios. Mean inelastic deformation demands on bilinear systems (simulating moment resisting frames) are considered as the basis for comparative purposes. Additional SDOF models representing partially‐restrained and concentrically‐braced (CB) frames are introduced and employed to assess the influence of different force‐displacement relationships on peak inelastic displacement ratios. The studies presented in this paper illustrate that the ratio between the overall yield strength and the strength during pinching intervals is the main factor governing the inelastic deformations of partially‐restrained models and leading to significant differences when compared with predictions based on bilinear structures, especially in the short‐period range. It is also shown that the response of CB systems can differ significantly from other pinching models when subjected to low or moderate levels of seismic demand, highlighting the necessity of employing dedicated models for studying the response of CB structures. Particular attention is also given to the influence of a number of scalar parameters that characterise the frequency content of the ground motion on the estimated peak displacement ratios. The relative merits of using the average spectral period T_aver, mean period T_m, predominant period T_g, characteristic period T_c and smoothed spectral predominant period T_o of the earthquake ground motion, are assessed. This paper demonstrates that the predominant period, defined as the period at which the input energy is maximum throughout the period range, is the most suitable frequency content scalar parameter for reducing the variability in displacement estimations. Finally, noniterative equivalent linearisation expressions based on the secant period and equivalent damping ratios are presented and verified for the prediction of peak deformation demands in steel structures. Copyright © 2011 John Wiley & Sons, Ltd.