Dynamic analysis by direct superposition of Ritz vectors

The solution of the eigenvalue problem for large structures is often the most costly phase of a dynamic response analysis. In this paper, the need for the exact solution of this large eigenvalue problem is eliminated. A new algorithm, based on error minimization, is presented for the generation of a sequence of Ritz vectors. These orthogonal vectors are used to reduce the size of the system. Only Ritz vectors with a large participation factor are used in the subsequent mode superposition analysis. In all examples studied, the superposition of Ritz vectors yields more accurate results, with fewer vectors, than if the exact eigenvectors are used. The proposed method not only reduces computer time requirements significantly but provides an error estimation for the dynamic analysis. The approach automatically includes the advantages of the proven numerical techniques of static condensation, Guyan reduction and static correction due to higher mode truncation.     

Solution of viscously damped linear systems using a set of load-dependent vectors

This paper considers a solution method for viscously damped linear structural systems which are subjected to transient loading. The equations of motion of such systems are written in a first-order form. A solution subspace is generated using the damped dynamic matrix and the static deflection from the first-order form of the equations of motion. Two convenient bases, Lanczos vectors and Ritz vectors, are constructed from this subspace. An approximate solution is then obtained by superposition of the Lanczos vectors or the Ritz vectors. In contrast to the traditional mode superposition method using complex eigenvectors, the Lanczos vectors or the Ritz vectors are less expensive to generate than the complex eigenvectors, yet yield comparable accuracy. In addition, there is no need for a static correction since the static deflection is already contained in our solution subspace. Numerical examples are presented to show the potential of using the Ritz vectors to compute responses of damped dynamic systems.     

Frequency dependent ritz vectors

In order to reduce the size of problems involving analysis of the dynamic response of structural systems, a transformatio based on appropriately selected Ritz shapes is commonly employed. The lower mode shapes may at times serve a effective Ritz shapes. However, the computation of mode shapes is a time consuming task; in addition, the mode shapes may not form the best basis for representing the spatial distribution of loads. The recently developed load dependent vectors, which are derived from a static solution for the applied loads, address some of the problems inherent in the use of mode shapes. However, both the natural mode shapes and the load dependent vectors fail to account for the frequency content of the loading, a parameter that may influence strongly the response, particularly for loading with a high frequency content. A procedure is presented here for the generation of frequency dependent vectors. A combination of load dependent and frequency dependent vectors will often form a very efficient basis for the representation of the response, as illustrated by several examples presented here.     

Method for solving dynamic equilibrium equations based on displacement excitation

Mode superposition is a widely used method for solving the dynamic equilibrium equation in structural dynamic analysis. However, the accuracy of this method may be reduced when the dynamic equilibrium equations are set up using displacement excitation. A new method for developing solutions for dynamic equilibrium equations based on displacement excitation is introduced. The dynamic equilibrium equation is decomposed into two parts, namely displacement excitation and velocity excitation, and precise integration and mode superposition methods are combined to solve the equation. Ritz vectors are then used to calculate the static response of the truncated modes of the structure, and a method for determining the number of participating modes is obtained. Using multi-degree-of-freedom systems as two computational examples, the differences in the structural responses obtained from the displacement excitation and acceleration excitation are compared and analyzed. It is shown that the new solution method generates consistent accuracy between the displacement excitation and acceleration excitation.     

Multiple‐support response spectrum analysis using load‐dependent Ritz vectors

Load‐dependent Ritz (LDR) vectors are used in conjunction with the multiple‐support response spectrum (MSRS) combination rule for analysis of structures subjected to spatially varying earthquake ground motions. The LDR vector approach for MSRS analysis is motivated by the fact that LDR vectors in general are more accurate and computationally simpler than eigenvectors in mode superposition analysis, and because many researchers and engineers are using LDR vectors in linear structural dynamic analysis. Mode truncation rules for the original MSRS method are modified to apply to LDR vectors. Two methods for selecting LDR vectors for multicomponent MSRS analysis are introduced. Idealized models of two real bridges with differing structural characteristics are used to investigate the accuracy and efficiency of the two LDR‐MSRS methods in comparison with results obtained by the original MSRS method as well as an extended version that accounts for the static contribution of truncated modes. The results show that the LDR‐MSRS method is generally more accurate than the original MSRS method and at least as accurate as the extended MSRS method. Copyright © 2014 John Wiley & Sons, Ltd.     

On the accuracy of mode superposition analysis in structural dynamics

