Mathematical Problems in Engineering

Volume 2008 (2008), Article ID 471080, 25 pages

http://dx.doi.org/10.1155/2008/471080

## The Effect of Infinitesimal Damping on the Dynamic Instability Mechanism of Conservative Systems

^{1}Department of Civil Engineering, University of Thessaly, Pedion Areos, 38 334 Volos, Greece^{2}Laboratory of Metal Structures, Department of Civil Engineering, National Technical University of Athens, Zografou Campus, 157 80 Athens, Greece^{3}Research Center Pure and Applied Mathematics, Academy of Athens, Soranou Efessiou 4, 115 27 Athens, Greece

Received 17 December 2007; Accepted 19 February 2008

Academic Editor: Jose Balthazar

Copyright © 2008 Dimitris S. Sophianopoulos et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

The local instability of 2 degrees of freedom (DOF) weakly damped systems is thoroughly discussed using the Liénard-Chipart stability criterion. The individual and coupling effect of mass and stiffness distribution on the dynamic asymptotic stability due to mainly infinitesimal damping is examined. These systems may be as follows: (a) unloaded (free motion) and (b) subjected to a suddenly applied load of constant magnitude and direction with infinite duration (forced motion). The aforementioned parameters combined with the algebraic structure of the damping matrix (being either positive semidefinite or indefinite) may have under certain conditions a tremendous effect on the Jacobian eigenvalues and then on the local stability of these autonomous systems. It was found that such systems when unloaded may exhibit periodic motions or a divergent motion, while when subjected to the above step load may experience either a degenerate Hopf bifurcation or periodic attractors due to a generic Hopf bifurcation. Conditions for the existence of purely imaginary eigenvalues leading to global asymptotic stability are fully assessed. The validity of the theoretical findings presented herein is verified via a nonlinear dynamic analysis.

#### 1. Introduction

In previous studies of the 3rd author, based on 2-DOF and 3-DOF cantilevered models [1] under partial follower loading
(nonconservative systems), it was shown that in a small region of *divergence **instability*,
flutter (dynamic instability) may occur before divergence (static instability),
if very small damping is included [2, 3]. Bolotin et al. [4] using an aeroelastic
model presented a similar result. Païdoussis et al. [5] and Païdoussis [6] have
shown that flutter may occur in an inherently conservative system but for large
damping. However, the effect of damping, being of paramount importance in
nonconservative autonomous systems, was in general ignored when these systems
are subjected to a step conservative (potential) loading. This is so because it
was widely accepted that dynamic stability in nondissipative conservative
systems, which are *stable*, does not
change by the inclusion of damping [7].

The local dynamic stability of discrete systems under
step conservative loading when small dissipative forces are included is
governed by the matrix-vector differential equation [8–11]: where the dot denotes a
derivative with respect to time *t*; is an *n*-dimensional *state* vector with coordinates
and are real symmetric matrices.
More specifically, matrix **M** associated with the total kinetic energy of the system is a function of the
concentrated masses being always *positive definite*; matrix **C** the elements of which are the damping coefficients
may be *positive definite*, positive *semidefinite* as in the case of pervasive
damping [12, 13], or *indefinite* [14–16]; **V** is a generalized *stiffness* matrix with coefficients whose elements are also linear functions of a *suddenly* applied external load with
constant direction and infinite duration [17], that is, .
Apparently, due to this type of loading, the system under discussion is *autonomous*. When the external loading
is applied *statically,* one can obtain
the static (divergence) instability or buckling loads by vanishing of the determinant of the stiffness
matrix ,
that is, Clearly, (1.2) yields an *n*th degree algebraic
equation in . Assuming *distinct* critical
states, the matrix is *positive **definite* for *positive
semidefinite* for and *indefinite* for

Kounadis in two very recent publications [10, 11] has
established the conditions under which the above autonomous dissipative systems
under step (conservative) loading may exhibit dynamic bifurcational modes of
instability *before* divergence (static)
instability, that is, for , when *infinitesimal* damping is included. These
bifurcational modes may occur through either a *degenerate* Hopf bifurcation (leading to periodic motion around
centers) or a *generic* Hopf
bifurcation (leading to periodic attractors or to flutter). These unexpected
findings (implying failure of Ziegler’s kinetic criterion and other singularity
phenomena) may occur for a certain combination of values of the mass
(primarily) and stiffness distribution of the system in connection with a *positive semidefinite* or an *indefinite* damping matrix.

