Tilbake
7.3

7.3 Luftmotstand og førsteordens separabel ODE

Drag-modeller, terminalfart, og den separable førsteordens ligningen `v(t) = v_T(1 − e^{−gt/v_T})` — med det som kobler rett videre til Python-sporet.

65 min
8 oppgaver
Luftmotstandførsteordens separabel ODE
Din fremgang i kapitlet
0 / 8 oppgaver

Forkunnskaper — sist du var her

Dette kapitlet bygger på kap. 7.1 og kap. 3.2. Det du trenger derfra, ferdig oppfrisket:

Bevegelseslikningen er Newtons 2. lov med akselerasjonen skrevet som derivert. Når ingen kraft avhenger av posisjonen, kan vi arbeide direkte med farten, og ordenen faller til én:

mdvdt=F(v)m\frac{dv}{dt} = \sum F(v)

Likevektsløsningen finner du ved å sette den deriverte lik null. For et fallende legeme kalles den terminalfarten, og den er en gratis poenghøster: den krever ikke at likningen løses.

Separabel betyr at alt med den ukjente kan samles på den ene siden og alt med tiden på den andre. Alle likningene i dette kapitlet er separable.

Fra matematikken trenger du to ting: Variabelskifte (substitusjon) og Integrasjon av 1/x og eksponentialfunksjoner. Det ene integralet som gjør hele jobben, er

duu=lnu+C\int \frac{du}{u} = \ln\lvert u\rvert + C

Fra kap. 1.2 trenger du dessuten ex1x+12x2e^{-x} \approx 1 - x + \tfrac12x^2 for små xx — den brukes til å vise at løsningen går over i fritt fall like etter slipp.

Notasjonsavtale for kapitlet. Vi regner positiv retning nedover i fallproblemene, slik at tyngden blir positiv og farten positiv når legemet synker. Farten skrives vv som kursiv skalar med fortegn. Luftmotstandens koeffisienter er kvk_v (lineær modell, enhet kg/s\text{kg/s}) og DD (kvadratisk modell, enhet kg/m\text{kg/m}) — de er ulike størrelser med ulike enheter og kan aldri byttes om. Terminalfarten er vTv_T; en periode ville vært TpT_p, og snordrag TT, men ingen av dem opptrer her. Tyngden er mgm\mathbf{g}, aldri G\mathbf{G}. Vi bruker g=9,81 m/s2g = 9{,}81\ \text{m/s}^2.

Løkke 1 — de to dragmodellene (~14 min)

Slipp et kaffefilter fra skulderhøyde. Det når full fart nesten med det samme og daler så jevnt ned. Slipp en stein fra samme høyde, og den akselererer hele veien.

Forskjellen er luftmotstanden, og den har to ansikter i dette faget. Hvilket du får, avhenger av hvor raskt og hvor stort legemet er — men i en eksamensoppgave står det alltid i teksten hvilken modell som gjelder. Din jobb er å kjenne dem igjen, sette dem opp med riktig fortegn, og vite hva de gjør med likningen.

Luftmotstand (drag)

En kraft fra lufta eller væsken på et legeme som beveger seg gjennom den, alltid rettet mot bevegelsen.

Kraften har ingen egen kilde i form av et fast legeme — den kommer av at mediet må skyves til side og av friksjon mot overflaten. Størrelsen vokser med farten, og det er nettopp den fartsavhengigheten som gjør bevegelseslikningen til en ekte differensiallikning.

Merk retningsregelen: «mot bevegelsen» er en betingelse som gjelder til enhver tid, også når legemet snur. Det er derfor modellen må skrives slik at kraften bytter fortegn samtidig med farten.

Lineær dragmodell
Modellen der motstanden er proporsjonal med farten:

FD=kvvF_D = -k_v v

Koeffisienten kvk_v har enhet kg/s\text{kg/s}. Minustegnet gjør at kraften peker mot bevegelsen: er vv positiv, er FDF_D negativ, og omvendt. Modellen håndterer retningsskifte helt av seg selv, siden vv inngår i første potens.

Lineær motstand gjelder for små, langsomme legemer i et seigt medium: en liten kule i olje, et støvkorn i luft, et lodd i glyserol. Den er dessuten den ene av de to modellene som gir en enkel lukket løsning, og derfor den som brukes når oppgaven vil ha et analytisk svar.

Kvadratisk dragmodell
Modellen der motstanden vokser med kvadratet av farten:

FD=DvvF_D = -D\lvert v\rvert v

Koeffisienten DD har enhet kg/m\text{kg/m} og avhenger av tverrsnittet, formen og luftas tetthet. Absoluttverdien er ikke pynt. Uten den ville leddet Dv2-Dv^2 pekt samme vei uansett fartsretning, og et legeme på vei oppover ville blitt drevet videre oppover i stedet for bremset.

