Tilbake
7.1

7.1 Differensiallikninger fra Newtons 2. lov

Hvordan `m d²y/dt² = ΣF` blir en differensialligning, hvilke typer som forekommer, og hvorfor «å integrere som om `y` er konstant» er meningsløst.

60 min
9 oppgaver
Differensiallikninger fra Newtons 2. lov
Din fremgang i kapitlet
0 / 9 oppgaver

Forkunnskaper — sist du var her

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

Newtons 2. lov, komponentvis, én gang per legeme:

Fx=max,Fy=may\sum F_x = ma_x, \qquad \sum F_y = ma_y

Akselerasjonen er den dobbelt deriverte av posisjonen:

v=dydt,a=dvdt=d2ydt2v = \frac{dy}{dt}, \qquad a = \frac{dv}{dt} = \frac{d^2y}{dt^2}

Det er nettopp denne sammenhengen som gjør Newtons 2. lov til en differensiallikning: venstre side inneholder krefter, høyre side inneholder en derivert av den ukjente funksjonen.

Kraftmodellene fra kap. 3.2: tyngden mgm\mathbf{g}, normalkraften N\mathbf{N}, dynamisk friksjon fd=μdN\lvert\mathbf{f}_d\rvert = \mu_d N mot bevegelsesretningen, og fjærkraften

F=k(xx0)\mathbf{F} = -k(x - x_0)

der kk er fjærkonstanten og x0x_0 er posisjonen der fjæra er avspent.

Konstant-aa-formlene v=v0+atv = v_0 + at og x=x0+v0t+12at2x = x_0 + v_0t + \tfrac12at^2 gjelder bare når akselerasjonen er konstant. Så snart en kraft avhenger av farten eller posisjonen, er de ubrukelige — og da er du i dette kapitlet.

Fra matematikken trenger du derivasjon og integrasjon: Derivasjon av vektorfunksjoner og Variabelskifte (substitusjon). R2 dekker ikke differensiallikninger som eget tema — hele typologien innføres her, fra bunnen.

Notasjonsavtale for kapitlet. Vi arbeider stort sett i én dimensjon og skriver da posisjonen yy (eller xx) og farten vv som kursive skalarer med fortegn: positiv verdi betyr «i den positive akseretningen». Krefter som opptrer som vektorer skrives fete: F\mathbf{F}, N\mathbf{N}. Tyngden skrives mgm\mathbf{g} eller Fg\mathbf{F}_galdri G\mathbf{G}, siden GG er reservert gravitasjonskonstanten senere i boka. Snordrag er TT; en svingeperiode skriver vi TpT_p for å unngå kollisjon. Vi bruker g=9,81 m/s2g = 9{,}81\ \text{m/s}^2.

Løkke 1 — når kraften ikke er konstant (~12 min)

Du slipper pedalene i en lang nedoverbakke. De første sekundene øker farten jevnt, som i fritt fall langs bakken. Men farten fortsetter ikke å øke jevnt: etter hvert flater den ut mot en verdi der du bare triller, uansett hvor lang bakken er.

Hvorfor? Fordi luftmotstanden vokser med farten. Kraften som bremser deg, avhenger av det du prøver å regne ut. Og i det øyeblikket er ikke lenger v=v0+atv = v_0 + at til noen hjelp: den formelen forutsetter at aa er konstant, og her er den ikke det.

Intuisjon: så lenge kreftene er konstante, er mekanikk en oppskrift. Så snart en kraft avhenger av fart eller posisjon, blir den en likning om en funksjon — og det er en differensiallikning.

Differensiallikning
En differensiallikning er en likning der den ukjente er en funksjon, og der likningen knytter funksjonen sammen med én eller flere av dens deriverte.

I mekanikken er den ukjente funksjonen posisjonen y(t)y(t) eller farten v(t)v(t), og likningen har formen

md2ydt2=Fm\frac{d^2y}{dt^2} = \sum F

Det avgjørende skillet: i en vanlig likning som 3x+2=113x + 2 = 11 leter du etter et tall. I en differensiallikning leter du etter en funksjon — en hel kurve y(t)y(t) som oppfyller likningen i hvert eneste tidspunkt. Å svare med et tall på et spørsmål om en differensiallikning er derfor alltid feil.

📜Bevegelseslikningen: Newtons 2. lov som differensiallikning
Sett opp Newtons 2. lov for legemet, sett inn kraftmodellene, og erstatt akselerasjonen med den andrederiverte av posisjonen:

md2ydt2=F(y, dydt, t)m\frac{d^2y}{dt^2} = \sum F\left(y,\ \frac{dy}{dt},\ t\right)

Høyresiden kan inneholde tre slags avhengighet, og det er dem du skal lete etter:

- avhengighet av posisjonen yy — fjærkraft, gravitasjon, potensialkrefter;
- avhengighet av farten dy/dtdy/dt — luftmotstand, dempere, viskøs motstand;
- eksplisitt avhengighet av tiden tt — en motor som skrus opp, en varierende ytre påkjenning.

Er høyresiden helt uavhengig av alle tre, er akselerasjonen konstant, og du kan bruke konstant-aa-formlene i stedet. Det er den ene sjekken som avgjør om du i det hele tatt har en differensiallikning å hanskes med.

Løsning av en differensiallikning

En løsning er en funksjon som gjør likningen sann for alle tidspunkt i det aktuelle intervallet. Du sjekker en påstått løsning ved å derivere den og sette den inn — akkurat som du sjekker en tallikning ved å sette inn tallet.

