9512.net

甜梦文库

甜梦文库

当前位置：首页 >> >> # PHYSICAL INTERPRETATION OF THE PROPER ORTHOGONAL MODES USING THE SINGULAR VALUE DECOMPOSITI

Journal of Sound and <ibration (2002) 000(0), 000}000 doi:10.1006/jsvi.2001.3930, available online at http://www.idealibrary.com on

PHYSICAL INTERPRETATION OF THE PROPER ORTHOGONAL MODES USING THE SINGULAR VALUE DECOMPOSITION

G. KERSCHEN AND J. C. GOLINVAL

??AS - <ibrations et Identi,cation des Structures, ;niversite de ?iege, Chemin des Chevreuils, H ` 1 (Batiment B52) 4000 ?iege, Belgium. E-mail: g.kerschen@ulg.ac.be ` (Received 12 November 2000, and in ,nal form 6 April 2001) Proper orthogonal decomposition is a statistical pattern analysis technique for "nding the dominant structures, called the proper orthogonal modes, in an ensemble of spatially distributed data. While the proper orthogonal modes are obtained through a statistical formulation, they can be physically interpreted in the "eld of structural dynamics. The purpose of this paper is thus to provide some insights into the physical interpretation of the proper orthogonal modes using the singular value decomposition 2002 Academic Press

1. INTRODUCTION

Proper orthogonal decomposition (POD) is a procedure for extracting a basis for a modal decomposition from an ensemble of signals. A very appealing property of the POD is its optimality. Among all possible decompositions of a random "eld, the POD is the most e$cient in the sense that for a given number of modes, the projection on the subspace used for modelling the random "eld will on average contain the most energy possible. Although POD has been regularly applied to non-linear problems, it is essential to underline that it is a linear technique and that it is optimal only with respect to other linear representations. The applications of this procedure are extensive in modelling of turbulence [1, 2] and image processing [3], and POD is now emerging as a useful tool in the "eld of structural dynamics. For instance, it has been applied to estimate the dimensionality of a system [4], to build reduced order models [5, 6], and to the identi"cation and updating of non-linear systems [7}9]. The purpose of this paper is to determine whether a physical interpretation can be attributed to the modes obtained from the decomposition, i.e., the proper orthogonal modes (POMs). Particularly, it is inquired when the POMs are related to the vibration eigenmodes. This work is closely related to the paper of Feeny and Kappagantu [10]. However, in the present paper, the emphasis is shifted towards the singular value decomposition of the displacement matrix rather than the eigenvalue problem of the covariance matrix. Furthermore, the case of linear systems under harmonic and white noise excitations is discussed in greater detail. The paper is organized as follows. In Section 2, the POD is brie#y introduced. Section 3 gives a brief review of the singular value decomposition and its properties that are relevant in the context of this paper. Sections 4, 5 and 6 study the physical interpretation of the POMs of discrete linear systems, respectively, for the free response in the undamped and damped cases, and for the harmonic response. Section 7 o!ers a geometric approach to the

0022-460X/02/000000#17 $35.00/0 2002 Academic Press

PROD.TYPE: TYP

ED: RAVI PAGN: TJ GSRS N SCAN:PREETHY

JSVI 20013930

2

G. KERSCHEN AND J. C. GOLINVAL

comparison between vibration eigenmodes and POMs. It also investigates the relationship between non-linear normal modes (NNMs) and POMs. Finally, the discussion of the stationary random response of a linear system to a white noise excitation is included in Appendix A.

2. PROPER ORTHOGONAL DECOMPOSITION

Proper orthogonal decomposition, also known as Karhunen}Loeve transform, was introduced by Kosambi [11]. It is also worth pointing out that POD is closely related to principal component analysis (PCA) introduced by Hotelling [12]. For a detailed historical review of POD or PCA, the reader is referred to references [10, 13]. Let v(x, t) be a zero mean random "eld on a domain . In practice, the "eld is sampled at a "nite number of points in time. Then, at time t , the system displays a snapshot v (x) which G G is a continuous function of x in . The aim of the POD is to "nd the most persistent structure

(x) among the ensemble of n snapshots. This is equivalent to minimizing the objective function : L " (

(x)!v (x)) ?x3 . G G

Minimize

(1)

Equation (1) can also be written in terms of a maximization problem [6]: (1/N) , (

(x)v (x) d ) L L "

(x)

(x) d ?x3 .

Maximize

(2)

Finally, the optimization problem can be reduced to the following integral eigenvalue problem [6]:

K(x, x)

(x) dx"

(x),

where K is the two-point correlation function 1 L K(x, x)" v (x)v (x). G G n G

(3)

(4)

Equation (3) has a "nite number of orthogonal solutions

G(x), called the proper orthogonal modes (POMs) with corresponding real and positive eigenvalues G. In practice, the snapshots are available at discrete measurement points x where k"1,2, m and the I integral eigenvalue problem (3) reduces to "nd the eigensolution of an (m;m) space correlation tensor K(x , x ) G" 2 2 2 K(x , x ) K 2 .

(5)

K(x , x ) 2 K(x , x ) K K K

JSVI 20013930

PROPER ORTHOGONAL MODE DECOMPOSITION

3