Kvadratisk motstand gjelder for større og raskere legemer i luft: en fallskjermhopper, en ball, en sykkel, en bil. Det er den realistiske modellen for det meste av dagliglivets bevegelser, og derfor den som dominerer i de numeriske oppgavene.

Hvorfor absoluttverdien er avgjørende. Tenk på en ball som kastes rett opp og faller ned igjen.

- På vei opp er v>0v > 0 (med positiv retning oppover). Motstanden skal peke nedover. Med Dvv=Dv2-D\lvert v\rvert v = -Dv^2 blir den negativ ✓.
- På vei ned er v<0v < 0. Motstanden skal nå peke oppover. Med Dvv-D\lvert v\rvert v blir leddet Dvv=+Dv2>0-D\lvert v\rvert v = +Dv^2 > 0 ✓, altså oppover.
- Med Dv2-Dv^2 ville leddet vært negativt i begge tilfeller, og på vei ned ville lufta ha dyttet ballen fortere nedover. Modellen ville brutt sammen fysisk.

Intuisjon: v\lvert v\rvert gir størrelsen, og den siste vv-en gir retningen. Produktet har derfor riktig størrelse og riktig fortegn i ett og samme uttrykk.

Dette er felle #4 i bokas feilregister, og den følger deg videre inn i Del 8: i koden skrives leddet -D*abs(v[i])*v[i], aldri -D*v[i]**2. Skriver du kvadratet, taper du poeng både i utledningen og i koden — to steder for samme feil.

Når kan du forenkle? Bare når du kan garantere at farten ikke skifter fortegn i det tidsrommet du regner på. En fallskjermhopper som faller nedover hele veien, kan trygt beskrives med Dv2-Dv^2 — men du må skrive ut betingelsen. Å forenkle uten å si hvorfor er å håpe at sensor antar det beste.

📝Oppgave 1

Forklar med egne ord, i to–tre setninger hver:

a) hvorfor luftmotstand gjør bevegelseslikningen til en differensiallikning som ikke kan løses med konstant-aa-formlene,

b) hvorfor kvadratisk luftmotstand må skrives Dvv-D\lvert v\rvert v og ikke Dv2-Dv^2,

c) hvorfor koeffisientene kvk_v og DD ikke kan ha samme enhet.

Løkke 2 — terminalfarten (~12 min)

Det første spørsmålet i en dragoppgave er nesten alltid det samme: hvor fort går det til slutt? Og det gode ved spørsmålet er at du kan svare uten å løse differensiallikningen i det hele tatt.

Intuisjon: når farten er så stor at motstanden akkurat balanserer tyngden, er resultantkraften null. Da er akselerasjonen null, og farten endrer seg ikke mer. Dette er nøyaktig likevektsløsningen fra kap. 7.1, bare med et fysikknavn på.

Terminalfart
Den konstante farten et fallende legeme nærmer seg når tyngden og luftmotstanden balanserer hverandre:

F=0mg=FD(vT)\sum F = 0 \quad \Longrightarrow \quad mg = \lvert F_D(v_T)\rvert

Terminalfarten er en likevektsløsning av bevegelseslikningen: den konstante funksjonen v(t)=vTv(t) = v_T oppfyller likningen for alle tt.

Legemet når aldri vTv_T i endelig tid — det nærmer seg verdien eksponentielt. I praksis er farten innenfor 5 %5\ \% av vTv_T etter et par tidskonstanter, og da er forskjellen uten betydning.

Terminalfarten er en gratis poenghøster. Den krever bare at du setter akselerasjonen lik null, og den er ofte en egen deloppgave som kan besvares selv om du står fast på resten.

📜Terminalfart for de to modellene
Lineær motstand. mg=kvvTmg = k_v v_T gir

vT=mgkvv_T = \frac{mg}{k_v}

Kvadratisk motstand. mg=DvT2mg = D v_T^2 gir

vT=mgDv_T = \sqrt{\frac{mg}{D}}

Legg merke til skaleringen — den er verdt en kommentar i en besvarelse. Med lineær motstand er terminalfarten proporsjonal med massen; med kvadratisk motstand vokser den bare som kvadratrota av massen. En fallskjermhopper som veier dobbelt så mye, faller derfor bare 2=1,41\sqrt{2} = 1{,}41 ganger så fort — ikke dobbelt så fort.

Benevningskontroll for den kvadratiske: [mg/D]=N/(kg/m)=(kgm/s2)m/kg=m2/s2[mg/D] = \text{N}/(\text{kg/m}) = (\text{kg}\cdot\text{m/s}^2)\cdot\text{m}/\text{kg} = \text{m}^2/\text{s}^2, og kvadratrota gir m/s\text{m/s} ✓.

En tredje sammenheng er verdt å merke seg, fordi den forenkler alle regnestykkene i løkke 3: for lineær motstand er

τmkv=vTg\tau \equiv \frac{m}{k_v} = \frac{v_T}{g}