It is pointed out that the number of modes which should be included in a mode superposition dynamic response analysis depends on both the frequency content and the distribution of the loading. If the loading frequency is low the effect of the higher modes can be approximated by a static analysis. A technique is described for calculating this static contribution from the higher modes; the total response is then represented by the sum of the lower mode dynamic response and the higher mode static effects. The effectiveness of the procedure is demonstrated by a numerical example.     

Ritz method for dynamic analysis of large discrete linear systems with non-proportional damping

Real and complex Ritz vector bases for dynamic analysis of large linear systems with non-proportional damping are presented and compared. Both vector bases are generated utilizing load dependent vector algorithms that employ recurrence equations analogous to the Lanczos algorithm. The choice of static response to fixed spatial loading distribution, as a starting vector in recurrence equations, is motivated by the static correction concept. Different phases of dynamic response analysis are compared with respect to computational efficiency and accuracy. It is concluded that the real vector basis approach is approximately eight times more efficient than the complex vector basis approach. The complex vector basis has some advantages with respect to accuracy, if the excitation is of piecewise linear form, since the exact solution can be utilized. In addition, it is demonstrated that both Ritz vector bases, real and complex, possess superior accuracy over the adequate eigenvector bases.     

Earthquake analysis of concrete gravity dams including dam-water-foundation rock interaction

A general procedure for analysis of the response of concrete gravity dams, including the dynamic effects of impounded water and flexible foundation rock, to the transverse (horizontal) and vertical components of earthquake ground motion is presented. The problem is reduced to one in two dimensions, considering the transverse vibration of a monolith of the dam. The system is analysed under the assumption of linear behaviour for the concrete, foundation rock and water. The complete system is considered as composed of three substructures—the dam, represented as a finite element system, the fluid domain, as a continuum of infinite length in the upstream direction, and the foundation rock region as a viscoelastic half-plane. The structural displacements of the dam are expressed as a linear combination of Ritz vectors, chosen as normal modes of an associated undamped dam-rock system. The effectiveness of this analytical formulation lies in its being able to produce excellent results by considering only a few Ritz vectors. The generalized displacements due to earthquake motion are computed by synthesizing their complex frequency responses using Fast Fourier Transform procedures. The stress responses are calculated from the displacements. An example analysis is presented to illustrate results obtained from this analytical procedure. Computation times for several analyses are presented to illustrate the effectiveness of the procedure.     

Generalized force vectors for multi‐mode pushover analysis

A generalized pushover analysis (GPA) procedure is developed for estimating the inelastic seismic response of structures under earthquake ground excitations. The procedure comprises applying different generalized force vectors separately to the structure in an incremental form with increasing amplitude until a prescribed seismic demand is attained for each generalized force vector. A generalized force vector is expressed as a combination of modal forces, and simulates the instantaneous force distribution acting on the system when a given response parameter reaches its maximum value during dynamic response to a seismic excitation. While any response parameter can be selected arbitrarily, generalized force vectors in the presented study are derived for maximum interstory drift parameters. The maximum value of any other response parameter is then obtained from the envelope of GPAs results. Each nonlinear static analysis under a generalized force vector activates the entire multi‐degree of freedom effects simultaneously. Accordingly, inelastic actions develop in members with the contribution of all ‘instantaneous modes’ in the nonlinear response range. Target seismic demands for interstory drifts at the selected stories are calculated from the associated drift expressions. The implementation of the proposed GPA is simpler compared with nonlinear response history analysis, whereas it is less demanding in computational effort when compared with several multi‐mode adaptive nonlinear static procedures. Moreover, it does not suffer from the statistical combination of inelastic modal responses obtained separately. The results obtained from building frames have demonstrated that GPA is successful in estimating maximum member deformations and member forces with reference to the response history analysis. When the response is linear elastic, GPA and response spectrum analysis produce identical results. Copyright © 2010 John Wiley & Sons, Ltd.     

Generalized force vectors for multi‐mode pushover analysis of torsionally coupled systems

A generalized multi‐mode pushover analysis procedure was developed for estimating the maximum inelastic seismic response of symmetrical plan structures under earthquake ground excitations. Pushover analyses are conducted with story‐specific generalized force vectors in this procedure, with contributions from all effective modes. Generalized pushover analysis procedure is extended to three‐dimensional torsionally coupled systems in the presented study. Generalized force distributions are expressed as the combination of modal forces to simulate the instantaneous force distribution acting on the system when the interstory drift at a story reaches its maximum value during seismic response. Modal contributions to the generalized force vectors are calculated by a modal scaling rule, which is based on the complete quadratic combination. Generalized forces are applied to the mass centers of each story incrementally for producing nonlinear static response. Maximum response quantities are obtained when the individual frames attain their own target interstory drift values in each story. The developed procedure is tested on an eight‐story frame under 15 ground motions, and assessed by comparing the results obtained from nonlinear time history analysis. The method is successful in predicting the torsionally coupled inelastic response of frames responding to large interstory drift demands. Copyright © 2014 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.     

