Prilog istraživanju toplotnog opterećenja amortizera poluaktivnog sistema oslanjanja motornog vozila Текст научной статьи по специальности «Медицинские технологии»

VOJNOTEHNIČKI GLASNIK/MILITARY TECHNICAL COURIER, 2013., Vol. LXI, No. 3

PRILOG ISTRAŽIVANJU TOPLOTNOG OPTEREĆENJA AMORTIZERA POLUAKTIVNOG SISTEMA OSLANJANJA MOTORNOG VOZILA

Miroslav D. Demić

Fakultet inženjerskih nauka, Kragujevac

DOI: 10.5937/vojtehg61-3572

OBLAST: motorna vozila

VRSTA ČLANKA: originalni naučni članak

Sažetak:

U razvoju motornih vozila značajnu ulogu ima dinamička simulaci-ja, koja je zasnovana na modeliranju. Uloga modeliranja je veoma zna-čajna u prvim fazama projektovanja, tokom definisanja opredeljujućih parametara vozila. Amortizer, kao izvršni organ poluaktivnog sistema za oslanjanje vozila, trpi toplotna opterećenja koja mogu da dovedu do njegovog oštećenja i degradacije karakteristika. Zbog toga je u ovom radu analizirano pretvaranje mehaničkog rada u toplotnu energiju, uz korišćenje metode dinamičke simulacije.

Ključne reči: vozilo, poluaktivni sistem elastičnog oslanjanja, amortizer, toplotno opterećenje amortizera.

Amortizeri su sastavni element poluaktivnog sistema za oslanjanje v ozila (u daljem tekstu PASeO), koji neposredno utiču na aktiv-nu bezbednost vozila. Pored osnovne uloge, uloga amortizera je da ap-sorbuje mehaničke oscilacije prenesene sa podloge i da obezbedi bezbednost putnika u toku vožnje.

Kao što je poznato, u amortizeru se kinetička energija sistema pretvara u mehanički rad, odnosno toplotnu energiju (Genta, 2003), (Miliken and Mili-ken, 1995), (Mohd, et al, 2008), (Mitschke, 1997), (Reisenburg, 1970), (Si-mid, 1998). U praksi se, u fazi projektovanja vozila, parametri amortizera bi-raju iz uslova prigušivanja oscilacija sistema vozila, ali se zbog izbegavanja negativnog uticaja na njihovu funkciju vodi računa i o termičkim opterećenji-ma amortizera (Mitschke, 1997). Cilj je da se što je moguće više mehanič-

* ZAHVALNICA: Autor se zahvaljuje Ministarstvu za obrazovanje, nauku i tehnološki razvoj Republike Srbije za finansijsku podršku (Naziv projekta: Istraživanje bezbednosti vozila kao dela kibernetskog sistema: vozač - vozilo - okruženje; Oznaka projekta: TR35041).

Uvod

[email protected]

kog rada primljenog sa podloge pretvori u toplotnu energiju koja se predaje okolini i time amortizer hladi. Pogrešan izbor amortizera sa stanovišta ter-mičkog opterećenja može izazvati ubrzanu degradaciju njegovih karakteristi-ka tokom eksploatacije. Prevelika količina toplote koja bi se eventualno za-državala „unutar" amortizera izazvala bi ubrzano propadanje zaptivnih ele-menata i gubitak funkcije prigušnog elementa.

Ispitivanja su pokazala da se mehanički rad jednim delom pretvara u toplotnu energiju koja se prenosi na radni fluid i telo amortizera, a preostali deo količine toplote predaje se okolini i time amortizer hladi. To se matema-tički može prikazati formulom (Atkins, 2010), (Bojić, 2011), (Ilić i dr, 1996), (Marić, 2002), (Mohd, et al, 2008), (Moran, et al, 2010), (Šelmić, 1986):

A = Qt + Qf + Q0 (1)

gde su:

• A - mehanički rad (jednak količini toplote),

• Qt - količina toplote koja se predaje telu amortizera,

• Qf - količina toplote koja se predaje radnon fluidu i

• Qo - količina toplote koja se predaje okolini.

Rad sile u amortizeru interesantan je zbog analize njegove transfor-macije u toplotnu energiju. Ovaj rad sile u amortizeru je merljiv i eksperi-mentalno, ali je teško merljiva količina toplote koja se oslobađa iz amortizera (Fermi, 2011). Ta pojava je kompleksna i teško merljiva, jer je pozna-to da se deo energije troši za zagrevanje elemenata amortizera, kao što su klip, cilindri, radni fluid, itd. Pored toga, i priroda prenosa toplote iz amortizera na okolinu je veoma složena. Prenos toplote vrši se konvekcijom koja je dominantna, ali i kondukcijom i zračenjem (Mitschke, 1997). Sa aspekta maksimalnog hlađenja, pravilan izbor amortizera zahteva sveobuhvatnu analizu transformacije mehaničke energije u toplotnu. Način transformacije mehaničke energije u toplotnu najvećim je delom određen samom kon-strukcijom amortizera. Na kondukcioni prelaz toplote i zračenje toplotne energije iz amortizera se ne može direktno uticati. Potrebno je povećati uti-caj na konvekcioni prelaz toplote sa amortizera na okolnu sredinu, kao do-minantne pojave. To se najbolje postiže konstruktivnim rešenjem karoseri-je (blatobrana) motornog vozila. Zahtev je da se iskoristi opstrujavanje va-zduha oko amortizera uz najmanje usložavanje konstrukcije.

Napominjemo da cilj ovog rada nije bio da se analizira hlađenje amortizera, već samo termičko opterećenje kojem je on izložen. Zbog toga je oce-njeno celishodnim da se izvrši analiza toplote koja se dobija pretvaranjem mehaničkog rada u toplotnu energiju u jedinici vremena. Mehanički rad u amortizerima poluaktivnog PASEO izračunavan je uz pomoć mehaničkog modela vozila, o čemu će biti reči u narednom tekstu.

CiD

Demić, M., Prilog istraživanju toplotnog opterećenja amortizera poluaktivnog sistema oslanjanja motornog vozila, pp. 26-50

VOJNOTEHNIČKI GLASNIK/MILITARY TECHNICAL COURIER, 2013., Vol. LXI, No. 3

Model vozila sa poluaktivnim sistemom za oslanjanje

Motorna vozila imaju složene dinamičke karakteristike, koje se mani-festuju pojavom prostornih kretanja, promenom karakteristika u toku eks-ploatacije, velikim brojem poremećajnih uticaja, pojavom zazora, trenja, hi-sterezisa i sl. (Belingardi and Demić, 2009, pp.45-53), (Demić, 1997a), (Demić, 2004), (Demić, 1999, pp.183-211), (Demić, 1994, pp.351-377), (Demić, 1996, pp.22-46), (Demić, et al, 2006, pp.145-158), (Demić and Belingardi, 2010, pp.6-17), (Demić and Diligenski, 2012, pp.185-190), (De-mić, 1997b, pp.135-146), (Genta, 2003), (Gillespie, 1992), (Miliken and Mi-liken, 1995), (Rajamani, 2005). Pomenute dinamičke pojave, a posebno vi-bracije, dovode do zamora korisnika i vozača, smanjuju vek trajanja vozila i njegovih sistema i sl. Osnovni cilj sistema za oslanjanje jeste da smanji pomenute negativne efekte, poboljša ponašanje vozila na putu i omogući eksploataciju vozila u širem spektru eksploatacionih uslova.