altså kan tidskonstanten leses rett av terminalfarten.

✏️Eksempel 1: To fallskjermfaser

En fallskjermhopper med utstyr har samlet masse m=80,0 kgm = 80{,}0\ \text{kg}. I fritt fall før utløsning er den kvadratiske dragkoeffisienten D1=0,250 kg/mD_1 = 0{,}250\ \text{kg/m}; med åpen skjerm er den D2=30,0 kg/mD_2 = 30{,}0\ \text{kg/m}.

a) Finn terminalfarten i begge fasene.

b) Hvor stor er akselerasjonen i det øyeblikket skjermen åpner seg, hvis hopperen da har nådd terminalfart i første fase?

c) Kommenter tallet i b).

Metodevalg: terminalfarten er likevektsløsningen, så vi setter F=0\sum F = 0. I b) bruker vi Newtons 2. lov direkte, med den farten hopperen faktisk har i det øyeblikket — ingen differensiallikning trengs.

a) Fritt fall før utløsning:

vT,1=mgD1=80,09,810,250=3139,2=56,0 m/sv_{T,1} = \sqrt{\frac{mg}{D_1}} = \sqrt{\frac{80{,}0 \cdot 9{,}81}{0{,}250}} = \sqrt{3139{,}2} = 56{,}0\ \text{m/s}

nedover, altså 202 km/h202\ \text{km/h}.

Med åpen skjerm:

vT,2=mgD2=784,830,0=26,16=5,11 m/sv_{T,2} = \sqrt{\frac{mg}{D_2}} = \sqrt{\frac{784{,}8}{30{,}0}} = \sqrt{26{,}16} = 5{,}11\ \text{m/s}

nedover, altså 18,4 km/h18{,}4\ \text{km/h} — omtrent farten du får ved å hoppe ned fra halvannen meters høyde, som er nettopp poenget med en fallskjerm.

b) I åpningsøyeblikket er farten fortsatt 56,0 m/s56{,}0\ \text{m/s}, men dragkoeffisienten er blitt D2D_2. Med positiv retning nedover:

ma=mgD2v2=784,830,03139,2=784,894176=93391 Nma = mg - D_2v^2 = 784{,}8 - 30{,}0 \cdot 3139{,}2 = 784{,}8 - 94\,176 = -93\,391\ \text{N}

a=9339180,0=1,17103 m/s2a = \frac{-93\,391}{80{,}0} = -1{,}17 \cdot 10^3\ \text{m/s}^2

altså om lag 1,17 km/s21{,}17\ \text{km/s}^2 rettet oppover (bremsing).

c) Det er 119119 ganger tyngdeakselerasjonen — helt urealistisk for et menneske, som tåler i størrelsesorden 10g10g i kort tid. Konklusjonen er at modellen med en brå overgang mellom to konstante DD-verdier er for grov: i virkeligheten folder skjermen seg ut over et par sekunder, og DD vokser gradvis. En bedre modell lar DD variere med tiden, og den kan bare løses numerisk.

Sensorblikk: en deloppgave som ber om en kommentar til et tall, gir full uttelling for et presist resonnement uten regning. Å regne ut 119g119g og la det stå ukommentert er å la poeng ligge — det er nettopp urimeligheten som er svaret.

📝Oppgave 2

Et kaffefilter med masse 2,50 g2{,}50\ \text{g} faller gjennom lufta. Ved de lave fartene det er snakk om, er motstanden godt beskrevet av den lineære modellen med kv=0,0150 kg/sk_v = 0{,}0150\ \text{kg/s}.

a) Finn terminalfarten.

b) Hva er akselerasjonen i slippøyeblikket?

c) Hva er akselerasjonen når farten er 1,00 m/s1{,}00\ \text{m/s}?

d) To like filtre stables oppi hverandre, slik at massen dobles mens formen og dermed kvk_v er uendret. Hva blir terminalfarten da?

📝Oppgave 3

En hagleskål og en tennisball har samme form og samme dragkoeffisient D=0,00120 kg/mD = 0{,}00120\ \text{kg/m}, men ulik masse. Motstanden er kvadratisk.

a) Finn terminalfarten for en ball med masse 0,145 kg0{,}145\ \text{kg}.

b) Hvor mange ganger så tung må en ball være for å få dobbelt så stor terminalfart?

c) En medstudent sier: «Da faller en dobbelt så tung ball dobbelt så fort.» Hva er galt?

d) Ved hvilken fart er luftmotstanden akkurat halvparten av tyngden for ballen i a)?

Løkke 3 — den separable løsningen (~20 min)

Nå skal vi finne hele forløpet, ikke bare sluttfarten. Vi tar den lineære modellen, som er den ene som gir et pent svar, og vi gjør det steg for steg — for det er utregningen sensor følger.

Bevegelseslikningen for et legeme som slippes fra ro, med positiv retning nedover:

mdvdt=mgkvv,v(0)=0m\frac{dv}{dt} = mg - k_v v, \qquad v(0) = 0

Klassifiseringen først, som alltid: første orden, lineær, inhomogen — og separabel. Det er separabiliteten som gir metoden.

📜Slik løses den separable likningen
Steg 1 — del på massen og faktoriser. Med τ=m/kv\tau = m/k_v:

dvdt=gvτ=1τ(gτv)\frac{dv}{dt} = g - \frac{v}{\tau}= \frac{1}{\tau}\left(g\tau - v\right)

og siden gτ=mg/kv=vTg\tau = mg/k_v = v_T:

dvdt=vTvτ\frac{dv}{dt} = \frac{v_T - v}{\tau}

Steg 2 — separer. Alt med vv til venstre, alt med tt til høyre:

dvvTv=dtτ\frac{dv}{v_T - v} = \frac{dt}{\tau}

Steg 3 — integrer begge sider fra start til tidspunktet tt, med v(0)=0v(0) = 0:

0vdvvTv=0tdtτln(vTvvT)=tτ\int_0^{v}\frac{dv'}{v_T - v'} = \int_0^t \frac{dt'}{\tau} \quad \Longrightarrow \quad -\ln\left(\frac{v_T - v}{v_T}\right) = \frac{t}{\tau}

Steg 4 — løs for vv. Ta eksponentialfunksjonen på begge sider:

vTvvT=et/τ v(t)=vT(1et/τ) \frac{v_T - v}{v_T} = e^{-t/\tau} \quad \Longrightarrow \quad \boxed{\ v(t) = v_T\left(1 - e^{-t/\tau}\right)\ }

og siden τ=vT/g\tau = v_T/g, kan det også skrives

v(t)=vT(1egt/vT)v(t) = v_T\left(1 - e^{-gt/v_T}\right)

Intuisjon steg for steg: minustegnet i steg 3 kommer av at vv står med negativt fortegn i nevneren — det er stedet flest mister et fortegn. Og legg merke til at integrasjonsgrensene bærer initialbetingelsen: skriver du dem inn med det samme, slipper du en ubestemt konstant helt.

Kontroll ved innsetting. Deriver løsningen:

dvdt=vTτet/τ=vTvτ\frac{dv}{dt} = \frac{v_T}{\tau}e^{-t/\tau} = \frac{v_T - v}{\tau} \quad ✓

siden vTv=vTet/τv_T - v = v_Te^{-t/\tau}. Likningen er oppfylt for alle tt, og v(0)=vT(11)=0v(0) = v_T(1-1) = 0 ✓.

Tidskonstant
Størrelsen τ=m/kv\tau = m/k_v, med enhet sekund, som setter tidsskalaen for hvor raskt farten nærmer seg terminalfarten.

v(τ)=vT(1e1)=0,632vTv(\tau) = v_T(1 - e^{-1}) = 0{,}632\,v_T

Etter én tidskonstant er altså 63 %63\ \% av terminalfarten nådd, etter to 86 %86\ \%, etter tre 95 %95\ \% og etter fem over 99 %99\ \%. Tommelfingerregelen «tre tidskonstanter, så er du praktisk talt framme» er verdt å ha.

For den lineære modellen gjelder dessuten τ=vT/g\tau = v_T/g: kjenner du terminalfarten, kjenner du tidsskalaen. En tung, langsom gjenstand har lang τ\tau; et lite, lett legeme når terminalfarten nesten umiddelbart.

✏️Eksempel 2: En isoporkule i luft

En liten isoporkule med masse m=12,0 gm = 12{,}0\ \text{g} slippes fra ro. Luftmotstanden er lineær med kv=0,0450 kg/sk_v = 0{,}0450\ \text{kg/s}.

a) Finn terminalfarten og tidskonstanten.

b) Skriv v(t)v(t).

c) Hvor lang tid tar det før kula har nådd 95 %95\ \% av terminalfarten?

d) Hvor langt har kula falt etter 1,00 s1{,}00\ \text{s}, og hvor langt ville den falt i fritt fall?

Metodevalg: vi bruker den separable løsningen som ble utledet over. Terminalfarten og tidskonstanten regnes ut først, siden begge inngår i uttrykket for v(t)v(t). Positiv retning er nedover.

a) vT=mgkv=0,01209,810,0450=0,117720,0450=2,62 m/snedoverv_T = \frac{mg}{k_v} = \frac{0{,}0120 \cdot 9{,}81}{0{,}0450} = \frac{0{,}11772}{0{,}0450} = 2{,}62\ \text{m/s} \quad \text{nedover}

τ=mkv=0,01200,0450=0,267 s\tau = \frac{m}{k_v} = \frac{0{,}0120}{0{,}0450} = 0{,}267\ \text{s}