Innsettingsprøven er både en kontroll og en poenghøster: sensor godtar en riktig løsning som er verifisert ved innsetting selv om veien fram var kort. Og på en «vis at y(t)=y(t) = \dots løser likningen»-deloppgave er innsetting hele metoden.

✏️Eksempel 1: Fritt fall som differensiallikning

En ball slippes fra en høyde. Vi regner yy positivt oppover og ser bort fra luftmotstanden.

a) Sett opp bevegelseslikningen.

b) Vis at y(t)=y0+v0t12gt2y(t) = y_0 + v_0t - \tfrac12 gt^2 er en løsning.

c) Hva betyr y0y_0 og v0v_0 fysisk?

Metodevalg først: vi bruker Newtons 2. lov på ballen alene. Med luftmotstanden neglisjert er tyngden den eneste kraften, så oppsettet er kort — men strukturen er den samme som i alle senere oppgaver.

a) Frilegemediagrammet for ballen har én pil: tyngden mgm\mathbf{g} rett nedover, altså i negativ yy-retning. Newtons 2. lov i yy-retningen gir

md2ydt2=mgd2ydt2=gm\frac{d^2y}{dt^2} = -mg \quad \Longrightarrow \quad \frac{d^2y}{dt^2} = -g

Massen faller ut. Høyresiden er en konstant — den avhenger verken av yy, av farten eller av tiden.

b) Vi deriverer den påståtte løsningen to ganger:

dydt=v0gt,d2ydt2=g\frac{dy}{dt} = v_0 - gt, \qquad \frac{d^2y}{dt^2} = -g

Den andrederiverte er g-g, som er nøyaktig høyresiden. Likningen er altså oppfylt for alle tt, og funksjonen er en løsning.

c) Setter vi t=0t = 0, får vi y(0)=y0y(0) = y_0 og dy/dtt=0=v0dy/dt\big|_{t=0} = v_0. Altså er y0y_0 starthøyden og v0v_0 startfarten (positiv oppover). De to konstantene er ikke noe likningen bestemmer — de kommer fra hvordan ballen ble sluppet.

Benevningskontroll: venstre side i d2y/dt2=gd^2y/dt^2 = -g har enhet m/s2\text{m/s}^2, og det har gg også. Det stemmer.

Merk hva vi gjorde her: vi verifiserte en løsning uten å utlede den. Det er en helt gyldig og ofte rask metode, og på eksamen er svaret ofte oppgitt i oppgaveteksten nettopp for at du skal kunne gjøre dette.

📝Oppgave 1

Skriv opp bevegelseslikningen for hvert av de tre tilfellene, og si i hver av dem hvilken funksjon som er den ukjente. Du skal ikke løse noe.

a) En kloss glir uten friksjon nedover et fast skråplan med helning θ\theta. La ss være posisjonen langs planet, positiv nedover.

b) En kloss på et vannrett, friksjonsfritt underlag er festet til en fjær med fjærkonstant kk. La xx være utslaget fra der fjæra er avspent.

c) En kasse skyves over et vannrett gulv av en konstant kraft FF, med dynamisk friksjonstall μd\mu_d mot underlaget. La xx være posisjonen langs gulvet, positiv i skyveretningen.

Løkke 2 — klassifiseringen som gir poeng (~16 min)

Når likningen står der, er neste steg alltid det samme: hva slags likning er dette? Sensor gir egne poeng for klassifiseringen, og grunnen er praktisk — det er klassifiseringen som bestemmer hvilken metode du kan bruke. Hopper du over den, gjetter du på metode.

Fire spørsmål, i denne rekkefølgen. De tar tjue sekunder til sammen når de sitter.

Ordenen til en differensiallikning
Ordenen er nummeret på den høyeste deriverte som opptrer i likningen.

mdvdt=mgkvver av første ordenm\frac{dv}{dt} = mg - k_v v \quad \text{er av \textbf{første} orden}

md2ydt2=kyer av andre ordenm\frac{d^2y}{dt^2} = -ky \quad \text{er av \textbf{andre} orden}

Newtons 2. lov skrevet for posisjonen gir alltid andre orden, fordi akselerasjonen er den andrederiverte. Skriver du den i stedet for farten, og posisjonen ikke inngår, faller ordenen til én. Det trikset er halve verktøykassa i luftmotstandsoppgaver.

Ordenen bestemmer hvor mange ubestemte konstanter løsningen får — og dermed hvor mange initialbetingelser du trenger.

Lineær differensiallikning

En likning er lineær dersom den ukjente funksjonen og alle dens deriverte opptrer i første potens, hver for seg, uten å multipliseres med hverandre og uten å stå inni en funksjon som sin\sin,  \sqrt{\ } eller en potens.

Lineært:  d2ydt2+3dydt+5y=2t\ \dfrac{d^2y}{dt^2} + 3\dfrac{dy}{dt} + 5y = 2t

Ikke lineært:  d2θdt2=gsinθ\ \dfrac{d^2\theta}{dt^2} = -\dfrac{g}{\ell}\sin\theta \quad (den ukjente står inni en sinus)

Ikke lineært:  mdvdt=mgDvv\ m\dfrac{dv}{dt} = mg - D\lvert v\rvert v \quad (den ukjente er kvadrert)

Hvorfor det betyr noe: lineære likninger med konstante koeffisienter har en fullstendig, mekanisk løsningsoppskrift, som du finner i Rottmann. Ikke-lineære har det som hovedregel ikke — de må enten lineariseres (kap. 7.2) eller løses numerisk.

