On the relevance of inertia related terms in the equations of motion of a ﬂexible body in the ﬂoating frame of reference formulation

The ﬂoating frame of reference formulation is an established method for the description of linear elastic bodies within multibody dynamics. An exact derivation leads to rather complex equations of motion. In order to reduce the computational burden, it is common to neglect certain terms. In the literature this is done by strict application of the small deformation assumption to the kinetic energy. This leads to a remarkably simpliﬁed set of equations. In this work, the signiﬁcance of all terms is investigated at the level of the equations of motion. It is shown that for a large number of applications the previously mentioned set of simple equations is sufﬁcient. Furthermore, scenarios are described in which this simple set is no longer accurate enough. Finally, guidelines are provided, so that engineers can decide which terms should be considered or not. The theoretical conclusions drawn in this work are underlined by qualitative numerical investigations.


Introduction
The floating frame of reference formulation (FFRF) is a widely used strategy for the inclusion of a flexible body into multibody system dynamics. The key idea is the separation of the overall motion into nonlinear rigid body motion superimposed by a small elastic deformation. The latter one is described by the superposition of weighted shape functions which are expressed in a body attached coordinate system which translates and rotates with the body. The final equations of motion for each flexible body contain both simple and complex terms. Hence, the numerical evaluation of some terms is cheap, yet expensive for the others. The numerical costly terms are associated with the body's flexibility, where the numerical effort increases with the number of considered shape functions. For many practical applications it has been observed that the flexibility related terms are of minor significance. Those terms are B W. Witteveen wolfgang.witteveen@fh-wels.at 1 Upper Austria University of Applied Sciences, Wels, Austria the inertia of the deformed body, inertia coupling and parts in the centrifugal and the Coriolis forces, which are related to elastic deformations. In the literature two approaches can be found that deal with the question of which terms need to be considered and which need not. The first one is the strict application of a small deformation assumption at the level of the kinetic energy; see, for example, [12] or [21]. With this assumption, a remarkably simple set of equations is obtained. This work demonstrates why the latter description is probably sufficient for many applications, but not for all. The second approach, examples of which can be found in the software package MSC.ADAMS [11], gives the user the possibility to activate or deactivate certain terms in the equations of motion. The final decision is up to the user, but the guidance is not very clear. However, for both approaches there is a final uncertainty which is addressed in this work.
All inertia related terms of a flexible body in the FFRF are investigated on the level of the equations of motion with respect to their relevance. It is shown why the simplest possible formulation, stemming from the strict application of the small deformation assumption at the level of the kinetic energy, is sufficient for many applications. Beyond that, a clear guidance will be given when, in addition, the deformation dependent inertia tensor, inertia coupling, or elastic deformation related parts of the centrifugal and Coriolis forces should be considered. This will be finally summarized in a "set of guidelines".
The paper is organized as follows: The first section briefly recaps the equations of motion of a free flexible body considering the FFRF. The assumptions, on which this paper is based, together with the resulting implications for the equations of motion, are documented in the next section. The Finite Element structures which are used for the sake of numerical illustration of the theoretical conclusions are introduced in the subsequent section. Following this, the small deformation assumption and its implications for the magnitude of the modal coordinates are discussed. In the next section, the maximum entries of significant matrices are estimated. The latter two estimations will be important for the determination of the significance of certain terms. The subsequent section is the key section of this contribution. All terms will be investigated with respect to their significance. It starts with the inertia of the deformable body followed by the inertia coupling of the rotational and flexible degrees of freedom. Finally, the terms of the quadratic velocity vector will be discussed. All the theoretical considerations in this chapter are accompanied by qualitative numerical investigations. In the subsequent chapter, a set of guidelines is presented in order to simplify the decision when and which term has to be considered or not. The paper will be concluded by summarizing the benefits when certain terms of the equations of motion can be neglected.

Brief review of the equations of motion
In order to avoid redundancy, an extensive derivation of the equations of motion is omitted and only relevant equations are reviewed. Details can be found in Shabana [1] and the work of Sherif and Nachbagauer [2]. For better readability with respect to the mentioned publications, the same notation is used in this work. In contrast to [1], a superscript referring to the number of the flexible body will be omitted.
For the description of the position vector r, pointing to an arbitrary point P on a deformable body, two coordinate systems are used (see Fig. 1). One is a space-fixed coordinate system, which is called the inertial frame (IF), and the second is a so-called floating frame (FF), which is a body attached coordinate system.
The vector r can be given as where R is the (3 × 1) position vector from the origin of the IF to the origin of the FF. The matrix A is the orthogonal rotation matrix transferring vectors expressed in the FF to their counterparts expressed in the IF. The vector from the origin of the FF to the point P is represented byū. The overbar denotes that the quantity is expressed with respect to the FF. The vectorū can be subdivided intoū 0 andū f whereū 0 holds the time invariant vector from the origin of the FF to the particle's position in the undeformed configuration. The (3 × 1) vectorū f represents the elastic deformations and is expressed by a superposition of weighted shape functions in the formū f = Sq f . The (3 × n f ) matrix S contains in its columns the time invariant shape functions, and the (n f × 1) vector q f holds the time variant weighting factors. In the context of the Finite Element Method, shape vectors are used instead of shape functions. Those shape vectors are commonly called "modes". The expressions "shape functions, "shape vectors" and "modes" are used in this work synonymously. The kinetic energy T of the deformable body is given as with the generalized coordinates The (n Θ × 1) vector Θ contains the generalized rotational coordinates where the number n Θ depends on the particular choice (Euler parameters, Euler angles, . . . ). All coordinates are collected in the (n × 1) vector q with n = 3 + n Θ + n f . According to [1, Chaps. 5.1 and 5.2], the submatrices can be computed using where the matrix I holds a (3 × 3) identity matrix andũ holds the (3 × 3) skew symmetric matrix representation of the vectorū. The (3 × n Θ ) matrixḠ transforms the generalized velocities into the (3 × 1) angular velocity vectorω, The parameters m and ρ represent the body's mass and density, respectively, and V is the volume of the body. The application of the Lagrange formalism, together with the consideration of the body's elastic forces, leads to the equations of motion in the form of where the (n × n) time invariant matrix K can be given as (see [1,Eq. (5.145)]). In contrast to the referenced equation, the constraints for a nonminimal set of rotational coordinates (e.g., Euler parameters) are omitted here, since they are not relevant for the considerations in this work. According to [2,Eqs. (25), (54) and (66)], the so-called (n × 1) quadratic velocity vector Q v can be written as Note that the different formulas for the quadratic velocity vector given in [2] and [1] are equal. This is discussed in detail in [2].

