振り子型動吸振器の模擬アニメーション 台湾の台北101 に設置されている振り子型動吸振器配置図 動吸振器 （どうきゅうしんき、dynamic vibration absorber、DVA）またはダイナミックダンパ (dynamic damper)とは、振動 する対象物に、補助的な質量体をばね などを介して付加することにより、対象物の固有振動数 周辺での共振 現象を抑制する装置のことである[1] [2]
振動抑制技術のうち最も基本的なものの一つであり、機械 をはじめ建築 、土木 などの分野でも用いられる。同調質量ダンパ （チューンドマスダンパ、Tuned Mass Damper、TMD）や質量ダンパ （マスダンパ、Mass Damper）などとも呼ばれる。
1883年 にP・ワッツ(P. Watts)に考案され[3] [4]
概要 機械や建造物に振動が発生するとき、多くの場合で振動が害をなすので振動を抑制したい。特に、振動を発生させる力の振動数と対象物の固有振動数 が近い場合、共振 が発生して大きな振動が発生する（強制振動 なども参照）。これを避けるためには対象物の固有振動数を変更するなどの適切な振動特性を対象物に与える必要がある。しかし実際の機械や建造物の設計では、種々の制約条件により対象物自体の特性を都合良く変更することができないことも多い。このようなときに、対象物に補助的な質量体を取り付け、この質量体に対象物の振動を吸収させて代わりに振動させることで対象物の振動抑制を図るのが動吸振器である。建造物における設計思想の「耐震 」・「制振 」・「免震 」の内、「制振」に分類される装置に該当する[5]
基礎理論 非減衰動吸振器 主系・従系ばね のみの最も単純なモデル、主系が加振力f (t )を受ける場合。図中において、k はばね定数 、m は質量、x はばねの自然長（つり合いの位置）からの変位、t は時刻をそれぞれ表す変数、m、aは主系、従系を表す添え字である。 最も単純な減衰の無い2自由度 系の線形ばね質量系について考える。ばね k m で支えられた質量m m からなる主系（制振対象）に、ばねk a と質量m a からなる従系（動吸振器）が取り付けられたモデルが、動吸振器の最も単純なモデルとなる。このモデルでは、2つの質量は質点（ある一点に質量が集中している）とし、ばねの重さは考えない。質点m m 、m a の変位（ばねk m 、k a がともにつり合っている位置からの移動量）をそれぞれx m 、x a と表すこととし、時刻をt とおくと、主系が力振幅f 0 、角振動数Ω の調和振動 形の加振力を受けるとき、運動方程式は以下のようになる[6]
m m x ¨ m + k m x m + k a ( x m − x a ) = f 0 cos ( Ω t ) {\displaystyle m_{m}{\ddot {x}}_{m}+k_{m}x_{m}+k_{a}(x_{m}-x_{a})=f_{0}\cos(\Omega t)} m a x ¨ a + k a ( x a − x m ) = 0 {\displaystyle m_{a}{\ddot {x}}_{a}+k_{a}(x_{a}-x_{m})=0} この系全体としての固有振動数ω 1 、ω 2 は次のように求まる。
( ω 1 ω 2 ) = ω m 2 + ( μ + 1 ) ω a 2 ∓ [ ω m 2 + ( μ − 1 ) ω a 2 ] 2 + 4 μ ω a 4 2 {\displaystyle {\dbinom {\omega _{1}}{\omega _{2}}}={\sqrt {\frac {\omega _{m}^{2}+(\mu +1)\omega _{a}^{2}\mp {\sqrt {\left[\omega _{m}^{2}+(\mu -1)\omega _{a}^{2}\right]^{2}+4\mu \omega _{a}^{4}}}}{2}}}} ここで、
ω a = k a m a , ω m = k m m m , μ = m a m m {\displaystyle \omega _{a}={\sqrt {\frac {k_{a}}{m_{a}}}},\qquad \omega _{m}={\sqrt {\frac {k_{m}}{m_{m}}}},\qquad \mu ={\frac {m_{a}}{m_{m}}}} である。このとき各質点の振幅倍率は以下のように得られる。
X m x s t = ω m 2 ( ω a 2 − Ω 2 ) ( ω 1 2 − Ω 2 ) ( ω 2 2 − Ω 2 ) {\displaystyle {\frac {X_{m}}{x_{st}}}={\frac {\omega _{m}^{2}(\omega _{a}^{2}-\Omega ^{2})}{(\omega _{1}^{2}-\Omega ^{2})(\omega _{2}^{2}-\Omega ^{2})}}} X a x s t = ω m 2 ω a 2 ( ω 1 2 − Ω 2 ) ( ω 2 2 − Ω 2 ) {\displaystyle {\frac {X_{a}}{x_{st}}}={\frac {\omega _{m}^{2}\omega _{a}^{2}}{(\omega _{1}^{2}-\Omega ^{2})(\omega _{2}^{2}-\Omega ^{2})}}} 非減衰主系非減衰動吸振器の共振曲線、質量比を変化させた場合 ここで、X m /x st の式に注目すると、ω a = Ω のとき、X m /x st = 0となる。すなわち加振力の振動数Ω が既知のとき、動吸振器の単体固有角振動数ω a をΩ と一致させるように設計することで、主系の振動を完全に消失させることができる[6] 同調 (tuning)とよぶ[6] 反共振（英語版） とよぶ。
また、ω a = Ω のときの従系質量体の振動変位の解[6]
x a = f 0 k a cos ( Ω t + π ) {\displaystyle x_{a}={\frac {f_{0}}{k_{a}}}\cos(\Omega t+\pi )} となり、振幅f 0 / k a で、位相 は加振力と180°ずれて振動する。さらに、このときの従系質量体がばねを通じて主系質量体へ及ぼす力F a は、x m = 0なので、
F a = − k a ( x m − x a ) = f 0 cos ( Ω t + π ) = − f 0 cos ( Ω t ) {\displaystyle F_{a}=-k_{a}(x_{m}-x_{a})=f_{0}\cos(\Omega t+\pi )=-f_{0}\cos(\Omega t)} となり、加振力を完全に打ち消すような力が、従系から主系へ加わっていることがわかる[7]
以上のように、理論上はω a をΩ と一致させるように設計すれば主系の振動を0にできるが、実際にはΩ が一定値に限定できる場合は少ない。共振曲線を見ると、Ω /ω a =1の反共振点のすぐそばにω 1 、ω 2 による共振点が存在する。すなわち、Ω が反共振点から変動すると振幅はすぐに大きくなる傾向がある。一方、主系・従系質量比μ = m a /m m に注目すると、μ が大きいほど（=従系質量が主系質量に近いほど）、反共振点から離れても振幅倍率の立ち上がりが緩やかである。振動抑制の観点からは、このように主系・従系質量比を大きく取る方が都合が良いが、実際の設計ではそのような大きな動吸振器を付けることは通常は制約がある[8]
減衰付動吸振器 主系ばねのみ、従系ばね・減衰有りのモデル、主系が加振力f (t )を受ける場合。図中において、k はばね定数、c は減衰定数、m は質量、x はばねの自然長（つり合いの位置）からの変位、t は時刻をそれぞれ表す変数、m、aは主系、従系を表す添え字である。 減衰の無いモデルの動吸振器では、加振力の振動数と動吸振器の単体固有角振動数が一致または狭い範囲で近くないと効果を発揮できない。減衰のある動吸振器では、比較的広い範囲に加振力の振動数が変動する場合でも、主系の振動を吸収することが可能となる。1909年にフラームにより考案された動吸振器は減衰が無い単純なものであった[4] [4]
ばねk m で支えられた質量m m からなる主系（制振対象）に、ばねk a 、減衰器c a と質量m a からなる従系（動吸振器）が取り付けられたモデルを考える。減衰器がないモデルと同様に、各質点の変位をx 、時刻をt とおくと、この運動方程式は以下のようになる[9]
m m x ¨ m + k m x m + c a ( x ˙ m − x ˙ a ) + k a ( x m − x a ) = f 0 e i Ω t {\displaystyle m_{m}{\ddot {x}}_{m}+k_{m}x_{m}+c_{a}({\dot {x}}_{m}-{\dot {x}}_{a})+k_{a}(x_{m}-x_{a})=f_{0}e^{i\Omega t}} m a x ¨ a + c a ( x ˙ a − x ˙ m ) + k a ( x a − x m ) = 0 {\displaystyle m_{a}{\ddot {x}}_{a}+c_{a}({\dot {x}}_{a}-{\dot {x}}_{m})+k_{a}(x_{a}-x_{m})=0} この運動方程式より、主系の変位倍率は次のように求まる[9]
| X m x s t | = ( α 2 − β 2 ) 2 + ( 2 ζ a α β ) 2 [ ( α 2 − β 2 ) ( 1 − β 2 ) − μ α 2 β 2 ] 2 + ( 2 ζ a α β ) 2 ( 1 − β 2 − μ β 2 ) 2 {\displaystyle \left|{\frac {X_{m}}{x_{st}}}\right\vert ={\sqrt {\frac {(\alpha ^{2}-\beta ^{2})^{2}+(2\zeta _{a}\alpha \beta )^{2}}{\left[(\alpha ^{2}-\beta ^{2})(1-\beta ^{2})-\mu \alpha ^{2}\beta ^{2}\right]^{2}+(2\zeta _{a}\alpha \beta )^{2}(1-\beta ^{2}-\mu \beta ^{2})^{2}}}}} ここで、
ω a = k a m a , ω m = k m m m , μ = m a m m , α = ω a ω m , β = Ω ω m , c c a = 2 m a k a , ζ a = c a c c a , x s t = f 0 k m {\displaystyle \omega _{a}={\sqrt {\frac {k_{a}}{m_{a}}}},\ \omega _{m}={\sqrt {\frac {k_{m}}{m_{m}}}},\ \mu ={\frac {m_{a}}{m_{m}}},\ \alpha ={\frac {\omega _{a}}{\omega _{m}}},\ \beta ={\frac {\Omega }{\omega _{m}}},\ c_{ca}=2{\sqrt {m_{a}k_{a}}},\ \zeta _{a}={\frac {c_{a}}{c_{ca}}},\ x_{st}={\frac {f_{0}}{k_{m}}}} である。
減衰比ζ a を変化させていくと、ζ a → 0のときは、上記の非減衰モデルに一致し、ζ a → ∞のときは、主系と従系は一体にふるまい、質量m = m a + m m 、ばね定数 k = k m の1自由度系のモデルに一致する[10] ζ a → 0でも、ζ a → ∞でも、共振点で振幅が無限大に発散することになる[11] [11]
定点理論 非減衰主系減衰付動吸振器の共振曲線 ブロックらの定点理論による最適結果の例 今、従系（動吸振器）の減衰特性を変化させていくことにより、主系の変位倍率がそれに連れてどのように変化するかに注目する。主系と従系の質量比μ 、主系と従系の単体固有角振動数比α を固定し、従系の減衰比ζ a を変化させて変位倍率X m /x st の変化を見ると、ある2つの主系単体固有角振動数と加振力振動数の比βの値で、ζa に無関係にX m /x st の値が定まる2つの点（P、Q）がある。これらの点を定点(fixed point)と呼ぶ[9] β の代わりにΩ /ωa の値で変化を見たときも同様である[10]
減衰の無い動吸振器では、加振力の振動数と動吸振器の単体固有角振動数が一致または狭い範囲で近くないと効果を発揮できないので、減衰を付与することで幅広い範囲で振動（振幅倍率）を抑えるようにしたい。上記の定点が存在する性質を利用して動吸振器特性の最適化を図るのが、動吸振器の定点理論 である。定点理論は、1932年のエーリッヒ・ハンカム (Erich Hahnkamm)の研究に始まり、1946年にJ・E・ブロック(J. E. Brock)によってほぼ完成された[4] [10]
共振曲線で、2つの定点が同じ値を取る。 共振曲線で、2つの定点が極大値を取る。 後者の操作は定点を共振点と一致させることと同義でもある。具体的には、主系と従系の質量比μ より以下のような最適値が求まる[12]
α o p t = 1 1 + μ {\displaystyle \alpha _{opt}={\frac {1}{1+\mu }}} ζ a o p t = 3 μ 8 ( 1 + μ ) {\displaystyle \zeta _{a\ opt}={\sqrt {\frac {3\mu }{8(1+\mu )}}}} 以上のような、最適な主系・従系質量比α opt を求めることを最適同調 、最適な従系減衰比ζ a opt を求めることを最適減衰 とよぶ[13] [4]
上記の最適減衰の式は厳密解ではなく、平均に基づく近似値である[14] μ ≪ 1と見なせる限り、実用上は特に問題ない。誤差が生じる原因は、P点の最適減衰によると、右側の共振点がQ点からわずかにずれ、同様に、Q点の最適減衰では、左側の共振点がP点からずれるためで、定点理論はP点、Q点での最適減衰を平均化して結果としている。西原らの厳密解との比較によると、μ = 0.1のとき相対差0.023%、μ = 1のとき相対差0.5%、μ = 10のとき相対差2.3%である[15] [16]
z o p t = α o p t 2 ( 1 − r o p t ) [ r o p t − ( 1 + μ ) α o p t 2 ] 2 r o p t [ ( 1 + μ ) 2 α o p t 2 − r o p t ] {\displaystyle z_{opt}={\sqrt {\frac {\alpha _{opt}^{2}(1-r_{opt})\left[r_{opt}-(1+\mu )\alpha _{opt}^{2}\right]}{2r_{opt}\left[(1+\mu )^{2}\alpha _{opt}^{2}-r_{opt}\right]}}}} α o p t = r o p t ( 1 + 4 + 3 μ ) 3 ( 1 + μ ) 2 {\displaystyle \alpha _{opt}={\sqrt {\frac {r_{opt}(1+{\sqrt {4+3\mu }})}{3(1+\mu )^{2}}}}} ここで、
z = c a 2 μ m m k m = α ζ a {\displaystyle z={\frac {c_{a}}{2\mu {\sqrt {m_{m}k_{m}}}}}=\alpha \zeta _{a}} r o p t = 8 [ ( 4 + 3 μ ) 3 / 2 − μ ] 64 + 80 μ + 27 μ 2 {\displaystyle r_{opt}={\frac {8\left[(4+3\mu )^{3/2}-\mu \right]}{64+80\mu +27\mu ^{2}}}} 最小分散規範 最小分散規範（赤）と定点理論（青）の比較例 非常に不規則な励振を受けるなどの場合は、特定の励振周波数近辺よりも、すべての周波数域で振動が最少となるように設計した方がよい[12] 最小分散規範 と呼ばれる最適化法がある。最小分散規範は、1963年にステファン・H・クランドル(Stephen H. Crandall)とウィリアム・D・マーク(William D. Mark)により発表された [17]
最小分散規範では、伝達される振動エネルギーに注目して、これが最少となるように設計する。すなわち、共振曲線を積分して得られる面積二乗値が最少となるようにする[12] ホワイトノイズ ランダム振動を行う場合は、主系と従系の質量比μ より以下のような最適値が求まる[18]
α o p t = 1 1 + μ 2 + μ 2 {\displaystyle \alpha _{opt}={\frac {1}{1+\mu }}{\sqrt {\frac {2+\mu }{2}}}} ζ a o p t = μ ( 4 + 3 μ ) 8 ( 1 + μ ) ( 2 + μ ) {\displaystyle \zeta _{a\ opt}={\sqrt {\frac {\mu (4+3\mu )}{8(1+\mu )(2+\mu )}}}} 上記の通り、定点理論と異なり2つの共振点の高さは一致せず、常に曲線上左側（周波数が低い側）の共振点が、曲線上右側（周波数が高い側）の共振点よりも大きくなる特徴がある[19]
減衰付主系・減衰付動吸振器 主系、従系ともにばね・減衰有りのモデル、主系が加振力を受ける場合
主系、従系ともにばね・減衰有りのモデル、基礎が振動変位する場合
より一般的な、主系にも減衰がある場合を考える。ばねk m 、減衰器c m で基礎に支えられた質量m m からなる主系（制振対象）に、ばねk a 、減衰器c a と質量m a からなる従系（動吸振器）が取り付けられたモデルの運動方程式は、主系に対して力励振f (t )が加わる場合と基礎に対して変位励振x 0 (t )が発生する場合、それぞれで以下のようになる。
主系に対して力励振f (t )が加わる場合：
m m x ¨ m + c m x ˙ m + k m x m + c a ( x ˙ m − x ˙ a ) + k a ( x m − x a ) = f ( t ) {\displaystyle m_{m}{\ddot {x}}_{m}+c_{m}{\dot {x}}_{m}+k_{m}x_{m}+c_{a}({\dot {x}}_{m}-{\dot {x}}_{a})+k_{a}(x_{m}-x_{a})=f(t)} m a x ¨ a + c a ( x ˙ a − x ˙ m ) + k a ( x a − x m ) = 0 {\displaystyle m_{a}{\ddot {x}}_{a}+c_{a}({\dot {x}}_{a}-{\dot {x}}_{m})+k_{a}(x_{a}-x_{m})=0} 基礎に対して変位励振x 0 (t )が発生する場合：
m m x ¨ m + c m ( x ˙ m − x ˙ 0 ( t ) ) + k m ( x m − x 0 ( t ) ) + c a ( x ˙ m − x ˙ a ) + k a ( x m − x a ) = 0 {\displaystyle m_{m}{\ddot {x}}_{m}+c_{m}({\dot {x}}_{m}-{\dot {x}}_{0}(t))+k_{m}(x_{m}-x_{0}(t))+c_{a}({\dot {x}}_{m}-{\dot {x}}_{a})+k_{a}(x_{m}-x_{a})=0} m a x ¨ a + c a ( x ˙ a − x ˙ m ) + k a ( x a − x m ) = 0 {\displaystyle m_{a}{\ddot {x}}_{a}+c_{a}({\dot {x}}_{a}-{\dot {x}}_{m})+k_{a}(x_{a}-x_{m})=0} 主系に対して力励振f = f 0 sin(Ωt )が加わる場合は、主系の変位倍率は次のように求まる[20]
| X m x s t | = ( α 2 − β 2 ) 2 + ( 2 ζ a α β ) 2 A 2 + 4 B 2 {\displaystyle \left|{\frac {X_{m}}{x_{st}}}\right\vert ={\sqrt {\frac {(\alpha ^{2}-\beta ^{2})^{2}+(2\zeta _{a}\alpha \beta )^{2}}{A^{2}+4B^{2}}}}} A = β 4 − { ( μ + 1 ) α 2 + 4 ζ a ζ m α + 1 } β 2 + α 2 {\displaystyle A=\beta ^{4}-\left\{(\mu +1)\alpha ^{2}+4\zeta _{a}\zeta _{m}\alpha +1\right\}\beta ^{2}+\alpha ^{2}} B = − { ζ m + ( μ + 1 ) ζ a α } β 3 + ( ζ m α + ζ a ) α β {\displaystyle B=-\left\{\zeta _{m}+(\mu +1)\zeta _{a}\alpha \right\}\beta ^{3}+(\zeta _{m}\alpha +\zeta _{a})\alpha \beta } ここで、
ω a = k a m a , ω m = k m m m , μ = m a m m , α = ω a ω m , β = Ω ω m {\displaystyle \omega _{a}={\sqrt {\frac {k_{a}}{m_{a}}}},\ \omega _{m}={\sqrt {\frac {k_{m}}{m_{m}}}},\ \mu ={\frac {m_{a}}{m_{m}}},\ \alpha ={\frac {\omega _{a}}{\omega _{m}}},\ \beta ={\frac {\Omega }{\omega _{m}}}} c c a = 2 m a k a , c c m = 2 m m k m , ζ a = c a c c a , ζ m = c m c c m , x s t = f 0 k m {\displaystyle c_{ca}=2{\sqrt {m_{a}k_{a}}},\ c_{cm}=2{\sqrt {m_{m}k_{m}}},\ \zeta _{a}={\frac {c_{a}}{c_{ca}}},\ \zeta _{m}={\frac {c_{m}}{c_{cm}}},\ x_{st}={\frac {f_{0}}{k_{m}}}} である。ζ m → 0のとき、上記の主系に減衰無しの場合の変位倍率と一致する。
一般に、主系に減衰要素が存在する場合は動吸振器の最適パラメータ（α opt 、ζ a opt ）の厳密解を得ることはできない[18] [21] 数値解析 により最適値を得る必要があり、多くの研究が行われてきている[3]
近似式 数値解析結果をもとにした、浅見らによる最適パラメータを求める近似式を以下に示す[22] [23] ζ m < 0.1程度の範囲まで実用的には十分な精度がある[24] [25] [21]
ここで、
ω a = k a m a , ω m = k m m m , μ = m a m m , α = ω a ω m {\displaystyle \omega _{a}={\sqrt {\frac {k_{a}}{m_{a}}}},\ \omega _{m}={\sqrt {\frac {k_{m}}{m_{m}}}},\ \mu ={\frac {m_{a}}{m_{m}}},\ \alpha ={\frac {\omega _{a}}{\omega _{m}}}} c c a = 2 m a k a , ζ a = c a c c a , c c m = 2 m m k m , ζ m = c m c c m {\displaystyle \ c_{ca}=2{\sqrt {m_{a}k_{a}}},\ \zeta _{a}={\frac {c_{a}}{c_{ca}}},\ c_{cm}=2{\sqrt {m_{m}k_{m}}},\ \zeta _{m}={\frac {c_{m}}{c_{cm}}}} とすれば、
最小分散規範による最適値： [22]
α o p t = 1 1 + μ 1 + μ 2 − ζ m ( 4 + μ ) μ 8 ( 1 + μ ) 3 ( 2 + μ ) ( 4 + 3 μ ) + ζ m 2 μ ( 192 + 304 μ + 132 μ 2 + 13 μ 3 ) 8 ( 1 + μ ) 2 ( 4 + 3 μ ) 2 2 ( 2 + μ ) 3 − ζ m 3 b 1 16 μ 3 2 ( 1 + μ ) 5 ( 2 + μ ) 5 ( 4 + 3 μ ) 7 {\displaystyle {\begin{alignedat}{3}\alpha _{opt}=&{\frac {1}{1+\mu }}{\sqrt {1+{\frac {\mu }{2}}}}-\zeta _{m}(4+\mu ){\sqrt {\frac {\mu }{8(1+\mu )^{3}(2+\mu )(4+3\mu )}}}\\&+\zeta _{m}^{2}{\frac {\mu (192+304\mu +132\mu ^{2}+13\mu ^{3})}{8(1+\mu )^{2}(4+3\mu )^{2}{\sqrt {2(2+\mu )^{3}}}}}\\&-\zeta _{m}^{3}{\frac {b_{1}}{16}}{\sqrt {\frac {\mu ^{3}}{2(1+\mu )^{5}(2+\mu )^{5}(4+3\mu )^{7}}}}\\\end{alignedat}}} ζ a o p t = μ ( 4 + 3 μ ) 8 ( 1 + μ ) ( 2 + μ ) − ζ m μ 3 4 ( 1 + μ ) ( 4 + 3 μ ) 2 ( 2 + μ ) 3 + ζ m 2 − 64 − 80 μ + 15 μ 3 32 2 μ 5 ( 1 + μ ) 3 ( 2 + μ ) 5 ( 4 + 3 μ ) 5 + ζ m 3 μ 3 b 2 32 ( 1 + μ ) 2 ( 4 + 3 μ ) 4 2 ( 2 + μ ) 7 {\displaystyle {\begin{alignedat}{3}\zeta _{a\ opt}=&{\sqrt {\frac {\mu (4+3\mu )}{8(1+\mu )(2+\mu )}}}-\zeta _{m}{\frac {\mu ^{3}}{4(1+\mu )(4+3\mu ){\sqrt {2(2+\mu )^{3}}}}}\\&+\zeta _{m}^{2}{\frac {-64-80\mu +15\mu ^{3}}{32}}{\sqrt {\frac {2\mu ^{5}}{(1+\mu )^{3}(2+\mu )^{5}(4+3\mu )^{5}}}}\\&+\zeta _{m}^{3}{\frac {\mu ^{3}b_{2}}{32(1+\mu )^{2}(4+3\mu )^{4}{\sqrt {2(2+\mu )^{7}}}}}\\\end{alignedat}}} ただし、
b 1 = 4096 + 13056 μ + 15360 μ 2 + 8080 μ 3 + 1780 μ 4 + 101 μ 5 {\displaystyle b_{1}=4096+13056\mu +15360\mu ^{2}+8080\mu ^{3}+1780\mu ^{4}+101\mu ^{5}} b 2 = 2048 + 6912 μ + 8064 μ 2 + 3616 μ 3 + 288 μ 4 − 125 μ 5 {\displaystyle b_{2}=2048+6912\mu +8064\mu ^{2}+3616\mu ^{3}+288\mu ^{4}-125\mu ^{5}} α o p t = 1 1 + μ 1 + μ 2 − ζ m ( 4 + μ ) μ 8 ( 1 + μ ) 3 ( 2 + μ ) ( 4 + 3 μ ) + ζ m 2 μ ( 704 + 1328 μ + 804 μ 2 + 157 μ 3 ) 8 ( 1 + μ ) 2 ( 4 + 3 μ ) 2 2 ( 2 + μ ) 3 + ζ m 3 b 1 16 μ 2 ( 1 + μ ) 5 ( 2 + μ ) 5 ( 4 + 3 μ ) 7 {\displaystyle {\begin{alignedat}{3}\alpha _{opt}=&{\frac {1}{1+\mu }}{\sqrt {1+{\frac {\mu }{2}}}}-\zeta _{m}(4+\mu ){\sqrt {\frac {\mu }{8(1+\mu )^{3}(2+\mu )(4+3\mu )}}}\\&+\zeta _{m}^{2}{\frac {\mu (704+1328\mu +804\mu ^{2}+157\mu ^{3})}{8(1+\mu )^{2}(4+3\mu )^{2}{\sqrt {2(2+\mu )^{3}}}}}\\&+\zeta _{m}^{3}{\frac {b_{1}}{16}}{\sqrt {\frac {\mu }{2(1+\mu )^{5}(2+\mu )^{5}(4+3\mu )^{7}}}}\\\end{alignedat}}} ζ a o p t = μ ( 4 + 3 μ ) 8 ( 1 + μ ) ( 2 + μ ) − ζ m μ 3 4 ( 1 + μ ) ( 4 + 3 μ ) 2 ( 2 + μ ) 3 + ζ m 2 4096 + 13760 μ + 18608 μ 2 + 12640 μ 3 + 4287 μ 4 + 576 μ 5 32 2 μ 3 ( 1 + μ ) 3 ( 2 + μ ) 5 ( 4 + 3 μ ) 5 + ζ m 3 μ 3 b 2 32 ( 1 + μ ) 2 ( 4 + 3 μ ) 4 2 ( 2 + μ ) 7 {\displaystyle {\begin{alignedat}{3}\zeta _{a\ opt}=&{\sqrt {\frac {\mu (4+3\mu )}{8(1+\mu )(2+\mu )}}}-\zeta _{m}{\frac {\mu ^{3}}{4(1+\mu )(4+3\mu ){\sqrt {2(2+\mu )^{3}}}}}\\&+\zeta _{m}^{2}{\frac {4096+13760\mu +18608\mu ^{2}+12640\mu ^{3}+4287\mu ^{4}+576\mu ^{5}}{32}}{\sqrt {\frac {2\mu ^{3}}{(1+\mu )^{3}(2+\mu )^{5}(4+3\mu )^{5}}}}\\&+\zeta _{m}^{3}{\frac {\mu ^{3}b_{2}}{32(1+\mu )^{2}(4+3\mu )^{4}{\sqrt {2(2+\mu )^{7}}}}}\\\end{alignedat}}} ただし、
b 1 = 65536 + 241664 μ + 369920 μ 2 + 305664 μ 3 + 148720 μ 4 + 43500 μ 5 + 7339 μ 6 + 576 μ 7 {\displaystyle b_{1}=65536+241664\mu +369920\mu ^{2}+305664\mu ^{3}+148720\mu ^{4}+43500\mu ^{5}+7339\mu ^{6}+576\mu ^{7}} b 2 = 524288 + 2818048 μ + 6621184 μ 2 + 8864512 μ 3 + 7377280 μ 4 + 3896224 μ 5 + 1271168 μ 6 + 233491 μ 7 + 18432 μ 8 {\displaystyle b_{2}=524288+2818048\mu +6621184\mu ^{2}+8864512\mu ^{3}+7377280\mu ^{4}+3896224\mu ^{5}+1271168\mu ^{6}+233491\mu ^{7}+18432\mu ^{8}} 定点理論による最適値： [23]
α o p t = 1 1 + μ − ζ m 1 1 + μ 1 2 ( 1 + μ ) ( 3 + 4 μ − 2 ( 2 + μ ) ( 9 + 4 μ ) 2 + μ ) + ζ m 2 b 0 − 4 ( 5 + 2 μ ) 2 ( 2 + μ ) ( 9 + 4 μ ) 4 ( 1 + μ ) 2 ( 2 + μ ) ( 9 + 4 μ ) {\displaystyle {\begin{alignedat}{2}\alpha _{opt}=&{\frac {1}{1+\mu }}-\zeta _{m}{\frac {1}{1+\mu }}{\sqrt {{\frac {1}{2(1+\mu )}}\left(3+4\mu -{\frac {\sqrt {2(2+\mu )(9+4\mu )}}{2+\mu }}\right)}}\\&+\zeta _{m}^{2}{\frac {b_{0}-4(5+2\mu ){\sqrt {2(2+\mu )(9+4\mu )}}}{4(1+\mu )^{2}(2+\mu )(9+4\mu )}}\\\end{alignedat}}} ζ a o p t = 3 μ 8 ( 1 + μ ) + ζ m 60 + 63 μ + 16 μ 2 − 2 ( 3 + 2 μ ) 2 ( 2 + μ ) ( 9 + 4 μ ) 8 ( 1 + μ ) ( 2 + μ ) ( 9 + 4 μ ) + ζ m 2 b 1 ( A + B ) 2 + μ + b 2 ( A − B ) μ 32 ( 1 + μ ) ( 9 + 4 μ ) μ ( 2 + μ ) {\displaystyle {\begin{alignedat}{2}\zeta _{a\ opt}=&{\sqrt {\frac {3\mu }{8(1+\mu )}}}+\zeta _{m}{\frac {60+63\mu +16\mu ^{2}-2(3+2\mu ){\sqrt {2(2+\mu )(9+4\mu )}}}{8(1+\mu )(2+\mu )(9+4\mu )}}\\&+\zeta _{m}^{2}{\frac {b_{1}(A+B){\sqrt {2+\mu }}+b_{2}(A-B){\sqrt {\mu }}}{32(1+\mu )(9+4\mu ){\sqrt {\mu (2+\mu )}}}}\\\end{alignedat}}} A = 3 ( 2 + μ ) − μ ( 2 + μ ) {\displaystyle A={\sqrt {3(2+\mu )-{\sqrt {\mu (2+\mu )}}}}} B = 3 ( 2 + μ ) + μ ( 2 + μ ) {\displaystyle B={\sqrt {3(2+\mu )+{\sqrt {\mu (2+\mu )}}}}} ただし、
b 0 = 52 + 41 μ + 8 μ 2 {\displaystyle b_{0}=52+41\mu +8\mu ^{2}} b 1 = − 1296 + 2124 μ + 6509 μ 2 + 5024 μ 3 + 1616 μ 4 + 192 μ 5 {\displaystyle b_{1}=-1296+2124\mu +6509\mu ^{2}+5024\mu ^{3}+1616\mu ^{4}+192\mu ^{5}} b 2 = 48168 + 112887 μ + 105907 μ 2 + 49664 μ 3 + 11632 μ 4 + 1088 μ 5 {\displaystyle b_{2}=48168+112887\mu +105907\mu ^{2}+49664\mu ^{3}+11632\mu ^{4}+1088\mu ^{5}} α opt とζ a opt の式は、上記の主系力加振系と同形式である。ただし係数は以下のように変わる。
b 0 = 52 + 113 μ + 76 μ 2 + 16 μ 3 {\displaystyle b_{0}=52+113\mu +76\mu ^{2}+16\mu ^{3}} b 1 = − 1296 + 2124 μ + 7157 μ 2 + 5924 μ 3 + 2032 μ 4 + 256 μ 5 {\displaystyle b_{1}=-1296+2124\mu +7157\mu ^{2}+5924\mu ^{3}+2032\mu ^{4}+256\mu ^{5}} b 2 = 48168 + 105111 μ + 91867 μ 2 + 40172 μ 3 + 8784 μ 4 + 768 μ 5 {\displaystyle b_{2}=48168+105111\mu +91867\mu ^{2}+40172\mu ^{3}+8784\mu ^{4}+768\mu ^{5}} 動吸振器の種類 イングランド、ストックトン＝オン＝ティーズ のインフィニティ橋(en:Infinity Bridge)に設置されている同調質量ダンパ ミレニアム・ブリッジ (ロンドン)に設置されている同調質量ダンパ 台湾の台北101に設置されている振り子型同調質量ダンパ 以下に、派生形も含めた動吸振器の種類を示す。ここでは、補助質量体により対象物の振動を抑制・吸収する装置を動吸振器の定義としているが、ばねのような復元力要素を持つもののみを動吸振器に分ける場合もある[26]
同調質量ダンパ 固体質量体を用いた最も基本的な動吸振器で、上記で基本原理を説明しているものと同じである。動吸振器と同義の場合もあるが、同調質量ダンパ 、チューンドマスダンパ (en:Tuned mass damper、TMD) 、質量ダンパ 、マスダンパ (Mass Damper) などと呼ばれる。特に土木 ・建築 分野ではこちらの名称で呼ばれる[27] 振り子 や倒立振子機構にすることで復元力を得るもの、従系質量体を1つと限らずに複数備えるもの、振動方向が上下あるいは水平に対応するものなどがある。
建築分野では、超高層建築物 に対して風による揺れが問題になることから制振装置として採用される[28] アメリカ のシティグループ・センターとジョン・ハンコック・センター で、どちらも構造エンジニア のウィリアム・ルメジャー (en:William LeMessurier) の設計によるものである[29]
2007年現在において、世界で最重量の同調質量ダンパは、台湾の台北101 に設置されている730tonの振り子型同調質量ダンパである[30] 千葉ポートタワー のものが高層建築物に対する初の採用である[31]
自動車にも採用されており、2006年のフォーミュラ1 では、エアロダイナミクスの効果を高めるために車体の前端にマスダンパを装備して車体の安定性を向上させたフォーミュラ1カー が使用された[32] ルノー・R26 なども参照。
同調液体ダンパ（スロッシングダンパ） 同調液体ダンパ 、あるいはスロッシングダンパ (tuned sloshing damper)とは、液体のスロッシング 現象を利用して振動を抑制するもので、建物などの制振に利用される。スロッシングの固有振動数と建物の固有振動数を同調させることで効果を発揮する。スロッシングの固有振動数の固有振動数は、液体を入れる水槽の大きさ、水深などにより決定される。日本では、神奈川県横浜市の新横浜プリンスホテル 、北海道函館市の五稜郭タワー などで利用されている[33]
フードダンパ フードダンパ (houdaille damper, houde damper)とは、質量と減衰器からなる動吸振器で、通常のものからばね要素を無くした構造が特徴である[34] [35] [34] [35] [36] [36]
ランチェスターダンパ（摩擦ダンパ） ランチェスターダンパ (Lanchester damper)とは、フードダンパにおける減衰要素を摩擦要素で置き換えたもの[37] フレデリック・ランチェスター により往復内燃機関 のクランクシャフト 用として発明されたもので、ランチェスターダンパとはそのような構造・用途に限って意味する場合もある[38] 摩擦ダンパ などとも呼ぶ[39] [40]
インパクトダンパ（衝撃ダンパ） インパクトダンパ (impact damper)とは、制振対象に衝突体を入れた容器を取り付け、さらに容器内の振動方向に隙間を設けることで、対象物振動時に衝突 が発生して対象物の振動を抑制するものである[41] [26] 運動量 を衝突により補助質量体へ移動させて対象物の運動を抑制し、移った補助質量体の運動量は摩擦などで散逸させる仕組みを取る[42] [40]
動力付き自動車模型のミニ四駆 では、コースアップダウンセクション通過によるジャンプから着地時の車体の上下動安定を目的とした「マスダンパー」と呼ばれる改造部品がある[43]
アクティブ動吸振器 一般的に指す動吸振器とは受動制振（パッシブ）による振動抑制機構である。すなわち、質量、ばね、減衰器などを取り付けるだけで、制振のために特別にエネルギを外部から入力するようなことは行わない機構である[44] [44] アクチュエータ などを利用して直接的に振動抑制力を与えるのが能動制振型で[1] アクティブ動吸振器 、アクティブマスダンパ (active mass damper, AMD)などとよぶ。例えば、質量体を直接アクチュエータで押し引きして建物の振動抑制を図るものなどがある[45]
アクティブ式の他に、アクチュエータは使用しない代わりに応答状態に応じて動吸振器の減衰特性を変化させるセミアクティブ式や、アクチュエータとばね・ダンパを併用するハイブリッド式などもある[46] [1]
脚注 [脚注の使い方 ]
参考文献 日本機械学会 編『機械工学辞典』（第2版）丸善、2007年1月20日。ISBN 978-4-88898-083-8。 制振工学ハンドブック編集委員会 編『制振工学ハンドブック』コロナ社、2008年5月13日。ISBN 978-4-339-04585-7。 日本機械学会 編『振動のダンピング技術』（第1版）養賢堂、1998年9月1日。ISBN 4-8425-9816-6。 瀬戸一登『動吸振器とその応用』（初版）コロナ社、2010年8月20日。ISBN 978-4-339-04607-6。 末岡淳男・金光陽一・近藤孝広『機械振動学』（初版）朝倉書店、2002年6月20日。ISBN 4-254-23706-5。 前澤成一郎『振動工学』（第1版）森北出版、1973年11月20日。 Jerome J. Connor (2003). Introduction to Structural Motion Control (1st ed.). Prentice Hall. ISBN 0-13-009138-3 西原修、松久寛「動吸振器の最大振幅倍率最小化設計 : 代数的な厳密解の導出」『日本機械学会論文集 C編』第63巻第614号、1997年10月25日、3438-3445頁、ISSN 03875024、NAID 130004232443。 浅見敏彦、桃瀬一成, 細川[ヨシ]延「主系の減衰を考慮した動吸振器の設計式について : 最小分散規範に基づく設計法」『日本機械学会論文集C編』第59巻第566号、1993年10月25日、2962-2967頁、ISSN 03875024、NAID 110002380702。 浅見敏彦、細川[ヨシ]延「主系の減衰を考慮した動吸振器の設計式について : 第2報,定点理論に基づく設計法」『日本機械学会論文集C編』第61巻第583号、1995年3月25日、915-921頁、ISSN 03875024、NAID 130004231037。 Kefu Liu; Gianmarc Coppola (2010). “Optimal design of damped dynamic vibration absorber for damped primary systems” (pdf). Trans. Can. Soc. Mech. Eng. 34 (1): 120-135. http://www.tcsme.org/Papers/Vol34/Vol34No1Paper8.pdf . 関連項目 ウィキメディア・コモンズには、同調質量ダンパー
外部リンク 「動吸振器」 - 機械工学事典（日本機械学会 ） アメリカ合衆国特許第 989,958号 （フラームによる1909年の特許）屋上の水槽によるスロッシングダンパーの模式アニメーション（アイディールブレーン社スーパー・スロッシング・ダンパー製品紹介ページ） Tuned Mass Damper - YouTube（Moog CSA Engineering社による同調質量ダンパ説明動画）典拠管理データベース: 国立図書館
![{\displaystyle {\dbinom {\omega _{1}}{\omega _{2}}}={\sqrt {\frac {\omega _{m}^{2}+(\mu +1)\omega _{a}^{2}\mp {\sqrt {\left[\omega _{m}^{2}+(\mu -1)\omega _{a}^{2}\right]^{2}+4\mu \omega _{a}^{4}}}}{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/25ca4ce321901331d495ccd2ffffdf9b11bcf643)
![{\displaystyle \left|{\frac {X_{m}}{x_{st}}}\right\vert ={\sqrt {\frac {(\alpha ^{2}-\beta ^{2})^{2}+(2\zeta _{a}\alpha \beta )^{2}}{\left[(\alpha ^{2}-\beta ^{2})(1-\beta ^{2})-\mu \alpha ^{2}\beta ^{2}\right]^{2}+(2\zeta _{a}\alpha \beta )^{2}(1-\beta ^{2}-\mu \beta ^{2})^{2}}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f5145a785ab8610995aab63a2094912a5a68531e)
![{\displaystyle z_{opt}={\sqrt {\frac {\alpha _{opt}^{2}(1-r_{opt})\left[r_{opt}-(1+\mu )\alpha _{opt}^{2}\right]}{2r_{opt}\left[(1+\mu )^{2}\alpha _{opt}^{2}-r_{opt}\right]}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d52ff4d9e0b59bcf1e1515a3416ed2ab077c4366)
![{\displaystyle r_{opt}={\frac {8\left[(4+3\mu )^{3/2}-\mu \right]}{64+80\mu +27\mu ^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bfd441d031b59f33506369af6cb1f58d01c64bd0)
