線形多自由度系の振動

振動工学における**線形多自由度系の振動**（せんけいたじゆうどけいのしんどう）は、線形な特性を持ち、さらに2以上の自由度を持つ系で起きる振動である。運動方程式は一般的に連立2階常微分方程式となり、行列およびベクトルで表現される。

線形多自由度系の振動では、固有モードという多自由度系特有の概念が現れ、自由度の数だけ固有モードと固有振動数の組が存在する。固有モードの直交性によって、減衰の無い系であれば固有モードごとの1自由度系の問題に帰着でき、振動解析を容易化できるのが特徴である。この手法を利用した振動解析手法はモード解析と呼ばれる。減衰のある系でも、比例粘性減衰という仮定を導入することによって、同様なことが可能となる。

モード解析手法は、振動実験結果から振動特性を同定するのにも使われる。有限要素法による連続体の振動計算においても、線形多自由度系の理論にもとづくモード解析手法が強力な効果を発揮し、振動解析を容易にする。 

## 背景
振動の問題を分類すると、あるいは振動現象を模した動力学モデルを分類すると、いくつかの視点が存在する。

物体の運動を表すのに必要な座標あるいは変位の数を、自由度という。自由度が1の系を**1自由度系**という。ある1つの質点がばねとダッシュポットを介して地面に固定され、質点が上下方向のみ動く例を考える。三角関数で表現される励振力がこの質点に加わるとき、この質点の運動方程式は次のような形で与えられる。