Assumptions and simplifications
In this section, some very common assumptions are discussed, since they lead to significant simplifications in the mass and stiffness matrix, as well as in the quadratic velocity vector.

Euler parameters
Euler parameters are used in this work for the parametrization of the rotation, see [1,Chap. 2.1]. Therefore, n Θ = 4 and the constraint Θ T Θ = 1 has to be considered in the final set of equations. An interesting discussion on Euler parameters can be found in [3]. However, due to its irrelevance for the subject of this paper, the constraint equation and its effect on the equations of motion are not considered. A consequence of the use of Euler parameters is the identityĠΘ = 0 (9) (see [1,Eqs. (2.69) and (2.87)]). This leads to simplifications in (8) since terms containing the productĠΘ evaluate to zero.

Linearized mean-axis conditions
As mention before, the position of a mass particle is described with respect to an FF attached to the body. A distinct definition of the FF requires the determination of its origin and orientation. The most relevant suggestions, together with interesting discussions on this issue, can be found, among other references, in [4][5][6] and [7]. However, in this work a linearized mean axis frame, which is sometimes referred to as Buckens or linearized Tisserand frame, is considered. The constraint equations for the origin and the orientation of the FF are In the next subsection, an advantageous choice of the shape functions in S is discussed. Those shape functions automatically fulfill the constraints (10), and hence they need not be considered explicitly in the final set of equations.

Use of (pseudo) free surface modes together with an FF origin fixed to the center of gravity of the undeformed body
Again, just a brief review is given on the implications of the choice of free surface modes and an FF origin which is initially attached to the center of gravity of the undeformed body. More details can be found in [5] and [8].
Attaching the origin of the FF to the center of gravity of the undeformed body implies Free surface modes are the eigenvectors of the unconstraint body with nonzero eigenvalues. The zero eigenvalue modes belong to the rigid body modes, see [8]. The properties of free surface modes can be given as which is explained in detail in [8]. Due to the fact thatū f = Sq f , and thereforeu f = Sq f , the constraints (10) are fulfilled automatically, when free surface modes are used. The application of (12) to the equations of motion leads to significant simplifications of the equations of motion. Details on all resulting simplifications can be found in [5] and [8].

Use of mass normalized modes
Free surface modes are stiffness and mass orthogonal. Mass normalized modes are characterized by a scaling so that where I nf is an (n f × n f ) identity matrix. The use of mass normalized modes leads to a diagonal stiffness matrix K ff , which is hereafter denoted as Ω 2 . Its ith main diagonal entry Ω 2 i holds the eigenvalue of mode number i, which is equal to the squared eigenfrequency of the mode. From now on, it is assumed that the mode numbering correlates with increasing eigenfrequencies.

Use of central principal axis of inertia
It is assumed that the initial orientation of the axis of the FF matches the principal axis of inertia of the undeformed body.

Final equations of motion
The application of the former assumptions to the mass matrix (2) and the quadratic velocity vector (8) lead to much simpler equations of motions of the free flexible body in the form of (see [1,Eqs. (8.15) and (8.23)] and [2, Eqs. (54) and (66), together with the subsequent discussion]). The slight difference to the referenced equations (ũ f instead ofũ) in the third term of (Q v ) Θ stems from a simplification due to (12). The same simplification leads toũ f instead ofũ in m Θf , which is well documented in [5] and [8]. In Shabana [1], so-called "inertia shape" integrals can be found (see [1,Eqs. (5.70), (5.155), (5.156) and (5.158)]). The scalarū 0k holds the kth component of the vectorū 0 (ū T 0 = ū 01ū02ū03 ). The (1 × n f ) row vector S k holds the kth row of S. In case free surface modes are used, it can be easily shown by the right equation of (12) thatĪ kl is equal toĪ lk . This interesting relationship has been found in unpublished records of Sherif ([2, 3, 8, 18] and [19]). By the use of (15), the integrals in Eq. (14) can be evaluated, and the equations of motion can be given in the form where the (3 × 3) matrix I 0 , the (3 × 3n f ) matrix W 1 , the (3n f × 3n f ) matrix W 2 and the (3n f × n f ) matrix W 3 hold invariant quantities and are defined as and the (3n f × 3) matrix Q f holds the modal coordinates in the form of This convenient nomenclature for the matrix W 1 and W 2 was introduced by Sherif, see [18] and [19]. The matrix W 3 can be found in the work of Shabana [1, Eqs. (5.74) to (5.76)] but is not denoted as W 3 there. It is worth mentioning that the (3 × 3) matrix I 0 represents the inertia tensor of the undeformed body. Note that (14) and (16) contain the same equations of motion, the notation is merely different. Both equations are used in the following investigations.

