Введение

Мы рассматриваем линейную цепочку N маятников связанных друг с другом упругими связями. Обозначим угол вращения [Maple Math] -го маятника как [Maple Math] и угловую скорость как [Maple Math] . На каждый маятник действует угловой момент из-за силы тяжести и угловой момент из-за двух упругих связей, который пропорционален разности углов вращений связанных маятников с коэффициентом [Maple Math] . При предположении, что все маятники имеют один момент инерции J, система уравнений движения, взятая из [1], имеет вид

[Maple Math]

где [Maple Math] есть угловые моменты от упругих связей и [Maple Math] есть гравитационный угловой момент.

Угловые моменты левой и правой упругих связей равны [Maple Math] и [Maple Math] , соответственно и следовательно

[Maple Math]

Угловой момент из-за силы тяжести может быть выражен как

>    Gamma2[i]:=-M*d*g*sin(phi[i]);

Gamma2[i] := -M*d*g*sin(phi[i])

где g есть гравитационная постоянная.

В виду выражений для Gamma1[i]  and Gamma2[i] , уравнения движения могут быть переписаны как

>    J*diff(phi[i](t),t$2)=kappa*(phi[i+1](t)-2*phi[i](t)+phi[i-1](t))-M*d*g*sin(phi[i](t));

J*diff(phi[i](t),`$`(t,2)) = kappa*(phi[i+1](t)-2*phi[i](t)+phi[i-1](t))-M*d*g*sin(phi[i](t))

Записанная выше система уравнений движения в континуальном пределе преобразуется к следующему уравнению в частных производных

>    J*diff(phi(x,t),t$2)=Kappa*diff(phi(x,t),x$2)-K[G]*sin(phi(x,t));

J*diff(phi(x,t),`$`(t,2)) = Kappa*diff(phi(x,t),`$`(x,2))-K[G]*sin(phi(x,t))

где K[G] = Md*g , K= kappa h ^2.

В новых пространственных и временных переменных, [Maple Math] , [Maple Math] , последнее уравнение приобретает стандартную форму уравнения синус-Гордон

[Maple Math]

Можно видеть, что система связанных маятников описывается дискретным аналогом уравнения синус-Гордон. Далее мы продемонстрируем некоторые решения уравнения синус-Гордон, используя описанную выше модель.