弹性应力波传播 返回
冲击 · 波动力学

弹性应力波传播模拟器

选择材料A和B,实时在x-t拉格朗日线图上显示应力波的反射与透射。基于声阻抗比计算反射系数R和透射系数T。

材料 A(入射侧)
材料
材料 B(透射侧)
材料
几何与荷载
杆长 L
m
冲击持续时间 τ
μs
冲击力 F
kN
计算结果
计算结果
波速 cA
m/s
波速 cB
m/s
反射系数 R
透射系数 T
应力脉冲传播与反射动画
波速 c [m/s]
脉冲位置 x [m]
峰值应力 [×σ₀]
通过时间 L/c [μs]
入射脉冲
反射脉冲
叠加(合成)
脉冲以 $c=\sqrt{E/\rho}$ 传播并在端部反射。固定端以相同符号反射,边界处应力加倍;自由端反射时符号反转(压缩↔拉伸)。钢(E=200 GPa, ρ=7850)的 c≈5050 m/s。
x-t 拉格朗日线图
入射波
反射波
透射波
应力波传播 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}$

什么是弹性应力波传播

🙋
“弹性应力波”是什么?听起来好抽象啊。
🎓
简单来说,就像你用手快速拍一下一根很长的铁杆一端,那个“拍”的力并不会瞬间传到另一端,而是像水波一样,以压缩变形的形式“跑”过去。这个“跑”的压缩变形就是弹性应力波。在实际工程中,比如用锤子敲打桩基,这个冲击力就是以应力波的形式传到地下的。
🙋
诶,真的吗?那如果这根杆子中间材料不一样了会怎样?
🎓
问得好!这正是我们这个模拟器的核心。当波跑到两种材料的交界处时,一部分会像撞到墙一样“弹回来”,这叫反射;另一部分会“穿过去”继续跑,这叫透射。你可以在模拟器里试试,把材料A从钢换成铝,把材料B保持为钢,然后点击运行,你会看到波在界面处分裂成两股的动态过程。
🙋
哦!那“弹回来”多少,“穿过去”多少,是怎么决定的呢?
🎓
这取决于两种材料的“声阻抗”,你可以把它理解为材料对波通过的“阻力”。模拟器右上角会根据你选择的材料自动计算反射系数R和透射系数T。试着把材料B的密度ρ调得特别大,你会发现R接近1,意味着波几乎全被反射回来了,就像撞到了固定端。

物理模型与关键公式

波在均匀杆中的传播速度由材料的弹性模量和密度决定,这是波传播的基本特性。

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

其中,$c$是波速(m/s),$E$是材料的弹性模量(Pa),$\rho$是材料密度(kg/m³)。波速越快,冲击荷载传播得越快。

当波传播到两种不同材料的界面时,其反射和透射行为由声阻抗决定。声阻抗是材料抵抗波传播能力的度量。

$$Z = \rho c A = \sqrt{E \rho}A$$

其中,$Z$是声阻抗,$A$是杆的横截面积。反射系数$R$和透射系数$T$由界面两侧材料的声阻抗比$Z_A$和$Z_B$决定:

$$R = \frac{Z_B - Z_A}{Z_B + Z_A}, \quad T = \frac{2Z_B}{Z_B + Z_A}$$

$R$和$T$分别表示反射波和透射波的应力幅值与入射波应力幅值的比值。注意,$R$可以是负数,表示反射波相位反转(如自由端反射)。

现实世界中的应用

桩基检测:工程现场常见的是“高应变动力试桩”。用重锤敲击桩顶,产生的应力波沿桩身向下传播,遇到桩底或缺陷(如缩颈、断桩)时会发生反射。通过分析反射回来的波信号,就能判断桩的完整性和承载力,而无需把桩挖出来。

霍普金森杆实验:这是测量材料在高应变率下力学性能的核心实验。子弹撞击入射杆产生应力波,波通过试件传到透射杆。通过测量杆上的应变,就能反推出试件在高速冲击下的应力-应变曲线,广泛应用于军工和航空航天材料测试。

超声波无损检测:比如检查飞机蒙皮、铁轨或大型铸件内部是否有裂纹。探头发出高频应力波(超声波),波在材料内部传播,遇到裂纹(即阻抗突变界面)会反射。通过接收反射波的时间和强度,就能精确定位缺陷的大小和位置。

汽车碰撞安全分析:在汽车碰撞试验的CAE模拟中,冲击载荷在车架结构中的传播本质就是复杂的应力波传播过程。工程师需要分析波在钣金件、焊点、不同材料连接处的反射与透射,以优化结构吸能和乘员保护区的刚度设计。

常见误解与注意事项

开始使用本模拟器时,有几个容易误解的地方。首先,人们常误以为“波形保持不变”。在现实世界中,波会因材料内摩擦和几何扩散(例如三维球面波)而衰减和色散,但在这个一维模型中这些因素被忽略了。因此,虽然模拟中完美的矩形波会永远往复传播,但这只是一种理想化的行为。

其次,参数设置的现实性。例如,将密度设置得极低以模拟“自由端反射”虽有助于学习,但实际材料的密度和杨氏模量范围是有限的。比如钢的密度约为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这类看似奇特的数值,不妨验证一下能量计算,结果会自洽。