Abstract

The development of reduced-order models remains an active research area, despite advances in computational resources. The present work develops a novel order-reduction approach that is designed to incorporate isolated regions that contain, for example, nonlinearitites or accumulating damage. The approach is designed to use global modes of the overall system response, which are then naturally coupled to the response in the isolated region of interest. Two examples are provided to demonstrate both the accuracy and the computational efficiency of the proposed approach. The performance of this approach is compared to the exact response corresponding to a finite element simulation for the chosen problems. In addition, the accuracy and computational efficiency are shown relative to a standard Galerkin reduction based on the linear normal modes. It is found that the proposed reduction offer computational efficiency comparable to a Galerkin reduction, but more accurately represents the response of the system when both are compared to the finite element simulation.

1 Introduction

Despite recent advances in computing power, large-scale finite element analysis (FEA) can still present computational difficulties simply arising from, for example, the overall size of the problems considered or large collections of analyses required for uncertainty quantification. Therefore, reduction techniques originating when computational resources were a fraction of current capabilities continue to retain their importance and usefulness. Component mode synthesis (CMS), based on the concept of substructuring, is one of the most widely used reduction methods in structural dynamics. Originally based on the the work of Hurty [1] and extended by Craig and Bampton [2], CMS decomposes the response of an individual component of a larger structure into component and constraint modes [3,4]. The dominant modes of the component can be identified and retained, thus significantly reducing the overall order of the system, while the constraint modes can be used to couple adjacent components together. Such an approach is particularly attractive when the response of the overall structure is localized to only a few of the components of the system, but can be used to develop reduced-order models of the response of complex structural systems.

However, there are numerous examples where the response of the overall system can be characterized as global in nature, so that the response of the system is more closely represented by the structural modes of the overall system rather than localized within only a few components. Jointed structures provide a canonical example of such systems [510], where the response of the overall system with joints is closely related to the modes of the overall structure in the absence of the interfaces, known as the monolithic system. In such a case, the structural modes of the monolithic system are significantly different from the vibration modes of the individual components. As identified by Segalman [1113], in such cases, it is desirable to describe the response and ultimately develop reduced-order models in terms of the structural modes of the overall system, rather than a collection of component-level modes. However, in such a system, the dynamics introduced by the presence of the interfaces, which are internal to the structure, nonetheless significantly influence the response of the overall system [6].

More generally, the response of system with isolated features, including nonlinearitites, could in principle be represented by CMS, so that the isolated effects are introduced through coupling between the constraint modes, while the exterior linear regions were decomposed into their individual component modes. However, the same issue remains—the response of the overall system may be better represented by the modes of the overall system rather than those of the individual (linear) components. In the present work, we develop a novel order-reduction method for systems with isolated nonlinearitites. In contrast to CMS, which utilizes the modes of the individual components to achieve the reduction in the overall model order, the present method employs the modes of the global system, obtained from the corresponding monolithic system in the absence of the nonlinearitites.

2 Global Model Reduction

We consider a general structure as illustrated in Fig. 1, composed of two exterior regions, identified as C1 and C2, which are described by their internal variables x_1RN1 and x_2RN2, respectively. In addition, there exists an internal isolated region Cr with variables x_rRNr. The equations of motion for the exterior regions C1 and C2 are assumed to be linear, while the localized interior region includes terms represented by N_RNr, which could describe nonlinearitites in the isolated region arising from, for example, the presence of joints and interfaces, or accumulating damage in the interior region. The overall linear mass, damping, and stiffness matrices are represented as M, C, and K, respectively, so that the overall equations of motion are
(1)
with
(2)
The mass, damping, and stiffness matrices are defined over the entire structure and in particular include the linear terms in the equations of motion within the region Cr. Moreover, in this formulation, it is assumed that the additional terms defined by N_ do not depend on the response of the exterior regions as described by x_1 and x_2 but instead only depend on the response in the isolated region, so that N_=N_(x_r). In the investigations below, this system is referred to as the reference system. A mixed formulation for the equations of motion, similar to that presented in Quinn [14], is developed here but extended to a form appropriate for the equations of motion arising from a finite element discretization of the equations of motion for a general structural system.
Fig. 1
Structural system with isolated nonlinearities
Fig. 1
Structural system with isolated nonlinearities
Close modal

2.1 Mixed Formulation.