To summarize, if the responses (e.g., the displacements) q (t) of a discrete dynamical I system with m degrees of freedom (d.o.f.) are sampled n times and if the (m;n) matrix q (t ) Q" 2 q (t ) L 2 2 "[q(t )2q(t )] K q (t ) 2 q (t ) K K L 2

(6)

is formed, then the POMs are merely the eigenvectors of G"(1/n)QQ2 and the corresponding eigenvalues are the proper orthogonal values (POVs). A POV measures the relative energy of the system dynamics contained in the associated POM.

3. SINGULAR VALUE DECOMPOSITION

The objective of this section is to review the singular value decomposition (SVD) and its features that are relevant in the context of POD. Particularly, it is pointed out that the POMs are optimal with respect to energy content. For a detailed description of SVD and its several possible applications in structural dynamics, the reader is referred to references [14, 15]. Since the matrices considered throughout the paper are built from system responses, e.g., displacements, the discussion is restricted to real matrices only. For any real (m;n) matrix A, there exists a real factorization A"U V2, (7)

where U is an (m;m) orthonormal matrix. Its columns form the left singular vectors. is an (m;n) pseudo-diagonal and semi-positive-de"nite matrix with diagonal entries containing the singular values . V is an (n;n) orthonormal matrix. Its columns form the right G singular vectors.

3.1.

GEOMETRIC INTERPRETATION

The SVD of a matrix, seen as a collection of column vectors, provides important insight into the oriented energy distribution of this set of vectors. It is worth recalling that 1. the energy of a vector sequence a building an (m;n) matrix A is de"ned via the I Frobenius norm K L N where p"min(m, n), E(A)"#A#" a " $ GH I G H I (8)

so that the energy of a vector sequence is equal to the energy in its singular spectrum; 2. the oriented energy of a vector sequence in some direction p with unit vector e of the N m-dimensional column space is the sum of squared projections of the vectors on to direction p L E (A)" (e2 a ). N N I I (9)

JSVI 20013930

4

G. KERSCHEN AND J. C. GOLINVAL

One essential property of SVD is that extrema in this oriented energy distribution occur at each left singular direction [15]. The oriented energy measured in the direction of the ith left singular vector is equal to the ith singular value squared. Since the POMs are directly related to the left singular vectors, it can be stated that they are optimal with respect to energy content in a least-square sense, i.e., they capture more energy per mode than any other set of basis functions.

3.2.

RELATION WITH THE EIGENVALUE PROBLEM

The SVD of a matrix can be calculated by means of solving two eigenvalue problems, or even one if only the left or the right singular vectors are required. Indeed, AA2"U U2 A2A"V V2. (10)

Consequently, the singular values of A are found to be the square roots of the eigenvalues of AA2 or A2A. The left and right singular vectors of A are the eigenvectors of AA2 and A2A respectively. Applying this reasoning to POD, it is now clear that the POMs, de"ned as the eigenvectors of the covariance matrix G"(1/n)AA2, are the left singular vectors of A. The POVs, de"ned as the eigenvalues of the covariance matrix, are the square of the singular values divided by the number of samples n. In conclusion, POD can be carried out directly by means of an SVD of matrix A. An interesting interpretation of the eigenvalue problem is that if a matrix is real, symmetric and positive de"nite, then the eigenvectors of the matrix are the principal axes of the associated quadratic form which is an n-dimensional ellipsoid centered at the origin of the Euclidean space [16]. Since AA2 is real, symmetric and positive de"nite, the POMs as eigenvectors of the covariance matrix are the principal axes of the family of ellipsoids de"ned by y2Gy"c where y is a real non-zero vector and c is a positive constant. It is worth pointing out that Feeny and Kappagantu showed that if each data has unit mass, then the POMs are the principal axes of inertia [10].

4. UNDAMPED AND UNFORCED LINEAR SYSTEMS

The aim of this section is to "nd the existing relationships between the POMs and the eigenmodes of an undamped and unforced linear system with m. d.o.f. The equation of motion may be written as follows: Mq#Kq"0, K (11)

where M and K are the mass and sti!ness matrices, respectively, and q is the vector of displacement co-ordinates. The system response due to initial conditions may be expressed as K K q(t)" (A cos t#B sin t)x " e (t)x , G G G G G G G G G

(12)

JSVI 20013930

PROPER ORTHOGONAL MODE DECOMPOSITION

5

where , x are the natural frequencies (in rad/s) and eigenmodes of the system; A and B G G G G are constants depending on the initial conditions; and e (t)"A cos t#B sin t G G G G G represents the time modulation of mode x . G The time discretization of the system response leads to n sampled values of the time functions which form an (m;n) matrix whose columns are the members of the data ensemble Q"[q(t )2q(t )] L " which can also be written as e (t ) 2 e (t ) K 2 e (t ) L 2 K e (t )x 2 e (t )x , G G G L G G G K (13)

Q"[x 2x ] K

2 e (t ) K L

"[x 2x K

e2 ] 2 e2 K

"[x 2x ] [e 2e ]2 K K "XE2 "X[I Z] [E R]2, (14)

