# Modal Mass and Length of Mode Shapes in Structural Dynamics

Authors: M. Aenlle; Martin Juul; R. Brincker
Journal: Shock and Vibration
 (2020)

---

## Abstract

The literature about the mass associated with a certain mode, usually denoted as the modal mass, is sparse. Moreover, the units of the modal mass depend on the technique which is used to normalize the mode shapes, and its magnitude depends on the number of degrees of freedom (DOFs) which is used to discretize the model. This has led to a situation where the meaning of the modal mass and the length of the associated mode shape is not well understood. As a result, normally, both the modal mass and the length measure have no meaning as individual quantities but only when they are combined in the frequency response function. In this paper, the problems of defining the modal mass and mode shape length are discussed, and solutions are found to define the quantities in such a way that they have individual physical meaning and can be estimated in an objective way.


---

## Body

## 1. Introduction

The classical equation of motion for a system withN degrees of freedom (DOFs) is [1–10](1)Myt+Cyt+Kyt=pt,where pt is the force input vector, yt is the response vector, and M,C, and K are the mass, the damping, and the stiffness matrices, respectively.Assuming proportional damping, it is well known that the solution is given as [1–10](2)yt=Ψqt=ψ1q1t+ψ2q2t+ψ3q3t+….where Ψ=ψ1ψ2ψ3… is the mode shape matrix containing the real-valued mode shapes representing the spatial solution and qt is the modal coordinate vector containing the modal coordinates qt=q1tq2t…T representing the time solution. Inserting the solution given by equation (2) into the equation of motion, equation (1), and multiplying from the right by the mode shape matrix transpose, it can be shown that the resulting matrices are diagonalizable, and we obtain a set of equations of independent 1-DOF systems, one for each of the modal coordinates, i.e.,(3)mrq˙t+crq˙t+krqt=prt,where mr, cr, kr, and prt=ψrTpt are the modal mass, the modal damping, the modal stiffness, and the modal load, respectively, corresponding to the r-th mode.In the frequency domain, the equation of motion is expressed as [1–10](4)Hω⋅Pω=Uω,where Hω is the frequency response function matrix given by(5)Hω=−ω2M+iωC+K−1.The FRF can also be expressed in terms of modal parameters, which for proportional damped models, it is given by [1–4](6)Hω=∑r=1NMψrψrTmrωr2−ω2+i2ζrωωr,where NM is the number of modes, ωr is the natural frequencies, ζr is the damping ratio, ψr is the mode shape, and mr is the modal mass, respectively, of the r-th mode. From equations (3) and (5), it is easily seen that the modal mass is needed in all applications where the frequency response function (FRF) (or the impulse response function (IRF)) has to be constructed from the modal parameters, such as structural modification, health-monitoring applications, and damage detection [1–4, 10].A mode shape is said to be mass-normalized when the modal mass is dimensionless, i.e.,m=1. On the contrary, a mode shape is said to be unscaled if it is not mass-normalized. The mass normalized ϕ and the unscaled ψ mode shapes are related by the following equation [1–10]:(7)ϕ=1mψ,where the modal mass m is a real scalar in undamped and proportionally damped models. If the mode shape ψ is dimensionless, the modal mass has units of mass kg.The concept of modal mass is addressed in the classical books of structural dynamics [5, 6], modal analysis [1–4, 10], and also in research papers [11, 12]. In some books of classical dynamics [5, 13], the concept of modal mass has also been named as generalized mass or effective modal mass.


### 1.1. The Concept of Effective Mass

Ewins [1] defined a new parameter denoted as effective mass which is different from the concept of modal mass. This effective mass is related to a certain mode and a certain DOF. If mass-normalized mode shapes are used, the effective mass at DOF j for mode r is defined as(8)mjjr=1ϕjr2.Due to the fact that mass-normalized mode shapes are unique and not subject to any arbitrary scaling factors, this new concept is also unique and represents a useful description of the behavior of the structure point by point and mode by mode [1].In base-excited systems, the term effective mass has also been used to define a different parameter [5, 6, 14], which must not be mistaken with the modal mass.For a structure subjected to an acceleration support excitation, with a mass matrixM, mode shapes ψ, and influence coefficient l, the modal participation factor for the r-th mode can be obtained as [5, 6, 14](9)Γr=ψrT⋅M⋅lψrT⋅M⋅ψr=ψrT⋅M⋅lmr,which provides a measure of the system mass participating in that particular mode. In engineering, the earthquake vector l represents the displacements of the story masses resulting from a static unit ground displacement in the direction of the seismic excitation [6].The effective modal massmeffr for the r-th mode is defined as(10)meffr=Γr2mr.The concept of effective mass was firstly proposed by Bamford et al. [14] in the early 70s, and it is based on the assumptions that the system is excited at the base and also that the base is rigid. A mode with a large effective mass is usually a significant contributor to the response of the system m, i.e., only the modes which have significant effective masses are needed to represent the response of the structure in a certain frequency band [15, 16]. Moreover, the sum of the effective masses for all modes in a given response direction must equal the total mass of the structure [5, 6]. If some modes of the system are truncated, the effective mass of the truncated modes is added directly to the base as a residual mass [15]. The equations of motion obtained for displacement and acceleration loadings are different from each other because of the construction of the effective external load [17, 18]. When multiple constraints need to be defined by different motions, displacement loading should be the preferred one.


### 1.2. Modal Mass and Modal Analysis

Most of the literature related to modal mass is devoted to methodologies and techniques to estimate the modal masses of structures and mechanical systems from experiments. In [19], it is mentioned that modal mass is the least reliable parameter when classical modal analysis (CMA, also known as experimental modal analysis) is being used. Furthermore, it is also very sensitive to response magnitude.In classical modal analysis, the modal masses are obtained from the experimental FRF using some of the several identification techniques, which are fully described in the classical books of modal analysis [1–3]. Although it is well known that a relatively high uncertainty is expected in the modal masses when using CMA, only few papers can be found in the literature discussing this subject.Zivanovic et al. [19] estimated the modal mass corresponding to the first vertical mode of a footbridge with CMA using two different excitations (random excitation and unaveraged chirp excitation) obtaining 53.200 kg and 48.800 kg, respectively, the difference being 8.3%.Allen and Sracic [20] applied CMA to determine the modal parameters in a free-free beam, excited with an impulse hammer measuring the response with a scanning laser vibrometer. The modal masses for all the modes were estimated with an error less than 20%. Much of the discrepancy was attributed to the variation in the hammer blows and uncertainty in the damping estimates.In operational modal analysis (OMA), the forces are unknown, and the modal masses cannot be estimated, i.e., only the unscaled mode shapes can be identified for each mode [4, 21]. In order to overcome this problem, several approaches have been proposed recently. One is based on modifying the dynamic behavior of the structure by changing the stiffness and/or the mass and then performing operational modal analysis on both the original and the modified structure [21–25]. Another one is using the mass matrix of a finite element model [26] or considering for the experimental mode shapes having the same modal masses as those of a finite element model [27].Several methods have also been proposed to estimate modal masses combining OMA and the response of the structure subject to a certain artificial excitation. Pavic and Brownjohn [28] developed a method for estimating the modal masses of footbridges known to respond to pedestrian excitation. Their method is based on quantifying the dynamic forces that generate the lively response, along with a characterization of the initial part of the resonant build-up due to these forces. In [29], the modal parameters of some British rock lighthouses are estimated with CMA. The most important modes have frequencies ranging between 4 Hz and 7 Hz, and the modal masses are of the order of 200 t.Brandt et al. [30] proposed the OMAH technique where operational modal analysis is combined with monoharmonic excitation applied by an actuator. Cara [31] developed the equations needed to compute the modal masses from a state space model when CMA and OMA are combined. Guillaume et al. [32] proposed the OMAX technique, which is an operational modal analysis approach where unknown forces as well as measurable input forces, which can be deterministic or stochastic, are available.Hwang et al. [33] proposed a method for estimating the modal masses of a structure using a mass-type damping device. Controllability and observability matrices are constructed using the identified system matrices and modal space system matrices. The modal masses are obtained considering that the product of controllability and observability matrices does not change with the type of system matrices.Füllekrug [34] developed the equations for estimating effective masses and modal masses from base forces measured in base-driven tests. The author obtained accurate results, and the accuracy only decreases in cases where the magnitude of the effective masses is low.


### 1.3. Tuned Mass Dampers (TMDs)

