Identification of Eigen-Frequencies and Mode-Shapes of Beams with Continuous Distribution of Mass and Elasticity and for Various Conditions at Supports

In the present article, an equivalent three degrees of freedom (DoF) system of two different cases of inverted pendulums is presented for each separated case. The first case of inverted pendulum refers to an amphi-hinge pendulum that possesses distributed mass and stiffness along its height, while the second case of inverted pendulum refers to an inverted pendulum with distributed mass and stiffness along its height. These vertical pendulums have infinity number of degree of freedoms. Based on the free vibration of the above-mentioned pendulums according to partial differential equation, a mathematically equivalent three-degree of freedom system is given for each case, where its equivalent mass matrix is analytically formulated with reference on specific mass locations along the pendulum height. Using the three DoF model, the first three fundamental frequencies of the real pendulum can be identified with very good accuracy. Furthermore, taking account the 3 (cid:1) 3 mass matrix, it is possible to estimate the possible pendulum damages using a known technique of identification mode-shapes via records of response accelerations. Moreover, the way of instrumentation with a local network by three accelerometers is given via the above-mentioned three degrees of freedom.


Introduction
An ideal three degrees of freedom system that is equivalent with the modal behavior of an infinity number of degree of freedom of two cases of pendulums is proposed for each case. The first pendulum has hinges at the two ends, while the second pendulum is a cantilever (inverted pendulum). Both pendulums are presented analytically in the present article. This equivalent three DoF system can be used in instrumentation of such pendulums, where the concept of the concentrated masses is not existing, with a local network of three accelerometers. This issue is a main problem that appears very common during the instrumentation of inverted pendulums or bridge beams or steel stairs [1][2][3] or wind energy powers [4,5] in order to identify the real vibration mode-shapes and the fundamental eigen-frequencies of the structure via records of response accelerograms at specific positions due to ambient excitation [6,7].