where X is the (m;m) modal matrix whose columns are the eigenmodes of the system; E is an (n;m) matrix whose columns are the functions e (t) at times t ,2, t ; I is an (m;m) G L identity matrix; Z is an (m;(n!m)) matrix full of zeros; R is an (n;(n!m)) matrix; and e "[e (t )2e (t )]2. G G G L Attention should be paid to the fact that R does not in#uence Q since it is multiplied by a matrix full of zeros. Equation (14) can be expressed in a more familiar form as Q"[X] [I Z] [E R]2"U V2. GHI GHI GHI

U V2

(15)

Accordingly, the above decomposition of Q may be thought of as the SVD of this matrix. However, this decomposition requires matrices U and V to be orthonormal as mentioned in Section 3. The aim now is to "nd the conditions when the columns of U(,X) and V(,[E R]) are orthogonal. 1. The columns of U are formed by the eigenmodes of the structure. The eigenmodes are orthogonal to each other in the metrics of the mass and sti!ness matrices. If the mass matrix is proportional to the identity matrix, it turns out that x 2 x " . G H GH Consequently, X is orthogonal if the mass matrix is proportional to the identity matrix.

JSVI 20013930

6

G. KERSCHEN AND J. C. GOLINVAL

2. It remains to determine when the columns of V are orthogonal. For this purpose, equation (14) may be rewritten as follows: Q"X[I Z] [E R]2 "X[diag(#e #) Z] [E diag(#e #\) R]2 G G #e # 0 2 0 #e # 2 "[x 2x ] K 2 2 2 0 0 0 0 2 0 0 0 2 2 0 0 0 e e 2 e 2 K R . #e # #e # #e # K (16)

2 2 2 2

2 #e # K

If the natural frequencies are distinct, it can be easily argued that the columns of G E diag(#e #\) are orthogonal if we consider an in"nite set of sampled values, i.e., G e e H P0 if nP iOj. G R, #e # #e # H G (17)

Since R does not have an in#uence on Q, its columns can be computed in order that they are orthogonal to those of E diag(#e #\). As can also be seen from equation (16), POD is G a bi-orthogonal decomposition that uncouples the spatial and temporal information contained in the data. To summarize, if the mass matrix is proportional to the identity matrix and if the number of samples is in"nite, the singular value decomposition of Q is such that (1) the columns of U are the eigenmodes; (2) the "rst n columns of V are the normalized time modulations of the modes. As stated in section 3.2, the POD basis vectors are just the columns of the matrix U in the singular value decomposition of the displacement matrix. Therefore, it can be concluded that the POMs converge to the eigenmodes of an undamped and unforced linear system whose mass matrix is proportional to identity if a su.cient number of samples is considered. Feeny and Kappagantu [10] previously obtained the same conclusion by a di!erent way. They based their demonstration on the fact that the POMs are the eigenvectors of the covariance matrix. In the case of a mass matrix not proportional to identity, the POMs no longer converge to the eigenmodes since the former are orthogonal to each other while the latter are orthogonal with respect to the mass matrix. However, knowing the mass matrix, it is still possible to retrieve the eigenmodes from the POMs. Equation (11) has to be rewritten through the co-ordinate transformation q"M\p as p#M\KM\p"0. K (18)

In equation (18), the system matrices are still symmetric while the e!ective mass matrix is equal to the identity. Thus, the left singular vectors of P"[p(t )2p(t )], i.e., the POMs, L converge to the eigenmodes y of this system. It is a simple matter to demonstrate that the G eigenmodes x of system (11) are related to those of system (18) by the following G relationship: x "M\y . G G (19)

JSVI 20013930

PROPER ORTHOGONAL MODE DECOMPOSITION

7

This section has investigated the discrete case. A detailed study of distributed systems can be found in reference [17]. This paper underlines that the conclusions are still valid if the distributed system is uniformly discretized.

5. DAMPED AND UNFORCED LINEAR SYSTEMS

Consider now a damped but still unforced linear system with m. d.o.f. for which the equation of motion is given as follows: Mq#Cq#Kq"0. K (20)

If the structure is lightly damped or with the assumption of modal damping, the system response can be readily written as K K e (t)x . q(t)" A exp\CG SGR cos((1! t# )x " G G G G G G G G G Using the same procedure as in the previous section yields Q"[q(t )2q(t )] L " K e (t )x 2 e (t )x G G G L G G G e2 ] 2 e2 K (22) K (21)

"[x 2x K

"XE2 "X[I Z] [E R]2 "X[diag(#e #) Z] [E diag(#e #\) R]2 G G #e # 0 2 0 #e # 2 "[x 2x ] K 2 2 2 0 "U V2, where e "[A exp\CG SG R cos((1! t # )2A exp\CG SG R cos((1! t # )]2. G G G G G G G G L G Again, the columns of U(,X) are orthogonal if the mass matrix is proportional to the identity matrix. The main di!erence with the undamped case is that the time modulations e (t)P0 if tP since the system returns to the equilibrium position in a "nite time. R G Consequently, it can no longer be a$rmed that #e #P if nP and that the columns of R R G E diag(#e #\) are orthogonal to each other. This causes a set of POMs di!erent from the G 0 0 0 2 0 0 0 2 2 0 0 0 e e 2 e 2 K R . #e # #e # #e # K

2 2 2 2

2 #e # K

JSVI 20013930

8

G. KERSCHEN AND J. C. GOLINVAL

eigenmodes to be obtained. However, if the damping is low and if a su$cient number of points are considered, E diag(#e #\) is almost orthogonal. In conclusion, the POMs of G a lightly damped and unforced linear system are a very good approximation of the eigenmodes of this system. This is in accordance with the result obtained in reference [10] using the eigensolution perspective.

6. HARMONIC AND FORCED HARMONIC RESPONSES OF A LINEAR SYSTEM

This section is divided into two parts. Firstly, the harmonic response of a linear system is considered. By harmonic response, we mean the combination of the free and forced responses. Secondly, attention is focused only on the forced response of the linear system.

6.1.

HARMONIC RESPONSE

The equation of motion of a linear system with m. d.o.f. excited by an harmonic force with a constant amplitude is Mq#Kq"f sin K C t. (23)

Equation (23) may be transformed by considering a new variable s"sin t that accounts C for the harmonic force Mq#Kq"fs K with s(0)"0, sK # s"0 C which yields

CFDFE

M* K*

sR (0)"

C

,

(24)

CFFDFFE

M 0 0 1

K q K # sK 0

!f C

q 0 " . s 0

(25)

For the sake of clarity, note ( , x ) the eigensolutions of the initial system (23) and G G ( *, x * ) the eigensolutions of the transformed system (24). This latter system may be G G viewed as an unforced system with m#1 d.o.f. (25). If the mass matrix is proportional to identity and if the number of samples is large enough, section 4 allows us to conclude that the POMs of the transformed system response converge to the eigenmodes of that system. Let us now compute the eigenmodes of the transformed system. These are the solution of (K*! * M*)x * "0 G G if * is a root of the algebraic equation G det(K*! M*"det K !f 0 C ! M 0 "0. 0 1 (27) (26)

JSVI 20013930

PROPER ORTHOGONAL MODE DECOMPOSITION

9

This equation becomes K! M 0 !f ! C "( ! ) det (K! M)"0. C

det

(28)

As can be seen from equation (28), the transformed system has m#1 eigenvalues. m eigenvalues are equal to those of the initial system (23) *" with i"1,2, m G G (29)

and the additional eigenvalue is equal to the square of the excitation frequency (in rad/s) * " . K> C (30)

The eigenmodes corresponding to these eigenvalues may now be calculated. As illustrated in equation (26), the eigenmodes are the solution of (K*! * M*)x * " G G For K! * M G 0 !f ! * C G x * "0. G

(31)

*" , an obvious solution of system (31) is G G x x* " G G 0 . (32)

Accordingly, the eigenmodes of the transformed system corresponding to have the G "rst m components equal to the eigenmodes of the initial system. The last component is equal to 0. . With this aim, It remains to evaluate the eigenmode corresponding to , i.e., x * C K> "nding the eigensolutions of the transformed system is also equivalent to "nding the eigensolutions of matrix M*\K*" M\ 0 0 1 K !f 0 C " 0 M\K !M\f C (33)

and M\K !M\f 0 C [x * 2x * x * ] K K> diag( *,2, * ) G K 0 0 * K> 0

"[x * 2x * x * ] K K> M\K 0

(34)

!M\f C

X 0

x* K> :K " X * 0 x K> K>

x* K> :K * x K> K>

0 , C

(35)

JSVI 20013930

10

G. KERSCHEN AND J. C. GOLINVAL

where X"[x 2x ] and "diag( ,2, ) are the eigensolutions of the initial K K system. It follows from equation (35) that M\KX"X , , M\Kx* !M\fx * "x * K> :K K> K> K> :K C 0"0, x* . "x * C K> K> K> K> C (36) (37) (38) (39)

Equation (37) allows us to calculate the "rst m components of the last eigenmode x * : K> x* "[M\K! I]\M\fx* "[K! M]\fx* . K> :K C K> K> C K> K> (40)

[K! M]\ is the dynamic in#uence coe$cient matrix and its spectral expansion is [18] C K x x2 G G [K! M]\" . C ( ! ) C G G G Let us now introduce the spectral expansion (41) in equation (40) x* " K> :K x x2 G G fx * . K> K> ( ! ) C G G G K (42) (41)

Since an eigenmode is de"ned as a scale factor and since x * is a scalar, the "nal K> K> expression for the eigenmode corresponding to is C x* " K> :K x x2 G G f. ( ! ) C G G G K (43)

To summarize, consider a matrix which contains the response of the transformed system (24), i.e., its "rst m rows contain the response of the initial system (23) and its (m#1)th row is the applied force Q*" q(t )2q(t ) L . s(t )2s(t ) L (44)

This matrix has m#1 POMs that have m#1 components. The dominant POM is related to the forced harmonic response of the system and its "rst m components are given by equation (43). Furthermore, if the mass matrix is proportional to density, the "rst m components of the remaining POMs are merely the eigenmodes of the linear system. This perspective should be useful in the context of modal analysis.

6.2.

FORCED HARMONIC RESPONSE

The forced response is de"ned as the part of the response synchronous to the excitation q(t)"q sin t. D C (45)

JSVI 20013930

PROPER ORTHOGONAL MODE DECOMPOSITION

11

The forced response amplitude is described by the following expression [18]: q "" D and q(t)" x x2 G G f sin t. C ! ) ( C G G G K (47) x x2 G G f ( ! ) C G G G K (46)

The displacement matrix Q becomes Q"[q(t )2q(t )] L " x x2 K x x2 G G G G f sin t 2 f sin t . C C L ( ! ) ( ! ) C G C G G G G G K (48)

Equation (48) may be expressed in the form x x2 G G f ( ! ) C G G G K sin sin t 2 C 2 t C

Q"

"q e2 D 1 "[q S ] D 0 $ 0 0 \ 2 0 0 2 0 $ $ $ $ [e R]2

0 2 2 0 0 2 0 #q # #e# D 0 $ 0 q D S #q # D

U

0 0

2 0 0 2 0 $ \ $ $ $ e 2 R #e#

"

q D S #q # D

2 2 0 0 2 0

"

[I ]

GFHFI GHI GFHFI,

V2

e 2 R #e#

(49)

where S is an (m;(m!1)) matrix, I is an (m;n) matrix containing only one non-zero element #q # #e#, and R is an (n;(n!1)) matrix. D Matrices S and R do not in#uence equation (49) since they are both multiplied by zero elements. If S and R are chosen in order that U and V are unitary matrices, then equation (49) is the singular value decomposition of the matrix Q.

JSVI 20013930

12

G. KERSCHEN AND J. C. GOLINVAL

In conclusion, the following points may be noted. 1. Since there is only one non-zero singular value, the forced harmonic response of a linear system is captured by a single POM whatever the number of d.o.f. is. Nevertheless, all the eigenmodes are necessary to reconstruct the response. This property is independent of the mass distribution and underlines the optimality of the POMs described in section 3.1. 2. An analytical expression of the POM is obtained: + K x x 2 /( ! G G G G POM" #+ K x x 2 /( ! G G G G ) C ) C ,f G . ,f# G (50)

Knowing the structural matrices and the spatial discretization of the excitation, the POM may be calculated without "rst simulating the system response as required in the de"nition of the POMs. 3. The expression of the POM (50) is equal, to the norm, to the last eigenmode of the transformed system for the harmonic response (43). This last eigenmode is thus related to the forced harmonic response. 4. The convergence of the dominant POM to an eigenmode is no longer guaranteed. The POM appears now as a combination of all the eigenmodes. However, if the excitation frequency tends to a resonant frequency of the system, for instance, then the C H denominator ! of the jth term of combination (50) tends to zero. It is thus H C observed that this term has a much larger amplitude than the others: + K x x 2 /( G G G POM" #+ K x x 2 /( G G G " x H if P . C H ! G ! G ) C ) C x x 2 /( K H H ,f# #x x 2 /( H H G G ,f ! H ! H ) C ) C f H f# H (51)

Since x 2 f represents a scalar product, the POM has the same direction as the H eigenmode x which means that the POM is equal to the resonating mode shape. This H is consistent with that obtained in reference [10] using the eigensolution perspective. It is worth pointing out that the non-resonating mode shapes should not be revealed by POD.

7. LINEAR NORMAL MODES, NON-LINEAR NORMAL MODES AND PROPER ORTHOGONAL MODES: A GEOMETRIC APPROACH

For the sake of clarity, the eigenmodes of a linear system are called here linear normal modes (LNMs). The determination of LNMs is reduced to the equivalent problem of computing the eigensolutions of linear transformations. Obviously, such an approach as well as the superposition principle is inadmissible for non-linear systems. The concept of synchronous non-linear normal mode (NNM) for discrete conservative oscillators was introduced for non-linear systems by Rosenberg [19]: &&A nonlinear system vibrates in normal modes when all masses execute periodic motions of the same period, when all of them pass through the equilibrium position at the same instant, and when, at any time t, the position of all the masses is uniquely de"ned by the position of any one of them.'' The objective of this section is to examine the geometric interpretation of LNMs, NNMs and POMs.

JSVI 20013930

PROPER ORTHOGONAL MODE DECOMPOSITION

13

Figure 1. Linear system consisting of masses and springs.

7.1.

LINEAR SYSTEMS

Consider a linear system consisting of masses and springs (Figure 1). If the displacement of the ith mass from its equilibrium position is denoted by q , then the equations of motion G of the system are m q "k (q !q )!k (q !q ) where i"1, 2,2, n, K G G G G\ G G> G G> If the co-ordinates are normalized using the transformations equation (52) becomes k K " kG G\ ! G ! G> G> ! G G m m m m m m G G G G\ G> G q "q ,0. L> (52)

"mq , i"1,2, n, G G G

where m "m , R. L>

(53)

The transformed equations of motion (53) may be regarded as those of a unit mass which moves in an n-dimensional space. The right-hand side of equation (53) derives from a potential function K "*; , with ;" L> kG G\ ! G . ! G * 2 m m G G G\ G

(54)

If no external force is present and if the motion is due to an initial displacement, the system occupies at time t"0 a position of maximum potential ;" !; . This latter equation de"nes an ellipsoid which is symmetric with respect to the origin. This ellipsoid is called the bounding ellipsoid because all solutions must lie in this domain. In its de"nition of a normal mode for a linear system, Rosenberg [19] stated that it is a straight line in the ( ,2, ) space which passes through the origin of that space and L which intersects the bounding ellipsoid orthogonally. It follows from the de"nition that the LNMs are the principal axes of the bounding ellipsoid in the ( ,2, ) space. This result L can also be obtained with the interpretation of the eigenvalue problem (section 3.2). Further discussion is given in Appendix B. If the mass matrix is proportional to identity, the LNMs are also the principal axes of the bounding ellipsoid in the (q ,2, q ) space whose expression is L L> k 1 G (q !q )" q2Kq. ;" G\ G 2 2 G (55)

JSVI 20013930

14

G. KERSCHEN AND J. C. GOLINVAL

Figure 2. LNMs and POMs. Principal axes of similar and similarly placed ellipsoids: } }, K; ))))); G; **, POM; , LNM.

