弾性波伝播 戻る
衝撃・波動力学

弾性波伝播シミュレーター

材料A・Bの棒の界面での応力波の反射・透過をx-t線図でリアルタイム可視化。反射係数R・透過係数Tをインピーダンス比から計算。

材料 A(入射側)
材料
材料 B(透過側)
材料
荷重・形状
棒長さ L2.0 m
衝撃持続時間 τ100 μs
衝撃力 F100 kN
計算結果
波速 cA
m/s
波速 cB
m/s
反射係数 R
透過係数 T

計算式

波速:$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}$
x-t 線図(ラグランジュ線図)
入射波
反射波
透過波
応力分布(t = L/c, 2L/c, 3L/c)

弾性波伝播シミュレーターとは

🧑‍🎓
このシミュレーターで「弾性波」って何がわかるんですか?衝撃を与えたら、棒の端まで一気に変形するんじゃないの?
🎓
ざっくり言うと、衝撃は「波」として伝わっていくんだ。例えばハンマーで鉄の棒の端を叩くと、その変形は一瞬で反対側には届かず、波速 $c$ で伝播する。このシミュレーターでは、材料AとBの境目で波がどう反射・透過するかがリアルタイムで見られるよ。まずは上のスライダーで「衝撃力 F」を変えて、波の大きさがどう変わるか試してみて。
🧑‍🎓
え、そうなんですか!「反射係数 R」がマイナスになることもあるみたいですけど、これはどういう意味ですか?
🎓
良いところに気が付いたね。反射係数がマイナスというのは、波の「位相」が反転する、つまり圧縮波が引張波として反射されることを意味するんだ。実務で多い例は、鋼(材料A)から空気(材料B)のような非常にインピーダンスの小さいものへの境界だね。シミュレーターで材料Bの密度 $\rho$ を極端に小さく設定すると、Rが-1に近づいて「自由端反射」の状態が再現できるよ。
🧑‍🎓
なるほど!でも、この「x-t線図」ってどう見ればいいんですか?斜めの線がたくさんあって…。
🎓
あの図は、横軸が位置、縦軸が時間で、波面の動きを追跡した「道のり」みたいなものだ。傾きが波速の逆数になるから、材料が変わると傾きも変わるんだ。界面で線が折れ曲がるのが反射や透過の瞬間。パラメータで材料AとBのヤング率 $E$ や密度 $\rho$ を変えて、線の傾きと反射の様子がどう連動するか観察してみよう。超音波探傷の原理が直感的に理解できるはずだよ。

物理モデルと主要な数式

一次元の棒を伝わる縦波(圧縮波)の伝播速度は、材料の剛性(ヤング率)と密度で決まります。

$$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ソフトが内部で何をしているのか、その本質が見えてくるはずだ。