Adaptive Reduced Basis Methods Applied to Structural Dynamic Analysis

Abstract

The reduced basis methods (RBM) have been demonstrated as a promising numerical technique for statics problems and are extended to structural dynamic problems in this paper. Direct step-by-step integration and mode superposition are the most widely used methods in the field of the finite element analysis of structural dynamic response and solid mechanics. Herein these two methods are both transformed into reduced forms according to the proposed reduced basis methods. To generate a reduced surrogate model with small size, a greedy algorithm is suggested to construct sample set and reduced basis space adaptively in a prescribed training parameter space. For mode superposition method, the reduced basis space comprises the truncated eigenvectors from generalized eigenvalue problem associated with selected sample parameters. The reduced generalized eigenvalue problem is obtained by the projection of original generalized eigenvalue problem onto the reduced basis space. In the situation of direct integration, the solutions of the original increment formulation corresponding to the sample set are extracted to construct the reduced basis space. The reduced increment formulation is formed by the same method as mode superposition method. Numerical example is given in Section 5 to validate the efficiency of the presented reduced basis methods for structural dynamic problems.

Share and Cite:

Huang, Y. and Huang, Y. (2015) Adaptive Reduced Basis Methods Applied to Structural Dynamic Analysis. American Journal of Computational Mathematics, 5, 317-328. doi: 10.4236/ajcm.2015.53029.

1. Introduction

Nowadays structural dynamic problems are usually solved by the finite element technique. Solution of displacement responses of all the nodes requires great effort. The scale and complexity of dynamics problems of practical engineering structure are ever increasing such that it requests more memory and computing time than before. Despite of the continuing advances in computer speeds and hardware capabilities, the dimension for numerical simulation is too large to provide real-time response in the design, optimization, control and characterization of engineering components or systems. Thus there are many motivations to develop methods that can not only reduce significantly the problem size and computational cost but also retain the accuracy of the solution and the physics of the structures.

Model order reduction techniques [1] -[11] have been proposed to reduce the size of a large-sized model before a detailed analysis performed. They are widely used in global-local analysis, reanalysis and structural dynamic optimization, eigenvalue problem, structural vibration and buckling, sensitivity studies and control parameter design, model update, and damage detection. A detailed review on model reduction techniques can be found in Noor [12] . These reduction methods usually include two steps. The first step is the classic finite element discretization; the second is the computation of some basis vectors in order to perform a Rayleigh-Ritz analysis. Clearly, the success of the method depends chiefly on the proper selection of the basis vectors.

However, order reduction has long been focused on control problems [4] [13] [14] ; most of the reduction methods in that field are designed for small or moderate-size systems and cannot be directly applied in the large- scale case. Nevertheless, the reduced model cannot retain all features of the full model due to the truncated errors. Even for features within an interested frequency range, they may not be exactly kept in the reduced model resulting from most of the model reduction techniques. In recent years, the requirement of reduction techniques for large-scale systems has triggered a revival of research activities in model order reduction [13] [15] [16] . Many powerful reduction techniques have been devised, in particular for linear time-invariant systems. Despite this progress, there are still many open problems.

Different from the traditional reduction methods, the reduced basis method (RBM) [17] - [20] is a very promising method which requires a projection onto the parameter-induced reduced basis space, as makes it very suitable for the analysis of large-scale system. The RBM has first been introduced for single-parameter problems in nonlinear structural analysis in the late 1970s and subsequently developed for multi-parameter problems. However, RBM rarely has been extended to perform model reduction in the structural dynamic problems yet.

In this paper we adopt the reduced basis method to perform the dynamic analysis of structures based on mode superposition method and direct integration method, respectively. A greedy algorithm is suggested to perform the adaptively selection of reduced basis vectors. Numerical example of a simplified one-dimensional seismic model is presented to demonstrate the feasible application of reduced basis method in structural dynamic problems. The error of the reduced system is evaluated numerically.

2. Theoretical Background

In structural dynamic analysis, the equations of motion are generally written as a set of linear second-order differential equations. The matrix form of these equations may be expressed by:

(1)

where, and are the acceleration, velocity, and displacement response vectors of the nodes, respectively, in the total Cartesian coordinate system. The upper dot means derivative with respect to time; is the