A tuned mass damper (TMD) is a dynamical device (it consists of a mass, a spring, and a damper) that is attached to a structure in order to reduce the dynamic’s response of a structure [35–39]. The mass and the stiffness of the TMD are tuned to the natural frequency that needs to be damped, typically the first mode.A TMD is characterized by three ratios: frequency ratio (natural frequency of the TMD to the natural frequency of the structure), mass ratio (mass of the TMD to the modal mass of the structure), and damping (damping of the TMD to the damping ratio of the structure) [35–37]. Thus, the modal mass of the structure plays an important role in the design of tuned mass dampers. As the mass of the damper depends on the modal mass of the structure, how the modal mass is determined is an important factor in the design of a TMD [40].The tuning of the properties of the TMD depends on the loading acting on the structure. Den Hartog and Ormondroyd [35] proposed a method to minimize the response to sinusoidal loading. The first successful analysis of the TMD for seismic loading was introduced by Wirsching and Yao [36].The effect of considering soil-structure interaction in the seismic response of reinforced concrete chimneys was investigated by Elias [41], where it is concluded that the increase in the mass ratio reduces the seismic response. An integrated damping system using both TMDs and double skin facade (DSF) damping system was investigated in [42]. The author concluded that this system requires a significantly reduced TMD mass ratio compared to traditional TMDs. The response of the wind turbine tower under wind loads was investigated by Gaur et al. [37]. They concluded that an optimum damping ratio exists for each mass ratio.


### 1.4. Mode Shape Normalization

Mode shapes can be normalized in many different ways, the most common techniques being mass normalization, normalization to the unit length of the mode shape (length scaling), and normalization to a component (usually to the largest component) equal to unity (DOF scaling). For this reason, the modal mass of a mode shape is not unique because it is directly related to the normalization method which has been used to define the mode shape [1], making the physical interpretation difficult.Due to the fact that the modal mass corresponding to the mass-normalized mode shapesϕ is m=1, the so-defined dimensionless modal mass introduces problems interpreting the modal equations given by equation (3) and the modal coordinates that now have the unit m/kg. This is misleading and unfortunate.In [1–3], the modal mass is defined as a scaling parameter for the mode shapes, i.e., it is used to convert the original unscaled mode shape vector ψ to the scaled (mass-normalized) mode shape vector ϕ. The concept of scaling factor α, which is related to the modal mass by(11)α=1mis also used in the field of structural dynamics [4].Expressions for calculating the modal mass in continuous beams can be found in the literature [5–9]. As an example, in continuous straight planar beams with length L, distributed mass density ρx, and cross section with area Ax, the modal mass (also denoted as generalized mass in some books of structural dynamics) corresponding to an arbitrary normalized continuous mode shape vector ψx is given by [5–9](12)mψ=∫0LρxAxψx2dx.In discrete systems, the mode shapes are vectors, and the modal mass corresponding to the mode shapeψ can be calculated with the following equation [1–8]:(13)m=ψTMψ,where M is the mass matrix.A mode shape is said to be normalized to the unit length if its length is unity. In one-dimensional continuous systems, the Euclidean lengthLE of a function ψx, also known as Euclidean norm or L2−norm, is defined as [43](14)LE=∫0Lψx2dx,where the subindex ‘E’ indicates Euclidean.In discrete systems, the length of the mode shape vectorψ (length of a vector in a Euclidean space) is defined as [43](15)LE=ψTψ.A mode shapeψ can be scaled to the unit length by(16)ψL=ψLE,where ψL indicates the mode shape normalized to the unit length (Euclidean length equals to unity).It is easy to realize that this procedure will define the modal mass estimated from equation (13) as something that is heavily dependent upon the number of DOFs in the model. In fact, it is easy to show that, for this case, the modal mass is approximately proportional to(17)m∝1N.Another approach to scaling is to scale the largest component (or some other component) to a certain value (for instance, equal to unity) [1–7]. This introduces the DOF scaled mode shape ψD:(18)ψD=ψψDn,where ψDn is the considered DOF used for scaling.If the mode shapes are normalized to the largest component equal to unity, the DOF used for normalization should be the same in the FE model and the test, i.e., the DOF must be shared. If this is not the case, then the modal mass is dependent upon the selected DOF.


### 1.5. Physical Interpretation of Modal Mass

Although the concept of modal mass is defined in the classical books of structural dynamics and modal analysis, the physical interpretation of this modal parameter is rarely addressed.The modal mass for the undamped and the proportionally damped cases is usually defined as a scaling factor for the mode shapes [1, 3, 11], i.e., it is used to convert the original unscaled mode shape vector {ψ} to the more useful mass-normalized mode shape vector {ϕ} [1–4].When using modal superposition (see equation (3)), the modal mass represents the mass of the single DOF systems, after decoupling the equations by transformation to modal coordinates [1–4].In [44], it is demonstrated that the modal mass of a simply-supported beam, when the mode shapes are normalized to the largest component equal to unity, is the half of the total mass of the system, whereas in fixed-fixed beams, the modal mass is equal to the total mass of the system.However, it seems like none of the classical ways of dealing with the scaling problem of the mode shapes leads to a modal mass that has the needed clarity in terms of physical meaning. Furthermore, since the results are dependent upon the number of DOFs in the model, the individual quantities do not seem to have much meaning at all.In this paper, a new and better definition of the modal mass is presented. This definition leads to significant improvements in the physical understanding of this modal parameter. Firstly, the most important requirements for a new definition of the modal mass are proposed. Then, a new expression for the length of a continuous mode shape is formulated introducing the concept of the volume matrix, which again leads to a new definition of the modal mass for the continuous case. The expressions derived for continuous systems are later extended to discrete cases. Finally, the concepts and equations formulated in the paper are illustrated by simple examples.


## 1.1. The Concept of Effective Mass

Ewins [1] defined a new parameter denoted as effective mass which is different from the concept of modal mass. This effective mass is related to a certain mode and a certain DOF. If mass-normalized mode shapes are used, the effective mass at DOF j for mode r is defined as(8)mjjr=1ϕjr2.Due to the fact that mass-normalized mode shapes are unique and not subject to any arbitrary scaling factors, this new concept is also unique and represents a useful description of the behavior of the structure point by point and mode by mode [1].In base-excited systems, the term effective mass has also been used to define a different parameter [5, 6, 14], which must not be mistaken with the modal mass.For a structure subjected to an acceleration support excitation, with a mass matrixM, mode shapes ψ, and influence coefficient l, the modal participation factor for the r-th mode can be obtained as [5, 6, 14](9)Γr=ψrT⋅M⋅lψrT⋅M⋅ψr=ψrT⋅M⋅lmr,which provides a measure of the system mass participating in that particular mode. In engineering, the earthquake vector l represents the displacements of the story masses resulting from a static unit ground displacement in the direction of the seismic excitation [6].The effective modal massmeffr for the r-th mode is defined as(10)meffr=Γr2mr.The concept of effective mass was firstly proposed by Bamford et al. [14] in the early 70s, and it is based on the assumptions that the system is excited at the base and also that the base is rigid. A mode with a large effective mass is usually a significant contributor to the response of the system m, i.e., only the modes which have significant effective masses are needed to represent the response of the structure in a certain frequency band [15, 16]. Moreover, the sum of the effective masses for all modes in a given response direction must equal the total mass of the structure [5, 6]. If some modes of the system are truncated, the effective mass of the truncated modes is added directly to the base as a residual mass [15]. The equations of motion obtained for displacement and acceleration loadings are different from each other because of the construction of the effective external load [17, 18]. When multiple constraints need to be defined by different motions, displacement loading should be the preferred one.


## 1.2. Modal Mass and Modal Analysis

