Published: 30 September 2019

Vibrations of the nonlinear system in which stationary harmonic excited multivalued regimes in the vicinities of resonances do not exist

K. Ragulskis1
L. Ragulskis2
1Kaunas University of Technology, K. Donelaičio Str. 73, LT-44249, Kaunas, Lithuania
2Vytautas Magnus University, Vileikos Str. 8, LT-44404, Kaunas, Lithuania
Corresponding Author:
K. Ragulskis
Views 41
Reads 10
Downloads 1175

Abstract

A nonlinear dynamical system is investigated which consists from a mass between two linear elastic connecting elements with different coefficients of stiffness. Laws of vibrations and characteristics of eigenvibrations of the system as well as of self-decaying vibrations of the system with damping and of the system with harmonic excitation are determined. Dynamical qualities of the system are revealed. It is shown that the system has infinite number of eigenfrequencies and that in the resonance zones multivalued stable and unstable motions do not exist in the system.

1. Introduction

Basic concepts of investigation of vibrating systems are presented in [1]. Stabilisation of periodic nonlinear systems is investigated in [2]. Frequency analysis of a typical system is performed in [3]. Pendulum mechanism is investigated in [4]. Nonlinear vibrations of a piecewise linear model are analysed in [5]. Multiple resonant zones are investigated in [6]. Sommerfeld effect is analysed in [7]. Isolated resonances are investigated in [8]. Dynamics of vibromotors is analysed in [9].

This paper is dedicated for the investigation of a nonlinear system which does not possess multivalued regimes. Moreover, amplitude-frequency characteristics of such kind of system are linear. Such types of systems are important from the point of view of different engineering applications where stationary regimes do posses the stability in respect to different law of motions and energy regimes. The governing system of equations describing the investigated class of systems reads:

1
x¨+2h1x˙+p12x=fsinωt, x0,
2
x¨+2h2x˙+p22x=fsinωt, x0.

It is assumed that:

3
p1<p2,hi,pi, i=1,2, f=const, ˙=ddt.

2. The conservative system

In this case in the Eqs. (1), (2) it is assumed that:

4
hi=0, i=1,2, f=0.

According to the Eqs. (1), (4) when the initial conditions of motion are:

5
t=0:x=0, x˙=-x˙-, x˙->0,t=T-1:x=0, x˙=+x˙-,

it is obtained in the interval t0,T-1:

6
x=-x˙+psinp1t,x˙=-x˙+cosp1t.

From where it is obtained:

7
T-1=πp1.

Motion in the interval:

8
tT-1,T-1+T2,

takes place continuously, that is at t=T-1:x=0, x˙=x˙-. According to the Eqs. (2), (4):

9
x=x˙-p2sinp2t-T-1,x˙=x˙-cosp2t-T-1.

From the Eqs. (8), (9):

10
T-2=πp2,

and:

11
T-=πp1+p2p1p2=2πω-,

where the frequency of eigenvibrations ω-:

12
ω-=2p1p2p1+p2.

The whole solution by performing general notation x is obtained:

13
x=-x˙-p1sinp1t t0,T-1+x˙-p2sinp2t-T-1 tT-1,T-,
14
x˙=-x˙-cosp1t t0,T-1+x˙-cosp2t-T-1 tT-1,T-,
15
x¨=p1x˙-sinp1t t0,T-1-p2x˙-sinp2t-T-1 tT-1,T-.

Further x, x˙, x¨ by taking into account the Eqs. (13)-(15) expansion into the Fourier series with respect to ω- is performed:

16
x=n=0xnnω-t, x˙=n=0x˙nnω-t, x¨=n=0x¨nnω-t,

and the following relationships are obtained:

17
xa=n=0xna, x˙a=n=0x˙na, x¨a=n=0x¨na,

where xa is the maximum deviation from the medium value of the whole process, xna are the amplitudes of the respective harmonics. In the same way the expansions of x˙ and x¨ are performed.

According to the Eqs. (13)-(17) by assuming:

18
p1=1, p2=2,

and the initial conditions of motion at t=0:

19
x=0, x˙=-x˙-<0,

