Gekoppelte Schwingung, Frequenzspaltung

Hi @all,

Wahrscheinlich haben einige von Euch von der Aufspaltung der
Energieniveaus gelesen, wenn 2 Atome näher aneinander rücken.
Es steht auch in vielen Bücher nachzulesen, dass dieser Effekt
sich "wie eine Frequenzaufspaltung bei gekoppelten Schwingungen"
vorzustellen sei. Leider halten sich alle diese Bücher mit dem
genauen Modell dieser gekoppelten Schwingung zurück.
Es ist sehr ärgerlich, wenn die Bücher nur so mit verschiedenen
Modellen um sich werfen ohne sie anständig zu erklären.

Aber gut, hier habe ich versucht so eine Kopplung zu realisieren.
Bei einfachem Feder-Masse System ist bekanntlich die Kreisfrequenz
Omega = sqrt(k/m).
Was ich aber machen will, ist Wand-Feder-Masse-Feder-Masse-Feder-Wand
oder kurz wkmkmw. Unten ist der Code und ein Bild der Simulation
habe ich angehängt.

1

model wkmkmw " |||/\/\/\/\O/\/\/\/\O/\/\/\/\/\/||| " 

2

  parameter Real m1(unit = "kg") = 1 "masse";

3

  parameter Real m2(unit = "kg") = 100 "masse";

4

  parameter Real k1(unit = "N/m") = 30 "federkonstante";

5

  parameter Real k2(unit = "N/m") = 30 "federkonstante";

6

  parameter Real k3(unit = "N/m") = 30 "federkonstante";

7

  parameter Real p(unit = "m") = 0 "ruhelage";

8

  parameter Real q(unit = "m") = 200 "ruhelage";

9

10

  Real x1(start = p, unit = "m") "ort";

11

  Real x2(start = q+100, unit = "m") "ort";

12

  Real v1(start = 0, unit = "m/s") "geschwindigkeit";

13

  Real v2(start = 0, unit = "m/s") "geschwindigkeit";

14

15

equation 

16

  v1 = der(x1);

17

  v2 = der(x2);

18

  -m1*der(v1) - k1*(x1-p) + k2*(x2-q-(x1-p)) = 0;

19

  -k2*(x2-q-(x1-p)) - m2*der(v2) - k3*(x2-q) = 0;

20

end wkmkmw;

Die 2-te Masse ist riesieg, sie wird um ihre Ruhelage mit
einer kleinen Frequenz oszillieren. Die 1-te dagegen vibriert
deutlich stärker. Im Bild fällt auf, dass beide reine Sinus-
schwingunen sind. Was noch auffällt(nach mehreren Simulationen),
dass die Amplitude der ersten Schwingung gleich der anfänglichen
Auslenkung ist, Amplitude der 2-ten nur die Hälfte davon.

Man könnte im Prinzip, denke ich, mit den Sinusansätzen in die
wild aussehende DGL. Interessanter und wilder wird der Fall, wenn
die Massen gleich sind. Im nachfolgenden Post kommt ein Bild dazu.