CC BY
21
DETERMINING THE DYNAMIC CHARACTERISTICS OF HYDRAULIC EXCAVATORS
Republic of Serbia,
e-mail: boro.lazarevic92@gmail.com,
ORCID iD: http://orcid.org/0000-0002-1563-3623
http://dx.doi.org/10.5937/vojtehg66-14400
FIELD: Construction Engineering ARTCLE TYPE: Original Scientific Article ARTICLE LANGUAGE: English
CP CP
>
Olgica Z. Lazarevic3, Boro M. Lazarevicb g
a University of Defense in Belgrade, Military Academy, 8
Military Mechanical Engineering Department, Belgrade, Republic of Serbia,
e-mail: olgica.l.novakovic@gmail.com, ORCID iD: http://orcid.org/0000-0001-5381-6888 b China Communication Construction Company, Belgrade,
T3
Summary:
This paper deals with the research into the dynamic q characteristics of hydraulic excavators i.e. natural frequencies and movements. The way hydraulic excavators work is described by
<u
ro
N <0
the mathematical model of six generalized coordinates based on o the Lagrange equations of the second kind. The obtained mathematical model of the excavator system represents the system of six nonlinear, nonhomogeneous differential equations of the second kind. By simulating real parameter values, we managed to obtain the functions of movement and their natural frequency changes within the allowed degrees of freedom, making it possible to examine the nature and intensity of those movements and the impacts of these parameteres on excavator operators as well as on excavator construction optimisation.
Key words: hydraulic excavator, dynamic characteristics, mathematical model, natural frequencies, movement.
<1> tr>
CD
Introduction
® Excavators are machines used for up to category four soil excavation
and for transporting the excavated soil to a landfill or to some means of transport. They belong to the group of mobile machines. By the way they work, they can be continuous and cyclic machines. Soil excavation can be ai done both on land and under water. The research deals with the hydraulic ¡2 wheel excavator BTH-600. Excavators are primarily used for excavating g the soil with cyclic work. ° The research based on both theory and experiments proved that
< hydraulic excavators, especially when in an open configuration position, behave as oscillation systems while transporting or unloading the soil. That o is why the stability of the balance system position is analysed by the m dynamic actuation of the system. £ In the available reference books (Janosevic, 1997), (Holclajtner et al,
¡i 2010, pp.5-25), (Lazarevic et al, 2008), (Frimpong et al, 2008, pp.20-29), (Vaha & Skibniewski, 1993, pp.148-158), (Lazarevic et al, 2014), the dynamic behaviour of excavators is analysed by maximum five degrees of freedom. However, from the construction itself, together with the way excavators work, it is quite natural for the model space movement to be crj expanded by one more degree of freedom, that is by revolving around the SZ transverse main central axis of the undercarriage inertia, which has been done in this work.
X LU I—
§ Mechanical-dynamic model of excavators
o
> A mechanical model of hydraulic excavators (Figure 1) consists of
undercarriage (1), rotary link (platform) (2) and manipulators (3, 4, 5). The lower part of the excavator, the undercarriage, may be placed on caterpillars, pneumatics or on a floating object. The elements of the subsystem for transferring the power and the elements of the subsystem for transferring the signal are placed on the rotary link (the platform). The plane of digging and the plane of unloading are decided by rotating the rotary link (platform) in relation to the basic undercarriage. The basic and rotary links are connected with a joint in the shape of a radial axial bearing and they form the basis or the excavator machine (Lazarevic, 2003).
Figure 1 - Physical model of a hydraulic BTH-600 excavator: 1 - undercarriage (pneumatics), 2 - rotary link (platform), 3 - boom, 4 - arm (stick), 5 - bucket Рис. 1 - Механическая модель гидравлического экскаватора BTH-600: - пневмоколесное ходовое устройство, 2 - поворотная платформа 3 - стрела, 4 - рукоять, 5 - ковш Слика 1 - Механички модел хидрауличког багера BTH-600: 1 - ослоно-кретни члан (пнеуматици), 2 - обртни члан (платформа), 3 - стрела, 4 - рука (држач), 5 - кашика
When describing the movement of the observed system (Figure1 and Figure 2) the following assumptions were taken into account:
- Dynamic model of the excavator is a non-conservative system with stationary and ideal links,
Vo - Small system oscillations around the stable balance position are
observed,
- Excavator backing has the elastic and damping properties, oc - Parts of the excavator kinematic chain are rigid bodies,
E - Due to oil viscosity and stiffness, hydraulic actuators of the drive
o mechanisms are elastic and damping elements,
" - Modulus of hydraulic oil stiffness is constant and independent from
pressure and temperature.
X
o
LU
H
>-
Q1 <
Generalized coordinates of the mechanical system of the hydraulic excavator
The description of a dynamic excavator model (Figure 2) is defined by the generalized coordinates:
M=k (pdaPy) (1)
0 where:
2 zc -is the vertical displacement of the mass centre of the
<
>o
>
undercarriage,
w p - is the angle of rotation around the longitudinal main central axis
o of inertia of the undercarriage,
o 9 - is the angle of rotation around the transversal main central axis
of inertia of the undercarriage,
a - is the angle of manipulator boom rotation around the axis O3y3 of the joint to which the rotary link is attached; P - is the angle of manipulator arm rotation around the axis O4 y4 of the joint in the end of the boom, and
y - is the angle of manipulator bucket rotation around the axis
O5 y5 of the joint in the end of the arm.
It is assumed that the generalized coordinates of small size are measured from the position of the stable balance of the system. The critical excavator positions are analysed when the manipulator plane forms the right angle with the undercarriage longitudinal plane (Figure 3).
Figure 2 - Dynamic model of the hydraulic excavator. The platform and the kinematic chain have been rotated by у = 900 Рис. 2 - Динамическая модель гидравлического экскаватора. Платформа и кинематическая цепь экскаватора под углом поворота у = 900 Слика 2 - Динамички модел хидрауличког багера. Платформа и кинематички ланац багера заокренути су за угао у = 900
CP CP
<л о
го
>
го о
X ф
о "5
"U
ф
о
го .с о о Е
тз
Е ей <и Q
"го
ц—' ф
о
■О >
си го
N ГО
Coordinate systems
The position of the excavator is defined in relation to the absolute coordinate system Oxyz. The coordinate beginning O of the absolute system is in the undercarriage mass centre for the whole system (undercarriage, rotary link, boom, arm and bucket). The Ox axis is directed
ф (Л
CD CD
"о
>
0
01
Q¿ Ш
O O
_l
<
o 2:
X
o
Ш
I—
>-
Q1 <
(Л <
-J
O >0
x
Ш I—
O z
—3
O >
towards the longitudinal axis, while the Oy axis is directed towards the main transversal central axis of inertia in the position of static balance of the whole system, and the Oz axis is directed vertically downwards. To determine the positions of the excavator kinematic chain parts, local coordinate systems have been set up (Figure 2).
Figure 3 - Accepted generalized mechanical system coordinates Рис. 3 - Принятые обобщенные координаты механической системы Слика 3 - Усво^ене генералисане координате механичког система
1. O1 x1 y1 z1 - mobile coordinate system accepted for the undercarriage mass center Oj. The Oxj axis is parallel with the longitudinal axis, the Oy1 axis is parallel with the transverse main central axis of the undercarriage inertia in the position of balance and the Ozj axis is directed vertically downwards.
2. O2x2y2z2 - mobile coordinate system connected with the center of the radial axial bearing O2.
3. O3x3y3z3 - coordinate system connected with the joint center O3 and redirected by the generalized coordinate a in relation to the coordinate system O2x2y2z2.
4. O4x4y4z4 coordinate system connected with the joint center O4 and redirected by the generalized coordinate p in relation to the coordinate system O3x3y3z3.
5. O5x5y5z5 coordinate system connected with the joint center O5 and redirected by the generalized coordinate y in relation to the coordinate system O4x4y4z4.
Transformation matrices
cp cp
ro
>
ro о X ф
To determine the speeds of the characteristic points with such complex mechanical systems, the method applied is the one based on setting up and transforming the adequate mobile coordinate systems.
To determine vectors of the characteristic points' position in relation to the absolute coordinate system, we determined the matrices of elementary rotations around the axes Ox, Oy and Oz and the transformation matrices of the immobile coordinate system Oxyz into the mobile coordinate systems Oixiyizi where /'=1-5 as well as the matrices of mobile coordinate
system transitions from one into another.
For calculating kinetic energy, it is necessary to calculate the speed of the mass centres of some hydraulic excavator parts. To be able to do this, we need the system characteristic points vector positions (Figure 4). The absolute coordinate system, the undercarriage coordinate system and the undercarriage mass center are placed together in one point, i.e. O = Oj.
T3
ф о
ro .с о о Е
тз
Е ф
ф а
"то ф
О
■о >
CD ГО N ГО
Figure 4 - Radius vectors of the characteristic system points Рис. 4 - Радиус-векторы характеристических точек системы Слика 4 - Радиус вектори карактеристичних тачака система
<1> tr>
CO
Z)
X
>
Mathematical excavator model
co~ Mathematical excavator model is based on general mechanical
theorems or on applying Lagrange equations of the second kind. These are the starting assumptions taken into account when creating the mathematical model:
1. The backing and the parts of the excavator kinematic chain are Si modelled with rigid bodies,
^ 2. During its manipulating tasks, the excavator work is stable with no
o possibility for the movements in the zero joint to appear,
3. Excavator kinematic chain is observed as an open configuration o chain with its last part - the bucket - being affected by the
technological resistance of digging, and
¡^ 4. During its manipulating tasks, the excavator kinematic chain is
affected by the forces of gravity, inertia and some external forces -& digging resistance.
The excavator model space is determined by the absolute coordinate system Oxyz . The excavator backing lies in the horizontal plane Oxy of the absolute coordinate system, while the vertical axis Oz is parallel with the kinematic pair seen as the undercarriage and the rotary part. Each d excavator part is determined in its local coordinate system Oixiyizi
o (Figures 1 and 2). The coordinate beginning of the local coordinate system is parallel with the center of the joint connecting that particular TE part with the previous one, while some axis are directed in such a way o that they are parallel with the joint axis, and as for manipulators (boom, g arm and bucket), they are parallel with the kinematic longitude of the part. For the bucket, the local coordinate system axis O5x5 goes through
the joint center O5 and through the center of the bucket cutting edge.
Differential movement equations are based on Lagrange equations of the second kind:
d fdEk} dEh dE
dt
k = QN, qr = zc,V,G,a,P,Y,r = 1,2,3,4,5,6 (2)
dqr) dqr dq, where:
Ek - is system kinetic energy, Ep - is system potential energy, QN - is generalized forces.
Resolving the systems of differential equations
The presented mathematical model represents the system of six nonlinear, nonhomogeneous differential equations of the second kind that was resolved by using discrete numerical methods and packets in the MATLAB program. The Runge-Kutta Method applies perfectly well to these systems of differential equations so it is necessary to write the equation in the matrix form:
A ■ q + C ■ q + B ■ q = D (3)
where the correspnding matrices A, B, C and D are the functions of the generalized coordinates and speeds. The matrix A is the mass (inertia) matrix,
aii ai2 ai3 a14 a15 ai6
a2l a22 a23 a 24 a25 a26
A = a3i a32 a33 a34 a35 a36
a4l a42 a43 a44 a45 a46
a5l a52 a53 a54 a55 a56
a6l a62 a63 a64 a65 a66 _
The elements of the mass (inertia) matrix are
condition has been fulfilled
a.. = a
J'
The matrix B is a rigidity matrix.
~bu b12 b13 0 0 0 "
b21 b22 b23 0 0 0
B = b3i b32 b33 0 0 0
0 0 0 b44 0 0
0 0 0 0 b55 0
0 0 0 0 0 b66 _
(4)
(5)
(6)
up cp
O
ro
>
ro o X <u
T3
<u o
sh
ro
-C
o o
e
T3
E
CD CD
Q
m
CD
o
■O >
CD
ro
N <0
The matrix D is an actuation forces matrix.
<D tí)
CD CD
"5
>
0
01
Q¿ UJ
ZD O o
_l
<
o 2:
X
o
LU
I—
>-
Q1 <
D =
dr
d2
d
31
d4
d
51
(7)
d6
The matrix C is a damping matrix and it is hard to form it because
zc,9,à,j3 and f are close-coupled and they are also close-coupled
with zc,ç,6,a, P and f making them almost impossible to split up. It is
possible to make the systems of differential equations linear, and being taken as given, it would be possible to find out the matrix C. However, this does not have to be done, because resolving nonlinear equations is provided by MATLAB. Because of applying the MATLAB programme package it is necessary to write the equation in the matrix form:
CO <
-j
(!) >Q
x
LU I—
O z
—3
O >
where it comes that:
A =
a,
a,,
a
a.
a
a
A • q = -(C • q + B • q - D)
A • q = - S
a
12
a
31
a
a
51
a
a
a
a
32
a
a
52
a
a
a
a
33
53
ai4 ai5
a
16
a
a34 a35
a44 a45
a54 a55
a64 a65
a
36
a
a
56
a
zC sii
(p s21
d s3i
á s41
P s51
Y s61
(8)
(9)
24
From the previous it comes that:
q = - A'1 • S
Or
ZC 5ii
p 521
9 = A-1 - 531
a 541
P 551
_y_ _561 _
(10)
where the matrix S is the column
- S =
21
41
It is possible to resolve accelerations and speeds prepared in this way with the MATLAB programme package with variable changes done beforehand as following:
Xi Zç, ^2 Z
C X4 = p, x5 =9, x6 =9
x7 = a, x8 = a, x9 = p, x10 = p, xu = y, x12 = y The form of the system of differential equations looks like this: zC
(11)
A(zc ,p, 9, a, P, y)-
S (zc, p, 9, a, P, y, Zc, P, 9, a, P, y)
(12)
'C P
9
a P
y
and is transformed into the system of twelve differential equations of the first kind making it suitable for numerical resolving by applying the MATLAB programme package.
p p
o
ro
>
ro o X CD O
"5
sh
T3
CD o
ro
-C
o o
'E ro
c ^
T3 CD -C
D) C
'c
'e
CD CD
Q
CD
o
■O >
CD
ro
N
ro
5
5
31
5
51
5
61
<D tr>
CD CD
>
Determining the dynamic behaviour of hydraulic excavators
o Resolving equation systems (3) used for describing the dynamic
behaviour of the hydraulic excavator can be divided into: • determinig natural frequencies and oscillation types □cT • determining the vectors of mobility qi (qi = zc Q, a,
o Determining natural values
<
Vibration issues are mostly connected with resonance, where one or more forms of vibration are activated by work power forces. Vibration forms lying among dynamic work forces' frequency levels are always a £ potential threat to cause problems. The modal analysis makes it possible < to determine frequencies that may cause resonance to appear. The modal analysis is the process of determining all modal parametres sufficient for formulating a mathematical-dynamic model. The modal parametres are: modal frequencies, modal dampings and modal shapes. One of the very w important mode properties is that any forced or free structure response
CD
>Q
may be brought together in one set of modes (Batinic, 2013), (Den Hartog, 1972).
Natural system values, i.e. natural frequencies and the corresponding EH vectors - modal forms - will be obtained through the free vibration system analysis while considering the homogeneous part of equation (3):
g A ■ q + B ■ q = 0 (13)
For the purpose of mathematical modelling of the free system vibrations, the aim of which is to determine their own values, it could be assumed that damping C within the system is negligible. According to the well-known procedure for determining the equation system own values (3), the dynamic matrix U is formed together with its characteristic determinant, i.e. a characteristic system equation:
det(A ■ I - U)= 0
where:
- U = A-1 ■ B - the dynamic system matrix where A and B are the inertia matrix and the system rigidity matrix,
- unit matrix,
л =
а
- i.e. а is the system natural frequency.
The observed hydraulic excavator is described by means of the system of nonlinear differential equations meaning that the system dynamic properties cannot be taken as constants. It applies in the same way to both rigidity and damping coefficients and their matrices coefficients. In such cases, it is not possible to make any separations based on their characteristic forms, by means of the vertical coordinates into a set of independent equations. A possibility of a step-by-step numerical integration is used for this kind of analysis - for certain time interval, the system is cosidered as linear with the characterisitcs it had at the beginning of that time interval. So, a nonlinear analysis is approximated by the series of variable gradual system analyses (Batinic, 2000).
The matrix U own values are zeros of the matrix U characteristic polynomial i.e. the roots of the characteristic equation (13). In order to obtain own values and their corresponding vectors, we considered the system parametres at a certain point of time, taken from the results of the numerical system analysis of differential equations at a certain point of time, for example at a point of time t = 0.100 s, i.e. in iteration n = 19999
at a small initiative system motion. This is a necessity because the mass matrices and the system rigidities are subordinate to the generalized coordinates and their speeds, and as such they cannot be used for obtaining dynamic system matrices. The matrices of rational values were obtained by this procedure (Jovanovic et al, 2004). The system state vector is shown in Table 1.
Table 1 - System State Vector Таблица 1 - Вектор состояния системы Табела 1 - Вектор стаъа система
zc (P P в в
-0.0022 0.0913 -0.0010 -0.9068 -0.0013 0.0105
Ol
со
CP CP
<л о
го
>
го о
X ф
тз
ф
о го
-С
о о Е
тз
Е ф
ф а
"то ф
о
■о ■>
ф го
N
го
a a ß ß У у
0.0049 -0.1762 -0.0903 -3.2058 0.0398 8.7507
The inertia matrix A and the rigidity matrix B of the system at the determined point of time are:
<u tn
CO CD
"o
>
o CM
QC UJ QC
O
o <
o 2:
X
o
LU
I—
>-
CC <
< -J
CD >Q
x
LU I—
O z
—3
O >
1.0633 0.9423 - 0.2444 1.5645 - 0.2090 0.0080 "
0.9423 6.6536 - 0.3062 5.8923 -1.8336 0.8307
_ 0.2444 - 0.3062 4.1787 0.4836 0.0661 - 0.0009
.0e + 0.04 ■
1.5645 - 5.8923 0.4836 9.1443 1.3128 - 0.3437
- 0.2090 _ 1.8336 0.0661 1.3128 0.8656 - 0.3928
- 0.0080 0.8307 - 0.0009 0.3437 - 0.3928 0.4432
" 0.600 0 - 0.6120 0 0 0 "
0 1.6869 -0.9890 0 0 0
- 0.6120 - 0.9890 2.4970 0 0 0
1.0e + 0.08 ■
0 0 0 0.1169 0 0
0 0 0 0 0.0056 0
0 0 0 0 0 0.0038
Applying the function [V, D\ = eig (U) from the MATLAB programme package will make it possible to obtain the diagonal matrix D containing own values on the main diagonal. We also obtained the full matrix V whose columns are the corresponding matrix U vectors.
D = 1.0e + 0.04 ■
1.4791 0 0 0 0 0
0
0.8468 0 0 0 0
0 0
0.2449 0 0 0
0 0 0
0.0281 0 0
0 0 0 0
0.0130 0
0 0 0 0 0
0.0041
V =
- 0.1588 -0.8613 - 0.6478 - 0.0059 - 0.0169 0.0071
- 0.5614 0.2806 - 0.3303 - 0.0034 - 0.0115 0.0088
0.1679 0.1327 - 0.4723 - 0.0027 - 0.0082 0.0051
- 0.2906 - 0.0266 - 0.2479 0.1243 0.2179 - 0.0655
-0.7119 0.2957 - 0.4279 - 0.6816 0.1570 -0.7386
0.2004 - 0.2715 0.0603 - 0.7210 0.9630 0.6708
The range of natural frequencies was obtained by finding the square root of the diagonal matrix D solution and dividing it by 2n (Table 2).
Table 2 - Natural frequencies Таблица 2 - Естественные частоты Табела 2 - Сопствене фреквенцще
Natural frequenci es ®Zc
rad/s 121.6164 92.0227 49.4916 16.7553 11.4087 6.4191
Hz 19.3558 14.6459 7.8768 2.6667 1.8158 1.0216
<л о
го
>
го о
X ф
тз
The natural frequencies for the state vector of equal system activation with the time of t = 0.025s and the iteration of 4999 were calculated as shown in Table 3.
Table 3 - Natural frequencies Таблица 3 - Естественные частоты Табела 3 - Сопствене фреквенци^е
Natural frequencies ®Zc ®a ФУ
rad/s 0+131.90/ 96.16 78.41 47.39 14.12 8.62
Hz 0+20.9918/ 15.3041 12.4792 7.5420 2.2474 1.3718
ф о
го
-С
о о Е
тз
Е
ей ф
а
"то ф
О
■CJ ■>
си
N ГО
Determining motion vectors on the generalized coordinates
Differential equations were resolved by applying the MATLAB packet programme, and the generalized coordinates of free oscillation models and their derivatives for the time t = 0.1s were obtained as well. The simulation of the observed hydraulic excavator behaviour was done at the beginning of moving the soil away, after the excavation of the soil had been finished with the angles of a0 =-35°, =-80° and y0 = 1600
together with the angle y/ = 900 regarded as the most unfavourable for excavator operations. The initial movements were very near zero. The simulation results are shown in Figures 5 - 10.
<u
(Л
CO CD
"o
>
О CM
DC Ш DC ZD О
о <
с;
Z X О ш
i>-
DC <
(Л <
-J
О >о
X ш н
о
О >
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1
<И
Figure 5 - Vertical motion of the
undercarriage mass center - zc
Рис. 5 - Вертикальное перемещение центра массы опорно-ходового
устройства - zc
Слика 5 - Вертикално помераше средишта масе ослоно-кретног
члана - zc
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.
Figure 6 - Angle of the undercarriage rotation around the longitudinal main central
axis of inertia oi x1 - p Рис. 6 - Угол поворота опорно-ходового устройства вокруг продольной главной центральной оси инерции о1 Xl -p
Слика 6 - Угао обрташа ослоно-кретног члана око уздужне главне централне осе
инерци]е O1 x1 -p
Figure 7-Angle of the undercarriage rotation around the transversal main cental axis of inertia O1 y1 -в
Рис. 7 - Угол поворота опорно-ходового устройства вокруг поперечной главной центральной оси инерции о1 y1 - в Слика 7 - Угао обрташа ослоно-кретног члана око попречне главне централне осе инерцц'е O1 y1 - в
Figure 8 - Angle of the boom rotation of the manipulator around the axis O3y3 - a
Рис. 8 - Угол поворота срелы манипулятора вокруг оси O3y3 - a
Слика 8 - Угао обрташа стреле манипулатора око осе озy3 - a
Figure 9-Angle of the arm rotation of the manipulator around the axis O4y4 of the
joint in the end of the boom - p
Рис. 9 - Угол вращения рукояти манипулятора вокруг оси O4y4
шарнира на конце стрелы - p
Слика 9 - Угао обртаъа руке манипулатора око осе O4 y4 зглоба
на кра]у стреле - p
Results analysis
¡н
Figure 10-Angle of the bucket rotation of the manipulator around the axis O5y5 of the
joint in the end of the arm - у
Рис. 10 - Угол вращения ковша манипулятора вокруг оси O5y5 шарнира
на конце рукояти Слика 10 - Угао обртаъа кашике манипулатора око осе O5y5 зглоба на
кра]у руке - у
The observed hydraulic excavator system was activated by a very small initial value of one of the generalized coordinates. This possibility is allowed by the simulation programme package. That way, the system is exposed to free oscillations. This state is described by homogeneous nonlinear differential equations.
The natural frequencies obtained by the mathematical model based on the generalized coordinate z have the values of o>Zc = 19.3558#z and
coZc = 0 + 20.9918/. These frequencies on the z axis have the biggest
impact on the excavator operator. When compared with some resonant frequencies [3-20Hz] it is possible to claim that, for some human organs, they partially overlap. This could cause problems to the operator in case of being in action over a longer period of time.
By analyzing the previously described numerical experiments, it is clear that some kinds of movement freedom, i.e. generalized coordinates, generate free oscillations of certain frequencies that do not allow the
CP CP
ГО
>
ro о X Ф
T3
Ф о
<0 -С
о о 'Е
тз
Е
ей <и
а
"го
<и
о
■о ■>
ф го
N
го
system to easily calm down after the initiative activation time. It was even proved that, with only one activation, within the observed time interval, at one iterration only, the system has free oscillations of one degree of freedom, and that at the second iterration it has free socillations of another > degree of freedom. All these facts clearly indicate that the system is dynamically unstable or at least on the very edge of the stability reference values. The described fenomenon of free oscillation is proved on the oc generalized coordinate zc.
E Figure 5 shows the nature of the vertical movement of the
o undercarriage mass center. It is clear that the oscillation amplitudes are " within the tolerable limits and that they gradually go down with the system tending to calm down.
Figure 6 shows the nature of the rotation angle change q> around the o undercarriage longitudinal axis. The oscillation amplitude is tolerable, but i- there is no evidence of its going up or down within the observed interval of & simulation.
^ Figure 7 illustrates the rotation angle change Q around the transversal
undercarriage axis O1 y1. The oscillation amplitude is within tolerable limits, but like in the previous case, there is no evidence of the amplitude tending to go up or down within the observed interval of simulation. fj The change of the rotation angle of the manipulator boom a around
u the axis O3y3 - a inside the joint connecting the manipulator and the
undercarriage is described in Figure 8. The oscillating nature of this angle's change within the observed time interval is obvious, but with the amplitude tending to go down.
Figure 9 shows the change in the rotation angle of the manipulator arm fi around the axis O4y4 of the joint in the end of the boom. The displayed change is not oscillatory, in the observed time, and it has a tendency to grow up to unacceptable values. This tendency is conditioned by free oscillation of the excavator in the generalized coordinates z.
Figure 10 shows the change in the bucket rotation angle of the manipulator y around the axis O5y5 of the joint in the end of the arm. The
displayed change has the oscillatory character with an increased amplitude value. In the observed simulation time, the tendency to a reduced or increased amplitude cannot be estimated.
The graphics of the generalized coordinates in Figures 9 and 10 do not show the system calm down within the observed time interval because the system is within the unstability area or on the very edge of the stability reference values, which was caused by free system oscillation on one of
CO
the generalized coordinates. Because of this, immediately after the time of =0.1 s the system gets into resonance with a significantly increased oscillation amplitude. In order to get the system calm down, i.e. to make the oscillation amplitudes on all generalized coordinates reach lower values, preferably zero, all real parts of the characteristic system equation square roots ought to be negative. This precondition for the observed excavator has not been met.
Conclusion
This paper deals with the analysis of the dynamic behaviour of hydraulic excavators. By applying the dynamic-mathematical model, we determined the level and nature of dynamic motions and natural frequencies. An excavator operator is a person who is seated in the rotating platform. System vibrations reach him via the excavator backing. This predominantly refers to the vibrations on the generalized coordinate zc, which is a longitudinal direction, i.e. the direction of the operator's backbone and that is why exposure to these vibrations was addressed to with special attention. Besides these vibrations, the operator is exposed to the vibrations on the generalized coordinates q> and Q which additionally makes the situation more complicated and increases possible harmful impacts on the operator's health. The hydraulic excavator natural frequency values, derived from analytical calculations, prove the intensity so high that it might cause serious problems to the operator's health. By comparing the natural frequency intensity values on the generalized coordinates Q, q> and Q with the resonance frequencies for some human organs, the possibility of serious harmful effects is obvious.
Describing the hydraulic excavator dynamic behaviour by the system of nonlinear, nonhomogenous differential equations made it possible to improve the analytical-mathematical model by which a simulated situation is almost like a real life situation of an excavator in operation.
By applying the modal analysis, it was possible to determine certain system frequencies and their corresponding vectors which might cause resonance to appear and thus spoil the operator's working condition comfort.
The numerical analysis was carried out for the BTH-600 excavator structure. The obtained results indicate that the structure was not composed at its best and that a dynamic analysis should be used to improve the dynamic characteristics of the observed system.
<u
(Л
CO
О
>
References
Batinic, V. 2013. Planetary gear dynamic response to mesh parametric excitation. Vojnotehnicki glasnik / Military Technical Courier, 61(1), pp.58-68 (in ° Serbian). Available at: http://dx.doi.org/10.5937/vojtehg61 -2006.
Batinic, V. 2000. Unutrasnja dinamika opterecenja planetarnih prenosnika, o MSc thesis, University of Belgrade: Faculty of Mechanical Engineering (in
0C
w Den Hartog, J.P., 1972. Vibracije u masinstvu. Beograd: Gradevinska
is knjiga (in Serbian).
§ Frimpong, S., Hu, Y. & Inyang, H. 2008. Dynamic Modeling of Hydraulic
< Shovel Excavators for Geomaterials. International Journal of
0 Geomechanics, 8(1), pp.20-29. Available at: http://doi.org/10.1061/(ASCE)1532-
1 3641(2008)8:1(20). ^ Holclajtner, M., Jeremic, O. & Milinovic, M. 2010. Simulacija rafala
visecevnog bacaca raketa. Vojnotehnicki glasnik/Military Technical oc Courier, 58(2), pp.5-25 (in Serbian). Available at: ■ http://doi.org/10.5937/vojtehg1002005H.
Janosevic, D. 1997. Optimalna sinteza pogonskih mehanizama hidraulickih bagera, Ph.D. thesis, University of Nis: Faculty of Mechanical Engineering (in Serbian).
Jovanovic, M., Janosevic, D. & Milic, P. 2004. Modalna fem analiza kostura (5 kretnog mehanizma hidraulickog bagera. In: Naucno-strucni skup Istrazivanje i
2 razvoj masinskih elemenata i sistema, IRMES, Kragujevac (in Serbian). Lazarevic, D. 2003. Inzinjerijske masine 1. Belgrade: Vojna akademija
x Generalstab Vojske Srbije i Crne Gore - Uprava za skolstvo i obuku (in Serbian).
Lazarevic, O., Batinic, V. & Randelovic, A. 2014. Odredivanje brzina
karakteristicnih tacaka hidraulickog bagera. In: 4. medunarodna konferencija za
о geometriju i grafiku moNGeometrija, Vlasina (in Serbian).
Lazarevic, O., Krsmanovic, M. & Peric, S. 2008. Analysis of dynamic behavior generators of hydraulic excavator. In: Kongres MVM, Kragujevac.
Vaha, P.K. & Skibniewski, M.J. 1993. Dynamic Model of Excavator. Journal of Aerospace Engineering, 6(2), pp.148-158. Available at: https://doi.org/10.1061/(ASCE)0893-1321(1993)6:2(148).
ОПРЕДЕЛЕНИЕ ДИНАМИЧЕСКИХ ПАРАМЕТРОВ ГИДРАВЛИЧЕСКИХ ЭКСКАВАТОРОВ
Олгица Ж. Лазаревич3, Боро М. Лазаревич6 а Университет обороны в г. Белград, Кафедра военно
машиностроительного инжиниринга, г. Белград, Республика Сербия 6 Китайская компания коммуникаций и строительства, г. Белград, Республика Сербия
ОБЛАСТЬ: машиностроение со ВИД СТАТЬИ: оригинальная научная статья
ЯЗЫК СТАТЬИ: английский а.
Резюме: _
го
В данной работе представлены исследования динамических д параметров гидравлического экскаватора: естественные ф частоты и движение. Работа экскаватора описана путем математических моделей с шестью обобщенными координатами, полученными на основании уравнений Лагранжа второго рода. Разработанная математическая модель системы экскаватора представляет собой систему нелинейных неоднородных дифференциальных уравнений второго рода. С помощью моделирования с реальными значениями параметров получены функции изменений движения и естественной частоты при допустимых степенях свободы, благодаря чему можно провести о анализ характера и интенсивности движения, воздействия этих | параметров на машиниста и оптимизировать конструкцию экскаватора (изменение массы, жесткости и демпфирования).
Ключевые слова: гидравлический экскаватор, динамические характеристики, математическая модель, естественные частоты, движение.
б
ОБЛАСТ: машинство
ВРСТА ЧЛАНКА: оригинални научни чланак JЕЗИК ЧЛАНКА: енглески
тз
ф '
о го
тз
Е ф
'
ф а
ОДРЕЪИВА^Е ДИНАМИЧКИХ КАРАКТЕРИСТИКА й
ХИДРАУЛИЧКИХ БАГЕРА ф
Олгица Ж. Лазаревич8, Боро М. Лазаревич6 а Универзитет одбране у Београду, Во]на академи]а, Катедра во]номашинског инженерства, Београд, Република Срби]а
го
СМпа Соттипюа^оп Соп81гисйоп Сотрапу, Београд, Република Срби]а ^
Сажетак:
У раду су приказана истраживаъа динамичких карактеристика хидрауличких багера: сопствене фреквенци]е и померак>а. Рад багера описан }е математичким моделом са шест генералисаних координата ко\и }е доби/'ен применом Лагранжевих ]едначина друге врсте. Доби}ени математички модел система багера представка систем од шест нелинеарних нехомогених диференцц'алних ]едначина другог
реда. Симулацирм са реалним вредностима параметара доби}ене су функци]е промене помераъа и сопствене фреквенци]е по дозвоъеним степенима слободе, што § омогуЬава испитиваъе карактера и интензитета тих
помераъа, утица]а параметара на руковаоца багера и оптимизаци}у конструкци}е багера (промене маса, крутости и пригушеъа).
Къучне речи: хидраулички багер, динамичке карактеристике, математички модел, сопствене фреквенци]е, помераъе.
Е
q Paper received on / Дата получения работы / Датум приема чланка: 30.06.2017. о Manuscript corrections submitted on / Дата получения исправленной версии работы / ^ Датум достав^а^а исправки рукописа: 26.11.2017. о Paper accepted for publishing on / Дата окончательного согласования работы / Датум х коначног прихвата^а чланка за об]ав^ива^е: 28.11.2017. о
¡^ © 2018 The Authors. Published by Vojnotehnicki glasnik / Military Technical Courier
(www.vtg.mod.gov.rs, втг.мо.упр.срб). This article is an open access article distributed under the CC terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/rs/)
<
© 2018 Авторы. Опубликовано в «Военно-технический вестник / Vojnotehnicki glasnik / Military Technical Courier» (www.vtg.mod.gov.rs, втг.мо.упр.срб). Данная статья в открытом доступе и распространяется в соответствии с лицензией «Creative Commons»
ад (http://creativecommons.org/licenses/by/3.0/rs/). <
О © 2018 Аутори. Обjавио Воjнотехнички гласник / Vojnotehnicki glasnik / Military Technical Courier
g; (www.vtg.mod.gov.rs, втг.мо.упр.срб). Ово jе чланак отвореног приступа и дистрибуира се у
>0 складу са Creative Commons licencom (http://creativecommons.org/licenses/by/3.0/rs/).
z x
Ш I—
О z
О >
