🔬 physics

Numerical Modeling of Surface Waves in the Block-Structured Soil Massifs

本文通过使用 3Dyn_Granular 软件进行三维数值研究，证明了具有顺应性及塑性变形特征的周期性土块系统能够有效地衰减表面波，并保护工程结构免受破坏性冲击。

原作者： Vladimir M. Sadovskii, Ilia S. Nikitin, Oxana V. Sadovskaya

发布于 2026-08-06

📖 1 分钟阅读☕ 轻松阅读

原作者： Vladimir M. Sadovskii, Ilia S. Nikitin, Oxana V. Sadovskaya

原始论文采用 CC BY 4.0 许可（https://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

大地的无形震颤与伟大的土壤护盾

想象一下，你脚下的地面不仅仅是坚硬、不变的岩石，而是一个巨大的、有生命的鼓。当地震发生或大规模爆炸发生时，它不仅仅是在大地中发出一次“砰”的撞击，而是会发出涟漪，就像把石头丢进池塘一样。虽然其中一些涟漪会深入地下并迅速消散，但最危险的是“表面波”。这些涟漪沿着地表的最顶端传播，像冲浪者乘着浪潮一样紧贴着表面移动。因为它们停留在表面，所以能量损失不像深层波那样快。当它们到达城市、桥梁或摩天大楼时，依然带着足以摧毁地基和折断管道的能量在剧烈摇晃。

科学家们早已知道，要阻止这些破坏性的波，你不能仅仅建造一堵墙；你必须诱骗这些波失去能量。一个聪明的想法是在波的路径中埋入一系列由松散颗粒材料（如沙子或粘土）制成的“软块”。这就像是尝试让一辆疾驰的汽车驶过一片深厚、蓬松的雪地，而不是撞向一面坚实的混凝土墙。雪并不仅仅是将汽车弹回，它会抓住汽车的能量，被挤压下去，并将这种运动转化为热量，从而显著减慢汽车的速度。本文深入探讨了数学模型和计算机模拟，以观察这种“软块”策略是否真的对地球表面波有效，以及它究竟能防止多少破坏。

伟大的波浪赛跑：松散的土块能挡住震动吗？

在这项研究中，来自俄罗斯的一个研究小组决定利用一台超级强大的计算机，进行一场高风险的“波浪捉迷藏”游戏。他们想看看当一个巨大的表面波在坚固的弹性景观上疾驰时，撞上一排由松散、柔软的土壤块组成的特殊障碍物时会发生什么。

设置：数字沙盒
研究人员在他们的计算机中构建了一个巨大的虚拟世界。想象一个巨大的矩形土壤箱，长45米，宽30米，深15米。在左侧，他们生成了一个完美的、强大的表面波，模拟自然地震或人造爆炸释放的那种能量。这种波被设计为沿表面传播，直奔一系列障碍物而去。

这些障碍物就是“屏障块”。它们不是坚固的混凝土；它们是由一种具有顺应性、颗粒状的材料（沙子和粘土的混合物）制成的，这种材料拥有一个秘密超能力：它可以改变其行为。当波撞击它们时，这些块体可以挤压、拉伸，甚至发生“塑性变形”。用通俗的话说，这意味着这些块体可以永久改变形状来吸收能量，而不是像橡胶球一样把波弹回去。这些块体呈周期性排列——就像埋在地下的一排栅栏——每个块体厚3米，间隔6米。

模拟：观看波浪撞击
使用名为 3Dyn Granular 的定制软件程序，团队在一个拥有96个处理器的超级计算机集群上运行了一个复杂的模拟。他们观察着波浪在数字景观上竞逐的过程。

首先，他们测试了一个“均匀”的情景，即整个地面只是坚固的弹性土壤。在这种情况下，波在传播过程中几乎没有能量损失，就像在完美且无尽的大海中冲浪的人一样。

然后，他们引入了障碍物。他们运行了两个不同版本的实验：

