Determination of the first natural frequency of an elastic rod of a discretecontinuous vibratory system
Oleksii Lanets^{1} , Oleksandr Kachur^{2} , Vitaliy Korendiy^{3} , Petro Dmyterko^{4} , Sergij Nikipchuk^{5} , Iryna Derevenko^{6}
^{1, 2, 3, 4}Department of Robotics and Integrated Mechanical Engineering Technologies, Lviv Polytechnic National University, Lviv, Ukraine
^{5}Department of Motor Vehicle Transport, Lviv Polytechnic National University, Lviv, Ukraine
^{6}Department of Strength of Materials and Structural Mechanics, Lviv Polytechnic National University, Lviv, Ukraine
^{3}Corresponding author
Vibroengineering PROCEDIA, Vol. 37, 2021, p. 712.
https://doi.org/10.21595/vp.2021.21981
Received 1 April 2021; received in revised form 10 April 2021; accepted 18 April 2021; published 7 May 2021
A continuous rodshaped member (a body with distributed inertia and rigidity parameters), which is the object of the investigations, is considered. To ensure the optimal natural frequency concerning the oscillations of the reactive (exciting) mass of the threemass discrete vibratory system, with the use of the KrylovDuncan functions, the mathematical model describing forced oscillations of the continuous member considered as a disturbing body of the threemass discretecontinuous vibratory system is established, and the corresponding frequency equation is analytically derived. The obtained theoretical results are verified using the Finite Element Method in SolidWorks software. The novelty of the present paper consists in substantiation of the possibilities of implementing the continuous rodshaped members with distributed inertia and rigidity parameters for exciting the oscillations of the threemass discretecontinuous interresonance vibratory systems.
 The mathematical model describing forced oscillations of the continuous member considered as a disturbing body of the threemass discretecontinuous vibratory system is established.
 The frequency equation and the analytical expression for calculating the first natural frequency of the continuous member (elastic rod) of the discretecontinuous vibratory system are derived.
 Frequency analysis of a rodshaped continuous body is performed using Finite Elements Method in SolidWorks software; the obtained results are in satisfactory agreement with analytical calculations.