The variables x_r within the interior region are further decomposed as
(3)
where the components described by x_α and x_β describe degrees-of-freedom (DOF) on the boundary between components C1 and Cr, and C2 and Cr, respectively. In particular, the coupling between x_1 and x_α is linear, as is that between x_2 and x_β. However, the degrees-of-freedom x_c within Cr are assumed to be nonlinearly coupled to these boundary coordinates x_α and x_β. Therefore, the general N degree-of-freedom mechanical system given in Eq. (1) can be expanded in the form
(4)
Recall that the nonlinearities described by N_() are localized to the interior region Cr. As a result they do not depend on the response of the exterior regions, as described by x_1 and x_2, but instead depend only on the response in the interior region, as indicated in Eq. (4).
Note that this representation is very similar to that taken in the development of the CMS approach. For example, the equations of motion restricted to the adjoining regions C1 and Cα can be expressed as
(5)
so that C1 describes a subcomponent and Cα represents the boundary interval, while the right-hand side describes the external forcing on the appropriate region.
Over the interior region let y_ be the response of an ideal system to be subsequently defined, so that the response of the reference system in the isolated region can be expressed as
(6)
Thus, z_ is the difference between the response of the ideal and actual system, referred to as the deviatoric response. Provided that N_ is small these are expected to be in some sense close. Using this expansion the equations of motion for C1 for example can be expressed as
(7)
This equation of motion exactly represents the response of the original structure but requires the response for the deviatoric displacement z_α. A similar expression holds in the exterior region C2. Nonetheless proceeding, a mixed displacement vector can be defined as
(8)
where (x_1,x_2) represent the exact response in the exterior regions and (y_α,y_c,y_β) described the response of the ideal system in the interior region. The resulting equations of motion for w_ can be written as
(9)
or
(10)
with Q_ defined as
(11)
Similar to the definition of z_, Q_ can be interpreted as the deviatoric force acting on the external regions from the internal region. The displacement vector w_ mixes the response of the original system in the exterior regions C1 and C2 with that of the ideal system in the interior region. Note that the ideal system is not the system in the absence of the nonlinearities as described by N_=0_. Instead, the coordinates y_ describe the response of the system with the additional forcing arising from Q_.

The closure of this system requires the deviatoric response (z_α,z_β) in the boundary regions be determined outside of this system of equations. In addition, Eq. (10) cannot be used to determine the response within the interior region defined by the boundary and isolated regions (Cα,Cc,Cβ). However, if (z_α,z_β) can be determined then the resulting response in the exterior regions is exact. In Ref. [14], the deviatoric force Q_ was specified as a constitutive equation, based on a reduced-order model for an interface as developed by Ref. [15].

In terms of y_ and z_, the equations of motion in the interior region can be reduced to
(12)
or
(13)
Thus, the deviatoric response as identified by z_ is driven through the nonlinearities N_ in the interior region by the response of the ideal displacements y_. The equations of motion in terms of w_ and z_ are referred to below as the unreduced isolated system.

In summary, the nonlinear equations of motion for the reference system given in Eq. (1) are written in terms of x_RN, with N = N1 + Nr + N2. These have been decomposed into a set of linear equations in terms of the mixed variable w_RN given in Eq. (10) together with a system of nonlinear equations given in Eq. (13), written in terms of z_RNr. At first glance, this appears to be a step backwards, as the original N degree-of-freedom system has been replaced by a system with N + Nr degrees-of-freedom. However, the linear equations for w_ may be able to be reduced using existing model reduction techniques, so that if NNr, a significant computational reduction may be achieved and, as illustrated below, a significant increase in accuracy as compared to traditional Galerkin reduction in terms of the eigenmodes of the linear system.

2.2 Modal Analysis.