equivalent force vector acting on the nodes; the total mass matrix, is the number of elements. The total damping matrix, and the total stiffness matrix. For a finite element:

and, where is the shape function, is the strain matrix,

and is the elasticity matrix. In the following analysis, the structure is subjected to initial conditions given by

(2)

3. The Reduced Basis Method Applied to Dynamic Problems

In the following analysis, the stiffness and mass matrices are assumed as parameter-decomposition forms

(3)

(4)

In Equation (3) and Equation (4), are the numbers of stiffness matrix and mass matrix that can be decomposed, respectively. They are determined by the problem itself.

The damping matrix is considered to be proportional.

(5)

3.1. Reduced Basis Method Based on Model Superposition Technique

3.1.1. Brief Introduction of Mode Superposition Technique

The mode superposition method can be used to perform a time history analysis to obtain the response of structure due to a transient loading as a function of time. It requires the solution of Equation (6) for the frequencies and mode shapes.

(6)

where mode shapes can be shown to be orthogonal to the mass and stiffness matrices, as permit the equations of motion to be uncoupled.

(7)

(8)

The accelerations, velocities, and displacements in Equation (1) are transformed to a different coordinate system:

(9)

Substituting Equation (9) into Equation (1) and premultiplying by yields

(10)

where

―modal stiffness matrix;

―modal damping matrix;

―modal mass matrix;

―modal load vector.

Equation (1) can be decoupled by substituting Equations (7) and (8) into Equation (10).

(11)

where is the loading of the ith order mode, is the damping ratio for the ith mode and is the frequency of the ith order mode as following.

Equation (11) can be solved by a procedure for solving single-degree-of-freedom dynamic problems.

It should be mentioned that the higher mode shapes of the system are unimportant for a practical engineering structure or component. Neglecting the higher frequencies and mode shapes of the system generally does not introduce significant errors. Thus modal truncation is often considered to reduce the computational effort when the number of DOF is large.

3.1.2. Reduced Basis Method to Generalized Eigenvalue Problem

Before the application of reduced basis method, a sample set of parameter domain is selected in a training space, which comprised of parameters spanning the parameter domain roughly.

(12)

The truncated eigenvectors corresponding to the parameters in the sample set are extracted to construct the reduced basis space

(13)

where m is the number of mode is retained in terms of the required accuracy.

It should be noted that the basis vectors are the solutions of the system equations at different parameters. They are perhaps nearly oriented in the same direction. Consequently, the resulted reduced system is very ill-posed especially for large, i.e. the condition number of the reduced stiffness matrix grows exponentially with. To guarantee the basis vectors’ linearly independence and make the reduced system well-posed, QR decomposition is used to generate a new basis which is orthogonal and able to retain all approximation properties of the original basis

(14)

The corresponding transform matrix is.

Then, the eigenvectors corresponding to a new parameter can be expressed as a linear combination of the basis vectors

(15)

The above equation also can be rewritten in a matrix form

(16)

To get the reduced system, the parameter-independent matrices are projected onto the reduced basis in terms of a Galerkin form.

(17)

(18)

From this parameter-decomposition expression, the reduced system can be easily obtained and explored in the whole parameter domain.

(19)

(20)

Obviously, the reduced eigenvalue problem can be solved more efficiently for each new parameter in test parameter-space.

(21)

The truncated eigenvectors can be regenerated approximately by

(22)

The approximation of eigenvalues can be demonstrated in terms of Rayleigh’s quotient.

(23)

3.1.3. The Adaptively Selection of Basis Vectors for Reduced Generalized Eigenvalue Problem

The basis vectors selection is critical for the efficiency and accuracy of the reduced basis method. Too many or too few vectors selected should be avoided. The former results in computational inefficiency, while the latter in unacceptable error. To obtain an appropriate basis space, a greedy algorithm is suggested to select the vectors adaptively.

At first, the error in approximated eigenvalues is presented.

(24)

The maximum error in the training space is definite as

(25)

The performing procedure of greedy algorithm is summarized as follows.

Step 1. One parameter in the training space is selected as the start point; the associated truncated eigenvectors are extracted as the vectors of the reduced basis space.

Step 2. QR decomposition is applied to perform orthogonalization of basis vectors.