As far as the POMs are concerned, they are the principal axes of the ellipsoid c"q2Gq where G is the covariance matrix (cf. Section 3.2). Since for an unforced system with a mass matrix proportional to identity, the POMs and the LNMs coincide, it can be concluded that ; " q2Kq and c"q2Gq are similar and similarly placed ellipsoids. This is illustrated in Figure 2 (two d.o.f. system with an initial displacement).

7.2.

NON-LINEAR SYSTEMS

If an LNM is a straight line in the co-ordinate space, an NNM is represented by a line, straight (similar NNM) or curved (non-similar NNM). But generally, NNMs are non-similar and the POMs, characterized by straight lines in the co-ordinate space, cannot be merged with NNMs. However, due to their optimality property if the motion is a single, synchronous NNM, the resonant POM minimizes the square of the distance with the NNM under the constraint that it passes through the origin of the co-ordinate system and as stated in reference [10], the POM can be considered as the best linear representation of the NNM. It is also worth pointing out that the LNMs are the tangent space to the NNMs [20].

8. CONCLUSION

This paper has presented a new way, based on the singular value decomposition, of interpreting the POMs in the "eld of structural dynamics. This work has underlined some features of POD which might be useful in the future. Since the POMs are related to the vibration eigenmodes in some cases, POD should be an alternative way of modal analysis for extracting the mode shapes of a structure. POMs could also be used to reconstruct a signal using a minimum number of modes.

