Stress Wave Propagation Back
Impact & Wave Mechanics

Elastic Stress Wave Propagation Simulator

Select materials A and B to see stress wave reflection and transmission at the interface on a live x-t Lagrange diagram. Computes wave speed, acoustic impedance, R and T in real time.

Material A (incident side)
Material
Material B (transmission side)
Material
Geometry & Load
Rod length L2.0 m
Impact duration τ100 μs
Impact force F100 kN
Results
Wave speed cA
m/s
Wave speed cB
m/s
Reflection R
Transmission T

Key Equations

Wave speed: $c = \sqrt{E/\rho}$ m/s
Impedance: $Z = \rho c A$
Reflection: $R = \dfrac{Z_B - Z_A}{Z_B + Z_A}$
Transmission: $T = \dfrac{2Z_B}{Z_B + Z_A}$
x-t Lagrange Diagram
Incident
Reflected
Transmitted
Stress Distribution (t = L/c, 2L/c, 3L/c)

What is Stress Wave Propagation?

🧑‍🎓
What exactly is a "stress wave"? I think of sound waves, but this simulator shows a rod getting hit. Are they related?
🎓
Basically, yes! A stress wave is a mechanical disturbance that travels through a material, just like sound. When you hit the end of this rod with the "Impact force F" slider, you create a local compression. That compression travels down the rod at a specific speed. Try setting the first material to "Aluminum" and hitting "Simulate"—you'll see the wavefront move.
🧑‍🎓
Wait, really? So the speed depends on the material? What happens when the wave reaches the other end, where the two materials meet?
🎓
Great question. The speed is $c = \sqrt{E/\rho}$, controlled by Young's modulus (E) and density (ρ). When the wave hits the interface, part of it reflects back and part transmits into the second material. This depends on their "acoustic impedance." In the simulator, try making Material A "Steel" and Material B "Rubber." You'll see a strong reflection because their impedances are very different.
🧑‍🎓
That makes sense. So the reflection and transmission coefficients (R & T) tell me how much force or stress comes back or goes through? Can they ever be negative?
🎓
Exactly. They determine the stress amplitude. R can be negative! If the wave moves from a low-impedance material to a high-impedance one (like Rubber to Steel), R is positive and the reflected wave is compressive. If it goes from high to low (Steel to Rubber), R is negative, meaning the reflected wave is tensile. Play with the E and ρ sliders to see R and T update in real-time. A common case is a crack forming if that tensile reflection is too strong.

Physical Model & Key Equations

The one-dimensional wave speed in a slender rod is governed by the material's stiffness and inertia. This is the fundamental speed at which any mechanical disturbance travels.

$$c = \sqrt{\frac{E}{\rho}}$$

$c$: Longitudinal wave speed (m/s)
$E$: Young's Modulus, material stiffness (Pa)
$\rho$: Density, material inertia (kg/m³)

When a stress wave encounters a change in material (an interface), its behavior is determined by the acoustic impedance mismatch. The coefficients define the amplitude of the reflected and transmitted stress waves relative to the incident wave.

$$R = \frac{Z_B - Z_A}{Z_B + Z_A}\quad,\quad T = \frac{2Z_B}{Z_B + Z_A} \quad,\quad Z = \rho c A$$

$R, T$: Reflection & Transmission coefficients for stress
$Z_A, Z_B$: Acoustic impedance of Material A & B (Pa·s/m)
$A$: Cross-sectional area of the rod (constant here). A higher impedance material is "stiffer" to wave motion.

Real-World Applications