Comment on the use of different mode basis
The literature offers several possibilities for the particular choice of shape vectors. The most famous mode base is probably the one of Craig [23]. Comparative studies of other possibilities can be found, among other references, in [24] and [25]. However, at this point of this work it is interesting to note that arbitrary mode shapes can be transformed into a mode base possessing the properties given before, see [8].

Comments on the complexity of the equations of motion and goal of this work
An obvious disadvantage of the equations of motion (16) are the necessary matrix vector operations related to the invariant matrices W 1 , W 2 and W 3 stemming from the body's flexibility. This increased numerical effort affects the residuum of (16) and the computation of the Jacobian matrix. In order to minimize this computational burden, some suggestions for the negligence of certain terms can be found in the literature and multibody dynamic software packages. As an example, [12] can be cited. This paper refers to the classical paper [13] where the hypothesis of small elastic deformations around the undeformed configuration is applied on the level of kinetic energy. As a consequence, the coupling between the rotational and flexible degrees of freedom in the mass matrix and in the quadratic velocity vector vanishes in the final equations of motion. Another example is the commercially available software package MSC.ADAMS [14], which offers the user a possibility to deactivate certain terms of the mass matrix. The explanation of the terms is quite vague and, hence, the decision is difficult to make for a user without extensive knowledge of the theory of flexible multibody systems. The goal of this work is a systematic investigation on the relevance of all the terms related to the invariant matrices on the level of the equations of motion. As far as we know, this kind of in depth analysis of the FFRF is a new contribution to the existing literature. Beside the theoretical insights a set of guidelines is provided as practical benefit. Thereby engineers can decide which terms should be considered or not.

Illustrative examples used in the paper
In the practical use of multibody simulation, the flexible bodies are mostly modeled via the Finite Element Method. In that context, the displacements are evaluated at certain grid points and instead of shape functions, shape vectors are used to describe the elastic deformation. More comments on the transition from a "continuous formulation" to a formulation based on Finite Element models can be found in [1, Chap. 5.2, Sect. "Lumped mass"]. Although the details for the computation of the invariant matrices I 0 , W 1 , W 2 and W 3 may change, the final equations of motion (16) remain the same. Those matrices are computed based on the result of the Finite Element analysis and are part of the procedure when a flexible body is imported into the multibody simulation software. In MSC.ADAMS [14], for example, those matrices can be found in the so-called *.mtx file.
In this work several theoretical considerations are presented. The drawn conclusions are of fundamental character for all flexible bodies independent of the overall multibody system where the flexible bodies are embedded. Therefore, it is not necessary to present results of time integrations of particular multibody systems. Instead, the presented examples deal with estimations of effects based on the former mentioned invariant matrices, which are flexible body data only.
The Finite Element (FE) structures used to this end are introduced in this section. Six very different FE structures have been selected in order to cover a broad band of applications. All the assumptions mentioned in the former subsection are valid for the following FE structures. For all the Finite Element models given below, the consistent set of units, ton, millimeter and Newton, have been used. The necessary matrices I 0 , W 1 , W 2 and W 3 have been computed using the flexible body interface of MSC.ADAMS [14]. Consequently, the mode base of Craig [23] is used, which consists of constraint normal modes (CNM) and  constraint interface modes (CIM). Screenshots of the FE models can be seen in Fig. 2 and additional information is collected in Table 1.
For reasons of confidentially, the picture of the car body can only be shown in lowresolution. The generic structure GenStruct1 has been created in order to underline the effect of a particular structural compliance with respect to centrifugal forces. The second generic structure GenStruct2 is similar as GenStruct1 without having this special sensitivity to centrifugal forces.
The eigenvalues of the final modes are later used to estimate an upper limit of the modal coordinates. Therefore, the eigenfrequencies of the first ten modes are given in Fig. 3.

Assumption of small deformations
All investigated terms are related to effects caused by the flexibility of the body. These effects are increasing with the elastic deformation of the body. Therefore, it is important to specify a realistic upper bound of expectable deformations, and hence an upper bound of the modal coordinates.

Small elastic deformations with respect to the body's dimension
Commonly, the mode shapes are obtained by a linear Finite Element analysis. Consequently, the attribute "linear" defines the range of validity. The deformations have to be small enough, so that no geometric nonlinearities and no material nonlinearities take place. In the simple case of a tensile bar, the longitudinal deformation is given as l = σ l 0 E , where σ is the stress, l 0 the length of the bar and E the Young modulus. The use of an unrealistic high yield stress limit σ = 10 3 N/mm 2 and E = 2 × 10 5 N/mm 2 (steel) leads to a maximum displacement of 0.5% of the original length. In cases where slender structures like beams or plates are loaded by bending, geometric nonlinearities are the limiting criterion for the validity of a linear computation. In the case of a two-sided clamped beam or a four-sided clamped plate, geometric nonlinearities significantly influence the result even if the deflection of the center point is very small. The regime of validity is larger when just a one-sided clamped beam or plate is considered. However, in the literature it is often stated that for linear analysis, the maximum deflection should be in the range of the cross-section dimensions. To give an example, for the Finite Element solver in SolidWorks [9], it is suggested to use nonlinear techniques in case of deformations larger than 1/20 times the largest dimension of the body [10]. However, in this paper it is supposed that the elastic deformations are at least one to two magnitudes smaller than the body's largest extension.

