Finite-element analysis of time-dependent large-deformation problems

An updated Lagrangian finite-element formulation has been developed for time-dependent problems of soil consolidation involving finite deformations. Large plastic strains as well as rotations occur in such problems and nominal stress measures are introduced in the formulation to redefine stresses. This leads to corrective terms for equilibrium and yield violations in addition to geometric stiffening terms in the governing integral equations. The soil is considered to be either a linear elastic or an elastoplastic, critical-state material. Some simple numerical examples are studied to validate the formulation, followed by a detailed analysis of the problem of penetration of a pile into soil. The results of this problem are viewed with emphasis on the physical interpretation and practical significance.     

Intensity map of Mw 6.9 2011 Sikkim–Nepal border earthquake and its relationships with PGA: distance and magnitude

We compiled available information of damages and other effects caused by the September 18, 2011, Sikkim–Nepal border earthquake from the print and electronic media, and interpreted them to obtain Modified Mercalli Intensity (MMI) at over 142 locations. These values are used to prepare the intensity map of the Sikkim earthquake. The map reveals several interesting features. Within the meizoseismal area, the most heavily damaged villages are concentrated toward the eastern edge of the inferred fault, consistent with eastern directivity. The intensities are amplified significantly in areas located along rivers, within deltas or on coastal alluvium such as mud flats and salt pans. We have also derived empirical relation between MMI and ground motion parameters using least square regression technique and compared it with the available relationships available for other regions of the world. Further, seismic intensity information available for historical earthquakes which have occurred in NE Himalayas along with present intensity has been utilized for developing attenuation relationship for NE India using two-step regression analyses. The derived attenuation relation is useful for assessing damage of a potential future earthquake (earthquake scenario-based planning purposes) for the northeast Himalaya region.

     

Three‐dimensional modal pushover analysis of buildings subjected to two components of ground motion,including its evaluation for tall buildings

The modal pushover analysis (MPA) procedure, presently restricted to one horizontal component of ground motion, is extended to three‐dimensional analysis of buildings—symmetric or unsymmetric in plan—subjected to two horizontal components of ground motion, simultaneously. Also presented is a variant of this method, called the practical modal pushover analysis (PMPA) procedure, which estimates seismic demands directly from the earthquake response (or design) spectrum. Its accuracy in estimating seismic demands for very tall buildings is evaluated, demonstrating that for nonlinear systems this procedure is almost as accurate as the response spectrum analysis procedure is for linear systems. Thus, for practical applications, the PMPA procedure offers an attractive alternative whereby seismic demands can be estimated directly from the (elastic) design spectrum, thus avoiding the complications of selecting and scaling ground motions for nonlinear response history analysis. Copyright © 2010 John Wiley & Sons, Ltd.     

A framework for the evaluation of ground motion selection and modification procedures

This study develops a framework to evaluate ground motion selection and modification (GMSM) procedures. The context is probabilistic seismic demand analysis, where response history analyses of a given structure, using ground motions determined by a GMSM procedure, are performed in order to estimate the seismic demand hazard curve (SDHC) for the structure at a given site. Currently, a GMSM procedure is evaluated in this context by comparing several resulting estimates of the SDHC, each derived from a different definition of the conditioning intensity measure (IM). Using a simple case study, we demonstrate that conclusions from such an approach are not always definitive; therefore, an alternative approach is desirable. In the alternative proposed herein, all estimates of the SDHC from GMSM procedures are compared against a benchmark SDHC, under a common set of ground motion information. This benchmark SDHC is determined by incorporating a prediction model for the seismic demand into the probabilistic seismic hazard analysis calculations. To develop an understanding of why one GMSM procedure may provide more accurate estimates of the SDHC than another procedure, we identify the role of ‘IM sufficiency’ in the relationship between (i) bias in the SDHC estimate and (ii) ‘hazard consistency’ of the corresponding ground motions obtained from a GMSM procedure. Finally, we provide examples of how misleading conclusions may potentially be obtained from erroneous implementations of the proposed framework. Copyright © 2014 John Wiley & Sons, Ltd.     

Evaluation of the exact conditional spectrum and generalized conditional intensity measure methods for ground motion selection

