CYBERNET

CAEを学ぶ

構造解析

はじめての振動解析(3)減衰の基礎

はじめての振動解析(3)減衰の基礎 の概要

はじめての振動解析(1)(2)(3)  完全保存版 PDF を "無料" 公開中!

はじめての振動解析(1)(2)(3) 完全保存版 PDF サンプル

これさえ読めば振動解析の基本がわかる!

全3回のシリーズとなっております「はじめての振動解析」について、3回分の解析講座を1冊でお読みいただける完全保存版 PDF を公開いたしました。どなたでも簡単なフォーム入力のみでダウンロードいただけますので、是非この機会にご利用ください。

ページ数:30ページ
ファイル形式: PDF

 

はじめての振動解析(1)(2)(3)  完全保存版 PDF の目次

(1) 振動の基礎とモーダル解析
 1. はじめに
 2. 振動のメカニズム
 3. 振動解析の基礎
 4. モーダル解析(固有値解析)
 5. モーダル解析でできること
 6. Ansys によるモーダル解析
 7. Ansysモーダル解析事例
 8. (補足) 実固有値解析と複素固有値解析
 9. おわりに
(2) 周波数応答解析モード重ね合わせ法の基礎
 1. はじめに
 2. 周波数応答解析
 3. 周波数応答解析でできること
 4. 周波数応答解析の解法:フル法(直接法)
 5. 周波数応答解析の解法:モード重ね合わせ法
 6. Ansysによる周波数応答解析
 7. Ansys周波数応答解析事例
 8. 付録
 9. おわりに
(3) 減衰の基礎
 1. はじめに
 2. 減衰の基礎
 3. 減衰のタイプ
 4. 粘性減衰
 5. 構造減衰
 6. モード減衰
 7. 各種減衰の使い分け
 8. 付録
 9. おわりに

1 はじめに

本シリーズ第2 回では周波数応答解析について取り上げました。周波数応答解析は加振力を考慮した振動応答を求めることができますが、得られる振幅は減衰の影響を大きく受けることをご説明いたしました。ところが減衰には様々な種類があり、解析タイプや解析手法によって使える減衰が異なるなど難解な部分が多いため、よく理解されずに利用されているケースも見受けられます。そこで本稿では、減衰の基礎とAnsysで利用できる各種減衰の特徴と使い分けをご説明いたします。

2 減衰の基礎

2.1 減衰は振動を抑える抵抗力

減衰とは何でしょうか?クルマのサスペンションについているダンパーを思い浮かべる方も多いと思いますが、そのように明確な形でなくても構造物には多かれ少なかれ減衰が含まれています。もし減衰がなければ、構造物はいったん振動すればずっと振動し続けるはずです。減衰によりエネルギーが消費され、振動エネルギーが小さくなっているのです。実際の振動現象に減衰が含まれている以上、CAEで解析する際にも減衰を考慮しなければ正しい解析になりません。

2.2 抵抗力は速度に比例

減衰は構造物に何らかの振動(動き)がないと働きませんので、減衰による抵抗力は速度に比例すると考えるのが簡単です。そこで減衰を速度に比例する粘性減衰として表現します。粘性減衰力fig1を(式1)で表し、比例定数fig4を(粘性)減衰係数と呼びます(図1)。

fig2

図1 粘性減衰

2.3 減衰するさまを臨界減衰比で評価

減衰fig4という物理量は、質量fig5や剛性fig6と違い非常にイメージしにくいです。そこで振動するかしないかの臨界減衰係数fig7と減衰係数fig4の比である臨界減衰比fig8を算出して評価するのが一般的です。簡単な1 自由度粘性減衰系振動モデル(図2)で臨界減衰比を算出してみます。

図2 1自由度粘性減衰系振動モデル

運動方程式は次式となります。

fig10

第1回のモーダル解析で行ったのと同じように、変位をfig11と仮定し代入します。

fig12

この式の解としてfig13単に振動しないことを意味するため除外すると次式となります。

fig14

fig15の部分に着目すると、fig4fig16より小さいか大きいかにより実数または虚数に分かれることがわかります。ここでfig17、固有振動数fig18と定義すると次式になります。

fig19

  • fig20のとき
    fig21は複素数となります。ここでfig22と定義するとfig23となり変位は(式6)で表現できます。fig24は振動状態を表し、fig25は時間が進むにつれ減衰する状態を表しますので減衰しつつ振動する結果となります。グラフで表すと(図3-b)です。

fig26

  • fig27のとき
    fig21は実数となります。fig28となり変位は(式7)で表現できます。今度は虚数成分がありませんので、時間が進むにつれ減衰する結果となります。グラフで表すと(図3-c)です。
fig29
 