“弹性”块： 这些块体由相同的沙粘土混合物组成，但表现得像完美的弹簧。它们会挤压并弹回，而不会发生永久性的形状改变。
“弹塑性”块： 块体材料相同，但这次允许它们发生永久变形。当波撞击时，这些块体会被挤压，产生热量并吸收能量，就像“蓬松的雪”那个类比一样。

结果：软块体的胜利
结果清晰且令人兴奋。当波撞击仅含弹性的块体时，应力波的振幅（高度或强度）下降了约 20%。这是一个好的开始，就像在路上放了一个减速带；它减慢了速度，但汽车仍然跑得很快。

然而，当研究人员切换到弹塑性块体时，游戏规则彻底改变了。应力波的振幅下降了超过 40%。为什么呢？因为这些块体不仅是弹开波，它们还“吃掉”了能量。随着波迫使这些松散的土壤块压缩和变形，材料将这种剧烈的摇晃运动转化为了热量和内部摩擦。这就像波穿过一片深厚、粘稠的泥沼，泥沼抓住了它的腿并将其拖慢。

意外的转折
不过，发生了一件奇怪的事情。虽然波在表面及其传播方向上的强度显著下降，但在质量块底部，向下挤压地面的压力（称为压应力， $\sigma_{11}$ ）实际上增加了。在弹性块体的情景中，这种压力增加了一倍多！看来，虽然这些块体擅长阻止摇晃运动，但它们也将部分能量重新定向到了更深的地下。

结论
研究得出结论，埋入一排松散的塑性土壤块是减缓表面波的一种根本有效的手段。模拟显示，通过选择合适的材料和排列方式，你可以削减表面波近一半的破坏力。研究人员建议，使用更具孔隙率的材料（具有更多气孔的材料）可以使这种效果更强，因为从松散土壤到紧实土壤的转变会产生消耗更多能量的冲击波。

虽然这只是计算机模拟而非现实世界的地震测试，但数学逻辑是成立的。论文指出，如果工程师能够正确建造这些“软屏障”，他们就可以保护桥梁、管道和建筑群免受最严重的震动影响，将破坏性的地震波转化为一种温和得多的隆隆声。这是一个充满前景的想法，它将地面本身变成了一个巨大的能量吸收盾牌。

技术摘要：块状结构土体中表面波的数值模拟

问题陈述
本文探讨了保护延展性基础设施（如桥梁、管道和跑道）免受破坏性地震表面波（特别是瑞利波和勒夫波）影响的挑战。与体波相比，表面波沿弹性半空间的自由边界传播，其衰减速度显著较慢，因此对近地表结构构成了严重威胁。本研究重点关注一种垂直防护地震屏障系统的数学建模与数值模拟。该系统由嵌入均匀弹性土体中的周期性排列块体组成，这些块体由顺应性颗粒材料（具体为松散的砂粘混合物）制成。研究目标是量化表面波振幅在与此类块状结构屏障相互作用时的衰减情况，并评估不同材料流变特性（弹性与弹塑性）在阻尼这些波方面的有效性。

研究方法
研究采用了一种结合解析初始数据与复杂屏障材料流变模型的三维数值方法。

初始条件与边界条件： 问题是在空间笛卡尔坐标系中构建的。入射表面波是在均匀弹性介质中使用 Shemyakin 的瑞利波传播广义解生成的。该解利用复变量的任意解析函数来定义具有特定振幅峰值的局部化波形，而非纯周期性波。定义域左侧的边界条件通过匹配该精确解来设定，而远端和底部边界则采用弱反射条件，以模拟无限质量体。
材料建模：
- 弹性介质： 周围土壤被建模为均匀、各向同性的弹性介质。
- 颗粒屏障： 屏障块采用预松散颗粒介质的流变模型进行描述。该模型考虑了以下因素：
  1. 弹塑性行为： 对拉伸和压缩具有不同的抗力，通过具有特定屈服点 ( $\sigma_s^+$ 和 $\sigma_s^-$ ) 的塑性铰进行建模。
  2. 孔隙率与压实： 一个带有初始间隙 ( $\varepsilon_0$ ) 的刚性接触单元代表初始孔隙率。当间隙闭合时发生压实，导致从低刚度状态向高模量硬化状态的转变。
  3. 变分不等式： 针对多孔介质的本构关系通过变分不等式进行公式化。这种基于 Prandtl-Reuss 塑性理论及特定压实锥的方法，能够对弹性、塑性和压实状态之间的转变进行严谨的数学处理。
