Topography effects in the 1999 Athens earthquake : engineering issues in seismology

Thesis (Sc. D.)--Massachusetts Institute of Technology, Dept. of Civil and Environmental Engineering, 2004.Includes bibliographical references.It is well known that irregular topography can substantially affect the amplitude and frequency characteristics of seismic motion. Macroseismic observations of destructive earthquakes often show higher damage intensity at the tops of hills, ridges and canyons than at lower elevations and on flat areas. Systematic seismic motion amplification over convex topographies has been confirmed by instrumental studies and also predicted by theoretical and numerical simulations of wave diffraction. Nonetheless, for the most part, the former have been limited to weak motion data and the later have treated topographic asperities as simple geometric irregularities on the surface of homogeneous, linearly elastic halfspaces. Despite the qualitative agreement between theory and observations on topography effects, there is still much uncertainty concerning the actual severity of amplification near topographic irregularities, inasmuch as predictive methods are still lacking on the quantitative aspects of seismic amplification near such features. Focusing of seismic rays by convex topographies does play a significant role as shown theoretically, yet it is not the only physical phenomenon involved. On the other hand, weak motion data may not be applicable to describe topography effects for strong shaking, and indeed there exist very few- if any- well documented case studies demonstrating the severity of topographic effects for strong ground motion. In this dissertation, we find that topography and local soil conditions need to be accounted for simultaneously for the prediction of site amplification factors, especially when earthquake motions are strong enough to elicit clear nonlinear soil behavior.(cont.) We examine how local stratigraphy, material heterogeneity and nonlinear soil response can alter the focusing mechanism at the vertex of cliff-type topographies, and how the free-field response is further modified on account of soil-structure interaction. By means of a case-study from the Athens 1999 earthquake, we validate the effects of local soil conditions by comparison with weak motion data, and illustrate the effects of nonlinear soil behavior and soil-structure interaction on strong motion amplification. Our finite-element, nonlinear simulations seem to explain the uneven distribution of severe damage in the community of Adàmes that borders the crest of the Kifissos river canyon at its deepest point. They also resolve in part previously unexplained discrepancies, often observed between strong amplification during actual earthquakes and moderate values predicted by simple theoretical models. Combining our findings with earlier published results, we propose a period- and space-dependent factor, referred to as Topographic Aggravation Factor (TAF), which can be used in engineering design to modify site-specific design spectra of seismic code provisions to account for topography effects.by Dominic Assimaki.Sc.D

User's guide to computer programs JET 5A and CIVM-JET 5B to calculate the large elastic-plastic dynamically-induced deformations of multilayer partial and/or complete structural rings

These structural ring deflections lie essentially in one plane and, hence, are called two-dimensional (2-d). The structural rings may be complete or partial; the former may be regarded as representing a fragment containment ring while the latter may be viewed as a 2-d fragment-deflector structure. These two types of rings may be either free or supported in various ways (pinned-fixed, locally clamped, elastic-foundation supported, mounting-bracket supported, etc.). The initial geometry of each ring may be circular or arbitrarily curved; uniform-thickness or variable-thickness rings may be analyzed. Strain-hardening and strain-rate effects of initially-isotropic material are taken into account. An approximate analysis utilizing kinetic energy and momentum conservation relations is used to predict the after-impact velocities of each fragment and of the impact-affected region of the ring; this procedure is termed the collision-imparted velocity method (CIVM) and is used in the CIVM-JET 5 B program. This imparted-velocity information is used in conjunction with a finite-element structural response computation code to predict the transient, large-deflection, elastic-plastic responses of the ring. Similarly, the equations of motion of each fragment are solved in small steps in time. Provisions are made in the CIVM-JET 5B code to analyze structural ring response to impact attack by from 1 to 3 fragments, each with its own size, mass, translational velocity components, and rotational velocity. The effects of friction between each fragment and the impacted ring are included