With N_0_, Eq. (1) has no simple representation and in particular no modal decomposition is assumed to be available2. However, such a eigenanalysis can be performed in the linear system with N_0_. More generally, various order reduction techniques could be applied to this linear system, including projection based methods [16]. Moreover, assuming proportional damping, the undamped modes ϕ_ defined as
(14)
can be used to decouple the mixed coordinates w_. In particular, these identified modes are global to the system and are not localized to a single exterior region, in contrast to a typical reduction based on component mode synthesis. Using the resulting eigenvectors, the modal decomposition of Eq. (10) reduces to
(15)
with qk(t) the kth modal coordinate, defined as
(16)
where δmn is the Kronecker delta function. Other expansions, including those based on Proper Orthogonal Decompositions [1719], could be employed as well. In this, ϕ_kTf_(t) describes the modal forcing while ϕ_kTQ_ describes the effect of the isolated nonlinearitites, both acting on the kth mode. However, as described above, the deviatoric forcing Q_ must be determined outside of Eq. (10). In particular, the deviatoric displacement z_ is driven by the mixed response of the interior region as defined by
(17)
which is expressed in terms of the modal coordinates through the mixed coordinates
(18)
If the response of the interior region in the mixed formulation can be well approximated by an M mode truncation then the deviatoric response in the interior region as defined by z_ can be determined from Eq. (13), which can then be used to determine the reduced-order modal response of the system as given in Eq. (15). With the modal reduction, the resulting system has M + Nα + Nc + Nβ degrees-of-freedom. While MN, typically the retained modes offer a significant reduction in the overall degree-of-freedom for the system, so that MN1 + N2. This formulation nonetheless couples together all of the retained modes through the reconstruction of the mixed response in the interior region from the modal expansion, as determined from Eq. (18). The resulting response of the system can then be reconstructed from the modal coordinates q_(t) and the interior coordinates z_(t) as
(19)

3 Example: Finite Element Analysis Rod Elements

As an illustrative example for this approach, consider the finite discretization of a one-dimensional rod with fixed-free boundary conditions. If overall mass and length of the rod are m and ℓ, respectively, while the axial rigidity is EA, then for the the uniform rod with N elements the corresponding mass and stiffness matrices are
(20)
where the elements are equally spaced, so that the distance between nodes is Δ = ℓ/N. The spatial position of the rod, denoted s lies in the interval 0 ≤ s ≤ ℓ, and if the interior region is located in the interval s ∈ (s1, s2), then the boundary nodes are assumed to correspond to α=s1N and β=s2N, so that
(21)
The nonlinearities between nodes i and j in the isolated region are assumed to depend on the relative displacement between node pairs, defined as vi,j(t) ≡ xi(t) − xj(t), so that individual forces are assumed to take the form
(22)
where h[v(t)] is a hysteretic operator acting on the relative displacement v(t). These additional terms used in the interior region are representative of nonlinearities that are often seen in structural systems, including linear detuning and cubic stiffness and damping terms, as well as a hysteretic force commonly arising from an interface [20]. These terms were chosen to highlight the ability of the proposed approach to incorporate such common terms. Therefore, the nonlinear force acting between nodes i and j is defined as
(23)
The resulting nonlinearities in the interior region are therefore identified as
(24)
In this description, h[v(t)] is a force component arising from hysteretic damping and is determined by a differential constitutive equation. In the case of a Bouc model for the hysteretic damping, the equations of motion are augmented by constitutive relations of the form
(25)
Note that any representation for the hysteretic damping could be used in the above formulation, including Bouc-Wen or Iwan models [20]. In the present work, Bouc models are chosen for simplicity [21].
For the corresponding undamped linear system, the energy in each mode is defined as
(26)
so that the total linearized energy in the system is defined as
(27)
The error between the reference system and the reduced order approximation can be defined as
(28)
where E_ is the vector of modal energies, so that [E_]i=Ei.
For a system with N degrees-of-freedom, the system is simulated with the nominal parameters
(29)
The underlying linear model for the rod could in principle be nondimensionalized so that the above choices for (ρ, ℓ, EA) do not restrict the system dynamics. In the internal region the parameters associated with the nonlinearities are scaled by the number of elements N, so that the equivalent stiffness and damping across the internal interval s ∈ (s1, s2) = (0.25, 0.35) remains (approximately) constant as the number of elements varies.
The original equations of motion are simulated with N = 80, and for this discretization, the number of elements in each region is
(30)
To illustrate the performance of the above model reduction approach, the model described in Sec. 2.2 is simulated and the response is compared against both the response of the original equations of motion and the response from a standard Galerkin reduction. The Galerkin reduction is obtained by assuming that the response is expanded as
(31)
analogous to the expansion given in Eq. (18), but substituted directly into the original coordinates without isolating the response of the interior region.

For this comparison, the initial conditions of the system are chosen so that the rod is statically displaced by a unit force applied at s = 0.30 and then released from rest. Note that this initial state excites every mode in the structure, and the response of the system at s = 0.50 is shown in Fig. 2. In Fig. 2(a), the response of the system is shown for the original equations of motion, while in Fig. 2(b), the response is shown at the same location for the modally reduced isolated model, retaining M = 10 structural modes. Note that this system has nine degrees-of-freedom associated with the interior region (c.f., Eq. (12)), so that the total degree-of-freedom of the reduced model is 19. Finally, for comparison, the response for a 19 degrees-of-freedom Galerkin reduction is shown in Fig. 2(c).