Homogen og inhomogen likning
Samle alle ledd som inneholder den ukjente funksjonen på venstre side. Står det null igjen på høyre side, er likningen homogen. Står det et ledd som ikke inneholder den ukjente — en konstant kraft, en ytre påkjenning — er den inhomogen.

d2udt2+ω02u=0homogen\frac{d^2u}{dt^2} + \omega_0^2 u = 0 \qquad \text{homogen}

dvdt+kvmv=ginhomogen (tyngden er drivleddet)\frac{dv}{dt} + \frac{k_v}{m}v = g \qquad \text{inhomogen (tyngden er drivleddet)}

Fysisk lesning: det inhomogene leddet er den kraften som holdes utenfra og som ikke forsvinner selv om systemet står stille. Ofte kan du gjøre en inhomogen likning homogen ved å måle utslaget fra likevekten i stedet for fra origo — det trikset er kjernen i kap. 7.2.

Separabel differensiallikning
En førsteordens likning er separabel dersom den kan skrives slik at alt som inneholder den ukjente står på den ene siden, og alt som inneholder den frie variabelen på den andre:

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

Da kan du integrere hver side for seg. Alle de førsteordens luftmotstandslikningene i kap. 7.3 er separable, også de ikke-lineære — separabel og lineær er to helt ulike egenskaper, og en likning kan være det ene uten å være det andre.

Andreordens likninger er ikke separable på denne måten; der er karakteristisk likning verktøyet.

Klassifiseringen peker rett på metoden. Denne tabellen er hele beslutningstreet for Del 7, og den er verdt å kunne uten å slå opp:

Likningen er …Typisk fysikkMetodeHvor i boka
1. orden, separabelfall med lineær eller kvadratisk luftmotstandseparér variablene og integrerkap. 7.3
2. orden, lineær, homogen, konstante koeffisienterfjær, pendel for små utslag, dempet svingningkarakteristisk likningkap. 7.2
2. orden, lineær, inhomogenhengende fjær med tyngdeflytt origo til likevekten, så som overkap. 7.2
ikke-lineær, men nær en likevektpendel, ikke-lineær fjærTaylor-utvikle og linearisérkap. 7.2
ikke-lineær, langt fra likevekt, eller 2Dkast med luftmotstand, planetbanenumerisk løsningDel 8

Den nederste raden er ikke en nødløsning. Numerisk løsning er en egen, fullverdig metode i dette emnet, og den er testet i alle sju settene. Del 8 bygger den fra grunnen. Frem til da skal du kunne si i klartekst at likningen må løses numerisk, og hvorfor — det er en vanlig deloppgave i seg selv.
✏️Eksempel 2: Klassifiser fire bevegelseslikninger

Klassifiser hver av likningene etter orden, lineær/ikke-lineær og homogen/inhomogen, og si hvilken metode klassifiseringen peker på. Alle konstantene er positive.

a)  mdvdt=mgkvv\ m\dfrac{dv}{dt} = mg - k_v v — lodd som faller med lineær luftmotstand

b)  md2ydt2=ky\ m\dfrac{d^2y}{dt^2} = -ky — kloss på vannrett fjær

c)  md2ydt2=kybdydt\ m\dfrac{d^2y}{dt^2} = -ky - b\dfrac{dy}{dt} — samme fjær, med en demper

d)  md2θdt2=mgsinθ\ m\ell\dfrac{d^2\theta}{dt^2} = -mg\sin\theta — pendel med lengde \ell

Metodevalg: vi kjører de fire spørsmålene i rekkefølge på hver likning. Klassifiseringen skrives alltid ut i klartekst — det er den som gir poeng, ikke bare svaret etterpå.

a) Høyeste deriverte er dv/dtdv/dt, altså første orden. Den ukjente vv opptrer i første potens alene, så likningen er lineær. Flytter vi kvvk_v v over, står det mdv/dt+kvv=mgm\,dv/dt + k_v v = mg, og høyresiden mgmg inneholder ikke vv — den er inhomogen. Den er også separabel: alt med vv kan samles på den ene siden.
Metode: separér og integrer.

b) Høyeste deriverte er d2y/dt2d^2y/dt^2, altså andre orden. Den ukjente står i første potens: lineær. Alt inneholder yy, høyresiden blir null: homogen, og med konstante koeffisienter.
Metode: karakteristisk likning.

c) Fortsatt andre orden, fortsatt lineær (både yy og dy/dtdy/dt er i første potens), fortsatt homogen. Dempeleddet endrer ikke klassifiseringen — det endrer bare hvilke røtter den karakteristiske likningen får.
Metode: karakteristisk likning, samme oppskrift som b).

d) Andre orden. Men θ\theta står inni en sinus, så likningen er ikke lineær. Den er homogen i den forstand at høyresiden er null når θ=0\theta = 0, men den vanlige oppskriften gjelder ikke.
Metode: Taylor-utvikle sinθθ\sin\theta \approx \theta for små utslag, og løs den lineariserte likningen; for store utslag må den løses numerisk.

Sensorblikk: legg merke til at c) og d) ser like «kompliserte» ut, men bare den ene har en lukket løsningsformel. Det er nettopp det klassifiseringen avslører — og derfor den skal skrives før du velger metode.

📝Oppgave 2

Klassifiser hver av likningene etter (i) orden, (ii) lineær eller ikke-lineær, (iii) homogen eller inhomogen. Alle konstanter er positive.

a)  dvdt=Dmvv\ \dfrac{dv}{dt} = -\dfrac{D}{m}\lvert v\rvert v