つまりfig8が1より小さいか大きいかにより、振動するかしないかが分かれることになります。この境界となる減衰係数fig30臨界減衰係数と呼び、臨界減衰係数との比であるfig8を臨界減衰比と呼びます(式8, 9)。
 
fig32
臨界減衰比による振動の違いをまとめると(図3)のようになります。


減衰の表し方としては、臨界減衰比以外にも次にご紹介する対数減衰率や損失係数もよく使われます。混同しないよう注意してください。

図3 臨界減衰比による振動の違い

2.4 対数減衰率

対数減衰率fig34は減衰自由振動の振幅比(fig35 )の自然対数をとったもので、臨界減衰比fig8と以下の関係式が成り立ちます(図4)。

図4 対数減衰率

2.5 損失係数

損失係数とは、減衰の定義の1 つで、ヒステリシス曲線における最大変位時の力と変位ゼロ時の力の比から求めます。記号fig37がよく使われます。複素平面による表現方法ではfig38とも呼ばれます。臨界減衰比fig8と以下の関係があります。
fig39

2.6 エネルギーと臨界減衰比

系に減衰が存在すると振動時にエネルギーが消費されます。持っていたエネルギーfig40と1周期間に消費されるエネルギーfig41の比は以下の関係式が成り立ちます。

fig42

3 減衰のタイプ

すべてのメカニズムが解明されているわけではありませんが、減衰は以下の3つのタイプに分けることができます。

  • 粘性減衰
    空気抵抗や構造内部の摩擦抵抗などによって生じる減衰力で速度に比例します。
    Ansysではレイリー減衰などで表現します。
  • 構造減衰(履歴減衰、複素減衰とも呼ばれる)
    構造内部に発生する減衰力でひずみ速度にほとんど依存しません。材料固有の減衰特性を持ちます。
    Ansysでは一定構造減衰などで表現します。
  • クーロン減衰
    固体接触面に生じる摩擦による減衰力です(本稿では紙面の都合上取り扱いません)。

4 粘性減衰

粘性減衰は速度に比例する減衰系で、数式上、速度に対数減衰率 比例した粘性減衰系で表現します。

4.1 要素粘性減衰

もっとも直接的に減衰を表現するのが要素粘性減衰です。バネマス系におけるダッシュポットの役割に相当するもので、系にダンパーが部品として組み込まれているような場合に適します。Ansysではバネを定義すると減衰係数も合わせて指定できます(図5)。

図5 要素粘性減衰

4.2 α 減衰

α減衰は、構造物を取り巻く流体などの粘性による減衰を表したものです。質量に比例する質量比例減衰で、低周波側で過減衰に働きます(図6)。

図6 α減衰

4.3 β減衰

β減衰は、構造物を内部で発生する微小な摩擦による減衰を表したものです。剛性に比例する剛性比例減衰で、高周波側で過減衰に働きます(図7)。

図7 β減衰

4.4 レイリー減衰

レイリー減衰は、質量と剛性に比例する比例粘性減衰で、α とβ 減衰の両方を考慮します。2つの主要な振動モードにおける減衰比からα、β を算出します(図8)。
レイリー減衰はα 減衰とβ 減衰をともに定義することで設定します。(図9)は周波数応答解析で系全体にレイリー減衰を設定するときの例です。(フル法では)材料ごとに指定することもできます。

図8 レイリー減衰

図9 レイリー減衰の設定

5 構造減衰

構造減衰は履歴減衰、複素減衰とも呼ばれ、速度に関係なく減衰比fig8が一定です。減衰として取扱いやすいのが特徴で、周波数応答解析でよく利用されます。計算上は周波数依存性を持つ等価粘性減衰で表現されます。

5.1 等価粘性減衰係数

等価粘性減衰係数は、粘性減衰力が1 周期でなす仕事量が消費エネルギーfig41に変換され描くループ面積と同等であることを利用して求められます(図10)(式12)。

図10 強制振動時の履歴曲線

fig49

また、最大変位fig50に達したときに蓄えられるひずみエネルギーfig52は次式で求められます。

fig53

これらの式より粘性減衰力が1周期でなす仕事量と等価な減衰fig54を求めると(式15)となります。このfig54が等価粘性減衰係数です。

fig55

5.2 調和加振における運動方程式

等価粘性減衰係数を用いて、調和加振の運動方程式を記述すると次式となります。
fig56
 
ここで変位fig57となるため、次式に書き換えることができます。
fig58
式17)を見ると、減衰項がfig59という複素数のバネ(複素剛性)で表現されていることがわかります(図11)。
式からも明らかなように、粘性減衰とは異なり速度(振動数)に関係なく減衰比fig8が定義できます。このfig60を構造減衰係数または損失係数と呼び、周波数応答解析でよく用いられます。

図11 1自由度構造減衰系振動モデル

