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.     

Ritz vectors and generation criteria for mode superposition analysis

A generation procedure of Ritz vectors to control the inclusion of static effect and the number of vectors in mode superposition dynamic analysis is presented. The original algorithm of the Ritz vectors¹⁵ is modified to improve stability in the generation procedure and to include the use of static residual. To reject unimportant Ritz vectors, cut-off criteria, which are based on the participation of mass distribution and spatial load distribution, are proposed. Numerical examples are presented to illustrate the effectiveness of the derived Ritz vectors over the eigenvectors and the performance of the cutoff criteria in the mode superposition dynamic analysis.     

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.     

非比例阻尼线性体系地震反应计算的振型分解反应谱法

以非比例阻尼线性体系地震反应计算实数形式的一般解答为基础,推导得到了非比例阻尼线性体系水平地震作用计算的多种形式,建立了非比例阻尼线性体系地震反应计算振型分解反应谱法的基本过程与步骤。最后,以一个五层剪切型结构为例,通过与各种常用直接积分方法计算结果的比较,证实了本文非比例阻尼线性体系地震反应计算实数形式的一般解答的高精度与可靠性。通过对多种形式地震作用所得地震效应的比较,证实了非比例阻尼线性体系地震反应振型分解反应谱方法的可靠性及可行性。     

A subspace modal superposition method for non-classically damped systems

For structures with non-proportional damping, complex eigenvectors or mode shapes must be used in order to decoe the equations of motion. The resulting equations can then be solved in a systematic way. The necessity of solvie complex eigenvalue problem of a large system remains an obstacle for the practical application of the method. This stres utilizes the fact that in practice only a small number of the complex modes are needed. Therefore, these complex modes be approximated by a linear combination of a small number of the undamped modes, which can be obtained by established methods with less cost. An additional eigenvalue problem is then solved in a subspace with a much sm dimension to provide the best combination coefficient for each complex mode. The method of solution for the decoue equations is then carried over, using the approximate complex modes expressed in undamped mode shapes, to resue simple formulas for the time- and frequency-domain solution. Thus, an efficient modal superposition method is develoe for non-proportionally damped systems. The accuracy of this approximate method is studied through an example. Comparing the frequency response result using the approximate method with that using the exact complex modes, found that the error is negligible.     

Lanczos method for dynamic analysis of damped structural systems

In this paper we extend the Lanczos algorithm for the dynamic analysis of structures⁷ to systems with general matrix coefficients. The equations of dynamic equilibrium are first transformed to a system of first order differential equations. Then the unsymmetric Lanczos method is used to generate two sets of vectors. These vectors are used in a method of weighted residuals to reduce the equations of motion to a small unsymmetric tridiagonal system. The algorithm is further simplified for systems of equations with symmetric matrices. By appropriate choice of the starting vectors we obtain an implementation of the Lanczos method that is remarkably close to that in Reference 7, but generalized to the case with indefinite matrix coefficients. This simplification eliminates one of the sets of vectors generated by the unsymmetric Lanczos method and results in a symmetric tridiagonal, but indefinite, system. We identify the difficulties that may arise when this implementation is applied to problems with symmetric indefinite matrices such as vibration of structures with velocity feedback control forces which lead to symmetric damping matrices. This approach is used to evaluate the vibration response of a damped beam problem and a space mast structure with symmetric damping matrix arising from velocity feedback control forces. In both problems, accurate solutions were obtained with as few as 20 Lanczos vectors.     

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.     

Response of non-classically damped structures in the modal subspace

The evaluation of the dynamic response of non-classically damped linear structures requires the solution of an eigenproblem with complex eigenvalues and modal shapes. Since in practice only a small number of complex modes are needed, the complex eigenvalue problem is solved in the modal subspace in which the generalized damping matrix is not uncoupled by classical real modes. It follows that the evaluation of the structural response requires in both cases the determination of complex modes by numerical techniques, which are not as robust as techniques currently used for the solution of the real eigenvalue problem, and the use of complex algebra. In the present paper an unconditionally stable step-by-step procedure is presented for the response of non-classically damped structures in the modal subspace without using complex quantities. The method is based on the evaluation of the fundamental operator in approximated form of the numerical procedure. In addition, the method can be easily modified to incorporate the modal superposition pseudo-static correction terms.     

On precise time integration method for non-classically damped MDOF systems

In the complex mode superposition method, the equations of motion for non-classically damped multiple-degree-of-freedom （MDOF） discrete systems can be transferred into a combination of some generalized SDOF complex oscillators. Based on the state space theory, a precise recurrence relationship for these complex oscillators is set up; then a delicate general solution of non-classically damped MDOF systems, completely in real value form, is presented in this paper. In the proposed method, no calculation of the matrix exponential function is needed and the algorithm is unconditionally stable. A numerical example is given to demonstrate the validity and efficiency of the proposed method.     

Dynamic analysis of structures using lanczos co-ordinates

A procedure for deriving the Lanczos vectors is explained and their use in structural dynamics analysis as an alternative to modal co-ordinates is discussed. The vectors are obtained by an inverse iteration procedure in which orthogonality is imposed between the vectors resulting from successive iteration cycles. Using these Lanczos vectors the equations of motion are transformed to tridiagonal form, which provides for a very efficient time-stepping solution. The effectiveness of the method is demonstrated by a numerical example.     

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.     