Step 3. The reduced generalized eigenvalue problem is solved in the training space to yield the approximated modes,.

Step 4. The maximum error is determined.

Step 5. The truncated eigenvectors corresponding to the maximum error will be selected as the next basis vectors and added to the reduced basis space. Then steps 2 to 4 are repeated. The greedy algorithm will terminate when the maximum error is less than a prescribed tolerance.

3.2. Reduced Basis Method Based on Direct Integration Technique

3.2.1. Brief Introduction of Direct Integration Technique

Direct integration provides a step-by-step numerical procedure to solve the equations of motion in Equation (1) directly without prior transformation of the equations to a different form. It can compute an approximate solution at discrete time intervals, where T is duration of the input motion or loading and is the time step. The widely used explicit methods are only conditionally stable such that some restrictions over the size of the selected time step. On the other hand implicit methods may be unconditionally stable, but the computational work and storage requirement per time step can be much greater than explicit methods because a solution of a coupled algebraic system is always involved. Generally, explicit methods are most preferred for wave propagation problems, while implicit methods are widely employed and advocated for structural dynamic problems.

Newmark method is considered as the example. It is a widely employed linear one-step implicit method with two basic assumptions

(26)

(27)

and is unconditionally stable under the following parameter limitation.

Given approximated values, , , for the displacement, velocity, and acceleration at time, the algorithmic values at time are the solution of the linear algebraic equations.

(28)

From the initial condition given by Equation (2), the initial acceleration given by

(29)

3.2.2. Reduced Basis Method Based on Direct Integration Technique

Just as the same in mode superposition, a sample set will be introduced from a training space comprised of a span of parameters and all time steps. The reduced basis space is defined as the span of finite element displacement response and doesn’t change with time for a new arbitrary parameter.

(30)

For the same reason mentioned in foregoing section, QR decomposition is applied to generate an orthogonal reduced basis space.

(31)

The transform matrix for projection can be written as:

(32)

The displacement response corresponding to new parameter and new time step can be approximated as the linear combination of the vectors in the reduced basis space.

(33)

It also can be expressed as a matrix form.

(34)

The approximated velocity and acceleration can be obtained by first order and second order derivatives of the approximated displacement response with respect to time, respectively.

(35)

(36)

The reduced Newmark formulation can be obtained by Galerkin projection of original space onto the reduced basis space.

(37)

The reduced stiffness, mass and damping matrices are respectively given by

(38)

(39)

(40)

where the reduced parameter-independent matrices are

(41)

(42)

The reduced load vector is

(43)

The initial condition corresponding to the reduced system is

(44)

3.2.3. The Adaptively Selection of Basis Vectors for Reduced Newmark Formulation

Similarly, the greedy algorithm is adopted to select the vectors adaptively and subsequently obtain an appropriate reduced basis space.

At first, the projection error is defined in the training space as

(45)

where is the projection displacement,.

The maximum norm of the projection error is defined as

(46)

The perform procedure of greedy algorithm is summarized as follows.

Step 1. To span the training space, the displacement of the last time step is selected as the first basis vector, corresponding to one source in the training space.

Step 2. QR decomposition is applied to perform orthogonalization of basis vectors.

Step 3. The reduced Newmark’s formula is solved in the training space to yield the reduced basis displacements.

Step 4. The maximum norm of the projection error is determined.

Step 5. The displacement corresponding to the maximum norm of the projection error will be selected as the next basis vector and steps 2 to 4 are repeated. The greedy algorithm will terminate when the maximum norm of projection error is less than a prescribed tolerance.

4. Numerical Example

4.1. Numerical Model

A simplified one-dimensional seismic model [21] is presented to numerically validate the application of reduced basis method to structural dynamic problem. The pressure variable during an earthquake is governed by dynamic equilibrium equation:

(47)

The earthquake source and occurring time are considered as system parameters, and denoted by, which vary within the domain. and, which

denote the spatial distribution and the temporal characteristics of earthquake source respectively, are showed in Figure 1 and Figure 2, respectively. The spatial domain is divided into linear elements and normalized to unit length,. The pressure changes with the occurring time and doesn’t change in spatial distribution such that it is fixed as 0.5. The pressure is zero in the earth’s crust and the pressure gradient is zero on the earth’s surface. The initial condition is given by

