動的エネルギー解析

動的エネルギー解析（DEA）^[1]は、複雑な構造物における構造伝播音と振動を数値的にモデル化する手法です。中〜高周波数範囲に適用可能で、この領域では従来の決定論的アプローチ（有限要素法や 境界要素法など）よりも計算効率が優れています。統計的エネルギー解析（SEA）^[2]などの従来の統計的アプローチと比較すると、 DEAはより多くの構造詳細を提供し、サブシステム分割に関する問題が少なくなります。DEA法は、（線形）輸送方程式の観点から複雑な構造物全体の振動波エネルギーの流れを予測します。これらの方程式は離散化され、メッシュ上で解かれます。

• 数値音響学における高周波法。
• エネルギーの流れはメッシュ全体で追跡されます。個々の光線ではなく光線の密度を使用するレイ トレーシングと考えることができます。
• 既存の FEM メッシュを使用できます。再モデリングは必要ありません。
• 計算時間は周波数に依存しません。
• 必要なメッシュ解像度は周波数に依存せず、FEM よりも粗く選択できます。ジオメトリを解決するだけで済みます。
• サブシステムごとに 1 つの数値のみを提供する SEA とは対照的に、微細な構造の詳細を解決できます。
• DEA で使用できるモデルの柔軟性が向上しました。SEA のような暗黙の仮定 (弱く結合されたサブシステムの平衡) はありません。

複雑な構造物（自動車、船舶、飛行機など）の振動音響特性のシミュレーションは、さまざまな設計段階で日常的に行われています。低周波数の場合、確立された選択肢は有限要素法（FEM）です。しかし、FEM を使用した高周波解析では、より短い波長を捉えるために車体構造の非常に細かいメッシュが必要になるため、計算コストが極めて高くなります。さらに、高周波での構造応答は、材料特性、形状、境界条件の小さな変化に非常に敏感です。このため、単一の FEM 計算の出力の信頼性が低くなり、アンサンブル平均が必要になるため、計算コストがさらに高くなります。したがって、高周波では、計算効率のよい他の数値手法が好まれます。

統計的エネルギー解析（SEA）^[2] は高周波問題を扱うために開発され、比較的小型で単純なモデルを生み出します。しかし、SEAは検証が難しい一連の仮定に基づいており、弱く結合した（弱く減衰した）サブシステム内で拡散波場と波エネルギーの準平衡が事実上必要とされます。

SEA の代替案の 1 つは、高周波限界における元の振動波の問題を考慮し、構造振動の光線追跡モデルを作成することです。^{[注 1]} 多重反射にわたる個々の光線の追跡は、軌道の増殖のため計算上実行可能ではありません。代わりに、転送演算子によって伝播される光線の密度を追跡する方がよい方法です。これは、参考文献で紹介されている動的エネルギー解析 (DEA)法の基礎となります。^[3] DEA は、拡散場と十分に分離されたサブシステムの仮定を取り除いた SEA の改良版と見なすことができます。位置と運動量の両方に依存するエネルギー密度を使用します。DEA は、エネルギーが隣接するメッシュ セル間を自由に流れる比較的細かいメッシュで動作できます。これにより、SEA によって課せられる制限と比較して、DEA で使用されるモデルの柔軟性が大幅に高まります。DEA は FE 解析用に作成されたメッシュを使用できるため、SEA のようなリモデリングは必要ありません。その結果、SEA よりも細かい構造の詳細を DEA で解決できます。