Use of ritz vectors in wave propagation and foundation response

The accuracy of a numerical method is demonstrated for the dynamic analysis of large complex finite element systems in which the spatial distribution of the loading is constant. The method is based on the use of a special class of Ritz vectors which were previously proposed and can be generated with minimum numerical effort. The purpose of this paper is to extend the use of these vectors to the solution of wave propagation and foundation response problems. The method is applied to one-, twoand three-dimensional problems in order to illustrate the efficiency and accuracy of the technique. Unless it is necessary to evaluate the very high-frequency behaviour of a structural system, it is shown that a small number of Ritz vectors will produce excellent results. Therefore, they can be very effective in the solution of three-dimensional soil-structure systems subjected to earthquake loading.     

Dynamic analysis of flexible rectangular fluid containers subjected to horizontal ground motion

In this paper, an analytical method is proposed to determine the dynamic response of 3‐D rectangular liquid storage tanks with four flexible walls, subjected to horizontal seismic ground motion. Fluid–structure interaction effects on the dynamic responses of partially filled fluid containers, incorporating wall flexibility, are accounted for in evaluating impulsive pressure. The velocity potential in which boundary conditions are satisfied is solved by the method of separation of variables using the principle of superposition. The impulsive pressure distribution is then computed. Solutions based on 3‐D modeling of the rectangular containers are obtained by applying the Rayleigh–Ritz method using the vibration modes of flexible plates with suitable boundary conditions. Trigonometrical functions that satisfy boundary conditions of the storage tank such that the flexibility of the wall is thoroughly considered are used to define the admissible vibration modes. The analysis is then performed in the time domain. Moreover, an analytical procedure is developed for deriving a simple formula that evaluates convective pressure and surface displacements in a similar rigid tank. The variation of dynamic response characteristics with respect to different tank parameters is investigated. A mechanical model, which takes into account the deformability of the tank wall, is developed. The parameters of such a model can be obtained from developed charts, and the maximum seismic loading can be predicted by means of a response spectrum characterizing the design earthquake. Accordingly, a simplified but sufficiently accurate design procedure is developed to improve code formulas for the seismic design of liquid storage tanks. Copyright © 2013 John Wiley & Sons, Ltd.     

Earthquake response analysis of multistorey buildings including foundation interaction

An efficient method, based on the Ritz concept, for dynamic analysis of response of multistorey buildings including foundation interaction to earthquake ground motion is presented. The system considered is a shear building on a rigid circular disc footing attached to the surface of a linearly elastic halfspace. In this method, the structural displacements are transformed to normal modes of vibration of the building on a rigid foundation. The analysis procedure is developed and numerical results are presented to demonstrate that excellent results can be obtained by considering only the first few modes of vibration. As the number of unknowns are reduced by transforming to generalized co-ordinates, the method presented is much more efficient than direct methods.     

Comparison between non-linear dynamic and static seismic analysis of structures according to European and US provisions

Several procedures for non-linear static and dynamic analysis of structures have been developed in recent years. This paper discusses those procedures that have been implemented into the latest European and US seismic provisions: non-linear dynamic time-history analysis; N2 non-linear static method (Eurocode 8); non-linear static procedure NSP (FEMA 356) and improved capacity spectrum method CSM (FEMA 440). The presented methods differ in respect to accuracy, simplicity, transparency and clarity of theoretical background. Non-linear static procedures were developed with the aim of overcoming the insufficiency and limitations of linear methods, whilst at the same time maintaining a relatively simple application. All procedures incorporate performance-based concepts paying more attention to damage control. Application of the presented procedures is illustrated by means of an example of an eight-storey reinforced concrete frame building. The results obtained by non-linear dynamic time-history analysis and non-linear static procedures are compared. It is concluded that these non-linear static procedures are sustainable for application. Additionally, this paper discusses a recommendation in the Eurocode 8/1 that the capacity curve should be determined by pushover analysis for values of the control displacement ranging between zero and 150% of the target displacement. Maximum top displacement of the analyzed structure obtained by using dynamic method with real time-history records corresponds to 145% of the target displacement obtained using the non-linear static N2 procedure.     