Most of the literature related to modal mass is devoted to methodologies and techniques to estimate the modal masses of structures and mechanical systems from experiments. In [19], it is mentioned that modal mass is the least reliable parameter when classical modal analysis (CMA, also known as experimental modal analysis) is being used. Furthermore, it is also very sensitive to response magnitude.In classical modal analysis, the modal masses are obtained from the experimental FRF using some of the several identification techniques, which are fully described in the classical books of modal analysis [1–3]. Although it is well known that a relatively high uncertainty is expected in the modal masses when using CMA, only few papers can be found in the literature discussing this subject.Zivanovic et al. [19] estimated the modal mass corresponding to the first vertical mode of a footbridge with CMA using two different excitations (random excitation and unaveraged chirp excitation) obtaining 53.200 kg and 48.800 kg, respectively, the difference being 8.3%.Allen and Sracic [20] applied CMA to determine the modal parameters in a free-free beam, excited with an impulse hammer measuring the response with a scanning laser vibrometer. The modal masses for all the modes were estimated with an error less than 20%. Much of the discrepancy was attributed to the variation in the hammer blows and uncertainty in the damping estimates.In operational modal analysis (OMA), the forces are unknown, and the modal masses cannot be estimated, i.e., only the unscaled mode shapes can be identified for each mode [4, 21]. In order to overcome this problem, several approaches have been proposed recently. One is based on modifying the dynamic behavior of the structure by changing the stiffness and/or the mass and then performing operational modal analysis on both the original and the modified structure [21–25]. Another one is using the mass matrix of a finite element model [26] or considering for the experimental mode shapes having the same modal masses as those of a finite element model [27].Several methods have also been proposed to estimate modal masses combining OMA and the response of the structure subject to a certain artificial excitation. Pavic and Brownjohn [28] developed a method for estimating the modal masses of footbridges known to respond to pedestrian excitation. Their method is based on quantifying the dynamic forces that generate the lively response, along with a characterization of the initial part of the resonant build-up due to these forces. In [29], the modal parameters of some British rock lighthouses are estimated with CMA. The most important modes have frequencies ranging between 4 Hz and 7 Hz, and the modal masses are of the order of 200 t.Brandt et al. [30] proposed the OMAH technique where operational modal analysis is combined with monoharmonic excitation applied by an actuator. Cara [31] developed the equations needed to compute the modal masses from a state space model when CMA and OMA are combined. Guillaume et al. [32] proposed the OMAX technique, which is an operational modal analysis approach where unknown forces as well as measurable input forces, which can be deterministic or stochastic, are available.Hwang et al. [33] proposed a method for estimating the modal masses of a structure using a mass-type damping device. Controllability and observability matrices are constructed using the identified system matrices and modal space system matrices. The modal masses are obtained considering that the product of controllability and observability matrices does not change with the type of system matrices.Füllekrug [34] developed the equations for estimating effective masses and modal masses from base forces measured in base-driven tests. The author obtained accurate results, and the accuracy only decreases in cases where the magnitude of the effective masses is low.


## 1.3. Tuned Mass Dampers (TMDs)

A tuned mass damper (TMD) is a dynamical device (it consists of a mass, a spring, and a damper) that is attached to a structure in order to reduce the dynamic’s response of a structure [35–39]. The mass and the stiffness of the TMD are tuned to the natural frequency that needs to be damped, typically the first mode.A TMD is characterized by three ratios: frequency ratio (natural frequency of the TMD to the natural frequency of the structure), mass ratio (mass of the TMD to the modal mass of the structure), and damping (damping of the TMD to the damping ratio of the structure) [35–37]. Thus, the modal mass of the structure plays an important role in the design of tuned mass dampers. As the mass of the damper depends on the modal mass of the structure, how the modal mass is determined is an important factor in the design of a TMD [40].The tuning of the properties of the TMD depends on the loading acting on the structure. Den Hartog and Ormondroyd [35] proposed a method to minimize the response to sinusoidal loading. The first successful analysis of the TMD for seismic loading was introduced by Wirsching and Yao [36].The effect of considering soil-structure interaction in the seismic response of reinforced concrete chimneys was investigated by Elias [41], where it is concluded that the increase in the mass ratio reduces the seismic response. An integrated damping system using both TMDs and double skin facade (DSF) damping system was investigated in [42]. The author concluded that this system requires a significantly reduced TMD mass ratio compared to traditional TMDs. The response of the wind turbine tower under wind loads was investigated by Gaur et al. [37]. They concluded that an optimum damping ratio exists for each mass ratio.


## 1.4. Mode Shape Normalization

Mode shapes can be normalized in many different ways, the most common techniques being mass normalization, normalization to the unit length of the mode shape (length scaling), and normalization to a component (usually to the largest component) equal to unity (DOF scaling). For this reason, the modal mass of a mode shape is not unique because it is directly related to the normalization method which has been used to define the mode shape [1], making the physical interpretation difficult.Due to the fact that the modal mass corresponding to the mass-normalized mode shapesϕ is m=1, the so-defined dimensionless modal mass introduces problems interpreting the modal equations given by equation (3) and the modal coordinates that now have the unit m/kg. This is misleading and unfortunate.In [1–3], the modal mass is defined as a scaling parameter for the mode shapes, i.e., it is used to convert the original unscaled mode shape vector ψ to the scaled (mass-normalized) mode shape vector ϕ. The concept of scaling factor α, which is related to the modal mass by(11)α=1mis also used in the field of structural dynamics [4].Expressions for calculating the modal mass in continuous beams can be found in the literature [5–9]. As an example, in continuous straight planar beams with length L, distributed mass density ρx, and cross section with area Ax, the modal mass (also denoted as generalized mass in some books of structural dynamics) corresponding to an arbitrary normalized continuous mode shape vector ψx is given by [5–9](12)mψ=∫0LρxAxψx2dx.In discrete systems, the mode shapes are vectors, and the modal mass corresponding to the mode shapeψ can be calculated with the following equation [1–8]:(13)m=ψTMψ,where M is the mass matrix.A mode shape is said to be normalized to the unit length if its length is unity. In one-dimensional continuous systems, the Euclidean lengthLE of a function ψx, also known as Euclidean norm or L2−norm, is defined as [43](14)LE=∫0Lψx2dx,where the subindex ‘E’ indicates Euclidean.In discrete systems, the length of the mode shape vectorψ (length of a vector in a Euclidean space) is defined as [43](15)LE=ψTψ.A mode shapeψ can be scaled to the unit length by(16)ψL=ψLE,where ψL indicates the mode shape normalized to the unit length (Euclidean length equals to unity).It is easy to realize that this procedure will define the modal mass estimated from equation (13) as something that is heavily dependent upon the number of DOFs in the model. In fact, it is easy to show that, for this case, the modal mass is approximately proportional to(17)m∝1N.Another approach to scaling is to scale the largest component (or some other component) to a certain value (for instance, equal to unity) [1–7]. This introduces the DOF scaled mode shape ψD:(18)ψD=ψψDn,where ψDn is the considered DOF used for scaling.If the mode shapes are normalized to the largest component equal to unity, the DOF used for normalization should be the same in the FE model and the test, i.e., the DOF must be shared. If this is not the case, then the modal mass is dependent upon the selected DOF.


## 1.5. Physical Interpretation of Modal Mass

Although the concept of modal mass is defined in the classical books of structural dynamics and modal analysis, the physical interpretation of this modal parameter is rarely addressed.The modal mass for the undamped and the proportionally damped cases is usually defined as a scaling factor for the mode shapes [1, 3, 11], i.e., it is used to convert the original unscaled mode shape vector {ψ} to the more useful mass-normalized mode shape vector {ϕ} [1–4].When using modal superposition (see equation (3)), the modal mass represents the mass of the single DOF systems, after decoupling the equations by transformation to modal coordinates [1–4].In [44], it is demonstrated that the modal mass of a simply-supported beam, when the mode shapes are normalized to the largest component equal to unity, is the half of the total mass of the system, whereas in fixed-fixed beams, the modal mass is equal to the total mass of the system.However, it seems like none of the classical ways of dealing with the scaling problem of the mode shapes leads to a modal mass that has the needed clarity in terms of physical meaning. Furthermore, since the results are dependent upon the number of DOFs in the model, the individual quantities do not seem to have much meaning at all.In this paper, a new and better definition of the modal mass is presented. This definition leads to significant improvements in the physical understanding of this modal parameter. Firstly, the most important requirements for a new definition of the modal mass are proposed. Then, a new expression for the length of a continuous mode shape is formulated introducing the concept of the volume matrix, which again leads to a new definition of the modal mass for the continuous case. The expressions derived for continuous systems are later extended to discrete cases. Finally, the concepts and equations formulated in the paper are illustrated by simple examples.


## 2. Basic Requirements

In order to achieve a better understanding of modal mass, one would expect that a new definition would lead to a better physical meaning of the modal mass, for instance, that the physical unit is indeedkg and that the size of the modal mass is a measure of the amount of mass moving in a given mode.Similarly, for the length of the mode shapes, since the mode shape is pure geometry, it is natural to expect that a good measure of the length is a pure geometrical measure, independent of the mass distribution in the considered body. It is also expected that since a “shape” normally means a continuum of dimensionless numbers, the length is dimensionless and a measure of the average movement over the mode shape.For practical applications (and for the physical meaning), it is important that the ways to estimate modal mass and mode shape length are objective so that different people will arrive at the same number for a given mode and a given structure. At least, they should be able to agree on some simple rules for obtaining these quantities so that they independently would arrive at the same estimates.At this point, we will conclude that the most important requirements for a new definition of the modal mass are as follows:(i)
The modal mass must be physically meaningful, having the unit ofkg and being a measure of the amount of mass moving in a given mode(ii)
The length of the mode shape must be a pure geometrical quantity, describing the average movement of the considered mode(iii)
Estimation of modal mass and mode shape length must be objective so that estimates obtained by different people are equal


## 3. Basic Concepts