(48)

To obtain parameter-decomposition forms of stiffness and mass matrices, namely, to express the stiffness and

Figure 1. The time history of earthquake source.

Figure 2. The spatial distribution of earthquake source.

mass matrices as the combination form of product of system parameter function and the matrix independent of

system parameters. The original x-domain is decomposed into the left zone, forcing zone, right zone and output zone. A standard y-domain is introduced as reference and decomposed into

A piecewise affine mapping from the standard y-domain to the original x-domain is given in Figure 3: from to; from to; from to; and the identity mapping from to. The resultant parameter-decomposition matrices are

(49)

(50)

where and are matrices, and are both independent of the system parameters, the parameter-dependent coefficients are:

4.2. The Numerical Results

As the reduced structural dynamic analysis performed by using mode superposition, the 12th truncation of mode is considered. It can be found from Figure 4 that the maximum error of approximated eigenvalue decreases rapidly with the increasing of the basis vectors. For Newmark integration case, the numerical parameters are selected as and the same convergence phenomenon as mode superposition can be found in Figure 5. However, it is obvious that the former converges sooner than the later.

The resulted reduced eigenvalue problem is 60 in dimensional, while the reduced Newmark formulation is 85 in dimensional for a prescribed error tolerance. Figure 6 and Figure 7 show that both reduced mode superposition and reduced Newmark integration approximate the original algorithms very well for engineering analysis. The CPU time for the reduced system and the original system are given in Table 1. The original Newmark costs expensively CPU time, while the reduced Newmark gives the lowest cost for test parameter

Figure 3. The affine mapping from y-domain to x-domain.

Figure 4. Maximum norm of projection error changing with the increasing of number of basis vectors.

Figure 5. Maximum norm of projection error changing with the increasing of number of basis vectors.

Figure 6. Maximum norm of projection error changing with the increasing of number of basis vectors.

Table 1. CPU time comparison of mode superposition and implicit Newmark method.

Figure 7. The pressure in the output domain corresponding to S = 0.5.

space. Despite the original mode superposition more effectively executed than the original Newmark method, the reduced form of the former costs more CPU time than the reduced form of the later. It can be concluded that the dynamic analysis have been performed much more effectively by either reduced mode superposition or reduced Newmark method.

It should be point out that the dimensional of the reduced system is determined by the reduced basis space and independent of the original system. For larger dynamic system, the efficiency of the reduced basis methods can be further enhanced.

5. Conclusion

Two kinds of reduced basis methods for dynamic problems are proposed in this paper. In the numerical example, the direct integration for the dynamic analysis is not numerically efficient as compared with the mode superposition method using eigenvectors due to the linear property of the seismic problem. However, it proves that the reduced basis method is available for structural dynamic analysis based on either mode superposition or direct integration. Though the undamped case studied, the reduced basis method can be applied to damped structures without any more effort. Furthermore, although the reduced Newmark method is only considered here, the reduced basis method can be very easily extended to other direct integration techniques.

Acknowledgements

This project is supported by National Natural Science Foundation of China (Grant No. 51305045), and by China Postdoctoral Science Foundation (No. 2014M562099).

Conflicts of Interest

The authors declare no conflicts of interest.

References

