kinematika-jednačine kretanja
Kinematika je dio teorijske mehanike koji proučava kretanje
tijela ne uzimajući u obzir njihovu masu i sile koje dejstvuju
na njih. Kretanje tijela se vrši tokom vremena u prostoru, pa
kinematika ima dvije osnovne veličine: put s
i vrijeme t. Njihove osnovne jedinice su metar i
sekunda. U klasičnoj mehanici vrijeme je pozitivna skalarna
veličina koja se stalno mijenja, pa se uzima za nezavisno promjenljivu.
Sve ostale veličine u kinematici se posmatraju kao funkcije vremena.
Razlikuju se početno vrijeme (početak kretanja) t0
, proteklo (od starta) vrijeme ti i interval
vremena Δt=i+1 - ti-1.
Za definiciju kretanja u kinematici se, osim osnovnih, koriste i
izvedene veličine: ubrzanje a i brzina v.
Ubrzanje je promjena brzine u intervalu vremena a=Δv/Δt.
Brzina je pređeni put u intervalu vrremena v=Δs/Δt.
U kinematici se posmatra kretanje krutih tijela, tj. tijela koja ne
mijenjaju svoj oblik (ne deformišu se). Kretanje nekog tijela je poznato
ako je tokom vremena poznat položaj svake tačke tog tijela. Zbog
toga se kinematika bavi kretanjem i tačke i tijela, pa se ona i
dijeli na: kinematiku tačke i
kinematiku krutog tijela
Položaj tačke u nekom trenutku zadaje se u Dekartovom
sistemu, sa koordinatnim početkom
O
i tri međusobno upravne ose
xyz. Putanja
(trajektorija) tijela je neprekidna linija koju opisuje težište
tijela pri kretanju u odnosu na Dekartov koordinatni sistem.
Težište je zamišljena tačka u kojoj djeluje
težina tijela.
Prema položaju trajektorije razlikuju se
- kretanje u ravni
- kretanje u prostoru
Prema obliku putanje kretanje se dijeli na:
- pravolinijsko kretanje
- krivolinijsko kretanje
Prema brzini kretanje može biti:
-
ravnomjerno ili jednoliko (brzina ne mijenja ni intenzitet
ni pravac)
-
ubrzano ili promjenljivo (brzina mijenja intenzitet i/ili
pravac)
Pravolinijsko kretanje može biti
jednoliko (brzina ne mijenja intenzitet)
i
promjenljivo (brzina mijenja intenzitet, postoji
ubrzanje).
Krivolinijsko kretanje je uvijek
promjenljivo, jer brzina mijenja
svoj pravac u prostoru tokom kretanja tijela.
Pravolinijsko kretanje koje ima konstantno ubrzanje
a
naziva se
jednako ubrzano pravolinijsko kretanje.
Za njega zakon promjene brzine glasi:
v = v0 +
a*t.
Ako je poćetna brzina
v0=0 jednaka nuli,
onda zakon izgleda
v = a*t.
Jednačina zavisnosti puta od brzine za jednako ubrzano pravolinijsko
kretanje je
s = v0*t + a*t2/2.,
a ako je početna brzina jednaka nuli
v0=0
zakon pređenog puta glasi
s = a*t2/2.
Vertikalni hitac naviše je poseban slučaj jednako ubrzanog
pravolinijskog kretanja tijela u polju sile zemljine teže, sa početnom
brzinom
v0 ≥ 0. Ubrzanje vertikalnog
hica jednako je ubrzanju zemljine teže
g = 9,81 m/s i
suprotnog je smjera od brzine, što znači da je smanjuje. Kod
vertikalnog hica naviše ubrzanje uslijed sile gravitacije djeluje
kao
usporenje.
Kosi hitac je ravansko krivolinijsko kretanje tijela u polju
gravitacije sa početnom brzinom v0 ≥ 0.
Tijelo se kreće pod uglom koji se posmatra u ondosu na horizontalni pravac.
Putanja tijela u slučaju kosog hitca je parabola. Odstojanje od koordinatnog
početka do padne tačke je domet XD.
Domet zavisi od početne brzine tijela v0 i od
ugla θ pod kojim hitac startuje(ugao elevacije).
Najviša tačka putanje naziva se tjeme putanje.
dinamika-sile, masa, inercija
Nije teško razlikovati objekat (tijelo) u kretanju i objekat koji miruje.
Nadzvučni avioni, sportski automobili ili biciklo mogu biti objekti u
kretanju ili mirovanju. Što uzrokuje da se jedan objekt kreće ili da
miruje? Što je razlog da se objekat u kretanju zaustavi?
Odgovor je jednostavan. Uzrok je sila (force). Sila je bilo kakvo
dejstvo koje može učiniti da neko tijelo promijeni svoj oblik ili
kretanje.
Dinamika je dio fizike (posebno klasične mehanike) koji se bavi
izučavanjem uticaja sila na kretanje materijalnih tijela. Temelje
dinamike postavio je Sir Isaak Newton definicijom
zakona kretanja. Dinamika, za razliku od kinematike uzima u obzir
tri fizička svojstva materijalnih tijela u dodiru, masu rakete m,
gustinu vazduha ρ kao i masu gasova koja ističe iz raketnog
motora ω.
Tijelo rakete, vazduh i produkti sagorijevanja barutnog goriva djeluju
uzajamno. Produkti sagorijevanja potiskuju raketu dajući joj brzinu,
težina usporava raketu, a vazduh svojom gustinom stvara aerodinamičku
silu otpora koja takođe usporavaju raketu.
Sile koje djeluju na raketu u letu su:
- Sila potiska - T (thrust)
- Sila teže - W (weight)
- Sila otpora - D (drag)
Dejstvo sila na raketu je kompleksan fenomen. Da bi se pojednostavio
proračun putanje, stvara se uprošćen dinamički model kretanja
rakete. Model predpostavlja da:
-
Sila potiska djeluje saosno na raketu u smjeru kretanja,
u centru izlaznog presjeka mlaznika.
-
Sila otpora djeluje saosno na raketu u centru otpora rakete,
suprotnog smjera od smjera kretanja.
-
Sila teže djeluje u težištu rakete, saosno sa raketom,
suprotnog smjera od smjera kretanja.
-
Raketa se kreće vertikalno naviše u polju sile gravitacije
(putanja je vertikalni hitac) , sa početnom brzinom
v0 ≥ 0 .
-
Raketa se lansira sa nivoa mora. Temperatura vazduha je 20°.
Kretanja vazduha (vjetra) nema. .
-
Raketa je dealno uravnotežena. Težište se poklapa sa
koordinatnim početkom rakete. Uzdužna osa spoljnje konture je
identična sa uzdužnom koordinatnom osom rakete. Napadna tačka
sile potiska je tqkođe na uzdužnoj osi rakete. .
Model kretanja rakete ne dovodi do grešaka u proračunu, samo
ga čini jednostavnijim. Dobijeni rezultati parametara putanje
omogućavaju potpunu i realnu ocjenu ostvarenja ciljeva projekta.
Model idealizuje uslove u kojima se odvija kretanje rakete.
Uticaj stvarnih uslova razmatra se u završnoj fazi projekta.
Tada se unose korekcije i dopune u finalnu verziju konstrukcije.
Vertikalni hitac je pogodan za primjenu iz dva razloga.
Taj oblik putanje omogućava lakšu analizu (izračunavanja), a
osim toga takvu putanju imaju modelarske i istraživačke rakete.
Masa rakete je funkcija vremena,
m=f(t), u aktivnoj
fazi leta (dok radi raketni motor-sagorijeva gorivo i ističu produkti
gorenja). Na ostalom dijelu putanje (pasivni dio) ona je nepromjenljiva
m=const.
Promjena sile potiska u vrementu
T=f(t) je poznata
ili se odredi preliminarnim proračunima na osnovu zahtjeva
projekta. Promjena se daje u obliku dijagrama zavisnosti sile
potiska od vremena rada motora. Za modelarski motor totalnog
impulsa
Itot = 5Ns dijagram je dat na
prijethodnoj slici.
Postavljanjem fizičko modela kretanja rakete stvoraju se uslovi
da se proračunaju parametri (veličine) u svakoj tački putanje,
odnosno da se postave dinamičke jednačine kretanja.
Proračunom se određuju maksimalna brzina, promjena brzine tokom leta,
promjena visine tokom leta i maksimalna visina leta
Hmax.
.
primjena jednačina dinamike i proračun putanje
u više približenja od projektnih zahtjeva do konstrukcije rakete
Svaka raketa se pravi sa ciljem da da se ostvare postavljeni projektni
zahtjeve. To mogu biti vrsta i količina korisnog tereta, visina
koja se treba dostići itd. U slučaju modelarske rakete projektni
zahtjevi su pravilan let i usporeno spuštanje na tlo. Modelarska
raketa zbog toga ima ima dobro projektovan stabilizator i padobran
kao korisni teret za bezbjedno spuštanje modela na tlo.
Realizaicja projekta rakete odvija se u nekoliko koraka.
-
Definisanje korisnog tereta je početni korak. Procjenjuju se
vrsta, količina, oblik i dimenzije, količina (težina i
zapremina) itd.
-
Na osnovu poznatog korisnog tereta procjenjuju se
kostrukcija rakete (prečnici, dužine, debljine, raketni
motor, stabilizator itd).
-
Podaci dobijeni u koraku 2 omogućavaju da se ocijeni
preliminarna vrijednost pasivne težine rakete i dobije
poćetni oblik njenih konstrukcionih komponenti. Tako se
dobija mogućnost procjene aerodinamički i balističkih
karakteristika rakete, a sa time i mogućnost proračuna
putanje (maksimalne visine).
-
Korak 4 je traženje granica u kojima treba da se kreće
maksimalna brzina rakete. Zadaju se neke vrijednosti i
uz procjenu dužine aktivnog dijela puta izračunava
odgovarajuća maksmalna visina Hmax,.
Na osnovu dobijenih rezultata procjenjuje se vrijednost
maksimalne brzine rakete, kojom se ostvaruje projektom
zadata visina.
-
Slijedeći potez je da se na osnovu jednačine Ciolkovskog
izračuna količina raketnog goriva potrebna da se ostvari
brzina iz koraka 4. Time su dobijeni podaci potrebni
za prvi potpuni preliminani dinamički proračun putanje.
Na osnovu vrijednosti
Hmax, dobijene
proračunom, ocjenjuje aw koliko preliminarni podaci zadovoljavaju
projektne zahtjeve. Ako je
Hmax
manja od projektom postavljene,
razmatra se povećanje količine raketnog goriva. Razmatraju se i drugi
konstrukcioni parametri. Zatim se sa novim podacima ponovo ulazi u
proraćun maksimalne visine i tako više puta.
Ovaj postupak se naziva
približavanje projektom
zadatim vrijednostima. Kad se za parametre putanje (visinu), dobije
zadovoljavajući rezultat iterativni dinamički proračun je završen.
Zatim se pristupa detaljnoj definiciji konstrukcije rakete i sve
provjerava konačnim proračunom i prikladnim eksperimentima.
faze kretanja rakete na putanji
U kretanju rakete razlikuju se tri faze:
kretanje kroz lanserni uređaj,
aktivni dio puta, i pasivni dio puta.
Pod aktivnim dijelom puta podrazumijeva se dio putanje na kojoj radi raketni
motor (dio puta na kojem sagorijeva raketno gorivo). Jedan dio aktivne
putanje obuhvata i kretanje rakete kroz lansirni uređaj. To je neslobodno,
vođeno kretanje.
Drugi dio aktivne putanje odvija se van lansera. To je slobodno
aktivno kretanje. Na aktivnom dijelu putanje raketa dobija ubrzanje pod
dejstvom sile potiska T raketnog motora.
Raketa mora napuštiti lanser sa brzinom koja je dovoljna za stabilan let
na nevođenom dijelu aktivnog puta. Brzina napuštanja lansera se naziva
početna brzina i označava sa v0.. Vrijednost
početne brzine zavisi od ubrzanja rakete a. Na osnovu
početne brzine određuje se dužina lansera. Preporuka je da raketa ima
ubrzanje a ≥ 15 m/s2 da bi bila stabilna
po izlasku iz lansera.
Na kraju aktivnog puta raketa ima najveću brzinu Vmax
od tog momenta leti po inerciji. To je pasivni dio putanje, zato što nema
uticaja nikakve pogonske sile. Postoje samo aerodinamičke sile i sila
gravitacije. One usporavaju raketu (negativno ubrzanje).
uslovi i predpostavke
sile koje utiču na kretanje
U svakoj fazi leta na raketu djeluju sila potiska T,
sila aerodinamičkog otpora D i
sila gravitacije Q. Sila potiska T
daje raketi ubrzanje, sile D i Q
usporavaju raketu. Predpostavka je da su sve sile saosne osi
simetrije rakete i da nema nikakvih bočnih dejstava, koji bi
remetili let. Radi jednostavnije analize predpostavlja se da je
putanja rakete vertikalni hitac naviše. Kretanje je promjenljivo
pravolinijsko, jer se tokom leta stalno mijenja brzina (ubrzanje i
usporenje). Formula a = T/mR
daje ubrzanje rakete. Sila potiska u toj formuli određuje se iz
dijagrama sile potiska u funkciji vremena za primijenjeni raketni motor.
Za modelarsku raketu to može biti dijagram promjene sile potiska
u funkciji vremena modelarskog raketnog motora B4.
Iz dijagrama sile potiska T = f(t) za
svako zadato vrijeme može se očitati vrijednost sile
T. Masa rakete mR se
mijenja na aktivnom dijelu putanje, zbog sagorijevanja raketnog
goriva. I sila potiska i masa rakete značajno se mijenjaju tokom
vremena. Da bi uticaj tih promjena na vrijednosti proračuna brzine
rakete i pređenog puta bio zanemarljiv, proračuni se vrše za jako
male intervale vremena. Ovakav pristup omogučava da se kretanje
rakete u postavljenom intervalu vremena posmatra kao jednako
ubrzano. Na osnovu jednačina koje važe za pravolinijsko jednako -
ubrzano kretanje mogu se izračunati brzina i put rakete u svakom
intervalu. Njihovim sabiranjem mogu se dobiti parcijalne, ali i
ukupna vrijednost. Ovakav pristup proračunu naziva se numerička
integracija diferencijalnih jednačina.
aktivni dio putanje - kretanje u lanseru
Kretanje rakete počinje u lanseru iz stanja mirovanja . To znači
da je startna brzina vs = 0.
I startno vrijeme je ts = 0.
Lanser je postavljen vertikalno jer se razmatra vertikalni hitac.
Zanemaruju se sile trenja i druga dejstva koja mogu nastati
tokom kretanju kroz lanser i biti remetilački faktor.
Parametri vs i ts
su početne vrijednosti u proračunu kretanja kroz lanser.
Numerička integracija jedančina leta rakete u lanseru obavlja se
za male vremenske intervale, Δt = 0,001sec.
Integracijom se dobija put L, na kojem brzina ima vrijednost dovoljnu
za stabilan let. To je brzina napuštanja lansera,odnosno početna brzina
v0 slobodnog leta rakete.
Put koji odgovara brzini napuštanja lansera, je dužina vodilica
lansitnog uređaja L0. Istovremeno se
dobija i vrijeme prolaska kroz lanser. To vrijeme se
naziva početno vrijeme t0.
slobodna aktivna putanja - kretanje van lansera
Kao što su vs i startno vrijeme
ts polazna tačka za proračun
kretanja kroz lanser tako su i v0
, t0 i L0
polazina tačka za proračun elemenata slobodne aktivne putanje.
Na tom dijelu putanje rakete djeluju iste sile kao u slačaju kretanja
kroz lanser: sila potiska T,
sila aerodinamičkog otpora D i
sila gravitacije Q. Jedino nema vodilica koje
usmjeravaju kretanje rakete. Pretpostavka ja tokom aktivnog leta
nema bočnih udara koji bi djelovali kao remetilački faktor
(udari vjetra itd).Numerička integracija jedančina leta rakete na
ovom dijelu putanje takođe se obavlja za male vremenske intervale,
a preporuka je Δt = 0,01sec.
Na kraju aktivnog dijela putanje sila potiska raketnog motora
je T = 0. Raketa stiče maksimalnu brzinu
vmax. Numeričkom integracijom izračunavaju
se dužina (za vertikalni hitac naviše visina) aktivnog dijela
puta rakete Ha. i vrijeme potrebno da
se on pređe ta..
Maksimalna brzina rakete vmax, dužina
aktivnog dijela puta Ha i njegovo trajanje
ta određuju početnu tačku pasivnog dijela
putanje rakete. Numeričkom integacijom dobijaju se parametri ostalih
tačaka. Proračun se završava maksimalnom visinom leta, a to je
tjeme putanje HT.
U tjemenu je brzina kretanja v = 0,
pa se numerička integracija obavlja samo do te tačke. Pored
tjemena putanje tada se dobija i ukupno vrijeme leta rakete do
tjemena, tT.
Na pasivnom dijelu putanje raketa se kreće na osnovu Prvog Njutnovog
zakona, po inerciji. Nema sile potiska raketnogo
motora, T = 0. Kretanje se odvija samo pod
dejstvom aerodinamičkih sila i sile gravitacije.
To bi bio KRAJ proračuna parametara leta rakete u
izloženom primjeru
>
proračun putanje
jednačine
Kretanje rakete se posmatra kao pravolinijsko, jednako
ubrzano. Za njega kinematske jednačine brzine i pređenog
puta glase: v = a*t i S = v*t/2,
(v-brzina, S-put, a-ubrzanje i t-vrijeme).
Na rješavanje jednačina kretanja rakete primjenjuje se numerička
integracija. Numerička integracija razmatra veoma male intervale
promjena brzine i puta sa vremenom. Dovoljno male da aproksimacija
o pravolinijskom jednako ubrzanom kretanju daje zadovoljavajuće
rezultate. Jednačine kretanja primijenjene na male vremenske
intervale su date na slikama 3 i 4.
ΔSi je promjena puta u proizvoljnom
intervalu Δti , promjena brzine je
Δvi, a Δti
je proizvoljni interval vremena tokom numeričkog integraljenja putanje.
Ako se jednačine primijene na vertikalno, pravolinijsko, promjenljivo kretanje
onda je jednačina promjene puta (visine) ΔHi
kao na slici 5. Kinematske jednačine promjene visine i brzine rakete tokom
vertikalnog leta daju odnose brzine, ubrzanja, puta i vremena. Ako se zna
ubrzanje u nekom intervalu vremena može se izračunati prmjena brzine u
u tom intervalu,
Δvi = a*Δti .
Ako se zna promjena brzine može se izračunati pređena visina u tom intervalu
ΔHi = Δvi*
Δti/2.
Sabiranjem vremenskih intervala Δti
od početka kretanja do nekog proizvoljnog trenutka u toku leta
dobija se vrijeme leta ti do tog
trenutka. Sabiranjem promjena brzina
Δvi
u svim intervalima do nekog proizvoljnog intrvala i,
dobija se ukupna brzina u trenutku leta ti.
Sabiranjem promjena visine
ΔHi
u svim intervalima do nekog proizvoljnog intrvala i,
dobija se ukupna visina u trenutku leta ti.
Međutim, da bi se određivali pređeni put, brzina i vreijeme leta
potrebno je poznavati ubrzanje a kojem je podvrgnuta
raketa.
Prema Drugom Njutnovom zakonu ubrzanje, promjena brzine i
pravca su posledica dejstva spoljašnjih sila na tijelo, u ovom
slučaju raketu.
U jednačini su F - rezultantna sila na tijelo u N,
m - masa tijela u kg i a - ubrzanje u m/s2.
Drugi Njutnov zakon može da se napiše i u drugom obliku (slika 6)
slijedi jednačina (slika 7), u kojoj se ubrzanje prikazuje kao
promjena brzina u intervalu vremena. Suština aproksimacije o
jednako ubrazanom kretanju ovdje je jasno prikazano. Ubrzanje
se mijenja tokom leta, ali u malom intervalu vremena uzima se
njegova srednja vrijednost i tako se može samo u njemu smatrati
konstantnim.
Ovo unosi grešku u proračun. Ukoliko je primijenjeni
interval vremena manji, greška je manja. Praksa je pokazala da se
dovoljna tačnost postiže ako je vremenski inerbal numeričke
integacije 0.001 - 0.01 sec za kretanje u lanseru
i 0.01 - 0.1 za slobodni dio putanje. Sa trenutnim
stanjem računarske tehnike vremenski interval numeričke integracije
može biti i manji, pa se postiže veća tačnost.
Sila F je rezultanta svih sila koje djeluju na raketu
tokom leta:
sile potiska T, sile otpora D i sile gravitacije Q.
Njutnov zakon sada može da se prikaže u razvijenom obliku:
, a iz nje slijedi formula za promjenu brzine u proizvoljnom
intervalu vremena.
Brzina u proizvoljnom trenutku dobija se dodavanjem promjene
brzine u odgovarajućem intervalu.
Brzina vi-1 je brzina iz prethodne
iteracije.