b)  d2xdt2+4dxdt+9x=0\ \dfrac{d^2x}{dt^2} + 4\dfrac{dx}{dt} + 9x = 0

c)  d2ydt2=gkvmdydt\ \dfrac{d^2y}{dt^2} = -g - \dfrac{k_v}{m}\dfrac{dy}{dt}

d)  d2xdt2=kmxcmx3\ \dfrac{d^2x}{dt^2} = -\dfrac{k}{m}x - \dfrac{c}{m}x^3

Løkke 3 — initialbetingelsene (~12 min)

En differensiallikning beskriver loven. Den sier hva som styrer bevegelsen, men ikke hvor bevegelsen startet. Derfor har en differensiallikning uendelig mange løsninger — én for hver måte du kunne ha satt systemet i gang på.

Intuisjon: likningen d2y/dt2=gd^2y/dt^2 = -g gjelder like godt for en ball du slipper som for en ball du kaster oppover. Fysikken er den samme; startbetingelsene er ulike.

Initialbetingelse
En initialbetingelse er en opplysning om systemets tilstand ved et bestemt tidspunkt, typisk t=0t = 0: hvor legemet var, og hvor fort det gikk.

y(0)=y0,dydtt=0=v(0)=v0y(0) = y_0, \qquad \frac{dy}{dt}\bigg|_{t=0} = v(0) = v_0

Initialbetingelsene kommer aldri fra likningen — de kommer fra oppgaveteksten, og de skal begrunnes fysisk: «kula slippes fra ro, altså v(0)=0v(0) = 0», «loddet henger i ro i likevekt når snora kuttes, altså u(0)=0u(0) = 0 og v(0)=0v(0) = 0 målt fra den nye likevekten».

En setning som begrunner initialbetingelsen fysisk er selvstendig poenggivende. Å bare skrive «v(0)=0v(0) = 0» uten å si hvorfor, er å la poeng ligge igjen.

📜Like mange initialbetingelser som ordenen

Løsningen av en differensiallikning av orden nn inneholder nn ubestemte konstanter, og du trenger nøyaktig nn initialbetingelser for å bestemme dem.

- Første orden (fart med luftmotstand): én konstant, én initialbetingelse — typisk v(0)v(0).
- Andre orden (fjær, pendel): to konstanter, to initialbetingelser — typisk y(0)y(0) og v(0)v(0).

Intuisjon: ordenen teller hvor mange ganger du må «integrere opp», og hver integrasjon føder én ukjent konstant.

Praktisk konsekvens: ser du en andreordens likning og bare én oppgitt startopplysning, mangler du noe — les oppgaveteksten en gang til. Ofte ligger den andre betingelsen skjult i en formulering som «slippes fra ro» eller «passerer likevektspunktet».

Generell og partikulær løsning
Den generelle løsningen er hele familien av funksjoner som oppfyller likningen, med de ubestemte konstantene stående. Den partikulære løsningen er det ene medlemmet av familien som i tillegg oppfyller initialbetingelsene.

y(t)=y0+v0t12gt2er generell sa˚ lenge y0,v0 er friey(t) = y_0 + v_0t - \tfrac12gt^2 \quad \text{er generell så lenge } y_0, v_0 \text{ er frie}

På eksamen forventes begge deler: først den generelle løsningen (den viser at du kan metoden), deretter innsettingen av initialbetingelsene (den viser at du kan svare på denne oppgaven). Poengene er delt mellom de to stegene.

✏️Eksempel 3: Fra generell til partikulær løsning

En stein kastes rett oppover fra y=0y = 0 med farten 6,0 m/s6{,}0\ \text{m/s}. Luftmotstanden neglisjeres, og yy regnes positivt oppover.

a) Skriv den generelle løsningen av bevegelseslikningen.

b) Bestem konstantene fra initialbetingelsene, og begrunn dem fysisk.

c) Når er steinen på det høyeste, og hvor høyt kommer den?

Metodevalg: likningen er d2y/dt2=gd^2y/dt^2 = -g med konstant høyreside, altså den enkleste typen. Vi integrerer to ganger og lar konstantene stå til initialbetingelsene settes inn.

a) Første integrasjon gir farten, andre gir posisjonen:

dydt=gt+C1,y(t)=12gt2+C1t+C2\frac{dy}{dt} = -gt + C_1, \qquad y(t) = -\tfrac12gt^2 + C_1t + C_2

To ubestemte konstanter — som seg hør og bør for en andreordens likning.

b) Steinen forlater hånden i origo, så y(0)=0y(0) = 0, og den kastes oppover med 6,0 m/s6{,}0\ \text{m/s}, altså i positiv retning: v(0)=+6,0 m/sv(0) = +6{,}0\ \text{m/s}. Innsatt i uttrykkene over:

y(0)=C2=0,dydtt=0=C1=6,0 m/sy(0) = C_2 = 0, \qquad \frac{dy}{dt}\bigg|_{t=0} = C_1 = 6{,}0\ \text{m/s}

y(t)=6,0t4,905t2(meter, med t i sekunder)y(t) = 6{,}0\,t - 4{,}905\,t^2 \quad (\text{meter, med } t \text{ i sekunder})

c) På toppen er farten null:

dydt=6,09,81t=0t=0,6116 s\frac{dy}{dt} = 6{,}0 - 9{,}81t = 0 \quad \Longrightarrow \quad t = 0{,}6116\ \text{s}

