振動解析の数値計算

第六回：オイラーの公式と振動の表現方法

今回は第二回でお話ししたオイラーの公式について解説していきます。公式そのものは式①で表されます。

振動解析のコラムとしては余談となりますが、この式に x=π を代入して式を変形すると式②となり、これは「世界一美しい等式」と呼ばれるものになります。

この式の証明についてはテイラー展開によるものが有名ですが、本コラムでは証明そのものではなく、この式の工学的な有用性について解説してみようと思います。
結論を先に述べると、この式の最も大きな意味は「振動現象を三角関数ではなく複素指数関数で表現できる」ことにあります。本コラムでは「第三回：１自由度減衰系の自由振動」から複素指数関数を用いて解を表現していますが、これは三角関数での表現が非常に煩雑になるからです。

第三回で登場したこの運動方程式を満たす解は複素指数関数に限定されると述べましたが、実際には複素数を用いず、下記の式④という実関数解も考えられます。(１自由度不減衰系の解に時間とともに減衰するe^-αtを掛けることで解を表現しています)

この式の１階微分、２階微分をそれぞれ求めると

また元の運動方程式を不減衰固有各振動数ω₀と減衰比ζを用いて変形しておき、

先ほどの１階微分、２階微分を代入して、[⋯] cos ⁡ωt+[⋯] sin ⁡ωt=0の形に整理すると、この式が任意の時刻tで成立するためには、cos ⁡ωtとsin ⁡ωtは独立な関数であるから、それぞれの係数が０である必要があり

式⑪を式⑨に代入すると

ωは実数なのでζ≦1であり、また

となるため、C₂の符号で吸収可能であり、-ω_dは省略します。
以上より、三角関数を用いて解を表現すると以下のようになります。

この式に初期条件とζに具体的な値を入れてグラフ化すると、第三回で表示したグラフと同じものになります。
ここまで記述したことで既にお分かりかと思いますが、三角関数のみで表現しようとすると計算が非常に複雑になるのに対して、複素指数関数を用いると計算を圧倒的に簡略化することができます。またζ＞1の場合（過減衰=振動しない現象）は平方根の内部が負の値となり、実関数である三角関数での表現が出来ません。そのため、実務上では「振動現象＝複素指数関数」が前提となっています。
ここで動画をご覧ください。この動画はフーリエ変換とオイラーの公式を用いて、複雑な振動波形を複素指数関数で表現する流れを表したものになります。

全体の流れとしては、
①複雑な振動波形をフーリエ変換によって三角関数の級数に分解
②一つの周波数に注目し、加法定理を用いて振幅と位相を算出
③らせん運動を経由して、単振動を円運動に変換
④円運動を複素数平面に投影し、円運動を複素指数関数で表現
となります。
動画で示した通りですが、複雑な振動波形から複素指数関数に至るまで情報が失われておらず、可逆的な変換であることが分かると思います。すなわち振動現象の表現方法として紹介したa cosωt+b sinωt 、Acos(ωt+θ) 、Ae^i(ωt+θ) はいずれも等価であり、複素指数関数を用いる理由は「運動方程式の計算（微積分の処理）が容易になるから」につきます。

また式③のような2階線形常微分方程式の形で表される物理現象は物体の振動問題に限らず、
・電気回路（RLC直列回路）
・音響（ヘルムホルツ共鳴器）
・制御工学（２次遅れ系）
など、さまざまな分野に共通して登場し、いずれも複素指数関数によって表現されます。このようにオイラーの公式そのものは三角関数と複素指数関数をつなぐシンプルなものですが、今日の工学において、振動・波動・信号といった周期現象を扱う上で不可欠な基礎となっています。

今回はオイラーの公式について説明しました。次回は多自由度系への拡張と振動解析について説明する予定です。