First case: modal analysis of undamped amphi-hinge vertical pendulum with distributed mass and stiffness
According to the theory of continuous systems [8,9], consider a straight amphi-hinge vertical pendulum that is loaded by an external continuous dynamic loading p z x, t ð Þ, with reference to a Cartesian three-dimensional reference system oxyz, (Figure 1).
The vertical pendulum possesses a distributed mass m x ð Þ per unit height, which in the special case of uniform distribution is given as m x ð Þ ¼ m in tons per meter (tn/m). Furthermore, according to Bernoulli Technical Bending Theory, the beam has section flexural stiffness ΕΙ y x ð Þ, where in the special case of an uniform distribution of the stiffness it is given as ΕΙ y x ð Þ ¼ ΕΙ y , where E is the material modulus of elasticity and Ι y is the section moment of inertia about y-axis ( Figure 1). Next, we are examining a such amphi-hinge vertical pendulum that possesses constant value of distributed mass along its height, as well as constant value of distributed section flexural stiffness ΕΙ y . Due to fact that the vertical pendulum mass is continuously distributed, this pendulum/beam has infinity number of degrees of freedom for vibration along the horizontal oz-axis. In order to formulate of the motion equation of this beam, we consider an infinitesimal part of the vertical beam, at location x from the origin o, that has isolated by two very nearest parallel sections. The infinitesimal length of this part is the dx. On this infinitesimal length, we notice the flexural moment M x, t ð Þ, the shear force Q x, t ð Þ with their differential increments, while the axial force Ν x, t ð Þ is ignored, because it does not affect the horizontal beam vibration along z-axis. Moreover, noted the resulting force P z x, t ð Þ of the external dynamic loading. Therefore, we can write: where the resulting force P z x, t ð Þ acts at the total beam infinitesimal part. Furthermore, according to D'Alembert principle, the resulting inertia force F a x, t ð Þ is noted, where: Here, we agree that the time derivatives of the displacements are going to symbolize with full stops, while the spatial derivatives of the displacements are going to symbolize with accent. Next, the damping and the second order differential are ignored, so the force equilibrium on the infinitesimal part of the beam along z-axis gives: Moreover, the moment equilibrium with reference to weight center (w.c.) of the infinitesimal part of the beam (see Figure 1) gives: According to Euler-Bernoulli Bending theory (where the shear deformations are ignored), it is well-known that the following basic equation is true: Eqs. (4) and (5) are inserted into Eq. (3), so the motion equation without damping for the examined vertical beam is given: Eq. (6) is a partial differential equation that describes the motion u z x, t ð Þ of the vertical beam that is loaded with the external dynamic loading p z x, t ð Þ. In order to arise a unique solution from Eq. (6), the support conditions must be used at the two beam ends. It is worthy to note that the classical case of a beam with distributed mass and section flexural stiffness, under external horizontal excitation (Figure 2) on the two supports is mathematically equivalent with the vibration that is described by Eq. (6). Indeed, in the case of Figure 2, the total displacement u tot z x, t ð Þ of the beam at x-location is given: where u g t ð Þ is the displacement at the base, same for the two supports. But, it is known that the inertia forces of the beam are depended by the total displacement u tot z x, t ð Þ, while the distributed dynamic loading is null, p z x, t ð Þ ¼ 0. Thus, Eq. (3) is transformed into: Amphi-hinge vertical pendulum subjected with the same horizontal ground motion u g t ð Þ on the two end-supports.
Following, Eqs. (4) and (5) are inserting into Eq. (8), thus we are taken: By the comparison of Eqs. (6) and (9), we notice that the undamped beam vibration due to horizontal motion of the two supports is mathematically equivalent with the undamped vibration of the same beam where the two supports are fixed and the beam is loaded with the equivalent distributed dynamic loading p eq x, t ð Þ: In the case of the horizontal pendulum/beam free vibration without damping, we consider the first part of Eq. (9) that must be null: Furthermore, we ask the unknown spatial time-function u z x, t ð Þ, which is the solution of Eq. (11), must have the form of separated variants: where φ x ð Þ is an unknown spatial function and q t ð Þ is an unknown timefunction. Eq. (12) has been derived two times with reference to time-dimension t and more two times with reference to spatial-dimension x, so: Eqs. (13) are inserted into Eq. (11), giving: and, next, divided with the number m Á φ x ð Þ Á q t ð Þ, thus we are getting: The left part of Eq. (15) is a time-function, but the right part is a spatialfunction. In order to true Eq. (15) for all time values as well as for all spatial positions, the two parts of Eq. (15) must be equal with a constant λ. Thus, Eq. (15) is separated at two following differential equations: However, the time equation (16) indicates a free vibration of an ideal single degree of freedom system that has eigen-frequency ¼ ffiffi ffi λ p . Inserting the eigenfrequency ω into Eq. (17) gives: Next, we set the positive parameter β such as to be equal: because the parameters ω 2 , m, ΕΙ y are always positive. By the mathematic theory it is known that the general solution of Eq. (18) has the following form: where the four unknown parameters C 1 , C 2 , C 3 , C 4 must be calculated. In order to achieve this, four support conditions of the beam have to used. Indeed, for x ¼ 0 and x ¼ L the displacement u z 0, t ð Þ of the amphi-hinge vertical pendulum as well as the flexural moment M 0, t ð Þ, both are equal zero. The spatial function φ x ð Þ, which is the solution of Eq. (20), gives the modal elastic line of the beam. Having as known data that the following equation is true: The spatial function of the modal elastic line for x ¼ 0 is: and also for x ¼ 0, the function of the flexural moment due to examined modal elastic line of the beam is given by Eq. (5): Eq. (20) has been derived two times with reference to spatial-dimension x, thus arise: and Therefore, Eq. (23) is transformed: By Eqs. (22) and (26) directly arise C 2 ¼ 0 and C 4 ¼ 0, thus the general solution of Eq. (20) is the following: In addition, the parameters C 1 , C 3 are calculated by the support conditions of the second support of the beam. Therefore, for x ¼ L the vertical displacement u L, t ð Þ ¼ 0 be true. Thus, from Eq. (12) arises that φ L ð Þ ¼ 0 and Eq. (27) gives: In continuous, Eq. (23) gives: where φ 00 L ð Þ is directly getting from Eq. (25) that is equivalent with zero: However, when Eqs. (28) and (30) are re-written again, we get: And added part to part these two above-mentioned equations arise: But, the term sinh βL is not equal with zero, because then vibration is not existing. Therefore, C 3 has to equal with zero, so Eq. (28) is formed: Moreover, by Eq. (34) arise that either C 1 ¼ 0 that is mpossibility because φ x ð Þ 6 ¼ 0 by Eq. (27), either sin βL ¼ 0 that means the following equation must be true: However, Eq. (35) is transformed to Eq. (36): By the definition of parameter β, we can calculate the eigen-frequency ω: Thus, inserting Eq. (36) into Eq. (37), the eigen-frequency ω n directly arises for each n-value.
Therefore, the vibration mode-shape of the examined vertical pendulum arises by Eq. (27)-since previous inserting Eq. (35)-thus: The value of C 1 is arbitrary, and we usually get it equal to unit. Thus, for each value of parameter n, a mode-shape with its eigen-frequency is resulted. The fundamental (first) mode-shape is resulted for n ¼ 1, which shows a half sinusoidal wave, the second mode-shape shows a full sinusoidal wave, etc. (Figure 3). The order of the eigen-frequencies is ω 1 , ω 2 ¼ 4ω 1 , ω 3 ¼ 9ω 1 , ω 4 ¼ 16ω 1 , etc.