y(0,6116)=6,00,61164,9050,61162=3,6701,835=1,83 my(0{,}6116) = 6{,}0 \cdot 0{,}6116 - 4{,}905 \cdot 0{,}6116^2 = 3{,}670 - 1{,}835 = 1{,}83\ \text{m}

Steinen er på det høyeste etter 0,61 s0{,}61\ \text{s} og kommer 1,83 m1{,}83\ \text{m} over utgangspunktet.

Kontroll med en annen metode: energibetraktningen gir 12mv02=mgh\tfrac12mv_0^2 = mgh, altså h=v02/(2g)=36/19,62=1,83 mh = v_0^2/(2g) = 36/19{,}62 = 1{,}83\ \text{m} ✓. To uavhengige veier til samme tall er et godt tegn, og sensor godtar begge.

Merk: her kunne vi integrere rett fram, fordi høyresiden var en ren konstant. Det er nettopp den situasjonen som ikke gjelder i resten av kapitlet.

📝Oppgave 3
En kloss på en fjær beveger seg etter likningen

d2xdt2=25x\frac{d^2x}{dt^2} = -25\,x

med xx i meter og tt i sekunder.

a) Vis ved innsetting at x(t)=Acos(5t)+Bsin(5t)x(t) = A\cos(5t) + B\sin(5t) er en løsning for alle verdier av AA og BB.

b) Klossen slippes fra ro i x=0,12 mx = 0{,}12\ \text{m}. Bestem AA og BB, og begrunn begge initialbetingelsene fysisk.

c) Hvor mange initialbetingelser trengs, og hvorfor akkurat så mange?

Løkke 4 — metodefellen som gir null poeng (~12 min)

Nå kommer den viktigste enkeltbiten i kapitlet. Se på likningen for et lodd som faller med lineær luftmotstand:

mdvdt=mgkvvm\frac{dv}{dt} = mg - k_v v

En fristende, men feil, tanke er: «jeg integrerer begge sider med hensyn på tiden». Gjør du det og behandler vv på høyresiden som om den var en konstant, får du

mv=(mgkvv)tv=gt1+kvt/mm\,v = (mg - k_v v)\,t \quad \Longrightarrow \quad v = \frac{gt}{1 + k_v t/m}

Det ser ut som et svar. Det har riktig enhet. Det starter til og med i null og flater ut mot riktig sluttverdi. Og det gir null poeng, fordi funksjonen ikke løser likningen.

Hvorfor: når du integrerer kvvdtk_v v\,dt, kan du ikke sette vdt=vt\int v\,dt = v\,t. Det ville krevd at vv var konstant — og hele poenget med likningen er at vv endrer seg. Du bruker altså svaret ditt som forutsetning for å finne svaret.

Hvor galt blir det? Med et lodd der terminalfarten er 2,62 m/s2{,}62\ \text{m/s}, gir den gale metoden 1,38 m/s1{,}38\ \text{m/s} etter 0,300{,}30 sekunder, mens den riktige løsningen gir 1,77 m/s1{,}77\ \text{m/s}. Avviket er 22 %22\ \% — midt i det tidsrommet oppgaven som regel spør om. At de to kurvene ender på samme verdi når tt \to \infty, er en tilfeldighet ved akkurat denne modellen, og redder ingenting: sensor retter metoden, ikke sluttverdien.

📝Oppgave 4
En student skal finne farten til et lodd som faller i en væske. Bevegelseslikningen er

mdvdt=mgkvvm\frac{dv}{dt} = mg - k_v v

Studenten skriver: «Jeg integrerer begge sider fra 00 til tt og får mv=mgtkvvtmv = mgt - k_v vt, altså v=gt1+kvt/mv = \dfrac{gt}{1 + k_vt/m}

a) Forklar presist hvilket steg som er ugyldig, og hvorfor.

b) Test svaret ved å sette det inn i differensiallikningen for et lite tidspunkt. Hva ser du?

c) Studenten innvender: «Men uttrykket mitt gir jo riktig terminalfart når tt blir stor!» Er det et gyldig forsvar?

Likevektsløsning
En likevektsløsning er en konstant funksjon som oppfyller differensiallikningen. Du finner den ved å sette alle deriverte lik null og løse den vanlige likningen som blir igjen.

mdvdt=mgkvv dv/dt=0 v=mgkvm\frac{dv}{dt} = mg - k_v v \quad \xrightarrow{\ dv/dt = 0\ } \quad v = \frac{mg}{k_v}

Fysisk er dette den tilstanden der kreftene balanserer og ingenting endrer seg mer. For et fallende legeme kalles den terminalfarten, og den behandles grundig i kap. 7.3.

Dette er en gratis poenghøster: likevektsløsningen kan du finne uten å løse likningen, og den er ofte et eget spørsmål i deloppgave a) eller b). Den fungerer også som kontroll: den fulle løsningen skal nærme seg likevektsløsningen når tiden går, hvis likevekten er stabil.

📝Oppgave 5
Et lodd faller i olje. Bevegelseslikningen er

mdvdt=mgkvvm\frac{dv}{dt} = mg - k_v v

med m=0,0200 kgm = 0{,}0200\ \text{kg} og kv=0,0800 kg/sk_v = 0{,}0800\ \text{kg/s}. Positiv retning er nedover.

a) Finn likevektsløsningen, og forklar fysisk hva den betyr.

b) Hva er akselerasjonen i det øyeblikket loddet slippes fra ro?

c) Hva er akselerasjonen når farten er halvparten av likevektsverdien?

d) Uten å løse likningen: skisser med ord hvordan v(t)v(t) ser ut, og begrunn formen ut fra a)–c).

