音響-構造相互作用
概要
先生、音響-構造相互作用って、つまり音と構造が互いに影響し合うってことですか？
理論と物理
その通りだ。構造物が振動すると周囲の空気や流体を押して音波が発生する。逆に音圧が構造に荷重として作用して振動を引き起こす。この双方向のカップリングが音響-構造相互作用（ASI）だよ。
具体的にはどんな場面で問題になるんですか？
自動車の車室内騒音が代表例だね。エンジンの振動がボディパネルを励振し、パネルが車室内の空気を圧縮・膨張させて乗員が聞く音になる。他には航空機の胴体透過騒音、建築物の遮音設計、潜水艦のソナードーム設計などがある。
支配方程式
構造側と音響側、それぞれどんな方程式で記述するんですか？
構造側は通常の弾性体の運動方程式だ。
右辺の $[A]\{p\}$ が音圧から構造への荷重項で、$[A]$ は連成マトリクスだ。音響側はヘルムホルツ方程式の時間領域版、つまり波動方程式になる。
この2つはどこで繋がるんですか？
連成界面での境界条件で繋がる。音響側から見ると、界面の法線方向速度が構造の変位速度と一致する必要がある。
これを行列形式で書くと、連成系全体は非対称な方程式系になるんだ。
非対称なのが厄介そうですね。
そうなんだ。NastranではFLUID要素（CAABSF, CHEXAなど）を使うと内部でこの非対称連成系を組んでくれる。AbaqusではAcoustic-Structural連成のための*TIE制約が用意されているよ。
周波数領域の定式化
振動騒音って周波数で議論することが多いですよね？
調和応答を仮定すると $\{u\} = \{\hat{u}\}e^{i\omega t}$、$\{p\} = \{\hat{p}\}e^{i\omega t}$ として周波数領域に変換できる。
周波数ごとに解くので、周波数応答関数（FRF）を効率よく求められるんだ。
この定式化だと、どのモードがどの周波数で問題になるか見やすいですね。
まさにそう。特に構造の固有振動数と音響キャビティの固有振動数が近い場合、連成共振が起きて音圧が増幅される。車室内のブーミング問題がこの典型だよ。
実務上のポイント
解析する上で気をつけるべきことは何ですか？
まず音響メッシュの要素サイズだ。最高解析周波数 $f_{max}$ に対して、要素あたり最低6節点が波長内に入るようにする。
空気中で $c = 343$ m/s、$f_{max} = 1000$ Hz なら $h \leq 57$ mm になる。二次要素ならもう少し粗くてもよいが、一次要素では守った方がいい。
なるほど、低次要素だとかなり細かいメッシュが必要になりますね。全体像がだいぶ見えてきました。
「音が構造を変形させる」——音響放射圧という現象
音響-構造相互作用の理論を学ぶと、「音圧が構造に力を与える」という事実に改めて驚く人が多い。実は超音波洗浄機はまさにこの原理で動いていて、28〜40 kHzの音波が金属部品の表面に微細な気泡（キャビテーション）を生じさせ、汚れを物理的に剥がす。音圧レベルは150 dB以上にもなる。人間の耳には聞こえないのに、金属を変形させるほどのエネルギーを持つ——音響連成の理論はそんな「見えない力」の正体を方程式で記述するところから始まる。
各項の物理的意味
- 構造-熱連成項：温度変化による熱膨張が構造変形を誘発し、変形が温度場に影響する。$\sigma = D(\varepsilon - \alpha \Delta T)$。【日常の例】夏に線路のレールが伸びて隙間が狭くなる——温度上昇→熱膨張→応力発生の典型例。電子基板がはんだ付け後に反るのも、異なる材料の熱膨張率差による。エンジンのシリンダーブロックは高温部と低温部の温度差で熱応力が発生し、最悪の場合亀裂に至る。
- 流体-構造連成（FSI）項：流体圧力・せん断力が構造を変形させ、構造変形が流体領域を変化させる双方向の相互作用。【日常の例】強風で吊り橋のケーブルが振動する（渦励振）——風の力が構造を揺らし、揺れた構造が風の流れを変え、さらに振動が増幅する。心臓の血流と血管壁の弾性変形、航空機の翼のフラッタ（空力弾性不安定性）も典型的なFSI問題。片方向のみの連成で済む場合もあるが、変形が大きい場合は双方向連成が必須。
- 電磁-熱連成項：ジュール発熱 $Q = J^2/\sigma$ が温度上昇を引き起こし、温度変化が電気抵抗を変化させるフィードバックループ。【日常の例】電気ストーブのニクロム線は電流が流れると発熱（ジュール熱）して赤くなる——温度が上がると抵抗が変わり、電流分布も変化する。IHクッキングヒーターの渦電流発熱、送電線の温度上昇による弛み増加もこの連成の例。
- データ転写項：異なる物理場間のメッシュ不一致を補間で解決。【日常の例】天気予報で「気温のデータ」と「風のデータ」を合わせて体感温度を計算するとき、それぞれの観測地点が異なれば補間が必要——CAEの連成解析でも、構造メッシュとCFDメッシュは一般に一致しないため、界面でのデータ転写（補間）精度が結果の信頼性に直結する。
仮定条件と適用限界
- 弱連成仮定（片方向連成）：一方の物理場が他方に影響するが逆は無視可能な場合に有効
- 強連成が必要なケース：FSIでの大変形、電磁-熱連成での温度依存性が強い場合
- 時間スケールの分離：各物理場の特性時間が大きく異なる場合、サブサイクリングで効率化可能
- 界面条件の整合性：連成界面でのエネルギー・運動量保存が数値的に満たされることを確認
- 適用外ケース：3つ以上の物理場が同時に強く連成する場合、モノリシック手法が必要になることがある
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 熱膨張係数 $\alpha$ | 1/K | 鋼: 約12×10⁻⁶、アルミ: 約23×10⁻⁶ |
| 連成界面力 | N/m²（圧力）またはN（集中力） | 流体側と構造側で力の釣り合いを確認 |
| データ転写誤差 | 無次元（%） | 補間精度はメッシュ密度比に依存。5%以下が目安 |
数値解法と実装
FEMによる音響-構造連成の離散化
音響-構造連成って、FEMでどう扱うんですか？ 構造と音響で要素が違いますよね。
いい質問だ。構造側は通常のシェル要素やソリッド要素、音響側は圧力自由度を持つ音響流体要素を使う。連成界面では構造の変位自由度と音響の圧力自由度を結合させるんだ。
Nastranだとどうやるんですか？
NastranではFLUID要素（CHEXA, CPENTA, CTETRAの音響版）を使い、構造面にACMODAL面を定義する。SOL 108（周波数応答）やSOL 111（モーダル周波数応答）で連成系を解くんだ。
FEM-BEM連成法
外部放射問題はFEMだと無限領域が困りますよね？
その通り。内部問題（車室内など閉空間）ならFEMだけでいいが、外部放射問題ではBEM（境界要素法）と組み合わせる。構造振動をFEMで解き、放射面の速度をBEMの境界条件にする。
外部放射問題はFEMだと無限領域が困りますよね？
その通り。内部問題（車室内など閉空間）ならFEMだけでいいが、外部放射問題ではBEM（境界要素法）と組み合わせる。構造振動をFEMで解き、放射面の速度をBEMの境界条件にする。
BEMのKirchhoff-Helmholtz積分方程式は次の通りだ。
ここで $G = \frac{e^{ikR}}{4\pi R}$ が3次元自由空間のグリーン関数だ。
BEMだと体積メッシュが不要だから、外部問題に向いてるわけですね。
そう。Siemens Simcenter 3D（旧LMS Virtual.Lab）やFFT ASTRAがFEM-BEM連成の代表的なツールだ。HEAD acousticsのARTEMISなども境界要素ベースだね。
統計的エネルギー解析（SEA）
高周波になるとFEMのメッシュが膨大になりませんか？
鋭い。高周波域（自動車なら概ね500 Hz以上）ではモード密度が極めて高くなり、FEMの決定論的アプローチは現実的でなくなる。そこで統計的エネルギー解析（SEA）を使う。
SEAの基本方程式はパワーバランスだ。
ここで $\eta_{ij}$ は結合損失係数、$n_i$ はサブシステム $i$ のモード密度、$E_i$ はエネルギーだ。
じゃあ周波数帯域によって使い分けるんですね。
そう。低周波はFEM、高周波はSEA、中間帯域はハイブリッドFE-SEA（Simcenter 3DのHybrid FE-SEAやWaveが対応）。この使い分けが振動音響解析の実務の基本だよ。
| 周波数帯域 | 手法 | 代表ツール |
|---|---|---|
| 低周波（~500 Hz） | FEM / FEM-BEM | Nastran, Abaqus, COMSOL |
| 中周波 | ハイブリッドFE-SEA | Simcenter 3D, Wave6 |
| 高周波（500 Hz~） | SEA | VA One (ESI), AutoSEA |
時間領域の陽解法
衝撃音とか過渡的な音はどうするんですか？
時間領域で解くことになる。LS-DYNAは音響-構造連成の陽解法をサポートしていて、エアバッグ展開音やドア閉め音の過渡解析に使われるよ。Abaqus/Explicitでも音響要素との連成が可能だ。
手法ごとに得意分野がはっきりしてるんですね。振動音響の世界、面白いです。
FEMと境界要素法——音響連成解析の「相棒」選び
音響-構造相互作用の数値解法を選ぶとき、現場エンジニアは必ず悩む。「FEMだけで閉じた系（キャビティ）を解くか、外部音場にBEMを組み合わせるか」という問題だ。自動車の車室内騒音のように空間が閉じている場合はFEMが得意だが、車外への放射音を求めようとすると無限遠の境界処理が必要になる。そこでBEMの出番。ただしBEMは係数行列が密行列になるため、モデルが大きくなると計算コストが一気に跳ね上がる。近年のFMM（高速多重極法）の登場でこの問題はかなり緩和されたが、ソルバー選定は今でも「一長一短」の世界だ。
モノリシック法
全物理場を1つの連立方程式系として同時に解く。強い連成に対して安定だが、実装が複雑でメモリ消費が大きい。
パーティション法（分離反復法）
各物理場を独立に解き、界面でデータ交換。実装が容易で既存ソルバーを活用可能。弱い連成に適する。
界面データ転写
最近傍法（最も簡単だが精度低い）、射影法（保存的）、RBF補間（メッシュ非一致に強い）。保存性と精度のバランスが重要。
サブイタレーション
各連成ステップ内で十分な反復を行い、界面条件の整合性を確保。残差基準は各物理場の典型値に基づいてスケーリング。
Aitken緩和
連成反復の緩和係数を自動調整。過緩和による発散を防止し、収束を加速する適応的手法。
安定性条件
added mass効果（流体-構造連成で構造密度≈流体密度の場合）に注意。不安定な場合はロビン型界面条件やIQN-ILS法を適用。
Aitken緩和のたとえ
Aitken緩和は「シーソーのバランス取り」に似ている。一方が強く押しすぎると反対側が跳ね上がり、その反動でまた強く押しすぎる——この振動を抑えるために、押す力を自動的に調整するのがAitken緩和。連成反復が振動して収束しないとき、前回の修正量を見て次の修正量を自動調整する適応的手法。
実践ガイド
一般的なフローはこうだ。
1. 構造モデル作成 -- CADから構造FEMモデルを構築。パネル結合部やシール部の接続を正確にモデル化
2. 音響キャビティモデル作成 -- 車室内などの閉空間を音響要素でメッシュ化
3. 連成界面定義 -- 構造面と音響面のインターフェースを設定（NastranならACMODAL、AbaqusならTIE）
4. 加振条件設定 -- エンジンマウント力、路面加振力などの入力
5. 周波数応答解析実行 -- 対象周波数範囲でFRFを計算
6. 後処理 -- 音圧レベル分布、寄与解析、パネル寄与分析
ステップ3の連成界面定義が肝だと聞いたんですが。
その通り。界面の法線方向が一致していないと、エネルギーの授受が正しく計算されない。Nastranではwet surfaceの法線方向を揃えるのが定番のチェックポイントだ。
メッシュ設計のコツ
構造と音響でメッシュサイズの考え方が違いますよね？
音響メッシュは波長基準で $h \leq c/(6 f_{max})$。構造メッシュは曲げ波長基準で $h \leq \lambda_b / 6$ だ。曲げ波長は次の式で求まる。
ここで $D$ は曲げ剛性、$t$ は板厚だ。一般に構造側の方が細かいメッシュが必要になる。
構造と音響でメッシュが一致しない場合はどうするんですか？
非適合メッシュの連成になる。NastranのACMODALやAbaqusのTIEは非適合メッシュに対応している。ただしメッシュ差が大きすぎると連成精度が落ちるので、界面付近では2:1程度の比率に抑えたい。
パネル寄与解析
どのパネルがうるさいか特定する方法ってありますか？
パネル寄与解析（Panel Contribution Analysis）だね。受音点の音圧をパネルごとの寄与に分解する。パネル $j$ の寄与は次のように書ける。
位相と振幅の両方を含むので、あるパネルの寄与が他を打ち消している場合もある。Siemens Simcenter 3Dの寄与解析機能が有名だよ。
対策を打つパネルの優先順位がつけられるんですね。
よくある失敗と対策
| 症状 | 原因 | 対策 |
|---|---|---|
| 特定周波数で異常な音圧ピーク | 連成共振（構造と音響の固有振動数一致） | 構造側の固有振動数シフト（補強追加）または制振材追加 |
| 実測と解析の周波数がずれる | 材料定数（特にヤング率）の誤差 | 実験モーダルで同定した値を使う |
| 全周波数帯で音圧が過大 | 減衰の設定不足 | 構造減衰（$\eta_s$=0.01~0.03）と音響吸音率を適切に設定 |
| 高周波で解析結果が振動的 | メッシュが粗い（波長分解能不足） | 6要素/波長のルールを確認 |
減衰の設定って、意外と難しそうですね。
そう、実務では構造減衰係数をパネルごとに設定し、シートやカーペットの吸音率をインピーダンス境界条件として入れる。この調整が解析精度に大きく影響するよ。
ドア閉まり音は「音響-構造連成」の塊だった
高級車のドアを閉めるとき「バタン」ではなく「コトン」と鳴る——あの音質を設計する仕事がある。実はドア閉まり音は、金属パネルの振動（構造）→室内キャビティの空気振動（音響）→シール材の減衰（材料）が絡み合った、まさに音響-構造相互作用の応用例だ。実務では、パネルの板厚を0.1 mm変えるだけで音質が変わるため、試作なしに狙った音を出すために音響連成解析が使われる。「ドア音エンジニア」という専門職が存在するほど、現場では重要視されている。
解析フローのたとえ
風船を膨らませたことがありますか？ あの瞬間、実は高度な流体-構造連成が起きています。内部の空気圧（流体）がゴム壁（構造）を押し広げ→広がった壁が内部の圧力分布を変え→変わった圧力がさらに壁を変形させる…このキャッチボールを計算ステップごとに繰り返すのがFSI解析です。
初心者が陥りやすい落とし穴
「片方向連成で十分でしょ？」——この判断ミスが連成解析で最も危険です。構造の変形が微小なら確かに片方向で足りますが、心臓弁の開閉のように変形が流路を大きく変える場合、片方向では全く話になりません。目安は「変形量が代表長さの1%を超えるか」。超えるなら双方向連成は必須です。片方向で済ませてしまった場合、結果が「もっともらしいけど実は大間違い」になる——これが最も怖いパターンです。
境界条件の考え方
連成界面のデータ交換は「国境の出入国管理」と同じです。各国（物理場）には独自の法律（支配方程式）がありますが、国境（界面）で人や物（力・温度・変位）のやり取りを正確に管理しないと、両国の経済（エネルギーバランス）が崩壊します。メッシュが一致していない場合の補間は「通訳」のようなもの——誤訳（補間誤差）が小さいほど良い結果が得られます。
ソフトウェア比較
主要なものを整理すると以下の通りだ。
| ツール名 | 開発元 | 得意分野 | 音響手法 |
|---|---|---|---|
| MSC Nastran | Hexagon (旧MSC) | 自動車NVH | FEM音響、ACMODAL |
| Simcenter 3D | Siemens | NVH統合環境 | FEM-BEM, ハイブリッドFE-SEA |
| Actran | FFT (Siemens子会社) | 音響専門 | FEM, BEM, インフィニットエレメント |
| COMSOL Multiphysics | COMSOL AB | マルチフィジックス | FEM音響, PML |
| Abaqus | Dassault Systemes | 汎用非線形 | FEM音響-構造連成 |
| VA One | ESI Group | 中高周波SEA | SEA, FE-SEA |
Actranって初めて聞きました。
ActranはベルギーのFFT（Free Field Technologies）が開発した音響専門ソルバーだ。2019年にSiemensが買収した。無限要素（Infinite Element）や適応PML（Perfectly Matched Layer）が充実していて、外部放射問題に非常に強いよ。航空エンジンのファン騒音解析で広く使われている。
NVH業界でのツール棲み分け
自動車メーカーだとどういう使い方が多いんですか？
典型的な棲み分けはこうだ。
- 低周波ブーミング: Nastran SOL 111（モーダル周波数応答） + 音響キャビティFEM
- ロードノイズ: Nastran + 伝達パス解析（TPA）
- 風切り音: CFD（Fluent/STAR-CCM+） → Actran or Simcenter BEM
- 中高周波: VA One（SEA） or Simcenter Hybrid FE-SEA
NastranがNVHのデファクトスタンダードなんですね。
そう。DMAP（Direct Matrix Abstraction Program）で独自の処理も書ける柔軟性がある。ただし最近はSimcenter 3Dがpre/postを含めた統合環境として勢いがあるね。
オープンソース選択肢
商用ツール以外の選択肢はありますか？
OpenFOAMのacousticFoamソルバーでLES+音響アナロジーができる。Code_Aster+Salome-Mecaでも音響FEMは可能だ。学術分野ではpyFEMという音響FEMのPythonライブラリもあるよ。ただし商用ツールのような自動連成インターフェースは自前実装になる。
やはり商用ツールの方が生産性は高そうですね。ツール選定の参考になりました。
NastranとAbaqus——音響連成ツール選びの現実
音響-構造連成解析のツール選定で「NastranかAbaqusか」という議論は根深い。NastranはDMAPによるカスタマイズ性と軽量なモデルが強みで、自動車・航空宇宙の構造部門では圧倒的なシェアを持つ。一方Abaqusは材料非線形や複雑な接触条件の扱いに強く、複合材や樹脂部品が多い製品では選ばれやすい。音響解析単体ならActranやVirtualLabが専業ツールとして精度面で有利な場面もある。実務での選定基準は「既存の社内ライセンス」と「エンジニアのスキル資産」が9割を占める——ツールの性能差より、使いこなせるかどうかの方がはるかに重要だ。
選定で最も重要な3つの問い
- 「何を解くか」：音響-構造相互作用に必要な物理モデル・要素タイプが対応しているか。例えば、流体ではLES対応の有無、構造では接触・大変形の対応能力が差になる。
- 「誰が使うか」：初心者チームならGUIが充実したツール、経験者ならスクリプト駆動の柔軟なツールが適する。自動車のAT車（GUI）とMT車（スクリプト）の違いに似ている。
- 「どこまで拡張するか」：将来の解析規模拡大（HPC対応）、他部門への展開、他ツールとの連携を見据えた選択が長期的なコスト削減につながる。
先端技術
ハイブリッドFE-SEA
低周波はFEM、高周波はSEAと言ってましたが、中間帯域のハイブリッド手法って具体的にどういう仕組みですか？
ハイブリッドFE-SEAは、決定論的に記述すべき部品（エンジンブロックなど剛性が高く少数モードの構造）をFEサブシステムとして扱い、統計的に記述すべき部品（パネルなどモード密度の高い構造）をSEAサブシステムとして扱う。
理論的にはLangleyらのケンブリッジ大学グループが2005年前後に体系化した。ランダムな点質量を仮定して「拡散場相反定理（diffuse field reciprocity）」を使うことで、FEとSEAの整合的な結合を実現しているんだ。
使えるツールは？
Simcenter 3DのHybrid FE-SEA機能と、ESI GroupのVA Oneが二大ツールだ。Wave6（旧Wave、Dassault Systemes）も対応している。
トリム音響モデリング
車の内装材（トリム）って音響にどう影響するんですか？
フロアカーペット、ダッシュインシュレーターなどの多層吸遮音材はBiotのポーロ弾性体理論で記述する。支配方程式が3つの場（固体フレーム変位、流体変位、音圧）の連成になるので非常に重い。
ActranのTrim Acoustic Module（TAM）やCOMSOLのPoroacoustics Moduleがこれに対応している。実務ではJohnson-Champoux-Allardモデルなどの等価流体モデルで近似することが多いよ。
なるほど、多孔質材料の流れ抵抗率 $\sigma$ とかを測定して入れるわけですね。
機械学習との融合
最近のAI活用はどうなっていますか？
いくつかの方向がある。一つは音響伝達関数のサロゲートモデル化だ。大量のFEM解析結果からニューラルネットワークで伝達関数を学習し、設計変更時の応答をリアルタイムに予測する。
二つ目はPINN（Physics-Informed Neural Network）で波動方程式の解を直接近似する研究。まだ精度面で課題があるが、逆問題（音源同定）には有望とされている。
三つ目はトポロジー最適化との組み合わせ。音響パワーを最小化する構造形状をAdjoint法で求める研究が活発だ。COMSOLやOptiStructが対応しているよ。
CAEの世界もどんどん変わっていくんですね。先端動向を追いかけ続けないと。
音響メタマテリアル——「負の質量密度」という奇妙な概念
音響-構造相互作用の先端研究として注目されているのが「音響メタマテリアル」だ。特定の周波数帯で等価的に「負の質量密度」や「負の弾性率」を示す構造体で、音の遮断や屈折をこれまでとはまったく異なる原理で制御できる。例えば10 cm四方の板状メタマテリアルが、同面積の通常の防音壁を大幅に上回る遮音性能を発揮する事例も報告されている。「負の密度」という概念は直感に反するが、局所的な共振構造（スプリング-マス系）の集合体として設計すると実現できる。次世代の静粛性設計を変える技術として研究が加速している。
トラブルシューティング
順番に見ていこう。
1. 連成界面の法線方向不一致
症状: 音圧が全くのゼロになる、または期待値の数倍になる。
原因: 構造面の法線が音響キャビティの内側を向いていない。Nastranでは法線が外向き（音響側）でないとconnectivityが正しく取れない。
対策: プリプロセッサ（HyperMesh, ANSA等）で法線方向を可視化して統一する。NastranのBGPDTカードでチェック可能。
これ、何時間も原因が分からないやつですね。
2. モード抜け
症状: 特定の周波数帯で実測と大きく乖離。
原因: モーダル周波数応答（SOL 111等）で使用するモード数が不足。音響キャビティと構造の両方のモードを解析周波数上限の1.5倍まで含める必要がある。
対策: 残留ベクトル（residual vector）の使用。NastranのRESTVECパラメータ、AbaqusのRESIDUAL MODES。
3. 高周波での精度劣化
症状: 解析周波数上限付近でFRFが振動的になり信頼性が低い。
原因: 音響メッシュの波長分解能不足。6要素/波長のルールを満たしていない。
対策: メッシュ密度を見直すか、対象周波数範囲を狭める。二次要素を使えば同じ節点数で分解能が向上する。
4. 減衰パラメータの不整合
構造減衰と音響吸音って別々に設定するんですよね？
症状: 共振ピークが実測よりシャープすぎる、またはブロードすぎる。
原因: 構造の損失係数 $\eta_s$ と音響の壁面インピーダンスが実態と合っていない。
対策: 構造減衰は実験モーダルから同定する（半値幅法）。音響吸音率は残響室法やインピーダンス管で測定した値を使う。
| パラメータ | 典型値 | 備考 |
|---|---|---|
| 鋼板の損失係数 | 0.001~0.003 | 制振材なし |
| 制振材付き鋼板 | 0.01~0.05 | CLD処理後 |
| 車室内吸音率（低周波） | 0.05~0.15 | シート・カーペット含む |
| 車室内吸音率（高周波） | 0.2~0.5 | 同上 |
5. Nastranでの典型的エラー
連成解析特有のハマりポイントが分かって助かりました。実務でエラーが出たらこのリストを見返します。
「解析は収束したのに実測と全然違う」——よくある音響連成の落とし穴
音響-構造相互作用の解析で「計算は無事終わったのに実験値と10 dB以上ずれる」という悩みは現場で頻繁に聞く。原因として最も多いのは「境界条件の設定ミス」だ。実際の部品はボルトで固定されているが、解析では完全固定（固定端）として扱うと、実際の締結剛性との差が固有値をずらし、共鳴ピークの周波数がズレる。さらに、内部損失係数（減衰）の設定を実測値ではなく教科書値で入れると、ピーク付近で大きな誤差が出る。「モデルが合わない＝メッシュが粗い」と思いがちだが、多くの場合は境界条件と減衰値の見直しで解決する。
「解析が合わない」と思ったら
- まず深呼吸——焦って設定をランダムに変えると、問題がさらに複雑になる
- 最小再現ケースを作る——音響-構造相互作用の問題を最も単純な形で再現する。「引き算のデバッグ」が最も効率的
- 1つだけ変えて再実行——複数の変更を同時に行うと、何が効いたか分からなくなる。科学実験と同じ「対照実験」の原則
- 物理に立ち返る——計算結果が「重力に逆らって物が浮く」ような非物理的な結果なら、入力データの根本的な間違いを疑う
なった
詳しく
報告