Kontroll: τ=vT/g=2,616/9,81=0,267 s\tau = v_T/g = 2{,}616/9{,}81 = 0{,}267\ \text{s} ✓ — de to uttrykkene må gi samme tall.

b) v(t)=2,62(1et/0,267) m/snedoverv(t) = 2{,}62\left(1 - e^{-t/0{,}267}\right)\ \text{m/s} \quad \text{nedover}

c) Vi krever v=0,95vTv = 0{,}95\,v_T, altså et/τ=0,05e^{-t/\tau} = 0{,}05:

t=τln(0,05)=0,26672,996=0,799 st = -\tau\ln(0{,}05) = 0{,}2667 \cdot 2{,}996 = 0{,}799\ \text{s}

altså om lag tre tidskonstanter, som tommelfingerregelen sier.

d) Strekningen er integralet av farten:

y(t)=0tv(t)dt=vT[tτ(1et/τ)]y(t) = \int_0^t v(t')\,dt' = v_T\left[t - \tau\left(1 - e^{-t/\tau}\right)\right]

y(1,00)=2,616[1,000,2667(1e3,75)]=2,616[1,000,2605]=1,93 my(1{,}00) = 2{,}616\left[1{,}00 - 0{,}2667\left(1 - e^{-3{,}75}\right)\right] = 2{,}616\left[1{,}00 - 0{,}2605\right] = 1{,}93\ \text{m}

I fritt fall ville kula falt 12gt2=4,91 m\tfrac12gt^2 = 4{,}91\ \text{m} på samme tid. Luftmotstanden har altså mer enn halvert fallstrekningen på ett sekund.

Benevningskontroll: vTτv_T\tau har enhet (m/s)s=m(\text{m/s})\cdot\text{s} = \text{m} ✓, som seg hør og bør for en strekning.

Rimelighetskontroll: en isoporkule som daler ned med under 3 m/s3\ \text{m/s} etter et kvart sekund, stemmer godt med det du ser når du slipper en. En stein av samme størrelse ville hatt hundre ganger større masse og dermed hundre ganger større terminalfart — og fortsatt akselerert.

📝Oppgave 4

Et lodd med masse 35,0 g35{,}0\ \text{g} slippes fra ro ned i en beholder med olje. Motstanden er lineær med kv=0,0700 kg/sk_v = 0{,}0700\ \text{kg/s}.

a) Finn terminalfarten og tidskonstanten.

b) Skriv v(t)v(t) og finn farten etter 0,500 s0{,}500\ \text{s}.

c) Når er farten halvparten av terminalfarten?

d) Vis at løsningen din oppfyller differensiallikningen.

📝Oppgave 5

En liten kule slippes fra ro i en væske. Farten måles til 2,00 m/s2{,}00\ \text{m/s} etter 0,300 s0{,}300\ \text{s} og til 3,00 m/s3{,}00\ \text{m/s} etter 0,800 s0{,}800\ \text{s}. Motstanden antas lineær, og kula slippes fra ro.

a) Sett opp de to likningene målingene gir.

b) Bestem terminalfarten og tidskonstanten.

c) Hva er kulas masse hvis kv=0,0250 kg/sk_v = 0{,}0250\ \text{kg/s}?

d) Hvorfor er det ikke nok med én måling for å bestemme både vTv_T og τ\tau?

Løkke 4 — grensene, og der den lukkede løsningen ryker (~15 min)

En løsning er ikke ferdig kontrollert før du har sett på grensene. For v(t)=vT(1et/τ)v(t) = v_T(1-e^{-t/\tau}) er det to av dem, og begge er faste deloppgaver på eksamen.

Grensen tt \to \infty. Da går et/τ0e^{-t/\tau} \to 0, og

vvTv \to v_T

Løsningen nærmer seg terminalfarten, som den skal — og det bekrefter at likevektsløsningen vi fant i løkke 2, faktisk er det bevegelsen ender i.

Grensen t0t \to 0. Her bruker vi Taylor-utviklingen fra kap. 1.2, ex1x+12x2e^{-x} \approx 1 - x + \tfrac12x^2:

v(t)=vT(1et/τ)vT(tτt22τ2)=gtg2t22vTv(t) = v_T\left(1 - e^{-t/\tau}\right) \approx v_T\left(\frac{t}{\tau} - \frac{t^2}{2\tau^2}\right) = gt - \frac{g^2t^2}{2v_T}

Det første leddet er nøyaktig fritt fall: like etter slipp er farten null, motstanden derfor null, og legemet faller som om lufta ikke fantes. Det andre leddet er den første korreksjonen, og den er negativ — motstanden begynner å bite.

Intuisjon: modellen må gi fritt fall i starten og terminalfart til slutt. Gjør den ikke det, er noe galt. Dette er den billigste kvalitetskontrollen som finnes på et slikt uttrykk, og på eksamen er den ofte en egen deloppgave med egne poeng.