Løkke 5 — når analytisk løsning ryker, og hva som da skjer (~14 min)

Det finnes ingen generell metode for å løse differensiallikninger. Det som finnes, er noen få familier vi kan løse eksakt — og de dekker heldigvis mesteparten av det som testes på eksamen:

- førsteordens separable likninger, som luftmotstandslikningene i kap. 7.3;
- andreordens lineære med konstante koeffisienter, som fjæra og den lineariserte pendelen i kap. 7.2.

Alt annet — kvadratisk luftmotstand i to dimensjoner, pendelen med store utslag, en planet i bane rundt sola — må løses numerisk.

Det numeriske sporet er en egen del av denne boka (Del 8), og det er ikke en nødløsning: numerisk Python er testet i 7 av 7 sett (100 %), og det er nesten alltid siste deloppgave i nettopp ODE-oppgaven. Metoden heter Euler–Cromer, og den bygger på at du deler tiden i små steg og oppdaterer fart og posisjon steg for steg. Del 8 bygger den fra grunnen av. Det du skal ta med deg herfra, er å skrive likningen på en form som er klar for numerisk behandling: akselerasjonen isolert på venstre side, alt annet uttrykt ved posisjon, fart og tid.

Tilstandsform: én andreordens likning blir til to førsteordens
Enhver andreordens bevegelseslikning kan skrives om til to koblede førsteordens likninger ved å innføre farten som egen ukjent:

d2ydt2=f(y,dydt){dydt=vdvdt=f(y,v)\frac{d^2y}{dt^2} = f\left(y, \frac{dy}{dt}\right) \quad \Longleftrightarrow \quad \begin{cases} \dfrac{dy}{dt} = v \\[6pt] \dfrac{dv}{dt} = f(y, v)\end{cases}

Paret (y,v)(y, v) kalles systemets tilstand: kjenner du begge to i ett øyeblikk, er hele framtiden bestemt av likningen.

Denne omskrivingen er ikke bare kosmetisk. Den er selve formen enhver numerisk metode arbeider i — også Euler–Cromer i Del 8, som oppdaterer først vv, så yy. Kan du gjøre omskrivingen raskt, har du allerede halve den numeriske deloppgaven i boks.

✏️Eksempel 4 (eksamensnivå): fjær med luftmotstand

En kule med masse mm henger i en fjær med fjærkonstant kk og svinger loddrett i lufta. Fjæra er avspent ved y=0y = 0, og yy regnes positivt oppover. I tillegg til tyngden og fjærkraften virker en luftmotstand med tallverdi DvvD\lvert v\rvert v, rettet mot bevegelsen.

a) Sett opp bevegelseslikningen.

b) Klassifiser den.

c) Finn likevektsposisjonen.

d) Kan likningen løses analytisk? Begrunn, og si hva neste steg ville vært.

e) Skriv likningen på tilstandsform.

Metodevalg: ett frilegemediagram for kula, Newtons 2. lov loddrett, deretter klassifisering før noe annet. Kraften fra luftmotstanden skrives med absoluttverdi nettopp fordi kula skifter retning under svingningen.

a) Tre krefter virker på kula: fjærkraften ky-ky (peker mot avspent lengde), tyngden mg-mg (nedover) og luftmotstanden Dvv-D\lvert v\rvert v (alltid mot bevegelsen). Newtons 2. lov loddrett:

md2ydt2=kymgDdydtdydtm\frac{d^2y}{dt^2} = -ky - mg - D\left\lvert\frac{dy}{dt}\right\rvert\frac{dy}{dt}

Hvorfor absoluttverdi? Skriver du Dv2-Dv^2, blir motstandsleddet negativt uansett retning, og kula ville bli bremset på vei ned og skjøvet nedover på vei opp. Med Dvv-D\lvert v\rvert v snur leddet fortegn samtidig med farten, og kraften peker alltid mot bevegelsen. Dette er felle #4 i feilregisteret, og den koster poeng både i utledningen og senere i koden.

b) Andre orden (høyeste deriverte er d2y/dt2d^2y/dt^2). Ikke lineær, på grunn av vv\lvert v\rvert v. Inhomogen, fordi tyngdeleddet mg-mg ikke inneholder den ukjente.

c) I likevekt står kula i ro: både farten og akselerasjonen er null, så luftmotstanden faller bort:

0=kyLmgyL=mgk0 = -ky_L - mg \quad \Longrightarrow \quad y_L = -\frac{mg}{k}

Altså henger kula en strekk mg/kmg/k under avspent lengde — negativ yy, som seg hør og bør når yy er positiv oppover.

d) Nei. Likningen er andreordens og ikke-lineær, og den faller utenfor begge de familiene vi kan løse eksakt: den er ikke separabel (andre orden), og karakteristisk likning gjelder bare for lineære likninger. Neste steg er numerisk løsning: del tiden i små steg og oppdater fart og posisjon steg for steg med Euler–Cromer-metoden, som bygges i Del 8.

Delvis analytisk innsikt får vi likevel: for svært små svingninger er farten liten, og vv\lvert v\rvert v er da mye mindre enn vv selv. Da nærmer likningen seg den udempede fjærlikningen fra kap. 7.2, som vi kan løse. Det er alltid verdt å si i en slik deloppgave — det viser at du kjenner grensetilfellet.

e) Med v=dy/dtv = dy/dt:

