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.

     

$k=$ 20 [N/m]

$m=$ 0.4 [kg]

$x(t_0)=$ 0.02 [m]

Die Periode der Schwingungen hängen von der Anfangsamplitude $x(t_0)$ ab.

 Numerical 2nd order differential equation solver 

$ \large \frac{dx}{dt}=$

$v_x$

$ \large a_x=\frac{F_x}{m}=\frac{dv_x}{dt}=$

Intitial conditions:

$x(t_0)=$

$\Delta t=$

$v_x(t_0)=$

$N_{steps}$

$t_0=$

Plot:

vs.

 

 $t$       $x$      $v_x$