graphical relationships are obtained x=xt,x˙=x˙t,x¨=x¨t,x˙x¨=x˙x¨t,x˙=x˙x,x¨=x¨x˙, x˙x¨=x˙x¨x (see Fig. 1).

Fig. 1Dynamics of the system for the initial conditions of motion t=0, x0=0, x˙0= –1 (thin line), t=0, x0=0, x˙0= –2/3 (line of medium thickness) and t=0, x0=0, x˙0= –1/3 (thick line)

Dynamics of the system for the initial conditions of motion t=0, x0=0, x˙0= –1 (thin line), t=0, x0=0, x˙0= –2/3 (line of medium thickness) and t=0, x0=0, x˙0= –1/3 (thick line)

a) Displacement as function of time

Dynamics of the system for the initial conditions of motion t=0, x0=0, x˙0= –1 (thin line), t=0, x0=0, x˙0= –2/3 (line of medium thickness) and t=0, x0=0, x˙0= –1/3 (thick line)

b) Velocity as function of time

Dynamics of the system for the initial conditions of motion t=0, x0=0, x˙0= –1 (thin line), t=0, x0=0, x˙0= –2/3 (line of medium thickness) and t=0, x0=0, x˙0= –1/3 (thick line)

c) Acceleration as function of time

Dynamics of the system for the initial conditions of motion t=0, x0=0, x˙0= –1 (thin line), t=0, x0=0, x˙0= –2/3 (line of medium thickness) and t=0, x0=0, x˙0= –1/3 (thick line)

d) Velocity multiplied by acceleration as function of time

Dynamics of the system for the initial conditions of motion t=0, x0=0, x˙0= –1 (thin line), t=0, x0=0, x˙0= –2/3 (line of medium thickness) and t=0, x0=0, x˙0= –1/3 (thick line)

e) Phase trajectory: velocity as function of displacement

Dynamics of the system for the initial conditions of motion t=0, x0=0, x˙0= –1 (thin line), t=0, x0=0, x˙0= –2/3 (line of medium thickness) and t=0, x0=0, x˙0= –1/3 (thick line)

f) Phase trajectory: acceleration as function of velocity

Dynamics of the system for the initial conditions of motion t=0, x0=0, x˙0= –1 (thin line), t=0, x0=0, x˙0= –2/3 (line of medium thickness) and t=0, x0=0, x˙0= –1/3 (thick line)

g) Phase trajectory: velocity multiplied by acceleration as function of displacement

According to the Eq. (12) eigenfrequency of the system:

20
ω-=2121+2=43.

In the expansions into the Fourier series the first four members n=0,1,2,3 are taken into account.

By taking into account the Eqs. (11), (12) Table 1 is obtained.

According to the Eqs. (13)-(17) by taking into account Eqs. (18)-(20) graphical relationships are obtained (see Fig. 2).

Table 1Analysis of harmonics

n
0
1
2
3
n
nω-
0
43
243
343
n43
T-n=2πnω-
32π
322π
332π
3n2π

Fig. 2Amplitude frequency characteristics (constant part and first three harmonics)

Amplitude frequency characteristics (constant part and first three harmonics)

a) Displacement frequency characteristic

Amplitude frequency characteristics (constant part and first three harmonics)

b) Velocity frequency characteristic

Amplitude frequency characteristics (constant part and first three harmonics)

c) Acceleration frequency characteristic

Amplitude frequency characteristics (constant part and first three harmonics)

d) Velocity multiplied by acceleration frequency characteristic

3. Forced vibrations

According to the Eqs. (1), (2) by assuming:

21
p1=1,p2=2,h1=0.1,h2=0.2,fvar,

and at:

22
ω=34ω-,ω-,43ω-,53ω-,2ω-,52ω-,

steady state solutions are found and their spectral analysis ω, 2ω, 3ω, 4ω, ... is performed.

Steady state solutions are presented in Fig. 3, Fig. 4, Fig. 5, Fig. 6, Fig. 7, Fig. 8.

Fig. 3Steady state motions of the system for the first value of frequency of excitation and f= –1 (thin line), f= –2/3 (line of medium thickness), f= –1/3 (thick line)