The question which now arises is whether there are
combinations of values of the above-mentioned parameters (mass and stiffness
distribution) which in connection with the algebraic structure of damping matrices
may lead to dynamic bifurcational modes of instability when the system under
discussion is *unloaded*. Such local
(due to unforced motion) dynamic instability will be sought through the set of
asymptotic stability criteria of Liénard-Chipart [8, 18] which are elegant and more
readily employed than the well-known Routh-Hurwitz stability
criteria.

As another main objective of this work, some new dynamic bifurcations related to the algebraic structure of the damping matrix when the systems are loaded by the above type of step conservative load will be also discussed, using the Liénard-Chipart criterion by analyzing 2-DOF models for which a lot of numerical results are available. Finally, the conditions of a double purely imaginary root leading to a new dynamic bifurcation, whose response is similar to that of a generic Hopf bifurcation, are properly established.

#### 2. Basic Equations

Solution of (1.1) can be sought in
the form where is in general a complex number and **r** is a complex vector independent of
time *t*.

Introducing **q** from
(2.1) into (1.1), we
get

For given stiffness coefficients , the
generalized stiffness matrix **V** is a linear function of . Thus, if matrices **M**, **C**, **V** are given, solutions of (2.2) are intimately related to the algebraic
properties of the matrix-valued function , and more specifically to the Jacobian
eigenvalues obtained through the vanishing of the determinant: whose expansion gives the characteristic (secular) equation
for an N-DOF system: where the real
coefficients are determined by means of Bôcher formula [19].
The eigenvalues (roots) of (2.4) are, in general, *complex conjugate* pairs (where
and are real numbers and ) with corresponding *complex conjugate* eigenvectors and . Since ,
clearly and .
Thus, the solutions of (1.1) are of the
form where *A* and *B* are *constants* which are determined from the *initial* conditions. Solutions in (2.5) are *bounded*,
tending to zero as ,
if all eigenvalues of (2.4) have *negative* real parts, that is, when for all *j*. In this case, the algebraic
polynomial (2.4) is called a *Hurwitz
polynomial* (since all its roots have negative real parts) and the *origin * of the system is *asymptotically stable*.

##### 2.1. Criteria for Asymptotic Stability

The *necessary and sufficient* conditions which assure that all roots of (2.4) have negative real parts (i.e., for all *j*) which means that the
corresponding polynomial is a Hurwitz
polynomial are of great practical importance.

Consider the more general case of a polynomial in *z*
with real coefficients : for which we will seek the *necessary and sufficient* conditions so that all its roots have *negative* real parts.

Denoting
by the real roots and by the complex roots of (2.6), we may assure that all these roots in the
complex plane lie to the *left* of the
imaginary axis, that is, Then one can write

Since due to inequality (2.7), each term in the last
part of (2.8) has *positive* coefficients, it is deduced that *all* coefficients
of (2.6) are also *positive*. However,
this (i.e.,
for all *i* with )
is a *necessary* but by no means *sufficient* condition for all roots of (2.6) to lie in the left half-plane (i.e., ).

According to *Routh-Hurwitz* criterion [18]
of asymptotic stability for all roots of (2.6) to
have *negative* real parts, the *necessary* and *sufficient* conditions are where (with for ). Note the last equality

It should be noted that when the above necessary
conditions (for all *i*) hold, inequalities (2.9) are *not* independent. For instance, for , the Routh-Hurwitz
conditions reduce to the single inequality ; for , they reduce to and ; while for , they reduce also to two inequalities, .
This case was discussed by Liénard and Chipart who established the following
elegant criterion for asymptotic stability [8].

*The Liénard-Chipart Stability Criterion*

For a polynomial with real
coefficients to have all roots with *negative* real parts, it is *necessary **and sufficient* that

(1)all coefficients of *f*(*z*) be positive, that is,
(2)the determinant inequalities be also
positive, that is,
where denotes as before the Hurwitz
determinant of th order.

