Dynamic substructuring of damped structures using singular value decomposition Roger Ohayon, R. Sampaio, Christian Soize
To cite this version:
Roger Ohayon, R. Sampaio, Christian Soize. Dynamic substructuring of damped structures using sin- gular value decomposition. Journal of Applied Mechanics, American Society of Mechanical Engineers, 1997, 64 (2), pp.292-298. 10.1115/1.2787306. hal-00770023
HAL Id: hal-00770023 https://hal-upec-upem.archives-ouvertes.fr/hal-00770023 Submitted on 4 Jan 2013
HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Dynamic Substructuring of Damped Structures Using Singular Value Decomposition
R. Ohayon Chair of Mechanics, CNAM, 2 rue Conte,´ F-75003 Paris, France (mem. ASME)
R. Sampaio Department of Mechanical Engineering, Pontificia Universidade Catolica´ do Rio de Janeiro, 22453-900 Rio de Janeiro - RJ - Brazil
C. Soize Structures Department, ONERA, BP 72, F-92322 Chatillon Cedex, France Abstract
This paper deals with the theoretical aspects concerning linear elastodynamic of damped continuum medium in the frequency domain. Eigenvalue analysis and frequency response function are studied. The methods discussed here use a dynamic substructuring approach. The first method is based on a mixed variational formulation in which Lagrange multipliers are introduced to impose the linear constraints on the coupling interfaces. A modal reduction of each substructure is obtained using its free-interface modes. A practical construction of a unique solution is carried out using the Singular Value Decomposition (SVD) related only to the frequency-independent Lagrange multiplier terms. The second method is similar to the first one replacing the free-interface modes by the fixed-interface modes and elastostatic operator on the interface of each substructure.
1. Introduction field. As a consequence, the final system for the mixed formulation has a rank-deficiency in the matrix that describes the constraints. This leads In this paper, we are interested in eigenvalue and frequency response to non-uniqueness of the solution. In order to avoid this difficulty, a new function calculations of a linear dynamic three-dimensional bounded constructive approach is proposed consisting in using a Singular Value damped elastic structures subjected to prescribed forces. Recall that the Decomposition (SVD) of the frequency-independent constraint matrix frequency response functions allow deterministic and stationary random and chose a "least-square" solution that is in fact the solution of the orig- analyses to be performed (Kree and Soize, 1986; Argyris and Meljnek, inal problem. Due to a relatively small number of degrees of freedom 1991). More precisely, this paper is devoted to theoretical aspects of in the reduced model, the use of SVD is particularly efficient. Since the structure-structure coupling by dynamic substructuring methods using problem under consideration is linear, SVD is used only once. Conse- modal reduction procedures. The proposed methodology can be ap- quently, the SVD appears as an efficient and reliable tool to solve this plied to general linear coupled systems such as fluid-structure interaction rank-deficiency problem. It should be noted that SVD has been used for problems (Morand and Ohayon, 1995; Soize, Desanti and David, 1992). undamped linear vibration analysis of plates using dynamic substructur- For linear structural vibrations, dynamic substructuring techniques based ing by analytical methods (Jen, Johnson and Dubois, 1995). Let us recall on the use of the fixed-interface modes or free-interface modes (com- that SVD has also been used in the area of the nonlinear dynamical anal- pleted by static boundary functions, attachment modes, residual flexibil- ysis of multibody systems with nonlinear constraints (Singh and Likins, ity, etc.) of each substructure have been widely developed in the litter- 1985; Shabana , 1991; Schmidt and Muller,¨ 1993). ature: for conservative structures see for example (Hurty, 1965; Craig Now we give a short description of the content of each section. and Bampton, 1968; MacNeal, 1971; Rubin, 1975; Flashner, 1986; Min, Igusa and Achenbach, 1992; Farhat and Geradin, 1994) and for damped Section 2 deals with the displacement and mixed variational formulations structures (Klein and Dowell, 1974; Hale and Meirovitch, 1980; Leung, for the coupled linear structure-structure problem, Lagrange multiplier 1993; Farstad and Singh, 1995; Rook and Singh, 1995). field being introduced in the mixed problem. In Section 3, we present a dynamic substructuring method using the Some papers are based on a mixed formulation using a Lagrange multi- free-interface modes of each linear substructure. The modal reduction plier in order to impose the linear constraints on the coupling interfaces procedure is carried out using a new explicit construction of the La- (see Klein and Dowell, 1974; Min, Igusa and Achenbach, 1992; Farstad grange multiplier admissible space. Two practical constructions of the and Singh, 1995; Rook and Singh, 1995). Within the context of finite frequency response function of the global linear damped structure and element discretization of linear structural dynamic problems, Farhat and the eigenvalues of the associated conservative structure are performed Geradin (1994) have also introduced a Lagrange multiplier to take into using SVD once on a part of the linear system to be solved, namely on account incompatible meshes on the interface (their analysis is devoted the frequency-independent Lagrange multiplier terms. to undamped structures using a component mode method based on fixed- interface modes and static boundary functions). Section 4 is devoted to a dynamic substructuring method using the classi- cal Craig and Bampton fixed-interface modes and boundary static func- Below, we present an original general approach for damped structures tions of each linear substructure, presented in an original general frame- using continuum-based variational formulations and Ritz-Galerkin pro- work allowing various other decomposition procedures to be obtained. jection methods using free-interface modes and fixed-interface modes of After having constructed the reduced matrix model of each substructure, each substructure (in this paper we do not consider mathematical aspects we explain two procedures for the assemblage of the substructures and of error estimates connected to the truncation of the modal series). For the construction of a solution, (1) in a classical manner and (2) as in this purpose, various rigorous algebraic decompositions of admissible Section 3 using Lagrange multiplier field and SVD. classes of the unknown fields are introduced and leads to several linear Finally, in Section 5, some conclusions are presented. dynamic substructuring methods, the continuity of the displacement field on the interface being imposed through the use of a Lagrange multiplier
Journal of Applied Mechanics 1 Ohayon, Sampaio, Soize 2. Displacement and Mixed Variational For- sesquilinear forms on CΩr ×CΩr corresponding to the mass, stiffness and mulations for the Coupled Structure-Structure damping operators of substructure Ωr, are introduced as follows Problem mr(ur δur)= ρr ur δur dx (3) 2.1 General Mechanical Hypotheses Ωr In this section, the following hypotheses are introduced: r r r r r r - One considers the linear vibrations of a three-dimensional structure k (u δu )= σij (u ) εij ( δu ) dx (4) Ω about a static equilibrium configuration which is considered here as a r natural state (for the sake of brevity, prestress are not considered but r r r r r r d (u δu )= sij (u ) εij ( δu ) dx (5) could be added without changing the theory). Ωr - The structure is only submitted to prescribed external forces (no pre- It should be noted that the hermitian form mr is positive definite on scribed displacement). r r CΩr ×CΩr . The hermitian forms k and d are semi-definite positive With the above hypotheses, there are two cases. (degenerated forms) since rigid body displacement fields are allowed in (1)- The first one, which is the only case considered in this paper, cor- r 3 the present case. The set Rrig of Ê -valued rigid body displacement fields responds to prescribed external forces which are in equilibrium at each r (of dimension 6) is a subset of CΩr . Consequently, for all δu in CΩr , instant. Consequently, the displacement field of the structure is defined r r r r r r r r k (u δu ) and d (u δu ) are equal to zero for any u in Rrig. up to an additive rigid body displacement field. In this case, we are r We then define the following sesquilinear form z on CΩr ×CΩr only interested in the part of the displacement field due to the structural deformation. We will see below how the rigid body displacement field zr(ur δur)= −ω2 mr(ur δur)+iω dr(ur δur)+kr(ur δur) (6) can be disregarded. (2)- The second case corresponds to prescribed external forces which are Finally, we define fr by the relation not in equilibrium at some instants. To solve this problem, the method consists in transforming this case to the first case by adding an additional r r r r ≪f δu ≫= gΩr δu dx + gΓr δu ds (7) external force related to rigid body field. For the sake of brevity, this case Ωr Γr will not be considered in the present paper. 2.3 Continuum-Based Variational Formulations for Two Coupled One presents a variational formulation of the problem (first case), taking Substructures Ω1 and Ω2 into account an additional small structural damping based on a linear viscoelastic model with an instantaneous memory. A frequency do- We consider a structure composed of two substructures Ω1 and Ω2 that main formulation is used, the convention for the Fourier transform being interact through a common boundary Γ (the extension to the case of more u(ω)= e−iωt u(t) dt where ω denotes the circular frequency, u(ω) is than two substructures is straightforward). The notations introduced in
Ê 3
Ê C a vector in C and u(ω) its conjugate ( and denote the set of real and Section 2.2 are used with r =1 and r =2. The linear coupling conditions complex numbers respectively). on Γ are written as u1 = u2 on Γ (8) 2.2 Notation for a Substructure Ω r 1 1 2 2 σtot n = −σtot n on Γ (9) We consider a structure formed by substructures that will be denoted where nr is the unit normal to , external to r. by an index r. Let Ωr be the 3D-bounded domain occupied at static Γ Ω 1 2 equilibrium by the substructure labelled by index r. Let ∂Ωr = Γr ∪ Γ 2.3.1 Basic (u u ) Variational Formulation P0 with be the boundary of (assumed to be smooth). The Γr ∩ Γ = ∅ Ωr 1 2 1 2 boundary Γ will be the interaction surface with another substructure. The For all real ω in Ê and prescribed (f f ), find (u u ) in CΩ1 ×CΩ2 verifying the linear constraint u1 = u2 on Γ, such that, for all (δu1 δu2) external prescribed volumetric and surface force fields applied to Ωr and 1 2 r r r r in CΩ1 ×CΩ2 verifying the linear constraint δu = δu on Γ, one has Γr are denoted by gΩr and gΓr respectively. Let u =(u1 u2 u3) be the displacement field at each point x =(x1 x2 x3) in cartesian coordinates. 1 1 1 2 2 2 1 1 2 2
3 3 z (u δu )+ z (u δu )=≪f δu ≫ + ≪f δu ≫ (10) Ê The set of admissible displacement fields with values in C (resp. in ) is denoted by CΩr (resp. RΩr ) and is used for dissipative problems (resp. From the mathematical point of view (see Dautray and Lions, 1992), by associated conservative problems). For substructure Ωr, one denotes the 1 r 3 taking Sobolev space H C as admissible space , the existence r r (Ω ) CΩr test function (weighted function) associated with u as δu ∈CΩ (or in r and uniqueness of a solution of P0 can be proved. RΩr ) . The strain tensor is defined by 1 2
2.3.2 Mixed (u u Ð) Variational Formulation P1 1 ε (ur)= (ur + ur ) (1) This formulation consists in relaxing the linear constraint (defined by ij 2 i j j i Eq. (8)) used in P0 by the introduction of a Lagrange multiplier field Ð defined on Γ. Let ΛΓ be the admissible set of Lagrange multiplier fields in which v j denotes the partial derivative of v with respect to xj. The 3
defined on Γ with values in C . total stress tensor is defined by 1 2 1 2
Formulation P1 . For all real ω in Ê and prescribed (f f ), find (u u ) 1 2 r r r in CΩ1 ×CΩ2 and Ð in ΛΓ such that, for all (δu δu ) in CΩ1 ×CΩ2 and σtot = σ + iωs (2)
for all δ Ð in ΛΓ, one has r r r r
where σ is the elastic stress tensor defined by σ (u )= aijkh εkh(u ) 1 1 1 2 2 2 1 2 1 2 1 1 Ð ij z (u δu )+z (u δu )+b(Ð δu −δu )+b(δ u −u )=≪f δu ≫ + ≪f r r r and iωs is the viscous part of the total stress tensor such that sij(u )= (11) r bijkh εkh(u ) (using summation over repeated indices). The mechanical where b( ) is defined by coefficients aijkh and bijkh are independent of ω and verify the usual properties of symmetry and positivity (see Marsden and Hughes, 1983).
r r Ð r b(Ð u )= u ds (12) The mass density is denoted by ρ . For the dissipative problem, three Γ
Journal of Applied Mechanics 2 Ohayon, Sampaio, Soize Space of Traces on Γ. The set of the traces related to the boundary Γ, The present approach is based on the fact that any Ð in ΛΓ can be r r is denoted by CΓ. Therefore, if u ∈ CΩr , then the trace of u on Γ is expanded on a complete orthonormal set in CΓ and consequently, the r N denoted by u|Γ and belongs to CΓ. In Eq. (11), ΛΓ is the dual space of projection of the Lagrange multiplier Ð is done on the subspace WΓ N CΓ. of CΓ ⊂ ΛΓ. A characterization of WΓ requires the construction of a N basis of W denoted by {w1 w }. One possible method consists Remark. From the mathematical point of view (see Dautray and Li- Γ N
1 r 3 1 2 3 in extracting an independent system of N functions from the family C ons, 1992), by taking CΩr = H (Ω C )), CΓ = H (Γ ) and 1 1 2 2 N −1 2 3 {u1|Γ uN1 Γ u1|Γ uN2 Γ}. Consequently, for all Ð in WΓ , ΛΓ = H (Γ C ), the existence and uniqueness of a solution of | | one has formulation P1 can be proved using the so called LBB condition related N to the sesquilinear form b (see Brezzi and Fortin, 1991). It should be Ð = p w (18)
1 2 3 −1 2 3 γ γ C noted that H C is dense in H . (Γ ) (Γ ) γ =1
red 3. Dynamic Substructuring Using the Free-Inter- The Reduced Problem P1 . We use the Ritz-Galerkin method consist- face Modes of Each Substructure ing in substituting Eqs. (16) and (18) into Eq. (11). Using the orthogonal- ity conditions defined by Eqs. (14) and (15) and introducing the vectors The method is based on the use of the mixed variational formulation 1 1 1 2 2 2 of generalized coordinates q = (q1 qN1 ), q = (q1 qN2 ) and defined by P1. Then, a modal reduction is carried out using the Ritz- p = (p1 pN ), one deduces the following finite-dimension reduced Galerkin projection on the free-interface modes of each substructure. problem from P1 Finally, the Singular Value Decomposition (SVD) is used for the con- struction of the solution. 1 T 1 1 Z (ω) 0 B1 q F 2 T 2 F 2 3.1 Free-Interface Modes of a Substructure Ωr 0 Z (ω) B2 q = (19) B1 B2 0 p 0 A free-interface mode of a substructure Ωr (for r=1 or r=2) is defined as an eigenmode of the conservative problem associated with the substructure r 2 in which, for all real ω and for r = 1 and r = 2, [Z (ω)] is an (Nr × Ωr, subject to zero forces on Γ. The real eigenvalues ω ≥ 0 and the r Nr) complex symmetric matrix, [Br] a (N × Nr) real matrix which is eigenmodes u in RΩr are solutions of the following spectral problem: r Nr r independent of ω and F a C -valued vector. Matrix [Z (ω)] is defined find ω2 , ur ur 0 such that for all δur , one has ≥ 0 ∈RΩr ( = ) ∈RΩr by r 2 r r r 2 ω ω iω kr(ur δur)= ω mr(ur δur) (13) [Z ( )] = − [M ]+ [ D ]+[ K ] (20) r r r 2 where [M ] and [ K ] are diagonal positive-definite matrices such that It can be shown that there exist six zero eigenvalues 0=(ω−5) = [Mr] = r δ and [ Kr] = r ωr 2 δ , [ Dr] is a full symmetric ωr 2 (associated with the rigid body displacement fields) and αβ α αβ αβ α α αβ = ( 0) positive-definite matrix, such that [ Dr] = dr(ur ur ). Consequently, that the strictly positive eigenvalues (associated with the displacement αβ β α for all real ω, matrix [Zr(ω)] is invertible. Matrix [B ] is such that for field due to structural deformation) constitute the increasing sequence r r 2 r 2 r r all α in {1 Nr} and γ in {1 N}, one has 0 < (ω1) ≤ (ω2) . The six eigenvectors {u−5 u0} associated with zero eigenvalues span Rrig (space of the rigid body displacement r r r r [Br]γα = b(wγ uα) (21) fields). The family {u−5 u0; u1 } of all the eigenvectors forms a complete set in RΩ . For α and β in {−5 0; 1 }, we have the r Finally, vector F r is such that, for all α in {1 N }, one has orthogonality conditions r r r r r r r r Fα =≪f uα ≫ (22) m (uα uβ)= δαβ α (14)
r r r r r 2 k (u u )= δαβ ω (15) α β α α 3.3 Practical Construction of the Frequency Response Function of r red in which α > 0 is the generalized mass of mode α depending on the the Global Structure Using Reduced Problem P1 and SVD normalization of the eigenmodes. First, we introduce the (N × M) real matrix [ B ] such that 3.2 Modal Reduction of P1
Nr M = N1 +N2 [ B ]=[ B1 B2] (23) We introduce the subspace CΩr of CΩr , of dimension Nr, spanned by r r r Nr {u1 u } with Nr ≥ 1. For all u in CΩ , one has Nr r and write Eq. (19) as
Nr ur qr ur Z(ω) BT q F = α α (16) = (24) α =1 B 0 p 0
r in which qα are complex-valued generalized coordinates. Concerning the In order to solve Eq. (24), we use a Singular Value Decomposition trace of the displacement field (including rigid body displacement field) (SVD) of [ B ]. It is know that there exist algorithms (see Golub and Van on , the subspace spanned by the family ur ur ur is Γ { −5|Γ 0|Γ; 1|Γ } Loan, 1989) which are very efficient for the construction of the SVD a complete set in CΓ (for the two domains r =1 and r =2). Consequently, of reasonable size matrices. This is the case for the reduced problems r the family {u1|Γ } forms a complete set of the displacement field on obtained by modal projection as Eq. (24). In the proposed approach, Nr Γ due only to the structural deformation. Let CΓ be the subspace of CΓ it should be noted that SVD will only be applied to the submatrix [ B ] spanned by the finite family {ur ur }. Let WN be the subspace in Eq. (24). The SVD of (N × M) real matrix [ B ] with M ≥ N (see 1|Γ Nr |Γ Γ of CΓ of finite dimension N ≤ N1 + N2 defined by Section 3.2) consists in constructing the following decomposition
N N1 N2 T WΓ = CΓ ∪CΓ (17) [ B ]=[U ][Σ][V ] (25)
Journal of Applied Mechanics 3 Ohayon, Sampaio, Soize −1 k where [U ] is an (N × N) orthogonal real matrix, [V ] is an (M × M) ek =< [ Z(ω)] F V >. Then, the projection of Eq. (33) on the orthogonal real matrix and [Σ] is a (N× M) real matrix which is written remaining {Vn+1 VM } yields for all k in {n+1 M}, in block form as [Σ]= Σ+ 0 (26) n −1 k′ k −1 k ξ = − y ′ <[ Z(ω)] V V > + <[ Z(ω)] F V > + k k in which [0] is the (N× (M − N)) null matrix and [ Σ ] is the (N× N) k ′=1 diagonal matrix of positive or null singular values σk such that σ1 ≥ (35) σ2 ≥ ≥ σN ≥ 0. Let n be the integer such that 1 ≤ n ≤ N such that The corresponding algorithm is summarized below. Step 0: calculating the SVD of [ B ] in order to obtain its rank n and σ1 ≥ σ2 ≥ ≥ σn > σn+1 = = σN = 0 (27) V1 VM . Consequently, the rank of [ B ] is equal to n and Eq. (25) yields the SVD Then, for each real ω, expansion Step 1: solving the linear equation of dimension n with n+1 right-hand n side members {F; V1 Vn} k kT [ B ]= σk U V (28) k =1 [ Z(ω)] X0 = F ; [ Z(ω)] Xk = Vk k ∈ {1 n} ; (36) in which the vectors Uk and Vk are the columns of [U ] and [V ] and such that Step 2: constructing (n× n) complex symmetric matrix [ E(ω)] such j k j k k k′ ′