<1>
LO CD
O >
CONTRIBUTION TO RESEARCH OF THERMAL LOADS IN COMMERCIAL VEHICLES POWERTRAIN MOUNTS
Miroslav D. Demie3, Jasna D. Glisovicb
University of Kragujevac, Faculty of Engineering, cc Department for motor vehicles and motors,
Kragujevac, Republic of Serbia, 3 a e-mail: demic@kg.ac.rs,
0 ORCID iD: http://orcid.org/0000-0003-2168-1370, ^ b e-mail: jaca@kg.ac.rs,
| ORCID iD: http://orcid.org/0000-0002-8251-1722
¡^ http://dx.doi.org/10.5937/vojtehg65-14028
i-
FIELD: Mechanical Engineering, Motor Vehicles and Motors ¡i ARTICLE TYPE: Original Scientific Paper
ARTICLE LANGUAGE: English
Abstract:
Due to the effects of excitations of road micro-roughness and the ¡3 powertrain as well as due to the conditions of movement (acceleration,
braking, curvilinear movement), the weights of commercial vehicles ^ perform complex spatial oscillations. The negative impact of these
1 movements or dynamic loads caused by them and transmitted to the supporting system of the vehicle can be reduced by a proper selection of
NO the position and characteristics of the powertrain mounts. Elastic and
g damping forces in mounts perform mechanical work that turns into heat
V energy. Theoretical, experimental and combined methods can be used for
the analysis of thermal loads of mounts. In this paper, research was carried out with theoretical methods using a mechanical model of the powertrain. Vibrations were observed on a freight vehicle powertrain of FAP 1314 middle-class vehicles. Bearing in mind the presence of classical (rubbermetal) and hydraulic mounts in modern commercial vehicles, the authors analysed the thermal loads of the mentioned vehicle type using a mathematical model.
Key words: excitation, commercial vehicles, mounts, modeling, powertrain, heat.
ACKNOWLEDGMENT: The authors are grateful for the financial support of TR35041 Project named The research of vehicle safety as part of a cybernetic system: Driver-Vehicle-Environment funded by the Ministry of Science and Technological Development of the Republic of Serbia.
CO
o
CD
Introduction s
od !± CP
The powertrain vibrations are transmitted over the mounts to the vehicle frame, while noise is transferred to the vehicle frame (structural) through the environment (air) as well. In order to reduce this effect to a sa- e tisfactory measure, the powertrain is elastically bound to the vehicle carrier § system. At the same time, this system absorbs vibration excitations which -iE are transferred to the powertrain from the road surface through the vehicle suspension system. Due to the kinetic energy of the powertrain resulting from vibrations, mechanical work is converted to heat in the mounts Jü (Mitschke & Wallentowitz, 2004).
In practice, the characteristics of the powertrain system are chosen from the conditions of minimizing noise and vibration. The task can be solved using analytical (theoretical) methods, using mathematical models,
o
o
<u >
ro 'o <D
or experimentally. The progress of computer technology has contributed to e greater application of mathematical models, while experimental methods, 8
due to high costs, are rarely used. In doing so, the thermal loads which powertrain mounts are exposed to are usually ignored, which, in cases of poorly solved heat removal, can cause distortion of mount characteristics. Therefore, it is considered appropriate to analyze their thermal loads. j=
Research has shown that mechanical work due to powertrain vibrations turns into heat energy that is transmitted, in a wider sense, to the environment (Bojic, 2011), (Ilic et al, 1996), (Moran, 2010), (Fermi, 2011), (Demic & Diligenski, 2016), (Mitschke & Wallentowitz, 2004): <g
A = Qt + Qf + Qc, (1)
where A is a mechanical operation (equivalent to the amount of heat), J; Qt o part of the heat surrendering to the mount body, J; Qf part of the heat surrendering to the fluid (for hydraulic mounts), J; and Qo part of the heat surrendering to the environment, J.
The work of elastic and damping forces in mounts is important for the analysis of heat loads, and it can be measured and experimentally measured (Mitschke & Wallentowitz, 2004) with more difficulty. Due to the distribution of heat inside the mount, as well as around it, the problem becomes significantly more complex.
It is noted that most of the heat is transferred by convection (the Q phenomenon that occurs when the fluid convection material exchanges heat within itself). A slightly smaller part of the heat is transferred by conduction (molecules of solids that vibrate at a higher temperature vibrate more rapidly over the equilibrium position and in collision with molecules having a lower temperature cause an energy impulse that accelerates the vibrations of these molecules and increases the temperature) and radiation
.Q
<D
LO CD
"o >
O CM
of
UJ
a.
Z) O O
_l <
o
X
o
LU
H ^
a. <
H
<
CD >o
X LU H O
O >
(occurs between any two objects at different temperatures, whether vacuum is between them or not) (Bojic, 2011), (Mitschke & Wallentowitz, 2004).
Converting mechanical work into heat largely depends on the construction of mounts. Therefore, during the design of the powertrain system, consideration should be given to creating the conditions for the good transfer of heat to the environment. It should be noted that the aim of this paper is not to investigate the transfer of heat to the environment from powertrain mounts, but to determine the amount of heat generated by vibrations in mounts.
As it is known, classic mounts consist of rubber-metal elements, with a dampening of the hysteresis of the tire, while hydraulic mounts have additional damping due to oil which they contain. Bearing in mind the presence of classical (rubber-metal) and hydraulic mounts in modern freight vehicles, the authors analysed the thermal loads of a FAP 1314 vehicle.
Motor vehicles represent complex oscillation systems with a large number of masses of interconnected elasto-damping elements, of most often non-linear characteristics. When traveling on uneven roads, the exitations from road roughness are transferred through the wheels and the elastic suspension system to the frame and further on to the powertrain (through the powertrain suspension system). In addition, the drive unit additionally oscillates the movable mass of the vehicle. The negative impact of vibrations on the characteristics of motor vehicles can be reduced by a correct selection of the position and characteristics of the elastodamping elements. This can be done using theoretical-experimental, theoretical or experimental methods (Frolov & Furman, 1990).
In the literature (Demic, 1990), (Demic, 1997a), (Demic, 1997b), (Demic, 2002), (Demic, 2013), (Simic, 1988), (Simic & Demic, 1992), (Ben-dat, 1998), (Igami et al, 2008), (Genta, 2003), (Shanguan et al, 2016) there are a number of mathematical-experimental methods for modeling the behavior of the powertrain. In this paper, a freight vehicle from the FAP 1314 production program, (whose scheme of the suspension of the powertrain is shown in Figure 1) was observed and the modeling process was applied based on the principles of classical mechanics (Atkins, 2010), (Pars, 1981).
The aim of the research is to begin theoretical analyzes of the problem of thermal loads of powertrain mounts, and the verification of the obtained results will be performed at a later stage in experimental research. It is expected that the results obtained will be used in designing new vehicles from the FAP production program.
Model of a powertrain
Depending on the task being worked out, various mechanical models of powertrains and vehicles can be found. From (Demic, 1990), (Demic,
1997), (Rotenberg, 1972), (Genta, 2003), (Gillespie, 1992), (Milliken & Milliken, 1995), (Simic & Demic, 1992), (Kennings et al, 2014), (Shanguanet al, 2016), it is known that during the analysis of the problem of transferring dynamic loads from the powertrain to the frame (chassis) of the vehicle, the oscillations of the cabin and the carriage compartment can be ignored. More precisely, the analysis includes only the oscillatory movements of the powertrain of the vehicle and the corresponding exitations from the vehicle frame.
The structure of the vehicle model and the powertrain should therefore be chosen to enable the analysis of the desired parameters. This indicates the need to use as simple mathematical models as possible, because due to the absence of precise inertial and geometric parameters of the powertrain and the vehicle itself, as well as the characteristics of the mounts, more complex models can lead to major errors (Demic, 1990), (Frolov & Furman, 1990).
The powertrain, as a rigid body in space, has six degrees of freedom (three translations and three rotations), so, to describe its movement, it is necessary to use six independent generalized coordinates (Demic, 1990), (Frolov & Furman, 1990), (Genta, 2003), (Pars, 1981), (Gillespie, 1992), (Milliken & Milliken, 1995). Here are adopted geometric gravity axes for observing the movement of the powertrains, starting at an equilibrium position, which will be discussed later. It is worth mentioning that the use of the geometric gravity axis leads to the need of using centrifugal moments of inertia, but, for the sake of the ease of analysis, the assumption that they are also the principal axes of inertia is introduced. The powertrain performs spatial vibrations under the action of excitations it receives from the frame (originating from the road roughness, vibrations of the frame as an elastic system, inertial forces of braking, inertial force of acceleration, centrifugal forces in curvilinear movement and Coriolis forces due to the complex movement of the vehicle), as well as from the inertial forces and torque of the engine and the rotating masses.
Six generalized coordinates should be introduced to describe the spatial movement of the powertrain: three translations and three rotations (Demic, 1990), (Frolov & Furman, 1990), (Genta, 2003), (Pars, 1981), (Gillespie, 1992), (Milliken & Milliken, 1995). Two coordinate systems will be adopted (Demic, 1990), Figure 2:
• immobile (OXYZ), with the origin at the center of gravity of the powertrain (defined in relation to its equilibrium position). The axes of this system coincide with the geometric axes of the powertrain symmetry (it is at the same time the global coordinate system), • movable (C^rO, which is firmly attached to the powertrain and rotates with it. The transformation of the coordinate system from (CIC) to (OXYZ) and vice versa can be made using three independent angles (Euler angles - which lead to coupled differential equations), or angles of rolling, galloping and winding that lead to a
CO
o
CT>
.8
p. p
s,
o E
ira
rtr e
o p
s le cl
CD >
"¡0 O CD
E E o o
</) T3 CO
o
"co E
CD
O
o CO C1J </) CD
O
O
.Q
O
O
CO
t
CD
•O
'E
CD
Q
■■¡г ф
ю CD
О >
О
см
ОС LLJ
ОС ZD О
О _|
< с;
X О ш
н
>-
ОС <
(Л <
-J
О >о
X ш н
о
о >
scattered system of differential equations. Since the model used had only one angle- the rotation around the X axis for the angle it was not necessary to use three angles.
Figure 1 - Scheme of the powertrain mounts of a freight motor vehicle Рис. 1 - Схема опоры силового привода грузового автомобиля Слика 1 - Шема ослашаша погонске групе посматраног теретног моторног возила
Figure 2 - Introduced coordinate systems to describe the movement of the powertrain Рис. 2 - Внедрены системы координат для описания движения силового привода Слика - 2 Уведени координатни системи за описиваше креташа погонске групе
Since the aim of the research is to determine the thermal loads of the powertrain mounts, it was considered necessary to analyze only the bouncing and rolling of the powertrain, and in this case the translatory movement of the powertrain gravity center in space is defined by the coordinate Z and the rotation of the powertrain (as a rigid body ) around the axis % is defined by the rolling angle ф (these generalized coordinates create the conditions to include the vertical exitation from the vehicle frame and the resulting inertial force, as well as the engine torque into the analysis, while respecting the recommendation to apply as simple mathematical models as possible). The elastic-viscous forces perform mechanical work in the mounts (which is converted into heat energy) due to the relative movement of the powertrain relative to the vehicle frame (deformation of the mounts). The adopted mathematical model enables the analysis of the influence of the powertrain mounts on their thermal loads.
Fin - the resulting inertial force of the piston group engine torque
an attack point of the resulting inertial force and its coordinates in relation to the moving coordinate system the angle of installation of the powertrain (in this case it is 0°)
Me E
У
Figure 3 - Excitation due to inertia force and engine torque Рис. 3 - Возбуждение от инерционной силы и крутящего момента двигателя Слика 3 - Побуде од инерцц'алне силе и обртног момента мотора
Bearing in mind the above and using the Newton-Euler principle, differential equations can be written to describe the observed oscillatory movements of the powertrain in the form:
m • Z с = £ Z
Iu9 =
S Mf
(2)
(3)
In the literature (Igami et al, 2008), (Kennings et al, 2014), (Shanguan et al, 2016), (Simic & Demie, 1992) there are a large number of mathematical models of powertrain mounts, more or less complex. Since the goal of this paper was to compare the thermal loads of the mounts, it was considered necessary to simplify the expressions for the approximation of force in the mounts (for the projections in the direction of X, Y and Z axes). As far as reactive forces are concerned, they are opposed to movement.
The forces in elastic mounts are assumed in the form (Demie, 1990):
F = c, - Л. + c 2 -A2
cz zl г г 2
+ c.3 -Л3
(4)
where c-b On, ci2, and ci3 are the stiffness coefficients, and A is the relative deformation of the mount.
The viscous forces in the mounts are assumed in the form (Demic, 1990):
CO
о
CO
C\l oo
03 !± Ci
о E
2 -tE eu
о
Ci <л
о
eu >
«
'о eu E E о о
(Л
тз го
"rô Е
го eu
(Л
eu
о .а
с о О
"rô
eu
■о
'Е
eu
а
Faz= kl -A+ k2 -A • Slgn(A)
(5)
where k1, k2 are the damping coefficients, the relative deformation speed of the mount, and sign the corresponding mathematical function.
и
<1>
LO CD
O >
o CM
DC LLJ
DC ZD
o
o <
o 2:
X
o
LU
I— >-
QC <
CO <
-J
CD >0
X LU I— O z
—3
O >
The deformation and deformation projections of the powertrain mount in the Z direction are defined by the terms (Figure 2):
AZI = Z -Z0 -^(ç-ça),
A Zi = Z- Z0 - b (Ç-Ço ).
(6)
(7)
The vibrations of the powertrain also depend on the unbalance of the engine (torque and inertial forces).
In this particular case, a four-stroke four-cylinder regular diesel engine with a crankshaft, whose knees were at the same level (angle 180°), was used. The piston mechanisms present forces (Mihalec et al, 2015), which are driven by a piston (the forces of gases and the inertial forces of the piston group), and on the movable bearing of the crankshaft of the centrifugal and tangential force of the crankshaft. In balancing the inertial forces of the mass of the piston group (ideally, if the force develops into the Fourier series), there remain unbalanced inertial forces of the second and higher orders. It is noted that when there are differences in the masses of piston groups per cylinders, unbalanced first-order forces occur (in this case, the case where the masses have been equal to each other).
Assuming that high-order harmonics can be neglected, the unbalanced inertial force of the observed engine can be expressed in the form (Frolov & Furman, 1990), (Milliken & Milliken, 1995), (Rotenberg, 1972), (Genta, 2003), (Simic, 1988):
F. = 4 mr -o -A- coslat
(8)
where mr is the reduced mass of the piston group, r is the radius of the crankshaft knee, a is the angular speed of the engine crankshaft, A is the ratio of the radius of the crankshaft and the length of the piston engine, and t is time.
Based on elementary knowledge from the vector theory, statics, and Figures 2 and 3, the moment from the resulting inertial force defined by term (8) (Demic, 1990) is:
Mf,„ =
aL
-F,n sin y 0 F,n cos7
(9)
where the expressions in (9) are in accordance with Figure 3.
Centrifugal forces are partially counterbalanced by counter-strategies, or by other methods, for which readers are referred to (Mahalec et al, 2015).
Heat loads of the powertrain mounts
x T
A = {Fm-dzrel = { Fm-zrei-dt, (12)
The tangential force causes the torque of the engine, which due to its change has a variable value (unevenness is partially reduced by the flywheel) (Mahalec et al, 2015). S
In the absence of precise data, we assume that the torque that affects the powertrain can be described by the term (Demie, 1990):
M = -Me- i0- im (0,95+0,1-rnd), (10) §
where Me is the engine torque, and io is the transmission ratio in the differential, im is the transmission ratio in the transmission, and md random numbers are evenly distributed in the interval [0,1].
The vibrations of the powertrain are also influenced by the vibrations of the vehicle frame, which are incidental (Demic, 1997a), (Demic, 1997b), (Shanguan et al, 2016). Since the complexity of the vehicle spatial model exceeds the scope of this paper, it was considered appropriate not to use frame excitations obtained based on a model but to adopt broadband functions of the exitation in the following form:
exitation = max(md - 0.5), (11)
where max = 0.01 m, is work, and md has the same meaning as with the torque of the engine.
The projections of the generalized forces include all the components | of the force and momentum of the corresponding mounts in the direction of the observed axis (for mounts 1 to 4), the inertial forces and torque moments of the engine, or the unbalanced inertial forces, as well as the force of the active suspension which, due to the assumption that it operates at the powertrain center, has no moment.
o
o
.Q
Due to the relative movement of the powertrain and the vehicle frame, ° work is carried out in the mounts, which is equivalent to the amount of heat Q (Bojic, 2011), (Ilic et al, 1996), (Moran, 2010), (Fermi, 2011), (Demie, n 2001), (Demie, 2013), (Demie & Diligenski, 2016), (Mitschke & Wallento-witz, 2004). Mechanical work in mounts is defined by the expression (Mitschke & Wallentowitz, 2004):
<1J
where Fm is the elasto-damping force in the mount, zrel , ze is the relative
deformation and the deformation rate of the mount (respectively), and t is the time of deformation. > The mean mechanical strength is given by the expression:
T
o CM
IT LLI
o >
n A k f ,
P = — = — J zrel • dt. (13)
T T
^ As the effective value of the relative deformation speed of the mount is:
O 2 i T
o
. 2 \ .
_ z relef = — J zrel • dt . (14)
o T 0
x then the mean strength can be written in the form:
o
LU 2
^ P = k • z relef (15)
a: J ■
^ It turns into heat, with dominant convection:
P = «• S •Ar, (16)
where P is the power (heat flux), W; a is the coefficient of heat transfer, < W/(m2-°C); S is the convection surface in m2; and Ar is the temperature d difference between the mounts and the ambient air in °C. g Since a and S are of unknown values and require very complicated
research in order to be defined, they are not discussed in this paper. J Bearing in mind that all four mounts had the same characteristics, we
¡5 considered appropriate to carry out an analysis of total heat loads.
Numerical simulation and analysis of results
From expressions (2-11), it can be established that the differential equations are non-linear and of second order. Since the Runge-Kutta method requires them to be reduced to differential equations of the first order before solving them, they are transformed into this form. In differential equations 20 parameters were used to describe the characteristics of the mounts, two inertial parameters, three parameters describing translatory motion and three parameters describing angular motion (movement, velocity, acceleration), 12 coordinates of the mounts and three coordinates of the attack point of the resulting inertial force. Differential equations were solved using the software that Demic, M. had developed in Pascal. The integration was performed with a time step of 0.0001 s at 524288 points. This enabled a reliable analysis in the range
0.019-5000 Hz. Obviously, this range is considerably broader than the exi-tation range of engine motions and the complete powertrain of vehicles. The integration of differential equations was performed for the case of classical and hydraulic mounts.
The parameters of the vehicle and its powertrain are given in Table 1, and the coordinates of the connecting points (mounts) in Table 2 (FAP, 2017).
The mechanical and hydraulic mounts had identical stiffness values in the direction of the X, Y and Z axes, which were determined in FAP (FAP, 2017), and shown in Table 3.
Table 1 - Typical parameters of FAP1314 and its powertrain Таблица 1 - Типичные параметры ФАП1314 и его силового привода Табела 1 - Карактеристични параметри возила ФАП1314 и иегове погонске групе
Maximum engine power, kW 100
Maximum engine speed, min-1 2600
Maximum engine torque, Nm 428
Speed at maximum torque, min-1 1300
Vehicle mass, kg 12000
Mass of the powertrain, kg 1680
Moments of inertia Ix/Iy/Iz, kgm2 85/35/72
Transmission ratio in the main gearbox, - 3.83
Transmission ratio in direct transmission, - 1
Table 2 - The coordinates of the connection points of mounts Таблица 2 - Координаты точек соединения опоры Табела 2 - Координате везних тачака ослонаца
Position of mounts ^ a, m b, m c, m
Mount 1 0.5 0.4 0.1
Mount2 0.5 -0.4 0.1
Mount3 -0.5 0.4 0.1
Mount4 -0.5 -0.4 0.1
СО
о со
С\1 оо оо !± ü
о Е
2 -tE ф
о
<Л
о
Jz ф
>
о ф
Е Е о о
(Л
тз го
"то Е
го ф
(Л ф
о
.а +-<
с о О
"то ф
'Я 'Е ф Q
Table 3 - Stiffness of the used conventional and hydraulic mounts Таблица 3 - Жесткость используемых стандартных и гидравлических опор Табела 3 - Крутости коришПених класичних и хидрауличких ослонаца
ci1, N/m ci2, N/m2 ci3, N/m3
X 1200000 250000 60000
Y 1200000 250000 60000
Z 1200000 250000 60000
■■¡г ф
ю CD
О >
О
см
ОС LLJ
ОС ZD
о
о _|
< с;
X О ш
I—
>-
ОС <
In the absence of more precise data, the damping characteristics of the mounts are approximately defined based on the stiffness of the mounts and the weight they carry (Genta, 2003), (Simic, 1988), and given in Table 4.
Table 4 - Assumed damping characteristics of the mounts Таблица 4 - Предполагаемые демпфирующие характеристики опоры Табела 4 - Претпостав^ене пригушне карактеристике ослонаца
kt1 (mechanical/hydraulic), Ns/m 2 2 ki2(mechanical/hydraulic), Ns /m
X 620/62000 1/100
Y 620/62000 1/100
Z 620/62000 1/100
For the sake of illustration, the vertical oscillatory exitations of the powertrain derived from the frame exitations are shown in Figure 4. It is obvious that the exitations are time-varying and high thermal loads of the powertrain mounts are to be expected. A similar commentary also applies when it comes to frame rolling exitations that are not shown here.
w <
-J
CD >o
X Ш I—
о z
о >
о
— O.DOO
S
г
®
> -0004
Time, s
Figure 4 - Illustration of the vertical excitations of the frame Рис. 4 - Иллюстрация вертикального возбуждения рамы Слика 4 - Илустрацща вертикалних побуда оквира
The total heat loads of the four-point system possessing elastic and damping properties (for all components of force and momentum) are calculated using expressions (2-15), and the results are shown in Figure 5, where the values for the heat amount are given in the logarithmic proportions.
By analyzing the numerical data and Figure 5, it can be found that the mechanical mounts suffer about 49 times less thermal load than the hydraulic ones (for a time period of 52 s, the classical is about 4.85 x 108 J, and the hydraulic one is 2.37 x 1010 J). This is understandable if one takes
into account that the hydraulic mount, in addition to the hysteresis in the tire, also has additional movement of the fluid within the mounts. Figure 5 shows that the amount of heat grows over time, and that, if thermal energy were not conducted from the mounts, they would face distortion of the shape and characteristics.
Figure 5 - Heat amount (mechanical energy) depending on the type of mounts Рис. 5 - Количество тепла (механическойэнергии) в зависимости от типа опоры Слика 5 - Промена количине топлоте (механичког рада) у ослонцима, зависно од типа ослонаца
СО
о со
С\1 оо оо !±
о Е
2 -tE ф
о
<Л
о
Jz ф
>
"го
о ф
Е Е о о
(Л
тз го
"го Е
го ф
(Л ф
о
.О +-'
С О
О 15 ф
■о
'Е ф Q
Figure 6 - The heat flux for hydraulic mounts Рис. 6 - Тепловой поток для гидравлических опор Слика 6 - Промена топлотног флукса код хидрауличких ослонаца
As noted, the heat flux for both types of mounts is calculated, and the flux for the hydraulic mounts is shown in Figure 6 for illustration. It is
ф obvious that it changes stochastically with time, so it is useful to calculate the effective values for both types of mounts, given in Table 5.
ю
со Table 5 - Statistical data on heat flux
Ъ Таблица 5 - Статистические данные о тепловом потоке
> Табела 5 - Статистички подаци о топлотном флуксу
Mounts Effective value of heat flux, W
Classic 8.261 106
Hydraulic 3.639 108
o CM
QC LLJ
QC
O o
-j The analysis of the data in Table 5 indicates that the heat flux is much
o higher in hydraulic mounts than in the conventional (rubber-metal) mounts _ of the powertrain. This indicates a higher possibility of disrupting the o characteristics of hydraulic mounts than those of conventional ones. That is why they are larger than the classic ones. It points to the fact that such £ high values are the product of rigorous exitations from the vehicle frame " used in the simulation. In practice, they are significantly lower and so are the actual heat loads of the powertrain mounts. However, the research carried out aimed to determine the ratio of the thermal loads of classic and hydraulic mounts, with the help of models, so the results can be adopted w as orientational, which is acceptable in the freight vehicle design phase. It ^ should be noted that the amplitudes of vibrations of the powertrains in the case of hydraulic mounts are significantly lower than in the case of
<
CD
o conventional mounts (Demic, 1990), (Simic & Demic, 1992), (FAP, 2017).
As the aim of the research was to determine which component of the w exitation (including the exitation forces of the inertial force and the torque o of the engine) had the greatest impact on the thermal loads of the powertrain mount, it was judged appropriate to clarify this problem further > on. It is known that a dynamic system with multiple input quantities and one output quantity can be shown as in Figure 7. If there is a connection between the input quantities, the system must first be decoupled (Bendat & Piersol, 1993), (Bendat & Piersol, 2010). As the oscillatory model of the powertrain had 4 oscillatory exitations (2 shifts and two exitations from resulting inertial forces and the torque of the engine), we will summarize the theory that enables the analysis of the influence of each individual exitation on the thermal load of the powertrain mounts,heat fluxin particular. Figure 8 shows the distribution pattern (Bendat & Piersol, 1993), (Bendat & Piersol, 2010).
Powertrain
suspension
system
Q - amount of heat or heat flux,
z0 - vertical exitation from the vehicle frame (x1),
(p - exitation from rolling the vehicle frame (x2),
Fin - second-order inertial force (x3), Me - engine torque (x4)
Figure 7 - Block diagram of the suspension system of the powertrain
Рис. 7 - Блок-схема системы подвески силового агрегата Слика 7 - Блок-дщаграм система за ослаъаъе погонске групе
Figure 8 - Decoupling scheme Рис. 8 - Развязывающая схема Слика 8 - Шема распрезаъа
Bearing in mind the theory from (Demic, 1990), (Bendat & Piersol, 1993), (Bendat & Piersol, 2010) and Figure 8, it can be written (all quantities are given in the complex plane-Fourier or Laplace transformation):
X, — X ,
x21! xi
-Ll2 •xl,
(17)
x3 2! x3 L13 • xi L23 • x21!
x43! x4 A-4 ^xi I21-x24! L34 • x3-2!
where x1, x2, x3, x4 are the input quantities,
-second conditional input with the eliminated influence of the first input; the third conditional input with the elimination effect of the
CO
о
CO
C\l oo oo ci
CP
о E
2 -tE eu
о
CP
(Л
о
cu >
« О
eu E E о о
<Л
тз го
lö Е
го си
(Л
си
о
.а +-<
с о О
"то
eu
■О
'Е
eu
а
21! ' 32! ' 43!
<1>
LO CD
O >
O
>o
first and second conditional inputs; the fourth conditional input with the elimination of the influence of the first and the conditional second and the third conditional inputs, respectively.
The transfer functions between the input ports can be displayed as follows:
O t x2 T X3 T X4
CN ' - 2 ' - ' - 4
iy Xi Xi Xi
LLI 111
s
XX
n T — 3'2! T — 43!
L\2 _ L\3 _ L\4 _'
X
(18)
X2\! X2\!
<
where L12, L13, L14, L23, L24 are the corresponding transfer functions i between the conditional input quantities. uj Note that portable functions indicate how much the output quantity is
£ increased or decreased in relation to the input quantity and how big is the phase delay between them.
On the basis of expressions (17) and (18), partial coherence functions can be calculated, and they indicate the coupling between the input and output signals (Demic, 1990), (Bendat, 1998), (Bendat & Piersol, 1993), (Bendat & Piersol, 2010).
<
« v 2 __[_SiQ(i-\)!
^ 'iQ (i\)!
v 2 __(19)
'iQ (i\)! „ n > V '
xx (i-\)^ QQ (i-\)!
where SiQ(i.1)! is the cross-spectra (the power cross-spectra is a complex OT quantity that represents the boundary value of the mathematical expectation of the products of the conjugated-complex amplitude of one signal > and the complex amplitude of the second signal - it enables the analysis of the spectral content of the interdependencies of the two signals) and Sxx(i-1)! and SQQ(i-1)! are the corresponding autospectra (the power autospectrum is a real quantity that makes the limit value of the mathematical expectation of the products of the conjugated-complex and complex amplitude of the observed signal- it enables the analysis of the spectral content of a single signal). The function of ordinary coherence shows whether there is a connection between the input and output quantities of the system, and the partial coherence functions are used to analyze the interconnections of the input and output quantities in the case when there is a connection between the input quantities (Demic, 2001), (Bendat & Piersol, 1993), (Bendat & Piersol, 2010).
Figure 9 - Partial coherence functions for hydraulic mounts: B-vertical excitation from the frame, C-excitation from the rolling of the frame with excluded influence of the vertical excitation, D-inertial force with excluded influence of the vertical excitation and rolling of the vehicle frame and E-motor torque with excluded influence of vertical excitation and rolling of the vehicle frame and the inertial forces of the engine Рис. 9 - Частичные функции когерентности гидравлических опор: B-вертикальное возбуждение от рамы, C-возбуждение от качения рамы без возможности влияния вертикального возбуждения, D-инерционная сила без
возможности влияния вертикального возбуждения и качения рамы транспортного средства и E- момент двигателя без возможности влияния вертикального возбуждения и качения рамы транспортного средства и сил инерции двигателя Слика 9 - Парцц'алне функци]е кохеренци за хидрауличке ослонце: B - вертикалне побуде од оквира, C - побуда од ва^аша оквира са иск^ученим утица]ем вертикалних побуда, D - инерци]ална сила са иск^ученим утица^ем побуда од вертикалних и ва^аша оквира возила и Е - момент мотора са иск^ученим утица^ем вертикалних побуда и ва^аша оквира возила и инерци]алне силе мотора
For a more detailed analysis, the partial coherence functions (expression 19) are calculated using the programs developed by the author in Pascal (Demic, 2001). For the illustration purposes, Figure 9 shows the calculated values for hydraulic mounts. The data analysis for the mechanical hydraulic mounts shown in Figure 9 shows that the coherent partial functions depend on the frequency and type of exitations. It is considered appropriate to carry out a more detailed analysis of the data of the partial coherence functions. In that sense, the minimum and maximum values of the coherent partial functions are calculated and given in Tables 6 and 7.
CO
о
CO
C\l oo oo !±
о E
2 -tE cu
о
(Л
о
cu >
ro
О
cu E E о о
<Л
тз го
го Е
го си
(Л
си
о
.а +-<
с о О
го си
■о
'Е
си
а
■■¡г ф
ю CD
О >
О
см
ОС LLJ
ОС ZD
о
о _|
< с;
X О ш
н >-
ОС <
(Л <
-J
О >о
X ш н о
о >
Table 6 - Minimum and maximum values of the partial coherence functions for mechanical mounts Таблица 6 -Минимальные и максимальные значения парциальных функций когерентности для механических опор Табела 6 - Минималне и максималне вредности парци]алних функци]а кохеренци за
механичке ослонце
min max
B 0.308 0.843
C 0.230 0.869
D 0.661 1.000
E 0.541 0.956
Table 7 - Minimum and maximum values of the partial coherence functions for hydraulic mounts Таблица 7 - Минимальные и максимальные значения парциальных функций когерентности для гидравлических опор Табела 7 - Минималне и максималне вредности парцц'алних функци^а кохеренци
min max
B 0.252 0.741
C 0.216 0.897
D 0.721 1.000
E 0.600 0.839
The analysis of the data from the given tables shows that the type of mounts has an effect on the partial coherence functions, and that their values are in the interval of 0.216 to 1, indicating that there is a connection between the input quantities (the exitation from the vehicle frame and the resulting inertial force and the engine torque) and heat flux (Bendat & Piersol, 1993), (Bendat & Piersol, 2010). The resulting inertial force and the engine torque have a greater impact on the thermal loads of the powertrain mounts.
The realized research was aimed at determining the relation between the thermal loads of the classical (rubber-metal) mounts and hydraulic mounts, with the help of a model, so that the obtained results can be adopted as a guide, which is necessary during the design phase of the conceptual design of a freight motor vehicle. It should be noted that hydraulic mounts have significantly better powertrain vibration damping parameters than the classical ones (Demic, 1990), (Simic & Demic, 1992).
Conclusion
Based on the research carried out, it can be concluded that with the help of mechanical models, the influence of the characteristics of the powertrain mounts on their thermal loads can be analyzed. The analyzes
have shown that hydraulic mounts are more thermally loaded than the conventional ones. Due to their heat loads and a complex construction, it is understandable why they are more rarely applied in freight motor vehicles. In order to verify the results,experimental research should also be carried out.
References
Atkins, P., 2010. The laws of thermodynamics: A very short introduction.OUP Oxford.
Bendat, J.S., 1998. Nonlinear Systems-Techniques and Applications.London: Wiley.
Bendat, J.S., & Piersol, A.G., 1993. Engineering applications of correlation and spectral analysis.New York: Wiley-Interscience.
Bendat, J.S., & Piersol, A.G., 2010. Random Data: Analysis and measurement procedures.London: John Wiley and Sons.
Bojic, M., 2011. Termodinamika.Kragujevac: University of Kragujevac -Faculty of Mechanical Engineering (in Serbian).
Demic, M., 1990. A contribution to the optimization of the position and the characteristics of passenger car powertrain mounts. International Journal of Vehicle Design, 11(1), pp.87-103.
Demic, M., 1997a. Identification of vibration parameters for motor vehicles. Vehicle System Dynamics, 27(2), pp.65-88. Available at:http://dx.doi.org/10.1080/00423119708969323.
Demic, M., 1997b. Optimizacija oscilatornih sistema motornih vozila.Kragujevac: University of Kragujevac - Faculty of Mechanical Engineering, monograph (in Serbian).
Demic, M., 2001. Computer programs "Demparkoh", "Analsigdem".
Demic, M., 2002. Neki aspekti istrazivanja „leprsanja" upravljackih tockova teretnog motornog vozila. Vojnotehnicki glasnik/Military Technical Courier, 50(3), pp.300-308. Available at: http://dx.doi.org/10.5937/vojtehg0203300D (in Serbian).
Demic, M.D., 2013. Prilog istrazivanju toplotnog opterecenja amortizera poluaktivnog sistema oslanjanja motornog vozila. Vojnotehnicki glasnik/Military Technical Courier, 61(3), pp.26-50. Available at: http://dx.doi.org/10.5937/vojtehg61-3572 (in Serbian).
Demic, M.D., & Diligenski, B.M., 2016. Numerical simulation of shock absorbers heat load for semi-active suspension system. Thermal Science, 20(5), pp.1725-1739. Available at: http://dx.doi.org/10.2298/TSCI150624005D.
- FAP, 2017.FAP Informacije. Priboj: FAP (in Serbian).
Fermi, E., 2011. Thermodynamics.New York: Dover Books on Physics.
Frolov, K.V., & Furman, F.A., 1990. Applied theory of vibration protection systems.Moscow: Mashinostroenie.
CO
o
CD
co
U)
O
E
co -tE cu
o
to
o
cu >
"co o cu E E o o
CO T3 CO
o
"co E cu
o
o CO
cu
CO
cu o
o
.Q
o
o
CO
cu
•o
'E cu
0
ю
"J Genta, G., 2003. Motor Vehicle Dynamics.Singapore: World Scientific
Publishing Company.
Gillespie, T.D., 1992. Fundamentals of Vehicle Dynamics.Warrendale, PA: Society of Automotive Engineers. £ Igami, H., Thompson, M., & Osodo, K., 2008. Automotive Powertrain
Mounting System Design Optimization. SAE Technical Paper, 2008-01-0879. o Available at: http://dx.doi.org/10.4271/2008-01-0879.
Ilic, G., Radojkovic, N. & Stojanovic, I., 1996. Termodinamika II - Osnove w prostiranja toplote.Nis: University of Nis - Faculty of Mechanical Engineering (in is Serbian).
§ Kennings, P., Layfield, J., Tarabra, M., Fothergill, D., Syred, F., & Franks,
^ G., 2014. Developing Powertrain Mounting Systems in the Virtual Engineering ICA World Using a Full Vehicle NVH Simulator. In: INTER-NOISE and NOISE-CON HN Congress and Conference Proceedings, Melbourne, Australia, pp.2263-2272. ^ Mahalec, I., Kozarac, D. & Lulic, Z., 2015. Konstrukcije motora.Zagreb:
FSB (in Serbian).
oc Milliken, W.F., & Milliken, D.L., 1995. Race Car Dynamics.Warrendale,
Pennsylvania: Society of Automotive Engineers.
Mitschke, M., & Wallentowitz, H., 2004. Dynamik der Kraftfahrzeuge.Berlin - Heidelberg: Springer (in German).
Moran, M.J., & et al, 2010. Fundamentals of Engineering AS Thermodynamics.Hoboken: John Wiley and Sons. (5 Pars, L.A., 1981. A Treatise on Analytical Dynamics.Oxford: Ox Bow Press.
2 Rotenberg, R., 1972. Podveska avtomobilja. Moscow: Mashinostroenie (in
Russian).
>o
x Shanguan, W.B., & et al, 2016. Design method of automotive powertrain
E mounting system based on vibration and noise limitations of vehicle level.
о
>
Mechanical Systems and Signal Processing, 76, pp.677-695. Available at:
О http://dx.doi.Org/10.1016/j.ymssp.2016.01.009.
SimiC, D., 1988. Dinamika motornih vozila.Belgrade: Naucna knjiga (in Serbian).
Simic, D., & Demic, M., 1990. Elasticno oslanjanje pogonske grupe. In: Kragujevac: MVM, monograph (in Serbian).
ПРИЛОЖЕНИЕ К ИССЛЕДОВАНИЮ ТЕПЛОВОЙ НАГРУЗКИОПОРЫ СИЛОВОГО ПРИВОДА ГРУЗОВОГО АВТОМОБИЛЯ
Мирослав Д. Демич, Ясна Д. Глишович Университет в г. Крагуевац, Инженерный факультет, Департамент автотранспортных средств и двигателей, г. Крагуевац, Республика Сербия
ОБЛАСТЬ: машиностроение, автомобили и двигатели о
ВИД СТАТЬИ: оригинальная научная статья
ЯЗЫК СТАТЬИ: английский <§
Резюме: ^
Под действием возбуждений, возникающих в результате микрошероховатости дорог, трансмиссии, а также условий движения е (ускорение, торможение, криволинейное движение), массы коммерческих транспортных средств совершают сложные колебательные движения. Негативное воздействие этих перемещений, а соответственно и динамических нагрузок на опорную систему транспортного средства, может быть снижено при соответствующем выборе положения и характеристик опорыI силового агрегата. Как показывает практика, при согласовании характеристик опоры силового привода чаще всего учитываются такие факторыI, как минимальный шум и вибрации ^ транспортного средства, в то время как проблеме тепловых нагрузок, которым они подвергаются при эксплуатации, практически не
о
2 -е ф
о
(Л О
ф >
о ф
Е
уделяется никакого внимания. Необходимо учитывать тот факт, что о повышенные тепловые нагрузки неизбежно приводят к повреждению опоры/ и ухудшению ее характеристик. Следует отметить, что силыI упруго-демпферных опор выполняют механическую работу, преобразуя о механическую энергию в тепловую. Анализ тепловых нагрузок на опорыI го может быть выполнен с помощью теоретических, экспериментальных и комбинированных методов. В настоящей работе представленыI результаты исследования, полученные на основании теоретических методов и упрощенной механической модели трансмиссии. В ходе исследования испытаны вибрации силового привода грузового автомобиля среднего класса производственной линейки ФАП 1314. С учетом того, что современные грузовые автомобили оснащены как стандартными механическими, так и гидравлическими опорами был проведен анализ тепловой нагрузки с помощью математического моделирования, применяемого в обоих случаях. Исследования проведеныI с целью определения соотношения тепловой нагрузки стандартных 5 механических и гидравлических опор с помощью математических моделей, таким образом, полученные результаты могут считаться ориентировочными, что является важным фактором на этапе проектирования коммерческого автотранспортного средства. Анализ показал, что количество тепла, вырабатываемого при применении гидравлических опор, значительно выше, чем в случае применения стандартных механических опор. Кроме того, расчет парциальных функций когерентности теплового потока показал, что их значения находятся в диапазоне от 0,216 до 1, что указывает на взаимозависимость между входными значениями и тепловым потоком. Кроме того, доказано, что инерционная сила и крутящий момент
Е
.о
•о
■■¡г ф
ю
оказывают большее влияние на тепловую нагрузку силовых агрегатов, чем вибрационное возбуждение рамы транспортного средства. Ключевые слова: возбуждение, коммерческий автотранспорт,
ю опоры, моделирование, силовой привод, теплота.
"о >
ПРИЛОГ ИСТРАЖИВА^У ТОПЛОТНИХ ОПТЕРЕЪЕ^А
° ОСЛОНАЦА ПОГОНСКЕ ГРУПЕ ТЕРЕТНОГ МОТОРНОГ ВОЗИЛА
О!
ш Мирослав Д. Деми^, Jасна Д. Глишови^
к Универзитет у Крагу]евцу, Факултет инженерских наука, Катедра за
О моторна возила и моторе, Крагу]евац, Република Срби]а о
< ОБЛАСТ: машинство, моторна возила и мотори
— ВРСТАЧЛАНКА: оригинални научни чланак
х JЕЗИКЧЛАНКА: енглески
0 ш
Сажетак:
>
Под деjством побуда од микронеравнина путева, погонске групе, као и од услова кретаъа (убрзаваъе, кочете, криволинц'ско кретак>е) масе теретних возила врше сложена просторна осцилаторна кретаъа. Негативан утицаj тих кретаъа, односно динамичких оптереЯеъа ще она изазиваjу и преносе на носеЯи ^ систем возила, може се смаъити правилним избором положаjа и
щ карактеристика ослонаца погонске групе.
2 У пракси се усаглашаваъе карактеристика ослонаца погонске групе,
на1чешАе, врши из услова минималне буке и вибрацца возила, при чему
1 се, углавном, занемаруу топлотна оптереЯеъа щима су они изложени
н у експлоатацщи, а зна се да она неминовно доводе до оштеЯеъа самих
~>
о Треба напоменути да еласто-пригушне силе у ослонцима врше
механички рад щи се претвара у топлотну енерг^у. Анализа то-плотних оптереЯеъа ослонаца може се вршити коришЯеъем тео-ретских, експерименталних и комбинованих метода. У овом раду истраживаъа су вршена применом теор^ских метода уз коришЯе-к>е упрошЯеног механичког модела погонске групе. При томе су по-сматране вибрац^е погонске групе теретног моторног возила средне класе из производног програма ФАП 1314. ИмаjуЯи у виду заступл>еност класичних (гумено-металних) и хидрауличких осло-наца погонске групе код савремених теретних моторних возила, у раду су извршене анализе топлотних оптереЯеъа поменутог возила, уз коришЯеъе математичког модела.
Реализована истраживаъа имала су за цил> да се утврди однос топлотних оптереЯеъа класичних и хидрауличких ослонаца, уз помоЯ модела, па се доб^ени резултати могу усвоjити као
ори]ентациони, што }е неопходно у фази израде иде]ног про]екта о
см
личина топлоте код хидрауличких ослонаца вишеструко веЬа него га у случа}у коришПеша класичних ослонаца. £
Поред тога, израчунате парци]алне функци]е кохеренци топлот-ног флукса су показале да су шихове вредности у интервалу 0,216 о до 1, што указу}е на то да посто]и веза измену улазних величина и топлотног флукса. При томе, веЬи утица] на топлотна оптере- | Пеша ослонаца погонске групе има]у инерци]ална сила и обртни момент мотора него осцилаторне побуде оквира возила. Къучне речи: побуда, комерци]ална возила, ослонци,
Е
с го -tE <u
о
(Л
моделираъе, погонска група, топлота. ъ
ф
>
Paper received on / Дата получения работы / Датум приема чланка: 15.05.2017. -¡g
Manuscript corrections submitted on / Дата получения исправленной версии работы / Датум о
доставлажа исправки рукописа: 07.08.2017. ^
Paper accepted for publishingon / Дата окончательного согласования работы / Датум коначног ||
прихватажа чланка за об]авливаже: 09.08.2017. о
(Л
тз го
го Е
© 2017 The Authors. Published by Vojnotehnicki glasnik / Military Technical Courier (www.vtg.mod.gov.rs, втг.мо.упр.срб). This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license о
(http://creativecommons.org/licenses/by/3.0/rs/).
© 2017 Авторы. Опубликовано в «Военно-технический вестник / Vojnotehnicki glasnik / Military Technical Courier» (www.vtg.mod.gov.rs, втг.мо.упр.срб). Данная статья в открытом доступе и распространяется в соответствии с лицензией «CreativeCommons» (http://creativecommons.org/licenses/by/3.0/rs/).
го
© 2017 Аутори. Обjавио Воjнотехнички гласник / Vojnotehnicki glasnik / Military Technical Courier <л (www.vtg.mod.gov.rs, втг.мо.упр.срб). Ово jе чланак отвореног приступа и дистрибуира се у складу са Creative Commons лиценцом (http://creativecommons.org/licenses/by/3.0/rs/).
о
.а +-<
с о О
"го ф
■о
'Е ф
а