材料 A（入射側）
材料
材料 B（透過側）
材料
荷重・形状
計算結果
応力パルスの伝播・反射アニメーション
パルスは波速 $c=\sqrt{E/\rho}$ で進み、端で反射します。固定端では同符号で反射し境界で応力が2倍に重なり、自由端では符号が反転（圧縮↔引張）して反射します。鋼（E=200 GPa, ρ=7850）では c≈5050 m/s。
x-t 線図（ラグランジュ線図）
応力分布（t = L/c, 2L/c, 3L/c）
理論・主要公式
波速：$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}$
弾性波伝播シミュレーターとは
🙋
このシミュレーターで「弾性波」って何がわかるんですか？衝撃を与えたら、棒の端まで一気に変形するのではないの？
🎓
大まかに言うと、衝撃は「波」として伝わっていくんだ。例えばハンマーで鉄の棒の端を叩くと、その変形は一瞬で反対側には届かず、波速 $c$ で伝播する。このシミュレーターでは、材料AとBの境目で波がどう反射・透過するかがリアルタイムで見られるよ。まずは上のスライダーで「衝撃力 F」を変えて、波の大きさがどう変わるか確認してみて。
🙋
え、そうなんですか！「反射係数 R」がマイナスになることもあるみたいですけど、これはどういう意味ですか？
🎓
良いところに気が付いたね。反射係数がマイナスというのは、波の「位相」が反転する、つまり圧縮波が引張波として反射されることを意味するんだ。実務で多い例は、鋼（材料A）から空気（材料B）のような非常にインピーダンスの小さいものへの境界だね。シミュレーターで材料Bの密度 $\rho$ を極端に小さく設定すると、Rが-1に近づいて「自由端反射」の状態が再現できるよ。
🙋
なるほど！でも、この「x-t線図」ってどう見ればいいんですか？斜めの線がたくさんあって…。
🎓
あの図は、横軸が位置、縦軸が時間で、波面の動きを追跡した「道のり」みたいなものだ。傾きが波速の逆数になるから、材料が変わると傾きも変わるんだ。界面で線が折れ曲がるのが反射や透過の瞬間。パラメータで材料AとBのヤング率 $E$ や密度 $\rho$ を変えて、線の傾きと反射の様子がどう連動するか観察してみよう。超音波探傷の原理が直感的に理解できるはずだよ。
よくある質問
Rが負の値は、反射波の応力の符号が入射波と反転することを意味します。例えば、圧縮波が入射した場合、反射波は引張波として戻ります。これは材料Bのインピーダンスが材料Aより小さい場合に発生し、界面剥離などの評価で重要な指標となります。
波の伝播速度cは√(E/ρ)で決まるため、ヤング率と密度を同じ比率で変更すると速度は変わりません。例えば両方を2倍にすると速度は√(2/2)=1倍で不変です。速度を変えたい場合は、Eとρの比を変えてください。
理論上は音響インピーダンスZ=ρcAに断面積が含まれますが、本シミュレーターは材料A・Bの断面積を同一(等断面棒)と仮定しているため、AはZ_B/Z_A比で相殺され、反射係数R=(Z_B−Z_A)/(Z_B+Z_A)・透過係数T=2Z_B/(Z_B+Z_A)はAに依存しません。Rを変えるにはρcの異なる材料を選んでください。断面積が界面で変わる場合（段付き棒）はZにAを含めて再計算が必要です。
x-t線図の傾きは波の伝播速度を表します。界面で材料が変わると速度cが変化するため、波の軌跡の傾きが変わります。また、反射波は進行方向が逆になるため、傾きの符号が反転します。これにより折れ曲がった線が描かれます。
実世界での応用
材料試験（ホプキンソン棒試験）：高速衝撃による材料の動的応力-ひずみ関係を評価する際の基礎理論です。試験片を挟んだ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のような一見不思議な値が出ても、エネルギー計算をしてみると辻褄が合うことを確認してみて。