Non-linear seismic response analysis using vector superposition methods

Computational algorithms based on the pseudo-force method (PFM) and the tangent spectrum method (TSM) for the seismic analysis of elasto-plastic MDOF structures by mode superposition are presented. The emphasis is put on the effect of the truncation of higher modes on the convergence of the ductility demand and energy dissipated during the earthquake. Eigenvectors and load-dependent vectors have been used in comparative analyses. Applications on a flexible 25-storey building and a stiff, 5-storey shear building indicate that, for a flexible structure, the ductility demand computed from the PFM is more sensitive to basis truncation than that computed from the TSM. However, for the stiff structure, the opposite behaviour is observed. The results indicate that the use of load-dependent vectors in inelastic analyses maintains the computational advantages found for elastic analyses in previous investigations. Although the PFM is more stable and computationally more effective than the TSM, it does not provide any information on the evolution of tangent modal properties in time that reflects dynamic response modification as the structure becomes inelastic.     

Towards a simple dynamic process conceptualization in rainfall–runoff models using multi‐criteria calibration and tracers in temperate,upland catchments

Empirically based understanding of streamflow generation dynamics in a montane headwater catchment formed the basis for the development of simple, low‐parameterized, rainfall–runoff models. This study was based in the Girnock catchment in the Cairngorm Mountains of Scotland, where runoff generation is dominated by overland flow from peaty soils in valley bottom areas that are characterized by dynamic expansion and contraction of saturation zones. A stepwise procedure was used to select the level of model complexity that could be supported by field data. This facilitated the assessment of the way the dynamic process representation improved model performance. Model performance was evaluated using a multi‐criteria calibration procedure which applied a time series of hydrochemical tracers as an additional objective function. Flow simulations comparing a static against the dynamic saturation area model (SAM) substantially improved several evaluation criteria. Multi‐criteria evaluation using ensembles of performance measures provided a much more comprehensive assessment of the model performance than single efficiency statistics, which alone, could be misleading. Simulation of conservative source area tracers (Gran alkalinity) as part of the calibration procedure showed that a simple two‐storage model is the minimum complexity needed to capture the dominant processes governing catchment response. Additionally, calibration was improved by the integration of tracers into the flow model, which constrained model uncertainty and improved the hydrodynamics of simulations in a way that plausibly captured the contribution of different source areas to streamflow. This approach contributes to the quest for low‐parameter models that can achieve process‐based simulation of hydrological response. Copyright © 2009 John Wiley & Sons, Ltd.     

Seismic fracture analysis of concrete gravity dams

A numerical procedure for evaluation of the fracture process of gravity dams during strong earthquakes is presented. The BEM is used to discretize the dam reservoir system including the crack surfaces, and stress intensity factors at the crack tip are employed in a stage by stage procedure which simulates the crack extension. For each stage of constant crack length the mode superposition technique is applied; this is made possible by simulating the impact process of crack closing by a load pulse applied at the contact points which permits the structural stiffness to be assumed unchanged. To verify the proposed procedure, a cantilever beam model structure made of gypsum was tested on a shaking table. Good correlation with the numerical results was obtained, from which it is concluded that the procedure can be employed for evaluation of the crack propagation process in concrete structures subjected to dynamic loadings.     

A fundamental solution due to a periodic point force in the interior of an elastic half-space

The fundamental solution for a periodic point force in the interior of a three-dimensional, homogeneous, isotropic, elastic half-space is derived. The method of synthesis and superposition is employed to obtain the solution in the Laplace transform as well as the frequency domain. These correspond to the dynamic equivalent of Mindlin's static half-space point force solutions. It is reduced, for certain limiting conditions, to the dynamic equivalent of Boussinesq's and Cerruti's problems of a normal and tangential periodic point force respectively, on the boundary of a half-space. Also, static solutions of Mindlin, Boussinesq and Cerruti are recovered for small frequency parameters. Finally, results are presented and compared with other available solutions.     

Stability and accuracy analysis of direct integration methods

A systematic procedure is presented for the stability and accuracy analysis of direct integration methods in structural dynamics. Amplitude decay and period elongation are used as the basic parameters in order to compare various integration methods. The specific methods studied are the Newmark generalized acceleration scheme, the Houbolt method and the Wilson θ-method. The advantages of each of these methods are discussed. In addition, it is shown how the direct integration of the equations of motion is related to the mode superposition analysis.