Riktningsfält och lösningskurvor 3

1:a ordningens differentialekvationer

ma5eulersmetodv1
ma5eulersmetodv1

Eulers metod demonstreras för ett antal givna 1: a ordningens differentialekvationer av formen
dy/dt=f(t,y)

Du väljer stegländen  (från ett antal givna), och sedan kan du stegvis få fram lösningskurvan samt tabell. Riktningsfältet är upprittad, och den exakta lösningen kan visas så att du kan se noggranheten på metoden.

 

Flera simuleringar från samma sida 

Dämpad svängning

2:a ordningens differentialekvationer

ma5secondorddifequv1
ma5secondorddifequv1

En kropp med massan m svänger kring sitt jämviktsläge. Den påverkas dels av en återförande kraft från fjädern:
-k·x, där k är fjäderkonstanten och x avståndet från jämviktsläget; och dels av en dämpad (friktions) kraft som är proportionell mot hastigheten: -b·x’.

Newtons andra lag (F=ma=mx’’) ger differentialekvationen:
m·x’’+b·x’+k·x=0

Simuleringen undersöker denna ekvation. Du kan välja olika värden på variablerna, och begynnelsevillkoren (x och x’’ vid t=0) väljs i fönstret till vänster. När du trycker på "»", så visas en animering av rörelsen.  

 

Flera simuleringar från samma sida

Påtvingat svängning I

2:a ordningens differentialekvationer

ma5secondorddifeqforcedv1
ma5secondorddifeqforcedv1

En kropp med massan m svänger kring sitt jämviktsläge. Den påverkas av tre krafter. En återförande kraft från fjädern: -k·x, där k är fjäderkonstanten och x avståndet från jämviktsläget; en dämpad (friktions) kraft som är proportionell mot hastigheten: -b·x’ ; och en yttre påtvingad kraft: A·cos(ωt).

Newtons andra lag (F=ma=mx’’) ger differentialekvationen:
m·x’’+b·x’+k·x=A·cos(ωt)

Simuleringen undersöker denna ekvation. Du kan välja olika värden på variablerna, och begynnelsevillkoren  (x och x’’ vid t=0) väljs i fönstret till vänster. När du trycker på "»", så visas en animering av rörelsen.

Du kan välja att visa den homogena lösningen ("Transient" =övergående; eftersom den homogena lösningen har en dämpningsfaktor e-r·t där r >0 och därmed bli liten när t är stort), den partikulära ("Steady State" eftersom det är den som blir kvar efter lång tid) eller den exakta lösningen. 

 

Flera simuleringar från samma sida

Påtvingat svängning II

2:a ordningens differentialekvationer

ma5secondorddifeqforcedv2
ma5secondorddifeqforcedv2

En kropp med massan m=1kg svänger kring sitt jämviktsläge. Den påverkas av tre krafter. En återförande kraft från fjädern: -k·x, där k är fjäderkonstanten och x avståndet från jämviktsläget; en dämpad (friktions) kraft som är proportionell mot hastigheten: -b·x’ (nedre delen av uppställningen); och en yttre påtvingad kraft genom att fjäderns hängpunkt svänger enligt h(t)=cos(ωt). Detta resulterar till en påtvingad kraft: k·h(t)= k·cos(ωt). 

Newtons andra lag (F=ma=mx’’) ger differentialekvationen:
m·x’’+b·x’+k·x=k·cos(ωt)

Simuleringen undersöker den partikulära lösningen
P(t)=c·cos(ωt)+d·sin(ωt)= A·cos(ωt+φ)
(Efter lång tid dör den homogena lösningen ut eftersom den har en dämpningsfaktor:
e-rt, och kvar blir bara den partkulära lösningen). 

I fönstren längst upp till höger visas A och φ som funktion av ω. Det intressanta är att vid ett visst värde på ωA=√(k/m-b2/2m2) så ökar amplituden A kraftigt. Vi har fått en resonans. Ju mindre dämpningsfaktorn b är desto mera påtagligt blir resonansfenomenet.

Du kan välja olika värden på variablerna b, k och ω. När du trycker på "»", så visas en animering av rörelsen.  

 

Flera simuleringar från samma sida

×
×