JSVI 20013930

PROPER ORTHOGONAL MODE DECOMPOSITION

15

ACKNOWLEDGMENTS

Mr Kerschen is supported by a grant from the Belgian National Fund for Scienti"c Research which is gratefully acknowledged. This work presents research results of the Belgian programme on Inter-University Poles of Attraction initiated by the Belgian state, Prime Minister's o$ce, Science Policy Programming. The scienti"c responsibility for this paper is assumed by its authors.

REFERENCES 1. P. HOLMES, J. L. LUMLEY and G. BERKOOZ 1996 ?urbulence, Coherent Structures, Dynamical Systems and Symmetry. Cambridge: New York. 2, W. CAZEMIER 1997 Ph.D. ?hesis, Rijksuniversiteit, Groningen. Proper orthogonal decomposition and low dimensional models for turbulent #ows. 3. G. UYTTERHOEVEN 1999 Ph.D. ?hesis, Katholieke ;niversiteit, ?euven. Wavelets: software and applications. 4. J. P. CUSUMANO, M. T. SHARKADY and B. W. KIMBLE 1993 Aerospace Structures: Nonlinear Dynamics and System Response, American Society of Mechanical Engineers AD-33, 13}22. Spatial coherence measurements of a chaotic #exible-beam impact oscillator. 5. R. KAPPAGANTU and B. F. FEENY 1999 Journal of Sound and <ibration 224, 863}877. An optimal modal reduction of a system with frictional excitation. 6. M. F. A. AZEEZ and A. F. VAKAKIS 1998 ?echnical Report, ;niversity of Illinois at ;rbana Champaign. Proper orthogonal decomposition of a class of vibroimpact oscillations. 7. T. K. HASSELMAN, M. C. ANDERSON and W. G. GAN 1998 Proceedings of the 16th International Modal Analysis Conference, Santa Barbara ;.S.A., 644}651. Principal component analysis for nonlinear model correlation, updating and uncertainty evaluation. 8. V. LENAERTS, G. KERSCHEN and J. C. GOLINVAL 2000 Proceedings of the 18th International Modal Analysis Conference, San Antonio, ;.S.A. Parameter identi"cation of nonlinear mechanical systems using proper orthogonal decomposition. 9. V. LENAERTS, G. KERSCHEN and J. C. GOLINVAL 2001 Mechanical Systems and Signal Processing 15, 31}43. Proper orthogonal decomposition for model updating of non-linear mechanical systems. 10. B. F. FEENY and R. KAPPAGANTU 1998 Journal of Sound and <ibration 211, 607}616. On the physical interpretation of proper orthogonal modes in vibrations. 11. D. KOSAMBI 1943 Journal of Indian Mathematical Society 7, 76}88. Statistics in function space. 12. H. HOTELLING 1933 Journal of Educational Psychology 24, 417}441 and 498}520. Analysis of a complex of statistical variables into principal components. 13. B. RAVINDRA 1999 Journal of Sound and <ibration 219, 189}192. Comments on &&On the physical interpretation of proper orthogonal modes in vibrations''. 14. D. OTTE 1994 Ph.D. ?hesis, Katholieke ;niversiteit, ?euven. Development and evaluation of singular value analysis methodologies for studying multivariate noise and vibration problems. 15. J. STAAR 1982 Ph.D. ?hesis, Katholieke ;niversiteit, ?euven. Concepts for reliable modelling of linear systems with application to on-line identi"cation of multivariate state space descriptions. 16. L. MEIROVITCH 1980 Computational Methods in Structural Dynamics. Alphen a/d Rijn: Sijtho! and Noordho!. 17. B. F. FEENY 1997 Proceedings of ASME Design Engineering ?echnical Conferences, Sacramento, ;.S.A. Interpreting proper orthogonal modes in vibrations. 18. M. GERADIN and D. RIXEN 1994 Mechanical <ibrations, ?heory and Application to Structural Dynamics. Paris: Masson. 19. R. M. ROSENBERG 1962 Journal of Applied Mechanics 29. The normal modes of nonlinear n-degree-of-freedom systems. 20. S. W. SHAW and C. PIERRE 1993 Journal of Sound and <ibration 164, 85}124. Normal modes for non-linear vibratory systems. 21. A. E. BRYSON and Y. C. HO 1975 Applied Optimal Control (Optimization, Estimation and Control). New York: Wiley. 22. D. F. MORRISON 1967 Multivariate Statistical Methods, McGraw-Hill Series in Probability and Statistics. New York: McGraw-Hill.