Steady state motions of the system for the first value of frequency of excitation and  f= –1 (thin line), f= –2/3 (line of medium thickness), f= –1/3 (thick line)

a) Displacement as function of time

Steady state motions of the system for the first value of frequency of excitation and  f= –1 (thin line), f= –2/3 (line of medium thickness), f= –1/3 (thick line)

b) Velocity as function of time

Fig. 4Steady state motions of the system for the second value of frequency of excitation and f= –1 (thin line), f= –2/3 (line of medium thickness), f= –1/3 (thick line)

Steady state motions of the system for the second value of frequency of excitation and  f= –1 (thin line), f= –2/3 (line of medium thickness), f= –1/3 (thick line)

a) Displacement as function of time

Steady state motions of the system for the second value of frequency of excitation and  f= –1 (thin line), f= –2/3 (line of medium thickness), f= –1/3 (thick line)

b) Velocity as function of time

Fig. 5Steady state motions of the system for the third value of frequency of excitation and f= –1 (thin line), f= –2/3 (line of medium thickness), f= –1/3 (thick line)

Steady state motions of the system for the third value of frequency of excitation and  f= –1 (thin line), f= –2/3 (line of medium thickness), f= –1/3 (thick line)

a) Displacement as function of time

Steady state motions of the system for the third value of frequency of excitation and  f= –1 (thin line), f= –2/3 (line of medium thickness), f= –1/3 (thick line)

b) Velocity as function of time

Fig. 6Steady state motions of the system for the fourth value of frequency of excitation and f= –1 (thin line), f= –2/3 (line of medium thickness), f= –1/3 (thick line)

Steady state motions of the system for the fourth value of frequency of excitation and  f= –1 (thin line), f= –2/3 (line of medium thickness), f= –1/3 (thick line)

a) Displacement as function of time

Steady state motions of the system for the fourth value of frequency of excitation and  f= –1 (thin line), f= –2/3 (line of medium thickness), f= –1/3 (thick line)

b) Velocity as function of time

Fig. 7Steady state motions of the system for the fifth value of frequency of excitation and f= –1 (thin line), f= –2/3 (line of medium thickness), f= –1/3 (thick line)

Steady state motions of the system for the fifth value of frequency of excitation and  f= –1 (thin line), f= –2/3 (line of medium thickness), f= –1/3 (thick line)

a) Displacement as function of time

Steady state motions of the system for the fifth value of frequency of excitation and  f= –1 (thin line), f= –2/3 (line of medium thickness), f= –1/3 (thick line)

b) Velocity as function of time

Fig. 8Steady state motions of the system for the sixth value of frequency of excitation and f= –1 (thin line), f= –2/3 (line of medium thickness), f= –1/3 (thick line)

Steady state motions of the system for the sixth value of frequency of excitation and  f= –1 (thin line), f= –2/3 (line of medium thickness), f= –1/3 (thick line)

a) Displacement as function of time

Steady state motions of the system for the sixth value of frequency of excitation and  f= –1 (thin line), f= –2/3 (line of medium thickness), f= –1/3 (thick line)

b) Velocity as function of time

Amplitude frequency characteristics are presented in Fig. 9, Fig. 10, Fig. 11, Fig. 12, Fig. 13, Fig. 14.

Fig. 9Amplitude frequency characteristics (constant part and first three harmonics) for the first value of frequency of excitation

Amplitude frequency characteristics (constant part and first three harmonics)  for the first value of frequency of excitation

a) Displacement frequency characteristic

Amplitude frequency characteristics (constant part and first three harmonics)  for the first value of frequency of excitation

b) Velocity frequency characteristic

Fig. 10Amplitude frequency characteristics (constant part and first three harmonics) for the second value of frequency of excitation

Amplitude frequency characteristics (constant part and first three harmonics)  for the second value of frequency of excitation

a) Displacement frequency characteristic

Amplitude frequency characteristics (constant part and first three harmonics)  for the second value of frequency of excitation

b) Velocity frequency characteristic

Fig. 11Amplitude frequency characteristics (constant part and first three harmonics) for the third value of frequency of excitation

