弹性波传播特性学习笔记（一）：一维梁单元的弹性波仿真

前言

        应力波理论作为固体力学中极为活跃的研究课题，长期受到学者们的广泛关注。该理论主要研究应力波在介质中传播的规律，涉及波动产生、传播、反射、折射及其与材料相互作用的动力学过程，其核心在于揭示材料在动态载荷下的力学行为与能量传递机制。应力波理论的研究重点集中在地震、高速撞击、爆炸及冲击荷载等极端动态条件下材料的动态性能与结构的动态响应。这些研究不仅推动了对材料损伤、破坏机理的认识，也为工程防护与安全设计提供了理论依据。

        在应用方面，应力波技术已深入多个重要领域：在国防与航空航天领域，它被用于评估武器效应、航天器防护及结构抗冲击能力；在矿山与交通机械中，用于分析机械动态载荷与疲劳损伤；在安全防护与地震工程中，则协助设计抗震结构及监测地壳活动。近年来，随着检测技术的发展，应力波理论的应用进一步拓展至无损检测与健康监测领域，尤其在损伤检测与缺陷表征方面表现突出——例如，通过分析应力波在材料中的传播信号，可识别混凝土结构中的裂缝、分层等内部缺陷，评估复合材料的分层损伤，以及监测压力容器、管道等设备的疲劳裂纹。此外，在超声传感器优化、声发射技术及结构健康监测系统中，应力波理论也发挥着关键作用，为实现工程结构的安全评估与寿命预测提供了重要手段。

一维波动理论

         以一维圆柱杆元为例，引入波动理论解析方法。以位移为未知函数，建立杆件微元运动方程如下：

其中 ，u 为杆件横向位移值，γ 为泊松比，rg 为截面对圆柱中心水平轴的回转半径，ρ 为杆密度，c0 为无频散的纵波波速。可以将下式谐波形式解代入微元方程：

        记位移u对时间t的二阶导为ω，且 c=ω/k 则可以得到下面新的化简后微元方程。

        此时c为纵波波速，c0为无频散纵波波速，波动理论中众所周知有解析式c0=（E/ρ）^0.5，此时当圆柱杆半径较小时，则有c≈c0，可以近似用解析式计算理论纵波波速。

一维弹性波仿真

        对基于Abaqus的一维弹性波仿真，初学者可能存在多种疑问或忽视某些注意事项，这里就个人仿真心得和研究过程，对几个关键点进行总结说明，个人能力有限可能存在疏漏，欢迎交流指正。

关键点1：

        一维梁单元在Abaqus中可以建立wire然后指派梁单元截面类型实现，也可以通过二维草图拉伸为实体梁实现，那对于弹性波仿真问题应该选哪一种？

答：看需求；如果只是为了简单验证理论、研究物理参数比如弹模或线密度对波传播特性影响的话，其实用wire方式建立的一维梁单元就可以很方便地实现，也不需要纠结实体梁可能存在的网格划分问题；但是如果需要进行梁单元耦合缺陷、弹性支撑、附加质量、预应力等进阶研究课题，还是需要用实体梁单元。并且就本人简单尝试两种单元分别验证波速解析解的实例来说，精度都比较高，都是吻合波动理论的。

关键点2-1：

        怎么在后处理输出xy图？

答：在mesh板块预定义好要输出的节点到set集后，在后处理板块按上图所示步骤即可输出xy图。

关键点2-2：

        为什么进行显式动力学分析后，输出的xy图中x轴步长很小，曲线非常不平滑，抑或是识别弹性波波速时识别解和理论值根本对不上？

答：曲线不平滑很可能是field output的interval步数设置太少，一般abaqus默认是20，可以先把他改成100不行再逐步增加即可。如果改完interval发现弹性波波速识别还是和理论值对不上，那很可能是网格划分密度不够，可以加密网格试试。

关键点2-3：

        xy图画出来了之后要识别波速，要怎么导出这个曲线图？

答：在后处理界面输出玩xy图后，点击左侧的XYDATA，右键对应要输出的图表文件，点击Edit，进入后点击X或Y即可选择整列数据，直接提取到Excel中即可完成导出。

综上，设置合理的话即可输出较为平滑的XY图如下：