線形多自由度系の振動のソースを表示
←
線形多自由度系の振動
ナビゲーションに移動
検索に移動
あなたには「このページの編集」を行う権限がありません。理由は以下の通りです:
この操作は、次のグループに属する利用者のみが実行できます:
登録利用者
。
このページのソースの閲覧やコピーができます。
[[振動工学]]における'''線形多自由度系の振動'''（せんけいたじゆうどけいのしんどう）は、[[線型性|線形]]な特性を持ち、さらに2以上の[[自由度]]を持つ系で起きる[[振動]]である。[[運動方程式]]は一般的に連立2階常微分方程式となり、[[行列]]および[[ベクトル]]で表現される。 線形多自由度系の振動では、固有モードという多自由度系特有の概念が現れ、自由度の数だけ固有モードと固有振動数の組が存在する。固有モードの直交性によって、[[減衰振動|減衰]]の無い系であれば固有モードごとの1自由度系の問題に帰着でき、振動解析を容易化できるのが特徴である。この手法を利用した振動解析手法は[[モード解析]]と呼ばれる。減衰のある系でも、比例粘性減衰という仮定を導入することによって、同様なことが可能となる。 モード解析手法は、振動実験結果から振動特性を[[システム同定|同定]]するのにも使われる。[[有限要素法]]による[[連続体力学|連続体]]の振動計算においても、線形多自由度系の理論にもとづくモード解析手法が強力な効果を発揮し、振動解析を容易にする。 ==背景== [[振動]]の問題を分類すると、あるいは振動現象を模した動力学モデルを分類すると、いくつかの視点が存在する{{Sfnm|末岡・金光・近藤|2000|1p=4|下郷・田島|2002|2p=4|横山・日野・芳村|2015|3p=3}}。 [[File:Mass spring damper.svg|thumb|190px|[[質点]]が[[ばね]]と[[ダッシュポット]]を介して基礎に固定されたモデル]] 物体の運動を表すのに必要な[[座標]]あるいは[[変位]]の数を、[[自由度]]という{{Sfnm|砂子田・伊藤・鄭・平元|2012|1p=4|平井・水田|2018|2p=156}}。自由度が1の系を'''1自由度系'''という{{Sfn|下郷・田島|2002|p=20}}。ある1つの[[質点]]が[[ばね]]と[[ダッシュポット]]を介して地面に固定され、質点が上下方向のみ動く例を考える。[[三角関数]]で表現される励振力がこの質点に加わるとき、この質点の[[運動方程式]]は次のような形で与えられる{{Sfn|末岡・金光・近藤|2000|pp=12–13}}。 {{NumBlk|:|<math> m \ddot{x} + c \dot{x} + k x = f_0 \cos \omega t </math>|{{EquationRef|1.1}}}} ここで、''m'' は[[質量]]、''c'' は[[粘性減衰係数]]、''k'' は[[ばね定数]]、''f''<sub>0</sub> は外力の[[振幅]]、''ω'' は外力の[[角振動数]]、''t'' は[[時間]]である{{Sfn|末岡・金光・近藤|2000|p=13}}。''x'' は静的釣り合い位置からの[[変位]]で、上部の "˙" は[[時間微分]] (''d''/''dt'') を示す{{Sfn|末岡・金光・近藤|2000|p=13}}。式{{EqNoteN|1.1}}は、振動問題を考える際の最も基礎的な式となる{{Sfn|末岡・金光・近藤|2000|pp=12–13}}。1自由度系の振動の問題は、振動現象を理解する上で不可欠な様々な概念を内包しており、1自由度系のモデルが振動問題を扱う基本モデルといえる{{Sfnm|末岡・金光・近藤|2000|1pp=4, 8|下郷・田島|2002|2p=20}}。 一方、自由度が2以上の系は'''多自由度系'''と呼ばれる{{Sfn|吉川・松井・石井|1987|p=31}}。振動問題をモデル化する際には1自由度や2自由度にモデル化すれば解析が容易であり、モデル化時の自由度はできるだけ最小限とするのが大切である{{Sfnm|吉川・松井・石井|1987|1p=46|砂子田・伊藤・鄭・平元|2012|2p=5}}。しかし、現実の構造物は複雑で、すべてを1自由度系として扱うことができない{{Sfn|砂子田・伊藤・鄭・平元|2012|p=57}}。実際の機械や建物などでは、1自由度系のモデルではその特性を説明しきれず、多自由度系としての取り扱いが必要になることも多い{{Sfnm|宮本・永野・藤谷・吉村|2014|1p=83|入江・小林|2006|2p=82}}。 また、系の構成要素が全て[[線型性|線形]]であるとき、その系を'''線形系'''という{{Sfn|末岡・金光・近藤|2000|p=5}}。系が線形であれば、出力と入力の間には[[重ね合わせの原理]]が成り立ち、出力は入力に単純に[[比例]]する{{Sfn|下郷・田島|2002|p=5}}。振動系が線形であれば、式{{EqNoteN|1.1}}のように慣性力、減衰力、復元力が、それぞれ[[加速度]]、[[速度]]、[[変位]]に比例する{{Sfnm|砂子田・伊藤・鄭・平元|2012|1p=5|藤田|2016|2p=15|末岡・金光・近藤|2000|3pp=12–13}}。比例係数が時間変化する場合を除き、線形系の振動問題については理論的手法がほぼ確立できている{{Sfn|下郷・田島|2002|p=5}}。 一方で、線形ではない系を'''非線形系'''という{{Sfn|末岡・金光・近藤|2000|p=5}}。摩擦、がたつき、大変位時の材料特性などが非線形の要因となる{{Sfn|日本機械学会（編）|2004|p=53}}。厳密に考えると、実際の機械や構造物はなんらかの非線形特性を持ち、実際の振動現象のほとんどは非線形系といえる{{Sfnm|砂子田・伊藤・鄭・平元|2012|1p=135|入江・小林|2006|2p=150}}。非線形系の厳密解を得ることはほとんどできず、一般的には非線形系の問題に対しては近似解法や数値計算に頼らざるを得ない{{Sfnm|末岡・金光・近藤|2000|1p=5|下郷・田島|2002|2p=5}}。しかし、線形系の振動として取り扱うことによって、十分に多くの問題を解決することもできる{{Sfn|入江・小林|2006|p=150}}。非線形系であっても、安定な平衡状態周りの微小変位運動であれば線形の理論が当てはまる{{Sfn|入江・小林|2006|p=150}}。特に線形系であれば、多自由度系であっても後述のように1自由度系の問題に帰着できるという顕著な特性がある{{Sfnm|末岡・金光・近藤|2000|1p=4|下郷・田島|2002|2p=5}}。 ==運動方程式== ===一般の多自由度系=== 多自由度系では多数の変数や係数を扱うため、[[行列]]と[[ベクトル]]を使って[[運動方程式]]を記す{{Sfn|背戸・丸山|2002|p=201}}。自由度が ''n'' の線形多自由度系の一般的な運動方程式は、次のように書き表される{{Sfnm|下郷・田島|2002|1p=82|日本機械学会（編）|2004|2pp=41–42|末岡・金光・近藤|2000|3pp=43–44|長松・長松|2018|4p=96|入江・小林|2006|5p=105}}。 {{NumBlk|:|<math> \boldsymbol{M \ddot{x}} + \boldsymbol{C \dot{x}} + \boldsymbol{K x} = \boldsymbol{f} </math>|{{EquationRef|2.1}}}} ここで、'''''ẍ''''', '''''ẋ''''', '''''x''''', '''''f''''' は ''n'' 次元縦ベクトル、'''''M''''', '''''C''''', '''''K''''' は ''n'' 次[[正方行列]]で、以下のように表される{{Sfnm|日本機械学会（編）|2004|1pp=41–42|入江・小林|2006|2pp=105–106}}。 :<math> \boldsymbol{\ddot{x}} = \begin{pmatrix} \ddot{x}_1 \\ \ddot{x}_2 \\ \vdots \\ \ddot{x}_n \end{pmatrix} </math>, <math> \boldsymbol{\dot{x}} = \begin{pmatrix} \dot{x}_1 \\ \dot{x}_2 \\ \vdots \\ \dot{x}_n \end{pmatrix} </math>, <math> \boldsymbol{x} = \begin{pmatrix} x_1 \\ x_2 \\ \vdots \\ x_n \end{pmatrix} </math>, <math> \boldsymbol{f} = \begin{pmatrix} f_1 \\ f_2 \\ \vdots \\ f_n \end{pmatrix} </math> {{NumBlk|:|<math> \boldsymbol{M} = \begin{pmatrix} m_{11} & m_{12} & \cdots & m_{1n} \\ m_{21} & m_{22} & \cdots & m_{2n} \\ \vdots & \vdots & \ddots & \vdots \\ m_{n1} & m_{n2} & \cdots & m_{nn} \end{pmatrix} </math>|{{EquationRef|2.2}}}} {{NumBlk|:|<math> \boldsymbol{C} = \begin{pmatrix} c_{11} & c_{12} & \cdots & c_{1n} \\ c_{21} & c_{22} & \cdots & c_{2n} \\ \vdots & \vdots & \ddots & \vdots \\ c_{n1} & c_{n2} & \cdots & c_{nn} \end{pmatrix} </math>|{{EquationRef|2.3}}}} {{NumBlk|:|<math> \boldsymbol{K} = \begin{pmatrix} k_{11} & k_{12} & \cdots & k_{1n} \\ k_{21} & k_{22} & \cdots & k_{2n} \\ \vdots & \vdots & \ddots & \vdots \\ k_{n1} & k_{n2} & \cdots & k_{nn} \end{pmatrix} </math>|{{EquationRef|2.4}}}} '''''x''''' は'''変位ベクトル'''と呼ばれる{{Sfn|長松・長松|2018|p=70}}。その成分の ''x'' は、静的な釣り合いの位置を原点とした各自由度の[[変位]]あるいは[[一般化座標]]を表している{{Sfnm|長松・長松|2018|1p=70|日本機械学会（編）|2004|2p=40}}。'''''ẋ''''' は'''速度ベクトル'''、'''''ẍ''''' は'''加速度ベクトル'''と呼ばれる{{Sfn|長松・長松|2018|pp=70, 91}}。成分の ''ẋ'' は ''x'' の時間1階微分すなわち各自由度の[[速度]]を意味し、''ẍ'' は ''x'' の時間2階微分すなわち[[各自由度]]の加速度を意味する{{Sfnm|宮本・永野・藤谷・吉村|2014|1p=23|長松・長松|2018|2pp=70, 91|宮本・永野・藤谷・吉村|2014|3p=23}}。 '''''f''''' は'''外力ベクトル'''と呼ばれる{{Sfn|末岡・金光・近藤|2000|p=43}}。成分の ''f'' は各一般化座標に対応して作用する外力で、一般的には時間の関数である{{Sfn|日本機械学会（編）|2004|p=42}}。 '''''M''''' は'''質量行列'''や'''慣性行列'''、'''質量マトリックス'''や'''慣性マトリックス'''と呼ばれる{{Sfnm|下郷・田島|2002|1p=81|宮本・永野・藤谷・吉村|2014|2p=83|吉川・松井・石井|1987|3p=96}}。'''''K''''' は'''剛性行列'''や'''剛性マトリックス'''と呼ばれる{{Sfnm|末岡・金光・近藤|2000|1p=43|宮本・永野・藤谷・吉村|2014|2p=83}}。'''''C''''' は'''減衰行列'''や'''減衰マトリックス'''と呼ばれる{{Sfnm|末岡・金光・近藤|2000|1p=43|吉川・松井・石井|1987|2p=97}}。各成分の ''m'', ''c'', ''k'' は質量、粘性減衰係数、[[剛性]]（ばね定数）を表す{{Sfnm|長松・長松|2018|1pp=70, 91|下郷・田島|2002|2pp=81–82}}。ここでの ''m'' は'''慣性係数'''ともいい、いわゆる[[質量]]だけでなく[[慣性モーメント]]なども含む{{Sfn|下郷・田島|2002|pp=88–89}}。また、''k'' は'''復元係数'''ともいい、通常のばね定数の他に回転ばね定数なども含む{{Sfn|下郷・田島|2002|pp=88–89}}。 線形多自由度系の一般基礎式{{EqNoteN|2.1}}は、線形1自由度系の基礎式{{EqNoteN|1.1}}と形式は同じで、変位と外力がベクトルに、質量とばね定数と粘性減衰係数が行列に置き換わった式となる{{Sfn|末岡・金光・近藤|2000|p=44}}。ただし、対象の系の規模が大きくなり、複雑なものとなると、[[運動の第2法則|ニュートンの運動法則]]から運動方程式を導出するのは容易ではない{{Sfn|末岡・金光・近藤|2000|p=207}}。実際に多自由度系の運動方程式を立てる際は、エネルギーのスカラー量から形式的に運動方程式を導ける[[ラグランジュの運動方程式]]が便利で、間違いを犯しにくく、多用されている{{Sfnm|吉川・松井・石井|1987|1p=91|背戸・丸山|2002|2p=201|藤田|2016|3p=119}}。 ===2自由度系の例=== [[File:Two-story shear building model.svg|thumb|300px|2層構造物の揺れのモデリング]] 建築物の例では、複数の階を持つ多層構造物が多自由度系の問題となる{{Sfn|宮本・永野・藤谷・吉村|2014|p=83}}。柱と床から構成される[[ラーメン構造]]の建物が振動する場合を考える。柱は床に比べて柔らかいので、構造物が揺れるとき各階の床は水平方向に揺れ、柱は水平方向のばねとして働くと見なせる{{Sfnm|宮本・永野・藤谷・吉村|2014|1p=27|平井・水田|2018|2p=39}}。これは、水平方向のみに動く2つの質点をばねで連結した2自由度系モデルと等価となる{{Sfn|平井・水田|2018|pp=83–84}}。2層構造物における、1層目の床質量を ''m''<sub>1</sub>、2層目の床質量を ''m''<sub>2</sub>、1層目の水平変位を ''x''<sub>1</sub>、2層目の水平変位を ''x''<sub>2</sub>、基礎と1層目の間の等価ばね定数を ''k''<sub>1</sub>、1層目と2層目の間の等価ばね定数を ''k''<sub>2</sub> とすると、このモデルの運動方程式は {{NumBlk|:| <math>\begin{cases} m_1 \ddot{x}_1 + (k_1 + k_2)x_1 - k_2 x_2 = 0 \\ m_2 \ddot{x}_2 - k_2 x_1 +k_2 x_2 = 0 \end{cases}</math>|{{EquationRef|2.5}}}} となる{{Sfn|平井・水田|2018|p=87}}。行列表示すれば、 {{NumBlk|:|<math> \begin{pmatrix} m_{1} & 0 \\ 0 & m_{2} \end{pmatrix} \begin{pmatrix} \ddot{x}_{1} \\ \ddot{x}_{2} \end{pmatrix} + \begin{pmatrix} k_{1}+k_{2} & -k_{2} \\ -k_{2} & k_{2} \end{pmatrix} \begin{pmatrix} x_{1} \\ x_{2} \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \end{pmatrix} </math>|{{EquationRef|2.6}}}} となり、この式では、<math> \begin{pmatrix} m_{1} & 0 \\ 0 & m_{2} \end{pmatrix} </math>が質量行列、<math> \begin{pmatrix} k_{1}+k_{2} & -k_{2} \\ -k_{2} & k_{2} \end{pmatrix} </math>が剛性行列である{{Sfn|平井・水田|2018|p=87}}。 [[File:Coupled vibration model of vertical for pitching motion of automobile.svg|thumb|250px|自動車の上下振動とピッチング振動の簡易モデル]] 並進運動と回転運動が組み合わさった2自由度系の例として、自動車の上下振動とピッチング振動の簡易モデルがある{{Sfn|藤田|2016|p=98}}。自動車を側面から見て、上下方向変位とピッチング回転だけの動きを考える。前位側のタイヤとサスペンションを1つのばねと見なして、後位側のタイヤとサスペンションも同様に1つのばねと見なし、[[剛体]]を前位と後位を2つのばねが支えているモデルを考える{{Sfn|入江・小林|2006|p=88}}。前位側と後位側のばねのばね定数をそれぞれ ''k''<sub>1</sub>、''k''<sub>2</sub> とする。剛体の重心位置から前位側のばねまでの距離を ''l''<sub>1</sub>、前位側のばねまでの距離を ''l''<sub>2</sub> とする。剛体の質量を ''m''、重心位置周りのピッチング方向の[[慣性モーメント]]を ''I<sub>G</sub>'' とする。この剛体の重心の上下運動 ''x'' とその周りのピッチング運動 ''θ'' の運動方程式は次のようになる{{Sfn|横山・日野・芳村|2015|p=83}}。 {{NumBlk|:|<math> \begin{cases} m \ddot{x} + (k_1 + k_2)x + (k_2 l_2 - k_1 l_1)\theta = 0 \\ I_G \ddot{\theta} + (k_2 l_2 - k_1 l_1)x + (k_2 l^{2}_{2} + k_1 l^{2}_{1})\theta = 0 \end{cases} </math>|{{EquationRef|2.7}}}} ===連成項=== 自動車の上下振動とピッチング振動の簡易モデルである式{{EqNoteN|2.7}}において、もし ''k''<sub>1</sub>''l''<sub>1</sub> = ''k''<sub>2</sub>''l''<sub>2</sub> であれば、''x'' と ''θ'' は互いに独立した運動方程式となる{{Sfn|藤田|2016|p=100}}。この条件では運動方程式は {{NumBlk|:|<math> \begin{cases} m \ddot{x} + (k_1 + k_2)x = 0 \\ I_G \ddot{\theta} + (k_2 l^{2}_{2} + k_1 l^{2}_{1})\theta = 0 \end{cases} </math>|{{EquationRef|2.8}}}} となり、2つの式はそれぞれ単独で解くことができる{{Sfn|砂子田・伊藤・鄭・平元|2012|p=68}}。すなわち、''x'' の振動と ''θ'' の振動が、互いに干渉すること無しに独立して起こる状態になっている{{Sfn|横山・日野・芳村|2015|p=85}}。このような状態を'''非連成'''という{{Sfn|藤田|2016|p=100}}。運動方程式が単独では解けない連立方程式になっているとき、方程式は'''連成'''しているといい、式{{EqNoteN|2.7}}の (''k''<sub>1</sub>''l''<sub>1</sub> − ''k''<sub>2</sub>''l''<sub>2</sub>)''x'' と (''k''<sub>1</sub>''l''<sub>1</sub> − ''k''<sub>2</sub>''l''<sub>2</sub>)''θ'' は'''連成項'''と呼ばれる{{Sfn|砂子田・伊藤・鄭・平元|2012|p=58}}<ref name="小寺・矢野2009">{{Cite book ja-jp |author = 小寺 忠・矢野 澄雄 |title = 例題で学ぶ機械振動学 |publisher = 森北出版 |url = https://www.morikita.co.jp/books/book/1441 |edition = 第1版 |year = 2009 |isbn = 978-4-627-66711-2 }} p. 55</ref>。''k''<sub>1</sub>''l''<sub>1</sub> = ''k''<sub>2</sub>''l''<sub>2</sub> のときに、これら連成項が 0 になり、''x'' と ''θ'' が独立した運動になる{{Sfn|入江・小林|2006|p=89}}。連成状態にある系を'''連成系'''、逆に各自由度の運動が完全に独立している系を'''非連成系'''と呼ぶ{{Sfn|末岡・金光・近藤|2000|p=43}}。 一般に、質量行列 '''''M'''''（式{{EqNoteN|2.2}}）、減衰行列 '''''C'''''（式{{EqNoteN|2.3}}）、剛性行列 '''''K'''''（式{{EqNoteN|2.4}}）における非対角成分の存在は、振動が連成していることを示している{{Sfn|吉川・松井・石井|1987|p=97}}。もし系が非連成系であれば、これら全ての行列は、非対角成分が全て 0 の[[対角行列]]となる{{Sfn|末岡・金光・近藤|2000|p=44}}。特に、質量行列が非対角成分によって振動が連成していることを'''動連成'''や'''動的連成'''といい、質量行列の非対角成分を'''動連成項'''という<ref name="振動工学ハンドブック">{{Cite book ja-jp |editor = 振動工学ハンドブック編集委員会 |title = 振動工学ハンドブック |publisher = 養賢堂 |edition = 第3版 |year = 1985 }} p. 38</ref>{{Sfn|横山・日野・芳村|2015|pp=100–101}}。また、剛性行列が非対角行列であれば'''静連成'''や'''静的連成'''といい、剛性行列の非対角成分を'''静連成項'''という<ref name="振動工学ハンドブック"/>{{Sfn|横山・日野・芳村|2015|pp=100–101}}。多自由度系の問題は概して複雑で、個々のパラメータが結果に与える影響の見通しを立てることが難しい{{Sfn|下郷・田島|2002|p=85}}。多自由度系では一般的に連成が存在し、これが多自由度系の解析を難しいものにしている{{Sfn|背戸・丸山|2002|p=201}}。動力学的な設計を行う上では、系を非連成化して影響をわかりやすくすることが有効となる{{Sfn|下郷・田島|2002|p=85}}。 ==不減衰自由振動== ===振動数方程式=== 継続的な外力が作用せず、外乱だけが与えられて起こる振動を[[自由振動]]という{{Sfn|藤田|2016|p=13}}。自由振動の振動は、系自体が持つ特性によって決まり、系の動特性を知る上で重要な振動の形態である{{Sfn|長松・長松|2018|pp=4–5}}。自由振動は減衰要素が存在しない場合と存在する場合に分かれる{{Sfn|吉川・松井・石井|1987|p=30}}。減衰要素と励振力が存在しない場合の式{{EqNoteN|2.1}}は {{NumBlk|:|<math> \boldsymbol{M \ddot{x}} + \boldsymbol{K x} = \mathbf{0} </math>|{{EquationRef|3.1}}}} となる{{Sfn|藤田|2016|p=124}}。ここで、'''0''' は下記のような[[ゼロベクトル]]である。 :<math> \mathbf{0} = \begin{pmatrix} 0 \\ 0 \\ \vdots \\ 0 \end{pmatrix} </math> さらに、'''''M''''' と '''''K''''' は[[行列の定値性|正定値]]であると仮定する。すなわち、'''''M''''' と '''''K''''' は成分が全て[[実数]]の[[対称行列]]で、それらの二次形式は常に正となる{{Sfn|安田|2012|p=119}}。実際に、質量行列と剛性行列が正定値行列であることは線形多自由度系の一般的な特徴であり、多くの振動系でこの仮定は満たされる{{Sfnm|末岡・金光・近藤|2000|1p=46|安田|2012|2p=119}}。この条件を満たす ''M'' と ''K'' を前提にすれば、振動系の具体的な構成に依存しない一般性の高い議論を展開できる{{Sfn|末岡・金光・近藤|2000|p=47}}。 式{{EqNoteN|3.1}}は定数係数の[[線型微分方程式|線形常微分方程式]]であるため、その解法に従って解を {{NumBlk|:|<math> \boldsymbol{x} = \boldsymbol{u} e^{\lambda t} </math>|{{EquationRef|3.2}}}} という形式で仮定できる{{Sfn|安田|2012|p=120}}。ここで ''e'' は[[ネイピア数]]、''λ'' は未知定数、 '''''u''''' は ''n'' 個の未知定数 ''u'' から成る下記のような縦ベクトルである。 :<math> \boldsymbol{u} = \begin{pmatrix} u_1 \\ u_2 \\ \vdots \\ u_n \end{pmatrix} </math> 式{{EqNoteN|3.1}}に対して仮定として与えられる解には、他に、最初から単振動を仮定して[[複素指数関数]]や[[三角関数]]の形式もある{{Sfnm|横山・日野・芳村|2015|1p=108|日本機械学会（編）|2004|2p=41}}。 式{{EqNoteN|3.2}}を式{{EqNoteN|3.1}}に代入して整理すると下記のような式になる{{Sfn|藤田|2016|p=124}}。 {{NumBlk|:|<math> (\lambda^2 \boldsymbol{M} + \boldsymbol{K}) \boldsymbol{u} e^{\lambda t} = \mathbf{0} </math>|{{EquationRef|3.3}}}} この式が恒等的に成り立つには下記の条件が満たされる必要がある{{Sfn|藤田|2016|p=124}}。 {{NumBlk|:|<math> (\lambda^2 \boldsymbol{M} + \boldsymbol{K}) \boldsymbol{u} = \mathbf{0} </math>|{{EquationRef|3.4}}}} 数学的には、この形式の方程式は[[一般化固有値問題]]として知られ、''λ''<sup>2</sup> は[[固有値]]、'''''u''''' は固有ベクトルと呼ばれる{{Sfnm|末岡・金光・近藤|2000|1p=48|横山・日野・芳村|2015|2p=109}}。'''''u''''' = '''0''' であれば式{{EqNoteN|3.4}}の条件は満たされるが、これは最初の釣り合いの位置からそのまま静止しているだけ状態を意味する解である{{Sfn|砂子田・伊藤・鄭・平元|2012|p=82}}。よって、ここでは '''''u''''' = '''0''' は興味の対象外で、それ以外の式{{EqNoteN|3.4}}を満たす ''λ'' と '''''u''''' について知りたい{{Sfn|砂子田・伊藤・鄭・平元|2012|p=82}}。'''''u''''' ≠ '''0''' で、なおかつ式{{EqNoteN|3.4}}が満たされる条件は、''U'' の係数行列 ''λ''<sup>2</sup> '''''M''''' + '''''K''''' の[[逆行列]]が存在しないことである{{Sfn|安田|2012|p=120}}。したがって、係数行列の[[行列式]]が 0 であればこの条件が満たされる{{Sfn|安田|2012|p=120}}。したがって、''λ'' が、 {{NumBlk|:|<math> \mbox{det}(\lambda^2 \boldsymbol{M} + \boldsymbol{K}) = 0 </math>|{{EquationRef|3.5}}}} を満たすとき、式{{EqNoteN|3.4}}が '''''u''''' = '''0''' 以外の解を持つ{{Sfn|安田|2012|p=120}}。ここで、det( ) は行列式を表す。式{{EqNoteN|3.5}}の行列式を展開すると、 ''λ''<sup>2</sup> についての ''n'' 次多項式になる{{Sfn|安田|2012|p=120}}。したがって、原理的には ''λ''<sup>2</sup> の値を求めることができる{{Sfn|安田|2012|p=120}}。ただし、この多項式を解析的に解くことができるのはせいぜい2自由度あるいは3自由度までで、それ以上の自由度の系になると[[数値解析]]で固有値を計算する{{Sfn|宮本・永野・藤谷・吉村|2014}}。[[ヤコビ法 (固有値問題)|ヤコビ法]]などを使って数値計算するときは、[[コレスキー分解]]を使い、一般化固有値問題形式の{{EqNoteN|3.4}}を標準的な[[固有値問題]]の形へ変換する{{Sfn|藤田|2016|pp=205–207}}。 上記のとおり、'''''M''''' と '''''K''''' は正定値であると仮定した。このとき、式{{EqNoteN|3.5}}の解は全て負の実数となる{{Sfn|安田|2012|p=120}}（解を複素指数関数と仮定した場合は全て正の実数{{Sfn|末岡・金光・近藤|2000|p=109}}）。したがって、''n'' 個の ''λ''<sup>2</sup> の解を {{NumBlk|:|<math> \lambda^2 = -\omega_1^2 ,\ -\omega_2^2 ,\ \cdots,\ -\omega_n^2 </math>|{{EquationRef|3.6}}}} とおくことができ、平方根を取って 2''n'' 個の''λ'' の値 {{NumBlk|:|<math> \lambda = \pm j \omega_1 ,\ \pm j \omega_2 ,\ \cdots,\ \pm j \omega_n </math>|{{EquationRef|3.7}}}} が得られる{{Sfn|安田|2012|p=121}}。ここで ''j'' は[[虚数単位]]である。''n'' 個の ''ω'' は'''固有角振動数'''または'''固有円振動数'''と呼ばれ、値が小さいものから順に'''1次'''、'''2次'''、…、'''''n'' 次'''の固有角振動数と呼ぶ{{Sfnm|横山・日野・芳村|2015|1p=109|宮本・永野・藤谷・吉村|2014|2p=87|吉川・松井・石井|1987|3pp=109–110}}。特に、最も値が小さい1次の固有角振動数は'''基本振動数'''と呼ばれる{{Sfn|長松・長松|2018|p=75}}。式{{EqNoteN|3.5}}は角振動数を求める式であるため'''振動数方程式'''と呼ばれる{{Sfn|長松・長松|2018|p=75}}。あるいは、式{{EqNoteN|3.5}}は '''''M''''' や '''''K''''' といった系の特性によって構成される式であることから'''特性方程式'''とも呼ぶ{{Sfn|長松・長松|2018|p=75}}。 ===固有モード=== 得られた固有値あるいは固有角振動数の値を式{{EqNoteN|3.4}}に代入すれば、'''0''' 以外の '''''u''''' の解が得られる{{Sfnm|安田|2012|1p=121|下郷・田島|2002|2p=89}}。ここで得られる '''''u''''' の各成分は一意な値を持たず、定まるのは各成分の互いの比 ''u''<sub>1</sub> : ''u''<sub>2</sub> : … : ''u<sub>n</sub>'' だけである{{Sfnm|安田|2012|1p=121|砂子田・伊藤・鄭・平元|2012|2p=109}}。一つの固有値あるいは固有角振動数に対応して一つの '''''u''''' が定まり、'''''u''''' は ''n'' 個存在する{{Sfn|末岡・金光・近藤|2000|pp=48–49}}。このような固有値と固有ベクトルの組は固有ペアと呼ばれる{{Sfn|末岡・金光・近藤|2000|p=48}}。''r'' 次の固有角振動数 ''ω<sub>r</sub>'' に対応する '''''u''''' を '''''u'''<sub>r</sub>'' を表現すれば、 '''''u'''<sub>r</sub>'' は {{NumBlk|:|<math> \boldsymbol{u}_r = \begin{pmatrix} u_{1r} \\ u_{2r} \\ \vdots \\ u_{nr} \end{pmatrix} </math>|{{EquationRef|3.8}}}} というベクトルである{{Sfn|下郷・田島|2002|p=89}}。ただし、上記のとおり、各成分 ''u''<sub>1''r''</sub>, ''u''<sub>2''r''</sub>, …, ''u<sub>nr</sub>'' は互いの比の値を表している{{Sfn|下郷・田島|2002|p=89}}。 [[File:Nomal modes of Two-story shear building model.svg|thumb|320px|式{{EqNoteN|2.5}}で表現される2層構造物振動モデルの固有モードの例{{Sfn|平井・水田|2018|pp=90–91}}。''k''<sub>1</sub> = ''k''<sub>2</sub>, ''m''<sub>1</sub> = ''m''<sub>2</sub> の場合を示しており、各固有モードの成分の内の一番大きな絶対値が1となるようにして値を定めている。]] 式{{EqNoteN|3.8}}で表されるベクトルが式{{EqNoteN|3.4}}における固有ベクトルであり、振動工学では'''固有モード'''、'''振動モード'''、'''固有振動モード'''、'''基準振動モード'''、'''モードベクトル'''などと呼ぶ{{Sfnm|下郷・田島|2002|1p=89|横山・日野・芳村|2015|2p=110|吉川・松井・石井|1987|3p=110}}。''n'' 自由度系には ''n'' 個の固有角振動数があり、固有角振動数それぞれに対応する形で ''n'' 個の固有モードが存在している{{Sfn|吉川・松井・石井|1987|p=110}}。各自由度の振幅比を決める固有モードは、固有角振動数が「振動の速さ」を表しているのに対して、「振動の形」を表していると言える{{Sfn|長松・長松|2018|p=76}}。自由度の数の分だけ固有角振動数が存在し、それら固有角振動数に対応して固有モードが存在していることが、線形多自由度系の特有な性質といえる{{Sfn|末岡・金光・近藤|2000|p=42}}。 各々の固有モードを、対応する固有角振動数が小さい順に1次固有モード、2次固有モード、…、''r'' 次固有モード、…、''n'' 次固有モードと呼ぶ{{Sfn|藤田|2016|p=126}}。特に、基本振動数（1次の固有角振動数）に対応する固有モードは'''基本モード'''と呼ばれる{{Sfn|長松・長松|2018|p=75}}。固有モード '''''u'''<sub>r</sub>'' を次数が低い順に並べて作る、下記のような行列を'''モード行列'''や'''モードマトリックス'''という{{Sfnm|日本機械学会（編）|2004|1p=92|入江・小林|2006|2p=106|吉川・松井・石井|1987|3pp=111–112}}。 {{NumBlk|:|<math> \boldsymbol{U} = (\boldsymbol{u}_1,\ \boldsymbol{u}_2,\ \cdots,\ \boldsymbol{u}_n) = \begin{pmatrix} u_{11} & u_{21} & \cdots & u_{1n} \\ u_{21} & u_{22} & \cdots & u_{2n} \\ \vdots & \vdots & \ddots & \vdots \\ u_{n1} & u_{n2} & \cdots & u_{nn} \end{pmatrix} </math>|{{EquationRef|3.9}}}} ===固有モードの直交性と正規化=== 多自由度系の固有モードには'''直交性'''という重要な性質がある{{Sfnm|平井・水田|2018|1p=93|藤田|2016|2p=128}}。''r'' 次の固有モード '''''u'''<sub>r</sub>'' と ''s'' 次の固有モード '''''u'''<sub>s</sub>'' について考える。ここで、''r'' と ''s'' は任意だが、''r'' ≠ ''s'' である。'''''M''''' と '''''K''''' は上記のとおり対称行列とする。このとき、 '''''u'''<sub>r</sub>'', '''''u'''<sub>s</sub>'', '''''M''''', '''''K''''' には次のような関係がある{{Sfn|横山・日野・芳村|2015|p=110}}。 {{NumBlk|:|<math> \boldsymbol{u}_r^\top \boldsymbol{M} \boldsymbol{u}_s = 0 </math>|{{EquationRef|3.10}}}} {{NumBlk|:|<math> \boldsymbol{u}_r^\top \boldsymbol{K} \boldsymbol{u}_s = 0 </math>|{{EquationRef|3.11}}}} ここで、'''''u'''<sub>r</sub>''<sup>⊤</sup> は '''''u'''<sub>r</sub>'' の[[転置行列]]を表している。これらの式は '''''u'''<sub>r</sub>'' と '''''Mu'''<sub>s</sub>'' の[[内積]]および '''''u'''<sub>r</sub>'' と '''''Mu'''<sub>s</sub>'' の内積が零であることを示しており、'''''M''''' と '''''K''''' に関して '''''u'''<sub>r</sub>'' と '''''u'''<sub>s</sub>'' が[[直交]]していることを意味している{{Sfn|日本機械学会（編）|2004|p=41}}。これが質量行列および剛性行列を介した固有モードの直交性である{{Sfnm|横山・日野・芳村|2015|1pp=110–111|日本機械学会（編）|2004|2p=41}}。別の言い回しでは '''''u'''<sub>r</sub>'' と '''''u'''<sub>s</sub>'' が一般直交性を有しているともいう{{Sfn|長松・長松|2018|p=79}}。 固有モードの直交性の物理的な意味は、次数の異なる固有モードの振動の間で[[力学的エネルギー]]の移動が全く起きないことを表している{{Sfnm|長松・長松|2018|1p=81|末岡・金光・近藤|2000|2p=61}}。ある固有モードの振動において、[[運動エネルギー]]と復元力による[[ポテンシャルエネルギー]]は時間的に変化しているが、それらの和の力学的エネルギーは時間に依らず一定に保たれている{{Sfn|安田|2012|p=130}}。そのため、自由振動中にある固有モードの振動が他の固有モードに移り変わったり、他の固有モードが新たに誘起されたりすることはない{{Sfn|長松・長松|2018|p=81}}。 一方で、''r'' = ''s''、すなわち同じ次数同士の固有モードの場合は、上記の式でも左辺は 0 とはならず、ある定数となる{{Sfn|砂子田・伊藤・鄭・平元|2012|p=85}}。これらの定数を ''M<sub>r</sub>'' と ''K<sub>r</sub>'' と表せば、 {{NumBlk|:|<math> \boldsymbol{u}_r^\top \boldsymbol{M} \boldsymbol{u}_r = M_r </math>|{{EquationRef|3.12}}}} {{NumBlk|:|<math> \boldsymbol{u}_r^\top \boldsymbol{K} \boldsymbol{u}_r = K_r </math>|{{EquationRef|3.13}}}} と表される{{Sfn|砂子田・伊藤・鄭・平元|2012|p=85}}。固有モード '''''u'''<sub>r</sub>'' は成分間の比だけを持ち、定まった値を持たないので、''M<sub>r</sub>'' と ''K<sub>r</sub>'' の値もこの段階では定まらない{{Sfn|長松・長松|2018|p=82}}。固有モードの絶対値が決まった後に、''M<sub>r</sub>'' と ''K<sub>r</sub>'' の値も定まる{{Sfn|長松・長松|2018|p=82}}。 定数 ''M<sub>r</sub>'' と ''K<sub>r</sub>'' は正の値であり、それぞれを（''r'' 次の）'''モード質量'''および'''モード剛性'''という{{Sfn|横山・日野・芳村|2015|p=111}}。モード行列を使って式{{EqNoteN|3.12}}と式{{EqNoteN|3.13}}を表せば、 {{NumBlk|:|<math> \boldsymbol{U}^\top \boldsymbol{M} \boldsymbol{U} = \begin{pmatrix} M_1 & 0 & \cdots & 0 \\ 0 & M_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & M_n \end{pmatrix} </math>|{{EquationRef|3.14}}}} {{NumBlk|:|<math> \boldsymbol{U}^\top \boldsymbol{K} \boldsymbol{U} = \begin{pmatrix} K_1 & 0 & \cdots & 0 \\ 0 & K_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & K_n \end{pmatrix} </math>|{{EquationRef|3.15}}}} となる{{Sfn|入江・小林|2006|p=106}}。モード質量とモード剛性は、同じ次数の固有角振動数 ''ω<sub>r</sub>'' と下記のような関係がある{{Sfn|砂子田・伊藤・鄭・平元|2012|p=85}}。 {{NumBlk|:|<math> \omega_r = \sqrt{ \frac{K_r}{M_r}} </math>|{{EquationRef|3.16}}}} 固有モードは比の関係を定めるベクトルであった。絶対的な大きさを持たない固有モードに大きさを定める方法として、次のような方法がある。 #下記のように、式{{EqNoteN|3.12}}の右辺 ''M<sub>r</sub>'' の値が 1 となるように定める{{Sfnm|砂子田・伊藤・鄭・平元|2012|1p=86|吉川・松井・石井|1987|2p=110}}。ここで、'''''ū'''<sub>r</sub>'' は、等式を満たすに大きさが決定された固有モードを意味している。 #:{{NumBlk||<math> \boldsymbol{\overline{u}}_{r}^\top \boldsymbol{M} \boldsymbol{\overline{u}}_{r} = 1 </math>|{{EquationRef|3.17}}}} #下記のように、同じ次数の固有モード同士の[[内積]]が 1 となるように定める{{Sfnm|藤田|2016|1p=125|砂子田・伊藤・鄭・平元|2012|2p=86}}。 #:{{NumBlk||<math> \boldsymbol{\overline{u}}_{r}^\top \boldsymbol{\overline{u}}_{r} = 1 </math>|{{EquationRef|3.18}}}} #固有モードの成分の内、[[絶対値]]が最大のものを 1 とおいて定める{{Sfn|背戸・丸山|2002|p=209}}。 #1番目の自由度に対応する固有モード成分 (''u''<sub>1''r''</sub>) を 1 とおいて定める{{Sfnm|藤田|2016|1p=125|吉川・松井・石井|1987|2p=110}}。 以上のように固有モードの大きさを一意に定めることを'''正規化'''と呼び、正規化された固有モード '''''ū'''<sub>r</sub>'' を'''正規固有モード'''や'''正規化モード'''と呼ぶ{{Sfnm|砂子田・伊藤・鄭・平元|2012|1p=87|横山・日野・芳村|2015|2pp=109–110|平井・水田|2018|3p=91}}。1番目の手法による正規固有モードは、特に'''質量正規固有モード'''や'''''M''-正規固有モード'''という名で呼ばれる{{Sfnm|長松・長松|2018|1p=84|末岡・金光・近藤|2000|2p=50}}。式{{EqNoteN|3.17}}を満たすようにしておくと利点が多く、質量正規固有モードはよく使われる{{Sfnm|吉川・松井・石井|1987|1p=110|平井・水田|2018|2p=91}}。正規固有モード '''''ū'''<sub>r</sub>'' を並べて作る下記のような行列を、'''正規モード行列'''という{{Sfn|横山・日野・芳村|2015|p=111}}。 {{NumBlk|:|<math> \boldsymbol{\overline{U}} = (\boldsymbol{\overline{u}}_{1},\ \boldsymbol{\overline{u}}_{2},\ \cdots,\ \boldsymbol{\overline{u}}_{n}) = \begin{pmatrix} \overline{u}_{1,1} & \overline{u}_{1,2} & \cdots & \overline{u}_{1,n} \\ \overline{u}_{1,2} & \overline{u}_{2,2} & \cdots & \overline{u}_{2,n} \\ \vdots & \vdots & \ddots & \vdots \\ \overline{u}_{1,n} & \overline{u}_{n,2} & \cdots & \overline{u}_{n,n} \end{pmatrix} </math>|{{EquationRef|3.19}}}} 質量正規固有モードを採用すると、式{{EqNoteN|3.17}}に対して剛性行列は固有角振動数と {{NumBlk|:|<math> \boldsymbol{\overline{u}}_r^\top \boldsymbol{K} \boldsymbol{\overline{u}}_r = \omega_r^2 </math>|{{EquationRef|3.20}}}} という関係を持つ{{Sfn|末岡・金光・近藤|2000|p=50}}。 ===モード座標と一般解=== 固有モード（固有ベクトル）は互いに直交なので、これらを使った[[線形結合]]で変位ベクトル '''''x''''' を表すことができる{{Sfnm|長松・長松|2018|1p=85|横山・日野・芳村|2015|2p=112}}。つまり、固有モード '''''U''''' を用いて、'''''x''''' を {{NumBlk|:|<math> \boldsymbol{x} = q_1 \boldsymbol{\overline{u}}_{1} + q_2 \boldsymbol{\overline{u}}_{2}+ \dots + q_n \boldsymbol{\overline{u}}_{n} = \boldsymbol{\overline{U} q} </math>|{{EquationRef|3.21}}}} と表すことができる{{Sfn|藤田|2016|p=130}}。ここで、''q'' は時間の未知関数で、'''モード座標'''、'''基準座標'''、'''正規座標'''、'''規準座標'''、'''主座標'''などと呼ばれる{{Sfnm|藤田|2016|1p=131|入江・小林|2006|2p=106|横山・日野・芳村|2015|3p=112}}。特に正規化した固有モードを使っているモード座標を指して、正規座標と呼ぶこともある{{Sfn|砂子田・伊藤・鄭・平元|2012|p=92}}。'''''q''''' は ''q<sub>r</sub>'' を 1 次から ''n'' 次まで並べた次のような縦ベクトルである。 :<math> \boldsymbol{q} = \begin{pmatrix} q_1 \\ q_2 \\ \vdots \\ q_n \end{pmatrix} </math> モード座標 '''''q''''' に対して、元の座標 '''''x''''' を'''物理座標'''と呼ぶ{{Sfn|末岡・金光・近藤|2000|p=53}}。{{EqNoteN|3.21}}は、正規固有モードを介して物理座標をモード座標に変換していることを意味する{{Sfn|末岡・金光・近藤|2000|p=53}}。モード座標 '''''q'''<sub>r</sub>'' は、'''''x''''' に対して ''r'' 次固有モード'''''ū'''<sub>r</sub>'' が寄与する度合いを表しているとも言える{{Sfn|長松・長松|2018|p=85}}。式{{EqNoteN|3.21}}は展開定理とも呼ばれる{{Sfn|宮本・永野・藤谷・吉村|2014|p=89}}。 式{{EqNoteN|3.21}}を運動方程式{{EqNoteN|3.1}}へ代入して、左から転置したモード行列 '''''U'''''<sup>⊤</sup> を掛けると次のようになる{{Sfn|藤田|2016|p=131}}。 {{NumBlk|:|<math> \boldsymbol{U^\top M U \ddot{q}} + \boldsymbol{U^\top K U q} = \mathbf{0} </math>|{{EquationRef|3.22}}}} 固有モードの直交性（式{{EqNoteN|3.14}}と式{{EqNoteN|3.15}}）を上式に当てはめると、 {{NumBlk|:|<math> \begin{cases} M_1 \ddot{q}_1 + K_1 q_1 = 0\\ M_2 \ddot{q}_2 + K_2 q_2 = 0\\ \vdots \\ M_n \ddot{q}_n + K_n q_n = 0\\ \end{cases} </math>|{{EquationRef|3.23}}}} というような ''q'' に関する ''n'' 個の運動方程式が得られる{{Sfn|藤田|2016|p=131}}。それぞれの式の両辺を ''M<sub>r</sub>'' で割ると、下記のような形になる{{Sfn|藤田|2016|p=131}}。 {{NumBlk|:|<math> \begin{cases} \ddot{q}_1 + \omega_1^2 q_1 = 0\\ \ddot{q}_2 + \omega_2^2 q_2 = 0\\ \vdots \\ \ddot{q}_n + \omega_n^2 q_n = 0\\ \end{cases} </math>|{{EquationRef|3.24}}}} [[File:2自由度系振動時系列グラフ.svg|thumb|280px|式{{EqNoteN|2.7}}で表現される自動車の上下・ピッチング振動モデルの振動の例。図はピッチング ''θ'' の振動を示している。パラメータは ''m'' = 1600 kg, ''I<sub>G</sub>'' = 2500 kg-m<sup>2</sup>, ''k''<sub>1</sub> = 3500 N/m, ''k''<sub>2</sub> = 4100 N/m, ''l''<sub>1</sub> = 1.4 m, ''l''<sub>2</sub> = 1.6 m で{{Sfn|入江・小林|2006|p=89}}、初期条件は上下変位速度が ''ẋ'' = 1 m/s で他は全て 0 という条件の例。]] 式{{EqNoteN|3.24}}は非連成化されており、各式はそれぞれ独立している{{Sfn|横山・日野・芳村|2015|p=112}}。そのため、一つ一つの式は1自由度系の不減衰自由振動と同じであるから、''q'' の各解は以下のようになる{{Sfn|横山・日野・芳村|2015|p=112}}。 {{NumBlk|:|<math> \begin{cases} q_1 = c_1 \sin (\omega_1 t+ \theta_1) \\ q_2 = c_n \sin (\omega_2 t+ \theta_2) \\ \vdots \\ q_n = c_n \sin (\omega_n t+ \theta_n) \\ \end{cases} </math>|{{EquationRef|3.25}}}} ここで、''c<sub>r</sub>'' と ''θ<sub>r</sub>'' (''r'' = 1, 2, …, ''n'') は初期条件で決まる定数で、''c<sub>r</sub>'' が[[振幅]]、''θ<sub>r</sub>'' が初期[[位相]]を意味する{{Sfn|横山・日野・芳村|2015|p=112}}。'''''x''''' の一般解は、式{{EqNoteN|3.25}}を式{{EqNoteN|3.21}}に代入して {{NumBlk|:|<math> \boldsymbol{x} = \sum_{r=1}^n \boldsymbol{\overline{u}}_{r} C_r \sin (\omega_r t+ \theta_r) </math>|{{EquationRef|3.26}}}} と得られる{{Sfn|横山・日野・芳村|2015|p=112}}。すなわち、''n'' 自由度系における各自由度の自由振動は、''n'' 個の[[調和振動]]の重ね合わせ（和）となっており、それら調和振動のそれぞれの振動数は 1 次から ''n'' 次までの固有角振動数となっている{{Sfnm|横山・日野・芳村|2015|1p=113|安田|2012|2p=122}}。固有モードは、対応する固有角振動数を持つ調和振動成分の各自由度間の振幅比を定めている{{Sfn|安田|2012|p=122}}。位相角を陽に表さずに、余弦関数と正弦関数の和 {{NumBlk|:|<math> \boldsymbol{x} = \sum_{r=1}^n \boldsymbol{\overline{u}}_{r} (a_r \cos \omega_r t + b_r \sin \omega_r t) </math>|{{EquationRef|3.27}}}} や、指数の正負が異なる複素指数関数の和 {{NumBlk|:|<math> \boldsymbol{x} = \sum_{r=1}^n \boldsymbol{\overline{u}}_{r} (d_r e^{i \omega_r t} + g_r e^{- i \omega_r t}) </math>|{{EquationRef|3.28}}}} といった形で一般解を表すこともできる{{Sfn|砂子田・伊藤・鄭・平元|2012|pp=88–89}}。''a<sub>r</sub>'', ''b<sub>r</sub>'', ''d<sub>r</sub>'', ''g<sub>r</sub>'' も初期条件で決まる定数で、式{{EqNoteN|3.27}}と式{{EqNoteN|3.27}}は式{{EqNoteN|3.26}}へ式変形可能な同値な式である{{Sfn|砂子田・伊藤・鄭・平元|2012|pp=88–89}}。 以上のような、物理座標をモード座標へ変換し、固有モードの直交性を利用して多自由度系の問題を1自由度系の重ね合わせの問題に帰着させ、解析を行う手法を[[モード解析]]（特にモード重畳法）という{{Sfn|末岡・金光・近藤|2000|p=54}}<ref name="小松2009"/>。後述のように、モード解析は特に[[有限要素法]]へ適用することで有効性を発揮する{{Sfn|藤田|2016|p=199}}。実物の振動特性を求める上でもモード解析の理論を適用することで各特性を[[システム同定|同定]]でき、この手法は[[実験モード解析]]と呼ばれる{{Sfn|長松・長松|2018|pp=12–13}}。 ==減衰自由振動== ===一般の粘性減衰=== [[File:Three point masses vibration system.svg|thumb|x260px|基礎の上でばねと減衰器が一組となって質点と連結し、直列に連なった3自由度減衰系の例]] 減衰行列 '''''C''''' が存在する線形多自由度系の振動について考える。外力が無い場合の ''n'' 自由度系の運動方程式は以下のようになる{{Sfn|長松・長松|2018|p=91}}。 {{NumBlk|:|<math> \boldsymbol{M \ddot{x}} + \boldsymbol{C \dot{x}} + \boldsymbol{K x} = \boldsymbol{0} </math>|{{EquationRef|4.1}}}} 基礎の上でばねと減衰器が一組となって質点と連結し、直列に連なった3自由度減衰系の典型的な例では、'''''M''''', '''''C''''', '''''K''''' は次のような行列となる{{Sfn|背戸・丸山|2002|p=222}}。 {{NumBlk|:|<math> \boldsymbol{M} = \begin{pmatrix} m_1 & 0 & 0 \\ 0 & m_2 & 0 \\ 0 & 0 & m_3 \end{pmatrix} </math>|{{EquationRef|4.2}}}} {{NumBlk|:|<math> \boldsymbol{C} = \begin{pmatrix} c_1 & -c_1 & 0 \\ -c_1 & c_1 + c_2 & -c_2 \\ 0 & -c_2 & c_2 + c_3 \end{pmatrix} </math>|{{EquationRef|4.3}}}} {{NumBlk|:|<math> \boldsymbol{K} = \begin{pmatrix} k_1 & -k_1 & 0 \\ -k_1 & k_1 + k_2 & -k_2 \\ 0 & -k_2 & k_2 + k_3 \end{pmatrix} </math>|{{EquationRef|4.4}}}} 式{{EqNoteN|4.1}}の解を式{{EqNoteN|3.2}}と同じように仮定し、式{{EqNoteN|3.2}}を式{{EqNoteN|4.1}}に代入して整理すると、 {{NumBlk|:|<math> (\lambda^2 \boldsymbol{M} + \lambda \boldsymbol{C} + \boldsymbol{K})\boldsymbol{u} = \boldsymbol{0} </math>|{{EquationRef|4.5}}}} となり、'''''u''''' = '''0''' 以外の解を持つという条件から {{NumBlk|:|<math> \mbox{det}(\lambda^2 \boldsymbol{M} + \lambda \boldsymbol{C} + \boldsymbol{K}) </math>|{{EquationRef|4.6}}}} という特性方程式が得られる<ref name="安田2001">{{Cite book ja-jp |author = 安田 仁彦 |title = 振動工学 ―応用編 |publisher = コロナ社 |url = https://www.coronasha.co.jp/np/isbn/9784339045581/ |edition = 初版 |year = 2001 |isbn = 978-4-339-04558-1 }} p. 77–78</ref>。式{{EqNoteN|4.6}}を展開すると、''λ'' に関する 2''n'' 次多項式となる<ref name="安田2001"/>。この多項式を解いて得た解 ''λ''<sub>1</sub>, ''λ''<sub>2</sub>, …, ''λ''<sub>2''n''</sub> を式{{EqNoteN|4.5}}に代入し、'''''u''''' を定め、'''''u''''' に適当な正規化を行う<ref name="安田2001"/>。正規化された各ベクトルを '''''ū'''''<sub>1</sub>, '''''ū'''''<sub>2</sub>, …, '''''ū'''''<sub>2''n''</sub> と表すと、式{{EqNoteN|4.1}}の一般解は {{NumBlk|:|<math> \boldsymbol{x} = \boldsymbol{\overline{u}}_1 a_1 e^{\lambda_1 t} + \boldsymbol{\overline{u}}_2 a_2 e^{\lambda_2 t} + \cdots + \boldsymbol{\overline{u}}_{2n} a_{2n} e^{\lambda_{2n} t} </math>|{{EquationRef|4.7}}}} となる<ref name="安田2001"/>。ここで、''a''<sub>1</sub>, ''a''<sub>2</sub>, …, ''a''<sub>2''n''</sub> は初期条件で決まる定数である<ref name="安田2001"/>。2''n'' 個の ''λ'' は、解が[[減衰振動]]であれば ''n'' 組の[[複素共役|互いに共役な複素数]]になる<ref name="安田2001"/>{{Sfn|長松・長松|2018|p=96}}。このとき、固有モードも複素数となり、'''複素固有モード'''と呼ばれる{{Sfn|長松・長松|2018|p=96}}。 ===比例粘性減衰=== 式{{EqNoteN|4.1}}のように一般的な減衰行列 '''''C''''' が運動方程式に存在する場合、正規モード行列によって '''''M''''' と '''''K''''' は対角化できるが、'''''C''''' も同時に対角化することはできない{{Sfn|末岡・金光・近藤|2000|p=61}}。そのため、不減衰自由振動で可能だったモード座標に変換しての非連成化が不可能となり、モード解析の利点を活かすことができなくなる{{Sfn|末岡・金光・近藤|2000|p=61}}。そこで、'''''C''''' が下記のような'''比例粘性減衰'''として与えられると仮定し、実際の振動解析が行われることも多い{{Sfn|末岡・金光・近藤|2000|p=61}}。 {{NumBlk|:|<math> \boldsymbol{C} = \alpha \boldsymbol{M} + \beta \boldsymbol{K} </math>|{{EquationRef|4.8}}}} ここで、''α'' と ''β'' は定数で、'''比例減衰定数'''と呼ばれる<ref>{{Cite book ja-jp |editor = 日本機械学会 | title = 振動のダンピング技術 | publisher = 養賢堂 | year = 1998 | edition = 第1版 | isbn = 4-8425-9816-6 }} p.60</ref>。比例粘性減衰を有する系を'''比例粘性減衰系'''などと呼ぶ{{Sfn|長松・長松|2018|p=92}}。比例粘性減衰が成り立つと仮定すれば減衰行列を対角化できる{{Sfn|入江・小林|2006|p=107}}。 式{{EqNoteN|4.8}}の形で与えられる比例粘性減衰は、特に'''レイリー減衰'''や'''レイリー型減衰'''と呼ばれる{{Sfnm|藤田|2016|1p=134|宮本・永野・藤谷・吉村|2014|2p=92}}。''α'''M''''' のみで仮定されるものは'''質量比例型減衰'''、''β'''K''''' のみで仮定されるものは'''剛性比例型減衰'''などと呼ぶ{{Sfn|宮本・永野・藤谷・吉村|2014|p=92}}。 実際の減衰が比例粘性減衰になっていることはまれであり、比例粘性減衰はあくまでも近似的なものである{{Sfnm|末岡・金光・近藤|2000|1p=61|平井・水田|2018|2p=118|入江・小林|2006|3p=107}}。しかし、比例粘性減衰の仮定を導入することで、不減衰系と同じ取り扱いが可能となり、モード解析の手法が適用可能になる{{Sfnm|藤田|2016|1p=134|末岡・金光・近藤|2000|2p=61}}。もし減衰が全体に分布しているような構造であれば、適当な比例減衰定数を設定すれば、実際の現象を実用問題ないレベルで再現できるという一定の妥当性もある{{Sfn|長松・長松|2018|p=95}}。減衰は摩擦・材料減衰・流体粘性など様々な要因で起こるため、そもそも減衰の適切な定式化自体が難しいといった事情もある{{Sfn|長松・長松|2018|p=94}}。式{{EqNoteN|2.3}}のように速度の比例定数として与えられる一般の粘性減衰も、多種多様な発生機構によって減衰が起きるという実情に起因して厳密には成立しない{{Sfn|長松・長松|2018|p=96}}。 式{{EqNoteN|4.8}}を式{{EqNoteN|4.1}}に代入した場合、 {{NumBlk|:|<math> \boldsymbol{M \ddot{x}} + (\alpha \boldsymbol{M} + \beta \boldsymbol{K}) \boldsymbol{\dot{x}} + \boldsymbol{K x} = \boldsymbol{0} </math>|{{EquationRef|4.9}}}} という運動方程式になる{{Sfnm|長松・長松|2018|1p=92|入江・小林|2006|2p=107}}。解を式{{EqNoteN|3.2}}のように仮定して上式に代入し、整理すると、 {{NumBlk|:|<math> \{ (\lambda^2 + \alpha \lambda) \boldsymbol{M} + (\beta \lambda + 1) \boldsymbol{K} \} \boldsymbol{u} = \boldsymbol{0} </math>|{{EquationRef|4.10}}}} となる{{Sfnm|長松・長松|2018|1p=92|入江・小林|2006|2p=107}}。さらに、 {{NumBlk|:|<math> \gamma^2 = \frac{\lambda^2 + \alpha \lambda}{\beta \lambda + 1} </math>|{{EquationRef|4.11}}}} とおけば、式{{EqNoteN|4.10}}は {{NumBlk|:|<math> ( \gamma^2 \boldsymbol{M} + \boldsymbol{K} ) \boldsymbol{u} = \boldsymbol{0} </math>|{{EquationRef|4.12}}}} となる{{Sfnm|長松・長松|2018|1pp=92–93|入江・小林|2006|2pp=107–108}}。式{{EqNoteN|4.12}}は、（複素指数関数や三角関数で解を仮定したときの）不減衰振動の固有値問題における ''λ'' を ''γ'' に置き換えただけの式になる{{Sfnm|長松・長松|2018|1pp=92–93|入江・小林|2006|2pp=107–108}}。したがって、比例粘性減衰を仮定した減衰振動の固有モードは、同一の質量行列と剛性行列を有する不減衰振動の固有モードと同じである{{Sfnm|長松・長松|2018|1p=93|入江・小林|2006|2p=108}}。 一方、比例粘性減衰を仮定した減衰振動の固有角振動数は、同一の質量行列と剛性行列の不減衰振動の固有角振動数よりも小さくなる{{Sfn|長松・長松|2018|p=94}}。モード行列 '''''U''''' で減衰行列を対角化すると {{NumBlk|:|<math> \boldsymbol{U}^\top \boldsymbol{C} \boldsymbol{U} = \begin{pmatrix} \alpha M_1 + \beta K_1 & 0 & \cdots & 0 \\ 0 & \alpha M_2 + \beta K_2 & \cdots & 0 \\ \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & \cdots & \alpha M_n + \beta K_n \end{pmatrix} </math>|{{EquationRef|4.13}}}} となる{{Sfn|入江・小林|2006|p=107}}。この行列の成分 ''αM<sub>r</sub> + βK<sub>r</sub>'' (''r'' = 1, 2, …, ''n'') を'''モード減衰'''や'''モード減衰係数'''と呼び、 ''C<sub>r</sub>'' などで表す{{Sfnm|長松・長松|2018|1p=93|入江・小林|2006|2p=107}}。ここで、''M<sub>r</sub>'' はモード質量、''K<sub>r</sub>'' はモード剛性である。さらに、 {{NumBlk|:|<math> C_{c,r} = 2 \sqrt{ M_r K_r} </math>|{{EquationRef|4.14}}}} という量を導入して、これでモード減衰を割った量 {{NumBlk|:|<math> \zeta_r = \frac{ C_r }{ C_{c,r} } = \frac{ \alpha M_r + \beta K_r }{ 2 \sqrt{ M_r K_r} } </math>|{{EquationRef|4.15}}}} を'''モード減衰比'''と呼ぶ{{Sfnm|長松・長松|2018|1p=94|入江・小林|2006|2p=108}}。多自由度系の減衰系では固有モードごとに減衰の効果が異なっており、モード減衰比が固有モードごとの減衰の効果の程度を表している<ref name="jsme_モード減衰比">{{Cite web|和書|author= 日本機械学会 | work = 機械工学事典 |url= https://www.jsme.or.jp/jsme-medwiki/13:1012759 |title= モード減衰比 |accessdate= 2021-02-13}}</ref>。''ζ<sub>r</sub>'' ≥ 1 ならば過減衰の状態であり、その固有モードの振動は起こらない<ref name="jsme_モード減衰比"/>{{Sfn|長松・長松|2018|p=94}}。''ζ<sub>r</sub>'' < 1 ならば減衰振動となり、その固有角振動数は {{NumBlk|:|<math> \omega_{d,r} = \omega_r \sqrt{ 1 - \zeta_r^2 } </math>|{{EquationRef|4.16}}}} で与えられる{{Sfn|長松・長松|2018|p=94}}。''ω<sub>d,r</sub>'' を減衰固有角振動数と呼ぶ{{Sfn|末岡・金光・近藤|2000|p=62}}。以上のように、線形1自由度系の減衰振動の考え方が固有モードごとの振動にも当てはまる{{Sfn|長松・長松|2018|p=94}}。正規座標へ変換を行い、モード減衰比と固有角振動数を用いて運動方程式を表すと、下記のように表現できる{{Sfn|末岡・金光・近藤|2000|p=62}}。 {{NumBlk|:|<math> \begin{cases} \ddot{q}_1 + 2 \zeta_{1} \omega_{1} \dot{q}_1 + \omega_{1}^2 q_1 = 0 \\ \ddot{q}_2 + 2 \zeta_{2} \omega_{2} \dot{q}_2 + \omega_{2}^2 q_2 = 0 \\ \vdots \\ \ddot{q}_n + 2 \zeta_{r} \omega_{r} \dot{q}_r + \omega_{n}^2 q_n = 0 \\ \end{cases} </math>|{{EquationRef|4.17}}}} ''ζ<sub>r</sub>'' < 1 であれば、''q<sub>r</sub>'' の一般解は積分定数を ''a<sub>r</sub>'' と ''b<sub>r</sub>'' として、 {{NumBlk|:|<math> q_r = e^{-\zeta_{r} \omega_{r} t}(a_r \cos \omega_{d,r} t + b_r \sin \omega_{d,r} t) </math>|{{EquationRef|4.18}}}} となる{{Sfn|末岡・金光・近藤|2000|p=62}}。物理座標 '''''x''''' へ逆変換すれば、比例粘性減衰系の自由振動の一般解は下記のようになる{{Sfn|末岡・金光・近藤|2000|p=62}}。 {{NumBlk|:|<math> \boldsymbol{x} = \sum_{r=1}^n \boldsymbol{\overline{u}}_{r} e^{-\zeta_{r} \omega_{r} t}(a_r \cos \omega_{d,r} t + b_r \sin \omega_{d,r} t) </math>|{{EquationRef|4.19}}}} ==強制振動== ===不減衰系=== 任意の励振力を受ける減衰の無い系の運動方程式は以下のように表される{{Sfn|横山・日野・芳村|2015|p=115}}。 {{NumBlk|:|<math> \boldsymbol{M \ddot{x}} + \boldsymbol{K x} = \boldsymbol{f} </math>|{{EquationRef|5.1}}}} 例えば、励振力が角振動数 ''Ω'' の余弦関数で与えられるとすれば、 {{NumBlk|:|<math> \boldsymbol{M \ddot{x}} + \boldsymbol{K x} = \boldsymbol{f_a} \cos \mathit{\Omega} t </math>|{{EquationRef|5.2}}}} となる{{Sfn|末岡・金光・近藤|2000|pp=54–55}}。ここで、'''''f<sub>a</sub>''''' は、以下のような各自由度に対する励振力の振幅値の縦ベクトルである。 :<math> \boldsymbol{f_a} = \begin{pmatrix} f_{a,1} \\ f_{a,2} \\ \vdots \\ f_{a,n} \end{pmatrix} </math> 式{{EqNoteN|5.2}}に対して特解を '''''x''''' = '''''u''''' cos''Ωt'' と仮定し、式{{EqNoteN|5.2}}に代入すれば {{NumBlk|:|<math> (\boldsymbol{M} - \mathit{\Omega}^2 \boldsymbol{K})\boldsymbol{u} = \boldsymbol{f_a} </math>|{{EquationRef|5.3}}}} となるので、'''''u''''' および '''''x''''' は、 {{NumBlk|:|<math> \boldsymbol{u} = (\boldsymbol{M} - \mathit{\Omega}^2 \boldsymbol{K})^{-1} \boldsymbol{f_a} </math>|{{EquationRef|5.4}}}} {{NumBlk|:|<math> \boldsymbol{x} = (\boldsymbol{M} - \mathit{\Omega}^2 \boldsymbol{K})^{-1} \boldsymbol{f_a} \cos \mathit{\Omega} t </math>|{{EquationRef|5.5}}}} となる{{Sfnm|末岡・金光・近藤|2000|1p=55|安田|2012|2p=131}}。ここで "<sup>−1</sup>" は[[逆行列]]を意味する。したがって、逆行列 ('''''M''''' − ''Ω''<sup>2</sup>'''''K''''')<sup>−1</sup> を計算すれば式{{EqNoteN|5.5}}から '''''x''''' の値が分かる。しかし、この逆行列を解析的に解くことは困難で、数値計算を行うにしても自由度の数が増えると膨大な計算量になる{{Sfn|末岡・金光・近藤|2000|p=55}}。もし励振力の角振動数 ''Ω'' を変えると、そのたびに逆行列を計算する必要がある{{Sfn|安田|2012|p=132}}{{Sfn|末岡・金光・近藤|2000|p=55}}。そのため、実際に '''''x''''' の解を求めるために行われるのは、下記のようなモード解析による手法である{{Sfnm|安田|2012|1p=132|末岡・金光・近藤|2000|2p=55}}。 解を得るために、'''''x''''' を複素数に拡張し、励振力 '''''f''''' を[[複素指数関数]] '''''f<sub>a</sub>'''e<sup> jΩt</sup>'' の形で与えるとする{{Sfnm|安田|2012|1pp=130–131|藤田|2016|2p=132}}。この場合、計算して解が得られた後に実部あるいは虚部を取ることで、実際の解が得られる{{Sfnm|安田|2012|1pp=130–131|藤田|2016|2p=132}}。運動方程式{{EqNoteN|5.1}}の右辺を '''0''' としたときのモード行列 '''''U''''' が、事前に求められているとする。モード座標への変換式{{EqNoteN|3.21}}を運動方程式{{EqNoteN|5.1}}へ適用して、左から '''''U'''''<sup>⊤</sup> を掛け、式{{EqNoteN|3.14}}と式{{EqNoteN|3.15}}の対角化を適用する。すると、励振力を受ける不減衰系の運動方程式は {{NumBlk|:|<math> \begin{cases} M_1 \ddot{q}_1 + K_1 q_1 = F_{1} e^{j \omega t} \\ M_2 \ddot{q}_2 + K_2 q_2 = F_{2} e^{j \omega t} \\ \vdots \\ M_n \ddot{q}_n + K_n q_n = F_{n} e^{j \omega t} \\ \end{cases} </math>|{{EquationRef|5.6}}}} という独立・非連成の ''n'' 個の運動方程式に帰着する{{Sfnm|安田|2012|1pp=132–133|藤田|2016|2pp=132–133}}。ただし、右辺の ''F<sub>r</sub>'' は次のような値である{{Sfnm|安田|2012|1p=133|藤田|2016|2p=133}}。 {{NumBlk|:|<math> F_r = \boldsymbol{u}_r^\top \boldsymbol{f_a} </math>|{{EquationRef|5.7}}}} 式{{EqNoteN|5.6}}は線形1自由度系と同じなので、[[強制振動]]を表す特解は、 {{NumBlk|:|<math> q_r = \frac{ F_r }{ K_r - M_r \mathit{\Omega}^2 } e^{j \mathit{\Omega} t} </math>|{{EquationRef|5.8}}}} となる{{Sfnm|安田|2012|1p=134|藤田|2016|2p=133}}。式{{EqNoteN|3.21}}を使ってモード座標を物理座標へ逆変換し、固有角振動数を使って整理すれば強制振動の解は次のようになる{{Sfnm|安田|2012|1p=134|藤田|2016|2p=133}}。 {{NumBlk|:|<math> \boldsymbol{x} = \sum_{r=1}^n \frac{ F_r }{ K_r - M_r \mathit{\Omega}^2 } \boldsymbol{u}_r e^{j \mathit{\Omega} t} = \sum_{r=1}^n \frac{ F_r }{ M_r (\omega_r^2 - \mathit{\Omega}^2) } \boldsymbol{u}_r e^{j \mathit{\Omega} t} </math>|{{EquationRef|5.9}}}} 式{{EqNoteN|5.9}}から分かることは、励振力の角振動数 ''Ω'' が系の固有角振動数 ''ω''<sub>1</sub>, ''ω''<sub>2</sub>, …, ''ω<sub>n</sub>'' のどれかに近いとき、係数の[[極限]]が 0 となってその振動成分が極めて大きくなり、[[共振]]が起こるという点である{{Sfn|安田|2012|p=134}}。''Ω'' が固有角振動数のいずれかに一致するとき、'''''x''''' の振幅は無限大に発散する{{Sfn|末岡・金光・近藤|2000|p=56}}。 ===比例粘性減衰系=== 系が比例粘性減衰系であれば、励振力を受ける場合でもモード座標変換によって独立した1自由度系に分解でき、大自由度の問題も扱いが容易になる{{Sfn|長松・長松|2018|p=97}}。一般の運動方程式{{EqNoteN|2.1}}に、レイリー減衰の式{{EqNoteN|4.8}}およびモード座標への変換式{{EqNoteN|3.21}}を代入する。左から転置した質量正規固有モードの正規モード行列 '''''Ū'''''<sup>⊤</sup> を掛けて対角化する。さらに、2''ζ<sub>r</sub>ω<sub>r</sub>'' = ''α'' + ''βω<sub>r</sub>''<sup>2</sup> という関係を用いれば、運動方程式は下記のような独立・非連成の ''n'' 個の方程式に帰着する{{Sfn|下郷・田島|2002|pp=109–110}}。 {{NumBlk|:|<math> \begin{cases} \ddot{q}_1 + 2 \zeta_1 \omega_1 \dot{q}_1 + \omega_1^2 q_1 = \boldsymbol{\overline{u}}_1^\top \boldsymbol{f} \\ \ddot{q}_2 + 2 \zeta_2 \omega_2 \dot{q}_2 + \omega_2^2 q_2 = \boldsymbol{\overline{u}}_2^\top \boldsymbol{f} \\ \vdots \\ \ddot{q}_n + 2 \zeta_n \omega_n \dot{q}_n + \omega_n^2 q_n = \boldsymbol{\overline{u}}_n^\top \boldsymbol{f} \\ \end{cases} </math>|{{EquationRef|5.10}}}} 励振力 '''''f''''' が '''''f<sub>a </sub>'''''cos''Ωt'' で与えられるとする。線形1自由度系と同様の解法によって、正規座標 ''q<sub>r</sub>'' の強制振動の解が下記のように求まる{{Sfn|日本機械学会（編）|2004|p=42}}。 {{NumBlk|:|<math> q_r = \dfrac{ \boldsymbol{\overline{u}}_{r}^\top \boldsymbol{f_a} }{ \omega_r^2 \sqrt{ \left \{ 1 - \left ( \dfrac{ \mathit{\Omega} }{ \omega_r } \right ) ^2 \right \}^2 + \left \{ 2 \zeta_r \dfrac{ \mathit{\Omega} }{ \omega_r } \right \} ^2 } } \cos (\mathit{\Omega} t - \phi_r ) </math>|{{EquationRef|5.11}}}} {{NumBlk|:|<math> \phi_r = \tan^{-1} \left ( \frac{2 \zeta_r \omega_r \mathit{\Omega} }{\omega_r^2 - \mathit{\Omega}^2} \right ) </math>|{{EquationRef|5.12}}}} 解を複素数で表現した場合は、 {{NumBlk|:|<math> q_r = \dfrac{ \boldsymbol{\overline{u}}_{r}^\top \boldsymbol{f_a} }{ \omega_r^2 \left \{ 1 - \left ( \dfrac{ \omega_r }{ \mathit{\Omega} } \right ) ^2 + j 2 \zeta_r \dfrac{ \omega_r }{ \mathit{\Omega} } \right \} } e^{j \mathit{\Omega} t} </math>|{{EquationRef|5.13}}}} となり、物理座標の解は {{NumBlk|:|<math> \boldsymbol{x} = \sum_{r=1}^n \frac{ \boldsymbol{ \overline{u}_r^\top f_a \overline{u}_r }}{ \omega_r^2 \left \{ 1 - \left ( \dfrac{ \mathit{\Omega} }{ \omega_r } \right ) ^2 + j 2 \zeta_r \dfrac{ \mathit{\Omega} }{ \omega_r } \right \} } e^{j \mathit{\Omega} t} </math>|{{EquationRef|5.14}}}} となる{{Sfn|日本機械学会（編）|2004|p=42}}。''Ω''/''ω<sub>r</sub>'' は強制振動比などと呼ばれる{{Sfn|背戸・丸山|2002|p=220}}。通常の固有モードの場合は次式となる{{Sfn|入江・小林|2006|p=108}}。 {{NumBlk|:|<math> \boldsymbol{x} = \sum_{r=1}^n \frac{ \boldsymbol{ u_r^\top f_a u_r }}{ K_r \left \{ 1 - \left ( \dfrac{ \mathit{\Omega} }{ \omega_r } \right ) ^2 + j 2 \zeta_r \dfrac{ \mathit{\Omega} }{ \omega_r } \right \} } e^{j \mathit{\Omega} t} </math>|{{EquationRef|5.15}}}} ===周波数応答関数=== 強制振動の特解から、[[周波数応答関数]]あるいは[[伝達関数]]が求められる{{Sfnm|末岡・金光・近藤|2000|1p=57|入江・小林|2006|2p=109}}。周波数応答関数は、振動系への入力が調和振動の形で与えられるときの、出力と入力の振幅比および位相差を周波数の関数として表す{{Sfn|下郷・田島|2002|pp=56–57}}。比例減衰系において、角振動数 ''Ω'' の調和外力が ''q'' 番目の自由度のみに加わるとする。このときの ''p'' 番目の自由度の振動応答は {{NumBlk|:|<math> G_{pq} ( \Omega ) = \sum_{r=1}^n \frac{ u_{pr} u_{qr} }{ K_r \left \{ 1 - \left ( \dfrac{ \mathit{\Omega} }{ \omega_r } \right ) ^2 + j 2 \zeta_r \dfrac{ \mathit{\Omega} }{ \omega_r } \right \} } </math>|{{EquationRef|5.16}}}} と表される{{Sfn|長松・長松|2018|p=98}}。この場合の周波数応答関数 ''G''(''Ω'') は[変位]/[力]の単位を持ち、[[機械的コンプライアンス|コンプライアンス]]とも呼ばれる{{Sfnm|長松・長松|2018|1p=98|背戸・丸山|2002|2p=226}}。''q'' 番目加振・''p'' 番目応答の周波数応答関数と ''p'' 番目加振・''q'' 番目応答の周波数応答関数は、[[マクスウェル・ベティの相反作用の定理|相反定理]]によって互いに等しく、''G<sub>pq</sub>''(''Ω'') = ''G<sub>qp</sub>''(''Ω'') の関係である{{Sfn|長松・長松|2018|pp=98–99}}。不減衰系の場合は、周波数応答関数は次式となる{{Sfn|日本機械学会（編）|2004|p=214}}。 {{NumBlk|:|<math> G_{pq}(\omega) = \sum_{r=1}^n \frac{ u_{pr} u_{qr} }{ K_r - M_r \mathit{\Omega}^2 } </math>|{{EquationRef|5.17}}}} [[File:Frequency response of 3 degrees of freedom system.png|thumb|300px|式{{EqNoteN|4.2}},{{EqNoteN|4.3}},{{EqNoteN|4.4}}によって与えられた3自由度減衰系の周波数応答関数の応答曲線の例。図は質点 ''m''<sub>3</sub> が加振されたときの質量 ''m''<sub>2</sub> の応答の例。山になっている点が共振点で、谷になっている点が反共振点。]] これらの周波数応答関数によって、系の周波数特性が把握できる{{Sfn|末岡・金光・近藤|2000|p=57}}。線形多自由度系の場合は、固有角振動数と固有モードを求めれば、調和外力による強制振動に対する周波数特性も同時に理解することができる{{Sfn|末岡・金光・近藤|2000|p=57}}。横軸に振動数や強制振動比を取り、縦軸に調和外力に対する振幅比や位相差を図示したものを応答曲線や共振曲線と呼ぶ{{Sfnm|末岡・金光・近藤|2000|1pp=24–25|吉川・松井・石井|1987|2p=70}}。縦軸を[[デシベル]]にして図示したものは特に[[ボード線図]]と呼ばれる{{Sfn|吉川・松井・石井|1987|p=72}}。系の振動特性を振動実験から同定する実験モード解析では、実験値に上記の周波数応答関数で[[カーブフィッティング]]して系の各種パラメータを推定する{{Sfn|日本機械学会（編）|2004|p=216}}。 上述のように、励振力の振動数が系の固有角振動数に一致すると、振動は極大化する。一方、多自由度系では共振して振幅が発散する現象だけでなく、励振力が特定の角振動数のときに周波数応答関数が極小になる現象もある{{Sfn|長松・長松|2018|p=97}}。このような現象を[[反共振]]と呼ぶ{{Sfn|長松・長松|2018|p=97}}。ただし、共振が系全体で起こるのに対して、反共振は一組の加振点・応答点ごとにしか起こらない{{Sfn|長松・長松|2018|p=100}}。周波数応答関数の絶対値を縦軸にした応答曲線上では、共振点は曲線の鋭い山のように現れ、反共振は曲線の鋭い谷のように現れる<ref>{{Cite web|和書|author= 日本機械学会 | work = 機械工学事典 |url= https://www.jsme.or.jp/jsme-medwiki/13:1010422 |title= 反共振 |accessdate= 2021-03-14 }}</ref>。反共振は、代表的な[[制振|制振器]]である[[動吸振器]]の原理として活用される{{Sfnm|藤田|2016|1p=104|背戸・丸山|2002|2p=178}}。振動を抑えたい対象にばね・ダンパを介して付加質量を取り付けることによって振動を抑制でき、地震や強風に対する建築構造物の防振や回転機械の防振などに使われる{{Sfnm|下郷・田島|2002|1p=101|背戸・丸山|2002|2pp=178–179}}。 ===刺激係数=== '''刺激係数'''は、[[地震]]のように構造物の基礎が揺れている場合に、その励振が各モードに対してどのぐらい強く寄与するのかを表す{{Sfn|宮本・永野・藤谷・吉村|2014|p=94}}。構造物を多質点系でモデル化して、基礎の変位を ''y'' として、各質点の基礎からの相対変位を ''x'' とする。各質点は ''ẍ<sub>i</sub>'' + ''ÿ'' の加速度を受けるため、運動方程式は {{NumBlk|:|<math> \boldsymbol{M \ddot{x}} + \boldsymbol{C \dot{x}} + \boldsymbol{K x} = - \boldsymbol{M 1} \ddot{y} </math>|{{EquationRef|5.18}}}} となる{{Sfnm|宮本・永野・藤谷・吉村|2014|1p=94|砂子田・伊藤・鄭・平元|2012|2p=98}}。ここで '''1''' は成分が全て 1 の縦ベクトルである。式{{EqNoteN|5.18}}において、'''''C''''' をレイリー減衰で表し、さらに '''''x''''' をモード座標に置き換え、左から転置したモード行列を掛けて整理すれば、下記のような各モード座標についての式になる{{Sfn|宮本・永野・藤谷・吉村|2014|p=94}}。 {{NumBlk|:|<math> \begin{cases} \ddot{q}_1 + 2 \zeta_1 \omega_1 \dot{q}_1 + \omega_1^2 q_1 = - \dfrac{ \boldsymbol{u_1^\top M 1} }{ M_1 } \ddot{y} \\ \ddot{q}_2 + 2 \zeta_2 \omega_2 \dot{q}_2 + \omega_2^2 q_2 = - \dfrac{ \boldsymbol{u_2^\top M 1} }{ M_2 } \ddot{y} \\ \vdots \\ \ddot{q}_n + 2 \zeta_n \omega_n \dot{q}_n + \omega_n^2 q_n = - \dfrac{ \boldsymbol{u_n^\top M 1} }{ M_n } \ddot{y} \\ \end{cases} </math>|{{EquationRef|5.19}}}} さらに、上式の ''ÿ'' の係数を次のように変形する{{Sfn|宮本・永野・藤谷・吉村|2014|p=94}}。 {{NumBlk|:|<math> \dfrac{ \boldsymbol{u_r^\top M 1} }{ M_r } = \dfrac{ \boldsymbol{u_r^\top M 1} }{ \boldsymbol{ u_r^\top M u_r } } = \beta_r </math>|{{EquationRef|5.20}}}} この ''β<sub>r</sub>'' を ''r'' 次の刺激係数と呼び、''r'' 次モードの運動方程式は次のようになる{{Sfnm|宮本・永野・藤谷・吉村|2014|1p=94|平井・水田|2018|2pp=151–152}}。 {{NumBlk|:|<math> \ddot{q}_r + 2 \zeta_r \omega_r \dot{q}_r + \omega_r^2 q_n = - \beta_r \ddot{y} </math>|{{EquationRef|5.21}}}} つまり、''β<sub>r</sub>'' は励振加速度が各モードに対してどのぐらい寄与しているか、あるいは各モードは基礎励振に対してどのぐらい影響を受けやすいかを表している{{Sfnm|宮本・永野・藤谷・吉村|2014|1p=94|砂子田・伊藤・鄭・平元|2012|2p=98}}。さらに、''r'' 次の刺激係数と固有モードの積を ''r'' 次の'''刺激関数'''と呼び、刺激関数の総和は次のように '''1''' に等しい{{Sfn|宮本・永野・藤谷・吉村|2014|p=94}}。 {{NumBlk|:|<math> \beta_1 \boldsymbol{u}_1 + \beta_2 \boldsymbol{u}_2 + \ \cdots \ + \beta_n \boldsymbol{u}_n =\boldsymbol{1} </math>|{{EquationRef|5.22}}}} すなわち、対称の振動系に '''1''' という外力ベクトルが加わったときに各モードが受ける度合いを、刺激関数は表している{{Sfn|宮本・永野・藤谷・吉村|2014|p=94}}。 ==有限要素法による連続体の振動への応用== 多自由度系では、質点または剛体から成る系を想定し、剛性、減衰、慣性などの振動特性も局所的に集中して系内で点在しているというモデルを考えていた{{Sfn|末岡・金光・近藤|2000|pp=5, 103}}。一方で、実際の物体は、物体自体が変形する{{Sfn|下郷・田島|2002|p=115}}。実際の物体は[[連続体力学|連続体]]としての性質を有しており、質量、剛性などは連続的に系内に分布しているモデルとなる{{Sfnm|末岡・金光・近藤|2000|1p=103|下郷・田島|2002|2p=115}}。実際の問題では、多自由度系に近似して取り扱っても十分な場合も多いが、振動時の機械・構造物の各部の変形や応力といったものを知るには連続体として取り扱う必要ある{{Sfn|入江・小林|2006|p=121}}。 連続体の振動の運動方程式は時間と空間に関する[[偏微分方程式]]で記述され、厳密解が得られることは限られる{{Sfn|下郷・田島|2002|p=115}}。実際の複雑な形状の構造物で連続体の振動を扱うには、実用的には[[有限要素法]]という手法が用いられる{{Sfnm|藤田|2016|1p=199|安田|2012|2p=184}}。有限要素法では、対象の連続体を小さな有限要素に分割し、連続体を多自由度系に置き換えて解を計算する{{Sfn|末岡・金光・近藤|2000|p=151}}。有限要素法においても、線形多自由度系の理論にもとづくモード解析手法が強力な効果を発揮し、振動解析が精度良くかつ容易にできるようになる{{Sfnm|長松・長松|2018|1pp=11–12|藤田|2016|2p=199}}。 有限要素法による振動解析では、立てられた振動数方程式の数値計算を行い、まず固有振動数と固有モードを得る<ref name="小松2009">{{Cite book ja-jp |author = 小松 敬治 |title = 機械構造振動学 ―MATLABによる有限要素法と応答解析 |publisher = 森北出版 |url = https://www.morikita.co.jp/books/book/1428 |edition = 第1版 |year = 2009 |isbn = 978-4-627-66611-5 }} pp. 43, 105, 192</ref>。次いで、固有モードの直交性を利用して周波数応答関数を得て、応答解析を行う<ref name="小松2009"/>。もし、定式化された運動方程式が1万自由度だとしたら、解くべき方程式は1万元の2次連立微分方程式となり、コンピュータを用いても計算に長時間を要する{{Sfn|長松・長松|2018|pp=11–12}}。しかし、モード解析手法を用いれば、1自由度系の解を1万回解いて、重ね合わせるだけで解が得られる{{Sfn|長松・長松|2018|pp=12, 89}}。さらに、対象が大規模自由度になったとしても、自由度の分だけ現れるモードを全て計算する必要性もない{{Sfn|末岡・金光・近藤|2000|pp=59–60}}。実用的に興味のある外力振動数を含む次数までモードの重ね合わせでも、十分な精度の振動応答解析が可能となる{{Sfn|末岡・金光・近藤|2000|pp=59–60}}。上記の1万自由度の例えで言えば、1自由度系の解を1万回解く必要もなく、もっと少ない回数の計算で事足りるようになる{{Sfn|長松・長松|2018|p=89}}。これらの長所によって、モード解析手法は有限要素法による振動解析で絶大な威力を発揮し、数十万規模の自由度を扱うような有限要素法計算であっても、モード解析手法の適用によって特段の支障なく計算が可能となる{{Sfn|長松・長松|2018|p=12}}。 ==係数行列が非対称行列の場合== 質量行列 '''''M'''''（式{{EqNoteN|2.2}}）、減衰行列 '''''C'''''（式{{EqNoteN|2.3}}）、あるいは剛性行列 '''''K''''' （式{{EqNoteN|2.4}}）が[[行列の定値性|正定値]]の条件を満たさない場合、すなわち実対称行列ではなく、非対称行列であるとき、その系では不安定振動が起こることがある{{Sfn|日本機械学会（編）|2004|p=70}}。このような条件では、式{{EqNoteN|4.6}}で表される特性方程式の固有値 ''λ'' に、実部が正の固有値が含まれることがありえる{{Sfn|末岡・金光・近藤|2000|p=175}}。固有値に実部を正とする複素数が含まれるとき、時間とともに振幅が大きくなっていく振動が起こる<ref name="安田2001"/>。このようなメカニズムは、[[自励振動]]が起こりえる系で平衡点から振動が成長するか否かを考察する上で基礎となる{{Sfn|日本機械学会（編）|2004|pp=70–71}}。自励振動は1自由度系でも起きる現象だが、係数行列が非対称であることによって引きこされる種類の自励振動は、多自由度系特有のものである{{Sfn|日本機械学会（編）|2004|p=76}}。 例として、次のような2自由度不減衰系を考える<ref name="小野・岡田・山崎">{{Cite journal ja-jp |author = 小野 京右・岡田 徹・山崎 聡 |year = 1994 |title = 自励振動アクチュエータに関する研究 : 第1報,2自由度自励振動系の特性解析 |journal = 日本機械学会論文集 C編 |volume = 60 |issue = 577 |publisher = 日本機械学会 |doi = 10.1299/kikaic.60.2984 |pages = 2985–2986 }}</ref>。 {{NumBlk|:|<math> \begin{pmatrix} m_{1} & 0 \\ 0 & m_{2} \end{pmatrix} \begin{pmatrix} \ddot{x}_{1} \\ \ddot{x}_{2} \end{pmatrix} + \begin{pmatrix} k_{11} & k_{12} \\ k_{21} & k_{22} \end{pmatrix} \begin{pmatrix} x_{1} \\ x_{2} \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \end{pmatrix} </math>|{{EquationRef|7.1}}}} ただし、''k''<sub>12</sub> ≠ ''k''<sub>21</sub> で、剛性行列は非対称行列である<ref name="小野・岡田・山崎"/>。さらに、''k''<sub>12</sub> と ''k''<sub>21</sub> のどちらかが正でどちらかが負であるような異符号の関係にあるとき、固有値は {{NumBlk|:|<math> \lambda = \pm \lambda_r \pm i \lambda_i </math>|{{EquationRef|7.2}}}} という形の複素数となる<ref name="小野・岡田・山崎"/>。''λ<sub>r</sub>'' と''λ<sub>i</sub>'' は {{NumBlk|:|<math> \lambda_r = \frac{1}{2}\sqrt{-\omega_{tr} + \sqrt{\omega_{tr}^2 - \omega_{diff}^2-4\omega_{sk} } } </math>|{{EquationRef|7.3}}}} {{NumBlk|:|<math> \lambda_i = \frac{1}{2}\sqrt{\omega_{tr} + \sqrt{\omega_{tr}^2 - \omega_{diff}^2-4\omega_{sk} } } </math>|{{EquationRef|7.4}}}} で与えられ、ここで、''ω<sub>tr</sub>'' = ''k''<sub>11</sub>/''m''<sub>1</sub> + ''k''<sub>22</sub>/''m''<sub>2</sub>, ''ω<sub>diff</sub>'' = ''k''<sub>11</sub>/''m''<sub>1</sub> − ''k''<sub>22</sub>/''m''<sub>2</sub>, ''ω<sub>sk</sub>'' = (''k''<sub>12</sub>/''m''<sub>1</sub>)(''k''<sub>21</sub>/''m''<sub>2</sub>) である<ref name="小野・岡田・山崎"/>。''λ<sub>r</sub>'' は発散率と呼ばれ、自励振動の強さを表す<ref name="小野・岡田・山崎"/>。 このような係数行列の非対称性によって起きる自励振動の事例は機械振動の中で多く見られ、[[クーロン摩擦]]による摩擦振動や[[滑り軸受]]で起こるオイルホイップなどがある{{Sfn|日本機械学会（編）|2004|pp=70, 74–76}}。 ==出典== {{Reflist|2}} ==参照文献== *{{Cite book ja-jp |author = 末岡 淳男・金光 陽一・近藤 孝広 |title = 機械振動学 |series = 基礎機械工学シリーズ 6 |publisher = [[朝倉書店]] |url = http://www.asakura.co.jp/books/isbn/978-4-254-23706-1/ |edition = 初版 |year = 2000 |isbn = 4-254-23706-5 |ref = {{Sfnref|末岡・金光・近藤|2000}} }} *{{Cite book ja-jp |author = 横山 隆・日野 順市・芳村 敏夫 |title = 基礎振動工学 |publisher = [[共立出版]] |url = https://www.kyoritsu-pub.co.jp/bookdetail/9784320082113 |edition = 第2版 |year = 2015 |isbn = 978-4-320-08211-3 |ref = {{Sfnref|横山・日野・芳村|2015}} }} *{{Cite book ja-jp |author = 長松 昌男・長松 昭男 |title = 実用モード解析入門 |publisher = [[コロナ社]] |url = https://www.coronasha.co.jp/np/isbn/9784339082272/ |edition = 初版 |year = 2018 |isbn = 978-4-339-08227-2 |ref = {{Sfnref|長松・長松|2018}} }} *{{Cite book ja-jp |author = 下郷 太郎・田島 清灝 |title = 振動学 |series = 機械系 大学講義シリーズ 11 |publisher = コロナ社 |url = https://www.coronasha.co.jp/np/isbn/9784339040456/ |edition = 初版 |year = 2002 |isbn = 4-339-04045-2 |ref = {{Sfnref|下郷・田島|2002}} }} *{{Cite book ja-jp |author = 入江 敏博・小林 幸徳 |title = 機械振動学通論 |publisher = 朝倉書店 |url = http://www.asakura.co.jp/books/isbn/978-4-254-23116-8/ |edition = 第3版 |year = 2006 |isbn = 978-4-254-23116-8 |ref = {{Sfnref|入江・小林|2006}} }} *{{Cite book ja-jp |author = 吉川 孝雄・松井 剛一・石井 徳章 |title = 機械の力学 |publisher = コロナ社 |url = https://www.coronasha.co.jp/np/isbn/9784339042733/ |edition = 初版 |year = 1987 |isbn = 978-4-339-04273-3 |ref = {{Sfnref|吉川・松井・石井|1987}} }} *{{Cite book ja-jp |author = 安田 仁彦 |title = 振動工学 ―基礎編 |publisher = コロナ社 |url = https://www.coronasha.co.jp/np/isbn/9784339046243/ |edition = 改訂版 |year = 2012 |isbn = 978-4-339-04624-3 |ref = {{Sfnref|安田|2012}} }} *{{Cite book ja-jp |author = 藤田 勝久 |title = 振動工学 ―振動の基礎から実用解析入門まで |publisher = [[森北出版]] |url = https://www.morikita.co.jp/books/book/3059 |edition = 新装版 |year = 2016 |isbn = 978-4-627-66542-2 |ref = {{Sfnref|藤田|2016}} }} *{{Cite book ja-jp |editor = 日本機械学会 |title = 機械力学 |series = 機械工学便覧 基礎編 α2 |publisher = [[丸善]] |url = https://www.maruzen-publishing.co.jp/item/b297444.html |edition = 第1版 |year = 2004 |isbn = 4-88898-116-7 |ref = {{Sfnref|日本機械学会（編）|2004}} }} *{{Cite book ja-jp |author = 砂子田 勝昭・伊藤 智博・鄭 萬溶・平元 和彦 |title = わかりやすい振動工学 |publisher = 共立出版 |url = https://www.kyoritsu-pub.co.jp/kenpon/bookDetail/9784320081857 |edition = 初版 |year = 2012 |isbn = 978-4-320-08185-7 |ref = {{Sfnref|砂子田・伊藤・鄭・平元|2012}} }} *{{Cite book ja-jp |editor = 宮本 裕司 |author = 宮本 裕司・永野 正行・藤谷 秀雄・吉村 智昭 |title = 建築振動を学ぶ ―地震から免震・制震まで |publisher = 理工図書 |url = https://www.rikohtosho.co.jp/book/1763/ |edition = 初版 |year = 2014 |isbn = 978-4-8446-0824-0 |ref = {{Sfnref|宮本・永野・藤谷・吉村|2014}} }} *{{Cite book ja-jp |author = 平井 一男・水田 洋司 |title = 耐震工学入門 |publisher = 森北出版 |url = https://www.morikita.co.jp/books/book/3318 |edition = 第3版・補訂版 |year = 2018 |isbn = 978-4-627-46454-4 |ref = {{Sfnref|平井・水田|2018}} }} *{{Cite book ja-jp |author = 背戸 一登・丸山 晃市 |title = 振動工学 ―解析から設計まで |publisher = 森北出版 |url = https://www.morikita.co.jp/books/book/1411 |edition = 第1版 |year = 2002 |isbn = 4-627-66451-6 |ref = {{Sfnref|背戸・丸山|2002}} }} {{Good article}} {{DEFAULTSORT:せんけいたしゆうとけいのしんとう}} [[Category:振動工学]]
このページで使用されているテンプレート:
テンプレート:Cite book ja-jp
(
ソースを閲覧
)
テンプレート:Cite journal ja-jp
(
ソースを閲覧
)
テンプレート:Cite web
(
ソースを閲覧
)
テンプレート:EqNoteN
(
ソースを閲覧
)
テンプレート:Good article
(
ソースを閲覧
)
テンプレート:NumBlk
(
ソースを閲覧
)
テンプレート:Reflist
(
ソースを閲覧
)
テンプレート:Sfn
(
ソースを閲覧
)
テンプレート:Sfnm
(
ソースを閲覧
)
線形多自由度系の振動
に戻る。
ナビゲーション メニュー
個人用ツール
ログイン
名前空間
ページ
議論
日本語
表示
閲覧
ソースを閲覧
履歴表示
その他
検索
案内
メインページ
最近の更新
おまかせ表示
MediaWiki についてのヘルプ
特別ページ
ツール
リンク元
関連ページの更新状況
ページ情報