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.     

Electron impact polarization of resonance lines from Lithium-like ions

The Oppenheimer-Penney theory, as developed by Percival and Seaton (1958), is applied to calculate the polarization of resonance lines from Li-like ions. Two laws for the pitch-angle distribution of electrons around the magnetic field are accounted. The degrees of polarization are averaged over the energy of non-thermal electrons generated during the initial phase of solar flares. It is found that for the full space pitch-angle distribution, as adopted by Chandra and Joshi (1984), the degrees of polarization are nearly independent of the atomic number of ion. Whereas for the forward-come distribution used by Haug (1981), they depend on the choice of the free parameterE ₀. The polarization of the resonance lines from Li-like ions is two times larger than that of the L radiations from H-like ions. Hence, under favourable conditions, it may be detected during solar flares.     

Characteristics of spectral aerosol optical depths over India during ICARB

Spectral aerosol optical depth (AOD) measurements, carried out regularly from a network of observatories spread over the Indian mainland and adjoining islands in the Bay of Bengal and Arabian Sea, are used to examine the spatio-temporal and spectral variations during the period of ICARB (March to May 2006). The AODs and the derived Ångström parameters showed considerable variations across India during the above period. While at the southern peninsular stations the AODs decreased towards May after a peak in April, in the north Indian regions they increased continuously from March to May. The Ångström coefficients suggested enhanced coarse mode loading in the north Indian regions, compared to southern India. Nevertheless, as months progressed from March to May, the dominance of coarse mode aerosols increased in the columnar aerosol size spectrum over the entire Indian mainland, maintaining the regional distinctiveness. Compared to the above, the island stations showed considerably low AODs, so too the northeastern station Dibrugarh, indicating the prevalence of cleaner environment. Long-range transport of aerosols from tshe adjoining regions leads to remarkable changes in the magnitude of the AODs and their wavelength dependencies during March to May. HYSPLIT back-trajectory analysis shows that enhanced long-range transport of aerosols, particularly from the west Asia and northwest coastal India, contributed significantly to the enhancement of AOD and in the flattening of the spectra over entire regions; if it is the peninsular regions and the island Minicoy are more impacted in April, the north Indian regions including the Indo Gangetic Plain get affected the most during May, with the AODs soaring as high as 1.0 at 500 nm. Over the islands, the Ångström exponent (α) remained significantly lower (～1) over the Arabian Sea compared to Bay of Bengal (BoB) (～1.4) as revealed by the data respectively from Minicoy and Port Blair. Occurrences of higher values of α, showing dominance of accumulation mode aerosols, over BoB are associated well with the advection, above the boundary layer, of fine particles from the east Asian region during March and April. The change in the airmass to marine in May results in a rapid decrease in α over the BoB.     

Urban effects of Chennai on sea breeze induced convection and precipitation

Doppler radar derived wind speed and direction profiles showed a well developed sea breeze circulation over the Chennai, India region on 28 June, 2003. Rainfall totals in excess of 100 mm resulted from convection along the sea breeze front. Inland propagation of the sea breeze front was observed in radar reflectivity imagery. High-resolution MM5 simulations were used to investigate the influence of Chennai urban land use on sea breeze initiated convection and precipitation. A comparison of observed and simulated 10m wind speed and direction over Chennai showed that the model was able to simulate the timing and strength of the sea breeze. Urban effects are shown to increase the near surface air temperature over Chennai by 3.0K during the early morning hours. The larger surface temperature gradient along the coast due to urban effects increased onshore flow by 4.0m s⁻¹. Model sensitivity study revealed that precipitation totals were enhanced by 25mm over a large region 150 km west of Chennai due to urban effects. Deficiency in model physics related to night-time forecasts are addressed.     

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.