5.3 一定構造減衰

構造減衰は周波数によらず一定なため、一定構造減衰とも呼ばれます。(図12)はフル法周波数応答解析で系全体に一定構造減衰係数を設定するときの例です。
図13)は材料ごとに一定構造減衰を定義する例です。上段の減衰比(臨界減衰比)を入力すると、自動的に2 倍した値が下段の一定構造減衰係数(損失係数)に表示されます。

図12 一定構造減衰係数の設定

図13 材料別の一定構造減衰係数の設定

6 モード減衰

モード重ね合わせ法では、減衰マトリクス[C]を作成する代わりに、モードごとの臨界減衰比であるモード減衰比を設定します(図14)。モード減衰比は速度に依存しない減衰で、構造減衰とよく似ています。

図14 モード減衰比

図15)はモード重ね合わせ法周波数応答解析で系全体に一定のモード減衰比を設定するときの例です。材料ごとに設定することもできます。この例では全モードで同じ減衰比を使用していますが、モードごとに異なる減衰比を設定することもできます(コマンドが必要)。

図15 モード減衰比の設定

7 各種減衰の使い分け

ここまで様々な減衰をご紹介してきましたが、これだけ種類があると、どの減衰を用いればよいか迷ってしまうと思います。そこで目的や解析タイプ別の一般的な指針をご案内します。

7.1 目的による使い分け

  • とりあえず何らかの減衰が入ればよい場合
    一定構造減衰(モード重ね合わせ法では一定のモード減衰比)を設定するのが簡単です。減衰はゼロにしないことを推奨します。ほぼ減衰が無いと想定される場合は、小さな値(0.001 など)を設定してください。
  • 各部品や材料の減衰を個別に設定したい場合
    フル法の場合:材料ごとに臨界減衰比を設定するのが簡単です。(図13)を参照してください。
    モード重ね合わせ法の場合:(図15)の[モーダルからのEQV減衰比]をYesに設定してください。
  • 実験モーダル結果からモードごとの減衰が見積もれる場合
    MDAMPコマンドを利用して、モードごとに減衰を設定することができます。詳細はAnsysのコマンドリファレンスをご覧ください。
  • 材料の周波数ごとの減衰データをお持ちの場合
    これは主に粘弾性材料が該当すると思われます。粘弾性材料では周波数ごとの減衰比を定義するのが簡単です。詳細は弊社サポートへお問合せください。

7.2 解析タイプによる使い分け

  • フル法周波数応答解析の場合
    構造減衰がよく利用されますが、粘性減衰も使用できます。両者をともに使用すると効果が累積されます。
    粘弾性材料はフル法周波数応答解析のみ対応します。
  • モード重ね合わせ法周波数応答解析の場合
    モード減衰がよく利用されますが、レイリー減衰も使用できます。両者をともに使用すると効果が累積されます。
    要素減衰は使用できません(QR減衰法に変えることで対応可)。

8 補足

本稿ではフル法およびモード重ね合わせ法の周波数応答解析で用いられる減衰に絞って説明しています。時刻歴応答解析や他の解析タイプで用いられる減衰については、Ansysヘルプをご覧ください。
なお本稿でご紹介したAnsysのGUIはAnsys 2019 R2をもとに記載しております。バージョンにより多少異なる場合がありますので、ご注意ください。

9 おわりに

今回は減衰をご説明いたしました。CAEで利用する減衰には様々な種類があり場面に応じて適切な減衰を選択するのは大変ですが、振動解析を実施する上で避けて通れないため少々難解な部分まで解説いたしました。皆様の業務に少しでもお役立ていただけましたら幸いです。
全3 回にわたって振動解析の基礎をご紹介してまいりましたが、振動解析は非常に奥深い世界で、本シリーズで触れていない機能やテクニックが数多く存在します。弊社では豊富な知見をベースにセミナーや受託解析、コンサルティングなどのサービスを提供しております。もし振動解析でお困りの点がございましたら、弊社営業部にお問い合わせください。

参考文献

[ 1] Ansysオンラインヘルプ > 構造解析ガイド > 構造解析の概要 > 減衰

Ansys、ならびにANSYS, Inc. のすべてのブランド名、製品名、サービス名、機能名、ロゴ、標語は、米国およびその他の国におけるANSYS, Inc. またはその子会社の商標または登録商標です。その他すべてのブランド名、製品名、サービス名、機能名、または商標は、それぞれの所有者に帰属します。本ウェブサイトに記載されているシステム名、製品名等には、必ずしも商標表示((R)、TM)を付記していません。 CFX is a trademark of Sony Corporation in Japan. ICEM CFD is a trademark used by Ansys under license. LS-DYNA is a registered trademark of Livermore Software Technology Corporation. nCode is a trademark of HBM nCode.