計算式
波速:$c = \sqrt{E/\rho}$ m/sインピーダンス:$Z = \rho c A$
反射係数:$R = \dfrac{Z_B - Z_A}{Z_B + Z_A}$
透過係数:$T = \dfrac{2Z_B}{Z_B + Z_A}$
材料A・Bの棒の界面での応力波の反射・透過をx-t線図でリアルタイム可視化。反射係数R・透過係数Tをインピーダンス比から計算。
一次元の棒を伝わる縦波(圧縮波)の伝播速度は、材料の剛性(ヤング率)と密度で決まります。
$$c = \sqrt{\frac{E}{\rho}}$$$c$: 弾性波の伝播速度 [m/s], $E$: ヤング率 [Pa], $\rho$: 密度 [kg/m³]
異なる材料の界面での波の挙動は、音響インピーダンスの比によって決まり、反射係数と透過係数で表されます。
$$Z = \rho c A, \quad R = \frac{Z_B - Z_A}{Z_B + Z_A}, \quad T = \frac{2Z_B}{Z_B + Z_A}$$$Z$: 音響インピーダンス [kg/s], $A$: 断面積 [m²], $R$: 反射係数, $T$: 透過係数。$R$が負の時は応力の符号が反転(圧縮→引張)します。
材料試験(ホプキンソン棒試験):高速衝撃による材料の動的応力-ひずみ関係を評価する際の基礎理論です。試験片を挟んだ2本の棒を伝わる入射波・反射波・透過波を計測し、ここで学んだ反射・透過の関係式を用いて材料特性を求めます。
非破壊検査(超音波探傷):材料内部のきず(割れや空洞)を検出する技術です。きずはインピーダンスが異なる界面とみなせ、探触子から発信された超音波パルスがそこで反射して戻ってくる時間と強さを分析することで、きずの位置や大きさを推定します。
地震工学:地震波が地下の異なる地層(岩盤、土壌など)の境界で反射・透過・屈折する挙動を理解する基礎となります。地盤のインピーダンスの違いが、地表に到達する震動の大きさに大きく影響します。
音響デバイス設計:超音波振動子や音響フィルタなど、異種材料を接合して使うデバイスでは、境界での効率的なエネルギー(波)の透過や、不要な反射の抑制が重要です。インピーダンス整合の設計に本理論が応用されます。
このシミュレーターを使い始めるとき、いくつか勘違いしやすいポイントがあるよ。まず、「波の形が変わらない」と思いがちなこと。実際の世界では、材料の内部摩擦や幾何学的な広がり(例えば3次元の球面波)で波は減衰・分散するけど、この1Dモデルではそれが無視されている。だから、シミュレーション上では永遠にきれいな矩形波が往復するけど、それはあくまで理想化された挙動なんだ。
次に、パラメータ設定の現実性。例えば、密度を極端に小さくして「自由端反射」を再現するのは学習上は良いけど、実在する材料の密度やヤング率の範囲は限られている。例えば、鋼の密度は約7800 kg/m³、アルミニウムは約2700 kg/m³だ。材料Bに「空気」相当の値を入れると(密度約1.2 kg/m³)、インピーダンス比が極端になりR≒-1になる。これは理解を深める良い実験だけど、実務ではここまで極端な差の界面は「材料-真空」ぐらいしかないことに注意しよう。
最後に、「反射係数Rと透過係数Tの関係」。エネルギー保存則から、反射波と透過波のエネルギーの和は入射波のエネルギーに等しくなるはずだ。でも、ここで表示されているRとTは応力(または粒子速度)の比だから、単純にR+T=1にはならないんだ。エネルギーの反射率と透過率は、$R^2$と$(Z_A/Z_B)T^2$で計算される。シミュレーターでR=0.5, T=1.2のような一見不思議な値が出ても、エネルギー計算をしてみると辻褄が合うことを確認してみて。
この1D弾性波の基礎は、思ったよりずっと幅広い分野の根底にあるんだ。まず挙げるのは「音響工学と超音波加工」。超音波洗浄機や超音波溶着機では、振動子(圧電素子)から発生した弾性波をホーン(工具)に効率よく伝えないといけない。異なる材料を接合する部分では、まさにこのシミュレーターで学んだインピーダンス整合が設計の鍵になる。インピーダンスが合わないとエネルギーが反射されてしまい、加工部に十分な振幅が得られないんだ。
もう一つは「地盤工学と耐震設計」。地震のS波やP波は地中を伝わる弾性波そのものだ。硬い岩盤から軟らかい沖積層へ波が伝わるとき、インピーダンスの違いで反射・透過が起きる。ここで透過係数が大きい(=多くのエネルギーが地表に伝わる)と震動が増幅され、建物への影響が大きくなる。このシミュレーターで材料Bを材料Aよりずっと軟らかく(Eを小さく)設定すると、その状況をシンプルにモデル化できるよ。
さらに「生体医工学」への応用も面白い。超音波診断装置(エコー)では、探触子から体内へ超音波パルスを送る。皮膚、脂肪、筋肉、骨といった各組織の音響インピーダンスは全て異なる。骨のようにインピーダンスが極端に高い組織の境界では強い反射(エコー)が返ってくる。この原理を理解するのに、この1Dモデルは優れた第一歩になるんだ。
このシミュレーターに慣れたら、次のステップとして「波動方程式の数学」に触れてみよう。このシミュレーションの背後には、一次元の波動方程式 $\frac{\partial^2 u}{\partial t^2} = c^2 \frac{\partial^2 u}{\partial x^2}$ がある。ここで$u$は変位だ。この偏微分方程式を「ダランベールの解」という方法で解くと、$u(x,t) = f(x-ct) + g(x+ct)$ という「進行波」の形の解が得られる。シミュレーターで見ている右向きと左向きの波が、まさにこの $f$ と $g$ に対応しているんだ。
実務に近づくなら、「3次元への拡張とモードの理解」が重要だ。実際の構造物では、棒の軸方向の縦波(このシミュレーターで扱っている波)だけでなく、曲げ波(横波)やねじり波など、さまざまなモードが同時に発生し、互いに結合することもある。例えば、板に衝撃を与えたときに発生する「ラム波」や、パイプを伝わる「ガイド波」の解析は、非破壊検査の最先端で活用されている。
学習の具体的な次のトピックとしては、「数値解法(FDM、FEM)の初歩」を薦めるよ。このシミュレーターはおそらく解析解に近い方法で計算しているけど、複雑な形状や非線形材料を扱うには有限要素法(FEM)などの数値シミュレーションが必須だ。まずは一次元の波動方程式を、単純な有限差分法(FDM)で自分でプログラミングして解いてみると、時間ステップや空間分割の重要性、安定性条件(例えば、クーラン条件 $\Delta t \le \Delta x / c$)といった数値解析の基礎概念が体で理解できるようになる。そこから市販のCAEソフトが内部で何をしているのか、その本質が見えてくるはずだ。