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.     

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.     

Approximate analysis methods for asymmetric plan base‐isolated buildings

An approximate method for linear analysis of asymmetric‐plan, multistorey buildings is specialized for a single‐storey, base‐isolated structure. To find the mode shapes of the torsionally coupled system, the Rayleigh–Ritz procedure is applied using the torsionally uncoupled modes as Ritz vectors. This approach reduces to analysis of two single‐storey systems, each with vibration properties and eccentricities (labelled ‘effective eccentricities’) similar to corresponding properties of the isolation system or the fixed‐base structure. With certain assumptions, the vibration properties of the coupled system can be expressed explicitly in terms of these single‐storey system properties. Three different methods are developed: the first is a direct application of the Rayleigh–Ritz procedure; the second and third use simplifications for the effective eccentricities, assuming a relatively stiff superstructure. The accuracy of these proposed methods and the rigid structure method in determining responses are assessed for a range of system parameters including eccentricity and structure flexibility. For a subset of systems with equal isolation and structural eccentricities, two of the methods are exact and the third is sufficiently accurate; all three are preferred to the rigid structure method. For systems with zero isolation eccentricity, however, all approximate methods considered are inconsistent and should be applied with caution, only to systems with small structural eccentricities or stiff structures. Copyright © 2001 John Wiley & Sons, Ltd.     

Asymmetric one‐storey elastic systems with non‐linear viscous and viscoelastic dampers: Simplified analysis and supplemental damping system design

Investigated is the accuracy in estimating the response of asymmetric one‐storey systems with non‐linear viscoelastic (VE) dampers by analysing the corresponding linear viscous system wherein all non‐linear VE dampers are replaced by their energy‐equivalent linear viscous dampers. The response of the corresponding linear viscous system is determined by response history analysis (RHA) and by response spectrum analysis (RSA) extended for non‐classically damped systems. The flexible and stiff edge deformations and plan rotation of the corresponding linear viscous system determined by the extended RSA procedure is shown to be sufficiently accurate for design applications with errors generally between 10 and 20%. Although similar accuracy is also shown for the ‘pseudo‐velocity’ of non‐linear VE dampers, the peak force of the non‐linear VE damper cannot be estimated directly from the peak damper force of the corresponding linear viscous system. A simple correction for damper force is proposed and shown to be accurate (with errors not exceeding 15%). For practical applications, an iterative linear analysis procedure is developed for determining the amplitude‐ and frequency‐dependent supplemental damping properties of the corresponding linear viscous system and for estimating the response of asymmetric one‐storey systems with non‐linear VE dampers from the earthquake design (or response) spectrum. Finally, a procedure is developed for designing non‐linear supplemental damping systems that satisfy given design criteria for a given design spectrum. Copyright © 2003 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.     

A modal pushover analysis procedure to estimate seismic demands for unsymmetric‐plan buildings

An Erratum has been published for this article in Earthquake Engng. Struct. Dyn. 2004; 33:1429. Based on structural dynamics theory, the modal pushover analysis (MPA) procedure retains the conceptual simplicity of current procedures with invariant force distribution, now common in structural engineering practice. The MPA procedure for estimating seismic demands is extended to unsymmetric‐plan buildings. In the MPA procedure, the seismic demand due to individual terms in the modal expansion of the effective earthquake forces is determined by non‐linear static analysis using the inertia force distribution for each mode, which for unsymmetric buildings includes two lateral forces and torque at each floor level. These ‘modal’ demands due to the first few terms of the modal expansion are then combined by the CQC rule to obtain an estimate of the total seismic demand for inelastic systems. When applied to elastic systems, the MPA procedure is equivalent to standard response spectrum analysis (RSA). The MPA estimates of seismic demand for torsionally‐stiff and torsionally‐flexible unsymmetric systems are shown to be similarly accurate as they are for the symmetric building; however, the results deteriorate for a torsionally‐similarly‐stiff unsymmetric‐plan system and the ground motion considered because (a) elastic modes are strongly coupled, and (b) roof displacement is underestimated by the CQC modal combination rule (which would also limit accuracy of RSA for linearly elastic systems). Copyright © 2004 John Wiley & Sons, Ltd.