Short communication block lanczos method for dynamic analysis of structures

The simple Lanczos method presented in a recent paper by the writers, with application to single vector loads, is extended to include a more general dynamic loading represented as a linear combination of k vectors (load patterns). The result is a set of orthogonal vectors that is used to transform the equations of motion to a banded form, the half-bandwidth of which becomes k + 1. When k is small relative to the number of equations, this approach provides for a very efficient time-stepping solution.     

Complex complete quadratic combination method for damped system with repeated eigenvalues

A new response-spectrum mode superposition method, entirely in real value form, is developed to analyze the maximum structural response under earthquake ground motion for generally damped linear systems with repeated eigenvalues and defective eigenvectors. This algorithm has clear physical concepts and is similar to the complex complete quadratic combination (CCQC) method previously established. Since it can consider the effect of repeated eigenvalues, it is called the CCQC-R method, in which the correlation coefficients of high-order modal responses are enclosed in addition to the correlation coefficients in the normal CCQC method. As a result, the formulas for calculating the correlation coefficients of high-order modal responses are deduced in this study, including displacement, velocity and velocity-displacement correlation coefficients. Furthermore, the relationship between high-order displacement and velocity covariance is derived to make the CCQC-R algorithm only relevant to the high-order displacement response spectrum. Finally, a practical step-by-step integration procedure for calculating high-order displacement response spectrum is obtained by changing the earthquake ground motion input, which is evaluated by comparing it to the theory solution under the sine-wave input. The method derived here is suitable for generally linear systems with classical or non-classical damping.     

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.     

Earthquake response analysis considering non-proportional damping

Non-proportional damping may be defined as a form of linear viscous damping which introduces coupling between the undamped modal co-ordinate equations of motion. The standard mode superposition method of earthquake response analysis therefore cannot be employed with non-proportionally damped structures. In this paper, several methods for analysing the dynamic response of non-proportional damped structures are outlined. It is concluded that the most efficient procedure is to express the response in terms of a truncated set of undamped modal coordinates and to integrate directly the resulting coupled equations. The effectiveness of the method is demonstrated by a numerical example.     

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.     

Partial eigensolution of damped structural systems by Arnoldi's method

An efficient numerical algorithm is developed to solve the quadratic eigenvalue problems arising in the dynamic analysis of damped structural systems. The algorithm can even be applied to structural systems with non-symmetric matrices. The algorithm is based on the use of Arnoldi's method to generate a Krylov subspace of trial vectors, which is then used to reduce a large eigenvalue problem to a much smaller one. The reduced eigenvalue problem is solved and the solutions are used to construct approximate solutions to the original large system. In the process, the algorithm takes full advantage of the sparseness and symmetry of the system matrices and requires no complex arithmetic, therefore, making it very economical for use in solving large problems. The numerical results from test examples are presented to demonstrate that a large fraction of the approximate solutions calculated are very accurate, indicating that the algorithm is highly effective for extracting a number of vibration modes for a large dynamic system, whether it is lightly or heavily damped.     

Spectral analysis of systems with non-classical damping using classical mode superposition technique

A spectral method for random vibration analysis of a structural system with non-proportional damping is presented using classical (undamped) mode superposition technique. The method obtains the frequency response function of the system by solving the dynamic equilibrium equations in generalized co-ordinates through an iterative process. The iterative solution is written in closed form and the proof for convergence of the iterative process is given. Numerical examples show the convergence characteristics of the process and an excellent accuracy of the obtained results. The method turns out to be computationally more efficient than the conventional methods of spectral analysis using damped mode shapes and frequencies.     

Closed‐form solution for seismic response of adjacent buildings with linear quadratic Gaussian controllers

Closed‐form solution for seismic response of adjacent buildings connected by hydraulic actuators with linear quadratic Gaussian (LQG) controllers is presented in this paper. The equations of motion of actively controlled adjacent buildings against earthquake are first established. The complex modal superposition method is then used to determine dynamic characteristics, including modal damping ratio, of actively controlled adjacent buildings. The closed‐form solution for seismic response of the system is finally derived in terms of the complex dynamic characteristics, the pseudo‐excitation method and the residue theorem. By using the closed‐form solution, extensive parametric studies can be carried out for the system of many degrees of freedom. The beneficial parameters of LQG controllers for achieving the maximum response reduction of both buildings using reasonable control forces can be identified. The effectiveness of LQG controllers for this particular application is evaluated in this study. The results show that for the adjacent buildings of different dynamic properties, if the parameters of LQG controllers are selected appropriately, the modal damping ratios of the system can be significantly increased and the seismic responses of both buildings can be considerably reduced. Copyright © 2001 John Wiley & Sons, Ltd.     

三维感应测井响应计算的交错网格有限差分法

应用交错网格有限差分法计算三维复杂环境中的感应测井响应. 其中，利用Krylov子空间不变性求解离散得到的大型稀疏复对称线性方程组. 在构造Krylov子空间时使用其系数矩阵的伪逆以改善迭代的收敛性. 迭代中，使用不完全Cholesky分解共轭梯度法求解4个三维Poisson方程以得到新的Lanczos向量. 通常迭代不超过20次可得到理想结果. 另外，提出一种新的物质平均公式以计算电导率平均值，可保证电流守恒.