It can be shown that if the Hurwitz determinants of *odd* order are *positive*, then
those of *even* order are also *positive*, and vice versa. This holds
even when *only part* of the
coefficients of *f*(*z*) (with are positive. According to this, the Liénard-Chipart criterion
is defined as follows.*Necessary* and *sufficient* conditions for all roots of the real polynomial to have *negative
real parts* can be given in any one of the following forms [18]:

(1)(2)
This stability criterion was rediscovered by Fuller [20].

For instance, for a 2-DOF cantilevered
model, the characteristic (secular) (2.4) is

According to the last criterion,
all roots of (2.13) have negative real parts provided that
and . Clearly, from the last inequality, it follows that .
Hence, the positivity of and was assured via the above conditions (i.e., ).

#### 3. Mathematical Analysis

Subsequently, using the spring
cantilevered dynamical model of 2-DOF shown in Figure 1, we will examine in
detail the *effect of violation* of one
or more of the conditions of Liénard-Chipart criterion
on its asymptotic stability. The response of this dynamic model carrying two
concentrated masses is studied when it is either *unloaded* or *loaded* by a
suddenly applied load of constant magnitude and direction with infinite
duration. Such autonomous dissipative systems with infinitesimal damping (including
the case of zero loading) are properly discussed. If at least one root of the
secular equation (2.13) has a positive real part, the corresponding solution (2.5)
will contain an exponentially increasing function and the system will become
unstable.

The seeking of an *imaginary root* of the secular equation (2.13)
which represents a border line between dynamic stability and instability is a
first but important step in our discussion. Clearly, an imaginary root gives
rise to an oscillatory motion of the form around the trivial state. However, the existence
of at least one *multiple imaginary root* of
the th order of multiplicity leads to a solution containing functions of
the form which increases with time. Hence, the *multiple* imaginary root on the imaginary
axis of the complex plane denotes *local dynamic
instability*. The discussion of such a situation is also another objective
of this study.

The nonlinear equations of motion for the 2-DOF model of Figure 1 with rigid links of equal length are given by [11] where .

Linearization of (3.1) after setting gives where

Note that in case of a Rayleigh’s dissipative
function the damping coefficients are, , and ,
where is dimensionless coefficient for the *i*th rigid link. This case (for which ) is a specific situation of the damping matrix **C** which is not discussed herein.

The *static* buckling (divergence) equation is given by whose lowest root is the *first* buckling load equal to Clearly, for the entire
interval of values of , (3.6)
yields .

The characteristic (secular) equation is where

Let us first examine the effect
of violation of Liénard-Chipartcriterion on the system stabilityin
the case of *zero loading* (i.e., ). Then expressions in (3.8) due to relations (3.4) are written as
follows:
According to Liénard-Chipart criterion,
inequalities (2.12a) imply where and .
Clearly, from the last inequality, it follows that .

For , and (implying ), it is deduced that this criterion is *violated* if either one of is zero or is zero. These three cases will be discussed
separately in connection with the algebraic structure of the damping matrix

*Case 1 (). *If (yielding ), then Equation (3.11),
being independent of and
*k*, is satisfied only when the damping matrix **C** is *indefinite*, that is, Indeed, the last inequality due to relation
(3.11) implies which is always satisfied,
regardless of the value of *c _{22}*/

*c*, since for , the discriminant of (3.13) (equal to ) is

_{11}*always*

*negative*.

Thus, we have explored the

*unexpected finding*that an

*unloaded*(stable) system becomes

*dynamically unstable*at any small disturbance in case of an

*indefinite*damping matrix even when

*infinitesimal*damping is included.

Since all coefficients of (3.7) are positive from the

*theory of algebraic equations*it follows that this equation

*cannot*have positive root. Also the case of existence of a pair of pure imaginary roots associated with is ruled out, since (due to ). Hence, (3.7) has either two

*negative*roots combined with a pair of complex conjugate roots with

*positive*real part or two pairs of

*complex conjugate*roots with

*opposite*real parts. Both cases imply local dynamic instability.

*Case 2 (). *If (implying also ), then Namely, the damping matrix
is *indefinite* but with large negative
determinant (rather unrealistic case). Since the Liénard-Chipart criterion
is violated, the model is again locally dynamically unstable.

Since all coefficients of (3.7) are *positive*, from
the *theory of algebraic equations*, it is deduced that this equation
cannot have *positive* root. Also the
case of existence of a pair of pure imaginary roots associated with is ruled out, since (due to ). Hence, (3.7) has either *two negative* roots combined with a pair of complex conjugate roots
with *positive* real part or two pairs
of *complex conjugate* roots with *opposite* real parts. Both cases imply local
dynamic instability.

*Case 3 (). *In
this case, stability conditions in (3.10) are satisfied except for the last one, since which yields [11]
Note that implies (i.e., ).

This is a necessary condition for
the secular (3.7) to have one pair of *pure
imaginary* roots .
Indeed, this can be readily established by inserting into (3.7) and then equating to zero real and
imaginary parts.

Consider
now the more general case of *nonzero
loading* (i.e., ). Using
relations in (3.8),
(3.15) can be written as follows: where
where
For real ,
the discriminant **D** of
(3.16) must be greater than
or equal to zero, that is, Subsequently, attention is focused
on the following: (a)
matrix **C** is *positive semidefinite* (i.e., with )
and (b) matrix **C** is *indefinite* ( with ).

Using the symbolic manipulation of *Mathematica* [21], one can find that
where is an algebraic polynomial of 5th degree in

*Case 4 (). *For , (3.16) implying
admits a *double root*, which due to (3.17a), (3.17b), (3.17d) is given
by Using the *Reduce* command
embedded in *Mathematica*, one can find
the conditions under which ,
given in the appendix, relation (A.1).

*Case 5 (). *Moreover, it was found symbolically that the 5th degree polynomial possesses *three
real* roots (one double and one single), and two *pure imaginary* ones. Discussing their nature, one can find that the
double root of ,
being equal to , yields Then, the double root of (3.16) becomes which is always greater than and hence
of minor importance for the present analysis.

The third real root of ,
if substituted in (3.16), yields again a double root in ,
always less than zero, which is rejected. Thus, only the case of a *positive semidefinite* damping matrix may
lead to an acceptable value of the corresponding load (i.e., ) associated with a *degenerate* Hopf bifurcation, as theoretically was shown by Kounadis
[10, 11].

*Case 6 (). * If ,
(3.16) implies ,
which after symbolic manipulation of (3.17c) can be written in the following
form: where are given in the appendix, relations (A.2).
It is evident that and , a fact implying that (3.23)
can be satisfied only for if also ; otherwise (i.e.,
if ) the system may be locally stable or unstable. For ,
one can find the corresponding values of , given in (A.3)
and (A.4) in the appendix, which are always positive. This special case, for
which the trivial (unloaded) state is associated with a pair of *pure imaginary* eigenvalues (necessary condition
for a Hopf bifurcation), implies a local dynamic instability.

*Conditions for A Double Imaginary Root*

For a *double imaginary* root, the first
derivative of the secular equation (3.7) must be also zero, which yields Inserting into (3.24) we obtain and thus Introducing this expression of into (3.15), it follows that ,
which also implies that and hence (3.24) becomes .
If is inserted into the secular equation ,
for a *double imaginary* root, it
follows that which due to relations (3.8) yields
For real ,
the discriminant **D** of (3.25) must satisfy the inequality which yields Using the conditions found above relations (3.9) yield Adding the last two equations, we
obtain Since and it follows that both coefficients of *c _{11}* and

*c*in (3.30) are positive. Hence,

_{22}*c*and

_{11}*c*are of

_{22}*opposite*sign (i.e., ) and consequently ; thus the 2nd of inequalities (3.27) is excluded.

Solving simultaneously the system
of equations , in *k*, *m*, , two
ternaries of values for *k*, *m*, and are obtained, given in the appendix, expressions (A.5). For all these values to be greater than zero, the *Reduce* command embedded in Mathematica
[21] yields two sets of conditions, given also in the appendix, relations
(A.6). Further symbolic computations are needed for establishing the
conditions for *a double pure imaginary* root for loading values less than .
Nevertheless, suitable combinations of values of , *k*, and *m* may
be found. This will be demonstrated in Section 4. The corresponding dynamic
response, since the system is associated with a codimension-2 local bifurcation, is anticipated to
be related to isolated periodic orbits
which will be established via a straightforward complete nonlinear dynamic
analysis.

#### 4. Numerical Results

In this section, numerical results corresponding to all the above cases of violation of the Liénard-Chipart stability criterion are given below in tabular and graphical forms.

*Case 1 (). *(a) For an unloaded system with , choosing and ,
(3.11) yields and as expected the damping matrix is *indefinite* with determinant The two pairs of corresponding eigenvalues are and = 0.00335877
,
implying *local* dynamic instability.
Solving numerically the system of nonlinear equations (3.1), we find that the
dynamic response of the system is associated with a divergent motion, as
depicted in Figure 2, with the aid of the time series
, time velocities and
phase-plane portraits
.

(b)
. For a system with , and and for ,
(2.13) yields implying The trivial state is locally dynamically unstable,
since and The corresponding dynamic response is again
related to a *divergent* (unbounded) motion,
as shown in the phase-plane portraits of Figures 3(a) and 3(b).

*Case 2 (). *(a) : If
relation (3.14) is satisfied, for example, for the damping coefficients
, yielding and .
For this, rather unrealistic, subcase, the corresponding eigenvalues are equal
to and (local instability). Hence, the response of
the system is also related to a *divergent* (unbounded) motion, presented graphically in the phase-plane curves of
Figures 4(a) and 4(b).

(b) .
Similarly, for a system with (for which ) in order that we must choose an indefinite damping matrix
with .
Setting, for example, , and the trivial state is locally dynamically unstable
with and The system exhibits a *divergent* (unbounded) motion, as shown in Figures 5(a) and 5(b).

*Case 3 (). *(a) *Positive
semidefinite damping matrix *. Choosing
(and thus ), the 1st requirement of the 2nd set of conditions given in the appendix, relation (A.1), is satisfied (i.e., ).
The 2nd requirement, that is, yields , and hence one can choose . The 3rd requirement in (A.1) implies that and thus
one can take . Solving numerically (3.15) with respect to , we
obtain the value of , associated with a pair of
pure imaginary eigenvalues, while the other pair has negative real parts. The
evolution of both pairs of eigenvalues in the complex plane as
varies is presented in Figure 6 for .
For , *a degenerate Hopf bifurcation
occurs* and the system exhibits a periodic motion, whose amplitude depends
on the initial conditions. Relevant results in graphical form can be found in
recent publications of the 3rd author [10, 11].

(b) * Indefinite damping matrix *.
It has been proven by Kounadis [10, 11] that in this subcase (for ) all the necessary and sufficient conditions
for a *generic Hopf bifurcation* are fulfilled and hence the
system experiences a periodic attractor response (stable limit cycles) with
constant final amplitudes regardless of the initial conditions. Numerical
results are given in the aforementioned papers by Kounadis.

(c) * and *. If at the same time , one can find the values of
through (A.3) and (A.4) in the appendix, which are
always positive. A further investigation of this case as well as of the case
for the global stability of the system can be performed through a
nonlinear dynamic analysis.

(d)
*Double pure imaginary eigenvalues*. For this
special case, three combinations of damping matrix coefficients are examined. These, along with the corresponding critical values of k, and
*m*, satisfying relations (A.5) of the appendix, are given in Table 1. Note
that Cases 3(d)1 and 3(d)2 are the outcome of the 1st set from relations
(A.5), while Case 3(d)3 is the outcome from the 2nd set. Clearly, in all cases, .

In the three above subcases, the evolution of both pairs of -dependent eigenvalues
in the complex plane is depicted in Figures 7, 8, 9(a), and 9(b), from which it
is evident that for all ,
except for (where a codimension-2 bifurcation
occurs), the pairs of eigenvalues remain always in *opposite* planes of the Im axis, but symmetric with respect to the Re
axis. This symmetry is always present for the pair with negative real parts,
while for the other pair (with positive real parts), this feature remains until
their imaginary part vanishes simultaneously at a certain value of the loading
less than .

The dynamic response of the system for all these subcases is associated with isolated periodic orbits (whose final amplitude is constant and independent of the initial conditions), as shown in the phase-plane trajectories of Figures 10, 11, and 12.

The corresponding dynamic bifurcations related to the above *double pure imaginary* eigenvalues behave like a *generic* Hopf bifurcation, whose basic
feature is the intersection of the -dependent path of
one eigenvalue with the imaginary axis. On the other hand, in all the above
subcases, the branches of two consecutive -dependent
eigenvalues meet the imaginary axis at the same point with .
Namely, the *transversality* condition
is satisfied through two intersected lines at the same point of the imaginary
axis, but whose branches in the left (negative) and right (positive) half planes
belong to the 1st and 2nd pairs of eigenvalues, respectively.

Finally, Figure 13 verifies
the unexpected phenomenon (Kounadis [11]) of *discontinuity* in the dynamic loading associated with either a *degenerate* or a *generic* Hopf bifurcation.

#### 5. Concluding Remarks

This
study discusses in detail the coupling effect of (infinitesimal mainly) damping
with the mass and stiffness distribution in a 2-DOF cantilevered model under
step potential loading. Such an autonomous system may be associated either with
a positive semidefinite or indefinite damping matrix (with positive or negative
diagonal elements). Attention is focused on the violation of the Liénard-Chipart stability criterion
when this system is either *unloaded* or *loaded* by a suddenly applied load
of constant magnitude and direction with infinite duration (step potential
loading). The most important findings of this study are the following.

(1) Usage of Liénard-Chipart, simple and readily employed, stability criterion compared to that of Routh-Hurwitz brought into light a variety of new types of dynamic bifurcations reported below.

(2) The mass and stiffness distribution combined
either with a positive semidefinite or an indefinite damping matrix may have a
considerable effect on the *asymptotic
stability*, *prior* to divergence
instability, even in case of *infinitesimal* damping.

(3) The cantilevered model when *unloaded* (being statically stable),
strangely enough, under certain conditions becomes dynamically unstable to any
small disturbance leading to a *divergent* (unbounded) motion.

(4) The above model when *loaded* under analogous of the previous conditions exhibits also a *divergent* motion at a certain value of
the external load.

It is worth noting that both the above cases of divergent motion may occur for *negligibly small* negative determinant of
the damping (indefinite) matrix when , while for
(regardless of whether
or ), the determinant of the damping matrix is
negative but *finite*.

(5) When and four distinct responses may occur as follows.

(a)If (positive semidefinite damping
matrix), the system exhibits a periodic motion associated with a *degenerate* Hopf bifurcation. Then, the
final amplitude of motion depends on the initial conditions.(b)When the damping matrix is *indefinite* (even with infinitesimal but
negative determinant, ), the system may exhibit a periodic attractor
response associated with a *generic* Hopf bifurcation. It is of paramount practical importance the case where such
an unexpected dynamic instability occurs at a load less than the 1st buckling load.
In both the above cases,
it was confirmed the recently reported unexpected finding [11] of discontinuity
of the dynamic loading (associated with either a
degenerate or a generic Hopf bifurcation) occurring at a certain value of the *mass* distribution.(c)When at the same time , we
have a *local* dynamic instability for
, whose global stability can be established through a nonlinear dynamic analysis.(d)The case of a *double pure imaginary eigenvalue* may occur for an *indefinite* damping matrix with *finite* determinant and *negative ratio* of the corresponding
diagonal elements. In this special case, there are *two pairs* of eigenvalues in the complex plane which touch the
imaginary axis at the same point for a certain value . This
situation yields local instability leading to a motion with final constant
amplitude regardless of the initial conditions. Namely, such a dynamic
bifurcation behaves in a way similar to that of a *generic* Hopf bifurcation. This new type of dynamic bifurcation was
also verified via a nonlinear dynamic analysis.

#### Appendix

#### Results of Symbolic Computations

(i) Conditions for from (3.20), with are

(ii) Expressions of coefficients of (3.23) are and corresponding values of for are

(iii) Values of , *k*, and *m* for which the
secular equation (3.7) has a double pure imaginary pair of roots (eigenvalues) are
and the corresponding conditions for the above two sets of values to
be positive are as
follows.

For the 1st set (), For the 2nd set (), where

#### References

- H. Ziegler, “Die Stabilitätskriterien der Elastomechanik,”
*Ingenieur Archiv*, vol. 20, no. 1, pp. 49–56, 1952. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. N. Kounadis, “On the failure of static stability analyses of nonconservative systems in regions of divergence instability,”
*International Journal of Solids and Structures*, vol. 31, no. 15, pp. 2099–2120, 1994. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. N. Kounadis, “Non-potential dissipative systems exhibiting periodic attractors in regions of divergence,”
*Chaos, Solitons & Fractals*, vol. 8, no. 4, pp. 583–612, 1997. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - V. V. Bolotin, A. A. Grishko, and A. V. Petrovsky, “Secondary bifurcations and global instability of an aeroelastic non-linear system in the divergence domain,”
*Journal of Sound and Vibration*, vol. 191, no. 3, pp. 431–451, 1996. View at Publisher · View at Google Scholar - M. P. Païdoussis, D. Mateescu, and W.-G. Sim, “Dynamics and stability of a flexible cylinder in a narrow coaxial cylindrical duct subjected to annular flow,”
*Journal of Applied Mechanics*, vol. 57, no. 1, pp. 232–240, 1990. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - M. P. Païdoussis,
*Fluid-Structure Interactions: Slender Structures and Axial Flow*, vol. 2, Elsevier Academic Press, London, UK, 2003. View at Zentralblatt MATH · View at MathSciNet - K. Huseyin,
*Multiple Parameter Stability Theory and Its Applications*, vol. 18 of*Oxford Engineering Science Series*, Clarendon Press, Oxford, UK, 1986. View at Zentralblatt MATH · View at MathSciNet - F. R. Gantmacher,
*Lectures in Analytical Mechanics*, Mir, Moscow, Russia, 1970. View at Zentralblatt MATH · View at MathSciNet - K. Huseyin,
*Vibrations and Stability of Damped Mechanical Systems*, Sijthoff & Noordhoff, Alphen aan den Rijn, The Netherlands, 1978. View at Zentralblatt MATH · View at MathSciNet - A. N. Kounadis, “Hamiltonian weakly damped autonomous systems exhibiting periodic attractors,”
*Zeitschrift für Angewandte Mathematik und Physik*, vol. 57, no. 2, pp. 324–350, 2006. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. N. Kounadis, “Flutter instability and other singularity phenomena in symmetric systems via combination of mass distribution and weak damping,”
*International Journal of Non-Linear Mechanics*, vol. 42, no. 1, pp. 24–35, 2007. View at Publisher · View at Google Scholar · View at MathSciNet - E. E. Zajac, “The Kelvin-Tait-Chetaev theorem and extensions,”
*Journal of the Astronautical Sciences*, vol. 11, no. 2, pp. 46–49, 1964. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - E. E. Zajac, “Comments on “stability of damped mechanical systems and a further extension”,”
*AIAA Journal*, vol. 3, no. 9, pp. 1749–1750, 1965. View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - R. Sygulski, “Dynamic stability of pneumatic structures in wind: theory and experiment,”
*Journal of Fluids and Structures*, vol. 10, no. 8, pp. 945–963, 1996. View at Publisher · View at Google Scholar - A. Laneville and A. Mazouzi, “Wind-induced ovalling oscillations of cylindrical shells: critical onset velocity and mode prediction,”
*Journal of Fluids and Structures*, vol. 10, no. 7, pp. 691–704, 1996. View at Publisher · View at Google Scholar - A. K. Misra, S. S. T. Wong, and M. P. Païdoussis, “Dynamics and stability of pinned-clamped and clamped-pinned cylindrical shells conveying fluid,”
*Journal of Fluids and Structures*, vol. 15, no. 8, pp. 1153–1166, 2001. View at Publisher · View at Google Scholar - A. N. Kounadis, “A geometric approach for establishing dynamic buckling loads of autonomous potential two-degree-of-freedom systems,”
*Journal of Applied Mechanics*, vol. 66, no. 1, pp. 55–61, 1999. View at Publisher · View at Google Scholar - F. R. Gantmacher,
*The Theory of Matrices*, Chelsea, New York, NY, USA, 1959. View at Zentralblatt MATH · View at MathSciNet - L. A. Pipes and L. R. Harvill,
*Applied Mathematics for Engineers and Physicists*, International Student Edition, McGraw-Hill, Kogakusha, Tokyo, Japan, 3rd edition, 1970. View at Zentralblatt MATH · View at MathSciNet - A. T. Fuller, “Conditions for a matrix to have only characteristic roots with negative real parts,”
*Journal of Mathematical Analysis and Applications*, vol. 23, no. 1, pp. 71–98, 1968. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - T. B. Bahder,
*Mathematica for Scientists and Engineers*, Addison-Wesley, Reading, Mass, USA, 1995. View at Zentralblatt MATH · View at MathSciNet