2重振り子の振動①-微小な場合
おもりを2つ吊るした振り子の微小振動の動きをラグランジアンを使って解析してみましょう。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_default_2024_img1-2.png)
糸の重さは無視できるものとします。
微小振動としているので、
![](https://diff-eq.com/wp-content/uploads/2024/06/approximation_of_trigonometric_functions_sin_tan_theta.png)
このような場合は水平方向のみを考えればよいです。
まず運動エネルギーは、
![](https://diff-eq.com/wp-content/uploads/2024/06/kinetic_energy_double_pendulum1_img1.png)
ポテンシャルエネルギーは、
![](https://diff-eq.com/wp-content/uploads/2024/06/potential_energy_double_pendulum1_img-1024x34-1.png)
とした場合
となってしまうので、以下のように級数展開します。
![](https://diff-eq.com/wp-content/uploads/2024/06/cos_theta_taylor_expansion_series_img.png)
より、
![](https://diff-eq.com/wp-content/uploads/2024/06/cos_theta_taylor_expansion_series_img2.png)
とします。
これを使えば、
![](https://diff-eq.com/wp-content/uploads/2024/06/potential_energy_double_pendulum1_img2.png)
よってラグランジアンは次のようになります。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_lagrangian_function_img1.png)
より、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_lagrangian_function_calculation_img2-1-1024x252-1.png)
に関しても同様にして、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_lagrangian_function_calculation_img3.png)
ここで、
として見やすくすると、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_lagrangian_function_calculation_result_img1.png)
さらにとおけば、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_lagrangian_function_calculation_result_img2.png)
2つのおもりの角振動数を、初期位相を
とし、解を次のように置きます。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_lagrangian_solution_theta1_theta2_img.png)
それぞれ2度微分すれば、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_lagrangian_solution_theta1_theta2_img2-1024x136-1.png)
これらを、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_lagrangian_function_calculation_result_img3.png)
に代入します。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_lagrangian_function_calculation_result_img4.png)
これより、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_lagrangian_function_calculation_result_img5.png)
上記の式においてが同時に
にならない条件を求めます。
その条件は永年方程式を使って次のようになります。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_huckle_method_img.png)
この行列式を計算します。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_huckle_method_calculation_process_img1.png)
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_huckle_method_calculation_process_img2.png)
ここでとすると、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_huckle_method_calculation_process_img3.png)
この上記方程式を解いていきます。
解の公式を使って求めると、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_huckle_method_quadratic_formula_calculation_process_img.png)
置き換えたをもとに戻せば、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_huckle_method_calculation_process_img4.png)
以下のように求まります。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_huckle_method_calculation_result_img.png)
これが微小振動における2重振り子の基準振動数になります。