计算算法： 作者使用了专有的高性能计算软件程序包 3Dyn Granular。该算法采用带有物理分裂（physical splitting）的时间显式法求解方程组和变分不等式。
- 弹性步骤在空间变量上使用两循环分裂法进行求解，并采用显式单调 ENO（本质无振荡）预测-校正格式。
- 塑性步骤涉及通过对每个网格单元的变分不等式求精确解，来修正弹性解以满足屈服面约束。
- 并行化通过 MPI 库在集群架构上实现，将空间域分配给不同的处理器。

关键结果
数值实验在一个包含周期性屏障块系统（厚度 3 m，间距 6 m）的土体质量体（深 15 m，宽 30 m，长 45 m）上进行。研究对比了三种情景：均匀弹性质量体、具有纯弹性屏障块的质量体，以及具有弹塑性屏障块的质量体。

应力衰减：
- 在均匀弹性质量体中，由于数值耗散和边界效应，波在传播距离内自然衰减了约 12–20%。
- 在弹性屏障情景下，正应力 ( $\sigma_{22}$ ) 和切应力 ( $\sigma_{12}$ ) 的振幅比均匀情况降低了约 20%。
- 在弹塑性屏障情景下，应力振幅比均匀情况降低了超过 40%。对一维模型的分析进一步表明，利用高孔隙率材料可能会通过冲击波耗散机制带来更大的衰减，尽管这一具体的增强效果在三维数值结果中未被完全量化。
应力重分布： 屏障的存在引起了应力场的重分布。值得注意的是，在屏障后方的深处，压应力 ( $\sigma_{11}$ ) 显著增加（在弹性情况下增加了一倍以上），这表明存在复杂的相互作用，即波绕过块体，使应力集中在特定区域。
能量耗散机制： 结果证实，波振幅的降低是由三个因素驱动的：介质界面处的反射、颗粒材料内的塑性能量耗散，以及与松散材料压实相关的冲击波耗散。

意义与主张
本文声称证明了使用填充顺应性颗粒材料的块状土壤屏障来有效阻尼表面波的基本可能性。作者断言，其数值结果支持了利用此类屏障保护施工现场和工程结构免受破坏性地震影响的潜力。

具体而言，研究强调了：

塑性至关重要： 引入塑性变形和颗粒材料的特定流变特性（包括拉压不对称性和孔隙率）与纯弹性屏障相比，显著增强了阻尼能力。
可扩展性： 所开发的计算算法和软件程序包能够处理具有曲线界面的块状介质中的三维弹塑性问题，并可在高性能计算集群上运行。
未来方向： 虽然目前的工作确立了在三维环境下具有弹塑性块体的屏障概念的可行性，但作者指出，需要进一步的研究来完整模拟非一维设置下的材料孔隙率和冲击波转变效应，以实现更大的衰减。

这项工作通过提供用于模拟颗粒屏障的严谨数学框架（变分不等式），并利用高保真数值模拟证明了当考虑塑性耗散时，此类屏障可以降低超过 40% 的表面波振幅，从而为该领域做出了贡献。

大地的无形震颤与伟大的土壤护盾

伟大的波浪赛跑：松散的土块能挡住震动吗？

技术摘要：块状结构土体中表面波的数值模拟

类似论文