We will introduce basic concepts considering the formulation in the continuous case where a mode shapeψ=ψr is given as a function of the position r=x,y,zT. The mode shape at a given point r has the three mode shape components ψx,ψy,ψz , resulting in the magnitude(19)ψ=ψr=ψx2r+ψy2r+ψz2r.The modal mass for the continuous case is now naturally defined as(20)m=∫Vρψ2dV,where integration is performed over the volume V of the considered structure and ρ=ρr is the mass density. Similarly, the total mass is defined as(21)M=∫VρdV.


### 3.1. Mode Shape Length and Modal Mass

For the continuous case, we will define the squared lengthL2 of the mode shape ψ as the average of the squared length ψ2 of the mode shape over the considered volume V:(22)L2=1V∫Vψ2dV.This secures that the length definition has the same unit as the mode shape. Thus, if the mode shape is dimensionless, so is the length. If the mass densityρ is constant, then from equation (20), we have(23)m=ρV1V∫Vψ2dV=ML2,which is the central equation of this paper stating the following for the case of constant mass density.The modal mass is equal to the product between the total mass of the structure and the squared length.If the considered body is constituted by two parts with two volumes,V1 with the mass density ρ1 and V2 with the mass density ρ2, then the corresponding result is(24)m=V1V1∫V1ρ1ψV12dV+V2V2∫V2ρ2ψV22dV=M1L12+M2L22,where ψV2 and ψV2 indicate the DOFs of the mode shape ψ in volumes V1 and V2, respectively, and(25)L12=1V1∫V1ψV12dV,L22=1V2∫V2ψV22dVare the partial lengths now defined over the partial volumes V1 and V2, respectively.Equation (20) easily generalizes to cases with many volumes with constant mass density.As indicated above, the total mass can either be calculated as the integral of the mass density over the total volume, equation (21), or as the two similar integrals over the two partial volumes V1 and V2:(26)M=∫V1ρ1dV+∫V2ρ2dV=M1+M2.


### 3.2. Moving Mass

Now, assuming that the mode shapes vanish insideV2, ψr2=0, for r∈V2, and using equation (24) with L22=0, we get(27)m=M1L12,which is obviously the right answer.Using equation (23) on the total volume, we get(28)m=ML2=MVVL2=ρ1V1+ρ2V2V1+V2VL2.Due to the fact that(29)L2=1V∫Vψ2dV=1VV1V1∫V1ψV12dV+V2V2∫V2ψV22dV,the following relationship between the total length and the partial lengths is derived:(30)L2V=L12V1+L22V2.IfL22=0, equation (30) results in(31)L2V=L12V1,and equation (28) can also be expressed as(32)m=ρ1V1+ρ2V2V1+V2V1L12.Equation (32) is equal to equation (23) only if ρ1=ρ2 or V2=0. This illustrates the validity of equations (22)–(24).Ifρ1=ρ2=ρ, then from equation (32), we have that m=ML2=M1L12, i.e., we still get the same result for the modal mass. So, adding the dead mass M2 does not change the modal mass because L22=0.However, since we are looking for a unique relationship between the modal mass and the length, we want the total moving mass to be a well-defined quantity. We can find that moving massM1 from equation (23), by defining the mode shape ψU, as a unitary translation inside V1, that is, for ψr2=1,r∈V1 because then using equation (23),(33)m=ρV1V∫VψU2dV=ρ∫V1dV=M1.So, if the total moving mass is smaller than the total mass in the considered system, dead mass is present.The moving mass is a well-defined quantity, and equation (23) defines a unique relation between the modal moving mass and the mode shape length.


### 3.3. Apparent Mass

In the general case of a varying mass density, we can define an average mass density as(34)ρ=MV=1V∫VρdVso that equation (23) still holds, but where the length is now defined as(35)L2=1V∫Vψ2rdV,where r=ρ/ρ, and the total mass is given in terms of the average density M=ρV. We shall not use this definition of the length, the reason being that the quantity contains information about the mass distribution, and thus, as defined in equation (35), is not a pure geometrical measure.Instead, we will start from equation (24) formulated for an arbitrary number of constant mass distributions(36)m=∑nMnLn2.From equations (24) and (26), it is clear that the total mass is M=∑nMn and that the total length is(37)L2=1V∑nVnLn2so that from equation (36), we have(38)m=MaL2,where the apparent total mass Ma is given by(39)Ma=V∑nMnLn2∑nVnLn2.The advantage of using the definition of the length given by equation (37) is that the length remains a pure geometrical quantity. Unfortunately, this introduces a new mass quantity, namely, the apparent total mass Ma, which for nonuniform mass distribution is different from the total mass and changes value from mode shape to mode shape.From equation (39), it is clear that if and only if the mass density is uniform, then Ma=M. Thus, the deviation of the apparent total mass Ma from the total mass M is a measure of how much the mass distribution deviates from the uniform distribution.


### 3.4. Discrete Formulations

The discrete formulation follows directly from the continuous formulation replacing the integrals with proper similar summations. Similar to equation (22), the squared length L2 of the mode shape ψ over a discrete number N of nodal points in the 3D space is then(40)L2=1V∑nψn2ΔVn,where ψn2 is the length of the mode shape ψ=ψrn, where the vector rn is pointing to the nodal point n representing the volume ΔVn. Inserting the volumes ΔVn into the diagonal matrix V and, similarly, collecting the DOF of the mode shape in the mode shape vector, equation (36) can be written as(41)L2=1VψTVψ.We should be aware that we are, here, using the same symbol “ψ” for the mode shape as before, even though before ψ was vector field in the 3D space but now it is actually a discrete vector containing the DOFs of the mode shape. The discrete counterpart to equation (23) is now(42)m=ρV1VψTVψ=ML2,which, as before, is only valid for a constant mass density ρ. Comparing this equation to the classical equation for the modal mass given by equation (13) naturally extends the definition of the volume matrix to the general nondiagonal form(43)V=Mρ,where V and M are the volume and mass matrices of the system, respectively. In this more general understanding of the volume matrix, the matrix can be full and can content rotational degrees of freedom. This also means that the length definition given by equation (37) naturally extends to cases with rotational degrees of freedom.For two different volumes,V1 with the mass density ρ1 and V2 with the mass density ρ2, equation (24) still holds but now with the length definitions:(44)L12=1V1ψ1TVψ1,L22=1V2ψ2TVψ2.For the general case of a varying mass density, we find the discrete expression for the length from equation (24):(45)L2=1V∫Vψ2ρdV,which directly leads to(46)L2=1VψTMR−1ψ.If the mass matrixM is diagonal, then R is a diagonal matrix containing the elements ρn representative for the volume ΔVn. Since M is symmetric, MR−1 can be replaced by R−1M or by the symmetric form R−1/2MR1/2. At this moment, it is interesting to note that the symmetric form R−1/2MR1/2 naturally defines a new definition of a mass scaled mode shape ψm=R−1/2ψ, but since it leads to mode shapes with the unfortunate unit 1/kg, we will not follow this lead any further.Equation (46) naturally suggests the general volume matrix:(47)V=MR−1,which is not necessarily diagonal. For instance, using an FE model, with the value ρn of the mass density in a certain element with the local mass matrix Mn, the corresponding local volume matrix Vn is defined according to equation (43):(48)Vn=Mn/ρn,and the global volume matrix is then assembled by the same procedure as when assembling the local mass matrices to form the global mass matrix. Following this procedure, equations (32)–(35) still hold in the general case.Before we leave this section, let us conclude that the two important scalar quantities, the total massM and the corresponding total volume V, can be found from any discrete converged model as(49)M=ψUTMψU,V=ψUTVψU,where ψU is a unitary translation vector as mentioned earlier. We could also find the total mass and total volume adding all matrix elements corresponding to translational degrees of freedom (but then, we have to take into account if 2 or 3D formulations are being used dividing the total sum by 2 or 3, respectively).


## 3.1. Mode Shape Length and Modal Mass

For the continuous case, we will define the squared lengthL2 of the mode shape ψ as the average of the squared length ψ2 of the mode shape over the considered volume V:(22)L2=1V∫Vψ2dV.This secures that the length definition has the same unit as the mode shape. Thus, if the mode shape is dimensionless, so is the length. If the mass densityρ is constant, then from equation (20), we have(23)m=ρV1V∫Vψ2dV=ML2,which is the central equation of this paper stating the following for the case of constant mass density.The modal mass is equal to the product between the total mass of the structure and the squared length.If the considered body is constituted by two parts with two volumes,V1 with the mass density ρ1 and V2 with the mass density ρ2, then the corresponding result is(24)m=V1V1∫V1ρ1ψV12dV+V2V2∫V2ρ2ψV22dV=M1L12+M2L22,where ψV2 and ψV2 indicate the DOFs of the mode shape ψ in volumes V1 and V2, respectively, and(25)L12=1V1∫V1ψV12dV,L22=1V2∫V2ψV22dVare the partial lengths now defined over the partial volumes V1 and V2, respectively.Equation (20) easily generalizes to cases with many volumes with constant mass density.As indicated above, the total mass can either be calculated as the integral of the mass density over the total volume, equation (21), or as the two similar integrals over the two partial volumes V1 and V2:(26)M=∫V1ρ1dV+∫V2ρ2dV=M1+M2.