Klasični sistemi to ne mogu da zadovolje, pa se javila potreba za uvođenjem novih sistema za oslanjanje, sa regulisanim karakteristikama (ukratko nazvanim poluaktivni ili aktivni sistemi). Njihove podele i klasifi-kacije poznate su iz (Gillespie, 1992), (Popović, 2001), pa o tome ovde neće biti više reči.

Ukazuje se na činjenicu da klasični «pasivni» sistemi mogu ponuditi samo kompromis između oprečnih zahteva za oscilatornom udobnošću i ponašanjem vozila na putu (Belingardi and Demić, 2009, pp.45-53), (De-mić, 1997a), (Demić, 2004), (Demić, 1999, pp.183-211), (Demić, 1994, pp.351-377), (Demić, 1996, pp.22-46), (Demić et al, 2006, pp.145-158), (Demić and Belingardi, 2010, pp.6-17), (Demić and Diligenski, 2012, pp.185-190), (Demić, 1997b, pp.135-146), (Genta, 2003), (Gillespie, 1992), (Miliken and Miliken, 1995), (Rajamani, 2005), jer su im parametri krutosti i prigušenja fiksni. Imajući u vidu težnju za manjom potrošnjom goriva, posebno kod vozila srednje klase, a što se, pored ostalog, postiže i manjom sopstvenom masom vozila, očigledno je da se aspekt «optere-ćeno/neopterećeno» stanje manifestuje degradacijom parametara oscila-torne udobnosti i ponašanja vozila na putu.