Magnitude of modes and modal coordinates
Two very frequently used sets of consistent units are kilogram/meter and ton/millimeter, respectively. When those units are used, experience shows that the requirement of small deformations (no material and no geometric nonlinearity) leads to modal scaling factors (so-called "modal coordinates") typically much less than one. The following considerations show why this assumption mostly holds true.
A consequence of the mentioned mass normalization (Eq. (13)) is the relationship where m t holds the total mass of the structure andS p i represents a kind of averaged modal deflection of component i. It can be seen that the averaged modal deflection increases when the mass is decreasing, and vice versa. In the following, two extreme cases are investigated. For the first case, it is assumed that the modal deflections are uniformly distributed to all degrees of freedom. The averaged deflectionS p a1 can then be given as For the second case, it is assumed that the modal deflection is dominated by one component, like the first bending mode of a plate having only deflections normal to the plate plane. In that case the averaged deflectionS p a2 can be approximated as In Fig. 4 the averaged modal deflection is plotted vs total mass. In case when kilogram and meter are used as units, a mass of 1 kilogram leads to an averaged modal deflection of 1 meter. A structure of 100 kilograms still leads to an averaged deflection of around 0.1 meter. When ton and millimeter are in use, a mass of 0.001 tons (=1 kg) would lead to an averaged modal deflection of around 40 millimeters. A mass of 0.1 tons (=100 kg) implies an averaged deflection of 4 millimeters. Note once again that these are averaged deflections, the maximum deflection is probably at least a factor of 2 higher. In order to fulfill the restriction of small deformations, modal coordinates smaller than one seem to be a valid assumption. However, heavy structures may not fulfill that assumption.
The former conclusion that the modal coordinates are normally smaller than one depends on the units employed; in this case, kilogram and meter and ton and millimeter, respectively. This conclusion may not be valid for all combinations of length and mass units. Nevertheless, the investigations on the relevance of certain terms in the equations of motion are of general nature since these terms describe mechanical effects like the change of inertia due to deformation, coupling between rotation and deformation, centrifugal and Coriolis forces. These phenomena are of a fundamental nature and their relevance does not depend on the units.

Numerical examples
In order to estimate the potential maximum modal amplitudes of the investigated FE structures the maximum von Mises stress for the first 16 modes is given in Fig. 5. A modal coordinate of one would imply that the stress contribution of this mode is of the same magnitude, as shown in Fig. 5.
It can be seen that for almost all structures a modal amplitude lower than one can be expected, because otherwise material nonlinearities take place. Only the car and plate require a more detailed investigation.
Car The first mode of the car, scaled by 1, leads to a stress contribution of 258 N/mm 2 (von Mises). This is already a high stress level, but one may argue that a very good material could handle that. Figure 6 shows the deflection due to this mode for scaling factors 1 and −1. It can be seen that a vibration with such an amplitude is an unrealistic scenario for a passenger car during operation.
Plate The stresses inside the plate caused by the modes scaled with 1 are definitely in the linear range with respect to material nonlinearity. Figure 6 shows the deflection for the first mode when scaled with −1 and 1, together with the undeformed state (black). The maximum displacement of 42 mm is far beyond the plate thickness of 0.5 mm and occurs at the plates corners. It can be assumed that geometric nonlinear effects take place in case of such deformations. Therefore, the assumption that the modal coordinates remain smaller than 1 is probably a good estimate for that structure, too.
Based on the previous observations, the maximum modal coordinate for the first mode of the Finite Element structures given in this paper is 1 when the stress of the corresponding mode is lower than 500 N/mm 2 . Otherwise the maximum modal coordinate is chosen, so that the stress of the mode becomes 500 N/mm 2 . The maximum modal coordinate of all other modes is scaled by the ratio of the modes eigenvalue and the eigenvalue of the first mode. This is because the eigenvalue acts as stiffness (=Ω 2 ) in Eq. (14). The maximum modal coordinates of the first ten modes for the examples in this paper are depicted in Fig. 7. The latter figure reveals a general tendency: The stiffer and compacter a structure, the smaller the expected modal amplitudes. For very stiff and very compact structures, like the connecting rod, modal amplitudes much smaller than 1 can be expected.

Magnitude of entries in invariants W 1 , W 2 and W 3
In the key section of this paper, the significance of the terms holding the matrices W 1 , W 2 and W 3 is investigated. This chapter is devoted to the magnitude of the entries of the latter matrices in order to estimate whether a certain effect may be negligible with respect to another one.

Magnitude of entries in matrix W 1
Matrix W 1 is constructed using submatricesĪ kl , see (17). The components ofĪ kl are computed along (15), which leads to the evaluation of an integral in the form of Vū 0k S m l ρ dV . An application of the Cauchy-Schwarz inequality delivers According to (15), the terms ( V (ū 0k ) 2 ρ dV ) are defined as I kk , which are the main diagonal entries of I 0 . Equation (19) shows that the maximum value of V (S m l ) 2 ρ dV is 1. Due to inequality (23), it can be expected that the Euclidean norm of W 1 is normally significant less than the norm of I 0 . This assumption is confirmed by the numerical examples in the next section.

Numerical examples
The following quantities are computed for matrix W 1 and summarized in Table 2: • Maximum absolute entry.
• The average of the magnitude of all nonzero entries. This value gives an idea of whether most of the entries are around the maximum entry or much smaller. • Euclidean norms of W 1 and I 0 .
It can be seen in Table 2 that the Euclidean norm of I 0 is always greater than that of W 1 . Moreover, it is interesting to observe that the average value of the magnitude of all nonzero entries is always significantly less than the maximum entry. This leads to the conclusion that most of the entries in W 1 hold a magnitude which is significantly smaller than the maximum.