Fig. 2
Numerical simulation; x(t)|s=0.50; in panels (b) and (c) the response of the original system is also shown with the dashed line: (a) original formulation, (b) reduced isolated formulation (M = 10, DOF = 19), and (c) Galerkin reduction (M = 19)
Fig. 2
Numerical simulation; x(t)|s=0.50; in panels (b) and (c) the response of the original system is also shown with the dashed line: (a) original formulation, (b) reduced isolated formulation (M = 10, DOF = 19), and (c) Galerkin reduction (M = 19)
Close modal

In Fig. 3(a) the error as defined by Eq. (28) is shown for both the localized reduced-order model as well as the corresponding Galerkin reduction. Although not shown, when the system is simulated with the unreduced isolated formulation given in Eqs. (10) and (13), the response is identical to the simulation with the original equations up to numerical accuracy of the integrator used. The nonlinearities are localized in the interior region s ∈ (s1, s2) = (0.25, 0.35). As a result, accurate resolution of the dynamics in this region is critical to capture the effect of the nonlinearities, including the energy loss due to dissipation. In Fig. 3(b), the error in the response of the interior region is shown as a function of time for the modally reduced isolated model as well as for the Galerkin reduction, both retaining 19 total degrees-of-freedom. As illustrated, the isolated model is more accurate in the interior region as compared to the Galerkin reduction.

Fig. 3
Error (N = 80, DOF = 19); solid curves represent the reduced isolated system (M = 10), short-dashed curve represents Galerkin reduction: (a) modal energy error and (b) displacement error in the interior region; Δ ≡ ||Xexact − Xreduced||s∈(0.25,0.35)
Fig. 3
Error (N = 80, DOF = 19); solid curves represent the reduced isolated system (M = 10), short-dashed curve represents Galerkin reduction: (a) modal energy error and (b) displacement error in the interior region; Δ ≡ ||Xexact − Xreduced||s∈(0.25,0.35)
Close modal
The performance of this isolated reduction technique is further investigated for both the modally reduced isolated model and the Galerkin reduction as the total degree-of-freedom varies. In Fig. 4(a), the energy error given in Eq. (28) is shown as M, the number of structural modes retained, varies, while in Figs. 4(b), the computational efficiency is shown, defined as the ratio of the computational time of the reduced system relative to that of the original equations of motion, that is
(32)
In each simulation, the computational time T is measured using the matlab’s built-in toc function. As observed in Fig. 4(a), there is a minimum degree-of-freedom required for the modally reduced isolated model due to the fixed number of degrees-of-freedom associated with the interior region, which with N = 80 is 9. The error in the reduced isolated system quickly decreases as the number of retained spatial modes increases, while the error in the Galerkin reduced model decreases more slowly. However, because of the additional degrees-of-freedom associated with the interior region in the reduced isolated system, it requires a minimum number of spatial modes to match the accuracy of the Galerkin approximation, in this case with M = 6. Above this spatial resolution, however, the reduced isolated model with retaining M > 6 spatial modes is more accurate than the corresponding Galerkin approximation with M + 9 spatial modes. Moreover, as seen in Fig. 4(b), the computational effort associated with the reduced isolated model is comparable to that of the Galerkin reduction of the same size. Note that the simulation uses the matlab adaptive integrator ode23t, so that the computational time depends on a number of factors in additional to model size, including internal integrator tolerances and estimated numerical error.
Fig. 4
Error (N = 80); solid curves represent reduced isolated system (M = 10), short-dashed curves represent Galerkin reduction: (a) modal energy error and (b) displacement error in the interior region
Fig. 4
Error (N = 80); solid curves represent reduced isolated system (M = 10), short-dashed curves represent Galerkin reduction: (a) modal energy error and (b) displacement error in the interior region
Close modal

4 Example: Finite Element Analysis Beam Elements