JSVI 20013930

16

G. KERSCHEN AND J. C. GOLINVAL

APPENDIX A: STATIONARY RANDOM RESPONSE OF A LINEAR SYSTEM TO A WHITE NOISE EXCITATION

This study concerns linear systems subjected to white noise sequences. With this aim, the equation of motion is recast in the state variable from r "Ar#Bw, where A" Z I M\K M\C q"Dr, (A1, A2)

is the system matrix, B is the input matrix, D is the output matrix, and w(t) is a vector white noise process such that E[w(t)]"0 and E[w(t)w( )2]" (t! ). It is assumed that the system is stable and time invariant, and that all processes are Gaussian. In this context, it can be shown [21] that the covariance matrix of the steady state response G "E[r(t)r(t)2] satis"es the Lyapunov equation P AG #G A2#B B2"0. P P (A3)

It is worth pointing out that G also corresponds to the constant in the de"nition of the P controllability grammian W of the system. A If only the displacements are considered, then the covariance matrix of the system response becomes G "E[q(t)q(t)2]"DG D2. O P (A4)

Equation (59) means that the POMs may be evaluated without "rst simulating the system. Indeed, if the structural matrices are assumed to be known, the Lyapunov equation (58) may be solved in order to compute the covariance matrix G and consequently G . The P O POMs are then the eigenvectors of G . The analytical relationship between the POMs and O the eigenmodes is now obscured. If all states (displacement and velocity) are measured, the POMs are merely the eigenvectors of the controllability grammian W . A If the white noise excitation is Gaussian, then the POMs may be geometrically interpreted. In that case, the response of the linear system is also Gaussian and is characterized at each d.o.f. by a probability density function equal to 1 (q)" (2 where "E[q] is the mean and probability density function reads exp \ O\I N

, (A5)

"E[(q! )] is the standard deviation. The joint

1 (q ,2, q )" exp \ K (2 )K 2 K

K G

OG\IG NG

.

(A6)

JSVI 20013930

PROPER ORTHOGONAL MODE DECOMPOSITION

17

It can be demonstrated [22] that the contours of (q ,2, q ) consist of m-dimensional K ellipsoids and that the POMs are the principal axes of these ellipsoids.

APPENDIX B: LNMs AND PRINCIPAL AXES OF THE BOUNDING ELLIPSOID

The LNMs are the eigenvectors of the matrix M\K. In order that the LNMs be the principal axes of the ellipsoid q2M\Kq"1, the matrix M\K must be real, positive de"nite and symmetric [16]. This is the case if the mass matrix is proportional to identity, i.e., M" I. Accordingly, the LNMs are the principal axes of the ellipsoid (1/ )q2Kq"1. This latter expression is to a constant, the expression of the potential energy in the (q ,2, q ) space. Since it is assumed that the mass matrix is proportional to identity, this is L also the expression, to a constant, of the potential energy in the ( ,2, ) space. This is L another way to demonstrate that the LNMs are the principal axes of the bounding ellipsoid in the ( ,2, ) space. L

JSVI 20013930

AUTHOR QUERY FORM HARCOURT PUBLISHERS

JOURNAL TITLE: JSVI ARTICLE NO. : 20013930 Queries and / or remarks Manuscript Page/line 16 Details required

DATE: 22/11/2001

Author's response

Please provide numbers to [19].

page

- The proper orthogonal decomposition in the analysis of turbulent flows
- Modeling and Control of physical processes using proper orthogonal decomposition
- Using the orthogonal array with grey relational analysis to optimize the laser hybrid welding of a 6
- Combinatorial probability interpretation of certain modified orthogonal polynomials
- A new look at proper orthogonal decomposition
- Proper Orthogonal Decomposition and Particle Image Velocimetry
- Wavelet Neural Networks for Adaptive Equalization by Using the Orthogonal Least Square Algorithm
- Abstract Principal Rotation Representations of Proper NxN Orthogonal Matrices
- Some Orthogonal Designs and complex Hadamard matrices by using two Hadamard matrices
- Reduced Order Modeling of the Upper Tropical Pacific Ocean Model Using Proper Orthogonal De

更多相关文章：
**
反螺旋求解_戴建生.pdf
**

*of* *Physical* Sciences and Engineering King’s ...*orthogonal* to *the* *orthogonal* base.7 *The* one-to-...Journal *of* Robotic Systems2003 *Using* *the* ...**
2019-Discriminant Function Analysis - California St....ppt
**

*The* *interpretation* *of* discrim results are always ...*the* cell will become *singular* and cannot be ...interpreted *using* *the* loading matrix (later) ...**
Bayesian linear regression.pdf
**

Note that if XX is*singular* then we must ...*the* component *of* Y which is *orthogonal* to *the* ...Example *of* *using* *the* evidence to select model order...**
Adaptive Manifold Learning.pdf
**

illustrate its failure*modes* *using* two simple ...*orthogonal* projection *of* each xij to *the* tangent ...(k) (k ) *the* d largest *singular* vectors *of* ...**
Discriminant Function Analysis - California State U....ppt
**

*The* *interpretation* *of* discrim results are always ...*the* cell will become *singular* and cannot be ...interpreted *using* *the* loading matrix (later) ...**
Cost.pdf
**

A is N N , non*singular*, large, and sparse. ...We can make *the* residuals *orthogonal* by combining...*The* *interpretation* *of* *the* solution *of* a linear ...**
Robust HInfinity Filtering for A Class ***of* Uncertain.pdf

*the* robust stability and a Control *of* *singular* ...better describe *physical* systems than regular ones....then there exists two *orthogonal* matrices L1 and ...**
On Principal Angles Between Subspaces in R^n.pdf
**

[w” in terms*of* their *orthogonal* and oblique ...[3] used *the* *singular* value decomposition *the* ...*interpretation* submatrices is *the* following COROLLARY...**
Power Allocation Based on Fast Water-Filling for.pdf
**

In particular, OFDM (*Orthogonal* Frequency Division ...[8-9] *using* SVD (*Singular* Value Decomposition)....some circuit power to support *the* *proper* ...**
Multilinear Image Analysis for Facial Recognition.pdf
**

*of* PCA and *the* *singular* value decomposition (SVD...different *modes* underlying *the* formation *of* facial ...I t *orthogonal*izes these N spaces and decomposes...**
1978-More-***The* Levenberg-Marquardt algorithm-impleme....pdf

*the* main are *the* *proper* *use* *of* implicitly scaled...This produces an *orthogonal* matrix Q and a ...*the* *singular* value decomposition 2 2 ~? ° i ...**
Intorduction to HHT_图文.ppt
**

*Orthogonal* Functions aka *Singular* Spectral Analysis ...not to force *the* data to fit preconceived *modes*...*Using* *the* Empirical *Mode* Decomposition Method, Proc...**
...Partial Differential Equation Toolbox for ***use* wi....pdf

*The* MathWorks, Inc, Natick MA, 2006. [2] R....Part I: Branches *of* non*singular* solutions. ...*physical* processes *using* *proper* *orthogonal* ...**
CS 223-B L9 Optical Flow_图文.ppt
**

Rank (E) = 2, two equal nonzero*singular* ...*the* 2D projection *of* *the* *physical* movement *of* ...Warp I(t-1) towards I(t) *using* *the* estimated...**
Let A = PDP Gamma1.pdf
**

I + (D ? k I)?1 B must be*singular*. Therefore...A geometric *interpretation* *of* *the* bound on approximate...and close to zero if they are nearly *ort*...**
Discriminant Function Analysis - California State U....ppt
**

*The* *interpretation* *of* discrim results are always ...*the* cell will become *singular* and cannot be ...interpreted *using* *the* loading matrix (later) ...**
***The* Higher-Order *Singular* Value Decomposition --*The*....pdf

(3-array) Tijk , where*the* three *modes* ...*singular* values *of* A (*the* square roots *of* *the* ...augmented with columns to square and *orthogonal* (...**
Ellipse Detection Based on Principal Component Anal....pdf
**

*the* coordinate frame *of* principal axis *using* *the* ...*interpretation* *of* medical image, detection *of* ...,N) by *singular* value decomposition *of* ...**
Level spacing distribution ***of* pseudointegrable bill....pdf

prediction*of* GOE (Gaussian *Orthogonal* Ensemble). ...*the* numerical treatment lacking any *physical* ...*The* special feature in case *of* *the* *singular* ...**
A 70GOPS, 34mW Multi-Carrier MIMO Chip in 3.5mm 2.pdf
**

*The* chip implements a 4×4 adaptive *singular* ...*using* *the* architectural/circuit techniques that ...6. Measured tracking *of* eigen-*modes*. y3(k) s...

Note that if XX is

illustrate its failure

A is N N , non

[w” in terms

In particular, OFDM (

Rank (E) = 2, two equal nonzero

I + (D ? k I)?1 B must be

(3-array) Tijk , where

prediction