Magnitude of entries in matrices W 2 and W 3
Matrices W 2 and W 3 are constructed using the submatricesS kl (see (17)). The components ofS kl are computed as in (15), which leads to the evaluation of an integral in the form of Recalling Eq. (19), it is clear that the theoretical maximum value of V (S p k ) 2 ρ dV is 1. Therefore, the entries ofS kl have to be between −1 and 1. Consequently, the entries of W 2 and W 3 have to be between −2 and 2 (see (17)). Based on the evaluated integrals it can be stated that the norms of W 2 and W 3 are typically much lower than the norms of I 0 and W 1 . Due to the conservative assumption that the maximum value is 1, it can be expected that the magnitude of most of the entries is much smaller than 1. Both assumptions are confirmed by the numerical examples below.

Numerical examples
The following quantities are computed for matrices W 2 and W 3 and summarized in Table 3: • Maximum absolute entry.
• The average of the magnitude of all nonzero entries. This value gives an idea whether most of the entries are around the maximum entry or much smaller. • Euclidean norm.
It can be seen that the maximum values are always below 2, which is in accordance with the theoretical considerations. In comparison to the entries in Table 2, it turns out that the norms of W 2 and W 3 are smaller than those of I 0 and W 1 . It is notable that the average value of the magnitude of all nonzero entries is significant less than the maximum entry. This leads again to the conclusion that the entries of both matrices are mostly significantly smaller than the maximum entry. Finally, it can be stated that the entries of W 3 tend to be less than the entries of W 2 .

Investigations on the relevance of the single terms in the equations of motion
The following considerations are based on the full set of equations of motion (14), or when expressed in matrix notation (16). By setting all the terms related to W 1 , W 2 and W 3 to zero, the simplest possible representation of (16) is (25) An alternative way to obtain Eqs. (25) is by a strict application of the small deformation assumption to the kinetic energy (2). This means that in case of a sum of the body's extension and its deformation, the deformation is always neglected, as reported in [12] (and the reference therein to [13]). The obvious advantage of (25) is its simplicity and the fact that for the free body all coordinates are decoupled. The intention of this key chapter is to clarify whether and when all terms have to be considered (like (16)) or just the simplest possible form (like (25)), or something in-between has to be used.
For the following investigations, a free floating flexible body is assumed as described by (14) and (16). If the body is connected to other bodies via constraining or imposed forces, the situation is somehow different and the resulting consequences are discussed as well.

Theoretical considerations
In this section the change of the inertia due to deformation is investigated. The corresponding term in the mass matrix is the (n Θ ×n Θ ) matrix m ΘΘ . According to (14) and (16), this matrix can be given as The (3 × 3) matrix I ΘΘ holds the inertia of the deformed body. This matrix can be found in the first term of (Q v ) Θ in (14) as well. An evaluation of the integral in (26) gives where I ΘΘ,kk holds the kth main diagonal entry of I ΘΘ ,ū 0l holds the lth component of the vectorū 0 , andū f m holds the mth component of the vectorū f . The symbols k, l and m represent the three directions in space. The first integrand contains the inertia of the undeformed body and can be found in the matrix I 0 , which is a diagonal matrix due to the assumption of principal axis of inertia. The second integrand contains a linear change in inertia and is represented in matrix form by W 1 Q f . The third integrand contains a quadratic change of inertia and is given by Q T f W 2 Q f . Under the assumption that the deformation is one or two magnitudes smaller than the largest dimension and within the linear range, just two special scenarios can occur, for which the inertia tensor changes significantly due to the body's deformation: • The second and the third terms (W 1 Q f and Q T f W 2 Q f ) cannot be neglected if neither the body's dimension in the direction l nor m in (27) is significantly greater than the flexible deformations in the same directions. A long beam with a small cross-section would be such a structure. In this case, the deformation may have a similar magnitude as the body's dimension in that direction and the two state-dependent terms may be important for an accurate result. A close look at (27) gives no reason to assume that the quadratic term is less important than the linear part if such a beam-like structure is considered.
• The state-dependent terms may also be non-negligible in case of structures being particularly soft with respect to centrifugal forces, like GenStruct1. For such structures, small displacements may have a considerable impact on the moment of inertia. If such a structure is considered, it can be assumed that the linear term will dominate the quadratic term since the product of the body's extension and the elastic deformation will dominate the square of the elastic deformation, see Eq. (27). Note that the former two considerations are just meaningful when the body's state is dominated by its own inertia. Let us assume a body which fulfills one of the former two criteria for the consideration of the state-dependent terms. If such a body is stiffly connected to another body with remarkable inertia, it is sufficient to consider I 0 only. An example for such a situation is a rigid flywheel which is mounted on a flexible slender beam. As a general observation, it can be stated that the influence of W 1 Q f and Q T f W 2 Q f decreases when other bodies are stiffly connected to the flexible body under consideration.