## 3.2. Moving Mass

Now, assuming that the mode shapes vanish insideV2, ψr2=0, for r∈V2, and using equation (24) with L22=0, we get(27)m=M1L12,which is obviously the right answer.Using equation (23) on the total volume, we get(28)m=ML2=MVVL2=ρ1V1+ρ2V2V1+V2VL2.Due to the fact that(29)L2=1V∫Vψ2dV=1VV1V1∫V1ψV12dV+V2V2∫V2ψV22dV,the following relationship between the total length and the partial lengths is derived:(30)L2V=L12V1+L22V2.IfL22=0, equation (30) results in(31)L2V=L12V1,and equation (28) can also be expressed as(32)m=ρ1V1+ρ2V2V1+V2V1L12.Equation (32) is equal to equation (23) only if ρ1=ρ2 or V2=0. This illustrates the validity of equations (22)–(24).Ifρ1=ρ2=ρ, then from equation (32), we have that m=ML2=M1L12, i.e., we still get the same result for the modal mass. So, adding the dead mass M2 does not change the modal mass because L22=0.However, since we are looking for a unique relationship between the modal mass and the length, we want the total moving mass to be a well-defined quantity. We can find that moving massM1 from equation (23), by defining the mode shape ψU, as a unitary translation inside V1, that is, for ψr2=1,r∈V1 because then using equation (23),(33)m=ρV1V∫VψU2dV=ρ∫V1dV=M1.So, if the total moving mass is smaller than the total mass in the considered system, dead mass is present.The moving mass is a well-defined quantity, and equation (23) defines a unique relation between the modal moving mass and the mode shape length.


## 3.3. Apparent Mass

In the general case of a varying mass density, we can define an average mass density as(34)ρ=MV=1V∫VρdVso that equation (23) still holds, but where the length is now defined as(35)L2=1V∫Vψ2rdV,where r=ρ/ρ, and the total mass is given in terms of the average density M=ρV. We shall not use this definition of the length, the reason being that the quantity contains information about the mass distribution, and thus, as defined in equation (35), is not a pure geometrical measure.Instead, we will start from equation (24) formulated for an arbitrary number of constant mass distributions(36)m=∑nMnLn2.From equations (24) and (26), it is clear that the total mass is M=∑nMn and that the total length is(37)L2=1V∑nVnLn2so that from equation (36), we have(38)m=MaL2,where the apparent total mass Ma is given by(39)Ma=V∑nMnLn2∑nVnLn2.The advantage of using the definition of the length given by equation (37) is that the length remains a pure geometrical quantity. Unfortunately, this introduces a new mass quantity, namely, the apparent total mass Ma, which for nonuniform mass distribution is different from the total mass and changes value from mode shape to mode shape.From equation (39), it is clear that if and only if the mass density is uniform, then Ma=M. Thus, the deviation of the apparent total mass Ma from the total mass M is a measure of how much the mass distribution deviates from the uniform distribution.


## 3.4. Discrete Formulations

The discrete formulation follows directly from the continuous formulation replacing the integrals with proper similar summations. Similar to equation (22), the squared length L2 of the mode shape ψ over a discrete number N of nodal points in the 3D space is then(40)L2=1V∑nψn2ΔVn,where ψn2 is the length of the mode shape ψ=ψrn, where the vector rn is pointing to the nodal point n representing the volume ΔVn. Inserting the volumes ΔVn into the diagonal matrix V and, similarly, collecting the DOF of the mode shape in the mode shape vector, equation (36) can be written as(41)L2=1VψTVψ.We should be aware that we are, here, using the same symbol “ψ” for the mode shape as before, even though before ψ was vector field in the 3D space but now it is actually a discrete vector containing the DOFs of the mode shape. The discrete counterpart to equation (23) is now(42)m=ρV1VψTVψ=ML2,which, as before, is only valid for a constant mass density ρ. Comparing this equation to the classical equation for the modal mass given by equation (13) naturally extends the definition of the volume matrix to the general nondiagonal form(43)V=Mρ,where V and M are the volume and mass matrices of the system, respectively. In this more general understanding of the volume matrix, the matrix can be full and can content rotational degrees of freedom. This also means that the length definition given by equation (37) naturally extends to cases with rotational degrees of freedom.For two different volumes,V1 with the mass density ρ1 and V2 with the mass density ρ2, equation (24) still holds but now with the length definitions:(44)L12=1V1ψ1TVψ1,L22=1V2ψ2TVψ2.For the general case of a varying mass density, we find the discrete expression for the length from equation (24):(45)L2=1V∫Vψ2ρdV,which directly leads to(46)L2=1VψTMR−1ψ.If the mass matrixM is diagonal, then R is a diagonal matrix containing the elements ρn representative for the volume ΔVn. Since M is symmetric, MR−1 can be replaced by R−1M or by the symmetric form R−1/2MR1/2. At this moment, it is interesting to note that the symmetric form R−1/2MR1/2 naturally defines a new definition of a mass scaled mode shape ψm=R−1/2ψ, but since it leads to mode shapes with the unfortunate unit 1/kg, we will not follow this lead any further.Equation (46) naturally suggests the general volume matrix:(47)V=MR−1,which is not necessarily diagonal. For instance, using an FE model, with the value ρn of the mass density in a certain element with the local mass matrix Mn, the corresponding local volume matrix Vn is defined according to equation (43):(48)Vn=Mn/ρn,and the global volume matrix is then assembled by the same procedure as when assembling the local mass matrices to form the global mass matrix. Following this procedure, equations (32)–(35) still hold in the general case.Before we leave this section, let us conclude that the two important scalar quantities, the total massM and the corresponding total volume V, can be found from any discrete converged model as(49)M=ψUTMψU,V=ψUTVψU,where ψU is a unitary translation vector as mentioned earlier. We could also find the total mass and total volume adding all matrix elements corresponding to translational degrees of freedom (but then, we have to take into account if 2 or 3D formulations are being used dividing the total sum by 2 or 3, respectively).


## 4. SEREP Reduction

In this section, we will show that the so-defined modal mass and modal length are invariant to the SEREP (system equivalent reduction expansion process) reduction technique [45]. For this case, it is practical to consider a mode shape cluster Ψ instead of the single mode shape. We can then define the diagonal modal mass matrix(50)mn=ΨTMΨ.Using SEREP, we have the reduced mass matrix(51)Ma=TTMT,where T=ΨaΨa+ (and where we assume full rank, i.e., Ψa is Na×Na). The modal mass matrix for the reduced mass matrix is(52)man=ΨaTMaΨa=ΨaTΨaΨa+TMΨaΨa+Ψa.However, sinceΨaΨa+=I, from equation (52), it is derived that(53)man=mn.Using the volume matrixV from equation (41), we can define the modal length matrix(54)L2=1VΨTVΨ.This matrix is clearly nondiagonal but positive definite so that the notationL2 makes sense. Now, defining the SEREP-reduced volume matrix(55)Va=TTVT,we can calculate the length in the reduced system:(56)La2=1VΨaTVaΨa.And following the same procedure as for the modal mass,(57)La2=1VΨaTTTVTΨa=1VΨTVΨ=L2,or in words, Modal mass and length matrices are invariant to SEREP.We should note that since the modal mass and modal length are the same in all reduced systems when using SEREP [45], equation (19) still holds in all reduced systems if we use the simple rule that the total mass should be estimated from the original system. The same is true for the total volume.This means that the two key quantities, the total mass and the total volume, have to be estimated for the considered system once and for all. One can use an experimental approach and use a weight scale to measure the total weight and a volume scale, for instance, a fluid expansion method, to measure the volume. This is of course the best possible way to estimate these two key quantities.If we cannot use the experimental way to estimate these quantities, then we could use an analytical model or a numerical model. An analytical model is often the simplest and the best approach. For example, if we consider a beam with lengthL and constant area A made out of material with the constant mass density ρ, then the total mass is simply M=ρAL, and the total volume is simply V=AL.If we want to use a finite element model to estimate the quantities, then like for all other quantities, the model must have converged in order to make sure that the estimated quantities are reliable.


## 5. Mode Shape Scaling

We have now established the important linear relation between the total mass of the considered system, the modal mass of the considered mode shape, and the length of the considered mode shape in case of uniform mass density:(58)m=ML2.And we will now consider suitable rules for how to choose the value forL. When a meaningful value for L has been established, the modal mass is then totally defined by equation (58). As explained in Section 1, we have to at least consider the classical scaling, the length scaling, and the DOF scaling.


