●知識グラフ全体 ,●このページを中心とした周辺知識グラフ
応用分野: 強制振動:微分方程式の解法1 , 強制振動:微分方程式の解法3 , 強制振動:微分方程式の解法4 , 強制振動:共振 , 強制振動:微分方程式の解法2

強制振動:運動方程式 (equation of motion)

x 軸上を単振動する質量 m の質点には,その位置 x に比例した復元力 −cx  ( c :正定数)が作用し,運動方程式は

m d2x dt2 =−cx     ······ 

と表される.この質点に時間的に変動する外力 Fex(t) が作用する場合を考えると,その運動方程式は

m d2x dt2 =−cx +Fex(t)     ······ 

となる.これが強制振動の従う運動方程式であり,一般に,強制振動において外力 Fex(t) は周期的(振動的)な力を考える.

式において, ω0 = c/m を導入して整理すると,定数係数の2階非同次線形微分方程式

d2x dt2 + ω02x = Fex(t) m     ······ 

が得られる. ω0 は単振動の固有角振動数である.この微分方程式を解くと,一般解が

x(t)=A cos(ω0t+α) +xp(t)     ( A , α :任意定数)     ······ 

のように,単振動(右辺第1項目)と外力 Fex(t) により強制的に引き起こされる変動 xp(t) との和(重ね合わせ)として求まる.ここで, xp(t) は

xp(t)= 1mω0 { sinω0t ∫ Fex(t) cosω0t dt −cosω0t ∫ Fex(t) sinω0t dt }     ······ 

である( xp(t) は任意定数を含まないので,上式の積分では積分定数を考慮する必要はない).

ここで,外力が

Fex(t)= F0sinωt     ( F0 , ω :定数)     ······ 

のように角振動数 ω で周期的に変動する場合を考えると, 強制振動項は, ω≠ω0 であれば

xp(t)= F0 m ( ω02−ω2 ) sinωt     ······ 

となる.上式は, ω<ω0 のとき外力と同位相となり, ω>ω0 のとき逆位相となる(位相が π だけずれる).一方, ω=ω0 であれば,式は発散してしまうため,この場合は

xp(t)=− F0 2mω0 tcosω0t     ······ 

で表される.この式から分かるように,時間の経過とともに振動の振幅が t に比例して増大する(共振).

外力が式で与えられる場合について,外力の角振動数 ω を変化させたときの解のグラフを以下に示す.ここで,式の初期条件として x(0)=0 , v(0)= dxdt(0)=0 を考え,それを満たす解

x(t)= F0 m ( ω02−ω2 ) { − (ω/ω0) sinω0t +sinωt }     ( ω≠ω0 )     ······ 

について, F0 m | ω02−ω2 | を一定に保ってグラフを描画した.


図1  ω<ω0 のときの解(赤点線:単振動,緑点線:外力による強制変動)


図2  ω>ω0 のときの解(赤点線:単振動,緑点線:外力による強制変動)

ω が ω0 近いと,うなりのような現象が生じることがわかる.

ω=ω0 については,式を用いて,初期条件を満たす解は

x(t)= F0 2mω02 { sinω0t −ω0t cosω0t }     ······ 

となり,そのグラフは下図となる(中段のグラフ).上段と下段のグラフは ω と ω0 が近いときの状況を式で描いたグラフである.


図3  ω=ω0 のときの解(中段)

ω=ω0 のときは,時間の経過とともに振動の振幅が t に比例して無限に増大する.


ホーム>>カテゴリー分類>>力学>>質点の力学>>強制振動>>運動方程式

最終更新日:2026年1月24日

[ページトップ]

金沢工業大学

利用規約,誤植や誤りなどがございましたら、こちらからご連絡ください。