デジタルツインを実現するCAEの真価

振動低減の戦略 〜ばね−マス系の振動【1】〜CAEと計測技術を使った振動・騒音対策(8)(3/5 ページ)

» 2023年06月08日 09時00分 公開

伝達関数をグラフで表現する

 もう少しお付き合いください。限界減衰co、減衰比ζ、無減衰固有角振動数ωno、無減衰固有振動数fnoを導入します。次式で定義されます。

式17 式17
式18 式18
式19 式19
式20 式20

 限界減衰coの意味するところは、図5の自由運動の解が2通りあって、図6に示すように比例粘性係数cがcoより小さいと、つまり減衰力が小さいと自由運動の解は振動するのですが(図6左図)、cがcoより大きいと、つまり減衰力が大きいと自由運動の解は一度も往復運動せずにゼロに近づきます(図6右図)。co図6左図になるか図6右図になるかの境目の比例粘性抵抗です。

限界減衰の意味 図6 限界減衰の意味[クリックで拡大]

 式17、18、19式12に代入すると次式になります。

式21 式21
式22 式22

 共振点、つまり固有振動数は式20で計算されることが多いのですが、減衰があると少しずれます。式21の値が最大となる角振動数が共振点となり、式21のルートの中を微分することで共振点の角振動数が求まります。次式となります。

式23 式23
式24 式24

 伝達関数の最大値は、式23式20に代入することで求まり、次式となります。分母にばね定数kと減衰比ζがあることに注目してください。

式25 式25

 数値を代入した例を図7に示します。位相角ですが、アークタンジェントの戻り値は+90〜−90[deg]なのですが、外力に対する応答が外力よりも早い時間に起こることはあり得ず、応答は常に遅れるので、0〜+90[deg]のときは180[deg]だけ引き算しています。文献でよく見掛けるグラフになりました。

ばね−マス系の伝達関数 図7 ばね−マス系の伝達関数[クリックで拡大]

伝達関数の裾野について

 周波数分析について説明した際、「ピーク値近傍の裾野の値は単にsinとcosを掛算した総和であって、信号には含まれない」といった旨のお話をしました。しかし、ばね−マス系の伝達関数の裾野には意味があります。図8に示すように、このばね−マス系を120[rad/s]で加振したときの振幅の伝達関数は0.0002[m/N]であって、振幅1[N]で加振したら変位の振幅は0.0002[m]となります。

伝達関数の裾野 図8 伝達関数の裾野[クリックで拡大]

Copyright © ITmedia, Inc. All Rights Reserved.