メッシュ上での DEA の実装は、離散フロー マッピング (DFM)と呼ばれます。ここでは、DFM の背後にある考え方について簡単に説明します。詳細については、以下の参考文献^{[1] [3] [4] [5] [6] [7]} を参照してください。DFM を使用すると、マルチモーダル伝播や曲面を含む、高周波での複雑な構造における振動音響エネルギー密度を計算できます。DFM はメッシュ ベースの手法であり、転送演算子を使用して構造のサブシステムの境界を通るエネルギーの流れを記述します。エネルギーの流れは光線の密度で表されます。つまり、特定の表面を通るエネルギー フラックスは、方向 の 点で表面を通過する光線の密度によって与えられます。ここで、は表面をパラメーター化し、は表面に接する方向成分です。以下では、表面は、自動車の床を記述する FE メッシュのメッシュ セルのすべての境界の結合によって表されます。位相空間座標 の密度 は、境界積分演算子 を介して、1 つの境界から隣接する境界交差点に転送されます。 ${\displaystyle \rho}$${\displaystyle s}$${\displaystyle p}$${\displaystyle s}$${\displaystyle p}$${\displaystyle \rho (s,p)=\rho (X_{s})}$${\displaystyle X_{s}=(s,p)}$

ここで、は、方向を持つ点の境界セグメントから始まる光線が別の境界セグメントを通過する場所を決定するマップであり、は 減衰係数と反射/透過係数を含む係数です（SEAの結合損失係数に似ています）。また、面内波と曲げ波の両方の場合のモード変換確率を制御します。これは、波の散乱理論から導き出されます（^[8]を参照）。これにより、DEAは曲率と変化する材料パラメータを考慮に入れることができます。式（1）は、1つのメッシュセル全体の光線追跡を、1つの表面から隣接する表面にエネルギー密度を転送する積分方程式で記述する方法です。 ${\displaystyle \phi (X_{s})}$${\displaystyle s}$${\displaystyle p_{s}}$${\displaystyle w(X_{s})}$

次のステップでは、伝達演算子（１）は位相空間の基底関数の集合を用いて離散化される。行列が構築されると、各要素の境界位相空間上の最終エネルギー密度 は、初期密度に関して、 次の形式の線形システムの解によって 与えられる。${\displaystyle {\bf {B}}}$${\displaystyle \rho}$${\displaystyle \rho_{0}}$

初期密度は、船舶のエンジンなどの振動励起のいくつかのソース分布をモデル化します。最終密度(すべてのセル境界上のエネルギー密度を記述) が計算されると、構造内の任意の場所のエネルギー密度を後処理ステップとして計算できます。 ${\displaystyle \rho_{0}}$${\displaystyle \rho}$

用語に関して、「離散フロー マッピング (DFM)」と「動的エネルギー解析」という用語には曖昧さがあります。ある程度、一方の用語を他方の用語の代わりに使用できます。たとえば、プレートを考えてみます。DFM では、プレートを多数の小さな三角形に分割し、三角形から (隣接する) 三角形にエネルギーの流れを伝播します。DEA では、プレートを分割せずに、プレートの境界でいくつかの高次基底関数 (位置と運動量の両方) を使用します。ただし、原則として、両方の手順を DFM または DEA として記述することは許容されます。

応用例として、自動車床パネルのシミュレーション^{[9] [10]} を示します。0.04 ヒステリシス減衰の 2500 Hz での点励起が適用されました。周波数平均 FEM シミュレーションの結果を DEA シミュレーションと比較します (DEA の場合、周波数平均化は不要)。結果は定量的にも良好な一致を示しています。特に、エネルギー フローの方向依存性がわかります。これは、プロットされているように主に水平方向です。これは、水平方向に伸びた複数の面外膨らみによって発生します。FEM 予測と DFM 予測の偏差が見られるのは、エネルギー量が無視できるパネルの右下部分のみです。DFM 予測によって与えられた総運動エネルギーは、FEM 予測の 12% 以内です。詳細については、引用文献を参照してください。

より応用的な例として、ヤンマートラクターモデル（青色のボディ：シャーシ/キャビンのスチールフレームと窓） のDEAシミュレーション^[11]の結果を左に示します。引用された研究では、ギアケースの背面の励起について、数値DEA結果が400Hzから4000Hzの周波数での実験測定と比較されています。両方の結果は良好に一致しています。DEAシミュレーションを拡張して、運転者の耳の音圧レベルを予測することができます。

• ノッティンガム大学波動モデリング研究グループ 2017-04-04 にWayback Machineでアーカイブされています。この研究グループの焦点の 1 つは DEA です。