係数励振

係数励振（けいすうれいしん、英: parametric excitation）とは、系の係数（パラメータ）が周期的に変化することで起こる振動現象である。ばね-質量系の運動方程式でいえば、質量やばね定数、減衰係数などの通常では定数とされる係数が周期的に変化するような場合に発生する。パラメトリック励振、パラメータ励振、パラメトリック発振などとも呼ぶ。多数の応用があり、ごく一例としては、電気回路分野における発振器の原理としての応用がある。 遊具のブランコの一人乗りの揺らし方は、係数励振の例である。係数励振系の運動方程式や回路方程式は、マシュー方程式、ヒル方程式の形式に帰着できる場合が多い。 係数励振のメカニズムを利用した身近な例としては、遊具のブランコの1人乗りの揺らし方がある。（詳細は#ブランコの係数励振を参照） 電気回路の例としては、パラメトロンがある。これは、回路のインダクタンスを周期変化させて係数励振振動を生み出すものである。初期条件の差による逆位相の係数励振振動を合わせてを生みだし、2種類の振動を論理素子として利用する。 係数励振が害をもたらす事例としては、鉄道車両におけるパンタグラフの架線からの離線現象、フック形軸継手の不安定ねじり振動などがある。 古くからの係数励振の実験としては、1831年のマイケル・ファラデーによるファラデー波の実験や、1859年のメルデ(Franz Emile Melde)による音叉と弦の実験などが知られている。 公園などにある遊具のブランコは係数励振を原理として利用している。このブランコの動きについて、以下のような物理モデルが考えられる。まず、一般的な振り子の運動方程式を角運動量より導くと次のようになる。 d ( m l 2 θ ˙ ) d t = − m g l sin ⁡ θ {\displaystyle {{d}(ml^{2}{\dot {\theta }}) \over dt}=-mgl\sin \theta } ここで、t：時間、m：振り子質量（≒操作者質量）、θ：振り子角度、g：重力加速度である。微小振動を仮定して sin ⁡ θ ≈ θ {\displaystyle \sin \theta \approx \theta } とすれば、以下のように変形できる。 θ ¨ + 2 l ˙ l θ ˙ + g l θ = 0 {\displaystyle {\ddot {\theta }}+2{\frac {\dot {l}}{l}}{\dot {\theta }}+{\frac {g}{l}}\theta =0} ブランコを漕ぐ動作とは、立ち漕ぎの場合は上半身を上下移動させ、座り漕ぎの場合は足を上下させる動作を行う。これらは、漕ぎ手の重心を上下させていることに等しい。これを、ブランコとそれに乗る漕ぎ手を合わせた一体の系として考えると、振り子のロープの長さが短くなったり、長くなったりすることに等しい。振り子腕長さ（ロープの長さ）を l とし、ブランコを漕ぐことによる振り子腕長さの変化を正弦波による周期的変化として以下のように表す。 l = l 0 + a sin ⁡ ω t {\displaystyle l=l_{0}+a\sin {\omega t}} ここで、l0：平均腕長さ、a：腕長さ変動振幅、ω：腕長さ変動の角周波数である。 これを上記の運動方程式に代入して、かつ、 a ≪ l 0 {\displaystyle a\ll l_{0}} として 1 / ( 1 + a / l 0 ) ≈ 1 − a / l 0 {\displaystyle 1/(1+a/l_{0})\approx 1-a/l_{0}} を利用すれば、以下の運動方程式を得ることができる。 θ ¨ + 2 ω ( a cos ⁡ ω t l 0 − a 2 2 l 0 2 sin ⁡ 2 ω t ) θ ˙ + g l 0 ( 1 − a l 0 sin ⁡ ω t ) θ = 0 {\displaystyle {\ddot {\theta }}+2\omega \left({\frac {a\cos {\omega t}}{l_{0}}}-{\frac {a^{2}}{2l_{0}^{2}}}\sin {2\omega t}\right){\dot {\theta }}+{\frac {g}{l_{0}}}\left(1-{\frac {a}{l_{0}}}\sin {\omega t}\right)\theta =0} さらに、 a ≪ l 0 {\displaystyle a\ll l_{0}} より ( a / l 0 ) 2 ≈ 0 {\displaystyle (a/l_{0})^{2}\approx 0} として簡単化すれば、 θ ¨ + 2 ω a cos ⁡ ω t l 0 θ ˙ + g l 0 ( 1 − a l 0 sin ⁡ ω t ) θ = 0 {\displaystyle {\ddot {\theta }}+2\omega {\frac {a\cos {\omega t}}{l_{0}}}{\dot {\theta }}+{\frac {g}{l_{0}}}\left(1-{\frac {a}{l_{0}}}\sin {\omega t}\right)\theta =0} よってブランコの運動方程式は、 θ ˙ {\displaystyle {\dot {\theta }}} 、 θ {\displaystyle \theta } の係数が変動する係数励振系で表すことができる。 この系の周期関数係数を無視した固有振動数は、 ω 0 = g / l 0 {\displaystyle \omega _{0}=g/l_{0}} である。フロケ理論を基に近似的に上式の不安定領域を求めると、第1次係数共振域では、 a / l 0 ≈ 0 {\displaystyle a/l_{0}\approx 0} のとき ω / ω 0 ≈ 2 {\displaystyle \omega /\omega _{0}\approx 2} で不安定振動発生し、そこから a / l 0 {\displaystyle a/l_{0}} が大きくなるに連れて不安定振動（発散）が発生する ω / ω 0 {\displaystyle \omega /\omega _{0}} の領域が広がるような分布となる。よって漕ぐ動作の周期に着目すると、ブランコを大きく揺らす漕ぎ方としては、ブランコの固有周期の2倍周期で重心の上下運動を行うことが理想的となる。 ^ a b 日本機械学会 編『機械工学辞典』（第2版）丸善、2007年1月20日、358頁。ISBN 978-4-88898-083-8。  ^ 振動工学 p.241 ^ 片桐道男, 山崎弘郎「パラメトリック発振増幅器の開発とβ線厚さ計への応用」『計測自動制御学会論文集』第1巻第1号、計測自動制御学会、1965年、75-83頁、doi:10.9746/sicetr1965.1.75、ISSN 0453-4654。  ^ a b c d 機械振動学 p.183 ^ a b 振動工学 pp.248-249 ^ 振動工学 pp.247-248 ^ パラメータ励振 p.6 ^ 吉田雅昭『Mathieu方程式に基づく平面磁路形パラメトリック変圧器の動作特性と発振安定性に関する研究』 八戸工業大学〈博士 （工学） 甲第49号〉、2011年、177頁。http://id.nii.ac.jp/1078/00003416/。  ^ a b c パラメータ励振 pp.10-11 ^ パラメータ励振 p.106 ^ パラメータ励振 p.110 ^ 機械振動学 p.185 前澤成一郎『振動工学』（第1版）森北出版、1973年11月20日。  小寺忠『パラメータ励振』（第1版）森北出版、2010年7月7日。ISBN 978-4-627-66741-9。  末岡淳男・金光陽一・近藤孝広『機械振動学』（初版）朝倉書店、2002年6月20日。ISBN 4-254-23706-5。  強制振動 自励振動