Extending on the previous example, the bending deformation of a cantilever beam of unit length is considered. In addition, the equations of motion can be nondimensionalized so that both EI and ρA are considered to be unity and the system is modeled using Euler-Bernoulli beam elements. Each node contains both rotational and transverse degrees-of-freedom, so an P node discretization has N = 2P degrees-of-freedom. In addition, proportional damping is used with ζ = 0.02. The system is assumed to have contain cubic stiffness and damping nonlinearities over the interval s ∈ (0.40, 0.50). External transverse loading is applied to the beam locations at s = s1 = 0.33 and s = s2 = 0.67, of the form
(33a)
(33b)
This particular excitation was chosen somewhat arbitrarily, but contains two frequencies that are near the lowest modal frequency of ω = 3.516 of the beam. Finally, the beam is initially at rest in its undeformed position.

With P = 35 nodes for the system, the resulting system has N = 70 degrees-of-freedom. As illustrated in Fig. 5(a), the average error as defined in Eq. (28) depends on the number of modes retained in the expansion of the structural system, and as M increases the error decreases, while as shown in Fig. 5(b), the computational time ratio is reduced for M < N, and approaches unity as MN.

Fig. 5
Modally reduced isolated system (N = 70): (a) mean energy error, E, and (b) computational time, η
Fig. 5
Modally reduced isolated system (N = 70): (a) mean energy error, E, and (b) computational time, η
Close modal

5 Conclusions

A novel order-reduction approach has been developed specifically for systems that contain isolated features, including nonlinearitites. The approach assumes linear exterior regions connected by localized regions that could, for example, include nonlinearities such as those arising from joints and interfaces. In contrast to a familiar component mode synthesis, the overall system is reduced using the global modes of the system, which then couple to the deviatoric effects introduced by the isolated region. The result is a reduced-order modeling approach that accurately represents the response in the isolated region and therefore can accurately describe the nonlinearities, yet is nonetheless written in terms of the global modes of the system. Two examples are considered, the transient response of a longitudinal rod and the forced response of an Euler-Bernoulli beam. The former is compared against both a standard FEA analysis as well as a Galerkin expansion using the linear mode shapes; the latter is compared against FEA analysis. Note however, that the response within the isolated region as described by z_ has not been reduced. Future efforts can seek to develop appropriate reduction methods focused on this interior region, that can then be combined with the framework presented here.

Footnote

2

Note that N_ need not necessarily be nonlinear, but instead could, for example, include some linear loss of stiffness due to increasing damage. While in this instance a modal decomposition could in principle be undertaken in the damaged state, it is not utilized here. In the present work any additional terms in N_ will be referred to as “nonlinear” regardless of the specific form of the terms.

Acknowledgment

Sandia National Laboratories is a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International, Inc., for the U.S. Department of Energy’s National Nuclear Security Administration under contract DE–NA0003525. This paper describes objective technical results and analysis. Any subjective views or opinions that might be expressed in the paper do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

Conflict of Interest

There are no conflicts of interest.

Data Availability Statement

The datasets generated and supporting the findings of this article are obtainable from the corresponding author upon reasonable request. The authors attest that all data for this study are included in the paper. Data are provided by a third party listed in Acknowledgment.

References