Keywords: oscillations, discretecontinuous system, continuous member, elastic rod, KrylovDuncan functions, inertial parameters, Finite Element Method.
1. Introduction
Investigations of resonance and transresonance regimes of work and determination of dynamical characteristics of oscillatory systems are the most important issues in designing the vibratory technological equipment. As a rule, in such systems, the working members are bodies whose inertial parameters are distributed over planes of given configurations and geometric sizes. Plates, beams, etc. belong to them. In the work [1], there are considered linear and nonlinear parameters of the main frequency of coneshaped beams which had been calculated by means of the coupled displacement field method. Transverse vibrations of toughelastic beams are determined by means of Kelvin model with different parameters of toughness [2]. The governing equations and boundary conditions are formed according to generalized Hamilton principle.
The application of ANSYS software in order to calculate the first natural frequency of a plate and to determine the shapes of oscillations simplifies the body of mathematics to some extent [3, 4]. In the work [5], vibrational analysis of sandwichpanels was carried out; the equations of motion were there obtained according to Hamilton principle, and the solution was carried out by means of the method of Generalized Differential Quadratures. The modelling was conducted in SolidWorks software. In the paper [6], there were analyzed Tshaped plates, where the prior analytical calculations were carried out according to the RayleighRitz method, which allowed for determination of the natural frequency of the plate. In the work [7], a beam as a homogeneous rod with constant inertia and rigidity parameters along its length is considered. The paper [8] considers twomass resonance vibroimpact module with rodtype oscillatory system. By means of Finite Element Method, there were determined the frequencies of free oscillations, contact and equivalent stresses under different positions of the supports.
The analysis of numerous publications dedicated to oscillatory systems of various vibratory technological machines has shown that the problems of implementation of combined discretecontinuous threemass systems are not thoroughly investigated. Therefore, in the previous publications of the authors, the optimal inertia and stiffness parameters of the threemass discrete vibratory systems were substantiated [9], and the stressstrain state of the doublespan elastic beam was analyzed [10]. The novelty of the present paper consists in determination of natural frequencies and substantiation of the possibilities of implementing the continuous rodshaped members with distributed inertia and rigidity parameters for exciting the oscillations of the threemass discretecontinuous interresonance vibratory systems.
2. Modelling of forced oscillations of continuous member of vibratory system
As a continuous member, there is taken a straight rod the dimensions of crosssection of which are small in comparison with its length. The suggested design diagram of fixation of the rod is like that shown in Fig. 1. The advantage of such a fixation is the fact that there exists a cantilever with considerable inertness (as to its length); such a cantilever can be set in motion mechanically. This can be ensured if one of the supports is disturbed by an eccentric drive [9]. This enables us to a priori assert that the virtually separated mass and rigidity will be accepted by the twomass oscillatory system as a third independent mass, and the threemass discretecontinuous oscillatory system will be an interresonance one.
Fig. 1. Calculation diagram of a continuous rodshaped member (a body with distributed inertia and rigidity parameters) considered as a disturbing body of the threemass vibratory system
The determination of the natural frequency of a rod, as a beam which has three spans, is carried out in terms of the reactions ${R}_{1}$ and ${R}_{2}$ of the supports. As the reference point the left end of the rod is chosen; this left end is free (Fig. 1). The equations of deflections ${w}_{m}\left(x\right)$ and ${w}_{r}\left(x\right)$ in the middle section and in the right section, respectively, will be of the following form:
where ${w}_{l}\left(x\right)$ is the deflection equation of the left section of the beam; $E$ is the modulus of elasticity of elongation for the material of the rod; ${J}_{z}$ is the moment of inertia of the crosssection of the rod about the neutral line of the section; $L$, ${L}_{1}$, ${L}_{2}$ are the lengths of the rod’s sections:
where $\rho $ is the density of the material of the rod; $\omega $ is the forced circular frequency; $h,{b}_{n}$ are the width and the thickness of the elastic rod, respectively.
In general, the deflection of a beam’s section can be described by the equation [10]:
where $S\left(\xi x\right)$, $T\left(\xi x\right)$, $U\left(\xi x\right)$, $V\left(\xi x\right)$ are combinations of circular and hyperbolic functions (KrylovDuncan functions); $A$, $B$, $C$, $D$ are arbitrary constants.
For the middle and the right sections of the beam (see Eqs. (12)), the functions $V\left(\xi x\right)$ takes the following forms:
$V\left(\xi \left({x({L}_{1}+L}_{2})\right)\right)=\frac{1}{2}\left(\mathrm{s}\mathrm{h}\left(\xi \left({x({L}_{1}+L}_{2})\right)\right)\mathrm{s}\mathrm{i}\mathrm{n}\left(\xi \left({x({L}_{1}+L}_{2})\right)\right)\right).$
Since the left end of the rod is free, there take place lateral displacement $w\left(x\right)$ and the angle of flexion $\theta \left(x\right)=dw\left(x\right)/dx$; the bending moment $M\left(x\right)=E{J}_{z}{d}^{2}w\left(x\right)/d{x}^{2}$ and the shear force $Q\left(x\right)=E{J}_{z}{d}^{3}w\left(x\right)/d{x}^{3}$ are equal to zero.
The boundary conditions for the left free end of the beam are:
Thus, taking into account the boundary conditions Eq. (67), the general solution Eq. (4) in the left section $\left(0\le x\le {L}_{1}\right)$ takes the form:
Using Eq. (1), the deflection equation for the middle section $\left({L}_{1}\le x\le {L}_{2}\right)$ assumes the form:
According to Eq. (2), for the right section $\left({L}_{2}\le x\le L\right)$ the deflection equation is the following:
The Eq. (810) are written in general form. To determine the first natural frequency of the rod let us deduce a frequency equation. Let us form a system of four equations, which satisfy their corresponding boundary conditions at the ends of the rod and at the supports. Let us consider the free right end of the rod (Fig. 1). The boundary conditions of its fixation are:
Using Eq. (10) for $x=L$, and taking into account the condition Eq. (12), let us form the first equation of the aforesaid system:
$+\frac{{R}_{2}{\xi}^{2}}{{\xi}^{3}E{J}_{z}}\bullet T\left(\xi \left(L{L}_{1}{L}_{2}\right)\right),$
$T\left(\xi \left(L{L}_{1}\right)\right)=\frac{1}{2}\left(\mathrm{s}\mathrm{h}\left(\xi \left(L{L}_{1}\right)\right)+\mathrm{s}\mathrm{i}\mathrm{n}\left(\xi \left(L{L}_{1}\right)\right)\right),$
$T\left(\xi \left(L{L}_{1}{L}_{2}\right)\right)=\frac{1}{2}\left(\mathrm{s}\mathrm{h}\left(\xi \left(L{L}_{1}{L}_{2}\right)\right)+\mathrm{s}\mathrm{i}\mathrm{n}\left(\xi \left(L{L}_{1}{L}_{2}\right)\right)\right).$
Using Eq. (10) for $x=L$, and taking into account the condition Eq. (13), let us form the second equation of the aforesaid system:
$+\frac{{R}_{2}{\xi}^{3}}{{\xi}^{3}E{J}_{z}}\bullet S\left(\xi \left(L{L}_{1}{L}_{2}\right)\right),$
$S\left(\xi \left(L{L}_{1}{L}_{2}\right)\right)=\frac{1}{2}\left(\mathrm{c}\mathrm{h}\left(\xi \left(L{L}_{1}{L}_{2}\right)\right)+\mathrm{c}\mathrm{o}\mathrm{s}\left(\xi \left(L{L}_{1}{L}_{2}\right)\right)\right).$
Since there are four unknowns – constants $A$, $B$ and reactions of the supports ${R}_{1}$, ${R}_{2}$ – it is necessary to derive two more equations. Considering the vertical displacement ${Y}_{2}$ of the left support ($x={L}_{1}$), where the reaction is equal to ${R}_{1}$ (Fig. 1), Eq. (8) can be written as:
Taking into account the vertical displacement ${Y}_{2}+\epsilon $ of the right support ($x={L}_{1}+{L}_{2}$), where the reaction is equal to ${R}_{2}$ (Fig. 1), Eq. (9) can be written as:
$T\left(\xi \left({L}_{1}+{L}_{2}\right)\right)=\frac{1}{2}\left(\mathrm{s}\mathrm{h}\left(\xi \left({L}_{1}+{L}_{2}\right)\right)+\mathrm{s}\mathrm{i}\mathrm{n}\left(\xi \left({L}_{1}+{L}_{2}\right)\right)\right),$
$V\left(\xi {L}_{2}\right)=\frac{1}{2}\left(\mathrm{s}\mathrm{h}\left(\xi {L}_{2}\right)\mathrm{s}\mathrm{i}\mathrm{n}\left(\xi {L}_{2}\right)\right).$
From Eqs. (14, 16, 18, 20), the following system of four equations of forced oscillations of the rod is obtained:
3. Determination of first natural frequency of continuous member of vibratory system
From the system of Eq. (22), the determinant of the frequency matrix is formed:
A priory, the lengths $L=$ 0.8301 m, ${L}_{1}=$ 0.475 m, and ${L}_{2}=$ 0.355 m are set; then, the determinant of the frequency matrix Eq. (23) is equal to zero at $\xi L=$ 2.7472.
The circular frequency of free oscillations of the rod, according to Eq. (3), is the following:
where $E=$ 2.1∙10^{11} Pa; ${J}_{z}=h\bullet {b}_{n}/12=$ 7.092∙10^{10} m^{4}; $\rho =$ 7850 kg/m^{3}; $h=$ 0.045 m; ${b}_{n}=$5.74∙10^{3} m.
For a solidbody model of the rod, the first natural frequency of oscillations of the rod has been determined by means of Finite Elements Method in SolidWorks software (${v}_{frs}=$ 14.972 Hz) (Fig. 2).
The inaccuracy of the obtained results is scanty; therefore, it can be asserted that the rodshaped continuous member is adequately in compliance with the discrete model (with the mass ${m}_{3}=$ 0.313 kg and the stiffness coefficient ${k}_{23}=$ 2.865∙10^{3} N/m) ensuring the first natural frequency of the threemass vibratory system of approximately 15 Hz [9].
Fig. 2. Results of computer simulation of natural frequencies of the elastic rod in SolidWorks software
4. Conclusions
After considering a discrete model of the vibratory machine, the natural frequency concerning the oscillations of its reactive mass (disturbing body) can be determined; that frequency became an initial parameter which must be ensured by the continuous member. The simplified calculation diagram of the discretecontinuous vibratory system is shown in Fig. 1. In such a system, the elastic rodshaped continuous member is fixed by hinged supports to the oscillating body; this enables us to disturb one of the supports by an eccentrictype vibration exciter.
According to classical theory of oscillations of elastic rods, having used the KrylovDuncan functions, the mathematical model describing forced oscillations of the continuous member considered as a disturbing body of the threemass discretecontinuous vibratory system has been established, and the corresponding frequency equation has been analytically derived. According to this equation, the decisive parameter for synthesis of the interresonance discretecontinuous vibratory system – the value of the first natural frequency in compliance with the natural frequency concerning the oscillations of the reactive mass of the discrete model – is established. This indicates that the rodshaped continuous member is in compliance with the discrete model (with the mass ${m}_{3}=$ 0.313 kg and the stiffness coefficient ${k}_{23}=$ 2.865∙10^{3} N/m) ensuring the first natural frequency of the threemass vibratory system of approximately 15 Hz.
Computer simulation of natural frequencies of a rodshaped continuous body has been carried out using Finite Elements Method in SolidWorks software; the obtained results (Fig. 2) are in satisfactory agreement with the results of analytical calculations (Eq. (24)).
The implemented prototype of a vibratory machine with a continuous member has generalized the results of analytical calculations and those of simulation; it confirmed the adequacy of the suggested theoretical principles.
References
 Rajesh K., Saheb K. M. Large amplitude free vibration analysis of tapered Timoshenko beams using coupled displacement field method. International Journal of Applied Mechanics and Engineering, Vol. 23, Issue 3, 2018, p. 673688. [Publisher]
 Chen LiQun, Peng Li, Zhang AQiang, Ding Hu Transverse vibration of viscoelastic Timoshenko beamcolumns. Journal of Vibration and Control, Vol. 23, Issue 10, 2015, p. 15721584. [Publisher]
 Jaworski J. W., Dowell E. H. Free vibration of a cantilevered beam with multiple steps: comparison of several theoretical methods with experiment. Journal of Sound and Vibration, Vol. 312, Issues 45, 2008, p. 713725. [Publisher]
 Gharaibeh M. A., Obeidat A. M., Obaidat M. H. Numerical investigation of the free vibration of partially clamped rectangular plates. International Journal of Applied Mechanics and Engineering, Vol. 23, Issue 2, 2018, p. 385400. [Publisher]
 Joubaneh Eshagh F., Barry Oumar R., Tanbour Hesham E. Analytical and Experimental Vibration of Sandwich Beams Having Various Boundary Conditions. Shock and Vibration, Vol. 2018, 2018, p. 3682370. [Publisher]
 Shi Xianjie, Shi Dongyan Free and forced vibration analysis of Tshaped plates with general elastic boundary supports. Journal of Low Frequency Noise, Vibration and Active Control, Vol. 37, Issue 2, 2018, p. 355372. [Publisher]
 Buchacz A. The supply of formal notions to synthesis of the vibrating discretecontinuous mechatronic systems. Journal of Achievements in Materials and Manufacturing Engineering, International OCOSCO World Press, Vol. 44, Issue 2, 2011, p. 168178. [Search CrossRef]
 Gursky V., Kuzio I. Dynamic analysis of a rod vibroimpact system with intermediate supports. Acta Mechanica et Automatica, Vol. 12, Issue 2, 2018, p. 127134. [Publisher]
 Korendiy V., Lanets O., Kachur O., Dmyterko P., Kachmar R. Determination of inertiastiffness parameters and motion modelling of threemass vibratory system with crank excitation mechanism. Vibroengineering Procedia, Vol. 36, 2021, p. 712. [Publisher]
 Lanets O. S., Kachur O. Yu., Korendiy V. М. Classical approach to determining the natural frequency of continual subsystem of threemass interresonant vibratory machine. Ukrainian Journal of Mechanical Engineering and Materials Science, Vol. 5, Issues 34, 2019, p. 7787. [Publisher]
Cited By
Vibroengineering PROCEDIA
Oleksandr Kachur, Oleksii Lanets, Vitaliy Korendiy, Vasyl Lozynskyy, Oleh Kotsiumbas, Oleksandr Havrylchenko, Nadiia Maherus

2021