Hva med den kvadratiske modellen? Her må vi være presise, for her sitter en utbredt misforståelse.

For et legeme som slippes fra ro og faller nedover hele veien, er vv=v2\lvert v\rvert v = v^2, og likningen

mdvdt=mgDv2m\frac{dv}{dt} = mg - Dv^2

er fortsatt separabel. Den kan faktisk løses i lukket form:

v(t)=vTtanh(gtvT),vT=mgDv(t) = v_T\tanh\left(\frac{gt}{v_T}\right), \qquad v_T = \sqrt{\frac{mg}{D}}

Du kan verifisere den ved innsetting akkurat som før. Så påstanden «kvadratisk motstand kan aldri løses analytisk» er feil, og den skal du ikke skrive.

Det som derimot er sant, er hvor smalt gyldighetsområdet er. Løsningen over holder bare så lenge

1. farten ikke skifter fortegn — så snart legemet snur, må vv\lvert v\rvert v skrives ut, og likningen skifter form ved vendepunktet;
2. bevegelsen er i én dimensjon — i to dimensjoner kobles xx- og yy-likningene gjennom farten v=vx2+vy2v = \sqrt{v_x^2+v_y^2}, og systemet har ingen løsning i elementære funksjoner;
3. koeffisientene er konstante — varierer DD med høyden, som i tynnere luft, ryker den også.

Et skrått kast med luftmotstand bryter både punkt 1 og punkt 2, og det er nettopp det som gjør Del 8 nødvendig. Kastet er den vanligste innpakningen av den numeriske deloppgaven i eksamenssettene, og formelarket hjelper deg ikke der: den ferdigheten må sitte i fingrene.

📝Oppgave 6

Ta utgangspunkt i v(t)=vT(1egt/vT)v(t) = v_T\left(1 - e^{-gt/v_T}\right) for et legeme som slippes fra ro med lineær luftmotstand.

a) Vis at uttrykket gir fritt fall for små tt, ved å Taylor-utvikle til og med andreordens ledd.

b) Vis at vvTv \to v_T når tt \to \infty, og forklar fysisk hva det betyr.

c) For hvilken fart er akselerasjonen halvparten av gg?

d) En medstudent har regnet ut v(t)=vT(1+egt/vT)v(t) = v_T\left(1 + e^{-gt/v_T}\right). Hvordan ser du at det er feil uten å regne på nytt?

📝Oppgave 7
(Eksamensnivå, sjanger G med numerisk hale.) En kule med masse mm skytes rett opp med startfarten v0v_0 og påvirkes av tyngden og kvadratisk luftmotstand. Regn positiv retning oppover.

a) Sett opp bevegelseslikningen for hele bevegelsen, både opp og ned, og forklar hvorfor absoluttverdien er nødvendig.

b) Skriv likningen slik den gjelder på vei opp, og slik den gjelder på vei ned.

c) Kommer kula opp med samme fart som den kom ned igjen med, eller ikke? Begrunn fysisk uten å løse likningen.

d) Kan bevegelsen beskrives med ett enkelt lukket uttrykk for v(t)v(t) hele veien? Begrunn, og si hva neste steg ville vært.

📝Oppgave 8
(Kvalitativ.) Svar i to–fire setninger per punkt.

a) To fallskjermhoppere med samme utstyr og samme kroppsstilling, men ulik masse, hopper samtidig. Hvem når bakken først, og hvorfor?

b) En medstudent skriver: «Luftmotstand gir alltid en analytisk uløselig likning, så alt må gjøres numerisk.» Rett opp utsagnet så presist du kan.

c) Hvorfor er terminalfarten en gratis poenghøster i en eksamensoppgave, og hva må du likevel huske å skrive?

Begrepsbank

Begrepsbanken er flashcard- og repetisjonsstoff — den gjentar det du nettopp har lest. Hopp trygt over ved førstegangslesing; tidsanslaget for kapitlet gjelder kjernestoffet.

Separasjon av variable
Metoden der en førsteordens likning dv/dt=f(v)h(t)dv/dt = f(v)h(t) skrives om til

dvf(v)=h(t)dt\frac{dv}{f(v)} = h(t)\,dt

og hver side integreres for seg. Metoden virker fordi dvdv og dtdt kan behandles som differensialer — det er samme regnegrep som ved substitusjon i et integral.

Alle dragmodellene i dette kapitlet er separable. Legger du integrasjonsgrensene inn med det samme (fra 00 til vv på venstre side, fra 00 til tt på høyre), slipper du både integrasjonskonstanten og et vanlig fortegnsrot.

Grensen mot fritt fall

Egenskapen at v(t)gtv(t) \to gt når t0t \to 0: i det første øyeblikket etter slipp er farten null, motstanden derfor null, og legemet faller som i vakuum.