1.
Hurty
,
W. C.
,
1965
, “
Dynamic Analysis of Structural Systems Using Component Modes
,”
AIAA J.
,
3
(
4
), pp.
678
685
. 10.2514/3.2947
2.
Craig.
,
R. R.
, and
Bampton, Jr
,
M. C. C.
,
1968
, “
Coupling of Substructures for Dynamic Analyses
,”
AIAA J.
,
6
(
7
), pp.
1311
1319
. 10.2514/3.4741
3.
Craig, Jr
,
R. R.
, and
Kurdila
,
A. J.
,
2006
,
Fundamentals of Structural Dynamics
. 2nd ed.,
John Wiley & Sons, Inc.
,
Hoboken, NJ
.
4.
Farhat
,
C.
, and
Geradin
,
M.
,
1994
, “
On a Component Mode Synthesis Method and Its Application to Incompatible Substructures
,”
Comput. Struct.
,
51
(
5
), pp.
459
473
. 10.1016/0045-7949(94)90053-1
5.
Brake
,
M. R. W.
, editor,
2018
,
The Mechanics of Jointed Structures: Recent Research and Open Challenges for Developing Predictive Models for Structural Dynamics
,
Springer
,
Cham, Switzerland
.
6.
Gaul
,
L.
, and
Lenz
,
J.
,
1997
, “
Nonlinear Dynamics of Structures Assembled by Bolted Joints
,”
Acta Mech.
,
125
(
1
), pp.
169
181
. 10.1007/BF01177306
7.
Khattak
,
A. R.
,
Garvey
,
S.
, and
Popov
,
A.
,
2010
, “
Proper Orthogonal Decomposition of the Dynamics in Bolted Joints
,”
J. Sound. Vib.
,
329
, pp.
1480
1498
. 10.1016/j.jsv.2009.11.016
8.
Festjens
,
H.
,
Chevallier
,
G.
, and
Dion
,
J. L.
,
2014
, “
Nonlinear Model Order Reduction of Jointed Structures for Dynamic Analysis
,”
J. Sound. Vib.
,
333
, pp.
2100
2113
. 10.1016/j.jsv.2013.11.039
9.
Roettgen
,
D. R.
, and
Allen
,
M. S.
,
2017
, “
Nonlinear Characterization of a Bolted, Industrial Structure Using a Modal Framework
,”
Mech. Syst. Signal Proc.
,
84
, pp.
152
170
. 10.1016/j.ymssp.2015.11.010
10.
Lacayo
,
R. M.
, and
Allen
,
M. S.
,
2019
, “
Updating Structural Models Containing Nonlinear Iwan Joints Using Quasi-static Modal Analysis
,”
Mech. Syst. Signal Proc.
,
118
, pp.
133
157
. 10.1016/j.ymssp.2018.08.034
11.
Segalman
,
D. J.
,
A modal approach to modeling spatially distributed vibration energy dissipation
.
Technical Report SAND2010–4763
,
Sandia National Laboratories
,
August
2010
.
12.
Segalman
,
D. J.
,
2006
, “
Modelling Joint Friction in Structural Dynamics
,”
Struct. Control Health Monitoring
,
13
, pp.
430
453
. 10.1002/stc.119
13.
Segalman
,
D. J.
,
2007
, “
Model Reduction of Systems With Localized Nonlinearities
,”
ASME J. Comput. Nonlinear. Dyn.
,
2
(
3
), pp.
249
266
. 10.1115/1.2727495
14.
Quinn
,
D. D.
,
2012
, “
Modal Analysis of Jointed Structures
,”
J. Sound. Vib.
,
331
(
1
), pp.
81
93
. 10.1016/j.jsv.2011.08.017
15.
Miller
,
J. D.
, and
Quinn
,
D. D.
,
2009
, “
A Two-sided Interface Model for Dissipation in Structural Systems With Frictional Joints
,”
J. Sound. Vib.
,
321
(
1–2
), pp.
201
219
. 10.1016/j.jsv.2008.09.037
16.
Benner
,
P.
,
Gugercin
,
S.
, and
Willcox
,
K.
,
2015
, “
A Survey of Projection-based Model Reduction Methods for Parametric Dynamical Systems
,”
SIAM Rev.
,
57
(
4
), pp.
483
531
. 10.1137/130932715
17.
Feeny
,
B. F.
, and
Kappagantu
,
R.
,
1998
, “
On the Physical Interpretation of Proper Orthogonal Modes in Vibrations
,”
J. Sound. Vib.
,
211
(
4
), pp.
607
616
. 10.1006/jsvi.1997.1386
18.
Chatterjee
,
A.
,
2000
, “
An Introduction to the Proper Orthogonal Decomposition
,”
Curr. Sci.
,
78
(
7
), pp.
808
817
.
19.
Kerschen
,
G.
,
Golinval
,
J.-C.
,
Vakakis
,
A. F.
, and
Bergman
,
L. A.
,
2005
, “
The Method of Proper Orthogonal Decomposition for Dynamical Characterization and Order Reduction of Mechanical Systems: An Overview
,”
Nonlinear Dyn. Vol.
,
41
, pp.
147
169
. 10.1007/s11071-005-2803-2
20.
Mathis
,
A. T.
,
Balaji
,
N. N.
,
Kuether
,
R. J.
,
Brink
,
A. R.
,
Brake
,
M. R. W.
, and
Quinn
,
D. D.
,
2020
, “
A Review of Damping Models for Structures with Mechanical Joints
,”
ASME Appl. Mech. Rev.
,
72
, p.
040802
. 10.1115/1.4047707
21.
Leine
,
R.
, and
Nijmeijer
,
H.
,
2014
,
Dynamics and Bifurcations of Non-Smooth Mechanical Systems
.
Lecture Notes in Applied and Computational Mechanics
.
Springer
Berlin, Heidelberg
,
ISBN 9783642535796
.