Numerical examples
In order to quantify the influence of the deformation on the inertia, the following two quantities are defined: where I 0,kk holds the principle moment of inertia of the undeformed body around the axis k; I L,kk,p contains the result of W 1 Q f at the same position when the pth mode is scaled with a certain amplitude; I Q,kk,p gives the corresponding quantity due to Q T f W 2 Q f . Therefore I L,kk,max and I Q,kk,max are measures of the change in the inertia around axis k due to the linear and the quadratic part of Eq. (27). In Fig. 8, I L,kk,max and I Q,kk,max are given for the first ten modes of the investigated structures. Thereby the maximum modal amplitudes have been chosen as depicted in Fig. 7. It can be seen that all changes are below 5% except in GenStruct1 and the beam. This is in good accordance with the former theoretical considerations. The beam structure is characterized by a mass distribution which is dominated by one dimension in space, namely direction '3'. Therefore, I Q,33,max comes into the range of I 0,33 . It is interesting to note that the linear change is negligible in comparison to the quadratic one. As it was already assumed in the theoretical part, it is wrong to argue that the quadratic part is less important than the linear one in case of a beam-like structure.
The structure GenStruct1 is very soft around axis 1. Its first bending mode is depicted in Fig. 9 together with the undeformed structure. Due to the special mass distribution and the Fig. 9 First bending mode of GenStruct1 (scaled) particular softness, a remarkable effect on the inertia around axis 1 can be expected; even the deformation may remain in the linear range. Therefore I L,11,max and I Q,11,max come into the range of 10% of I 0,11 . A similar structure GenStruct2 does not show this effect because of its stiffer character.

Summary
The deformation-related change of a body's rotational inertia can be mostly neglected. Exceptions are only possible if the body has either a dominant extension in one direction or a very particular mass distribution in combination with particular (soft) stiffness properties. If one of those two situations takes place, an additional requirement needs to be fulfilled, which is that the body is a free body or at least very softly connected to other bodies, so that its state is mainly determined by its own inertia. Consequently, in nearly all industrial applications the terms W 1 Q f and Q T f W 2 Q f can be neglected, except when the previously mentioned criteria are fulfilled.

Theoretical considerations
In this section the coupling effect of m Θfq f on Θ (second line in equations of motion (14)) and m T ΘfΘ onq f (third line in equations of motion (14)) is investigated. For the subsequent investigations, the second line in (14) and (16) is pre-multiplied withḠ. This leads to and where the relationsḠḠ  (15), can be given as Note that the change of the inertia due to the body's deformation has already been investigated in the preceding section and is therefore neglected in the mass matrix and in the quadratic velocity vector. Another obvious observation is that the inertia coupling depends on the modal deflection but not on velocity, like the quadratic velocity vector does. As a consequence, it is not necessary to distinguish between moderate and high angular velocities, which is the case when the terms ofQ v are investigated. In the next step,q f of the first line in (29) is expressed by the second line, which leads to Expressingᾱ of the second line in (29) by the first line leads to In (34) the modal coordinates' inertia is reduced to the rotational coordinates while the rotational inertia is reduced to the flexible coordinates in (35). The reduced inertias can be given in integral and matrix form as and where I f →Θ has dimension (3 × 3) and I Θ→f has dimension (n f × n f ). The main diagonal entries of I f →Θ and I Θ→f can be given as
The scalar value I f →Θ,kk,p contains the result of Q T f W 3 W T 3 Q f at the main diagonal entry k due to a scaling of mode p. The scalar value I Θ→f,rr,p contains the result of 3 at the main diagonal entry r due to a scaling of mode p. A value of 1 indicates for both quantities an influence of the coupled inertia, which is of the same magnitude as the inertia of the corresponding degrees of freedom.
In Fig. 10, the quantities defined by (40) are given for the first ten modes of the investigated structures. Thereby the maximum modal amplitudes have been chosen as depicted in Fig. 7. It can be seen that only the structures GenStruct1 and "Beam" have the potential that the inertia coupling may significantly influence the equations of motion. The same two structures have been denoted as "critical" in the preceding subsection (see Fig. 8). This was predicted due to the theoretical considerations before. If we consider the conservative assumption for the maximum modal displacement, just the beam remains potentially sensitive to the inertia coupling.

Summary
The inertia coupling is relevant for the same type of structures, which are critical in terms of deformation induced inertia changes. The influence of the inertia coupling decreases with the number of stiff connections to other bodies. Similar to the preceding section, it can be concluded that inertia coupling is an effect of minor practical relevance.

Quadratic velocity vector
The quadratic velocity vector is formed by centrifugal, Coriolis and Euler forces. The latter evaluate to zero when Euler parameters are used. In such cases, along with the assumptions on which this paper is based, the quadratic velocity vector can be repeated as (see (14) and (30)). The centrifugal forces are described by the first terms in each line, which contain the square of the angular velocity. The Coriolis forces are given by the terms where a product of the angular and modal velocities can be found. A subdivision of (41) into centrifugal and Coriolis forces gives The quadratic velocity vector is responsible for the so-called rotordynamics effects, such as the gyroscopic effect and a meaningful coupling between rotational and elastic degrees of freedom. There is a lot of literature available on this subject (see Gasch [17]). From the literature it is known that rotor dynamic effects are of minor importance when the quotient of rotational speed and lowest eigenfrequencies of the assembled system is small enough. Note that the lowest eigenfrequencies of the assembled system are the relevant ones and not the first eigenfrequencies of a free single flexible body which is the square root of the first entries in the Ω 2 matrix. An example would be a multibody system where a rigid disk is constrained to an elastic beam. In that case, the lowest eigenfrequencies of the entire system (beam plus disk), depend on the inertia data of the disk, the position where the disk is fixed to the beam, and on the mounting of the beam.