Two existing, contemporary ground motion selection and modification procedures – (i) exact conditional spectrum (CS‐exact) and (ii) generalized conditional intensity measure (GCIM) – are evaluated in their ability to accurately estimate seismic demand hazard curves (SDHCs) of a given structure at a specified site. The amount of effort involved in implementing these procedures to compute a single SDHC is studied, and a case study is chosen where rigorous benchmark SDHCs can be determined for evaluation purposes. By comparing estimates from ground motion selection and modification procedures with the benchmark, we conclude that estimates from CS‐exact are unbiased in many of the cases considered. The estimates from GCIM are even more accurate, as they are unbiased for most – but not all – of the cases where estimates from CS‐exact are biased. We find that it is possible to obtain biased SDHCs from GCIM, even after employing a very diverse collection of intensity measures to select ground motions and implementing its bias‐checking feature, because it is usually difficult to identify intensity measures that are truly ‘sufficient’ for the response of a complex, multi‐degree‐of‐freedom system. Copyright © 2015 John Wiley & Sons, Ltd.     

Estimating the seismic displacement of friction pendulum isolators based on non‐linear response history analysis

A procedure for developing equations that estimate the isolator displacement due to strong ground motion is applied to buildings isolated with the friction pendulum system. The resulting design equations, based on rigorous non‐linear analysis, offer an alternative to the iterative equivalent‐linear methods used by current U.S. building codes. The governing equations of the system are reduced to a form such that the median normalized displacement of the system due to an ensemble of ground motions is found to depend on only the isolation period—a function of the curvature of the isolator—and the friction force at incipient slip normalized by peak ground velocity. The normalization is effective in minimizing the dispersion of the normalized displacement for an ensemble of ground motions, implying that the median normalized displacement is a reliable estimate of response. The design equations reflect the significant (20 to 38%) increase in displacement when the excitation includes two lateral components of ground motion instead of just one component. Equivalent‐linear methods are shown to underestimate by up to 30% the exact median displacement determined by non‐linear response history analysis for one component of ground motion, and building codes include at most a 4.4% increase for a second component. Copyright © 2003 John Wiley & Sons, Ltd.     

Three-dimensional analysis of spatially varying ground motions around a uniform canyon in a homogeneous half-space

Department of Civil Engineering, University of California, Berkeley, CA 94720, U.S.A. A direct boundary element method to determine the three-dimensional seismic response of an infinitely-long canyon of arbitrary but uniform cross-section cut in a homogeneous viscoelastic half-space is presented. The seismic excitation is represented by P, SV, SH or Rayleigh waves at arbitrary angles with respect to the axis of the canyon. The accuracy of the procedure and implementing computer program is demonstrated by comparison with previous solutions for the limiting case of two-dimensional response, recently obtained three-dimensional response results for infinitely-long canyons, and three-dimensional boundary method solutions presented in this paper for finite canyons.     

Numerical evaluation of the damping‐solvent extraction method in the frequency domain

The damping‐solvent extraction method for the analysis of unbounded visco‐elastic media is evaluated numerically in the frequency domain in order to investigate the influence of the computational parameters—domain size, amount of artificial damping, and mesh density—on the accuracy of results. An analytical estimate of this influence is presented, and specific questions regarding the influence of the parameters on the results are answered using the analytical estimate and numerical results for two classical problems: the rigid strip and rigid disc footings on a visco‐elastic half‐space with constant hysteretic material damping. As the domain size is increased, the results become more accurate only at lower frequencies, but are essentially unaffected at higher frequencies. Choosing the domain size to ensure that the static stiffness is computed accurately leads to an unnecessarily large domain for analysis at higher frequencies. The results improve by increasing artificial damping but at a slower rate as the total (material plus artificial) damping ratio ζ_t gets closer to 0.866. However, the results do not deteriorate significantly for the larger amounts of artificial damping, suggesting that ζ_t≈0.6 is appropriate; a larger value is not likely to influence the accuracy of results. Presented results do not support the earlier suggestion that similar accuracy can be achieved by a large bounded domain with small damping or by a small domain with larger damping. Copyright © 2002 John Wiley & Sons, Ltd.