Amplitude frequency characteristics (constant part and first three harmonics)  for the third value of frequency of excitation

a) Displacement frequency characteristic

Amplitude frequency characteristics (constant part and first three harmonics)  for the third value of frequency of excitation

b) Velocity frequency characteristic

Fig. 12Amplitude frequency characteristics (constant part and first three harmonics) for the fourth value of frequency of excitation

Amplitude frequency characteristics (constant part and first three harmonics)  for the fourth value of frequency of excitation

a) Displacement frequency characteristic

Amplitude frequency characteristics (constant part and first three harmonics)  for the fourth value of frequency of excitation

b) Velocity frequency characteristic

Fig. 13Amplitude frequency characteristics (constant part and first three harmonics) for the fifth value of frequency of excitation

Amplitude frequency characteristics (constant part and first three harmonics)  for the fifth value of frequency of excitation

a) Displacement frequency characteristic

Amplitude frequency characteristics (constant part and first three harmonics)  for the fifth value of frequency of excitation

b) Velocity frequency characteristic

Fig. 14Amplitude frequency characteristics (constant part and first three harmonics) for the sixth value of frequency of excitation

Amplitude frequency characteristics (constant part and first three harmonics)  for the sixth value of frequency of excitation

a) Displacement frequency characteristic

Amplitude frequency characteristics (constant part and first three harmonics)  for the sixth value of frequency of excitation

b) Velocity frequency characteristic

4. Conclusions

On the basis of the presented results the qualities of dynamic behavior of the nonlinear vibrating system which consists from a mass between two linear elastic connecting members with different coefficients of stiffness can be used in different engineering applications requiring the stability in a wide range of parameters.

The motions of the system and amplitude frequency characteristics are determined by analytical – numerical method relationships. It is determined that for the case of the conservative system eigenfrequencies do not depend on the value of the amplitude of excitation. It is shown that for the case of forced harmonic excitation stable and unstable multivalued regimes do not exist in the system. Also, separate classes of nonlinear dynamical systems are investigated which have qualities of similar type.

The presented results enable to perform the design of nonlinear vibrating systems of this type. This enables to create new systems which can be noted by high stability properties, because in the vicinities of resonances they do not have multivalued solutions.

References

  • Bolotin V. V. Vibrations in Engineering. Handbook, Vol. 1, Mashinostroienie, Moscow, 1978, (in Russian).
  • Zaitsev V. A. Global asymptotic stabilization of periodic nonlinear systems with stable free dynamics. Systems and Control Letters, Vol. 91, 2016, p. 7-13.
  • Salahshoor E., Ebrahimi S., Zhang Y. Frequency analysis of a typical planar flexible multibody system with joint clearances. Mechanism and Machine Theory, Vol. 126, 2018, p. 429-456.
  • Starossek U. Forced response of low-frequency pendulum mechanism. Mechanism and Machine Theory, Vol. 99, 2016, p. 207-216.
  • Wang S., Hua L., Yang C., Zhang Y., Tan X. Nonlinear vibrations of a piecewise-linear quarter-car truck model by incremental harmonic balance method. Nonlinear Dynamics, Vol. 92, 2018, p. 1719-1732.
  • Alevras P., Theodossiades S., Rahnejat H. On the dynamics of a nonlinear energy harvester with multiple resonant zones. Nonlinear Dynamics, Vol. 92, 2018, p. 1271-1286.
  • Sinha A., Bharti S. K., Samantaray A. K., Chakraborty G., Bhattacharyya R. Sommerfeld effect in an oscillator with a reciprocating mass. Nonlinear Dynamics, Vol. 93, 2018, p. 1719-1739.
  • Habib G., Cirillo G. I., Kerschen G. Isolated resonances and nonlinear damping. Nonlinear Dynamics, Vol. 93, 2018, p. 979-994.
  • Ragulskis K., Bansevičius R., Barauskas R., Kulvietis G. Vibromotors for Precision Microrobots. Hemisphere, New York, 1987.

About this article

Received
15 July 2019
Accepted
20 August 2019
Published
30 September 2019
Keywords
nonlinear system
coefficients of stiffness
amplitude-frequency characteristics
dynamical qualities