Grensen er en obligatorisk kontroll av ethvert dragresultat. Gir uttrykket ditt noe annet enn gtgt for små tt, er det galt — og kontrollen tar ti sekunder med Taylor-utviklingen ex1xe^{-x} \approx 1 - x.

Eksponentiell tilnærming
Måten farten nærmer seg terminalfarten på: avstanden vTvv_T - v krymper med en fast faktor for hvert like langt tidsintervall.

vTv(t)=vTet/τv_T - v(t) = v_T e^{-t/\tau}

Samme matematiske form dukker opp i radioaktiv nedbrytning, i avkjøling av en kopp kaffe og i utladning av en kondensator. Det er ikke tilfeldig: alle tre er systemer der endringsraten er proporsjonal med avstanden til likevekten.

Når gjelder hvilken modell

Grovregelen er at lineær motstand gjelder for små, langsomme legemer i et seigt medium, og kvadratisk motstand for større, raskere legemer i luft eller vann.

Skillet skyldes at strømningen rundt legemet skifter karakter: ved lave farter siger mediet jevnt forbi, ved høyere farter dannes virvler bak legemet. Overgangen beskrives av et dimensjonsløst forholdstall som ikke er pensum i dette emnet.

På eksamen står det alltid i oppgaveteksten hvilken modell som gjelder. Din jobb er å bruke den riktig — ikke å avgjøre hvilken som er fysisk mest treffende.

Dragkoeffisientene kvk_v og DD
To ulike koeffisienter med ulike enheter, som aldri må forveksles:

[kv]=Nm/s=kg/s,[D]=N(m/s)2=kg/m[k_v] = \frac{\text{N}}{\text{m/s}} = \text{kg/s}, \qquad [D] = \frac{\text{N}}{(\text{m/s})^2} = \text{kg/m}

DD vokser med tverrsnittet og med luftas tetthet, og synker med en strømlinjeformet form. En fallskjerm har en DD som er over hundre ganger større enn kroppens egen.

Enhetsforskjellen er den raskeste kontrollen på at du har brukt riktig modell: kommer du ut med en terminalfart som har feil enhet, har du blandet dem.

Tanh-løsningen for kvadratisk motstand
For et legeme som slippes fra ro og faller nedover hele veien, gir den kvadratiske modellen

v(t)=vTtanh(gtvT),vT=mgDv(t) = v_T\tanh\left(\frac{gt}{v_T}\right), \qquad v_T = \sqrt{\frac{mg}{D}}

Funksjonen tanh\tanh starter med stigningstall 1 i origo og flater ut mot 1, akkurat som eksponentialløsningen — kurvene ligner, men er ikke like.

Gyldigheten er smal: løsningen krever at farten ikke skifter fortegn og at bevegelsen er endimensjonal. Bryter du en av forutsetningene, må likningen løses numerisk.

Hvorfor 2D-drag må løses numerisk
I to dimensjoner er farten v=vx2+vy2v = \sqrt{v_x^2 + v_y^2}, og motstanden fordeler seg på komponentene som

FD,x=Dvvx,FD,y=DvvyF_{D,x} = -Dv\,v_x, \qquad F_{D,y} = -Dv\,v_y

Hver komponentlikning inneholder altså begge fartskomponentene, og de kan ikke løses hver for seg. Det er denne koblingen som gjør at et skrått kast med luftmotstand ikke har noen løsning i elementære funksjoner — og som gjør Del 8 til en nødvendig del av faget, ikke et vedheng.

Luftmotstand gjør alltid negativt arbeid

Fordi kraften alltid peker mot bevegelsen, er skalarproduktet FDdr\mathbf{F}_D \cdot d\mathbf{r} negativt i hvert eneste tidssteg.

Konsekvensen er at mekanisk energi ikke er bevart når luftmotstand er med, og at en kule som kastes opp alltid kommer ned med lavere fart enn den ble kastet opp med. Se kap. 5.1 for arbeidsbegrepet.

Dette er også grunnen til at luftmotstand aldri kan «hjelpe» et legeme framover: den er en ren tapskraft i mekanisk forstand.

Strekningen med luftmotstand
Posisjonen finner du ved å integrere farten en gang til. For lineær motstand og start fra ro:

y(t)=vT[tτ(1et/τ)]y(t) = v_T\left[t - \tau\left(1 - e^{-t/\tau}\right)\right]

For store tt nærmer dette seg vT(tτ)v_T(t - \tau): en rett linje med stigningstall vTv_T, forskjøvet vTτv_T\tau nedover i forhold til den rene terminalfart-linja. Forskyvningen er «etterslepet» fra oppstartsfasen, altså den strekningen legemet tapte mens det ennå ikke hadde full fart.

Den numeriske halen

Den faste siste deloppgaven i en dragoppgave: «skriv et program, eller pseudokode, som løser likningen numerisk». Den er med i 7 av 7 sett (100 %).

