2重振り子②-微小でない場合
2重振り子の振動②-微小でない場合
微小でない場合の2重振り子の振動
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_default_img.png)
上図のそれぞれのおもりはどちらも質量をとし、その重りをつないでいる糸の長さは
で曲がったりせずかつ重さは無視できるものとします。
軸に関しては下向きにとります。
こうした場合、ラグランジアンの式は次のような式になります。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_Lagrangian_functions_calculate_process_img1-1024x475-1.png)
途中の式では以下のような三角関数の公式、
![](https://diff-eq.com/wp-content/uploads/2024/06/trigonometric_formulas_img1.png)
を使っています。
についてのラグランジアンは、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_lagrangian_function_calculation_process_img1.png)
ここで三角関数の公式を使うと、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_lagrangian_function_calculation_process_img2.png)
ここで求まった上記式に対して以下の三角関数の公式を適用します。
![](https://diff-eq.com/wp-content/uploads/2024/06/trigonometric_formulas_img2.png)
そうすると以下のように求まります。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_lagrangian_function_calculation_process_img3.png)
この結果をについてのラグランジアンに使えば、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_lagrangian_function_calculation_process_img4.png)
次はに対しても同じようにすると、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_lagrangian_function_calculation_process_img5.png)
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_lagrangian_function_calculation_process_img6.png)
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_lagrangian_function_calculation_process_img7.png)
となるので以下のように求まります。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_lagrangian_function_calculation_process_img8.png)
上記の式に関してそれぞれ、
のそれぞれについて求めていきます。
まず上記式に関して、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_lagrangian_function_calculation_process_img9.png)
出てきたこの式を式に代入すると、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_lagrangian_function_calculation_process_img10-1024x74-1.png)
さらに計算していきます。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_lagrangian_function_calculation_process_img11_formula_correction20240412-1024x286-1.png)
これによりに関しては以下のように求まります。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_theta1_2dot_lagrangian_function_calculation_result_img.png)
また上記式導出においては三角関数におけるの性質である
を使っています。
次にを求めます。
まず、式において
について変形すると、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_theta1_2dot_lagrangian_function_calculation_img1.png)
この式を、先ほどの、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_theta_2dot_default_img.png)
の式の右辺のに代入し計算していきます。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_theta2_2dot_lagrangian_function_calculation_process_img1-1024x312-1.png)
これより、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_theta2_2dot_lagrangian_function_calculation_process_img2-1024x75-1.png)
この式の右辺をさらに以下のようにして計算していきます。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_theta2_2dot_lagrangian_function_calculation_process_img3-1024x303-1.png)
これによりは以下のように求まります。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_theta2_2dot_lagrangian_function_calculation_process_img4-1024x50-1.png)
まとめれば以下のようになります。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_theta1_2dot_theta2_2dot_lagrangian_function_calculation_result_img-1024x152-1.png)
出てきた式を見てわかるように、かなり複雑な運動を行うことが予想できると思います。
さらにおもりの重さの違いや糸の長さの違いなどが加わればもっと複雑な、なかばカオス的な状態を導くことになります。
近似式を用いた振動の考察
いったんここで、上式において次に示す近似式、
![](https://diff-eq.com/wp-content/uploads/2024/06/approximation_of_trigonometric_functions.png)
これらを代入して計算するとどうなるかを考察してみましょう。
まずより、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_theta1_2dot_approximate_value_solution_img1.png)
さらにより、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_theta2_2dot_approximate_value_solution_img1.png)
結果として次のような関係式が求まります。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_theta1_theta2_2dot_approximate_value_solution_img1.png)
さらに今度は上記関係式を連立させて計算していきます。
まず、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_theta1_theta2_2dot_approximate_value_solution_simultaneous_ode_img1.png)
次にの式に関してそれを2倍して辺々かけ合わせます。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_theta1_theta2_2dot_approximate_value_solution_simultaneous_ode_img2.png)
先ほどの複雑な2つの式が以下のように簡単な形式の連立微分方程式になります。
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum2_theta1_theta2_2dot_approximate_value_solution_simultaneous_ode_img3.png)
ここでとおいてみましょう。
すると、
![](https://diff-eq.com/wp-content/uploads/2024/06/double_pendulum1_lagrangian_function_calculation_result_img3-1.png)
これは2重振り子の振動①-微小な場合において求められた、
の連立式と同一なものであることがわかります。