ここで、*m* は質量、*c* は粘性減衰係数、*k* はばね定数、*f*0 は外力の振幅、*&omega;* は外力の角振動数、*t* は時間である。*x* は静的釣り合い位置からの変位で、上部の "˙" は時間微分 (*d*/*dt*) を示す。式は、振動問題を考える際の最も基礎的な式となる。1自由度系の振動の問題は、振動現象を理解する上で不可欠な様々な概念を内包しており、1自由度系のモデルが振動問題を扱う基本モデルといえる。

一方、自由度が2以上の系は**多自由度系**と呼ばれる。振動問題をモデル化する際には1自由度や2自由度にモデル化すれば解析が容易であり、モデル化時の自由度はできるだけ最小限とするのが大切である。しかし、現実の構造物は複雑で、すべてを1自由度系として扱うことができない。実際の機械や建物などでは、1自由度系のモデルではその特性を説明しきれず、多自由度系としての取り扱いが必要になることも多い。

また、系の構成要素が全て線形であるとき、その系を**線形系**という。系が線形であれば、出力と入力の間には重ね合わせの原理が成り立ち、出力は入力に単純に比例する。振動系が線形であれば、式のように慣性力、減衰力、復元力が、それぞれ加速度、速度、変位に比例する。比例係数が時間変化する場合を除き、線形系の振動問題については理論的手法がほぼ確立できている。

一方で、線形ではない系を**非線形系**という。摩擦、がたつき、大変位時の材料特性などが非線形の要因となる。厳密に考えると、実際の機械や構造物はなんらかの非線形特性を持ち、実際の振動現象のほとんどは非線形系といえる。非線形系の厳密解を得ることはほとんどできず、一般的には非線形系の問題に対しては近似解法や数値計算に頼らざるを得ない。しかし、線形系の振動として取り扱うことによって、十分に多くの問題を解決することもできる。非線形系であっても、安定な平衡状態周りの微小変位運動であれば線形の理論が当てはまる。特に線形系であれば、多自由度系であっても後述のように1自由度系の問題に帰着できるという顕著な特性がある。

## 運動方程式
### 一般の多自由度系
多自由度系では多数の変数や係数を扱うため、行列とベクトルを使って運動方程式を記す。自由度が *n* の線形多自由度系の一般的な運動方程式は、次のように書き表される。

 + \boldsymbol{C \dot{x}} + \boldsymbol{K  x} = \boldsymbol{f}
|}}

ここで、***&#7821;***, ***&#7819;***, ***x***, ***f*** は *n* 次元縦ベクトル、***M***, ***C***, ***K*** は *n* 次正方行列で、以下のように表される。

:
\boldsymbol{\ddot{x}} = 
\begin{pmatrix}
\ddot{x}_1 \\
\ddot{x}_2 \\
\vdots \\
\ddot{x}_n
\end{pmatrix}
, 
\boldsymbol{\dot{x}} = 
\begin{pmatrix}
\dot{x}_1 \\
\dot{x}_2 \\
\vdots \\
\dot{x}_n
\end{pmatrix}
, 
\boldsymbol{x} = 
\begin{pmatrix}
x_1 \\
x_2 \\
\vdots \\
x_n
\end{pmatrix}
, 
\boldsymbol{f} = 
\begin{pmatrix}
f_1 \\
f_2 \\
\vdots \\
f_n
\end{pmatrix}


***x*** は**変位ベクトル**と呼ばれる。その成分の *x* は、静的な釣り合いの位置を原点とした各自由度の変位あるいは一般化座標を表している。***&#7819;*** は**速度ベクトル**、***&#7821;*** は**加速度ベクトル**と呼ばれる。成分の *&#7819;* は *x* の時間1階微分すなわち各自由度の速度を意味し、*&#7821;* は *x* の時間2階微分すなわち各自由度の加速度を意味する。

***f*** は**外力ベクトル**と呼ばれる。成分の *f* は各一般化座標に対応して作用する外力で、一般的には時間の関数である。

***M*** は**質量行列**や**慣性行列**、**質量マトリックス**や**慣性マトリックス**と呼ばれる。***K*** は**剛性行列**や**剛性マトリックス**と呼ばれる。***C*** は**減衰行列**や**減衰マトリックス**と呼ばれる。各成分の *m*, *c*, *k* は質量、粘性減衰係数、剛性（ばね定数）を表す。ここでの *m* は**慣性係数**ともいい、いわゆる質量だけでなく慣性モーメントなども含む。また、*k* は**復元係数**ともいい、通常のばね定数の他に回転ばね定数なども含む。

線形多自由度系の一般基礎式は、線形1自由度系の基礎式と形式は同じで、変位と外力がベクトルに、質量とばね定数と粘性減衰係数が行列に置き換わった式となる。ただし、対象の系の規模が大きくなり、複雑なものとなると、ニュートンの運動法則から運動方程式を導出するのは容易ではない。実際に多自由度系の運動方程式を立てる際は、エネルギーのスカラー量から形式的に運動方程式を導けるラグランジュの運動方程式が便利で、間違いを犯しにくく、多用されている。

### 2自由度系の例

建築物の例では、複数の階を持つ多層構造物が多自由度系の問題となる。柱と床から構成されるラーメン構造の建物が振動する場合を考える。柱は床に比べて柔らかいので、構造物が揺れるとき各階の床は水平方向に揺れ、柱は水平方向のばねとして働くと見なせる。これは、水平方向のみに動く2つの質点をばねで連結した2自由度系モデルと等価となる。2層構造物における、1層目の床質量を *m*1、2層目の床質量を *m*2、1層目の水平変位を *x*1、2層目の水平変位を *x*2、基礎と1層目の間の等価ばね定数を *k*1、1層目と2層目の間の等価ばね定数を *k*2 とすると、このモデルの運動方程式は

となる。行列表示すれば、

となり、この式では、
\begin{pmatrix}
m_{1} & 0 \\
0 & m_{2}
\end{pmatrix}
が質量行列、
\begin{pmatrix}
k_{1}+k_{2} & -k_{2} \\
-k_{2} & k_{2}
\end{pmatrix}
が剛性行列である。

並進運動と回転運動が組み合わさった2自由度系の例として、自動車の上下振動とピッチング振動の簡易モデルがある。自動車を側面から見て、上下方向変位とピッチング回転だけの動きを考える。前位側のタイヤとサスペンションを1つのばねと見なして、後位側のタイヤとサスペンションも同様に1つのばねと見なし、剛体を前位と後位を2つのばねが支えているモデルを考える。前位側と後位側のばねのばね定数をそれぞれ *k*1、*k*2 とする。剛体の重心位置から前位側のばねまでの距離を *l*1、前位側のばねまでの距離を *l*2 とする。剛体の質量を *m*、重心位置周りのピッチング方向の慣性モーメントを *I<sub>G</sub>* とする。この剛体の重心の上下運動 *x* とその周りのピッチング運動 *&theta;* の運動方程式は次のようになる。

### 連成項
自動車の上下振動とピッチング振動の簡易モデルである式において、もし *k*1*l*1 = *k*2*l*2 であれば、*x* と *&theta;* は互いに独立した運動方程式となる。この条件では運動方程式は

となり、2つの式はそれぞれ単独で解くことができる。すなわち、*x* の振動と *&theta;* の振動が、互いに干渉すること無しに独立して起こる状態になっている。このような状態を**非連成**という。運動方程式が単独では解けない連立方程式になっているとき、方程式は**連成**しているといい、式の (*k*1*l*1 &minus; *k*2*l*2)*x* と (*k*1*l*1 &minus; *k*2*l*2)*&theta;* は**連成項**と呼ばれる。*k*1*l*1 = *k*2*l*2 のときに、これら連成項が 0 になり、*x* と *&theta;* が独立した運動になる。連成状態にある系を**連成系**、逆に各自由度の運動が完全に独立している系を**非連成系**と呼ぶ。

一般に、質量行列 ***M***（式）、減衰行列 ***C***（式）、剛性行列 ***K***（式）における非対角成分の存在は、振動が連成していることを示している。もし系が非連成系であれば、これら全ての行列は、非対角成分が全て 0 の対角行列となる。特に、質量行列が非対角成分によって振動が連成していることを**動連成**や**動的連成**といい、質量行列の非対角成分を**動連成項**という。また、剛性行列が非対角行列であれば**静連成**や**静的連成**といい、剛性行列の非対角成分を**静連成項**という。多自由度系の問題は概して複雑で、個々のパラメータが結果に与える影響の見通しを立てることが難しい。多自由度系では一般的に連成が存在し、これが多自由度系の解析を難しいものにしている。動力学的な設計を行う上では、系を非連成化して影響をわかりやすくすることが有効となる。

## 不減衰自由振動
### 振動数方程式
継続的な外力が作用せず、外乱だけが与えられて起こる振動を自由振動という。自由振動の振動は、系自体が持つ特性によって決まり、系の動特性を知る上で重要な振動の形態である。自由振動は減衰要素が存在しない場合と存在する場合に分かれる。減衰要素と励振力が存在しない場合の式は

 + \boldsymbol{K x} = \mathbf{0}
|}}

となる。ここで、**0** は下記のようなゼロベクトルである。

:
\mathbf{0} = 
\begin{pmatrix}
0 \\
0 \\
\vdots \\
0
\end{pmatrix}


さらに、***M*** と ***K*** は正定値であると仮定する。すなわち、***M*** と ***K*** は成分が全て実数の対称行列で、それらの二次形式は常に正となる。実際に、質量行列と剛性行列が正定値行列であることは線形多自由度系の一般的な特徴であり、多くの振動系でこの仮定は満たされる。この条件を満たす *M* と *K* を前提にすれば、振動系の具体的な構成に依存しない一般性の高い議論を展開できる。

式は定数係数の線形常微分方程式であるため、その解法に従って解を

という形式で仮定できる。ここで *e* はネイピア数、*&lambda;* は未知定数、 ***u*** は *n* 個の未知定数 *u* から成る下記のような縦ベクトルである。

:
\boldsymbol{u} = 
\begin{pmatrix}
u_1 \\
u_2 \\
\vdots \\
u_n
\end{pmatrix}


式に対して仮定として与えられる解には、他に、最初から単振動を仮定して複素指数関数や三角関数の形式もある。

式を式に代入して整理すると下記のような式になる。

この式が恒等的に成り立つには下記の条件が満たされる必要がある。

数学的には、この形式の方程式は一般化固有値問題として知られ、*&lambda;*2 は固有値、***u*** は固有ベクトルと呼ばれる。***u*** = **0** であれば式の条件は満たされるが、これは最初の釣り合いの位置からそのまま静止しているだけ状態を意味する解である。よって、ここでは ***u*** = **0** は興味の対象外で、それ以外の式を満たす *&lambda;* と ***u*** について知りたい。***u*** &ne; **0** で、なおかつ式が満たされる条件は、*U* の係数行列 *&lambda;*2 ***M*** + ***K*** の逆行列が存在しないことである。したがって、係数行列の行列式が 0 であればこの条件が満たされる。したがって、*&lambda;* が、

を満たすとき、式が ***u*** = **0** 以外の解を持つ。ここで、det( ) は行列式を表す。式の行列式を展開すると、 *&lambda;*2 についての *n* 次多項式になる。したがって、原理的には *&lambda;*2 の値を求めることができる。ただし、この多項式を解析的に解くことができるのはせいぜい2自由度あるいは3自由度までで、それ以上の自由度の系になると数値解析で固有値を計算する。ヤコビ法などを使って数値計算するときは、コレスキー分解を使い、一般化固有値問題形式のを標準的な固有値問題の形へ変換する。

上記のとおり、***M*** と ***K*** は正定値であると仮定した。このとき、式の解は全て負の実数となる（解を複素指数関数と仮定した場合は全て正の実数）。したがって、*n* 個の *&lambda;*2 の解を

とおくことができ、平方根を取って 2*n* 個の*&lambda;* の値 

が得られる。ここで *j* は虚数単位である。*n*　個の *&omega;* は**固有角振動数**または**固有円振動数**と呼ばれ、値が小さいものから順に**1次**、**2次**、…、***n'' 次'''の固有角振動数と呼ぶ。特に、最も値が小さい1次の固有角振動数は'''基本振動数'''と呼ばれる。式は角振動数を求める式であるため'''振動数方程式'''と呼ばれる。あるいは、式は ***M*** や ***K*** といった系の特性によって構成される式であることから'''特性方程式'''とも呼ぶ。

===固有モード===
得られた固有値あるいは固有角振動数の値を式に代入すれば、'''0''' 以外の ***u*** の解が得られる。ここで得られる ***u*** の各成分は一意な値を持たず、定まるのは各成分の互いの比 ''u''<sub>1</sub> : ''u''<sub>2</sub> : &hellip; : ''u<sub>n</sub>'' だけである。一つの固有値あるいは固有角振動数に対応して一つの ***u*** が定まり、***u*** は ''n'' 個存在する。このような固有値と固有ベクトルの組は固有ペアと呼ばれる。''r'' 次の固有角振動数 ''&omega;<sub>r</sub>'' に対応する ***u*** を ***u**<sub>r</sub>'' を表現すれば、 ***u*'r* は



というベクトルである。ただし、上記のとおり、各成分 *u*<sub>1*r*</sub>, *u*<sub>2*r*</sub>, &hellip;, *unr* は互いの比の値を表している。


式で表されるベクトルが式における固有ベクトルであり、振動工学では*'固有モード**、**振動モード**、**固有振動モード**、**基準振動モード**、**モードベクトル**などと呼ぶ。''n'' 自由度系には ''n'' 個の固有角振動数があり、固有角振動数それぞれに対応する形で ''n'' 個の固有モードが存在している。各自由度の振幅比を決める固有モードは、固有角振動数が「振動の速さ」を表しているのに対して、「振動の形」を表していると言える。自由度の数の分だけ固有角振動数が存在し、それら固有角振動数に対応して固有モードが存在していることが、線形多自由度系の特有な性質といえる。

各々の固有モードを、対応する固有角振動数が小さい順に1次固有モード、2次固有モード、…、''r'' 次固有モード、…、''n'' 次固有モードと呼ぶ。特に、基本振動数（1次の固有角振動数）に対応する固有モードは**基本モード**と呼ばれる。固有モード  ***u*'r* を次数が低い順に並べて作る、下記のような行列を*'モード行列**や**モードマトリックス**という。



===固有モードの直交性と正規化===
多自由度系の固有モードには**直交性**という重要な性質がある。''r'' 次の固有モード ***u*'r* と *s* 次の固有モード ***u**s* について考える。ここで、*r* と *s* は任意だが、*r* &ne; *s* である。***M*** と ***K*** は上記のとおり対称行列とする。このとき、  ***u**r*, ***u**s*, ***M***, ***K*** には次のような関係がある。





ここで、***u**r*<sup>&top;</sup> は ***u**r* の[[転置行列]]を表している。これらの式は ***u**r* と ***Mu**s* の[[内積]]および ***u**r* と ***Mu**s* の内積が零であることを示しており、***M*** と ***K*** に関して ***u**r* と ***u**s* が[[直交]]していることを意味している。これが質量行列および剛性行列を介した固有モードの直交性である。別の言い回しでは ***u**r* と ***u**s* が一般直交性を有しているともいう。

固有モードの直交性の物理的な意味は、次数の異なる固有モードの振動の間で[[力学的エネルギー]]の移動が全く起きないことを表している。ある固有モードの振動において、[[運動エネルギー]]と復元力による[[ポテンシャルエネルギー]]は時間的に変化しているが、それらの和の力学的エネルギーは時間に依らず一定に保たれている。そのため、自由振動中にある固有モードの振動が他の固有モードに移り変わったり、他の固有モードが新たに誘起されたりすることはない。

一方で、*r* = *s*、すなわち同じ次数同士の固有モードの場合は、上記の式でも左辺は 0 とはならず、ある定数となる。これらの定数を *Mr* と *Kr* と表せば、





と表される。固有モード ***u**r* は成分間の比だけを持ち、定まった値を持たないので、*Mr* と *Kr* の値もこの段階では定まらない。固有モードの絶対値が決まった後に、*Mr* と *Kr* の値も定まる。

定数 *Mr* と *Kr* は正の値であり、それぞれを（*r* 次の）*'モード質量**および**モード剛性**という。モード行列を使って式と式を表せば、





となる。モード質量とモード剛性は、同じ次数の固有角振動数 ''&omega;<sub>r</sub>'' と下記のような関係がある。


</math>|}}

固有モードは比の関係を定めるベクトルであった。絶対的な大きさを持たない固有モードに大きさを定める方法として、次のような方法がある。
#下記のように、式の右辺 ''M<sub>r</sub>'' の値が 1 となるように定める。ここで、***&#363;*'r* は、等式を満たすに大きさが決定された固有モードを意味している。
#:_{r}^\top \boldsymbol{M} \boldsymbol{\overline{u}}_{r} = 1 
</math>|}}
#下記のように、同じ次数の固有モード同士の[[内積]]が 1 となるように定める。
#:_{r}^\top \boldsymbol{\overline{u}}_{r} = 1 
</math>|}}
#固有モードの成分の内、[[絶対値]]が最大のものを 1 とおいて定める。
#1番目の自由度に対応する固有モード成分 (*u*<sub>1*r*</sub>) を 1 とおいて定める。

以上のように固有モードの大きさを一意に定めることを*'正規化**と呼び、正規化された固有モード ***&#363;*'r* を*'正規固有モード**や**正規化モード**と呼ぶ。1番目の手法による正規固有モードは、特に**質量正規固有モード**や***M*-正規固有モード**という名で呼ばれる。式を満たすようにしておくと利点が多く、質量正規固有モードはよく使われる。正規固有モード ***&#363;*'r* を並べて作る下記のような行列を、*'正規モード行列**という。

 = (\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>|}}

質量正規固有モードを採用すると、式に対して剛性行列は固有角振動数と

_r^\top \boldsymbol{K} \boldsymbol{\overline{u}}_r = \omega_r^2
</math>|}}

という関係を持つ。

===モード座標と一般解===
固有モード（固有ベクトル）は互いに直交なので、これらを使った[[線形結合]]で変位ベクトル ***x*** を表すことができる。つまり、固有モード ***U*** を用いて、***x*** を

_{1} + q_2 \boldsymbol{\overline{u}}_{2}+ \dots + q_n \boldsymbol{\overline{u}}_{n} 
=  \boldsymbol{\overline{U} q}
</math>|}}

と表すことができる。ここで、''q'' は時間の未知関数で、**モード座標**、**基準座標**、**正規座標**、**規準座標**、**主座標**などと呼ばれる。特に正規化した固有モードを使っているモード座標を指して、正規座標と呼ぶこともある。***q*** は ''q<sub>r</sub>'' を 1 次から ''n'' 次まで並べた次のような縦ベクトルである。

:<math>
\boldsymbol{q} = 
\begin{pmatrix}
q_1 \\
q_2 \\
\vdots \\
q_n
\end{pmatrix}
</math> 

モード座標 ***q*** に対して、元の座標 ***x*** を**物理座標**と呼ぶ。は、正規固有モードを介して物理座標をモード座標に変換していることを意味する。モード座標 ***q*'r* は、***x*** に対して *r* 次固有モード***&#363;**r* が寄与する度合いを表しているとも言える。式は展開定理とも呼ばれる。

式を運動方程式へ代入して、左から転置したモード行列 ***U***<sup>&top;</sup> を掛けると次のようになる。

 + \boldsymbol{U^\top K U q} = \mathbf{0}
</math>|}}

固有モードの直交性（式と式）を上式に当てはめると、



というような *q* に関する *n* 個の運動方程式が得られる。それぞれの式の両辺を *Mr* で割ると、下記のような形になる。




式は非連成化されており、各式はそれぞれ独立している。そのため、一つ一つの式は1自由度系の不減衰自由振動と同じであるから、*q* の各解は以下のようになる。



ここで、*cr* と *&theta;r* (*r* = 1, 2, &hellip;, *n*) は初期条件で決まる定数で、*cr* が[[振幅]]、*&theta;r* が初期[[位相]]を意味する。***x*** の一般解は、式を式に代入して

_{r} C_r \sin (\omega_r t+ \theta_r)
</math>|}}

と得られる。すなわち、*n* 自由度系における各自由度の自由振動は、*n* 個の[[調和振動]]の重ね合わせ（和）となっており、それら調和振動のそれぞれの振動数は 1 次から *n* 次までの固有角振動数となっている。固有モードは、対応する固有角振動数を持つ調和振動成分の各自由度間の振幅比を定めている。位相角を陽に表さずに、余弦関数と正弦関数の和

_{r} (a_r \cos \omega_r t + b_r \sin \omega_r t)
</math>|}}

や、指数の正負が異なる複素指数関数の和

_{r} (d_r e^{i \omega_r t} + g_r e^{- i \omega_r t})
</math>|}}

といった形で一般解を表すこともできる。*ar*, *br*, *dr*, *gr* も初期条件で決まる定数で、式と式は式へ式変形可能な同値な式である。

以上のような、物理座標をモード座標へ変換し、固有モードの直交性を利用して多自由度系の問題を1自由度系の重ね合わせの問題に帰着させ、解析を行う手法を[[モード解析]]（特にモード重畳法）という。後述のように、モード解析は特に[[有限要素法]]へ適用することで有効性を発揮する。実物の振動特性を求める上でもモード解析の理論を適用することで各特性を[[システム同定|同定]]でき、この手法は[[実験モード解析]]と呼ばれる。

==減衰自由振動==
===一般の粘性減衰===

減衰行列 ***C*** が存在する線形多自由度系の振動について考える。外力が無い場合の *n* 自由度系の運動方程式は以下のようになる。

 + \boldsymbol{C \dot{x}} + \boldsymbol{K  x} = \boldsymbol{0}
</math>|}}

基礎の上でばねと減衰器が一組となって質点と連結し、直列に連なった3自由度減衰系の典型的な例では、***M***, ***C***, ***K*** は次のような行列となる。





式の解を式と同じように仮定し、式を式に代入して整理すると、



となり、***u*** = *'0** 以外の解を持つという条件から



という特性方程式が得られる。式を展開すると、''&lambda;'' に関する 2''n'' 次多項式となる。この多項式を解いて得た解 ''&lambda;''<sub>1</sub>, ''&lambda;''<sub>2</sub>, &hellip;,  ''&lambda;''<sub>2''n''</sub> を式に代入し、***u*** を定め、***u*** に適当な正規化を行う。正規化された各ベクトルを ***&#363;***<sub>1</sub>, ***&#363;***<sub>2</sub>, &hellip;, ***&#363;***<sub>2''n''</sub> と表すと、式の一般解は

_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>|}}

となる。ここで、''a''<sub>1</sub>, ''a''<sub>2</sub>, &hellip;, ''a''<sub>2''n''</sub> は初期条件で決まる定数である。2''n'' 個の ''&lambda;'' は、解が[[減衰振動]]であれば ''n'' 組の[[複素共役|互いに共役な複素数]]になる。このとき、固有モードも複素数となり、**複素固有モード**と呼ばれる。

===比例粘性減衰===
式のように一般的な減衰行列 ***C*** が運動方程式に存在する場合、正規モード行列によって ***M*** と ***K*** は対角化できるが、***C*** も同時に対角化することはできない。そのため、不減衰自由振動で可能だったモード座標に変換しての非連成化が不可能となり、モード解析の利点を活かすことができなくなる。そこで、***C*** が下記のような**比例粘性減衰**として与えられると仮定し、実際の振動解析が行われることも多い。



ここで、''&alpha;'' と ''&beta;'' は定数で、**比例減衰定数**と呼ばれる。比例粘性減衰を有する系を**比例粘性減衰系**などと呼ぶ。比例粘性減衰が成り立つと仮定すれば減衰行列を対角化できる。

式の形で与えられる比例粘性減衰は、特に**レイリー減衰**や**レイリー型減衰**と呼ばれる。''&alpha;**M*** のみで仮定されるものは'''質量比例型減衰'''、''&beta;'''K*** のみで仮定されるものは**剛性比例型減衰**などと呼ぶ。

実際の減衰が比例粘性減衰になっていることはまれであり、比例粘性減衰はあくまでも近似的なものである。しかし、比例粘性減衰の仮定を導入することで、不減衰系と同じ取り扱いが可能となり、モード解析の手法が適用可能になる。もし減衰が全体に分布しているような構造であれば、適当な比例減衰定数を設定すれば、実際の現象を実用問題ないレベルで再現できるという一定の妥当性もある。減衰は摩擦・材料減衰・流体粘性など様々な要因で起こるため、そもそも減衰の適切な定式化自体が難しいといった事情もある。式のように速度の比例定数として与えられる一般の粘性減衰も、多種多様な発生機構によって減衰が起きるという実情に起因して厳密には成立しない。

式を式に代入した場合、

 + (\alpha \boldsymbol{M} + \beta \boldsymbol{K}) \boldsymbol{\dot{x}} + \boldsymbol{K  x} = \boldsymbol{0}
|}}

という運動方程式になる。解を式のように仮定して上式に代入し、整理すると、

となる。さらに、

とおけば、式は

となる。式は、（複素指数関数や三角関数で解を仮定したときの）不減衰振動の固有値問題における *&lambda;* を *&gamma;* に置き換えただけの式になる。したがって、比例粘性減衰を仮定した減衰振動の固有モードは、同一の質量行列と剛性行列を有する不減衰振動の固有モードと同じである。

一方、比例粘性減衰を仮定した減衰振動の固有角振動数は、同一の質量行列と剛性行列の不減衰振動の固有角振動数よりも小さくなる。モード行列 ***U*** で減衰行列を対角化すると

となる。この行列の成分 *&alpha;M<sub>r</sub> + &beta;K<sub>r</sub>* (*r* = 1, 2, &hellip;, *n*) を**モード減衰**や**モード減衰係数**と呼び、 *C<sub>r</sub>* などで表す。ここで、*M<sub>r</sub>* はモード質量、*K<sub>r</sub>* はモード剛性である。さらに、

という量を導入して、これでモード減衰を割った量

を**モード減衰比**と呼ぶ。多自由度系の減衰系では固有モードごとに減衰の効果が異なっており、モード減衰比が固有モードごとの減衰の効果の程度を表している。*&zeta;<sub>r</sub>* &ge; 1 ならば過減衰の状態であり、その固有モードの振動は起こらない。*&zeta;<sub>r</sub>* d,r* を減衰固有角振動数と呼ぶ。以上のように、線形1自由度系の減衰振動の考え方が固有モードごとの振動にも当てはまる。正規座標へ変換を行い、モード減衰比と固有角振動数を用いて運動方程式を表すと、下記のように表現できる。



*&zeta;r* < 1 であれば、*qr* の一般解は積分定数を *ar* と *br* として、



となる。物理座標 ***x*** へ逆変換すれば、比例粘性減衰系の自由振動の一般解は下記のようになる。

_{r} e^{-\zeta_{r} \omega_{r} t}(a_r \cos \omega_{d,r} t + b_r \sin \omega_{d,r} t)
</math>|}}

==強制振動==
===不減衰系===
任意の励振力を受ける減衰の無い系の運動方程式は以下のように表される。

 + \boldsymbol{K  x} = \boldsymbol{f}
</math>|}}

例えば、励振力が角振動数 *&Omega;* の余弦関数で与えられるとすれば、

 + \boldsymbol{K  x} = \boldsymbol{f_a} \cos \mathit{\Omega} t
</math>|}}

となる。ここで、***f<sub>a</sub>*** は、以下のような各自由度に対する励振力の振幅値の縦ベクトルである。

:<math>
\boldsymbol{f_a} = 
\begin{pmatrix}
f_{a,1} \\
f_{a,2} \\
\vdots \\
f_{a,n}
\end{pmatrix}
</math>

式に対して特解を ***x*** = ***u*** cos*&Omega;t* と仮定し、式に代入すれば



となるので、***u*** および ***x*** は、





となる。ここで "<sup>&minus;1</sup>" は[[逆行列]]を意味する。したがって、逆行列 (***M*** &minus; *&Omega;*<sup>2</sup>***K***)<sup>&minus;1</sup> を計算すれば式から ***x*** の値が分かる。しかし、この逆行列を解析的に解くことは困難で、数値計算を行うにしても自由度の数が増えると膨大な計算量になる。もし励振力の角振動数 *&Omega;* を変えると、そのたびに逆行列を計算する必要がある。そのため、実際に ***x*** の解を求めるために行われるのは、下記のようなモード解析による手法である。

解を得るために、***x*** を複素数に拡張し、励振力 ***f*** を[[複素指数関数]] ***f<sub>a</sub>**e j&Omega;t* の形で与えるとする。この場合、計算して解が得られた後に実部あるいは虚部を取ることで、実際の解が得られる。運動方程式の右辺を *'0** としたときのモード行列 ***U*** が、事前に求められているとする。モード座標への変換式を運動方程式へ適用して、左から ***U***<sup>&top;</sup> を掛け、式と式の対角化を適用する。すると、励振力を受ける不減衰系の運動方程式は



という独立・非連成の ''n'' 個の運動方程式に帰着する。ただし、右辺の ''F<sub>r</sub>'' は次のような値である。



式は線形1自由度系と同じなので、[[強制振動]]を表す特解は、



となる。式を使ってモード座標を物理座標へ逆変換し、固有角振動数を使って整理すれば強制振動の解は次のようになる。



式から分かることは、励振力の角振動数 ''&Omega;'' が系の固有角振動数 ''&omega;''<sub>1</sub>, ''&omega;''<sub>2</sub>, &hellip;, ''&omega;<sub>n</sub>'' のどれかに近いとき、係数の[[極限]]が 0 となってその振動成分が極めて大きくなり、[[共振]]が起こるという点である。''&Omega;'' が固有角振動数のいずれかに一致するとき、***x*** の振幅は無限大に発散する。

===比例粘性減衰系===
系が比例粘性減衰系であれば、励振力を受ける場合でもモード座標変換によって独立した1自由度系に分解でき、大自由度の問題も扱いが容易になる。一般の運動方程式に、レイリー減衰の式およびモード座標への変換式を代入する。左から転置した質量正規固有モードの正規モード行列 ***&#362;***<sup>&top;</sup> を掛けて対角化する。さらに、2''&zeta;<sub>r</sub>&omega;<sub>r</sub>'' = ''&alpha;'' + ''&beta;&omega;<sub>r</sub>''<sup>2</sup> という関係を用いれば、運動方程式は下記のような独立・非連成の ''n'' 個の方程式に帰着する。

_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>|}}

励振力 ***f*** が ***f<sub>a </sub>***cos''&Omega;t'' で与えられるとする。線形1自由度系と同様の解法によって、正規座標 ''q<sub>r</sub>'' の強制振動の解が下記のように求まる。 

_{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>|}}



解を複素数で表現した場合は、

_{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>|}}

となり、物理座標の解は

{ \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>|}}

となる。''&Omega;''/''&omega;<sub>r</sub>'' は強制振動比などと呼ばれる。通常の固有モードの場合は次式となる。

{ 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>|}}

===周波数応答関数===
強制振動の特解から、[[周波数応答関数]]あるいは[[伝達関数]]が求められる。周波数応答関数は、振動系への入力が調和振動の形で与えられるときの、出力と入力の振幅比および位相差を周波数の関数として表す。比例減衰系において、角振動数 ''&Omega;'' の調和外力が ''q'' 番目の自由度のみに加わるとする。このときの ''p'' 番目の自由度の振動応答は



と表される。この場合の周波数応答関数 ''G''(''&Omega;'') は[変位]/[力]の単位を持ち、[[機械的コンプライアンス|コンプライアンス]]とも呼ばれる。''q'' 番目加振・''p'' 番目応答の周波数応答関数と ''p'' 番目加振・''q'' 番目応答の周波数応答関数は、[[マクスウェル・ベティの相反作用の定理|相反定理]]によって互いに等しく、''G<sub>pq</sub>''(''&Omega;'') = ''G<sub>qp</sub>''(''&Omega;'') の関係である。不減衰系の場合は、周波数応答関数は次式となる。




これらの周波数応答関数によって、系の周波数特性が把握できる。線形多自由度系の場合は、固有角振動数と固有モードを求めれば、調和外力による強制振動に対する周波数特性も同時に理解することができる。横軸に振動数や強制振動比を取り、縦軸に調和外力に対する振幅比や位相差を図示したものを応答曲線や共振曲線と呼ぶ。縦軸を[[デシベル]]にして図示したものは特に[[ボード線図]]と呼ばれる。系の振動特性を振動実験から同定する実験モード解析では、実験値に上記の周波数応答関数で[[カーブフィッティング]]して系の各種パラメータを推定する。

上述のように、励振力の振動数が系の固有角振動数に一致すると、振動は極大化する。一方、多自由度系では共振して振幅が発散する現象だけでなく、励振力が特定の角振動数のときに周波数応答関数が極小になる現象もある。このような現象を[[反共振]]と呼ぶ。ただし、共振が系全体で起こるのに対して、反共振は一組の加振点・応答点ごとにしか起こらない。周波数応答関数の絶対値を縦軸にした応答曲線上では、共振点は曲線の鋭い山のように現れ、反共振は曲線の鋭い谷のように現れる。反共振は、代表的な[[制振|制振器]]である[[動吸振器]]の原理として活用される。振動を抑えたい対象にばね・ダンパを介して付加質量を取り付けることによって振動を抑制でき、地震や強風に対する建築構造物の防振や回転機械の防振などに使われる。

===刺激係数===
**刺激係数**は、[[地震]]のように構造物の基礎が揺れている場合に、その励振が各モードに対してどのぐらい強く寄与するのかを表す。構造物を多質点系でモデル化して、基礎の変位を ''y'' として、各質点の基礎からの相対変位を ''x'' とする。各質点は　''&#7821;<sub>i</sub>'' + ''&#255;'' の加速度を受けるため、運動方程式は

 + \boldsymbol{C \dot{x}} + \boldsymbol{K  x} = - \boldsymbol{M 1} \ddot{y}
</math>|}}

となる。ここで **1** は成分が全て 1 の縦ベクトルである。式において、***C*** をレイリー減衰で表し、さらに ***x*** をモード座標に置き換え、左から転置したモード行列を掛けて整理すれば、下記のような各モード座標についての式になる。



さらに、上式の ''&#255;'' の係数を次のように変形する。



この ''&beta;<sub>r</sub>'' を ''r'' 次の刺激係数と呼び、''r'' 次モードの運動方程式は次のようになる。



つまり、''&beta;<sub>r</sub>'' は励振加速度が各モードに対してどのぐらい寄与しているか、あるいは各モードは基礎励振に対してどのぐらい影響を受けやすいかを表している。さらに、''r'' 次の刺激係数と固有モードの積を ''r'' 次の**刺激関数**と呼び、刺激関数の総和は次のように **1** に等しい。



すなわち、対称の振動系に **1** という外力ベクトルが加わったときに各モードが受ける度合いを、刺激関数は表している。

==有限要素法による連続体の振動への応用==
多自由度系では、質点または剛体から成る系を想定し、剛性、減衰、慣性などの振動特性も局所的に集中して系内で点在しているというモデルを考えていた。一方で、実際の物体は、物体自体が変形する。実際の物体は[[連続体力学|連続体]]としての性質を有しており、質量、剛性などは連続的に系内に分布しているモデルとなる。実際の問題では、多自由度系に近似して取り扱っても十分な場合も多いが、振動時の機械・構造物の各部の変形や応力といったものを知るには連続体として取り扱う必要ある。

連続体の振動の運動方程式は時間と空間に関する[[偏微分方程式]]で記述され、厳密解が得られることは限られる。実際の複雑な形状の構造物で連続体の振動を扱うには、実用的には[[有限要素法]]という手法が用いられる。有限要素法では、対象の連続体を小さな有限要素に分割し、連続体を多自由度系に置き換えて解を計算する。有限要素法においても、線形多自由度系の理論にもとづくモード解析手法が強力な効果を発揮し、振動解析が精度良くかつ容易にできるようになる。

有限要素法による振動解析では、立てられた振動数方程式の数値計算を行い、まず固有振動数と固有モードを得る。次いで、固有モードの直交性を利用して周波数応答関数を得て、応答解析を行う。もし、定式化された運動方程式が1万自由度だとしたら、解くべき方程式は1万元の2次連立微分方程式となり、コンピュータを用いても計算に長時間を要する。しかし、モード解析手法を用いれば、1自由度系の解を1万回解いて、重ね合わせるだけで解が得られる。さらに、対象が大規模自由度になったとしても、自由度の分だけ現れるモードを全て計算する必要性もない。実用的に興味のある外力振動数を含む次数までモードの重ね合わせでも、十分な精度の振動応答解析が可能となる。上記の1万自由度の例えで言えば、1自由度系の解を1万回解く必要もなく、もっと少ない回数の計算で事足りるようになる。これらの長所によって、モード解析手法は有限要素法による振動解析で絶大な威力を発揮し、数十万規模の自由度を扱うような有限要素法計算であっても、モード解析手法の適用によって特段の支障なく計算が可能となる。

==係数行列が非対称行列の場合==
質量行列 ***M***（式）、減衰行列 ***C***（式）、あるいは剛性行列 ***K**' （式）が[[行列の定値性|正定値]]の条件を満たさない場合、すなわち実対称行列ではなく、非対称行列であるとき、その系では不安定振動が起こることがある。このような条件では、式で表される特性方程式の固有値 *&lambda;* に、実部が正の固有値が含まれることがありえる。固有値に実部を正とする複素数が含まれるとき、時間とともに振幅が大きくなっていく振動が起こる。このようなメカニズムは、[[自励振動]]が起こりえる系で平衡点から振動が成長するか否かを考察する上で基礎となる。自励振動は1自由度系でも起きる現象だが、係数行列が非対称であることによって引きこされる種類の自励振動は、多自由度系特有のものである。

例として、次のような2自由度不減衰系を考える。



ただし、*k*<sub>12</sub> &ne; *k*<sub>21</sub> で、剛性行列は非対称行列である。さらに、*k*<sub>12</sub> と *k*<sub>21</sub> のどちらかが正でどちらかが負であるような異符号の関係にあるとき、固有値は



という形の複素数となる。*&lambda;r* と*&lambda;i* は




で与えられ、ここで、*&omega;tr* = *k*<sub>11</sub>/*m*<sub>1</sub> + *k*<sub>22</sub>/*m*<sub>2</sub>, *&omega;diff* = *k*<sub>11</sub>/*m*<sub>1</sub> &minus; *k*<sub>22</sub>/*m*<sub>2</sub>, *&omega;sk* = (*k*<sub>12</sub>/*m*<sub>1</sub>)(*k*<sub>21</sub>/*m*<sub>2</sub>) である。*&lambda;r'' は発散率と呼ばれ、自励振動の強さを表す。

このような係数行列の非対称性によって起きる自励振動の事例は機械振動の中で多く見られ、クーロン摩擦による摩擦振動や滑り軸受で起こるオイルホイップなどがある。

## 出典

## 参照文献
*
*
*
*
*
*
*
*
*
*
*
*
*