Numerische Lösungen von Differentialgleichungen in drei Dimensionen
Ein Objekt, das sich in drei Dimensionen bewegt, wird durch sechs Variablen beschrieben: $x$, $y$, $z$, $v_x$, $v_y$ und $v_z$. Wenn die Kraft auf das Objekt bekannt ist, dann kann die Bewegung durch sechs Differentialgleichungen erster Ordnung beschrieben werden,
$\large \frac{dx}{dt}=v_x$ $\large \frac{dv_x}{dt}=F_x(x,y,z,v_x,v_y,v_z,t)/m$
$\large \frac{dy}{dt}=v_y$ $\large \frac{dv_y}{dt}=F_y(x,y,z,v_x,v_y,v_z,t)/m$
$\large \frac{dz}{dt}=v_z$ $\large \frac{dv_z}{dt}=F_z(x,y,z,v_x,v_y,v_z,t)/m$
Diese Gleichungen können Schritt für Schritt numerisch gelöst werden. Bei bekannten Anfangsbedingungen für $\vec{x}(t_0)$ und $\vec{v}(t_0)$ ergibt sich als Näherung, für $\vec{x}$ und $\vec{v}$ kurze Zeit $\Delta t$ später,
$\large \vec{x}(t_0+\Delta t) \approx \frac{d\vec{v}}{dt}|_{t_0}\Delta t$ and $\large \vec{v}(t_0+\Delta t) \approx \vec{F}(x(t_0),v_x(t_0),t_0)\Delta t/m.$
Sobald eine Näherung für die Position und die Geschwindigkeit des Objekts zum Zeitpunkt $t_0 + \Delta t$ berechnet ist, kann die Näherung der Position und der Geschwindigkeit zum Zeitpunkt $t_0 + 2\Delta t$ berechnet werden. Durch wiederholen dieses Schrittes kann eine Tabelle mit Zeiten $(t)$, Positionen $(x,y,z)$, der Geschwindigkeit $(v_x,v_y,v_z)$ und der Kräfte $(F_x,F_y,F_z)$ erstellt werden. Andere Größen wie die kinetische Energie, die geleistete Arbeit und die Momentanleistung können berechnet werden. Die kinetische Energie beträgt $E_{\text{kin}}=\frac{v_x^2+v_y^2+v_z^2}{2m}$ [J]. Wenn eine Kraft ein Teilchen um eine Strecke $d\vec{r}$ bewegt, führt sie ein kleinen Teil $dW=\int\vec{F}\cdot d\vec{r}$ [J] Arbeit aus. Diese Arbeit wird in einer Zeit $dt$ ausgeführt. Die dafür benötigte Leistung beträgt $P=dW/dt=\vec{F}\cdot\vec{v}$ [W].
In nachfolgenden Programm sind $F_x$, $F_y$ und $F_z$ die drei Komponenten der Kraft, $m$ die Masse und $t$ die Zeit. Ausgehend von den Anfangsbedingungen werden die Gleichungen für eine Gesamtzahl von $N_{Schritte}$ Schritten mit einer Schrittweite von $\Delta t$ integriert. Die Kraft kann in Abhängigkeit von den Positionen $x$, $y$, $z$, den Geschwindigkeiten $v_x$, $v_y$, $v_z$ und der Zeit $t$ angegeben werden.
$t$ [s] $x$ [m] $y$ [m] $z$ [m] $v_x$ [m/s] $v_y$ [m/s] $v_z$ [m/s] $F_x$ [N] $F_y$ [N] $F_z$ [N] $P$ [W] $E_{\text{kin}}$ [J] $W$ [J]
Fragen
Ein Objekt, das sich in drei Dimensionen bewegt, wird durch sechs Variablen beschrieben: $x$, $y$, $z$, $v_x$, $v_y$ und $v_z$. Wenn die Kraft auf das Objekt bekannt ist, dann kann die Bewegung durch sechs Differentialgleichungen erster Ordnung beschrieben werden,
$\large \frac{dx}{dt}=v_x$ $\large \frac{dv_x}{dt}=F_x(x,y,z,v_x,v_y,v_z,t)/m$
$\large \frac{dy}{dt}=v_y$ $\large \frac{dv_y}{dt}=F_y(x,y,z,v_x,v_y,v_z,t)/m$
$\large \frac{dz}{dt}=v_z$ $\large \frac{dv_z}{dt}=F_z(x,y,z,v_x,v_y,v_z,t)/m$
Diese Gleichungen können Schritt für Schritt numerisch gelöst werden. Bei bekannten Anfangsbedingungen für $\vec{x}(t_0)$ und $\vec{v}(t_0)$ ergibt sich als Näherung, für $\vec{x}$ und $\vec{v}$ kurze Zeit $\Delta t$ später,
$\large \vec{x}(t_0+\Delta t) \approx \frac{d\vec{v}}{dt}|_{t_0}\Delta t$ and $\large \vec{v}(t_0+\Delta t) \approx \vec{F}(x(t_0),v_x(t_0),t_0)\Delta t/m.$
Sobald eine Näherung für die Position und die Geschwindigkeit des Objekts zum Zeitpunkt $t_0 + \Delta t$ berechnet ist, kann die Näherung der Position und der Geschwindigkeit zum Zeitpunkt $t_0 + 2\Delta t$ berechnet werden. Durch wiederholen dieses Schrittes kann eine Tabelle mit Zeiten $(t)$, Positionen $(x,y,z)$, der Geschwindigkeit $(v_x,v_y,v_z)$ und der Kräfte $(F_x,F_y,F_z)$ erstellt werden. Andere Größen wie die kinetische Energie, die geleistete Arbeit und die Momentanleistung können berechnet werden. Die kinetische Energie beträgt $E_{\text{kin}}=\frac{v_x^2+v_y^2+v_z^2}{2m}$ [J]. Wenn eine Kraft ein Teilchen um eine Strecke $d\vec{r}$ bewegt, führt sie ein kleinen Teil $dW=\int\vec{F}\cdot d\vec{r}$ [J] Arbeit aus. Diese Arbeit wird in einer Zeit $dt$ ausgeführt. Die dafür benötigte Leistung beträgt $P=dW/dt=\vec{F}\cdot\vec{v}$ [W].
In nachfolgenden Programm sind $F_x$, $F_y$ und $F_z$ die drei Komponenten der Kraft, $m$ die Masse und $t$ die Zeit. Ausgehend von den Anfangsbedingungen werden die Gleichungen für eine Gesamtzahl von $N_{Schritte}$ Schritten mit einer Schrittweite von $\Delta t$ integriert. Die Kraft kann in Abhängigkeit von den Positionen $x$, $y$, $z$, den Geschwindigkeiten $v_x$, $v_y$, $v_z$ und der Zeit $t$ angegeben werden.
Wenn die Gesamtkraft auf ein Teilchen Null ist, bewegt es sich mit konstanter Geschwindigkeit geradlinig..
$\vec{F}=0$
Wenn die Gesamtkraft auf ein Teilchen konstant ist, bewegt es sich entlang einer Parabelbahn. Wenn die einzige Kraft auf ein Teilchen die Schwerkraft ist, dann ist die Kraft konstant und die Bewegung ist parabolisch.
$\vec{F}= -mg \hat{z}$
Hier ist $g = 9,807$ m/s² ist die Erdbeschleunigung an der Erdoberfläche.
Die Gesamtkraft auf ein Teilchen, welches senkrecht unter dem Einfluss der Gravitationskraft fällt und dabei eine lineare Reibungskraft erfährt, ist:
$F_z=-mg-bv_z,$
wobei $m$ die Masse, $g=9.81$ m/s² und $b$ der Reibungskoeffizient ist. Schließlich gleicht die Widerstandskraft die Gravitationskraft aus und das Teilchen bewegt sich mit einer Endgeschwindigkeit $v_T= mg/b$.
Wenn die Gesamtkraft linear proportional zur Verschiebung aus der Gleichgewichtsposition ist, ist die Bewegung harmonisch.
$\vec{F}= -kx \,\hat{x}$.
Hier ist $k$ die effektive Federkonstante.
Kreisbewegung ist eine Kombination aus harmonischer Bewegung in $x-$Richtung und harmonischer Bewegung in $y-$Richtung. Die Amplituden der Bewegung in $x-$- und $y-$-Richtung müssen gleich sein und die Anfangsbedingungen müssen so eingestellt sein, dass die maximale Verschiebung in $x-$-Richtung mit der Nullverschiebung in $y-Richtung zusammenfällt .
$\vec{F}= -kx \,\hat{x}-ky \,\hat{y}$.
Hier ist $k$ die effektive Federkonstante.
Satelliten Orbits
Ein Satellit umkreist die Erde. Die Schwerkraft auf den Satelliten ist,
$\large \vec{F} = -\frac{Gm_e m_s}{r^2} \hat{r}$,
mit $G= 6.6726 \times 10^{-11}$ N m²/kg² der Gravitationskonstanten, $m_e = 5.97219 \times 10^{24}$ kg der Masse der Erde, $m_s$ der Masse des Satellits, und $\vec{r}$ der Position des Satelliten, gemessen vom Mittelpunkt der Erde.
Solange die Erde weitaus schwerer ist, als der Satellit, hat die Masse des Satelliten keinen Einfluss. Die Masse des Satelliten taucht im Ausdruck für die Beschleunigung nicht auf. Einzig die Anfangsposition und Geschwindigkeit bestimmen die Bahn des Satelliten.
Die Anfangsbedingungen zur Zeit $t=0$ sind:
$x=$ m $y=$ m $z=$ m
$v_x=$ m/s $v_y=$ m/s $v_z=$ m/s $m_s=$ kg
Wenn der Orbit unter 6400000 m fällt, stürzt der Satellit auf die Erde. Es gibt eine Vielzahl an Umlaufbahnen, wie geosynchrone Umlaufbahnen, geostationäre Umlaufbahnen, Satellitenorbits, elliptische Umlaufbahnen, und Friedhofsorbits. Der Unterschied ist nur abhängig von den Anfangsbedingungen des Satelliten. Für die Berechnung der Umlaufbahn des Satelliten sollte ein großer Zeitschritt gewählt werden.
schiefe Ebene
Ein Block der Masse $m$ wird auf eine schiefe Ebene gestellt, die einen Winkel von $\theta$ Grad zur Horizontalen bildet.
Die Gravitationskraft $-mg\hat{y}$ lässt sich in eine Kraft senkrecht zur schiefen Ebene und eine Kraft parallel zur schiefen Ebene zerlegen.
Die Kraft parallel zur Ebene, $|F_{\parallel}|=mg\sin\theta$, lässt den Block entlang der Ebene gleiten. In Bezug auf ihre $x$- und $y$-Komponenten ist die Kraft entlang der Ebene,
$\vec{F}_{\parallel} = -mg\sin\theta\cos\theta \hat{x} - mg\sin\theta\sin\theta \hat{y}.$
Dies ist eine konstante Kraft; sie ist unabhängig von der Position, der Geschwindigkeit und der Zeit. Wir erwarten eine parabelförmige Bewegung der Masse und Formeln f&uulm;r den Ortsvektor und den Geschwindigkeitsvektor finden Sie auf der Seite Konstante Kraft = Parabelbewegung . Es ist natürlich auch möglich, diese Kraft und die Anfangsbedingungen in einen numerischen Differentialgleichungslöser zu stellen. Die Kraft und die Masse werden bei jeder Einstellung der Schieberegler in den Differentiallöser unten geladen. Dem Differentialgleichungslöser kann eine Widerstandskraft hinzugefügt werden, indem der Kraft ein Term hinzugefügt wird, der der Geschwindigkeit entgegengesetzt ist. Wenn diese Widerstandskraft linear $\vec{F}_{drag} = -a\vec{v}$ ist, dann ist es möglich Formeln zur Beschreibung der resultierenden Bewegung zu finden . Wenn die Widerstandskraft jedoch nichtlinear ist, wie z. B. $\vec{F}_{drag} = -a|\vec{v}|\vec{v}$, können keine einfachen Formeln und kein numerischer Differentialgleichungslöser gefunden werden ist die beste Option.
keine Widerstandskraft: $\vec{F} = -mg\sin\theta\cos\theta \,\hat{x} - mg\sin\theta\sin\theta \,\hat{y}$
lineare Widerstandskraft: $\vec{F} = -mg\sin\theta\cos\theta\, \hat{x} - mg\sin\theta\sin\theta \,\hat{y}-a\vec{v}$
nichtlineare Widerstandskraft: $\vec{F} = -mg\sin\theta\cos\theta \,\hat{x} - mg\sin\theta\sin\theta \,\hat{y}-a|\vec{v}|\vec{v}$
Luftreibung
Ein Stein der Masse 100 g wird mit einer Anfangsgeschwindigkeit $|\vec{v}|$ in einem Winkel $\theta$ von der Horizontalen geworfen. Der Stein wird aus seiner Position $\vec{r}=0$ gelöst. Wenn die Luftreibung vernachlässigt wird, lautet die Gleichung, die die Bewegung des Balls beschreibt,
$\vec{r} = vt\cos\theta \,\hat{x}+ \left( vt\sin\theta -\frac{gt^2}{2}\right)\,\hat{y}\,\text{m},$
wobei $g = 9,81$ m/s² ist die Erdbeschleunigung und $t$ wird in Sekunden gemessen. Der Stein folgt einer Parabel und kehrt zu der Höhe zurück, aus der er geworfen wurde, $y=0$, in einem Abstand von $x$ von wo er geworfen wurde. Für eine gegebene Anfangsgeschwindigkeit ist der Winkel, der den Abstand $x$ maximiert, $\theta = \frac{\pi}{4}\,\text{rad} = $ 45°.
Bezieht man in das Problem eine lineare Reibungskraft $\vec{F}_{drag} = -b\vec{v}$ ein, so ändert sich der optimale Winkel. Das folgende Form berechnet die Flugbahn des Steins für verschiedene Winkel, Reibungskonstanten und Anfangsgeschwindigkeiten.
$\vec{F} = m\frac{d^2\vec{r}}{dt^2} = -b\vec{v}-mg \,\hat{z}$
$b=$ N s/m
Die Anfangsbedingungen zum Zeitpunkt $t=0$ sind $\vec{r}=0$
$|\vec{v}| =$ m/s $\theta=$ deg
Ein Ball wird in den Wind geworfen
Ein Ball der Masse $m$ wird geworfen und erfährt einen Reibungswiderstand durch die Bewegung durch ein Gas oder eine Flüssigkeit. Die Kräfte, die auf den Ball wirken sind die Schwerkraft $-mg\hat{z}$ und die Reibungskraft. Weht ein Wind, kann dieser beschrieben werden durch,
$ \vec{F}_{fric} = -a(\vec{v}-\vec{v}_{\text{wind}}) - b(\vec{v}-\vec{v}_{\text{wind}})|(\vec{v}-\vec{v}_{\text{wind}})|,$
mit $a$ und $b$ Konstanten und $\vec{v}_{\text{wind}}$ der Windgeschwindgkeit, die von Ort und Zeit abhängig sein kann. Für niedrige Reynolds-Zahl , dominiert der lineare Term $-a(\vec{v}-\vec{v}_{\text{wind}})$ meistens, für eine hohe Reynolds-Zahl hingegen, dominiert der quadratische Term $- b(\vec{v}-\vec{v}_{\text{wind}})|(\vec{v}-\vec{v}_{\text{wind}})|$.
$\vec{F}= m\frac{d^2\vec{r}}{dt^2} = -a(\vec{v}-\vec{v}_{\text{wind}}) - b(\vec{v}-\vec{v}_{\text{wind}})|(\vec{v}-\vec{v}_{\text{wind}})|-mg\,\hat{z}$
$m=$ kg $a=$ N s/m $b=$ N s²/m²
Die drei Komponenten des Windvektors können Funktionen des Ortes und der Zeit sein.
$v_{\text{wind},x}=$ m/s $v_{\text{wind},y}=$ m/s $v_{\text{wind},z}=$ m/s
Die Anfangsbedingungen zur Zeit $t=0$ sind,
$x=$ m $y=$ m $z=$ m
$v_x=$ m/s $v_y=$ m/s $v_z=$ m/s
Raketenstart
Eine Modellrakete der Masse $m$ wird gestartet, mit einem Motor, welcher für 3 Sekunden eine aufwärts gerichtete Kraft $F_{\text{thrust}}$ erzeugt. Diese Kraft kann mathematisch beschrieben werden als $\vec{F}=F_{\text{thrust}}H(3-t)\hat{z}$ mit $H(x)$ der Heaviside Stufen Funktion . Da die Rakete Treibstoff verbraucht, nimmt ihre Masse ab. In dem Beispiel unten, verliert die Rakete die Hälfte ihrer Masse in den drei Sekunden ihrer Beschleunigung und behält dann eine konstante Masse, $m=0.1(2-H(3-t)t/3-H(t-3))$. Weitere Kräfte auf die Rakete wirken sind die Gravitation $-mg\hat{z}$ und eine Reibungskraft, die beschrieben werden kann durch,
$\large \vec{F}_{fric} = -a(\vec{v}-\vec{v}_{\text{wind}}) - b(\vec{v}-\vec{v}_{\text{wind}})|(\vec{v}-\vec{v}_{\text{wind}})|,$
mit $a$ und $b$ Konstanten und $\vec{v}_{\text{wind}}$ der Geschwindigkeit, die von Ort und Zeit abhängen kann.
$\vec{F}=0.1(2-H(3-t)t/3-H(t-3))\frac{d^2\vec{r}}{dt^2} = F_{\text{thrust}}H(3-t)\hat{z}-a(\vec{v}-\vec{v}_{\text{wind}}) - b(\vec{v}-\vec{v}_{\text{wind}})|(\vec{v}-\vec{v}_{\text{wind}})|-mg\,\hat{z}$
$m=$ kg $F_{\text{thrust}}=$ N $a=$ N s/m $b=$ N s²/m²
Die drei Komponenten des Windvektors können von Ort und Zeit abhängen .
$v_{\text{wind},x}=$ m/s $v_{\text{wind},y}=$ m/s $v_{\text{wind},z}=$ m/s
Die Anfangsbedingungen zur Zeit $t=0$ sind,
$x=$ m $y=$ m $z=$ m
$v_x=$ m/s $v_y=$ m/s $v_z=$ m/s
Ein springender Ball
Ein Ball der Masse $m$ wird geworfen und springt auf den Boden. Die auf diesen Ball wirkenden Kräfte sind die Schwerkraft $-mg\hat{z}$, eine der Geschwindigkeit proportionale Widerstandskraft $-a\vec{v}$ und die Kraft, die der Boden beim Aufprall ausübt. Der Boden kann als elastisches Material beschrieben werden, das mit einer Kraft $-kzH(-z)\hat{z}$ nach oben drückt. Hier ist $k$ die elastische Konstante und die Heaviside-Funktion $H(-z)$ sorgt dafür, dass diese Kraft nur dann wirkt, wenn die Ball $z=0$ unterschreitet.
Je größer die elastische Konstante $k$ ist, desto weniger sinkt der Ball unter $z=0$ und die Simulation prallt von einem harten Boden ab. Bei der Wahl des Zeitschritts $\Delta t$ ist jedoch Vorsicht geboten. Es sollte viel kleiner als $2\pi\sqrt{m/k}$ sein, damit der Solver die Abpraller vom Boden richtig beschreiben kann. Einige Differentialgleichungslöser passen den Zeitschritt automatisch an, um Probleme wie dieses zu berücksichtigen.
$\vec{F} = m\frac{d^2\vec{r}}{dt^2} = -a\vec{v}-mg\hat{z} -kzH(-z)\hat{z}$
$m=$ kg $a=$ N s/m $k=$ N/m
Die Anfangsbedingungen zur Zeit $t=0$ sind,
$x=$ m $y=$ m $z=$ m
$v_x=$ m/s $v_y=$ m/s $v_z=$ m/s
Curveball: Magnus effect
Ein sich drehender Ball, der durch die Luft geworfen wird, folgt aufgrund des Magnus-Effekts einer gekrümmten Bahn. Die Magnuskraft auf eine glatte Kugel ist,
$\vec{F}_{\text{Magnus}} = \frac{4}{3} \pi\rho r^3(\vec{\omega}\times\vec{v}),$
wobei $r$ der Radius der Ball, $\rho$ die Dichte der Luft, $\vec{\omega}$ die Kreisfrequenz der Drehung und $\vec{v}$ die Geschwindigkeit von ist der Ball. Die Luftwiderstandskraft wird angenähert als
$\vec{F}_{\text{drag}} = -\frac{\pi}{2} r^2 \rho C_d|\vec{v}|\vec{v},$
wobei $C_d$ die Luftwiderstandsbeiwert ist.
$\vec{F}= m\frac{d^2\vec{r}}{dt^2} = -\frac{\pi}{2} r^2 \rho C_d|\vec{v}|\vec{v}+ \frac{4}{3} \pi\rho r^3(\vec{\omega}\times\vec{v})-mg\,\hat{z}$
$m=$ kg $r=$ m
$\rho=$ kg/m³ $C_d=$
Die drei Komponenten des Vektors, der die Drehung beschreibt:
$\vec{\omega}=$ $\hat{x}$ + $\hat{y}$ + $\hat{z}$ rad/s
Im Allgemeinen ist der Luftwiderstandsbeiwert eine Funktion der Geschwindigkeit und kann in Form von vx, vy und vz ausgedrückt werden. Die Rotationsgeschwindigkeit nimmt mit der Zeit aufgrund von Reibung ab. Sie kann als Funktion der Zeit ausgedrückt werden.
Die Anfangsbedingungen zum Zeitpunkt $t=0$ sind,
$x=$ m $y=$ m $z=$ m
$v_x=$ m/s $v_y=$ m/s $v_z=$ m/s
Konstantes elektrisches Feld
Die Kraft auf ein Teilchen mit der Ladung $q$ in einem konstanten elektrischen Feld $\vec{E}$ beträgt:
$\vec{F} = q\vec{E}$.
Ein geladenes Teilchen folgt einer parabolischen Flugbahn in einem konstanten elektrischen Feld.
$m=$ kg $q=$ C
$\vec{E}=$ $\hat{x}$ + $\hat{y}$ + $\hat{z}$ V/m
Elektronen sind so leicht, dass ein kurzer Zeitschritt $\Delta t$ gewählt werden muss.
Coulomb-Potenzial
Ein geladenes Teilchen, das sich in einem Coulomb-Potenzial bewegt, folgt einer Umlaufbahn, die einer Satellitenumlaufbahn ähnelt. Die Kraft ist,
$\vec{F} = \frac{qQ}{4\pi \epsilon_0}\frac{\vec{r}-\vec{r}_0}{|\vec{r}-\vec{r}_0|^3} \hspace{1cm}\text{[N]}$.
Dabei ist $\vec{r}_0$ die Position der festen Ladung $Q$, die das Coulomb-Potenzial erzeugt, während $q$ die Ladung des beweglichen Teilchens ist. Der Einfachheit halber nehmen wir $\vec{r}_0=0$.
$m=$ kg $q=$ C kg $Q=$ C
Potenzial durch Ladungsleitung
Stellen Sie sich ein geladenes Teilchen vor, das sich im Potential bewegt, das durch eine Ladungslinie in $z-$Richtung verursacht wird. Das elektrische Feld ist,
$\vec{E}(\vec{r}) = \frac{\lambda}{2\pi\epsilon_0 (x^2+y^2)}\left(x\,\hat{x}+y\,\hat{y}\right)\quad\text{V/m},$
wobei $\lambda$ die lineare Ladungsdichte ist. Die Kraft auf ein Teilchen mit der Ladung $q$ und der Masse $m$ beträgt:
$\vec{F} = q\vec{E}\quad\text{N}.$
$m=$ kg $q=$ C kg $\lambda=$ [C/m]
Elektrische Ablenkung von Elektronen
Elektronen werden durch eine Spannung $V_x$ in Richtung einer positiv geladenen Platte beschleunigt. Einige der Elektronen passieren ein kleines Loch in der Platte und bilden einen Elektronenstrahl, der sich in einen Bereich bewegt, in dem ein elektrisches Feld aufgebaut wird, indem eine Spannung $V_y$ zwischen zwei Metallplatten angelegt wird, die einen Abstand $d$ voneinander haben.
Die Elektronen werden aus dem Ruhezustand mit einer Spannung von $V_x=5000$ V auf die positiv geladene Platte in einem Abstand von 5 cm beschleunigt. Die Elektronen werden abgelenkt, wenn sie das elektrische Feld zwischen den Platten passieren, wobei $V_y = 60$ V. Die Kraft beträgt in diesem Fall:
$\vec{F}=m\frac{d^2\vec{r}}{dt^2} = qE_xH(0.05-x)\,\hat{x} +qE_yH(0.1-|x-0.25|)\,\hat{y}$.
Eine Darstellung der elektrischen Ablenkung finden Sie unter Elektrische Ablenkung eines Elektronenstrahls .
Ein geladenes Teilchen in einem konstanten Magnetfeld
Die Kraft auf ein Teilchen mit der Ladung $q$ in einem Magnetfeld $\vec{B}$ beträgt:
$\vec{F} =q\vec{v}\times\vec{B}.$
Wenn die Geschwindigkeit $\vec{v}$ parallel zum Magnetfeld ist, ist die Kraft Null und das Teilchen bewegt sich geradlinig weiter. Liegt die Anfangsgeschwindigkeit nur in der Ebene senkrecht zum Magnetfeld, bewegt sich das Teilchen auf einem Kreis. Wenn die Anfangsgeschwindigkeit eine Komponente in der Ebene senkrecht zum Magnetfeld und eine Komponente parallel zum Feld aufweist, folgt das Teilchen einer Helix und dreht sich spiralförmig um die Richtung des Magnetfelds.
$m=$ kg $q=$ C
$\vec{B}=$ $\hat{x} + $ $\hat{y} + $ $\hat{z}$ [T]
Magnetische Ablenkung eines Elektronenstrahls
Ein Magnetfeld, das senkrecht zur Geschwindigkeit der Elektronen in einem Elektronenstrahl angelegt wird, beugt die Elektronen in eine Richtung, die sowohl zur Geschwindigkeit als auch zum Magnetfeld senkrecht ist. Im Diagramm unten werden Elektronen durch eine Spannung $V_x= 5000$ V in Richtung einer positiv geladenen Platte beschleunigt. Einige der Elektronen passieren ein kleines Loch in der Platte und bilden einen Elektronenstrahl, der sich zu einem kreisförmigen Bereich bewegt, in dem ein Magnetfeld $\vec{B} = 0,0025\,\hat{z}$ T von einem Elektromagneten erzeugt wird. Durch die Lorentzkraft werden die Elektronen abgelenkt und sie folgen einer Kreisbahn im Magnetfeld. Der Strahl wandert dann weiter und trifft auf einen Bildschirm, der 15 cm rechts von der Mitte des Magneten liegt.
Eine Simulation der magnetischen Ablenkung finden Sie unter Magnetische Ablenkung eines Elektronenstrahls .
J. J. Thomsons Experiment zur Bestimmung des Ladungs-zu-Masse-Verhältnisses von Elektronen
J. J. Thomson führte Experimente durch, um zu zeigen, dass Atome aus subatomaren Teilchen mit positiven und negativen Ladungen bestehen. Er stellte fest, dass die negativ geladenen Teilchen (Elektronen) viel leichter waren als die positiv geladenen Teilchen. Das Ladungs-zu-Masse-Verhältnis der Elektronen wurde gemessen, indem die Elektronen durch eine Spannung $V_x$ in Richtung einer positiv geladenen Platte beschleunigt wurden. Einige der Elektronen passieren ein kleines Loch in der Platte und bilden einen Elektronenstrahl, der in einen Bereich wandert, in dem ein elektrisches Feld und ein magnetisches Feld vorhanden waren. Die Elektronen wurden in $x$-Richtung beschleunigt, das elektrische Feld war in $y$-Richtung und das magnetische Feld war in $z$-Richtung.
Die Kraft auf ein Elektron aufgrund des elektrischen Feldes beträgt $\vec{F}_{\text{elec}}=-|e|E_y\hat{y}$ und die Kraft auf ein Elektron aufgrund des magnetischen Feldes Feld ist $\vec{F}_{\text{magn}}=|e|v_xB\hat{y}$. Die Geschwindigkeit der Elektronen kann bestimmt werden, da bei der Beschleunigung der Elektronen in $x$-Richtung die potentielle Energie $|e|V_x$ in kinetische Energie umgewandelt wird $\frac{1}{2}mv_x^2= |e|V_x$. Daher $v_x=\sqrt{2|e|V_x/m}$. Wenn das Experiment so eingestellt wird, dass sich die elektrischen und magnetischen Kräfte gegenseitig aufheben $|e|v_xB\hat{y}-|e|E_y\hat{y}=0$, dann ist das Ladungs-zu-Masse-Verhältnis gegeben von,
$\large \frac{e}{m}=\frac{V_y^2}{2V_x\mu_0^2n^2I^2d^2}.$
Die folgende Berechnung gilt für $V_x = 5000$ V, $V_y = 60$ V, $\vec{B}= 0,00025\,\hat{z}$ T. Anfangs haben die Elektronen eine Geschwindigkeit von Null. Die Metallplatten sind 10 cm lang und haben einen Abstand von $d=$ 2 cm. Weitere Details finden Sie unter J. J. Thomsons Experiment zur Bestimmung des Ladungs-zu-Masse-Verhältnisses von Elektronen .
Ein Elektron bewegt sich im Magnetfeld, das durch einen stromdurchflossenen Draht erzeugt wird.
Ein Magnetfeld wird durch einen Strom erzeugt, der durch einen langen geraden Draht in der +$z$-Richtung fließt
$\vec{B}(\vec{r})=\frac{\mu_0I}{2\pi(x^2+y^2)}\,\hat{z}\times (x\hat{x} + y\hat{y})\,\,\text{[T]}$
wobei $I$ der Strom ist. Für Details siehe: Das Magnetfeld um einen unendlich langen geraden Draht .
Die Kraft auf ein Teilchen mit der Ladung $q$ in diesem Magnetfeld beträgt:
$\vec{F} =q\vec{v}\times\vec{B}.$
Die Werte für an werden zunächst in das Formular geladen. Das Elektron bewegt sich auf einer komplizierten Bahn um den Draht. Beachten Sie, dass die geleistete Arbeit (letzte Spalte in der Tabelle unten) Null bleibt (innerhalb des numerischen Fehlers), da die Kraft immer senkrecht zur Geschwindigkeit ist.
$m=$ kg $q=$ C $I=$ A
Ein geladenes Teilchen in elektrischen und magnetischen Feldern
Wenn sich ein geladenes Teilchen mit einer Ladung $q$ und einer Masse $m$ in einem elektrischen Feld $\vec{E}$ und einem magnetischen Feld $\vec{B}$ bewegt, beträgt die Kraft auf das Teilchen:
$ \vec{F} = q(\vec{E} + \vec{v}\times \vec{B}).$
Ausgeschrieben in ihren drei Komponenten ist die Lorentzkraft:
$F_x = q(E_x+v_yB_z-v_zB_y)$, $F_y = q(E_y+v_zB_x-v_xB_z)$, $F_z = q(E_z+v_xB_y-v_yB_x).$
Die Komponenten der elektrischen und magnetischen Felder können in den Textfeldern unten angegeben werden. Diese können dann in den allgemeinen 3D-Bewegungsdifferentialgleichungslöser geladen werden, um die Flugbahn des geladenen Teilchens zu berechnen.
$m=$ kg $q=$ C
$E_x(x,y,z,t)=$ V/m $E_y(x,y,z,t)=$ V/m $E_z(x,y,z,t)=$ V/m
$B_x(x,y,z,t)=$ T $B_y(x,y,z,t)=$ T $B_z(x,y,z,t)=$ T
Wenn ein konstantes, gleichmäßiges elektrisches Feld senkrecht zu einem konstanten, gleichmäßigen Magnetfeld verläuft, verläuft die durchschnittliche Geschwindigkeit des geladenen Teilchens in der Richtung senkrecht sowohl zum elektrischen Feld als auch zum magnetischen Feld. Beachten Sie, dass sich Elektronen sehr schnell bewegen, daher muss ein kurzer Zeitschritt $\Delta t$ gewählt werden.
Masse-nichtlineare Feder System
Ein numerischer Differentialgleichungsrechner kann auch verwendet werden, wenn die Feder nonlinear ist.
Für eine lineare Feder ist die Kraft proportional zur Auslenkung der Feder, $F=-kx$.
Bei einer harten Feder, wächst die Kraft schneller als linear, wenn die Feder ausgelenkt wird.
Ein Beispiel für eine harte Feder ist $F=-kx|x|$.
Bei einer weichen Feder, wächst die Kraft langsamer als linear, wenn die Feder ausgelenkt wird.
Ein Beispiel für eine weiche Feder ist $F=-k\sin(x)$.
In der darunterliegenden Simulation wird die Feder von der Gleichgewichtsposition ausgelenkt, und die Masse $m$ aus dem Ruhezustand losgelassen.
Wenn Reibung vernachlässigt wird, dann oszilliert die Masse um die Gleichgewichtsposition der Feder.
Für nichtlineare Federn, ist die Frequenz der Oszillation von der Amplitude Abhängig.
Bei harten Federn erhöht sich die Frequenz mit der Amplitude, und bei weiche Federn verringert sich die Frequenz mit der Amplitude.
harte Feder $F=-kx|x|$
lineare Feder $F=-kx$
weiche Feder $F=-k\sin(x)$
$m=$ kg $k=$ N/m
Resonanz: Gedämpftes angetriebenes Masse-Feder-System
Eine Masse $m$ ist an einer linearen Feder mit der Federkonstante $k$ befestigt. Auf die Masse wirkt eine Reibungskraft in entgegengesetze Richtung der Geschwindigkeit $F_{\text{drag}}=-bv_x$ (sei $b$ der Reibungskoeffizient). Das System wird durch eine periodische Antriebskraft $F_0\cos(\Omega t)$ angetrieben. Die folgende Differentialgleichung beschreibt die Bewegung:
$\large m\frac{d^2x}{dt^2}+b\frac{dx}{dt}+kx=F_0\cos(\Omega t).$
Die Resonanzfrequenz ist $\omega=\sqrt{k/m-b^2/4m^2}=$ 0.943 rad/s. Die maximale Auswirkung folgt aus der Gleichheit der Treiber- und Resonanzfrequenz. Die Amplitude der daraus resultierenden Schwingungen (im Gleichgewichtszustand) ist $F_0/\sqrt{(k-m\Omega^2)^2+\Omega^2b^2}=$ 0.447 m.
Weitere Details finden Sie unter: Resonanz eines gedämpften harmonischen Oszillators .
Chaotische Lösungen für das angetriebene Pendel
Die Differentialgleichung, die die Bewegung eines angetriebenen Pendels beschreibt, lautet [Fitzpatrick 2006] a>,
$\large \frac{d^2\theta}{dt^2}+\frac{1}{q}\frac{d\theta}{dt}+\sin(\theta)=A\cos(\Omega t),$
wobei $q$ die Dämpfung beschreibt, $A$ das Drehmoment misst, das verwendet wird, um das Pendel mit der Frequenz $\Omega$ anzutreiben, und $\theta$ der von der Vertikalen gemessene Winkel ist. Bei $\theta=0$ hängt das Pendel nach unten und bei $\theta=\pi$ steht das Pendel auf. Es ist bekannt, dass diese Gleichung chaotische Lösungen aufweist.
Weitere Details finden Sie unter Chaotische Lösungen für das angetriebene Pendel .
Zwei gekoppelte Masse-Feder Schwinger
Die beiden Differentialgleichungen, die die reibungsfreien Bewegungen $x(t)$ beziehungsweise $y(t)$ zweier, über eine Feder gekoppelte Massen (der Masse $m$) beschreibt, lauten
$\large m \frac{d^2 x}{dt^2} = -kx + k_{12} (y-x)$
$\large m \frac{d^2 y}{dt^2} = -ky + k_{12} (x-y)$
wenn jede Masse zusätzlich mit einer Feder der Federkonstante $k$ gegenüber der koppelnden Feder $k_{12}$ verbunden ist.
Was geschieht, wenn beide Massen aus der Ruhelage mit der gleichen Anfangsauslenkung bzw. entgegengesetzter Anfangsauslenkungen zu schwingen beginnen?
Weitere Details finden Sie unter Schwebung bei zwei gekoppelten Schwingern .
Swinging spring
Die schwingende Feder oder das elastische Pendel ist eine Masse, die an einem Ende mit einem elastischen Stab verbunden ist, der sich dehnen, aber nicht biegen kann. Das andere Ende des Stabes ist an einem festen Punkt befestigt, so dass die Masse frei schwingen kann. Dieses System kann eine komplexe Bewegung ausführen, und für einige Anfangsbedingungen führt es eine chaotische Bewegung aus. Die Gleichungen lauten:
$\large m\frac{d^2x}{dt^2}=-m\omega^2\left(\frac{|\vec{r}| -l}{|\vec{r}|}\right)x$
$\large m\frac{d^2y}{dt^2}=-m\omega^2\left(\frac{|\vec{r}| -l}{|\vec{r}|}\right)y$
$\large m\frac{d^2z}{dt^2}=-m\omega^2\left(\frac{|\vec{r}| -l}{|\vec{r}|}\right)z-mg$
wobei $m$ die Masse des Stabes, $l$ die ungestreckte Länge des Stabes, $g = 9.807$ die Erdbeschleunigung an der Erdoberfläche, $\omega^2=k/m$, und $|\vec{r}| =\sqrt{x^2+y^2+z^2}$ ist.