Otuda se može zaključiti da samo sistemi oslanjanja sa regulisanim karakteristikama mogu da na zadovoljavajući način kompenzuju poreme-ćajne uticaje, a, istorijski posmatrano, prvi takav sistem je korišćen kod vasionskih letelica («Sputnik», između 1950. i 1960. godine, (Silani,

2004), dok kod motornih vozila njihova primena počinje 70-ih godina (Silani, 2004) i to najpre u «poluaktivnom» izvođenju.

U praksi postoje dve vrste sistema sa regulisanim karakteristikama, i to (Demić et al, 2006, pp.145-158), (Doule et al, 1990), (Florin et al,

2004), (Hac et al, 1999, pp.1-7), (Hrovat et al, 1982, pp.1-9), (Margolis, 1983, pp.317-330), (Nell et al, 1998, pp.25-40), (Popović, 2001), (Silani, 2004), (Slaski et al, 2004, pp.597-604), (Rajamani, 2005), (Tomović, 1974), (Yoshimura et al, 2003, pp.22-40):

1. «poluaktivni», kod kojih se vrši regulacija krutosti ili prigušenja; imaju povoljne karakteristike, a cena im je prihvatljiva i za vozila niže kategorije i

2. «aktivni», koji imaju generator sile (aktuator) koji može biti kombi-novan i sa klasičnim elastoprigušnim elementima; cena im je veo-ma visoka, pa primenu nalaze kod vozila veoma visoke klase ili kod specijalnih vozila.

Kako postoje različita konstruktivna izvođenja pomenutih sistema, u ovom radu će biti reči o poluaktivnom sistemu oslanjanja vozila sa regu-lacijom prigušenja u amortizerima i njihovom termičkom opterećenju.

Imajući u vidu da u literaturi postoji više algoritama za regulisanje prigušenja u amortizerima, ocenjeno je celishodnim da se neki od njih ov-de prikažu u najkraćim crtama.

U (Florin, 2004) kontrolisana sila u amortizerima definiše se izrazima:

Fa = k(t)(Z2 - Z1) (2)

pri čemu je njena veličina:

Fa

0 V- ksZ2(Z2 - Z1)

< kmax (Z2 - Z1 ) V - ksZ2(Z2 - Z1) <-kmax (Z2 - Z1)2

ksZ2 za ostale slučajeve

gde su:

• k(t) - kontrolisana vrednost koeficijenta prigušenja,

• kmax - maksimalna vrednost koeficijena prigušenja,

• ks - koeficijent u povratnoj sprezi (skyhook),

(3)

• zi, Z2 - brzine neoslonjene i oslonjene mase (za slučaj četvrtinskog modela vozila), respektivno.

U (Silani, 2004) kontrolni signal prigušenja u amortizeru definiše se:

\kd V Z2(Z2 - Z1) > 0

lkmm V Z2(Z2 - Z1 ) ^ 0

(4)

pri čemu je:

Demić, M., Prilog istraživanju toplotnog opterećenja amortizera poluaktivnog sistema oslanjanja motornog vozila, pp. 26-50

VOJNOTEHNIČKI GLASNIK/MILITARY TECHNICAL COURIER, 2013., Vol. LXI, No. 3

kd = i

kmax V ks-------— > k

k ——— V k < k ——— < k

s . . p s . . max

kp V ks

< km

(5)

gde je kp - nepromenljiva vrednost koeficijenta prigušenja, dok su ostale oznake identične onima iz (3).

Algoritam opisan u (Silani, 2004) definiše koeficijent prigušenja u obliku:

km

km

V Z2izm Z2max < Z‘2izm z2min

V z 2 ■ * z 2 ■ < z 2 ■ * z 2

2izm 2min 2izm 2max

(6)

gde su:

• z2izm ’ z2min ’ z2max - izmereno ubrzanje i ubrzanja oslonjene mase za mi-nimalno i maksimalno prigušenje u amortizeru, respektivno.

Napominje se da se umesto z2max,z2min u izrazu (6) često koristi od-govarajuća sila u amortizeru.

Regulacija prigušenja u amortizeru (Silani, 2004) definisana je algo-ritmom:

«0. 1 II (7)

pri čemu važi relacija :

0 < k{t) < kmax (8)

Pored toga, u (Silani, 2004) daje se i mogućnost linearne aproksima-cije izraza (8).

Neki autori (Nell et al, 1998, pp.25-40), (Silani, 2004) preporučuju korišćenje sistema sa preklopnicima ("uključeno”- “isključeno”), čija je ka-rakteristika jednostavnost konstrukcije i niska cena. Imajući u vidu eko-nomsku situaciju u našoj zemlji, kao i tehnološki nivo proizvođača kom-ponenti i vozila, ocenjeno je celishodnim da se u ovom radu ova koncep-cija regulacije sile u amortizeru detaljnije analizira.

Optimalno projektovanje sistema poluaktivnog oslanjanja biće de-monstrirano na primeru nelinearnog prostornog modela vozila.

<JD

Oscilatorni model vozila

Imajući u vidu cilj istraživanja, ocenjeno je celishodnim da se posma-traju sledeća sekundarna (oscilatorna) kretanja masa vozila (slika 1) .

- Oslonjena masa:

• z - vertikalne oscilacije oslonjene mase (poskakivanje),

• ф - valjanje i

• 0 - galopiranje

- Vertikalne vibracije neoslonjenih masa (točkova):

• zpl - prednji levi,

• zpd - prednji desni,

• zzl - zadnji levi i

• zzd - zadnji desni.

Diferencijalne jednačine oscilatornih kretanja prikazanog modela vozila napisane su uz primenu programskog paketa NEWEUL. Priprema ulaznih po-dataka izvršena je u skladu sa zahtevima pomenutog programa, a za defini-sanje dinamike modela vozila sa sedam stepeni slobode kretanja korišćeno je 24 koordinatnih sistema (uključujući i inercijalni-globalni koordinatni sistem OI1I2I3), u odnosu na koje su definisani inercijalni parametri, sile i sl. Imajući u vidu da je postupak pripreme ulaznih podataka detaljno opisan u (NEWEUL, 2000), (Sciehlen, 1986, pp.271-288), ovde o tome neće biti više reči.

Treba napomenuti da su generalisane koordinate oslonjene mase (z, ф i 0) generalisane koordinate poskakivanja masa točkova (zpl, zpd,zzl i zzd) defi-nisane od njihovog ravnotežnog položaja, pa sile teže posmatranih masa ne figurišu u diferencijalnim jednačinama kretanja (Demić, 1997a), (Demić, 2004), (Demić, 1999, pp.183-211), (Demić, 1994, pp.351-377), (Demić, 1996, pp.22-46), (Demić et al, 2006, pp.145-158), (Demić and Belingardi, 2010, pp.6-17), (Demić and Diligenski, 2012, pp.185-190), (Demić, 1997b, pp.135-146), (Genta, 2003), (Gillespie, 1992), (Miliken and Miliken, 1995), (Rajamani, 2005). Koristeći programski paket NEWEUL i sliku 1, napisane su nelinearne diferencijalne jednačine kretanja posmatranog modela vozila, (NEWEUL, 2000), (Sciehlen, 1986, pp.271-288):

Mq + K = QE (9)

gde su:

• M - matrica inercijalnih parametara (SSKxSSK - SSK broj stepeni slobo-

de kretanja),

• K - matrica Koriolisovih i centrifugalnih sila (1xSSK),

• QE - matrica generalisanih sila, •

• q,q - matrice generalisanih koordinata i njihovih ubrzanja, respektivno

(1xSSK).

Demić, M., Prilog istraživanju toplotnog opterećenja amortizera poluaktivnog sistema oslanjanja motornog vozila, pp. 26-50

VOJNOTEHNIČKI GLASNIK/MILITARY TECHNICAL COURIER, 2013., Vol. LXI, No. 3

Inercijalne parametre posmatranog modela vozila definišu veličine:

• Gs, ms, mp, mz, gravitaciona sila i oslonjena masa, neoslonjene mase sistema prednjih i zadnjih točkova, respektivno,

• Ix moment inercije oslonjene mase za «x» osu,

• Iy moment inercije oslonjene mase za «y» osu, a

• Iz moment inercije oslonjene mase za «z» osu.

U slučaju posmatranog prostornog modela vozila, diferencijalne jednačine koje opisuju oscilatorna kretanja su oblika:

z = (QE[1] - k[1 ])/ms 9 = (QE[2\ - k[2]) /(Ix cos2 в + 2Iz sin2 в) в = (QE[3\ - k[3])/Iy

zPi = (QE[4] - k[4])/ mp (10)

= (QE[5] - k[5])/ mp Zzi = (QE[6] - k[6])/ mz zzd = (QE[7] - k[7])/ mz

Matrice K, i QE date su izrazima:

K [1] = 0

K [2] = 2( Iz -1x )рв sin(e) cos(e)

2

K [3] = (Ix - I2 )9 sin(e)cos(e)

K [4] = 0 (11)

K[5] = 0 K [6] = 0 K [7] = 0

C3D

от = -(Fupi + Fupđ + FUZI + Fuzd )cos(p) cos(0) QE[2] = F -Fupđ)sx cos(p)cos(0) + F, -Fuzđ )s2 cos(p)cos(0)

QE[3] = (Fup, + Fupđ)acos(0) + (F^ -

Fuzi )si sinp) sin(0)- (Fup, + Fupđ)b cos(0)

+ (Fuzđ - Fuzl )S2 sin(P)sin(0)

QE[4] = -Fpp, + Fupl cos(p) cos(0)

QE[5] = ~Fppđ +Fupđ cos(p)cos(0)

QE[6] = -Fuzl + Fuzl Cos(P)Cos(0)

QE[7] = -Fpzđ +Fuzđ cos(p)cos(0)

Pomeranja karakterističnih tačaka koja su značajna za izračunava-nje zbirnih sila u oprugama i amortizerima (Fupl, Fupd, Fuzi i Fuzd) automat-ski su izračunate istim programom. Ona su data izrazima:

A„, Д„ = asin(0)±, sin(p)cos(0) + -. с°(р)со*(0 (13)

А,, A zd = -b sin(0) ± s1 sin(p) cos(0) + z cos(p) cos(0)

Sada možemo izračunati i relativna pomeranja klipova u amortizeri-ma:

A*

A*

A*

A*

pl A pi zpi

pđ = A pđ z pđ

zl =A zl - Zzl zđ = A zđ Zzđ

(14)

Brzine pomeranja (A pl, A pđ, A z,, A zđ) mogu se izračunati diferen-ciranjem izraza (14) po vremenu, što ovde neće biti učinjeno.

Radi praćenja daljeg teksta posmatraćemo sliku 2, na kojoj je prika-zana predložena šema upravljanja poluaktivnim sistemom oslanjanja za posmatrani model vozila. U konkretnom slučaju razrađen je koncept upravljanja sistemom oslanjanja koja se ostvaruje preko dve grupe po-vratnih sprega. Zadatak spoljašnjih povratnih sprega jeste da minimizira-

(®>

Demić, M., Prilog istraživanju toplotnog opterećenja amortizera poluaktivnog sistema oslanjanja motornog vozila, pp. 26-50

VOJNOTEHNIČKI GLASNIK/MILITARY TECHNICAL COURIER, 2013., Vol. LXI, No. 3

ju uticaj mikroneravnina puta, a unutrašnjih da upravljaju silama u amorti-zerima. Preciznije rečeno, spoljašnje povratne sprege vrše kontrolu osci-latorne udobnosti vozila, a unutrašnje omogućavaju povoljnu raspodelu opterećenja i njegovo nivelisanje, tokom manevara vozilom (aspekt po-našanja vozila na putu). Sa slike je očigledno da je ovaj koncept poluak-tivnog oslanjanja baziran na registrovanju vertikalnih ubrzanja poskakiva-nja, ubrzanja valjanja i galopiranja oslonjene mase (koja su integracijom omogućila izračunavanje odgovarajućih brzina), kao i hodova klipova ci-lindara i njihovih odgovarajućih brzina.

Kako su sile u elastoprigušnim elementima (zbir sile u opruzi i amor-tizeru) Fupl , Fupd , Fuzl i Fuzd zavisne od poskakivanja, valjanja i galopiranja oslonjene mase, neophodno je razdvojiti komponente koje izazivaju po-skakivanje - z od onih koje izazivaju valjanje- ф i galopiranje- 0 (videti blok «Rasprezanje» definisan na slici 2). Na osnovu generalisanih sila Fz=QE[1], F<p=QE[2] i F0=QE[3] (izraz 12), uz pomente generalisane koor-dinate možemo napisati relacije koje dovode u vezu rezultujuće signale fz, f<p i f0 sa delovima signala upravljanja silama u amortizerima fpl, fpd, fzl, fzd u obliku:

fz -c<pc6 - cqcQ - cqcQ - cqcQ

f„ = sxcqc6 - sxcqc6 s2 cqcQ - s2 cqcQ

fo_ ac6 - sxsqs6 ac6 + sxsqs6 - bcd- s2 sqsQ - bc6 + s2 sqsQ

fPi

fpd

fzl

(15)

gde su radi lakšeg pisanja uvedene skraćenice s(.)=sin(.) i c(.)=cos(.), a dimenzije matrice su 4x3.

Pri tome treba naglasiti da u izrazu (15) figurišu samo sile u amortizerima, jer se njihovi parametri kontrolišu, dok su parametri opruga fiksne veličine.

Izraz (15) možemo napisati i sažetije, u matričnom obliku:

F = AFf

(16)

odakle sledi:

Ff = a--f

(17)

C3D

Slika 1 - Oscilatorni model vozila Figure 1 - The oscillatory model of the vehicle

Iz linearne algebre poznat je postupak izračunavanja pseudoinverz-ne matrice A-1

A-1 = AT (AAT )-1 (18)

gde je AT transponovana matrica matrice A. Imajući u vidu složenost matrice A (izraz 14), pseudoinverzna matrica izračunata je uz korišćenje programa Mathematica (http://www.wolfram.com, nd), tako da je dobijena relacija:

b (s2 + s2 )spsd + s1Led 1

Loped 2L(s12 + s22)epe2d 2Led

4^ b (s12 + s2)spsd- s1Led 1

fpd Loped 2 L(s2 + s2)epe 2d 2 Led

fzi a (s12 + s22 )spsd - s2 Led 1

fzd _ Loped 2 L(s2 + s2)epe 2d 2 Led

a (s12 + s2 )spsd + s2Led 1

Loped 2L(s12 + s2)epe 2d 2Led

(19)

gde fzd, fzi, fpd i fpl (e7, e8, e9 i e10) predstavljaju delove upravljačkih sig-nala sila u amortizerima.

Demić, M., Prilog istraživanju toplotnog opterećenja amortizera poluaktivnog sistema oslanjanja motornog vozila, pp. 26-50

VOJNOTEHNIČKI GLASNIK/MILITARY TECHNICAL COURIER, 2013., Vol. LXI, No. 3

Slika 2 - Blok-dijagram posmatranog poluaktivnog sistema oslanjanja Figure 2 - Block sheme of the semi-active suspension system

Na osnovu usvojenog koncepta upravljanja sila u amortizerima i slike 2. mogu se napisati sledeće relacije:

e1 = ^želj -б

e2 = e1x[1]

e3 = Vželj -V

e4 = e3 x[2]

e5 = zžeij - ž

e6 = e3x[2]

fv = e2 - x[6](P

fe = e4 - x[5]e fz = e5 - x[4]ž

e7 = fd

e8 = fzi

e9 = fpd e10 = fpl

ell = e7 - x[7] e12 = e8 - x[8] e13 = e9 - x[9] e14 = e10 - x[10] e15 = step (e11) e16 = step (e12) e17 = step (e13) e18 = step (e14)

(20)

Ukazuje se na činjenicu da veličine e7-e18 predstavljaju upravljačke veličine za sile na zadnjem desnom, zadnjem levom, prednjem desnom i prednjem levom amortizeru (fazd,fazl,fapd,fapl), respektivno (videti sliku 2). U ovom radu su, imajući u vidu usvojeni koncept upravljanja prigušenjem i izraz (2), sile u amortizerima definisane relacijama:

Fazd = e15 • Aaz • (Zzd - Ztzd )

Fazl = e16 • Aaz • (Zzl - Ztzl )

Fapd = • Aap • (Z pd - Ztpd )

Fapl = e18 • Aap • (Zpl - Ztpl )

gde smo sa Aap i Aaz označili maksimalne površine otvora kroz koje struji ulje u prednjim i zadnjim amortizerima.

Imajući u vidu napomene u vezi s korišćenjem preklopnika "uključe-no” - "isključeno”, usvojena je zavisnost upravljanja koeficijentom prigu-šenja u amortizerima oblika:

k = Jkmax V Zz (- ^^t,z ) > 0

= |o V Zz (Zz - Ztz) < 0 (22)

z = pl, pd, zI , Zd

Na osnovu slike 2. očigledno je da veličine e15-e18 imaju vrednosti 0 ili 1, a u zavisnosti od odgovarajućeg proizvoda u izrazu (22) .

Nelinearne sile u oprugama koje nisu objekt upravljanja definišu iz-razi (Belingardi and Demić, 2009, pp.45-53), (Demić, 1997a), (Demić, 2004), (Demić, 1999, pp.183-211), (Demić, 1994, pp.351-377), (Demić, 1996, pp.22-46), (Demić et al, 2006, pp.145-158), (Demić and Belingardi, 2010, pp.6-17), (Demić and Diligenski, 2012, pp.185-190), (Demić, 1997b, pp. 135-146), (Genta, 2003), (Gillespie, 1992), (Miliken and Mili-ken, 1995), (Rajamani, 2005):

Fozd = 41 1] • (Azd - Zzd ) + 412] • (Azd - Zzd )3

Fozl = X[1 1] • (Azl - Zzl ) + ^[12] • (Azl - Zzl )3

3 (23)

Fopd = x[1 3] • (Apd - Zpd ) + x[1 4] • (Apd - Zpd )

Fopl = x[1 3] • (Apl - Zpl ) + x[14] • (Apl - Zpl )3

Demić, M., Prilog istraživanju toplotnog opterećenja amortizera poluaktivnog sistema oslanjanja motornog vozila, pp. 26-50

VOJNOTEHNIČKI GLASNIK/MILITARY TECHNICAL COURIER, 2013., Vol. LXI, No. 3

Radijalne sile u pneumaticima možemo izraziti na osnovu njihovih radijalnih deformacija (Demić, 1997a):

Ft = c 1 A + c 2 A 2 - c 2 A 3

t, s pi, о 4s p 2,s s p 2,s s

s = pl, pd, zl, zd.

gde su deformacije pneumatika date izrazima:

A = z — z

s s o,s

s = pl, pd, zl, zd.

(24)

(25)

Dinamičke procese sistema opisuju spregnute jednačine (8-25). Imajući u vidu njihov karakter, kao i slučajan ili udarni karakter pobude mikro- profila puta, evidentno je da ih moramo rešavati numerički, meto-dom Runge-Kuta. Numerička integracija izvršena je sa vremenskim kora-kom od 0,01, s u 4096 tačaka, što je dovelo do pouzdanosti rezultata o oblasti, približno 0,025-50 Hz (Bendat and Piersol, 2000), (Bendat and Piersol, 1980), (Bendat, 1998), a to je zadovoljavajuće sa aspekta oscila-torne udobnosti i ponašanja vozila na putu (Simić, 1988).

Na osnovu slike 1. može se videti da posmatrani model vozila osci-luje pod uticajem pobuda mikroneravnina puta na prednjim i zadnjim toč-kovima.

U literaturi postoji veliki broj podataka o neravninama mikroprofila puta (Barbarić, 1996, pp.699-710), (Demić, 1997a), ali su oni, pretežno, za-snovani na definisanju spektara snage. Ti spektri ne omogućavaju vršenje inverzne Furijeve transformacije, a koja je neophodna pri definisanju vre-menskih funkcija pobuđivanja (Demić, 1997a). Zbog toga je u ovom radu prihvaćen postupak definisanja poliharmonijske vremenske funkcije pobude iz (Demić, 1997a), a koji će ukratko biti opisan. Polazni podaci o mikro-neravninama puta preuzeti su iz (Demić, 1997a), a prikazani na slici 3.

Pretpostavićemo da funkciju pobude od mikroneravnina puta na le-vom prednjem točku vozila opisuje relacija:

NH

zopl(t) = z ) sm\_2nft + s{ f)] (26)

i

gde su:

• zopl - funkcija pobude od mikroneravnina na prednjem levom točku,

• A(f) - amplitude mikroneravnina puta,

• f - učestanost,

• 9(f) - fazni ugao,

• t - vreme,

• NH - željeni broj harmonika.

U ovom radu usvojena je funkcija pobude od 100 harmonika, koja omogućava pobuđivanje u oblasti učestanosti 0,015-35, Hz.

Na osnovu slike 3, možemo napisati:

A( f) pl = (A + B0 fp )/ v (27)

gde su:

• Ao, Bo - koeficijenti zavisni od vrste putnog pokrivača,

• l - talasna dužina.

• fp - talasna dužina neravnine (fp = 1/l),

• v - brzina vozila.

U nedostatku podataka za fazne uglove pretpostavićemo ih u obliku (Demić, 1997a):

s(f)pi = 2n(RND - 0.5) (28)

gde su RND ravnomerno raspoređeni slučajni brojevi u intervalu {0, 1}.

Funkcija pobude levog zadnjeg točka može se izračunati na osnovu izraza:

Zozl = (t - L / v) (29)

gde je L = a+b - osno rastojanje vozila.

Slika 3 - Granične vrednosti puteva (1- savremeni put, 2 - loš betonski put, 3 - reparirani asfaltni put, 4 - gradska ulica u dobrom stanju)

Figure 3 - Limit values of road roughness (1- for modern roads, 2 - for a road with a bad concrete cover, 3 - for a repaired asphalt road and 4 - for a city street in good condition)

Demić, M., Prilog istraživanju toplotnog opterećenja amortizera poluaktivnog sistema oslanjanja motornog vozila, pp. 26-50

VOJNOTEHNIČKI GLASNIK/MILITARY TECHNICAL COURIER, 2013., Vol. LXI, No. 3

Imajući u vidu slučajan karakter izraza (26), kao i (27-29), radi upro-šćenja usvojeno je da se i pobuda na desnoj strani vozila može opisati iz-razima (26-29).

Parametri analiziranog vozila prikazani su u tabeli 1 (Zastava infor-macije, 1977-2010).

Tabela 1 - Osnovni podaci o parametrima vozila Table 1 - Basic data on vehicle parameters

Osloniena masa, ms, kg 1500 **

Masa prednjeg točka, mB, kg 59

Masa zadnjeg točka, mz, kg 59

Moment inercije, Ix, kgm2 750 **

Moment inercije, Iy, kgm2 2160 **

Moment inercije, Iz, kgm2 2160 **

Osno rastojanje, L, m 3,1

Položaj težišta, a/b, - 1,4/1,7

*Krutost pneumatika uz linearni član, cp1,N/m 190000

*Krutost pneumatika uz kvadratni član, cp2,N/m2 1900000

*Krutost pneumatika uz linearni član, cp3, N/m3 19000000

Željena visina vozila, zžej, m 0,5

Željeno valjanje oslonjene mase, ф žej, rad 0

Željeno galopiranje oslonjene mase, 0 žej, rad 0

Karakteristična brzina vozila, v, m/s 30

* Za pneumatike prednjih i zadnjih točkova

** Inercijalni parametri za neopterećeno vozilo pomnoženi su faktorom 0,75

Termičko opterećenje amortizera poluaktivnog sistema oslanjanja

Usled relativnog kretanja oslonjene i neoslonjenih masa vrši se me-hanički rad u amortizeru koji je ekvivalentan količini toplote Q, J, (Atkins, 2010), (Mitschke, 1997), (Reisenburg, 1970). Mehanički rad (količinu toplote) definiše izraz (Mitschke, 1997):

^ T

Ai = j FaAr,i = j Fa,ivr,idt (30)

0 0

gde su:

• Fai - sila u amortizeru,

• sri - relativni put kretanja klipa u odnosu na telo amortizera (izraz 13),

• vri - relativna brzina kretanja klipa u odnosu na telo amortizera (izraz 14),

C4D

• t - vreme, a

• i - pl,pd,zl,zd.

Snaga koja je potrebna za izazivanje relativnog kretanja klipa amor-tizera u odnosu na telo data je izrazom:

i - pl,pd,zl,zd

Kao što je poznato (Atkins, 2010), (Bojić, 2011), (Ilić i dr, 1996), (Ma-lić, 1972), (Milinčić i dr, 1984), (Moran et al, 2010), (Reisenburg, 1970), (Šelmić, 1986), ova snaga identična je toplotnom protoku (fluksu) prosti-ranja toplote konvekcijom, a njena srednja vrednost je:

i - pl,pd,zl,zd

Odvođenje najvećeg dela količine toplote na okolni vazduh vrši se konvekcijom (Mitschke, 1997). Formalna zavisnost data je izrazom (Reisenburg, 1970):

gde su:

• a - koeficijent prelaza toplote,

• S - površina sa koje se toplota odvodi,

• At - razlika je temperatura vazduha i spoljašnjosti amortizera,

• i - pl, pd, zl, zd.

Kao što je već napomenuto, u ovom radu nisu vršene analize odvo-đenja toplote sa amortizera, jer nisu poznate veličine ai i Si, a očigledno je da su neophodna veoma opsežna eksperimentalna istraživanja radi od-ređivanja pomenutih veličina, što će svakako biti predmet posebne pa-žnje u narednom periodu.

Uz korišćenje programa razvijenog u Paskalu, izračunati su parame-tri toplotnog opterećenja svakog pojedinačnog amortizera. Preciznije re-čeno, za svaki amortizer izračunati su snaga (odgovara toplotnom fluksu) i rad (odgovara količini toplote). Ocenjeno je celishodnim da se u analizu

P(t). = F v .

\ /г a,г r,г

(31)

(32)

P = aSAt

(33)

Metod i analiza podataka

Demić, M., Prilog istraživanju toplotnog opterećenja amortizera poluaktivnog sistema oslanjanja motornog vozila, pp. 26-50

VOJNOTEHNIČKI GLASNIK/MILITARY TECHNICAL COURIER, 2013., Vol. LXI, No. 3

uključe brzina kretanja vozila (30 i 5, m/s), kao i opterećenost vozila (100%, 75%), a pri simulaciji su korišćeni izrazi (9-33).

Radi ilustracije, toplotni fluks za brzinu kretanja vozila od 30 m/s i potpuno opterećeno vozilo prikazan je na slici 4. Analizom podataka o to-plotnom fluksu, čiji je ilustrativan primer prikazan na slici 4, može se uoči-ti da termički procesi u amortizerima poluaktivnog sistema za oslanjanje, pri kretanju vozila po putu sa stohastičkim parametrima mikroprofila puta, predstavljaju stohastičke procese. Pri tome, brzina kretanja vozila i opte-rećenje vozila utiču na termički fluks. Takođe, postoji i uticaj položaja amortizera u odnosu na težište vozila, što je i razumljivo ako se imaju u vidu oscilatorni procesi vozila i njegovih sistema, izraz (14).

Q

N

Ш

_r

Q

o"

CL

6

_T

CL

m

(Л

3

0,000060

0,000055

0,000050

0,000045

0,000040

0,000035

0,000030

0,000025

0,000020

0,000015

0,000010

0,000005

0,000000

-0,000005

B

C

D

E

t---1----1---1---1---1----1---1---r

0 5 10 15 20

Vreme, s

Slika 4 - Toplotni fluks za brzinu vozila od 30 m/s i potpuno opterećeno vozilo Figure 4 - Heat flux for a speed of 30, m/s and a fully loaded vehicle

Kako je termički fluks slučajna veličina, ocenjeno je celishodnim da se izračunaju srednje snage - fluksevi (Bendat and Piersol, 1980), (Ben-dat and Piersol, 2000), (Bendat, 1998) za svaki amortizer posebno, a po-daci su dati u tabeli 2.

Analizom podataka iz tabele 2. može se utvrditi da termički fluksevi svih amortizera opadaju sa porastom brzine kretanja vozila, a što je lo-gično imajući u vidu da su funkcije pobude od mikroprofila puta manje pri većim brzinama (izraz 27).

Takođe, iz tabele 2 može se videti da toplotni fluksevi opadaju sa smanjenjem opterećenja vozila (inercijalnih parametara).

C4D

Tabela 2 - Srednji toplotni fluks za vreme kretanja od 20 s Table 2 - Medium heat flux for a movement of 20,s

Srednji fluks, W Prednji levi amortizer Prednji desni amortizer Zadnji levi amortizer Zadnji desni amortizer

Brzina 30 m/s, masa 100% 2.330434E-006 5.569178E-006 5.358273E-006 2.944404E-006

Brzina 5 m/s masa 100% 2.435122E-006 5.673833E-006 5.477176E-006 2.968751 E-006

Masa 75%, brzina 5 m/s 2.431798E-006 4.518553E-006 4.263119E-006 2.748116E-006

Radi ilustracije, na slici 5 prikazan je rad koji se pretvara u količinu toplote u amortizerima poluaktivnog sistema za oslanjanje, za brzinu kre-tanja vozila od 5 m/s.

Q

ч

Ш

_r

4

Q

o'

CL

6

CL

m

T3

&

0,00013-|

0,000120,00011 -0,000100,00009 -0,00008 -0,00007 -0,00006 -0,00005 -0,00004 -0,00003 -0,00002 -0,00001 -0,00000 -

-0,00001 A------1 1 1 1 1 1 1 1 10----------------------------------------------------------------5-10-15-20

Vreme, s

Slika 5- Mehanički rad (količina toplote) u amortizerima pri brzini kretanja vozila od 5 m/s pri

potpuno opterećenom vozilu

Figure 5- Mechanical work (the heat amount) in shock absorbers for a speed of 5, m/s and a fully

loaded vehicle

Analizom podataka o količini toplote za brzine kretanja od 5 i 30 m/s, čiji je ilustrativan primer za brzinu od 5 m/s prikazan na slici 5, može se utvrditi da produkcija količine toplote raste sa pređenim putem (vreme-nom) i da je veća za manje brzine kretanja vozila, što je u skladu sa tu-mačenjem toplotnog fluksa, čiji je to integral.

Demić, M., Prilog istraživanju toplotnog opterećenja amortizera poluaktivnog sistema oslanjanja motornog vozila, pp. 26-50

VOJNOTEHNIČKI GLASNIK/MILITARY TECHNICAL COURIER, 2013., Vol. LXI, No. 3

Da bi se utvrdio uticaj opterećenja vozila na količinu toplote koja se oslobađa u amortizerima poluaktivnog sistema vozila, posmatraće se ilu-strativan primer sa slike 6. Naime, na toj slici prikazan je mehanički rad (pretvara se u toplotu) koji se vrši u amortizerima delimično opterećenog vozila (75%) pri kretanju brzinom od 5 m/s.

Slika 6- Mehanički rad (količina toplote) u amortizerima delimično opterećenog vozila pri njegovom kretanju brzinom od 5, m/s

Figure 6- Mechanical work (heat amount) in shock absorbers for a partially loaded vehicle and a

speed of 5, m/s

Na osnovu podataka sa slika 5. i 6. može se utvrditi da smanjeno opterećenje vozila dovodi do manjeg termičkog opterećenja amortizera vozila, što je u skladu sa toplotnim fluksom.

Zaključak

Na osnovu izvršenih istraživanja može se zaključiti:

1. Korišćenjem modela vozila sa poluaktivnim sistemom za oslanja-nje može se izračunati termičko opterećenje njegovih amortizera,

2. Sa porastom brzine kretanja vozila pri njegovom kretanju po putu istovetnih karakteristika opada produkcija toplote u amortizerima.

3. Sa porastom opterećenja vozila raste produkcija toplote u amorti-zerima.

Literatura

Atkins, P. 2010. The Laws of Thermodynamics (A Very Short Introduction). Oxford: Oxford University Press.

Barbarić, Ž. 1996. Uprošćeni matematički model za formiranje slike terena optoelektronskim senzorima. Vojnotehnički glasnik, Vol. 44. No. 6, str. 699-710.

Belingardi, G., & Demić, M. 2009. A possible model for shock absorber by using the "Black box" method. Istraživanja iprojektovanja za privredu, 7(4), str. 45-53.

Bendat, J., & Piersol, A. 1980. Engineering applications of correlation and spectral analysis.New York: John Wiley and Sons.

Bendat, J. 1998. Nonlinear Systems - Techniques and Applications.Lon-don: John Wiley and Sons.

Bendat, J., & Piersol, A. 2000. Random data analysis and measure-ment.London: John Wiley and Sons.

Bojić, M. 2011. Termodinamika.Kragujevac: Mašinski fakultet u Kragujevcu.

Demić, M. 1994. Optimization of Vehicles Elasto-Damping Elements Characteristics from the Aspect of Ride Comfort. Vehicle System Dynamics, 23(1), pp. 351-377. doi:10.1080/00423119408969066

Demić, M. 1996. Optimization of Characteristics of Elasto-Damping Elements from the Aspect of Oscillatory Comfort and Vehicle Handling. Int. J. Vehicle Design, 17(1), pp. 22-46.

Demić, M. 1997a. Optimizacija oscilatornih sistema motornih vozila. Kragu-jevac: Mašinski fakultet u Kragujevcu.

Demić, M. 1997b. Prilog utvrđivanju pouzdanosti nove metode za identifikaciju oscilatornih parametara motornih vozila. Vojnotehnički glasnik, Vol. 45, No. 2. str. 135-146.

Demić, M. 1999. The Definition of the Tires Limit of Admissible Nonuniformity by Using the Vehicle Vibratory Model. Vehicle System Dynamics, 31(3), pp. 183-211. doi:10.1076/vesd.31.3.183.2012

Demić, M. 2004. Projektovanje putničkih automobila. Kragujevac: Mašinski fakultet u Kragujevcu.

Demić, M., Diligenski, Đ., Demić, I., & Demić, M. 2006. A method of vehicle active suspension design. Forsch. Ingenieurwes, 70(5), pp. 145-158. doi:10.1007/s10010-06-0025

Demić, M., & Belingardi, G. 2010. A Contribution to shock absorber modelling and Analysis of their Influence on Vehicle ride Charakteristics. Journal of Middle European Construction and Design of Cars (MECCA), 01, pp. 6-17.

Demić, M., & Diligenski, Đ. 2012. A Contribution to Research of Degradation of Characteristics of Vibration Parameters on Vibration Aspekt of Vehicle Comfort. Journal of Applied Engineering Science, 10(4), pp. 185-190.

Doule, J., & and other, 1990. Feedback Control Theory.Macmillan Publishing Co.

Dorf, R., & Bishop, R. 2004. Modern Control Systems.Prentice Hall.

Fermi, E. 2011. Thermodynamics.Dover Books on Physics.

Demić, M., Prilog istraživanju toplotnog opterećenja amortizera poluaktivnog sistema oslanjanja motornog vozila, pp. 26-50

VOJNOTEHNIČKI GLASNIK/MILITARY TECHNICAL COURIER, 2013., Vol. LXI, No. 3

Florin, M., & and other, 2004. U: CONAT2004, Brashov, pp. 20-22 CO-NAT200418.

Genta, A. 2003. Motor Vehicle Dynamics: Modeling and Simulation.World Scientific Publishing.

Gillespie, T. 1992. Fundamentals of Vehicle Dynamics.SAE International.

Hac, A., & and other, 1999. Elimination of Limit Cycles Due to Signal Estimation in Semi-active Suspensions. U: Steering and Suspension Technology Symposium 1999, Detroit, Michigan, pp. 1-7

Hrovat, D., & Hubbard, M. 1982. Optimum Vehicle Suspensions Minimizing RMS Ratllespace, Sprung mass Acceleration and Jerk. Journal of Dynamic Systems, Measurement, and Control, 103(3), pp. 228-236.

Ilić, G., & i dr., 1996. Termodinamika II - osnove prostiranja toplote.Niš: Ma-šinski fakultet u Nišu.

Malić, D. 1972. Termodinamika i termotehnika.Beograd: Građevinska knjiga.

Margolis, D.L. 1983. Semi-Active Control of Wheel Hop in Ground Vehicles. Vehicle System Dynamics, 12(6), pp. 317-330. doi:10.1080/00423118308968760

Marić, M. 2002. Nauka o toploti: Termodinamika, prenos toplote, sagoreva-nje.Novi Sad: Fakultet tehničkih nauka.

Merrit, H. 1967. Hydraulic Control Systems.New York, London, Sydney: John Willey & sons, Inc.

Miliken, W.F., & Miliken, D.L. 1995. Race Car Vehicle Dynamics.Warrenda-le: SAE.

Milinčić, D., & i dr., 1984. Problemi iz prostiranja toplote.Beograd: Građe-vinska knjiga.

Milinčić, D., & Voronjec, D. 1991. Termodinamika.Beograd: Mašinski fakultet u Beogradu.

Mitschke, M. 1997. Kraftfahrzeugkonstruktion, Teil D.Forlesung: TU Braunschweig.

Mohd Sh., at all., 2008. Experimental heat transfer study on the shock absorber operation. U: International Conference on Science & Technology: Aplica-tions in Industry & Education, Penang Malaysia, pp. 759-765

Moran, M.J., & at all., 2010. Fundamentals of Engineering Thermodyna-mics.John Wiley & Sons.

NEWEUL, Manual 2000. TU Stutgart.

Nell, S., & Steyn, J.L. 1998. An alternative control strategy for semi-active dampers on off road vehicles. Journal of Terramechanics, 35(1), pp. 25-40. doi:10.1016/S0022-4898(98)00010-X

Popović, V. 2001. Projektovanje i simulacija sistema aktivnog oslanja-nja.Beograd: Mašinski fakultet u Beogradu.

Rajamani, R. 2006. Vehicle Dynamics and Control.Springer Verlag Gmbh.

Riesenburg, O. 1970. Beitrag zur Klarung der Vorgange in einem hydrau-lischen Schwingungsdampfer.TU Braunschweig.

C4D

Schlelen, W. 1986. Modeling and Analysis of Nonlinear Multibody Systems. Vehicle System Dynamics, 15(5), pp. 271-288. doi:10.1080/00423118608968856

Silani, E. 2004. Active and semiactive suspensions control strategies in road vehicles.Politecniko di Milano, Dipartimento di electronica e informacione.

Simić, D. 1988. Dinamika motornih vozila.Beograd: Naučna knjiga.

Slaski, G., & Walerjanczyk, W. Possibilities of impruving active safety byusing semi-active suspension. U: KONMOT, Krakov. pp. 597-604

Šelmić, P. 1986. Tehnička termodinamika. Beograd: Naučna knjiga.

Tomović, R., & i dr., 1974. Uvod u nelinearne sisteme automatskog upra-vljanja.Beograd: Naučna knjiga.

Yoshimura, T., & Watanabe, K. 2003. Active suspension of a full car model using fuzzy reasoning based on single input rule modules with dynamic absorbers. Int. J. of Vehicle Design, 31(1), pp. 22-40.

Zastava informacije 1977-2010

Wolftram Research: Mathematica 5.2., Preuzeto sa http://www.wol-fram.com

A CONTRIBUTION TO THE INVESTIGATION OF THE HEAT LOAD OF SHOCK ABSORBERS OF SEMI-ACTIVE SUSPENSIONS IN MOTOR VEHICLES

FIELD: Motor Vehicles

ARTICLE TYPE: Original Scientific Paper

Summary:

Dynamic simulation, based on modeling, has a significant role during the process of vehicle development. It is especially important in the first stages of vehicle design, when relevant vehicle parameters are to be defined. Shock absorbers as executive parts of vehicle semiactive suspension systems suffer thermal loads, which may result in damage and degradation of ther characteristics. Therefore,this paper shows an attempt to analyze converting of mechanical work into heat by using the dynamic simulation method.

Introduction

Shock absorbers are integral elements of semi-active suspension systems for vehicles (hereinafter SASS). They directly affect the active vehicle safety. The role of shock absorbers is to absorb mechanical vibrations transferred from the road and to ensure the safety of passengers in a vehicle.

The kinetic energy of vehicle vibrations transforms into mechanical work or heat in shock absorbers. In practice, in the first stage of vehicle development, the shock absorber parameters are

C47>

Demić, M., Prilog istraživanju toplotnog opterećenja amortizera poluaktivnog sistema oslanjanja motornog vozila, pp. 26-50

VOJNOTEHNIČKI GLASNIK/MILITARY TECHNICAL COURIER, 2013., Vol. LXI, No. 3

chosen from the condition of damping vibrations of vehicles, but their thermal shock loads should be also taken into account.

Motor vehicles have complex dynamic characteristics manifested by spatial movement, parameters change during operation, a number of disturbing influences, backlash, friction, hysteresis, etc. The above-mentioned dynamic phenomena, especially vibration, lead to fatigue of driver and users, reduce the life of the vehicle and its systems, etc. The main objective of the system is to reduce the reliance of the above-mentioned negative effects, improving the vehicle behavior on the road and allow the exploitation of vehicles in a wide range of service conditions. Classical systems cannot satisfiy these conditions, so there was a need to introduce new suspension systems with controlled characteristics (briefly called "semi-active", or "active" systems).

Oscillatory model of vehicle

The differential equations of vibratory motion of the given vehicle model have been written using the "NEWEUL." software package. The data input has been done in accordance with the program requirements, and 24 coordinate systems (including inertial-OI1I2I3 global coordinate system) have been used to define the dynamic vehicle model with seven degrees of freedom.

Thermal load on the shock absorbers of the semi-active suspension system

Due to the relative motion between sprunged and unsprunged masses, mechanical work in the absorbers is equivalent to the amount of heat Q, J. As the heat flux is of a random size, it was assessed as appropriate to calculate the average power (flux) for each buffer separately, and the data are given in Table 1

Table 1 - Medium heat flux during a movement of 20, s

Medium flux, W Front left absorber Front right absorber Rear left absorber Rear right absorber

Speed 30,m/s, Mass 100% 2.330434E-006 5.569178E-006 5.358273E-006 2.944404E-006

Speed 5, m/s Mass 100% 2.435122E-006 5.673833E-006 5.477176E-006 2.968751 E-006

Mass 75%, Speed 5, m/s 2.431798E-006 4.518553E-006 4.263119E-006 2.748116E-006

For illustration, Fig. 1 shows the work converted to heat amount in the shock absorbers of the semi-active suspension system, for a vehicle speed of 5 m/s.

UJ

Q

4-

о

CO

I

re

о

'E

re

£

О

a>

0,00013

0,00012

0,00011

0,00010

0,00009

0,00008

0,00007

0,00006

0,00005

0,00004

0,00003

0,00002

0,00001

0,00000

-0,00001

B

C

D

E

T------'-----1-----1------1-----1-----1-----'------1-----'

0 5 10 15 20

Time, s

Figure 1 - Mechanical work (heat amount) in shock absorbers for a speed of vehicle 5,

m/s and a fully-loaded vehicle

The analysis of the heat amount data for a speed of 5 and 30, m/s,for an illustrative example of a speed of 5, m / s , given in Fig. 1 shows that the amount of heat production increases with the covered distance (time). If the vehicle speed increases, the amount of heat decreases.

The example in Fig. 2 shows the effect of the vehicle load on the heat amount generated in the semi-active suspension system absorbers. Namely, the figure shows mechanical work (converted into heat) in the shock absorbers of a partially loaded vehicle (75%) while moving at a speed of 5, m / s.

Demić, M., Prilog istraživanju toplotnog opterećenja amortizera poluaktivnog sistema oslanjanja motornog vozila, pp. 26-50

VOJNOTEHNIČKI GLASNIK/MILITARY TECHNICAL COURIER, 2013., Vol. LXI, No. 3

uj

“

о

lT

4-

m

I

TO

о

'E

то

■П

о

ф

0,00013

0,00012

0,00011

0,00010

0,00009

0,00008

0,00007

0,00006

0,00005

0,00004

0,00003

0,00002

0,00001

0,00000

-0,00001

B

C

D

E

т------1-----1-----1-----1-----1-----1-----1-----г

0 5 10 15 20

Time, s

Figure 2 - Mechanical work (heat amount) of the shock absorbers in a partially loaded

vehicle at a speed of 5, m/s

It can be determined that the reduced load on a vehicle leads to reduced thermal shock loading of vehicles, which is in accordance with heat flux.

Conclusion

Based on the study,the following can be concluded:

• Models of vehicles with semi-active suspension systems can be used for calculating thermal load on the suspension,

• speed increase during vehicle movement along the path of identical properties, decreases the production of heat in shock absorbers, and

• increase in load of the vehicle leads to the increase in the production of heat in shock absorbers.

Key words: Vehicle, Semi-active suspension system, Shock Absorber, Thermal load of shock absorbers

Datum prijema članka/Paper received on: 13. 03. 2013.

Datum dostavljanja ispravki rukopisa/Manuscript corrections submitted on: 04. 04. 2013.

Datum konačnog prihvatanja članka za objavljivanje/ Paper accepted for publishing on: 06. 04. 2013.