{dydt=vdvdt=kmygDmvv\begin{cases} \dfrac{dy}{dt} = v \\[6pt] \dfrac{dv}{dt} = -\dfrac{k}{m}y - g - \dfrac{D}{m}\lvert v\rvert v \end{cases}

Sensorblikk: legg merke til at DD deles på mm i akselerasjonslinja. Det er en av de vanligste feilene når et nytt kraftledd legges til — kraften må gjennom Newtons 2. lov før den blir en akselerasjon. To initialbetingelser trengs, y(0)y(0) og v(0)v(0), og de er samtidig startverdiene den numeriske løkka trenger.

📝Oppgave 6

En fallskjermhopper i fritt fall før skjermen åpnes påvirkes av tyngden og en kvadratisk luftmotstand med tallverdi Dv2Dv^2. Vi regner positiv retning nedover, og hopperen faller hele tiden nedover.

a) Sett opp bevegelseslikningen for farten v(t)v(t), og forklar hvorfor du kan skrive motstandsleddet som Dv2-Dv^2 i akkurat denne situasjonen.

b) Klassifiser likningen.

c) Finn terminalfarten uttrykt ved mm, gg og DD, og regn den ut for m=80,0 kgm = 80{,}0\ \text{kg} og D=0,250 kg/mD = 0{,}250\ \text{kg/m}.

d) Hvorfor må motstandsleddet skrives Dvv-D\lvert v\rvert v i det generelle tilfellet?

📝Oppgave 7

En kloss med masse mm ligger på et vannrett, friksjonsfritt bord og er festet til en fjær med fjærkonstant kk. Fjæra er avspent ved x=0x = 0. I tillegg virker en dempende kraft bdx/dt-b\,dx/dt fra en oljedemper, og en ytre motor drar i klossen med kraften F0cos(Ωt)F_0\cos(\Omega t).

a) Sett opp bevegelseslikningen.

b) Klassifiser den fullstendig.

c) Hvor mange initialbetingelser trengs, og hva ville de være fysisk?

d) Hva ville likningen bli hvis motoren skrus av, og hva ville den bli hvis både motoren og demperen fjernes?

e) Skriv den fulle likningen på tilstandsform.

📝Oppgave 8
(Eksamensnivå, sjanger G — differensiallikningsoppgaven.) En dykkerklokke med masse mm senkes i vann. Kreftene på den er tyngden mgmg nedover, en konstant oppdrift BB oppover, og en lineær vannmotstand med tallverdi kvvk_v\lvert v\rvert rettet mot bevegelsen. Regn positiv retning nedover, og anta mg>Bmg > B.

a) Sett opp bevegelseslikningen for farten v(t)v(t), gyldig så lenge klokka synker.

b) Klassifiser likningen.

c) Finn synkefarten i likevekt, og sett opp betingelsen for at klokka i det hele tatt synker.

d) Klokka slippes fra ro. Vil den nærme seg likevektsfarten ovenfra eller nedenfra? Begrunn uten å løse likningen.

e) Oppdriften økes langsomt til B>mgB > mg. Hva skjer med likevektsløsningen, og hva betyr fortegnet fysisk?

📝Oppgave 9
(Kvalitativ.) To studenter diskuterer utsagnet: «En differensiallikning har alltid uendelig mange løsninger, så den kan ikke forutsi noe.»

a) Forklar hva som er riktig og hva som er galt i utsagnet.

b) En kule svinger på en fjær. Hvorfor er det ikke nok å kjenne posisjonen ved t=0t = 0 for å vite hva som skjer videre?

c) For en førsteordens likning som mdv/dt=mgkvvm\,dv/dt = mg - k_v v holder det derimot med én opplysning. Forklar hvorfor forskjellen ligger i ordenen og ikke i fysikken.

Svar i to–fire setninger per punkt.

Begrepsbank

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

Startverdiproblem

En differensiallikning sammen med akkurat nok initialbetingelser til å bestemme alle konstantene. Det er startverdiproblemet, ikke likningen alene, som har ett entydig svar.

Nesten alle oppgaver i FYS1100 er startverdiproblemer: likningen kommer fra Newtons 2. lov, betingelsene fra formuleringer som «slippes fra ro» eller «kastes med farten v0v_0».

Randbetingelse (til forskjell fra initialbetingelse)

En initialbetingelse gir tilstanden i ett tidspunkt, typisk t=0t = 0. En randbetingelse gir opplysninger i to ulike punkter — for eksempel «kula er i x=0x = 0 ved t=0t = 0 og ved t=2,0 st = 2{,}0\ \text{s}».

I mekanikkens tidsutviklingsproblemer er det initialbetingelser som gjelder, og det er dem eksamen spør etter. Randbetingelser hører til i andre deler av fysikken (bølger på en streng, statiske bjelker) og er kjenne-stoff her.

Autonom differensiallikning
En likning der tiden ikke opptrer eksplisitt på høyresiden — bare gjennom den ukjente funksjonen og dens deriverte.

mdvdt=mgkvver autonomm\frac{dv}{dt} = mg - k_v v \quad \text{er autonom}

md2xdt2=kx+F0cos(Ωt)er det ikkem\frac{d^2x}{dt^2} = -kx + F_0\cos(\Omega t) \quad \text{er det ikke}

Autonome likninger har den praktiske egenskapen at det er likegyldig når du starter klokka: samme starttilstand gir samme forløp. Det er også de som har likevektsløsninger i vanlig forstand.

Konstante koeffisienter
En lineær likning har konstante koeffisienter når tallene foran den ukjente og dens deriverte ikke avhenger av tiden:

