Insights through Modal Minimal Models for Analysis of Linear and Nonlinear Dynamic Problems
Dynamic analysis of three-dimensional structures is common practice in industry to optimize products and in research to gain insights on the influence of parameters. The complexity differs based on the linearity or nonlinearity of the underlying problem, the type of excitation, e.g., forced or self-excitation, and the number of degrees of freedom that need to be examined. Reduced order models and optimized numerical methods are used to optimize the time and computational power needed to gain a certain insight. This article focusses on a specific class of problems where the modes of the structure do not or do not significantly change through the (damping) device or force that is added to the structure. Herein, lumped mass models are commonly used for analysis of the dynamic response of the system. In the article, it is highlighted that lumped mass models can give quantitative insight but modally reduced models allow a direct optimization of the problems with respect to the physical degree of freedom that, for example, is subject to self-excitation or dampened. The benefit of modal minimal models and its limitations are shown and discussed for different linear and nonlinear dynamic problems.
A dynamic analysis of complex three-dimensional structures is commonly done using numerical models with numerous degrees of freedom, whereas these models are efficient for a single linear analysis, e.g., to analyze the mechanical integrity of a structure, and they are inefficient for solving nonlinear problems or optimization problems, e.g., with contacting interfaces or the problem of the optimal placement of a nonlinear vibration mitigation device. Methods like the Harmonic Balance Method (HBM) , the Alternating Frequency-Time  (AFT) method, or conventional averaging methods  have been developed to determine the nonlinear response of the system, e.g., to periodical excitation. Herein, Reduced Order Models (ROMs) like the Component Mode Synthesis  are commonly used to represent the dynamic behavior of the structure for the specific problem and underlying parameters. These methods are especially useful if the nonlinear forces strongly influence the modes of the underlying structure or the nonlinearity cannot be approximated in a narrow frequency range. In this case, constraint modes [4, 5] can be used or a superposition of the modal ansatz functions/vectors that need to be able to picture the response of the nonlinear system.
The work done here is focusing on the special case that the modes of the structure are not significantly changed through the linear or nonlinear forces that are acting on the structure. In this case, lumped mass models are often used that represent a single degree of freedom (one mode) of the structure [6–10]. In contrast, this paper highlights the benefits of using modally reduced models [11–16] with one degree of freedom for this class of problems. For this purpose, different examples are shown that are usually optimized with lumped mass models or with complex three-dimensional numerical models. The examples include the efficient identification of the critical modes in a self-excited system and optimization of the placement of frictional and viscous damping devices and tuned-mass dampers. The examples are discussed and are benchmarked with the lumped mass models, and the suitability and limitations are discussed. The focus is on problems that can be solved analytically or semianalytically, and direct insight of parameter influences can be provided.
2. Comparing Modally Reduced and Lumped Mass Models
A mechanical system with several degrees of freedom is typically described using complex structural models based on finite elements or partial differential equations. The goal is to find minimal modal models that quantitatively deliver the same results and can be based on the complex models. These can be derived at the cost of a numerical or experimental modal analysis but can deliver analytical equations to predict loads or to optimize damper design and strategies in linear or nonlinear systems.
The approach using a modal single degree of freedom (SDOF) oscillator (modal DOF ) and a physical single-mass oscillator (physical DOF ) is analyzed, and the results are discussed. The equation of motion of a lumped mass model isor mass normalized as
The equation can be used to describe one mode with a modal damping value and an angular natural frequency . describes a force acting on the structure at a predefined point j. The ratio of the stiffness and the mass and the viscous damping need to be adjusted to account for a point of excitation j with respect to the mode shape of the structure and to adjust for the modal damping value and the natural frequency. The model is only valid for the point where it is observed which has to be the same point or node j where it is excited.
The modal single degree of freedom model with the modal DOF of any mode of the structure is determined directly from its mass normalized modal parameters:
Here, is the physical force that is acting at a position . The modal force at the right side of the equation is calculated by summation of the multiplication of the mass normalized mode shape amplitude with the physical forces for each position j, respectively. Complex excitation can be applied at multiple positions with different phase information.
Although both minimal models are leading to similar mathematical solutions, the modal approach offers certain advantages and more direct insights. Qualitative and quantitative insights are directly gained through the modal properties that are present in the analytical equations.
In case of the modal model, forces can be applied at every node and physical amplitudes can be calculated at every point of interest k by (assuming one contributing mode). The modal parameters are directly used and are not converted into substitute parameters, which is especially beneficial for investigating parameter influences such as frequency or mode shape changes as well as for parametric or statistical analysis. Furthermore, in the modally reduced model, the positions of the acting forces are directly considered in the differential equation via and multiple forces at multiple positions can be acting on the structure.
2.1. Critical Modes in Self-Excited Systems
The first example that highlights the benefits of a modal model is the identification of self-excited modes . There are different mechanisms for self-excitation of a system described in the literature. In the case where an elastic structure interacts with an interface, a common model for self-excitation is of a velocity-dependent contact force characteristic. Figure 1. Qualitatively a negative/positive slope of the nonlinear contact force characteristic corresponds to an energy input/energy dissipation for small monoharmonic vibrations with respect to a mean relative velocity or . The absolute value of the slope scales the net energy in one period of vibration linearly. The question answered here is which mode i is unstable or likely to be excited.
A modal model for each mode represented by the number is used that directly shows the influence of the point of excitation (the degree of freedom ) that is represented by the amplitude of the mode shape. The contact force characteristic is linearized with respect to a mean relative velocity .
This leads to the modal equation of motion:
Here, s denotes the slope of the contact force characteristic with respect to the mean relative velocity between the two interfaces . Assuming only one point of excitation j and one contributing mode i of the structure, all velocity-dependent parameters can be summarized resulting in an equivalent damping for the mode:
The sign of the equivalent damping is an indicator for the fix-point stability for the given parameters. The system is marginally stable if leading to two equivalent equations as follows:
That can be used to rank the stability of the modes with respect to the underlying self-excitation mechanism. is the (negative) slope of the velocity weakening contact force characteristic that needs to be exceeded (the absolute value) to get self-excitation of a mode with the modal properties , , and . The modal description accentuates the contribution of the amplitude of the mode shape at the contacting interface which is not visible in case of a lumped mass model.
The concept is visualized in Figure 1. An elastic structure is regarded that can be described by a specific number of modes. The degree of freedom j of the elastic structure is tangentially in contact with a rigid counterpart. The contacting interfaces at j have a relative velocity and a constant normal force orthogonal to j. The tangential force characteristic is assumed to be nonlinearly dependent on the relative velocity between the node j and linearly dependent on the normal force in the contact interface (Figure 1). As an example, the stability of a specific mode i is analyzed with reference to three points A, B, and C in Figure 1. In these points, the critical slope of the mode is related to the slope of the contact force characteristic. In point A, the relative velocity and the normal force are comparably low. The slope of the contact force characteristic is less negative compared with the mode specific value which is the slope of the black line through A. Therefore, following the logics of (5) and (6), not enough energy input is provided to overcome the energy dissipation through the modal damping (1% in the example). Therefore, point A is stable.
At point B, the normal force is doubled compared with point A. The contact force is twice as high. and the slope of the contact force characteristic is twice as steep. The critical slope of the mode , however, is constant and now less steep than the contact force characteristic. Again following the logics of equations (5) and (6), the system is unstable leading to self-excitation. A similar argumentation can be applied to point C that is stable.
2.2. Estimation of Amplitudes in Self-Excited Systems
The typical velocity-dependent contact force or torque characteristic is a superposition of frictional forces (sticking and sliding friction) and viscous forces. Generally, a positive slope is dissipative and a negative is exciting the system. The overall energy consumed or added in one period of the vibrations is given by the time-integral of the multiplication of the velocity-dependent contact force and the velocity and is zero in case of a limit cycle:
The slope is steep and positive for low relative velocity values close to zero (Figure 1). The reason is that the contact forces are dependent on the sign of the relative velocity. If the dissipative energy of this part is assumed to be dominant and the amplitude is approximated by a single harmonic , an upper amplitude  for the limit cyclecan be defined based on the assumption that the velocity fluctuation equals the mean velocity, . The loads and amplitudes are scaled by the ratio of the mode shape amplitudes to a point or node k within the structure. Again one contributing mode i is assumed. This assumption is valid because the dissipation of energy at low relative velocities is not only limiting the energy in the dominant mode i but also in all other modes that are typically completely suppressed. The approach leads to a simple and valid upper limit for the amplitude at any position of an elastic structure as a function of the modal properties of the mode.
2.3. Additional Damping Elements
The influence of different types of linear or nonlinear dampers such as fluid dampers, friction dampers, tuned-mass dampers can be investigated using modal models of the structure. The assumptions for the derivation of the analytical equations are a negligible influence of the forces and the damper device on the mode. Further, the interaction with other modes is negligible and no other modes are considerably excited due to other mechanisms or higher harmonic contents of nonlinear forces.
The model consists of two degrees of freedom (DOF) such as one modal DOF with the modal amplitude representing the structure and one physical DOF representing the physical displacement of the damper that is attached at the degree of freedom j. The physical displacement of the structure can be calculated by multiplication of the mode shape amplitude at the position or node j where the damper is forcedly connected to the structure. The equation of motion accounting for linear and nonlinear damping forces yields
The damper is characterized through a viscous damper constant d, a linear stiffness c, and an inertia J of the damper. The nonlinear forces of the damper can act on the position or node k of the structure and the inertia mass, respectively. The equation of motion is solved analytically or semianalytically for the different damper types. Herein, the equivalent modal damping ratio for a considered mode and the contributing critical parameters are derived and analyzed. In the following section, the analytical solutions of different damper principles are discussed. Figure 2 shows the qualitative curve of the damping provided by the various damper principles with respect to frequency (Figure 2, left) and amplitude (Figure 2, right). With respect to the frequency, the robustness of the damping effect is determined, whereas the change in amplitude can be used to make statements about the linearity. Thus, the fluid and tuned-mass dampers are linear while the inertia friction and friction bar dampers are nonlinear.
2.3.1. Inertia-Based Friction Damper
A rigid inertia is connected by frictional forces only to the structure in case of an inertia-based damper. Equation (9) is analyzed with the parameters and a frictional force . The derivation of the equivalent damping for the nonlinear system is described in detail the appendix. The general approach is representative for all damper principles discussed. An analytical solution  for the dissipated energy and the resulting equivalent modal damping ratio over a period are determined by assuming a harmonic motion of the th mode and using methods like averaging or harmonic linearization. Due to the nonlinearity, the damping coefficient that can be calculated for the friction damper is dependent on the modal amplitude of the considered mode (Figure 2, right). By determining the amplitude dependent optimum , the maximum provided damping and the corresponding amplitude can be determined as follows:
The solution shows that the damping is especially dependent on the amplitude of the mode shape and linearly dependent on the mass moment of inertia .
2.3.2. Stiffness-Based Friction Damper
The stiffness-based damper  bridges the amplitude at two points of the structure through a stiffness . To investigate a stiffness-based friction damper, equation (9) is used with and with and . In this case, the influence of the stiffness on the dynamic behavior is dominant in comparison to the inertia. The effect of the friction damper is thus dominated by the stiffness and the displacement between the points j and k. Again a monoharmonic motion is assumed. In addition, due to the comparably low inertia J, the inherent dynamics of the damper can be neglected for the derivation of the analytical solution . The nonlinearity of the friction force leads to an amplitude dependency of the damping coefficient (Figure 2, right). The maximum of the damping ratio and the corresponding amplitude are determined by
The maximum of the damping is increasing with the power of two with the difference of the amplitudes of the mode shape. The stiffness is contributing linearly to the damping. However, the stiffness may not be excessively increased to limit the effect on the mode shape and the amplitude difference. The stiffness-based damper is dependent on the angular natural frequency of the system and especially effective low frequencies whereas the inertia-based damper provides constant damping for the whole frequency range. The decisive difference is that the contact force of the stiffness-based damper is controlled by the displacement in comparison to the acceleration in case of the inertia-based damper.
2.3.3. Fluid Damper
The parameters of the ordinary fluid damper are , and . Assuming harmonic motion of the structure, an analytical solution for the dampingis obtained. With , an optimal damping ratio regarding an associated frequencyis derived. Similar to the friction damping, the equation incorporates one of the most important factors that is the placement with respect to the mode shape represented by the value of that is contributing to the damping with the power of two. The damping effect prevails over the entire amplitude range of the structural vibration. For a lumped mass model, similar results can be derived that do not contain any information regarding the placement of the damper .
2.3.4. Tuned Mass Damper
By adding a stiffness to the fluid damper, the natural frequency of the damper system can be tuned to the vibrational frequency of the underlying structure. The damping of the resulting tuned-mass damper (TMD) is strongly dependent on the frequency (Figure 2, left) . Based on the previously discussed assumptions, the relative movement between the inertia of the damper and the underlying elastic structure is mathematically similar to a system with a mass imbalance. Hence, an analytical approximation of the equivalent modal dampingof the modal damping ratio of the th mode provided by the TMD can be determined. The frequency ratio ( is the vibration frequency of the structure; is the natural frequency of the TMD) is particularly important. The damping provided can be approximated for different exceptional cases of the characteristic TMD frequency to
Equation (15) shows that for small natural frequencies of the TMD , i.e., for systems where , the resulting damping corresponds to the damping of the fluid damper. If the stiffness is comparably high , the damping effect of the TMD tends to 0. The frequency with the highest damping is at if the inertia of the damper is rather small . Besides the vibration frequency and the position of the TMD, the inertia and the damping constant are decisive. It should be noted that unrealistically high inertias or too low damping constants significantly influence the system behavior and thus render the approximated solution unsuitable.
In structural dynamics, several analytical and numerical methods exist to calculate the linear or nonlinear response of the system and to estimate system properties. Analytical equations have limited applicability but they are commonly used to gain direct insight. Single or two degrees of freedom lumped mass models are commonly used for this purpose. In the article, modal models are proposed that are based on exactly the same math.
The article discusses the benefits of a modal description for insights into various problems where the resulting linear or nonlinear equation of motion can be solved analytically or semianalytically. It is shown that the major benefit of the modal description is the visibility of the influence of the modal parameters such as the amplitude of the mode shape at a point of interest. Therefore, the point of observation of an amplitude or the placement of a damper can be derived at every position of a 3D-structure without any changes to the equations that are used. Exemplary results that oppose the results of a lumped mass models and modal models are summarized in Table 1.
The analytical solution for the equivalent damping coefficient of an inertia friction damper  with the equation of motion according to equation (9) with and a sliding friction force is derived using the harmonic balance method . Herein, a harmonic amplitude of the th mode is assumed. Based on this assumption, the relative displacement/velocity in the friction contact and subsequently the dissipated energy are estimated. For this purpose, the times , and that correspond to the transitions between sticking and sliding friction in one period of vibration are calculated. For modal amplitudes , the inertia force of the damper is smaller than the sticking friction force limit and only sticking occurs in one oscillation period with . Hence, no energy is dissipated . For higher amplitudes, the dissipated energy in one oscillation periodcan be calculated by four integrals that correspond to sticking and sliding regimes. The energy is only dissipated in the sliding regimes and . The transition times can be calculated analytically and semianalytically as described in . Generally, a modal dissipated energy of a linear system in one period and the corresponding linear modal damping ratio for one mode are determined by
By assuming the linear dissipated energy and the nonlinear dissipated energy are equal, the equivalent damping ratiocan be estimated. By solving equation (A1) (cf. ) and inserting the resulting equation into equation (A3), an analytical solution withis derived with an amplitude dependent maximum at at .
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
The authors would like to thank Baker Hughes for supporting the work and giving permission to publish this paper. The authors acknowledge support by the Open Access Publication Funds of the TU Braunschweig.
K. Magnus, K. Popp, and W. Sextro, Schwingungen: Physikalische Grundlagen und mathematische Behandlung von Schwingungen, Springer Fachmedien Wiesbaden, Wiesbaden, Germany, 2013.
T. M. Cameron and J. H. Griffin, “An alternating frequency/time domain method for calculating the steady-state response of nonlinear dynamic systems,” Journal of Applied Mechanics, vol. 56, no. 1, pp. 149–154, 1989.View at: Publisher Site | Google Scholar
J. A. Sanders, J. Murdock, and F. Verhulst, Averaging Methods in Nonlinear Dynamical Systems, Springer, New York, NY, USA, 2007.
R. R. Craig and M. C. C. Bampton, “Coupling of substructures for dynamic analyses,” AIAA Journal, vol. 6, no. 7, pp. 1313–1319, 1968.View at: Publisher Site | Google Scholar
M. P. Castanier, Y.-C. Tan, and C. Pierre, “Characteristic constraint modes for component mode Synthesis,” AIAA Journal, vol. 39, no. 6, pp. 1182–1187, 2001.View at: Publisher Site | Google Scholar
K. E. Hafner and H. Maass, Torsionsschwingungen in der Verbrennungskraftmaschine, Springer, Vienia, Australia, 1985.
I. López, J. M. Busturia, and H. Nijmeijer, “Energy dissipation of a friction damper,” Journal of Sound and Vibration, vol. 278, no. 3, pp. 539–561, 2004.View at: Publisher Site | Google Scholar
K. Popp, L. Panning, and W. Sextro, “Vibration damping by friction forces: theory and applications,” Journal of Vibration and Control, vol. 9, no. 3-4, pp. 419–448, 2003.View at: Publisher Site | Google Scholar
J. P. Den Hartog and G. Mesmer, Mechanische Schwingungen, Springer, Berlin, Germany, 1952.
S. V. Bakre and R. S. Jangid, “Optimum parameters of tuned mass damper for damped main system,” Structural Control and Health Monitoring, vol. 14, no. 3, pp. 448–470, 2007.View at: Publisher Site | Google Scholar
D. J. Ewins, Modal Testing: Theory, Practice and Application, Research Studies Press, Baldock, England, 2000.
A. Hohl, M. Tergeist, H. Oueslati et al., “Prediction and mitigation of torsional vibrations in drilling systems,” in Proceedings of the IADC/SPE Drilling Conference and Exhibition, Texas, TX, USA, September 2016.View at: Google Scholar
A. Hohl, M. Tergeist, H. Oueslati et al., “Derivation and experimental validation of an analytical criterion for the identification of self-excited modes in drilling systems,” Journal of Sound and Vibration, vol. 342, no. 1, pp. 290–302, 2015.View at: Publisher Site | Google Scholar
V. Kulke, G.-P. Ostermeyer, M. Tergeist, and A. Hohl, “Semi-analytical approach for derivation of an equivalent modal friction-damping ratio and its application in a self-excited drilling system,” Dynamics, vol. 4, 2019.View at: Publisher Site | Google Scholar
V. Kulke, P. Thunich, F. Schiefer, and G.-P. Ostermeyer, “A method for the design and optimization of nonlinear tuned damping concepts to mitigate self-excited drill string vibrations using multiple scales lindstedt-poincaré,” Applied Sciences, vol. 11, no. 4, p. 1559, 2021.View at: Publisher Site | Google Scholar
M. Neubauer and J. Wallaschek, “Vibration damping with shunted piezoceramics: fundamentals and technical applications,” Mechanical Systems and Signal Processing, vol. 36, no. 1, pp. 36–52, 2013.View at: Publisher Site | Google Scholar
L. P. Shilnikov, A. L. Shilnikov, D. V. Turaev, and L. O. Chua, Methods of Qualitative Theory in Nonlinear Dynamics: Part I, World Scientific Pub. Co, Singapore, 1998.