The equivalent three degrees of freedom system for amphi-hinges pendulum
At vertical pendulums or beams where the fundamental horizontal mode-shape does not activate the 90% of the total beam mass, we ask to consider the three first mode-shapes. Thus, for this purpose, we must define an ideal equivalent three degrees of freedom beam, which is going to give the first three frequencies and the first three mode-shapes of the examined beam. Therefore, which is the ideal three degrees of freedom system, where its three mode-shapes coincide with the real first three frequencies of the vertical pendulum with distributed mass and flexural stiffness?
In order to answer the above-mentioned question, consider a weightless vertical pendulum with height L and constant section in elevation, where carry three concentrated masses that each one has the same mass-value m eq , located per distance 0.25L, between one to one, and each one mass possesses an horizontal degree of freedom (Figure 4).
Furthermore, the pendulum flexibility matrix f can be calculated using a suitable method (Figure 4), and the inverse matrix gives the stiffness matrix k of the three degrees of freedom beam.
The equations of motion for the case of the free undamped vibration of the ideal beam is given as: The eigen-problem is written as: where the eigen-frequencies ω n and the three mode-shapes φ n are known by Eq. (38) and (39) and Figure 4. Therefore, the unique unknown parameter is the mass m eq . Thus, The numerical solution of Eq. (46) gives three roots for parameter m eq , where only the first root is acceptable, because the other two values rejected since do not have natural meaning (appear values greater from the total pendulum mass mL). Thus, the only one acceptable root is given as: Therefore, inserting the ideal equivalent mass m eq by Eq. (47) at three degrees of freedom system of Figure 4, the three eigen-frequencies coincide with the real values of the initial vertical pendulum that has distributed mass and flexural stiffness.

Discussion about the amphi-hinge vertical pendulums
By the previous mathematic analysis, an ideal three degrees of freedom system that is equivalent with the modal behavior of the amphi-hinge vertical pendulum with distributed mass and flexural stiffness along its height has been presented. This ideal three degrees of freedom system can be used in instrumentation of a such vertical tower, which does not possess concentrated masses. In the framework of the identification of mode-shapes of an amphi-hinge pendulum, the equivalent mass by Eq. (47) permits to locate accelerometers per 0.25L (as shown at Figure 4) and there measure the response horizontal acceleration histories, in order to calculate the real first three mode-shapes of the vertical pendulum in order to avoid the instability and resonance-vibrations between the examined vertical pendulum and the supported rockets or space bus before the last launched.

