CAEを学ぶ
構造解析
はじめての振動解析(3)減衰の基礎
はじめての振動解析(3)減衰の基礎 の概要
はじめての振動解析(1)(2)(3) 完全保存版 PDF を "無料" 公開中!
はじめての振動解析(1)(2)(3) 完全保存版 PDF サンプル
これさえ読めば振動解析の基本がわかる!
全3回のシリーズとなっております「はじめての振動解析」について、3回分の解析講座を1冊でお読みいただける完全保存版 PDF を公開いたしました。どなたでも簡単なフォーム入力のみでダウンロードいただけますので、是非この機会にご利用ください。
ページ数:30ページ
ファイル形式: PDF
はじめての振動解析(1)(2)(3) 完全保存版 PDF の目次
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 減衰の基礎
2.1 減衰は振動を抑える抵抗力
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の部分に着目すると、fig4がfig16より小さいか大きいかにより実数または虚数に分かれることがわかります。ここでfig17、固有振動数[画像:fig18]と定義すると次式になります。
[画像:fig19]
- fig20のとき
fig21は複素数となります。ここでfig22と定義するとfig23となり変位は(式6)で表現できます。fig24は振動状態を表し、fig25は時間が進むにつれ減衰する状態を表しますので減衰しつつ振動する結果となります。グラフで表すと(図3-b)です。
[画像:fig26]
- fig27のとき
fig21は実数となります。fig28となり変位は(式7)で表現できます。今度は虚数成分がありませんので、時間が進むにつれ減衰する結果となります。グラフで表すと(図3-c)です。
図3 臨界減衰比による振動の違い
2.4 対数減衰率
対数減衰率fig34は減衰自由振動の振幅比(fig35 )の自然対数をとったもので、臨界減衰比fig8と以下の関係式が成り立ちます(図4)。
図4 対数減衰率
2.5 損失係数
2.6 エネルギーと臨界減衰比
系に減衰が存在すると振動時にエネルギーが消費されます。持っていたエネルギーfig40と1周期間に消費されるエネルギーfig41の比は以下の関係式が成り立ちます。
[画像:fig42]
3 減衰のタイプ
すべてのメカニズムが解明されているわけではありませんが、減衰は以下の3つのタイプに分けることができます。
- 粘性減衰
空気抵抗や構造内部の摩擦抵抗などによって生じる減衰力で速度に比例します。
Ansysではレイリー減衰などで表現します。 - 構造減衰(履歴減衰、複素減衰とも呼ばれる)
構造内部に発生する減衰力でひずみ速度にほとんど依存しません。材料固有の減衰特性を持ちます。
Ansysでは一定構造減衰などで表現します。 - クーロン減衰
固体接触面に生じる摩擦による減衰力です(本稿では紙面の都合上取り扱いません)。
4 粘性減衰
4.1 要素粘性減衰
図5 要素粘性減衰
4.2 α 減衰
図6 α減衰
4.3 β減衰
図7 β減衰
4.4 レイリー減衰
図8 レイリー減衰
図9 レイリー減衰の設定
5 構造減衰
構造減衰は履歴減衰、複素減衰とも呼ばれ、速度に関係なく減衰比fig8が一定です。減衰として取扱いやすいのが特徴で、周波数応答解析でよく利用されます。計算上は周波数依存性を持つ等価粘性減衰で表現されます。
5.1 等価粘性減衰係数
等価粘性減衰係数は、粘性減衰力が1 周期でなす仕事量が消費エネルギーfig41に変換され描くループ面積と同等であることを利用して求められます(図10)(式12)。
図10 強制振動時の履歴曲線
5.2 調和加振における運動方程式
図11 1自由度構造減衰系振動モデル
5.3 一定構造減衰
図12 一定構造減衰係数の設定
図13 材料別の一定構造減衰係数の設定
6 モード減衰
モード重ね合わせ法では、減衰マトリクス[C]を作成する代わりに、モードごとの臨界減衰比であるモード減衰比を設定します(図14)。モード減衰比は速度に依存しない減衰で、構造減衰とよく似ています。
図14 モード減衰比
図15 モード減衰比の設定
7 各種減衰の使い分け
7.1 目的による使い分け
- とりあえず何らかの減衰が入ればよい場合
一定構造減衰(モード重ね合わせ法では一定のモード減衰比)を設定するのが簡単です。減衰はゼロにしないことを推奨します。ほぼ減衰が無いと想定される場合は、小さな値(0.001 など)を設定してください。 - 各部品や材料の減衰を個別に設定したい場合
フル法の場合:材料ごとに臨界減衰比を設定するのが簡単です。(図13)を参照してください。
モード重ね合わせ法の場合:(図15)の[モーダルからのEQV減衰比]をYesに設定してください。 - 実験モーダル結果からモードごとの減衰が見積もれる場合
MDAMPコマンドを利用して、モードごとに減衰を設定することができます。詳細はAnsysのコマンドリファレンスをご覧ください。 - 材料の周波数ごとの減衰データをお持ちの場合
これは主に粘弾性材料が該当すると思われます。粘弾性材料では周波数ごとの減衰比を定義するのが簡単です。詳細は弊社サポートへお問合せください。
7.2 解析タイプによる使い分け
- フル法周波数応答解析の場合
構造減衰がよく利用されますが、粘性減衰も使用できます。両者をともに使用すると効果が累積されます。
粘弾性材料はフル法周波数応答解析のみ対応します。 - モード重ね合わせ法周波数応答解析の場合
モード減衰がよく利用されますが、レイリー減衰も使用できます。両者をともに使用すると効果が累積されます。
要素減衰は使用できません(QR減衰法に変えることで対応可)。
8 補足
なお本稿でご紹介したAnsysのGUIはAnsys 2019 R2をもとに記載しております。バージョンにより多少異なる場合がありますので、ご注意ください。
9 おわりに
全3 回にわたって振動解析の基礎をご紹介してまいりましたが、振動解析は非常に奥深い世界で、本シリーズで触れていない機能やテクニックが数多く存在します。弊社では豊富な知見をベースにセミナーや受託解析、コンサルティングなどのサービスを提供しております。もし振動解析でお困りの点がございましたら、弊社営業部にお問い合わせください。
参考文献
[ 1] Ansysオンラインヘルプ > 構造解析ガイド > 構造解析の概要 > 減衰
関連情報
関連する解析事例
MORE電子銃内部の空間電荷効果の評価
Ansys Charge Plusを用いた電子銃内の電子温度、電子密度を可視化
パワーモジュールのさらなる信頼性向上へ、亀裂進展マクロを活用した寿命解析に挑む
計測技術で医療の未来を支える −磁気共鳴工学で実現する体に負担のない検査
北海道大学 大学院情報科学研究院 様
実験とシミュレーションの相互作用で現象の本質に迫る 〜教育現場で活きる Ansys Mechanical の事例
北海道大学 大学院工学研究院 様
そり変形を考慮した樹脂バックルのはめ合わせ接触解析事例
静止軌道のプラズマ環境下における衛星表面帯電
Ansys Charge Plus による帯電リスクの可視化
Ansys HFSSを用いた車載機器の放射イミュニティ規格試験解析
試験環境モデルを用いたISO 11452-2 解析
ラインを止めない!実機検証からシミュレーションへ、 デジタルツイン活用へ向けた取り組み事例
ロート製薬株式会社 様
関連する資料ダウンロード
MORE-
無料解析事例・技術資料
非線形解析の壁を突破するシミュレーション技術
〜Ansys LS-DYNAで解決!最大荷重評価のボトルネック〜
-
無料製品/サービス紹介
事例でご紹介!流体解析分野のエンジニアリングサービス 〜解析業務の委託・立ち上げ支援・カスタマイズによる効率化など〜
-
無料解析事例・技術資料
はんだ濡れ上がり形状予測解析
〜Ansys LS-DYNAで電子機器の信頼性向上に貢献〜
-
無料製品/サービス紹介CAE運用・活用
Ansys ユーザーのための PyAnsys 完全ガイド
Pythonで加速するCAEワークフロー
-
無料製品/サービス紹介
共振回避だけで終わらせない実レベルの振動解析
〜Ansys Mechanicalで実現する高度な製品開発〜
-
無料解析事例・技術資料
吸入器内の粒子挙動を可視化する
〜薬剤送達効率向上に向けた解析〜
-
無料解析事例・技術資料
医薬品バイアルの温度挙動解析
〜保管環境の影響把握と品質維持に向けた可視化アプローチ〜
-
無料製品/サービス紹介
そのFDTD計算、もっと速くできる!Lumerical+GPUでフォトニクス解析に革命を