Non-Destructive Testing (NDT): Engineers use stress waves to inspect structures for internal flaws without damaging them. A transducer sends a pulse into a pipeline or aircraft wing; the timing and amplitude of reflections (seen on the simulator's Lagrange diagram) reveal cracks or corrosion.

Seismic Engineering: Earthquake waves are large-scale stress waves traveling through the Earth's layers. Understanding how they reflect and transmit at soil-rock interfaces is critical for designing foundations and predicting ground motion in cities.

Ultrasonic Welding & Cleaning: High-frequency stress waves are used to fuse plastics or metals (welding) or to create cavitation for cleaning delicate parts. The efficiency depends on properly matching the tool's impedance to the workpiece, just like in the simulator.

Sports Equipment Design: The "feel" of a golf club or tennis racket hitting a ball is governed by stress wave propagation. Designers tweak materials and layers to control how much vibration (wave reflection) travels back to the player's hands, optimizing performance and comfort.

Common Misconceptions and Points to Note

When you start using this simulator, there are a few points that are easy to misunderstand. First, you might tend to think "the wave shape doesn't change." In the real world, waves attenuate and disperse due to internal friction in materials and geometric spreading (e.g., 3D spherical waves), but this 1D model ignores that. So, while a clean rectangular wave oscillates back and forth forever in the simulation, that's an idealized behavior.

Next, the realism of parameter settings. For example, setting an extremely low density to simulate "free-end reflection" is good for learning, but the range of densities and Young's moduli for real materials is limited. For instance, steel's density is about 7800 kg/m³, and aluminum is about 2700 kg/m³. If you input values equivalent to "air" for Material B (density ~1.2 kg/m³), the impedance ratio becomes extreme, resulting in R≈-1. This is a good experiment for deepening your understanding, but note that in practice, such an extreme interface difference only exists in cases like "material-vacuum."

Finally, "the relationship between the reflection coefficient R and the transmission coefficient T". From the law of energy conservation, the sum of the energies of the reflected and transmitted waves should equal the energy of the incident wave. However, the R and T displayed here are ratios of stress (or particle velocity), so they don't simply satisfy R+T=1. The energy reflection and transmission coefficients are calculated by $R^2$ and $(Z_A/Z_B)T^2$. If you see seemingly odd values like R=0.5, T=1.2 in the simulator, try checking the energy calculation to see that it balances out.

Related Engineering Fields

The fundamentals of this 1D elastic wave underlie a surprisingly wide range of fields. The first to mention is "Acoustical Engineering and Ultrasonic Machining." In ultrasonic cleaners and ultrasonic welders, elastic waves generated from a transducer (piezoelectric element) must be efficiently transmitted to the horn (tool). At junctions between different materials, the impedance matching you learned in this simulator becomes a key design factor. If the impedance doesn't match, energy is reflected, preventing sufficient amplitude from reaching the processing area.

Another is "Geotechnical Engineering and Seismic Design." Seismic S-waves and P-waves are essentially elastic waves propagating through the ground. When a wave travels from hard bedrock to soft alluvial soil, reflection and transmission occur due to impedance differences. If the transmission coefficient is large (= much energy reaches the surface), the ground motion is amplified, increasing the impact on buildings. You can simply model this situation in the simulator by setting Material B to be much softer (lower E) than Material A.

Furthermore, applications in "Biomedical Engineering" are also interesting. In ultrasonic diagnostic devices (echo), an ultrasonic pulse is sent from the probe into the body. The acoustic impedance differs for each tissue: skin, fat, muscle, bone, etc. At boundaries with extremely high-impedance tissues like bone, strong reflections (echoes) return. This 1D model is an excellent first step to understanding this principle.

For Further Learning

Once you're comfortable with this simulator, as a next step, try exploring "the mathematics of the wave equation." Behind this simulation lies the one-dimensional wave equation $\frac{\partial^2 u}{\partial t^2} = c^2 \frac{\partial^2 u}{\partial x^2}$, where $u$ is displacement. Solving this partial differential equation using a method called "d'Alembert's solution" yields a solution in the form of "traveling waves": $u(x,t) = f(x-ct) + g(x+ct)$. The rightward and leftward waves you see in the simulator correspond precisely to this $f$ and $g$.

To get closer to practical applications, understanding "extension to 3D and wave modes" is important. In actual structures, not just longitudinal waves along the axis of a rod (the waves handled in this simulator), but various modes like bending waves (transverse waves) and torsional waves can occur simultaneously and sometimes couple with each other. For example, the analysis of "Lamb waves" generated in plates upon impact or "guided waves" propagating through pipes is utilized at the forefront of non-destructive testing.

As a concrete next learning topic, I recommend "basics of numerical methods (FDM, FEM)." This simulator likely calculates using a method close to an analytical solution, but for complex shapes or nonlinear materials, numerical simulations like the Finite Element Method (FEM) are essential. Start by programming a simple Finite Difference Method (FDM) solution for the 1D wave equation yourself. This will help you physically grasp fundamental concepts of numerical analysis like the importance of time steps and spatial discretization, and stability conditions (e.g., the Courant condition $\Delta t \le \Delta x / c$). From there, you should start to see the essence of what commercial CAE software is doing internally.