### 5.1. Length Scaling

Considering the arbitrarily scaled mode shapeψ with the length L, then a classical length scaling of ψ would be defining the scaled mode shape ψL as(59)ψL=ψL.Using equation (41), we find the resulting length of ψL as(60)L2=1VψLTVψL=1L21VψTVψ=1.Thus, the classical length scaling leads to a version of equation (55), where(61)m=M,or in words,Using length scaling, the modal mass in constant mass density systems is always equal to the total mass of the considered system.One can say that this is a simple rule and at least easy to remember and to use and also practical because, now, different people can easily arrive to the same number (or at least similar estimates) because the art is only to estimate the total mass of the system.The drawback is of course that this modal mass does not give us any information about how much mass is actually moving in the different modes. If mass density is constant, all the modes have the same modal mass. Even a local mode, where nearly all movements are localized to a very small region, still has the same modal mass. This is unfortunate, but if one can live with this limitation, it might be a good way to deal with the modal mass.In nonconstant mass density systems, the modal mass is different for each mode and equal to the apparent mass.


### 5.2. DOF Scaling

The classical DOF scaling is based on choosing one of the DOFs in the system and then scale this value to a certain number, for instance, to unity [1–4].Before we go further into this scaling, let us consider an abstract example. Like in Section3, let us consider the case where the continuous body is constituted by two parts with two volumes, V1 and V2, and let us assume that we consider the abstract mode shape ψU as an unitary translation inside V1 and zero inside V2.(62)ψUr2=1;r∈V1,0;r∈V2.Just for simplicity, let us assume that the mass density is constant, so from equation (20),(63)m=∫V1ψU2ρdV=ρV1=M1.Thus, the modal mass is then equal to the total mass inside the volumeV1. It is also clear that if V1 extends to the total volume V=V1+V2, then the modal mass becomes equal to the total mass of the system.This case illustrates the DOF scaling, where the maximum value ofψUr2 is scaled to unity. Inspired by this case, we will define the DOF scaled version as(64)ψD=ψmaxψr.Using this scaling, the modal mass has the properties(65)m≤M,and thus is a direct measure of how much mass is moving in the considered mode shape. In words, we can express this asUsing largest component equal to unity, the modal mass in constant mass density systems is always smaller than or equal to the total mass of the considered system.The advantage of this scaling is just that the size of the modal mass now makes totally sense, but the drawback is that we need to scale the largest component to unity, which is the same as scaling:(66)maxψD=1.In a model with only a few components, where one of the DOFs is not directly equal to the DOF where the mode shape is maximum, it is not so easy to find the scaling. The solution is to know the ratio between one of the DOFs of the model and the valuemaxψD found in a refined model (ideally, a model where all quantities of importance have converged).


### 5.3. Scaling Using SEREP

We start out with an initial model where all important quantities have converged towards a stable value; then, it is easy to conclude from the earlier analysis that, for the length scaling, all reduced models satisfy equation (61), i.e., in all reduced models, the modal masses are the same and equal to the total mass of the system.Using DOF scaling, all common DOFs are still the same, and the length, which can be estimated in any of the systems reduced from the initial system, is the same; so, also, all modal masses are the same, and the property expressed in equation (65) is true in all the reduced systems. One should note that equation (66) only holds in systems which are reduced to a size where maxψ is the same as in the initial system. If this is not the case, then maxψ must be estimated in a more refined system.


## 5.1. Length Scaling

Considering the arbitrarily scaled mode shapeψ with the length L, then a classical length scaling of ψ would be defining the scaled mode shape ψL as(59)ψL=ψL.Using equation (41), we find the resulting length of ψL as(60)L2=1VψLTVψL=1L21VψTVψ=1.Thus, the classical length scaling leads to a version of equation (55), where(61)m=M,or in words,Using length scaling, the modal mass in constant mass density systems is always equal to the total mass of the considered system.One can say that this is a simple rule and at least easy to remember and to use and also practical because, now, different people can easily arrive to the same number (or at least similar estimates) because the art is only to estimate the total mass of the system.The drawback is of course that this modal mass does not give us any information about how much mass is actually moving in the different modes. If mass density is constant, all the modes have the same modal mass. Even a local mode, where nearly all movements are localized to a very small region, still has the same modal mass. This is unfortunate, but if one can live with this limitation, it might be a good way to deal with the modal mass.In nonconstant mass density systems, the modal mass is different for each mode and equal to the apparent mass.


## 5.2. DOF Scaling

The classical DOF scaling is based on choosing one of the DOFs in the system and then scale this value to a certain number, for instance, to unity [1–4].Before we go further into this scaling, let us consider an abstract example. Like in Section3, let us consider the case where the continuous body is constituted by two parts with two volumes, V1 and V2, and let us assume that we consider the abstract mode shape ψU as an unitary translation inside V1 and zero inside V2.(62)ψUr2=1;r∈V1,0;r∈V2.Just for simplicity, let us assume that the mass density is constant, so from equation (20),(63)m=∫V1ψU2ρdV=ρV1=M1.Thus, the modal mass is then equal to the total mass inside the volumeV1. It is also clear that if V1 extends to the total volume V=V1+V2, then the modal mass becomes equal to the total mass of the system.This case illustrates the DOF scaling, where the maximum value ofψUr2 is scaled to unity. Inspired by this case, we will define the DOF scaled version as(64)ψD=ψmaxψr.Using this scaling, the modal mass has the properties(65)m≤M,and thus is a direct measure of how much mass is moving in the considered mode shape. In words, we can express this asUsing largest component equal to unity, the modal mass in constant mass density systems is always smaller than or equal to the total mass of the considered system.The advantage of this scaling is just that the size of the modal mass now makes totally sense, but the drawback is that we need to scale the largest component to unity, which is the same as scaling:(66)maxψD=1.In a model with only a few components, where one of the DOFs is not directly equal to the DOF where the mode shape is maximum, it is not so easy to find the scaling. The solution is to know the ratio between one of the DOFs of the model and the valuemaxψD found in a refined model (ideally, a model where all quantities of importance have converged).


## 5.3. Scaling Using SEREP

We start out with an initial model where all important quantities have converged towards a stable value; then, it is easy to conclude from the earlier analysis that, for the length scaling, all reduced models satisfy equation (61), i.e., in all reduced models, the modal masses are the same and equal to the total mass of the system.Using DOF scaling, all common DOFs are still the same, and the length, which can be estimated in any of the systems reduced from the initial system, is the same; so, also, all modal masses are the same, and the property expressed in equation (65) is true in all the reduced systems. One should note that equation (66) only holds in systems which are reduced to a size where maxψ is the same as in the initial system. If this is not the case, then maxψ must be estimated in a more refined system.


## 6. Examples

### 6.1. A Cantilever Beam with Constant Density

Several simulations were carried out on a planar finite element model of a cantilever beam with constant mass density and constant cross section (Figure1). The geometric and material properties of the cross section are also shown in Figure 1.Figure 1
Cantilever beam with constant mass density.The analytical mode shapes corresponding to the first three modes and normalized to the largest component equal to unity are shown in Figure2 (solid lines). On the contrary, the analytical modal masses and the squared length of the mode shapes are presented in Table 1. In this particular case, the largest component is at the free border of the beam for all the modes, and the modal masses are equal for all the modes and so is the length of the mode shapes.Figure 2
First three mode shapes normalized to the largest component equal to unity. Constant mass density (solid lines). Nonconstant mass density (dashed lines).Table 1
Modal masses, effective masses, and length of mode shapes for the cantilever beams.
ModelConstant density (Figure1)Nonconstant density (Figure5)Total mass (kg)1616Apparent mass (kg)Equal to total mass,Ma1=Ma2=Ma3=16Ma1=14.10Ma2=15.93Ma3=16.06Squared length of mode shapesL12=L22=L32=0.25L12=0.2517L22=0.2049L32=0.2063Modal mass (kg)m1=m2=m3=4m1=3.549m2=3.265m3=3.313In order to study the convergence of the equations proposed in this paper, the finite element model was discretized with different number of linear bending finite elements (Figure1) ranging from 2 to 100. Moreover, both consistent and lumped mass matrices were considered in the simulations.The modal masses estimated with equation (13) for the first three modes using different number of elements are presented in Figure 3. It can be observed that the estimated modal masses converge very fast to the analytical solution when a consistent mass matrix is used, the error being less than 1% for all the modes when the model is discretized with 6 elements. When a lumped mass matrix is considered, the convergence is achieved with more elements (error less than 1% with 9, 23, and 36 elements for the first, second, and third modes, respectively). Due to the fact that the mass density is constant, the ratio m/L2 is constant for all the modes, and the same conclusions can be inferred for the length of the mode shapes (Figure 4).Figure 3
Convergence of the modal masses for the first three modes of the cantilever beam with constant mass density.Figure 4
Convergence of the length of the mode shapes for the first three modes of the cantilever beam with constant mass density.Table2 presents the results obtained with a reduced model using SEREP when only three DOFs are considered. As expected from the theory formulated in Section 6, the modal masses and the length of the mode shapes of the reduced model are exactly the same as those of the full model, independent of the number and the distribution of the active DOFs. The reduced mass and volume matrices are also shown in Table 2.Table 2
Results obtained with a reduced model using SEREP [45].
Active DOFs[33, 65, 99][49, 75, 99]Mass matrix (kg)18681−129414336−129418972−30094336−300910177083−86163735−861610493−45543735−45541985Volume matrixm32.3798−1.64850.5523−1.64851.1430−0.38340.5523−0.38340.12960.9023−1.09760.4578−1.09761.3367−0.58010.4578−0.58010.2528Length of mode shapesL12=L22=L32=0.25Modal mass (kg)m1=m2=m3=4Figure 5
Cantilever beam with nonconstant mass density.