ad2ydt2+bdydt+cy=0med a,b,c konstantea\frac{d^2y}{dt^2} + b\frac{dy}{dt} + cy = 0 \quad \text{med } a, b, c \text{ konstante}

Dette er den ene klassen der en fullstendig og mekanisk løsningsoppskrift finnes — den karakteristiske likningen. Fysisk svarer konstante koeffisienter til at massen, fjærkonstanten og dempingen ikke endrer seg underveis, noe som gjelder i praktisk talt alle oppgavene i dette emnet.

Analytisk løsning

En løsning skrevet som en formel i kjente funksjoner: potenser, eksponentialfunksjoner, sinus og cosinus. Den kan deriveres, grensebetraktes og settes inn i andre uttrykk.

Analytiske løsninger finnes bare for spesielle likninger. Når de finnes, er de å foretrekke: de viser hvordan svaret avhenger av parametrene, ikke bare hva tallet blir for én bestemt kalibrering.

Numerisk løsning

En tabell av tallverdier for den ukjente funksjonen i et endelig antall tidspunkt, beregnet steg for steg fra initialbetingelsene. Du får ingen formel, men du får kurven — og du får den for likninger som ikke kan løses analytisk.

Metoden emnet bygger på heter Euler–Cromer og utvikles i Del 8. Der lærer du også hvorfor akkurat den metoden foretrekkes framfor den enda enklere Forward Euler-metoden.

Diskretisering
Å erstatte den kontinuerlige tiden med et rutenett av tidspunkt t0,t1,t2,t_0, t_1, t_2, \dots med fast avstand Δt\Delta t, og de deriverte med differanser mellom nabopunkter:

dvdtvi+1viΔt\frac{dv}{dt} \approx \frac{v_{i+1} - v_i}{\Delta t}

Dette er det første steget i enhver numerisk metode, og forklarer hvorfor svaret blir en tallrekke og ikke en formel. Detaljene, og hva som skjer når Δt\Delta t velges for stort, hører til Del 8.

Superposisjonsprinsippet for lineære homogene likninger

Er y1(t)y_1(t) og y2(t)y_2(t) begge løsninger av samme lineære homogene likning, er også enhver sum Ay1(t)+By2(t)Ay_1(t) + By_2(t) en løsning.

Det er dette prinsippet som gjør at den generelle løsningen av en andreordens lineær homogen likning kan skrives som en sum av to grunnløsninger med hver sin konstant — for eksempel Acosω0t+Bsinω0tA\cos\omega_0 t + B\sin\omega_0 t. Prinsippet gjelder ikke for ikke-lineære likninger, og det er en av de viktigste praktiske forskjellene mellom de to klassene.

Kraftlov: hva høyresiden kan avhenge av

Høyresiden i en bevegelseslikning kan avhenge av posisjonen (fjærkraft, gravitasjon), av farten (luftmotstand, demping) og av tiden (en ytre påkjenning som varierer) — eller av flere av delene samtidig.

Å identifisere hvilken av de tre du har, er det som avgjør både klassifiseringen og metoden. Avhenger høyresiden av ingen av dem, er akselerasjonen konstant, og du er tilbake i kap. 2.1.

Verifisering ved innsetting

Å kontrollere en påstått løsning ved å derivere den så mange ganger som ordenen krever og sette resultatet inn i likningen. Går likningen opp for alle tt, er funksjonen en løsning.

Dette er både en kontroll og en fullverdig løsningsmetode når svaret er oppgitt: en «vis at y(t)=y(t) = \dots oppfyller likningen»-deloppgave besvares nøyaktig slik, og det er full uttelling for det.

Dimensjonskontroll av en differensiallikning

Hvert ledd i en differensiallikning må ha samme enhet. I md2y/dt2=kymgm\,d^2y/dt^2 = -ky - mg har alle tre ledd enheten newton.

Kontrollen avslører raskt en glemt masse eller en kraft som ikke er delt på mm i akselerasjonslinja. Den koster fem sekunder og er en av de billigste kvalitetssikringene som finnes — samme teknikk som benevningskontrollen i kap. 1.1.

Entydighet: én tilstand, én framtid

For de bevegelseslikningene som opptrer i dette emnet gjelder at én gitt starttilstand gir nøyaktig én løsning. To ulike bevegelser kan altså aldri ha samme posisjon og samme fart i samme øyeblikk.

Praktisk konsekvens: får du to ulike svar på samme startverdiproblem, er minst ett av dem galt. Dette er kjenne-stoff — de matematiske betingelsene bak er ikke pensum, men konklusjonen brukes hele tiden.

Eksplisitt form: akselerasjonen isolert
En bevegelseslikning står på eksplisitt form når den høyeste deriverte er alene på venstre side:

d2ydt2=1mF(y,dydt,t)\frac{d^2y}{dt^2} = \frac{1}{m}\sum F\left(y, \frac{dy}{dt}, t\right)

Dette er formen både klassifiseringen og den numeriske behandlingen tar utgangspunkt i. Første grep etter at Newtons 2. lov er satt opp: del på massen og få akselerasjonen alene.

Parameter kontra ukjent

I md2y/dt2=kymgm\,d^2y/dt^2 = -ky - mg er mm, kk og gg parametre — kjente tall som beskriver systemet. y(t)y(t) er den ukjente funksjonen.

Blandingen av de to er en klassisk kilde til rot: å «løse for kk» i en bevegelseslikning gir ingen mening med mindre oppgaven eksplisitt ber om å bestemme fjærkonstanten fra en observert bevegelse. Spør alltid: hva er det jeg leter etter — et tall eller en funksjon?

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.