[1] Guyan, R.J. (1965) Reduction of Stiffness and Mass Matrices. AIAA Journal, 3, 380-381.
http://dx.doi.org/10.2514/3.2874
[2] Wilson, E.L. and Bayo, E.P. (1967) Use of Special Ritz Vectors in Dynamic Substructure Analysis. AIAA Journal, 5, 1944-1954.
[3] Sirovich, L. and Kirby, M. (1987) Low-Dimensional Procedure for the Characterization of Human Faces. Journal of the Optical Society of America A, 4, 519-524.
http://dx.doi.org/10.1364/JOSAA.4.000519
[4] Moore, B.C. (1981) Principal Component Analysis in Linear Systems: Controllability, Observability, and Model Reduction. IEEE Transactions on Automatic Control, 26, 17-32.
http://dx.doi.org/10.1109/TAC.1981.1102568
[5] Friswell, M.I., Garvey, S.D. and Penny, J.E.T. (1995) Model Reduction Using Dynamic and Iterative IRS Technique. Journal of Sound and Vibration, 186, 311-323.
http://dx.doi.org/10.1006/jsvi.1995.0451
[6] Kammer, D.C. (1987) Test-Analysis Model Development Using an Exact Modal Reduction. The International Journal of Analytical and Experimental Modal Analysis, 2, 174-179.
[7] Leung, A.Y.T. (1978) An Accurate Method of Dynamic Condensation in Structural Analysis. International Journal for Numerical Methods in Engineering, 12, 1705-1715.
http://dx.doi.org/10.1002/nme.1620121108
[8] O’Callahan, J.C. (1989) A Procedure for an Improved Reduced System (IRS) Model. Proceedings of the 7th International Modal Analysis Conference, Las Vegas, 17-21.
[9] Singh, M.P. and Suarez, L.E. (1992) Dynamic Condensation with Synthesis of Substructure Eigenproperties. Journal of Sound and Vibration, 159, 581-606.
http://dx.doi.org/10.1016/0022-460X(92)90455-7
[10] Van Woerkom, P.Th.L.M. (1990) Mathematical Models of Flexible Spacecraft Dynamics: A Survey of Order Reduction Approaches. Control Theory and Advanced Technology, 4, 609-632.
[11] Qu, Z.Q. and Chang, W.J. (1997) An Extended Order Method (EOM) for Nonlinear Eigenvalue Problems. Mechanics and Practice, 19, 25-26.
[12] Noor, A.K. (1994) Recent Advances and Applications of Reduction Methods. Applied Mechanics Reviews, 47, 125-146.
http://dx.doi.org/10.1115/1.3111075
[13] Lall, S., Marsden, J.E. and Glavaski, S. (2002) A Subspace Approach to Balanced Truncation for Model Reduction of Nonlinear Control Systems. International Journal of Robust and Nonlinear Control, 12, 519-535.
http://dx.doi.org/10.1002/rnc.657
[14] Ito, K. and Ravindran, S.S. (1997) Reduced Order Methods for Nonlinear Infinite Dimensional Control Systems. Proceedings of the 36th conference on Decision & Control, San Diego, 10-12 December 1997, 2213-2218.
http://dx.doi.org/10.1109/CDC.1997.657101
[15] Willcox, K. and Peraire, J. (2002) Balanced Model Reduction via the Proper Orthogonal Decomposition. AIAA Journal, 40, 2323-2330.
http://dx.doi.org/10.2514/2.1570
[16] Qu, Z.Q. (2002) Model Reduction for Dynamical Systems with Local Nonlinearities. AIAA Journal, 40, 327-333.
http://dx.doi.org/10.2514/2.1649
[17] Almroth, B.O., Stern, P. and Brogan, F.A. (1978) Automatic Choice of Global Shape Functions in Structural Analysis. AIAA Journal, 16, 525-528.
http://dx.doi.org/10.2514/3.7539
[18] Noor, A.K. and Peters, J.M. (1980) Reduced Basis Technique for Nonlinear Analysis of Structures. AIAA Journal, 18, 455-462.
http://dx.doi.org/10.2514/3.50778
[19] Machiels, L., Maday, Y., Oliveira, I.B., Patera, A.T. and Rovas, D.V. (2003) Output Bounds for Reduced-Basis Approximations of Symmetric Positive Definite Eigenvalue Problems. Comptes Rendus de l’Académie des Sciences - Series I - Mathematics, 331, 152-158.
[20] Maday, Y., Patera, A.T. and Peraire, J. (1999) A General Formulation for a Posteriori Bounds for Output Functionals of Partial Differential Equations; Application to the Eigenvalue Problem. Comptes Rendus de l’Académie des Sciences— Series I—Mathematics, 328, 823-828.
[21] Akcelik, V., Bielak, J., Biros, G., Epanomeritakis, I., Fernandez, A., Ghattas, O., et al. (2003) High Resolution Forward and Inverse Earthquake Modeling on Terasacale Computers. Proceedings of SC2003, Phoenix, 15-21 November 2003.

Copyright © 2024 by authors and Scientific Research Publishing Inc.

Creative Commons License

This work and the related PDF file are licensed under a Creative Commons Attribution 4.0 International License.