### 6.2. A Cantilever Beam with Nonconstant Density

Simulations were also performed on the same cantilever beam but with nonconstant mass density. The geometric and material properties of the cross section are shown in Figure5.The analytical mode shapes corresponding to the first three modes and normalized to the largest component equal to unity are shown in Figure2 (dashed lines). It can be observed that the first mode shape is very similar to that corresponding to the constant mass density case, but significant differences can be seen in the second and third modes.The analytical modal masses and the squared length of the mode shapes are presented in Table1. As the mass density is not constant, equation (25) does not hold, and an apparent mass can be calculated for each mode (see Table 1).The finite element model was discretized with different number of linear bending elements ranging from 4 to 100. Both consistent and lumped mass matrices were considered in the simulations.The modal masses estimated with equation (14) using different number of elements are presented in Figure 6. It can be observed that, again, the estimated modal masses converge very fast to the analytical solution when a consistent mass matrix is used, the error being less than 1% for all the modes with 7 elements. When a lumped mass matrix is considered, the convergence is achieved with more elements (error less than 1% with 12, 24, and 36 elements for the first, second, and third modes, respectively).Figure 6
Convergence of the modal masses for the first three modes of the cantilever beam with nonconstant mass density.With respect to the length of the mode shapes (Figure7), the same speed of convergence as that observed for the constant mass density case was obtained for the nonconstant density case (error less than 1% with 6 elements and consistent mass matrix).Figure 7
Convergence of the length of the mode shapes for the first three modes of the cantilever beam with nonconstant mass density.


### 6.3. A Rigid Block on Springs

In order to illustrate how the equations proposed in this paper work on discrete systems modelled with lumped masses and lumped inertias, the modal masses and mode shape lengths of a rigid block on two springs (see Tables3–5) vibrating in the x-y plane (bouncing mode and pitch mode) have been calculated. We are illustrating the equations corresponding to both the continuous and the discrete case.Table 3
A rigid block on springs. Continuous formulas.
Mode 1Mode 2Mode shapesψ1x=0ψ1y=1ψ2y=2x/aψ2x=−2y/bLength of mode shapes,L2=1V∫Vψx2+ψy2dVL12=1L22=1/31+b2/a2Modal mass,m=∫Vρψx2+ψy2dVm1=Mm2=M/31+b2/a2Table 4
A rigid block on springs. Discretization with two translational DOFs.
Mode 1Mode 2Mode shapesψ1=11ψ2=−11Mass matrixMM/4+J/a2M/4−J/a2M/4−J/a2M/4+J/a2Volume matrixVV/4+V/121+b2/a2V/4−V/121+b2/a2V/4−V/121+b2/a2V/4+V/121+b2/a2Length of mode shapes,L2=1/VψTVψL12=1L22=1/31+b2/a2Modal mass,m=ψTMψm1=Mm2=M/31+b2/a2Table 5
A rigid block on springs. Discretization with one translation and one rotation.
Mode 1Mode shapesMode shapesψ1=10ψ2=02/aMass matrixMM00JVolume matrixVV00V/121+b2/a2Length of mode shapes,L2=1/VψTVψL12=1L22=1/31+b2/a2Modal mass,m=ψTMψm1=Mm2=M/31+b2/a2The block has constant densityρ, dimensions a×b×c, total mass M=ρabc, total volume V=M/ρ, and inertia J=M/12a2+b2 with respect to the center of gravity of the block.If we consider the block as a continuous system, the analytical expressions of the mode shapes in both directions,x and y, are needed (see Table 3). The mode shapes have been normalized to the largest vertical component equal to unity. The modal masses and the length of the mode shapes calculated with equations (16) and (18) are presented in Table 3.If the system is modelled with two translational DOFs (see Table4), the mass matrix, the volume matrices, and the mode shapes are presented in Table 4. The modal masses and the mode shape lengths calculated with equations (37) and (38), respectively, are shown in Table 4.The results obtained with the same system modelled with one translation and one rotational DOF are shown in Table5.In can be seen from Tables3–5 that the obtained values of the mode shape lengths and the modal masses are the same using different models. Moreover, as the density is constant, the modal mass is equal to the product between the total mass of the structure and the squared length. This example demonstrates that the concepts and equations formulated in this paper can also be easily extended to discrete systems modelled with lumped masses and/or lumped inertias.


## 6.1. A Cantilever Beam with Constant Density

Several simulations were carried out on a planar finite element model of a cantilever beam with constant mass density and constant cross section (Figure1). The geometric and material properties of the cross section are also shown in Figure 1.Figure 1
Cantilever beam with constant mass density.The analytical mode shapes corresponding to the first three modes and normalized to the largest component equal to unity are shown in Figure2 (solid lines). On the contrary, the analytical modal masses and the squared length of the mode shapes are presented in Table 1. In this particular case, the largest component is at the free border of the beam for all the modes, and the modal masses are equal for all the modes and so is the length of the mode shapes.Figure 2
First three mode shapes normalized to the largest component equal to unity. Constant mass density (solid lines). Nonconstant mass density (dashed lines).Table 1
Modal masses, effective masses, and length of mode shapes for the cantilever beams.
ModelConstant density (Figure1)Nonconstant density (Figure5)Total mass (kg)1616Apparent mass (kg)Equal to total mass,Ma1=Ma2=Ma3=16Ma1=14.10Ma2=15.93Ma3=16.06Squared length of mode shapesL12=L22=L32=0.25L12=0.2517L22=0.2049L32=0.2063Modal mass (kg)m1=m2=m3=4m1=3.549m2=3.265m3=3.313In order to study the convergence of the equations proposed in this paper, the finite element model was discretized with different number of linear bending finite elements (Figure1) ranging from 2 to 100. Moreover, both consistent and lumped mass matrices were considered in the simulations.The modal masses estimated with equation (13) for the first three modes using different number of elements are presented in Figure 3. It can be observed that the estimated modal masses converge very fast to the analytical solution when a consistent mass matrix is used, the error being less than 1% for all the modes when the model is discretized with 6 elements. When a lumped mass matrix is considered, the convergence is achieved with more elements (error less than 1% with 9, 23, and 36 elements for the first, second, and third modes, respectively). Due to the fact that the mass density is constant, the ratio m/L2 is constant for all the modes, and the same conclusions can be inferred for the length of the mode shapes (Figure 4).Figure 3
Convergence of the modal masses for the first three modes of the cantilever beam with constant mass density.Figure 4
Convergence of the length of the mode shapes for the first three modes of the cantilever beam with constant mass density.Table2 presents the results obtained with a reduced model using SEREP when only three DOFs are considered. As expected from the theory formulated in Section 6, the modal masses and the length of the mode shapes of the reduced model are exactly the same as those of the full model, independent of the number and the distribution of the active DOFs. The reduced mass and volume matrices are also shown in Table 2.Table 2
Results obtained with a reduced model using SEREP [45].
Active DOFs[33, 65, 99][49, 75, 99]Mass matrix (kg)18681−129414336−129418972−30094336−300910177083−86163735−861610493−45543735−45541985Volume matrixm32.3798−1.64850.5523−1.64851.1430−0.38340.5523−0.38340.12960.9023−1.09760.4578−1.09761.3367−0.58010.4578−0.58010.2528Length of mode shapesL12=L22=L32=0.25Modal mass (kg)m1=m2=m3=4Figure 5
Cantilever beam with nonconstant mass density.


## 6.2. A Cantilever Beam with Nonconstant Density