Theoretical considerations
Rotational part of quadratic velocity vector The rotational part of the quadratic velocity vector can be given in integral form as and in matrix form as The subdivision into the centrifugal and Coriolis forces can be written as The centrifugal force holds the inertia tensor of the deformed body and is responsible for the gyroscopic effects. The importance of the change in inertia due to the deformation of the body has already been discussed in the section "Rotational inertia of the deformed body". In this chapter it has been shown that the terms W 1 Q f and Q T f W 2 Q f are only important for beam-like structures and for structures which are soft with respect to centrifugal forces. Therefore, the terms W 1 Q f and Q T f W 2 Q f can almost always be neglected and the inertia tensor of the undeformed body is sufficient. In this work it is suggested to always consider the gyroscopic effects based on I 0 , since the angular velocity needs to be computed anyway and, hence, the computation of the centrifugal force is cheap. Furthermore, it is worth mentioning that such a consideration of the centrifugal force does not couple the rotational and flexible degrees of freedom of the free flexible body.
The Coriolis force couples the modal and the rotational degrees of freedom via a product of the angular and the modal velocities. It stems from the fact that the elastic body is characterized with respect to a co-rotating reference frame. It can be observed that all three terms of the Coriolis force consist of a matrix-matrix or a matrix-vector product of the angular and modal velocities together with a kind of inertia. One term holds a constant inertia (W 1 ) while the other terms hold an inertia which is scaled by the modal coordinates (Q T f W 2 and . From a former section it is known that the entries of W 1 are typically significantly higher as those of W 2 and W 3 . Moreover, the latter two matrices are scaled by modal coordinates smaller than 1, which leads to even smaller values. Based on the insights of the former sections, it can be expected that exceptions may occur in case of beam-like bodies and when a body is particularly soft with respect to the centrifugal force.

Flexible part of quadratic velocity vector
In the next subsection, the flexible part of the quadratic velocity vector is investigated in terms of its relevance for the overall solution.
The (n f × 1) vectorQ v,f contains the effect of centrifugal and Coriolis forces on the modal coordinates. This vector evaluates to zero, when the small deformation assumption is applied to the kinetic energy; see (25). The flexible part of the quadratic velocity vector can be given asQ or in matrix form A subdivision into the centrifugal and Coriolis forces can be written as A well-known effect of the centrifugal force is the elastic expansion (widening) of a flexible body when it rotates around one axis. The general question of considering or neglecting this kind of deformation cannot be answered in general and is application dependent. In case of compliant structure with respect to the centrifugal force, like GenStruct1, it seems obvious to consider the deformation induced by the centrifugal force. However, even small deformations of actually stiff structures can be important in some applications. Such an example would be the fatigue life prediction of a crankshaft based on modal stress recovery (see [15] and [16]). Although the crankshaft deformation due to the rotational speed is small, the thereby induced stress can have significant impact on the fatigue life prediction. This is because it acts as kind of mean stress which has considerable impact on fatigue lifetime. Elastic deformations caused by centrifugal forces are a fundamental effect of rotating elastic bodies. Their complete negligence only makes sense when either the rotational speed is low or the body is stiff and the neglected small deformations are non-critical. It can be seen that the centrifugal force consists of a part with constant inertia (=W 1 ) and a second one where the matrix W 2 is scaled by the modal coordinates. As in the section before, it can be stated that the entries of W 1 are much larger than those of W 2 . Moreover, W 2 is multiplied with modal coordinates, which are typically smaller than one. Therefore, the constant term normally dominates the state-dependent one. The numerical examples below underline the latter considerations. Exceptions may be special structures, which are particularly soft in terms of centrifugal forces or beam-like bodies. It may be interesting to note that the term holding W 2 represents a reduction of the structure's effective stiffness since the centrifugal force leads to an expansion of the body, which leads to a higher centrifugal force. If the rotational speed exceeds a certain limit, the effective stiffness becomes negative and the system is instable. This is a more academic scenario, however, because under such conditions the structure will not fulfill the small deformation assumption.
Once again, it is emphasized that the former considerations do not contain the centrifugal and Coriolis forces due to other bodies which are somehow stiffly connected to the flexible body under consideration. Even if the body is just connected to ground the mounting forces may influence the body's behavior much more than the latter discussed quadratic velocity vector. Again, a tendency can be assumed: The more the flexible body interacts with the ground or other bodies, the more the previously discussed terms become irrelevant for the overall solution.

Numerical examples
Rotational part of quadratic velocity vector In the discussion of the Coriolis force, it has been assumed that the term holding W 1 is normally dominating. Exceptions are expected in cases of slender structures and structures with a special mass distribution and an extraordinary sensitivity to centrifugal forces.
In Table 4 the ratio of the Euclidian norm of the deformation-dependent terms with respect to the non-deformation-dependent term can be seen. For the construction of Q T f the maximum modal deformation is used for the first ten modes, see Fig. 7. The modal coordinates of the remaining modes have been set to zero. For the construction of the modal velocities it has been assumed that those ten modes vibrate with the modes' eigenfrequencies. The angular velocity has been set toω T = [ 1 1 1]. Note that a scaling of the angular velocity does not change the ratio given in Table 4. It can be seen that the norm of the deformation-dependent terms is of two to three magnitudes smaller than that with the nondeformation-dependent term, except in the case of the Beam and GenStruct1. This is in accordance with the theoretical considerations mentioned before. It should be emphasized at this point that the assumptions made for the latter ratio are very conservative. Firstly, the maximum modal coordinates themselves are conservative (see the discussion around Fig. 7) and secondly, the fact that all maximum coordinates and velocities act (without phase shift) simultaneously is even more conservative. Therefore, it seems that in reality the ratio of the norms is expected to be even smaller.
Flexible part of quadratic velocity vector It has been assumed that the state-dependent part of the centrifugal force can normally be neglected in comparison to the constant one. Table 5 contains the ratio of the Euclidean norm of the two parts of the centrifugal force when the maximum modal coordinates are applied for the first ten modes simultaneously. It can be seen that the contribution of the deformation-dependent term is at least of two magnitudes smaller than that with the constant term, excluding GenStruct1. This was expected because of the special characteristics of this structure.  From rotor dynamics it is known that the Coriolis force can be neglected when the ratio of the assembled systems' first eigenfrequency and the rotational speed are considerable lower than 1. If this is not the case, the Coriolis forces should be taken into account, whereby, for the rotational degrees of freedom, the term W 1Qfω is sufficient. Exceptions are beamlike structures and, once again, rotational soft structures in the sense of GenStruct1, when the inertia forces of those structures are dominated by themselves and not by other bodies which are stiffly attached. which are obtained by the strict application of the small deformation assumption on the level of the kinetic energy can be used in the following cases: • Moderate angular velocities when the body is neither a beam nor an extraordinary rotational soft structure (in the sense of GenStruct1). Examples are ground vehicles, housings, aircrafts, and so on. • In cases of flexible bodies which are stiffly connected to other bodies, so that the effective inertia or effective stiffness of the assembled system is different to that of the free body. Examples are crankshafts, connecting rods, or a beam with a mounted flywheel.
In industrial operations, the latter cases will probably cover most the applications. II. If, in addition, the widening due to the centrifugal force is of interest, Eq. (50) should be considered. In this equation the centrifugal force acting on the flexible coordinates is additionally regarded. ⎡ III. In case of rotating structures, the following extraordinary situations can require additional terms: • If the body under consideration is very soft with respect to centrifugal forces (in the sense of GenStruct1) and this softness is not influenced by connections to other bodies (or ground), the linear change in inertia should be considered in the form of IV. The Coriolis force should be considered when the ratio of the assembled system's first eigenfrequency and the rotational speed is not considerably lower than 1 and when the inertia of the rotating body is determined by itself and not by other bodies which are stiffly connected. In such a case, the term −Ḡ T W 1Qfω is normally sufficient for the rotational degrees of freedom. The full expression (−Ḡ T [W 1 + 2Q T f W 2 ]Q fω − 2Ġ T Q T f W 3q f ) needs to be regarded only in cases with beam like structures or rotational soft structures, like GenStruc1. For the flexible coordinates the Coriolis force is then covered by the term −2W T 3Qfω .

Benefit
The well-founded negligence of certain terms in the equations of motion leads to the following benefits: The first benefit is obviously the reduction of the computational burden. The matrixes W 1 , W 2 and W 3 are involved in matrix-matrix and matrix-vector products which need to be computed at each iteration during the numerical time integration. This has to be done for the equations of motions and also for the computation of the Jacobian matrix. When terms including W 1 , W 2 and W 3 can be neglected, the simulation time decreases significantly in cases where models with a high number of flexible bodies and/or when flexible bodies with many modes are considered. Note that the latter becomes more and more relevant because there is a tendency to consider more and more local effects in flexible multibody dynamics (see, for example, [18][19][20][21] and [26]). In order to illustrate the potential, two examples are given. The first example deals with an elastic piston rod simulation where the mode base of Craig [23] is extended by so called "contact modes" so that in total 90 modes are in use. Those "contact modes" enable an efficient and accurate consideration of the nonlinear contact forces in the bearing cap, see Sect. 4.2. of [26] for more details. Two simulations are performed with MSC.ADAMS: one with the default setting "partial coupling" and another with "full coupling". In the first case just the terms holding W 1 are considered, while in the second case all terms of the equations of motion are considered. Even though the results are almost identical, the difference with respect to the CPU time is almost 50%. For the second example, the freely available code FreeDyn [27] is used. A 6 cylinder V-engine as in Fig. 11 has been constructed with rigid bodies and 6 flexible piston rods. For the piston rod, a Craig mode base [23] is used with 12 constraint interface modes and 15 constraint normal modes. The crankshaft is mounted to ground by a revolute joint and a moment accelerates the entire mechanism up to 8000 rpm. The same two kinds of simulation are performed as in the former example. It can be reported that the difference in the CPU time is around 15%, while the modal coordinates of both computations are almost identical.
A second benefit is that the computation of W 1 , W 2 and W 3 is not necessary at all when they are not needed. Consequently, the eigenvalues and mode shapes are enough for flexible multibody dynamics. This is standard output of many Finite Element codes, while W 1 , W 2 and W 3 are not. This simplifies the numerical implementation of the equations of motion enormously.
The third benefit is that negligence of the terms including W 1 , W 2 and W 3 leads to equations of motion where the mass matrix and the quadratic velocity vector are decoupled with respect to the translational, rotational and flexible degrees of freedom. This fact simplifies considerations with respect to separated time integration [19] or model reduction of multibody systems [22] tremendously.
Finally, the set of guidelines given above removes uncertainty concerning the question of which invariant needs to be considered and which not. In available software packages it is common to commit this decision to the user without clear guidance. The former suggestions can be used as such a guidance in order to activate or deactivate the proper invariants.

Conclusion
In this work the significance of all inertia related terms of a flexible body in the FFRF were investigated at the level of the equations of motion. It turned out that for a lot of applications a remarkably simple and decoupled set of equations are sufficient. This has already been suggested in the literature when the small deformation assumption is strictly applied at the level of the kinetic energy (see [12] or [21]). However, there are situations which require either the deformation-dependent inertia tensor, or the inertia coupling, or the centrifugal and Coriolis forces which are related to the elastic deformation. All of these terms were investigated with respect to their significance and, finally, these results were condensed in a set of guidelines which give simple advice about which term needs to be considered and which not. All theoretical considerations have been underlined by simple numerical investigations which have been applied to a couple of very different Finite Element structures.