弹性波传播仿真器 返回
Interactive Tool — Wave Physics

弹性波传播仿真器
(P波/S波)

实时可视化固体中纵波(P波)和横波(S波)的速度与波形。通过改变材料弹性模量和密度,体验超声波探伤和地震波分析的基础。

材料预设
物性参数
杨氏模量 E
GPa
GPa単位 — 鋼210, Al70, CFRP135, 混凝土30
泊松比 ν
密度 ρ
kg/m³
波形参数
频率 f
kHz
探伤超声波: 1〜10 MHz / 地震波: 1〜10 Hz
振幅
播放控制
波形比較
💡 画布点击: 设置波源位置。从点击位置产生波。
波形动画(实时) t = 0.00 μs
P波(纵波/压缩波) S波(横波/剪切波)
Wave
计算结果摘要
计算结果
cP [m/s]
cS [m/s]
λP [mm]
λS [mm]
Z [MRayl]

理论笔记 — 弹性波速度

固体中弹性波的速度由材料的弹性常数和密度决定。P波(Primary wave)是粒子振动方向与波传播方向相同的纵波,S波(Secondary wave)是垂直于传播方向的横波。

$$c_P = \sqrt{\frac{E(1-\nu)}{\rho(1+\nu)(1-2\nu)}}, \quad c_S = \sqrt{\frac{E}{2\rho(1+\nu)}}= \sqrt{\frac{G}{\rho}}$$

两者速度之比 $c_P/c_S = \sqrt{2(1-\nu)/(1-2\nu)}$ 仅取决于泊松比,与材料无关(钢约为1.73)。声学阻抗 $Z = \rho \cdot c_P$ 对超声波探伤中的反射率计算不可或缺。

波动方程: $\partial^2 u/\partial t^2 = c^2 \partial^2 u/\partial x^2$ 的解为 $u(x,t) = A\sin(kx - \omega t)$($k = \omega/c$)。

什么是弹性波(P波与S波)

🙋
P波和S波是什么?听起来好复杂。
🎓
简单来说,你可以把它们想象成固体材料里的两种“震动”。P波是推拉式的震动,像弹簧一样一缩一伸,粒子振动方向和波前进方向一致。S波是剪切式的震动,像抖动绳子,粒子振动方向和波前进方向垂直。在实际工程中,比如地震时,P波先到,S波后到,这就是为什么我们先感到上下颠簸,然后才是左右摇晃。
🙋
诶,真的吗?那为什么它们的速度不一样呢?
🎓
问得好!根本原因在于材料抵抗“压缩”和“剪切”变形的能力不同。P波速度主要取决于材料抵抗压缩的能力(体积模量),而S波速度取决于抵抗剪切变形的能力(剪切模量)。通常固体抵抗压缩比抵抗剪切更容易,所以P波比S波快。你可以在模拟器里试试看,把“杨氏模量”的滑块调高,你会发现两种波速都变快了,但增长的幅度不一样。
🙋
那旁边的“泊松比”滑块是干嘛的?调它好像也会影响速度。
🎓
泊松比描述材料被压缩时,横向会“胖”多少。它直接联系了材料的压缩和剪切特性。比如橡胶泊松比接近0.5(几乎不可压缩),钢铁约0.3。你试着把泊松比从0.1慢慢调到0.4,会发现P波速度在降低,而S波速度变化不大。这能帮你直观理解为什么像水这样的介质(泊松比≈0.5)只能传播P波(纵波),没有S波(横波)。

物理模型与关键公式

P波(纵波)的传播速度由材料的弹性模量和密度决定,它反映了介质对体积压缩的响应。

$$c_P = \sqrt{\frac{E(1-\nu)}{\rho(1+\nu)(1-2\nu)}}$$

其中,$c_P$是P波波速(m/s),$E$是杨氏模量(Pa),$\nu$是泊松比,$\rho$是密度(kg/m³)。分母中的$(1-2\nu)$项表明,当$\nu$接近0.5时(如橡胶、水),P波速度会变得非常大,因为材料几乎不可压缩。

S波(横波)的传播速度只取决于材料的剪切模量$G$和密度,它反映了介质对形状变化的抵抗能力。

$$c_S = \sqrt{\frac{E}{2\rho(1+\nu)}}= \sqrt{\frac{G}{\rho}}$$

其中,$c_S$是S波波速(m/s),$G$是剪切模量(Pa)。公式表明S波速度与泊松比$\nu$的关系相对简单。一个重要的工程参数是波速比$c_P/c_S$,它只与泊松比有关:$\frac{c_P}{c_S}= \sqrt{\frac{2(1-\nu)}{1-2\nu}}$,这在地震波分析中用于反推地下介质的性质。

现实世界中的应用

超声波无损探伤:在航空航天和制造业中,工程师向零件内部发射高频P波。通过分析反射回波的时间和强度,可以像“B超”一样检测出材料内部的裂纹、气孔等缺陷,而无需破坏零件本身。模拟器中的“频率”参数就对应了探伤超声波的频率。

地震波分析与地质勘探:地震发生时,全球地震台网通过记录P波和S波到达的时间差,可以精确定位震源。在石油勘探中,人工制造地震波,通过分析地下不同岩层反射和折射的P波,可以绘制出地下数千米的地质结构图,寻找油气储层。

医学超声成像:虽然人体软组织主要传播纵波(类似P波),但超声诊断设备(如B超、彩超)的核心原理同样是弹性波传播。通过计算声波在组织界面处的反射(与声阻抗$\rho c$相关),来构建人体内部的图像。

结构健康监测:在大型桥梁、大坝或风力发电机叶片上布设传感器网络,持续监测结构中传播的弹性波(声发射)。如果出现裂纹,会产生新的波源,通过分析波速和波形变化,可以实时预警结构损伤的位置和严重程度。

常见误解与注意事项

模型假设:本模拟器所用数学模型基于线性、均质、各向同性等简化假设。在将计算结果直接用于设计决策之前,务必确认实际系统是否满足这些假设。

单位与量纲:许多计算错误源于单位换算错误或数量级判断失误。请时刻注意各参数输入框旁标注的单位。

结果验证:始终将模拟器输出结果与物理直觉或手算结果进行核对。若结果出乎意料,请检查输入参数或采用独立方法进行验证。

进阶学习指引

深化理论:在本工具的简化模型基础上,进一步研究非线性效应、三维行为和时间依赖现象。阅读专业教材和学术论文,掌握严格的数学推导,是提升工程解题能力的关键。

数值方法:系统学习有限元法(FEM)、有限差分法(FDM)和有限体积法(FVM),理解商业CAE求解器的内部运行机制,这将显著提升您设置有效仿真的能力。

实验验证:理论和仿真结果必须通过实验数据加以验证。养成将计算结果与测量值进行对比的习惯,这正是V&V(验证与确认)的精髓所在。

CAE工具:准备好后,可进一步探索Ansys、Abaqus、OpenFOAM、COMSOL等业界主流工具。通过本模拟器培养的物理直觉,将帮助您更有效地配置和使用这些工具。

相关工具

截面惯性矩/惯性矩计算器
实时计算矩形、圆形、I形和T形截面的截面惯性矩和重心位置
SIMP相位最优化仿真器
2D有限元法(Q4单元)实时拓扑优化
应力集中系数 Kt
由形状参数计算孔、缺口和阶梯轴的Kt
欧拉屈曲载荷
由截面形状和端部条件(4种)计算临界载荷Pcr和屈曲模态形状