Simulations were also performed on the same cantilever beam but with nonconstant mass density. The geometric and material properties of the cross section are shown in Figure5.The analytical mode shapes corresponding to the first three modes and normalized to the largest component equal to unity are shown in Figure2 (dashed lines). It can be observed that the first mode shape is very similar to that corresponding to the constant mass density case, but significant differences can be seen in the second and third modes.The analytical modal masses and the squared length of the mode shapes are presented in Table1. As the mass density is not constant, equation (25) does not hold, and an apparent mass can be calculated for each mode (see Table 1).The finite element model was discretized with different number of linear bending elements ranging from 4 to 100. Both consistent and lumped mass matrices were considered in the simulations.The modal masses estimated with equation (14) using different number of elements are presented in Figure 6. It can be observed that, again, the estimated modal masses converge very fast to the analytical solution when a consistent mass matrix is used, the error being less than 1% for all the modes with 7 elements. When a lumped mass matrix is considered, the convergence is achieved with more elements (error less than 1% with 12, 24, and 36 elements for the first, second, and third modes, respectively).Figure 6
Convergence of the modal masses for the first three modes of the cantilever beam with nonconstant mass density.With respect to the length of the mode shapes (Figure7), the same speed of convergence as that observed for the constant mass density case was obtained for the nonconstant density case (error less than 1% with 6 elements and consistent mass matrix).Figure 7
Convergence of the length of the mode shapes for the first three modes of the cantilever beam with nonconstant mass density.


## 6.3. A Rigid Block on Springs

In order to illustrate how the equations proposed in this paper work on discrete systems modelled with lumped masses and lumped inertias, the modal masses and mode shape lengths of a rigid block on two springs (see Tables3–5) vibrating in the x-y plane (bouncing mode and pitch mode) have been calculated. We are illustrating the equations corresponding to both the continuous and the discrete case.Table 3
A rigid block on springs. Continuous formulas.
Mode 1Mode 2Mode shapesψ1x=0ψ1y=1ψ2y=2x/aψ2x=−2y/bLength of mode shapes,L2=1V∫Vψx2+ψy2dVL12=1L22=1/31+b2/a2Modal mass,m=∫Vρψx2+ψy2dVm1=Mm2=M/31+b2/a2Table 4
A rigid block on springs. Discretization with two translational DOFs.
Mode 1Mode 2Mode shapesψ1=11ψ2=−11Mass matrixMM/4+J/a2M/4−J/a2M/4−J/a2M/4+J/a2Volume matrixVV/4+V/121+b2/a2V/4−V/121+b2/a2V/4−V/121+b2/a2V/4+V/121+b2/a2Length of mode shapes,L2=1/VψTVψL12=1L22=1/31+b2/a2Modal mass,m=ψTMψm1=Mm2=M/31+b2/a2Table 5
A rigid block on springs. Discretization with one translation and one rotation.
Mode 1Mode shapesMode shapesψ1=10ψ2=02/aMass matrixMM00JVolume matrixVV00V/121+b2/a2Length of mode shapes,L2=1/VψTVψL12=1L22=1/31+b2/a2Modal mass,m=ψTMψm1=Mm2=M/31+b2/a2The block has constant densityρ, dimensions a×b×c, total mass M=ρabc, total volume V=M/ρ, and inertia J=M/12a2+b2 with respect to the center of gravity of the block.If we consider the block as a continuous system, the analytical expressions of the mode shapes in both directions,x and y, are needed (see Table 3). The mode shapes have been normalized to the largest vertical component equal to unity. The modal masses and the length of the mode shapes calculated with equations (16) and (18) are presented in Table 3.If the system is modelled with two translational DOFs (see Table4), the mass matrix, the volume matrices, and the mode shapes are presented in Table 4. The modal masses and the mode shape lengths calculated with equations (37) and (38), respectively, are shown in Table 4.The results obtained with the same system modelled with one translation and one rotational DOF are shown in Table5.In can be seen from Tables3–5 that the obtained values of the mode shape lengths and the modal masses are the same using different models. Moreover, as the density is constant, the modal mass is equal to the product between the total mass of the structure and the squared length. This example demonstrates that the concepts and equations formulated in this paper can also be easily extended to discrete systems modelled with lumped masses and/or lumped inertias.


## 7. Discussion

Considering a certain body, first, one has to decide which parts of the body have constant mass density. Let us first assume that we can identify a number of different bodies with the volumesVn, where this is the case. Then, for each of these partial bodies, the corresponding modal mass and modal length are calculated according equations (9) and (37):(67)mn=ψnTMnψn,Ln2=1VnψnTVnψn,where ψn is the part of the mode shape ψ corresponding to the considered partial body Vn, mn and Ln2 are the contributions to the modal mass and the mode shape length in the body, and Vn is the volume of the body. Mn is the mass matrix, and Vn is the volume matrix of the considered partial body, respectively, where Vn is found from equation (44). This means that we can calculate separately the modal mass mn and the length Ln2 of each partial body Vn (for example, in a finite element model, we can calculate the modal mass and the length of each of the finite elements of the model) and then calculate the total modal mass and the total length as the sum of the contributions of each body.The total mass and the total volume of each partial body are then(68)Mn=ψUTMnψU,Vn=ψUTVnψU,where ψU is a unitary translation. For each partial body, the so-calculated quantities should satisfy equation (47):(69)mn=MnLn2,for all modes. If this is not the case, the apparent total masses, as introduced in equation (39), are not equal to the total mass; then, it means that the assumption of constant mass density inside the considered partial body is wrong, and the body should be subdivided into smaller parts. However, if equation (65) holds for all partial bodies and all modes, then we can conclude the final results for the total body:(70)m=∑nmn,M=∑nMn,V=∑nVn,L2=∑nLn2.If we are not working with well-defined parts each having its own constant mass density but rather a density that is varying smoothly over the body, then it might be more convenient to use all elements as partial bodies. In this case, we can follow the same procedure as mentioned above, but where the volume matrixVn is now simply found from equation (49).The final question still remains—how should we choose the scaling? The answer is that it does not matter because whatever scaling we use, a linear relation exists between squared length and the modal mass given by equation (39), and we can accurately calculate the relation for any of the considered modes.At this point, it is worth noticing that we can in a way forget about the dead mass problem because any influence is removed using equation (39). However, of course, if the length contributions from all modes are zero for a certain partial body, then it might be useful to remove that body from the description/modelling of the problem.Finally, coming back to the initial requirements mentioned in Section2, we can conclude that most of them have been fulfilled with the present proposal. Both modal mass and mode shape length have now reasonable units and a well-defined meaning. Also, it is easy for anybody to achieve the same result no matter what scaling has been used because the final answer is not the modal mass, and the final answer is all quantities need to be defined by equation (49). However, one can say that since the DOF scaling maxψD=1 is more meaningful than the length scaling and because the modal masses using DOF scaling clearly illustrate the amount of mass moving under the considered mode shape, it might, in many cases, be more meaningful to use this scaling than any other possible scaling.


## 8. Conclusions

The concept of length of a continuous mode shape has been defined in this paper. The new definition depends on the mode shape and how the volume is distributed in the structure. Then, the concept of length has been extended to discrete systems by introducing the concept of a volume matrix. This length definition has the same unit as the mode shape, but it does not coincide with the length of a vector in an Euclidean space that is normally used in linear algebra. We omit the definition from linear algebra because in that case, the length depends on the number of components of the vector.A new and better definition of the modal mass, which is physically meaningful and does not depend on the number of DOFs of a discrete model, has been formulated. It is demonstrated that if the mass density of the system is constant, then the modal mass is always equal to the product between the total mass of the structure and the squared length.If the mass density is not constant, the concept of apparent mass is proposed. This apparent mass is different for each mode and depends on the mode shape and how the mass is distributed in the structure. In these cases, the modal mass is always equal to the product between the apparent mass and the length squared.If a model is reduced with the SEREP reduction technique, it has been demonstrated that the modal mass and the length of the mode shapes are invariant to SEREP, i.e., the modal mass and the modal length are the same in all reduced systems and equal to the modal mass and length of the full model.In constant mass density systems, when mode shapes are normalized to the length, the modal mass is always equal to the total mass of the considered system. The drawback is that this modal mass does not give us any information about how much mass is actually moving in different modes. For example, local and global modes have the same modal mass, in spite of the fact that global modes are moving much more mass. In nonconstant mass density systems, the modal mass is different for each mode and equal to an apparent mass.When mode shapes are normalized to the largest component equal to unity, the modal mass is a measure of how much mass is moving in the considered mode. The drawback is that we need to scale the maximum component equal to unity, and in a model with only a few DOFs, where the DOF with the largest component is not measured, it is not easy to find the right scaling because we have to find the largest component in the converged model.



---

*Source: 8648769-2020-06-23.xml*