Poengene gis for riktige initialbetingelser, et riktig akselerasjonsuttrykk og en riktig oppdateringsløkke — ikke for syntaks. Metoden emnet bygger på heter Euler–Cromer, og den utvikles fra grunnen i Del 8.

Det du kan gjøre allerede nå, er å skrive bevegelseslikningen på den formen koden trenger: akselerasjonen isolert på venstre side, motstandsleddet skrevet med absoluttverdi.

Halvveistiden mot terminalfart
Tiden det tar før farten har nådd halve terminalfarten, for lineær motstand og start fra ro:

t1/2=τln20,693τt_{1/2} = \tau\ln 2 \approx 0{,}693\,\tau

Den er uavhengig av vTv_T selv og bestemmes bare av tidskonstanten. Det gjør den til et praktisk mål i laboratorieforsøk: mål tiden til halv sluttfart, så har du τ\tau og dermed forholdet m/kvm/k_v uten å måtte kjenne noen av dem hver for seg.

Ensrettet bevegelse

Betingelsen som må være oppfylt for at Dvv-D\lvert v\rvert v kan forenkles til Dv2-Dv^2: farten skifter ikke fortegn i det tidsrommet du regner på.

En fallskjermhopper som synker hele veien, oppfyller den. En ball som kastes opp og faller ned igjen, gjør det ikke — der er toppunktet et skille mellom to ulike likninger.

Forenklingen er lov, men betingelsen skal skrives ut. En besvarelse som bruker Dv2-Dv^2 uten et ord om hvorfor, ser ut som felle #4 selv om regningen er riktig.

Startakselerasjonen i et fall med motstand
Ved v=0v = 0 er motstanden null uansett modell, så

a(0)=ga(0) = g

Legemet starter altså alltid som i fritt fall, og kurven v(t)v(t) har stigningstall gg i origo. Dette gjelder både for lineær og kvadratisk motstand.

Det er en av de raskeste kontrollene på et fartsuttrykk: deriver det og sett t=0t = 0. Får du noe annet enn gg, er svaret galt.

Sluttfart kontra terminalfart
Terminalfarten vTv_T er en egenskap ved systemet: den farten der kreftene balanserer. Sluttfarten er den farten legemet faktisk har når fallet er over — når det treffer bakken.

De to er bare like hvis fallhøyden er stor nok til at terminalfarten rekker å bli nådd, altså minst noen tidskonstanter. Faller et lodd i ett tiendedels sekund, er sluttfarten langt under vTv_T.

Å svare med vTv_T på et spørsmål om farten ved bakken er en klassisk feil, og den er lett å unngå: sjekk om falltiden er stor sammenliknet med τ\tau.

Skrått kast med luftmotstand

Den vanligste innpakningen av den numeriske deloppgaven: et legeme kastes skrått opp, og motstanden er kvadratisk.

Uten motstand er xx- og yy-bevegelsen uavhengige, og oppgaven løses med kastformlene fra kap. 2.2. Med motstand kobles de gjennom farten v=vx2+vy2v = \sqrt{v_x^2+v_y^2}, banen blir usymmetrisk, og rekkevidden blir kortere enn den klassiske formelen sier.

Det finnes ingen lukket løsning, og oppgaven ender derfor alltid i en numerisk deloppgave.

Parametertilpasning fra måledata

Å bestemme modellens ukjente størrelser — typisk vTv_T og τ\tau — ved å sette måledata inn i den kjente løsningsformen.

Regelen er at du trenger like mange uavhengige målinger som du har ukjente parametre. Én måling låser ikke to parametre, uansett hvor nøyaktig den er.

Et nyttig grep er å dele to måleuttrykk på hverandre: da faller den ene parameteren bort, og du står igjen med én likning i én ukjent.

Fortegnsregelen i separasjonen
Integralet der den ukjente står med negativt fortegn i nevneren:

dvvTv=lnvTv+C\int \frac{dv}{v_T - v} = -\ln\lvert v_T - v\rvert + C

Minustegnet kommer av kjerneregelen: den indre funksjonen vTvv_T - v har derivert 1-1. Dette er stedet flest mister et fortegn i hele Del 7.

Kontrollen er enkel: en løsning som vokser uten grense i stedet for å flate ut mot vTv_T, har mistet nettopp dette minustegnet.

Symbol- og formelliste
Repetisjonsoppgaver

Dette kapitlet er skrevet av Anthropics toppmodeller (Claude Opus og Claude Fable) og er foreløpig ikke manuelt gjennomgått — kvalitetskontrollen gjøres av uavhengige KI-agenter, og innmeldte feil rettes fortløpende. Funnet en feil? Meld fra, så retter vi den. Les mer om hvordan innholdet lages.

Skolesaga er en uavhengig læringsressurs og er ikke tilknyttet eller godkjent av den aktuelle utdanningsinstitusjonen. Dette er ikke offisielt studiemateriell. Les mer.