Second case: modal analysis of undamped inverted pendulum with distributed mass and stiffness
According to the theory of continuous systems [8,9], consider a straight vertical cantilever (inverted pendulum) that is loaded by an external continuous dynamic loading p z x, t ð Þ, with reference to a Cartesian three-dimensional reference system oxyz, (Figure 5).
The vertical inverted pendulum possesses a distributed mass m x ð Þ per unit height, which in the special case of uniform distribution is given as m x ð Þ ¼ m in tons per meter (tn/m). Furthermore, according to Bernoulli Technical Bending Theory, this cantilever has section flexural stiffness ΕΙ y x ð Þ, where in the special case of an uniform distribution of the stiffness it is given as ΕΙ y x ð Þ ¼ ΕΙ y , where E is the material modulus of elasticity and Ι y is the section moment of inertia about y-axis (Figure 1). Next, we are examining such an inverted pendulum that possesses constant value of distributed mass along its height, as well as constant value of distributed section flexural stiffness ΕΙ y . Due to fact that the cantilever mass is continuously distributed, this inverted pendulum/beam has infinity number of degrees of freedom for vibration along the horizontal oz-axis. In order to formulate of the motion equation of this beam, we consider an infinitesimal part of the vertical beam, at location x from the origin o, that has isolated by two very nearest parallel sections. The infinitesimal length of this part is the dx. On this infinitesimal length, we notice the flexural moment M x, t ð Þ, the shear force Q x, t ð Þ with their differential increments, while the axial force Ν x, t ð Þ is ignored, because it does not affect the horizontal cantilever vibration along z-axis. Moreover, noted the resulting force P z x, t ð Þ of the external dynamic loading. Eqs. (1)-(20) describe mathematically the modal analysis of the inverted pendulum. Afterward, the four unknown parameters C 1 , C 2 , C 3 , C 4 of Eq. (20) must be calculated taking account the conditions of the cantilever. In order to achieve this, four support conditions of the cantilever have to be used. Indeed, for the end at x ¼ 0 the displacement u z 0, t ð Þand the slope u 0 z 0, t ð Þ of displacement profile must be zero. So, Eq. (20) gives: and Moreover, for the free end at x ¼ L of the cantilever, the flexural moment M L, t ð Þ as well as the shear force Q L, t ð Þ must be both zero. So, from Eqs. (20) and (5) and afterward from Eq. (48) we take: However, re-writing Eqs. (50) and (51) again, we get the matrix form: Eq. (52) is a real eigenvalue problem. In order to calculate the eigenvalues, parameters C 1 and C 2 must not both equal zero. Thus, the determinant of the matrix by Eq. (52) must be zero, where it drives to the following frequency equation: However, Eq. (54) can be solved numerically only, where we obtain the first four roots (n ¼ 1, 2, 3, 4Þ: By the definition of parameter β, we can calculate the first four circular eigenfrequency ω using Eq. (36): Thus, inserting Eq. (54) into Eq. (55), the first four eigen-frequency ω n is directly arise for each n-value.
Therefore, the vibration mode-shape of the examined vertical inverted pendulum arises by Eq. (20)-since previous inserting Eq. (56)-thus: The value of C 1 is an arbitrary constant, and we usually get it equal to unit. Thus, for each value of parameter n, a mode-shape with its eigen-frequency is resulted. The fundamental (first) mode-shape is resulted for n ¼ 1, etc. (Figure 6).

The equivalent three degrees of freedom system of inverted pendulum
At inverted pendulums cantilevers, where the fundamental horizontal modeshape does not activate the 90% of the total cantilever mass, we ask to consider the three first mode-shapes. Thus, for this purpose, we must define an ideal equivalent three degrees of freedom beam, which is going to give the three mode-shapes of the examined beam. Therefore, which is the ideal three degrees of freedom system, where its three eigen-frequencies coincide with the real first three frequencies of the inverted pendulum with distributed mass and flexural stiffness?
In order to answer the above-mentioned question, consider a weightless vertical inverted pendulum with height L and constant section in elevation, where carry three concentrated masses that each one has the same mass-value m eq , located per distance 0.333L, between one to one, and each one mass possesses an horizontal degree of freedom (Figure 7).
The inverted pendulum displacement vector u of the three degrees of freedom, as well as the diagonal beam mass matrix m are written: Furthermore, the inverted pendulum flexibility matrix f can be calculated using a suitable method (Figure 7), and the inverse matrix gives the stiffness matrix k of the three degrees of freedom beam.  The equations of motion for the case of the free undamped vibration of the ideal beam is given as: The eigen-problem is written as: where the eigen-frequencies ω n and the three mode-shapes φ n are known by Eq. (56) and Figure 6. Therefore, the unique unknown parameter is the mass m eq . Thus, , B ¼ k 11 k 33 þ k 11 k 22 þ k 22 k 33 À k 2 12 À k 2 13 À k 2 23 ω 4 1 C ¼ À k 11 k 22 k 33 þ 2k 12 k 13 k 23 À k 11 k 2 23 À k 22 k 2 13 À k 33 k 2 12 ω 6 1 The numerical solution of Eq. (64) gives three roots for parameter m eq , where only the first root is acceptable, because the other two values rejected since do not have natural meaning (appear values greater from the total pendulum mass mL). Thus, the only one acceptable root is given: Therefore, inserting the ideal equivalent mass m eq by Eq. (65) at three degrees of freedom system of Figure 7, the three eigen-frequencies coincide with the real values of the initial vertical inverted pendulum that has distributed mass and flexural stiffness.
On the contrary to above-mentioned about the examined cantilever, it is worth noting that in the case where ask an ideal single degree of freedom system that has eigen-frequency equal to fundamental eigen-frequency of the real cantilever we write the following equations: • Eigen-frequency of single degree of freedom system with k the cantilever lateral stiffness and m eq, sdof the concentrated mass at the top of the cantilever: • Fundamental (first) eigen-frequency (see Eq. (56)) of the real cantilever with continuous distribution of mass and flexural stiffness: • Thus, inserting Eq. (67) into Eq. (66) the equivalent concentrated mass m eq,sdof at the top of the cantilever is given [9]: m eq,sdof ¼ 0:2426742 Á mL ð Þ (68)

Conclusions
The present article has presented a mathematic ideal three degrees of freedom system that is equivalent with the modal behavior of two cases of pendulums. First, an amphi-hinge vertical pendulum with distributed mass and flexural stiffness along its height has been examined. Second, an inverted pendulum that can be simulates a tower (or cantilever) has been examined too. For each case, an ideal three degrees of freedom system has been proposed that can be used in instrumentation of such a vertical tower, which does not possess concentrated masses. In the framework of the identification of mode-shapes of the above-mentioned pendulums, the equivalent mass by Eqs. (47) and (65) permits to locate accelerometers (as shown at Figures 4 and 7) and measures the response acceleration histories, in order to calculate the real first three frequencies.