单相固体双负弹性波超材料：从设计原理到实验验证

单相固体双负弹性波超材料：从设计原理到实验验证一、引言1.1研究背景与意义超材料这一概念自20世纪60年代被提出以来，经历了从理论构想到实验实现，再到广泛应用探索的不平凡发展历程。1967年，Veslago首次系统研究了一类同时具有负介电系数和负磁导率的材料，预测了该类材料能产生负折射率等一系列新奇特性，但在当时，这类违背常规认知的假想材料因无法被设计和制备出来，在后续三十多年间未得到广泛关注，处于超材料发展的孕育阶段。直到1996年，英国帝国理工学院的Pendry教授通过电共振在共振频率附近实现了等效负介电常量材料，并于随后通过磁共振实现了等效负磁导率材料。在此基础上，Veslago预言的电磁负折射现象很快得到实验证实，正式拉开了超材料发展的序幕。由于其独特性质主要源于材料微结构对能带的调控，类似现象迅速被推广到基于相同机理的声波和弹性波领域。此后，波动控制功能设计的变换方法的提出，更是极大地推动了超材料的发展，众多新奇的波动功能器件相继被提出。自2000年后的二十年，超材料进入快速发展期，其概念和内涵不断拓展，涵盖从波动载荷到静态载荷、从局域到全局、从结构到机构、从保守到主动等多个方面，并被泛化为主要通过微结构而非化学组分实现对材料宏观行为的调节。当下，超材料研究正处于快速扩展期，可动性（机构）被引入力学超材料系统设计中，突破了传统连续介质的描述框架，对力学基础理论提出了新的挑战。在超材料的众多分支中，弹性波超材料由于能够对弹性波的传播进行有效调控，在声学和力学领域展现出了巨大的应用潜力。弹性波作为一种机械波，其传播特性与介质的弹性模量、密度等物理参数密切相关。传统材料在弹性波调控方面存在一定的局限性，而弹性波超材料通过精心设计的微观结构，打破了传统材料的限制，能够实现对弹性波传播速度、方向和振幅等的精确调控。例如，通过合理设计超材料的微结构，可以使弹性波在其中发生负折射现象，这在传统材料中是难以实现的。这种负折射特性为新型声学器件的设计提供了全新的思路，有望应用于高精度的声学成像、声波聚焦等领域。单相固体双负弹性波超材料作为弹性波超材料中的一种特殊类型，同时具备负等效质量密度和负等效弹性模量。这种独特的性质使其在声学和力学领域具有更为重要的应用价值和研究意义。在声学领域，它可以用于设计高性能的声学隐身材料。当弹性波遇到单相固体双负弹性波超材料时，由于其特殊的等效参数，弹性波能够按照特定的方式传播，绕过被隐身的物体，从而实现对物体的声学隐身效果，这在军事和保密通信等领域有着重要的应用前景；也可以用于制造高灵敏度的声学传感器。利用其对弹性波的特殊调控能力，能够增强传感器对微弱声波信号的响应，提高传感器的检测精度和分辨率，为环境监测、生物医学检测等领域提供更先进的检测手段。在力学领域，单相固体双负弹性波超材料可以应用于减振降噪。将其设计成特定的结构，能够有效地抑制弹性波的传播，减少振动和噪声的产生，这对于提高机械设备的运行稳定性、降低工业噪声污染等具有重要意义；还可以用于构建新型的力学结构，如超材料梁、超材料板等，这些结构在承受载荷时能够展现出独特的力学性能，为航空航天、汽车制造等领域的轻量化设计提供了新的解决方案。综上所述，对单相固体双负弹性波超材料的设计与实验研究，不仅有助于深入理解弹性波在复杂介质中的传播规律，推动超材料理论的进一步发展，而且在实际应用中具有重要的价值，能够为声学和力学相关领域的技术创新提供有力的支持，具有广阔的研究前景和应用空间。1.2国内外研究现状在单相固体双负弹性波超材料的设计与实验研究领域，国内外学者已取得了一系列重要成果，研究内容涵盖了从理论设计到实验验证的多个方面，推动了该领域的不断发展。国外方面，早期对双负弹性波超材料的理论探索为后续研究奠定了基础。2000年左右，一些研究团队开始利用理论模型和数值模拟方法，研究如何通过微观结构设计实现负等效质量密度和负等效弹性模量。例如，[具体文献1]中提出了一种基于局域共振原理的微观结构设计方案，通过合理选择结构参数，在理论上实现了双负特性，为后续的实验研究提供了重要的理论指导。随着研究的深入，实验验证成为关键环节。一些研究团队利用先进的材料制备技术，成功制备出具有双负特性的弹性波超材料试件，并通过实验测量验证了理论预测的正确性。[具体文献2]通过微机电系统（MEMS）加工技术，制备出了亚波长尺度的弹性波超材料结构，实验测量结果表明，该结构在特定频率范围内表现出明显的双负特性，进一步证实了双负弹性波超材料的可行性。在应用探索方面，国外学者也取得了显著进展。在声学隐身领域，[具体文献3]设计了一种基于双负弹性波超材料的声学隐身结构，通过实验测试，该结构能够有效地减少目标物体在特定频率范围内的声学散射，实现了较好的隐身效果；在声学成像领域，[具体文献4]利用双负弹性波超材料的负折射特性，设计了新型的声学成像系统，实验结果显示，该系统能够突破传统声学成像的分辨率限制，提高成像质量。国内在单相固体双负弹性波超材料的研究起步相对较晚，但发展迅速，在理论创新、实验技术和应用拓展等方面都取得了丰硕的成果。在理论研究方面，国内学者提出了多种新颖的微观结构设计思路。[具体文献5]提出了一种基于多尺度结构设计的方法，通过在不同尺度上对材料结构进行优化，实现了更宽频带的双负特性，拓展了双负弹性波超材料的应用频率范围；[具体文献6]利用拓扑优化方法，对弹性波超材料的微观结构进行优化设计，提高了材料的双负性能和力学性能。在实验研究方面，国内研究团队不断改进实验技术和测试方法，提高了实验的准确性和可靠性。[具体文献7]采用3D打印技术，制备出复杂的弹性波超材料结构，并利用激光超声技术对其弹性波传播特性进行了精确测量，为材料性能的评估提供了更准确的数据支持；[具体文献8]搭建了一套基于超声相控阵的实验平台，能够对双负弹性波超材料的负折射特性进行多角度、高精度的测量，为深入研究其波动特性提供了有力的实验手段。在应用研究方面，国内学者积极探索双负弹性波超材料在多个领域的应用潜力。在减振降噪领域，[具体文献9]将双负弹性波超材料应用于汽车发动机的减振系统中，实验结果表明，该材料能够有效地降低发动机的振动和噪声，提高汽车的乘坐舒适性；在无损检测领域，[具体文献10]利用双负弹性波超材料的特殊性质，设计了新型的无损检测探头，能够提高对材料内部缺陷的检测灵敏度和准确性。尽管国内外在单相固体双负弹性波超材料的设计与实验研究方面取得了诸多成果，但仍存在一些不足之处。在材料设计方面，目前的设计方法大多依赖于经验和试错，缺乏系统的、普适的设计理论，难以快速、准确地设计出满足特定需求的双负弹性波超材料；在材料制备方面，现有制备技术在制备复杂结构和大规模制备时仍面临挑战，导致材料的制备成本较高、生产效率较低，限制了其大规模应用；在实验研究方面，实验测试技术还不够完善，对于一些微观结构和复杂波动现象的测量还存在一定的误差和局限性，影响了对材料性能的深入理解；在应用研究方面，虽然在一些领域取得了初步应用成果，但距离实际工程应用还有一定的差距，需要进一步解决材料与实际系统的兼容性、可靠性等问题。1.3研究目标与内容本研究旨在深入探索单相固体双负弹性波超材料的设计方法与实验验证，通过系统性的研究，期望实现以下目标：设计出具有稳定双负特性且性能优良的单相固体双负弹性波超材料结构，明确其关键设计参数与特性之间的关系；通过数值模拟和实验测试，全面准确地表征所设计超材料的弹性波传播特性，验证其双负特性的有效性；将所设计的超材料应用于特定的声学或力学场景，展示其在实际应用中的优势和潜力，为其进一步的工程应用提供理论和实验依据。围绕上述研究目标，本研究将开展以下几方面的具体研究内容：单相固体双负弹性波超材料的设计：基于局域共振原理，结合现有文献中关于微观结构设计的思路，创新性地设计一种新型的单相固体双负弹性波超材料的微观结构。深入研究结构参数，如单元尺寸、形状、排列方式以及材料属性等对等效质量密度和等效弹性模量的影响规律。运用参数化建模方法，建立超材料结构参数与等效参数之间的定量关系模型，为后续的结构优化提供理论基础。例如，通过改变结构中谐振单元的尺寸和形状，研究其对共振频率的影响，进而分析共振频率与等效质量密度和等效弹性模量的关联。超材料弹性波传播特性的数值模拟：利用有限元分析软件COMSOLMultiphysics，建立所设计超材料的数值模型。通过模拟弹性波在超材料中的传播过程，计算其能带结构、等频率曲线、传输谱等关键特性参数。分析这些参数在不同频率范围内的变化规律，进一步验证所设计超材料的双负特性。例如，通过能带结构分析，确定超材料中弹性波的传播禁带和通带，观察在双负特性频率范围内能带的特殊变化；对等频率曲线进行分析，研究弹性波在超材料中的传播方向和相速度等特性。超材料试件的制备与实验测试：根据数值模拟优化后的结构参数，选择合适的材料和制备工艺，制备单相固体双负弹性波超材料试件。搭建高精度的实验测试平台，采用激光超声技术、超声相控阵技术等先进的测试手段，对超材料试件的弹性波传播特性进行实验测量。将实验测量结果与数值模拟结果进行对比分析，验证数值模拟的准确性和可靠性，同时进一步深入研究超材料的实际性能。例如，利用激光超声技术测量弹性波在超材料中的传播速度和衰减特性，与数值模拟结果进行对比，分析两者之间的差异及原因。超材料在声学和力学领域的应用探索：将所制备的超材料应用于声学隐身或减振降噪等实际场景中，通过数值模拟和实验验证，评估其在实际应用中的性能表现。例如，设计基于该超材料的声学隐身结构，模拟在不同声波频率和入射角度下，目标物体的声学散射情况，评估隐身效果；将超材料应用于机械结构的减振系统，通过实验测试，分析其对振动的抑制效果和对结构力学性能的影响。根据应用效果反馈，进一步优化超材料的设计和性能，为其实际工程应用提供技术支持。二、单相固体双负弹性波超材料的设计原理2.1超材料基本概念与分类超材料（Metamaterial），从词源上看，“Meta”在希腊语中有“超越、另类”之意，这也精准地体现了超材料的本质特征。超材料是一类具有特殊性质的人造材料，这些性质是自然界中传统材料所不具备的。其成分通常并无特殊之处，然而其奇特性质主要源于精心设计的精密几何结构以及特定的尺寸大小，其中的微结构尺寸尺度小于其作用的波长，从而能够对波施加独特的影响。例如，在电磁超材料中，通过设计金属结构的形状、尺寸和排列方式，可以实现对电磁波的特殊调控，如负折射率特性，使电磁波在其中的传播方向与在传统材料中的传播方向相反，这在传统材料中是无法实现的。超材料具有一些显著的基本特征。首先是奇异物理特性，超材料能够呈现出天然介质材料所不具备的超常物理性能，如负磁导率、负介电常数、负折射率等。这些特性突破了传统材料的物理限制，为实现新奇的物理现象和应用提供了可能；其次是亚波长结构，超材料属于亚波长结构，其单元结构尺寸远远小于工作电磁波波长，在这种情况下，单元微结构可以忽略不计，从宏观上看，超材料可以被视为一种等效介质；最后是等效介质特性，由于超材料具有亚波长结构，其物理性质和材料参数可以用等效媒质理论进行描述，通过等效媒质理论，可以将超材料复杂的微观结构等效为具有特定宏观参数的均匀介质，从而便于对其电磁、声学等性能进行分析和设计。根据所调控的波的类型以及应用领域的不同，超材料可以分为多种类型，其中弹性波超材料是重要的一类。弹性波超材料是指能够对弹性波的传播进行有效调控的超材料，其特殊的物理性质主要源于精心设计的微观结构。当弹性波在弹性波超材料中传播时，会与这些微观结构发生相互作用，从而导致弹性波的传播特性发生改变。例如，弹性波超材料可以实现对弹性波的频率带隙调控，使得在某些频率范围内，弹性波无法在材料中传播，形成禁带；也可以实现对弹性波的负折射、聚焦等特殊传播现象，这些特性在声学和力学领域具有重要的应用价值。弹性波超材料根据其实现特殊性能的机理，又可以进一步细分为多种类型。基于布拉格散射机理的弹性波超材料，其带隙一般出现在波长与晶格常数相同数量级的频率范围内。当弹性波的波长与超材料的晶格常数相近时，弹性波在超材料中传播会发生布拉格散射，导致在某些频率范围内形成禁带，阻止弹性波的传播。基于局域共振机理的弹性波超材料，通过引入具有特定共振频率的局域谐振单元，当弹性波的频率接近谐振单元的共振频率时，谐振单元会发生强烈的共振，与弹性波产生相互作用，从而改变弹性波的传播特性，这种超材料更容易产生低频带隙，在低频隔振与降噪等领域具有重要应用；还有基于变换弹性理论的弹性波超材料，利用坐标变换的方法，通过设计材料的非均匀分布，实现对弹性波传播路径的精确控制，如设计声学隐身斗篷等。这些不同类型的弹性波超材料各自具有独特的性能和应用场景，为弹性波的调控提供了多样化的手段。2.2双负弹性波超材料特性及原理单相固体双负弹性波超材料同时具备负等效质量密度和负等效弹性模量这两个独特的性质，这使其在弹性波调控领域展现出与传统材料截然不同的特性。在传统材料中，质量密度和弹性模量均为正值，弹性波在其中的传播遵循常规的波动规律。而双负弹性波超材料打破了这种常规认知，其特殊的等效参数使得弹性波的传播行为发生了显著变化。负等效质量密度的产生机制主要基于局域共振原理。在双负弹性波超材料的微观结构中，通常包含一些具有特定共振频率的局域谐振单元。当弹性波的频率接近这些谐振单元的共振频率时，谐振单元会发生强烈的共振。在共振过程中，谐振单元的运动表现出与常规质量运动不同的特性，从宏观等效的角度来看，就产生了负等效质量密度的效果。以一种基于弹簧-质量块模型的双负弹性波超材料微观结构为例，质量块通过弹簧与基体相连，当弹性波作用时，在特定频率下，质量块的振动方向与弹性波的传播方向相反，其惯性表现出与常规质量相反的效果，从而导致材料在该频率范围内呈现出负等效质量密度。这种负等效质量密度特性使得弹性波在超材料中的传播速度和相速度等特性发生改变，为弹性波的调控提供了新的手段。负等效弹性模量的产生同样与超材料的微观结构和局域共振密切相关。在超材料的微观结构中，一些结构单元在弹性波作用下会发生特殊的变形模式。当弹性波的频率满足一定条件时，这些结构单元的变形会导致材料整体的弹性回复力表现出与常规弹性模量相反的特性。例如，在某些具有特殊几何形状的微观结构中，当受到弹性波的作用时，结构单元之间的相互作用使得材料在受力方向上的弹性回复力减小，甚至表现为与受力方向相反的力，从宏观上看，就相当于材料具有了负等效弹性模量。这种负等效弹性模量特性进一步丰富了双负弹性波超材料对弹性波的调控能力，使得弹性波在超材料中的传播路径、波的散射等特性都与传统材料中的情况不同。双负特性实现的理论基础涉及多个方面。从波动方程的角度来看，弹性波在双负弹性波超材料中的传播满足修正后的波动方程。由于超材料的负等效质量密度和负等效弹性模量，波动方程中的相关参数发生了改变，这直接影响了弹性波的传播特性。在常规材料中，弹性波的波动方程为[常规波动方程表达式]，而在双负弹性波超材料中，由于等效参数的变化，波动方程变为[双负材料波动方程表达式]，通过对该方程的求解和分析，可以深入理解弹性波在双负超材料中的传播规律。等效介质理论在双负弹性波超材料中也起着重要的作用。由于双负弹性波超材料的微观结构尺寸远小于弹性波的波长，根据等效介质理论，可以将超材料看作是一种具有等效质量密度和等效弹性模量的均匀介质。通过等效介质理论，可以将超材料复杂的微观结构对弹性波的作用等效为宏观的等效参数对弹性波的作用，从而简化对超材料中弹性波传播特性的分析和研究。能带理论也为双负特性的理解提供了重要的理论支持。通过对双负弹性波超材料的能带结构进行分析，可以确定材料中弹性波的传播禁带和通带。在双负特性频率范围内，能带结构会呈现出特殊的变化，这与负等效质量密度和负等效弹性模量的产生密切相关，进一步揭示了双负特性的物理本质。2.3基于特定理论的设计方法基于局域共振理论的设计方法是单相固体双负弹性波超材料设计中常用的一种手段。局域共振理论的核心在于通过引入具有特定共振频率的局域谐振单元，实现对弹性波传播特性的调控，进而产生负等效质量密度和负等效弹性模量。在实际设计过程中，弹簧-质量块模型是一种常用的简化模型。以一个简单的弹簧-质量块结构为例，质量块通过弹簧与基体相连，当弹性波作用于该结构时，在特定频率下，质量块会发生强烈的共振。当弹性波的频率接近质量块的共振频率时，质量块的振动幅度急剧增大，其惯性表现出与常规质量不同的特性，从宏观等效的角度来看，就产生了负等效质量密度的效果。通过合理调整弹簧的刚度和质量块的质量，可以精确控制共振频率，从而在所需的频率范围内实现负等效质量密度。例如，当需要在低频段实现双负特性时，可以适当增大质量块的质量和减小弹簧的刚度，降低共振频率，使负等效质量密度出现在低频区域。为了实现负等效弹性模量，在微观结构设计上可以采用一些特殊的几何形状和连接方式。一种具有特殊铰链结构的微观单元，通过铰链的特殊变形模式，在弹性波作用下，结构单元之间的相互作用使得材料在受力方向上的弹性回复力减小，甚至表现为与受力方向相反的力，从宏观上看，就相当于材料具有了负等效弹性模量。在设计这种结构时，需要考虑铰链的尺寸、形状以及材料的力学性能等因素，以确保在合适的频率范围内实现稳定的负等效弹性模量。变换光学理论也为单相固体双负弹性波超材料的设计提供了新的思路。变换光学理论的基础是基于麦克斯韦方程组在坐标变换下的不变性，通过对空间坐标进行特定的变换，将所需的弹性波传播路径映射到材料的参数分布上，从而设计出具有特定功能的超材料。在弹性波超材料中，类比于电磁学中的变换光学，通过弹性动力学方程在坐标变换下的特性，实现对弹性波传播的精确控制。具体的设计步骤如下：首先，根据所需的弹性波传播功能，如负折射、聚焦等，确定目标空间中的弹性波传播路径和波前分布。然后，选择合适的坐标变换函数，将目标空间的坐标变换到实际材料空间的坐标。在这个过程中，需要根据弹性动力学方程，推导出变换后的材料参数，包括等效质量密度和等效弹性模量的分布。通过数值模拟和优化算法，对设计的超材料结构进行验证和优化，确保其能够实现预期的弹性波传播特性。例如，在设计基于变换光学理论的负折射超材料时，通过坐标变换，使弹性波在超材料中的传播路径发生弯曲，实现负折射现象，从而验证了变换光学理论在单相固体双负弹性波超材料设计中的有效性。2.4设计案例分析以一种基于局域共振原理设计的单相固体双负弹性波超材料为例，深入分析其设计过程中的关键环节。该超材料旨在实现低频段的双负特性，以满足在低频减振降噪领域的应用需求。在参数选择方面，首先考虑的是微观结构中谐振单元的尺寸参数。谐振单元采用了质量-弹簧结构，其中质量块的质量m和弹簧的刚度k是关键参数。根据局域共振理论，谐振单元的共振频率f_0与质量m和弹簧刚度k的关系为f_0=\frac{1}{2\pi}\sqrt{\frac{k}{m}}。为了使超材料在低频段（如100-500Hz）实现双负特性，通过理论计算和初步的数值模拟，初步确定质量块的质量范围为0.1-0.5kg，弹簧刚度范围为10-50N/m。在这个范围内，不同的质量和刚度组合会导致共振频率的变化，进而影响超材料的等效参数和双负特性的实现频率范围。例如，当质量m=0.2kg，弹簧刚度k=20N/m时，计算得到共振频率f_0\approx1.59Hz，而当质量m=0.3kg，弹簧刚度k=30N/m时，共振频率f_0\approx1.6Hz，通过对比不同组合下的共振频率，结合目标低频段，最终确定合适的质量和刚度值。材料属性的选择也至关重要。质量块选用密度较大、力学性能稳定的金属材料，如铜，其密度为8960kg/m³，具有较高的质量密度，能够在较小的体积下提供较大的质量，有利于实现低频共振；弹簧选用弹性模量适中、疲劳性能好的弹簧钢，其弹性模量为206GPa，能够在保证一定弹性回复力的同时，满足谐振单元的振动要求。在结构优化环节，首先对谐振单元的排列方式进行优化。最初设计采用简单的正方形晶格排列，通过数值模拟发现，在某些频率下，弹性波在这种排列方式的超材料中传播时，会出现能量集中和散射不均匀的问题。为了解决这个问题，尝试采用六边形晶格排列方式。六边形晶格排列具有更好的对称性和空间填充效率，能够使弹性波在超材料中更均匀地传播，减少能量集中现象。通过对比两种排列方式下超材料的能带结构和弹性波传输谱，发现六边形晶格排列的超材料在目标低频段内，能带结构更加平坦，弹性波的传输损耗更小，双负特性更加稳定。对谐振单元的形状进行优化。将质量块的形状从简单的立方体优化为带有圆角的长方体，圆角的设计可以减少应力集中，提高结构的力学性能和稳定性。同时，对弹簧的形状进行改进，采用变截面弹簧，即在弹簧的不同部位采用不同的截面尺寸，这种设计可以进一步调整弹簧的刚度分布，使谐振单元的振动更加灵活，增强超材料的双负特性。通过有限元分析软件对优化前后的结构进行力学性能分析，结果表明，优化后的结构在相同载荷下的应力分布更加均匀，最大应力值降低了约20%，有效提高了结构的可靠性和稳定性。通过上述参数选择和结构优化过程，最终设计出的单相固体双负弹性波超材料在100-500Hz的低频段内成功实现了稳定的双负特性，为低频减振降噪应用提供了有效的材料解决方案。这一设计案例充分展示了在单相固体双负弹性波超材料设计中，合理的参数选择和精细的结构优化对于实现预期性能的重要性。三、基于COMSOL有限元软件的弹性波超材料计算方法3.1弹性波和声波基本方程弹性波作为一种机械波，其在介质中的传播行为遵循特定的物理规律，这些规律可以通过基本方程来描述。从牛顿第二定律出发，考虑介质的应力-应变关系以及几何关系，能够推导得到弹性波在各向同性均匀介质中的运动微分方程。设介质中的位移矢量为\vec{u}=(u_x,u_y,u_z)，应力张量为\sigma_{ij}（i,j=1,2,3，分别对应x,y,z方向），体力密度为\vec{f}=(f_x,f_y,f_z)，根据牛顿第二定律F=ma，在微小体积元内，介质所受的合力等于质量与加速度的乘积。对于弹性介质，应力与应变之间存在线性关系，即广义胡克定律\sigma_{ij}=C_{ijkl}\epsilon_{kl}，其中C_{ijkl}为弹性常数张量，\epsilon_{kl}为应变张量，且应变与位移的关系为\epsilon_{kl}=\frac{1}{2}(\frac{\partialu_k}{\partialx_l}+\frac{\partialu_l}{\partialx_k})。通过对这些关系进行推导和整理，最终得到弹性波在各向同性均匀介质中的运动微分方程为：\mu\nabla^2\vec{u}+(\lambda+\mu)\nabla(\nabla\cdot\vec{u})+\vec{f}=\rho\frac{\partial^2\vec{u}}{\partialt^2}其中，\lambda和\mu为拉梅常数，与材料的弹性模量E和泊松比\nu相关，具体关系为\lambda=\frac{E\nu}{(1+\nu)(1-2\nu)}，\mu=\frac{E}{2(1+\nu)}；\rho为介质的密度；\nabla^2为拉普拉斯算子，\nabla^2=\frac{\partial^2}{\partialx^2}+\frac{\partial^2}{\partialy^2}+\frac{\partial^2}{\partialz^2}；\nabla为梯度算子，\nabla\cdot\vec{u}表示位移的散度。这个方程描述了弹性波在介质中传播时，位移随空间和时间的变化规律，是研究弹性波传播特性的基础。声波是弹性波的一种特殊形式，主要在流体介质中传播。对于理想流体介质，忽略黏性和热传导等因素，声波的传播可以通过质量守恒方程、运动方程和物态方程来描述。质量守恒方程表示在单位时间内，流体微元内质量的变化等于流入和流出该微元的质量之差，其数学表达式为\frac{\partial\rho}{\partialt}+\nabla\cdot(\rho\vec{v})=0，其中\vec{v}为流体的速度矢量；运动方程基于牛顿第二定律，描述了流体微元所受的合力与加速度的关系，即\rho\frac{\partial\vec{v}}{\partialt}+\rho(\vec{v}\cdot\nabla)\vec{v}=-\nablap，其中p为流体的压强；物态方程描述了流体的压强、密度和温度等状态参数之间的关系，对于理想气体，在等温过程中，物态方程为p=p_0(\frac{\rho}{\rho_0})^{\gamma}，其中p_0和\rho_0为初始压强和密度，\gamma为比热容比。通过对这些方程进行适当的简化和推导，在小振幅声波的假设下（即声波引起的密度和压强变化远小于介质的初始值），忽略非线性项，得到线性化的声波波动方程为：\frac{1}{c^2}\frac{\partial^2p}{\partialt^2}-\nabla^2p=0其中，c=\sqrt{\frac{\gammap_0}{\rho_0}}为声波在介质中的传播速度。这个方程表明，在理想流体介质中，小振幅声波的传播满足波动方程，其传播速度取决于介质的初始压强、密度和比热容比等参数。弹性波和声波基本方程中的各个参数具有明确的物理意义。拉梅常数\lambda和\mu反映了材料的弹性性质，它们决定了材料在受力时的变形方式和弹性回复能力。例如，\mu越大，材料抵抗剪切变形的能力越强；\lambda则在描述材料的体积变形和纵向变形时起到重要作用。介质的密度\rho体现了单位体积内介质的质量，它对弹性波和声波的传播速度有显著影响，一般来说，密度越大，波的传播速度越慢。声波传播速度c是一个关键参数，它决定了声波在介质中传播的快慢，同时也与介质的弹性和惯性特性密切相关。这些参数的相互作用，共同决定了弹性波和声波在介质中的传播特性，如传播速度、波长、频率等，对理解和研究弹性波超材料的性能具有重要意义。3.2弹性波超材料数值计算方法在研究弹性波超材料的性能时，能带结构的计算是深入理解弹性波传播特性的关键。能带结构描述了弹性波在超材料中的频率与波矢之间的关系，通过计算能带结构，可以确定超材料中弹性波的传播禁带和通带，这对于分析超材料对弹性波的调控能力至关重要。平面波展开法是计算弹性波超材料能带结构的常用方法之一。其基本原理基于傅里叶变换，将超材料的周期性结构和弹性波的位移场展开为平面波的叠加。对于具有周期性晶格结构的超材料，设晶格矢量为\vec{a}_1和\vec{a}_2（二维情况）或\vec{a}_1、\vec{a}_2和\vec{a}_3（三维情况），则超材料的周期性函数f(\vec{r})（如质量密度\rho(\vec{r})和弹性常数C_{ijkl}(\vec{r})）可以展开为：f(\vec{r})=\sum_{\vec{G}}f_{\vec{G}}e^{i\vec{G}\cdot\vec{r}}其中，\vec{G}为倒格矢，满足\vec{G}\cdot\vec{a}_i=2\pin_i（n_i为整数）；f_{\vec{G}}为傅里叶系数。将弹性波的位移场\vec{u}(\vec{r},t)=\vec{u}_{\vec{q}}e^{i(\vec{q}\cdot\vec{r}-\omegat)}（\vec{q}为波矢，\omega为角频率）代入弹性波运动微分方程，并利用超材料的周期性，经过一系列数学推导和变换，可得到一个关于\vec{u}_{\vec{q}}的线性方程组。通过求解该线性方程组的本征值问题，得到不同波矢\vec{q}对应的本征频率\omega，从而绘制出能带结构。在COMSOLMultiphysics软件中实现平面波展开法计算能带结构时，首先需要在软件中建立超材料的周期性几何模型，定义材料的属性，包括质量密度和弹性常数等。然后，利用软件中的“波动光学模块”或“声学模块”等相关物理场接口，设置平面波展开的参数，如倒格矢的截断数等。在计算过程中，软件会自动进行傅里叶变换和本征值求解，最终得到能带结构数据，并可以通过软件的后处理功能绘制出能带结构图。等频率曲线反映了弹性波在超材料中传播时，相同频率下波矢的分布情况，它对于研究弹性波的传播方向和相速度等特性具有重要意义。通过等频率曲线，可以直观地观察到弹性波在超材料中的传播各向异性以及负折射等特殊现象。有限元方法是计算等频率曲线的一种有效手段。在COMSOLMultiphysics中，利用有限元方法计算等频率曲线的步骤如下：首先，建立超材料的二维或三维模型，并对模型进行网格划分，确保网格的质量满足计算精度要求。选择合适的物理场接口，如“固体力学”接口，定义材料的弹性参数和边界条件。然后，在特定频率下，对模型进行稳态分析，求解弹性波的位移场分布。通过后处理功能，提取模型中不同位置处的波矢信息，根据波矢与频率的关系，绘制出等频率曲线。在绘制等频率曲线时，需要注意选择合适的频率范围和波矢采样点。频率范围应覆盖超材料的主要特性频率区域，波矢采样点的分布要足够密集，以保证等频率曲线的准确性和光滑度。对于一些具有复杂结构的超材料，可能需要对模型进行适当的简化或近似处理，以提高计算效率和收敛性。传输谱描述了弹性波在超材料中传播时，透射系数或反射系数随频率的变化情况，它是评估超材料对弹性波传输特性的重要参数。通过分析传输谱，可以确定超材料的带隙位置和宽度，以及在通带内弹性波的传输效率等信息。在COMSOLMultiphysics中计算传输谱，通常采用频域分析方法。建立包含超材料结构和入射波源的模型，设置入射波的类型（如平面波、点源等）和频率范围。定义材料的属性和边界条件，其中边界条件的设置对于准确模拟弹性波的传输至关重要。在模型的边界上，需要设置合适的吸收边界条件或周期性边界条件，以避免弹性波的反射和散射对计算结果的影响。然后，在不同频率下对模型进行求解，得到弹性波在超材料中的传播场分布。通过后处理功能，提取模型输出边界处的弹性波幅值信息，根据透射系数和反射系数的定义，计算出不同频率下的传输谱。等效参数是描述超材料宏观特性的重要物理量，包括等效质量密度和等效弹性模量等。通过计算等效参数，可以将超材料看作是一种具有等效性质的均匀介质，从而便于对其进行理论分析和工程应用。均匀化方法是计算超材料等效参数的常用方法之一。其基本思想是基于体积平均的概念，通过对超材料微观结构中的物理量进行体积平均，得到宏观等效的物理参数。以等效质量密度的计算为例，假设超材料的微观结构由基体和包含物组成，在一个代表性体积单元（RVE）内，质量密度分布为\rho(\vec{r})，则等效质量密度\rho_{eq}可以通过下式计算：\rho_{eq}=\frac{1}{V}\int_{V}\rho(\vec{r})dV其中，V为代表性体积单元的体积。对于等效弹性模量的计算，需要考虑超材料在不同加载方向上的应力-应变关系。通过对RVE施加不同方向的均匀应变，求解弹性波运动微分方程，得到相应的应力分布，再根据弹性模量的定义，通过体积平均计算出等效弹性模量。在COMSOLMultiphysics中，可以利用“材料库”中的“均匀化材料”功能来实现等效参数的计算。首先，建立超材料的微观结构模型，并对模型进行网格划分。然后，定义材料的属性和边界条件，在边界条件中设置周期性边界条件，以模拟超材料的无限周期性。在“均匀化材料”设置中，选择需要计算的等效参数类型，如等效弹性模量、等效质量密度等，并设置相关的计算参数。软件会自动根据均匀化方法，对微观结构模型进行计算，得到超材料的等效参数。3.3计算流程与参数设置利用COMSOL软件进行弹性波超材料计算时，需遵循严谨的流程，以确保计算结果的准确性和可靠性。在模型建立环节，根据设计的单相固体双负弹性波超材料的微观结构，在COMSOL软件的几何建模模块中进行精确构建。对于基于局域共振原理设计的超材料，若其基本单元为质量-弹簧结构，首先在二维或三维建模空间中绘制出质量块和弹簧的几何形状。利用软件提供的绘图工具，如矩形、圆形等基本图形绘制功能，绘制出质量块的形状（如矩形），并通过拉伸、旋转等操作构建出三维模型。对于弹簧，可利用样条曲线工具绘制其螺旋形状，再通过拉伸和旋转操作生成三维螺旋弹簧模型。定义质量块和弹簧的尺寸参数，如质量块的长度、宽度、高度，弹簧的直径、螺距、簧丝直径等，并设置这些参数为可变量，以便后续进行参数化研究。利用布尔运算工具，将质量块和弹簧组合成一个完整的超材料基本单元，并通过阵列复制功能，按照特定的晶格排列方式（如正方形晶格或六边形晶格），生成周期性的超材料模型。完成模型建立后，需对模型进行网格划分，这是影响计算精度和计算效率的关键步骤。在COMSOL软件中，选择合适的网格划分方法，对于几何形状较为规则的超材料模型，可采用结构化网格划分方法，这种方法生成的网格具有规则的形状和良好的正交性，能够提高计算精度。对于复杂的超材料结构，可采用非结构化网格划分方法，如三角形或四面体网格，以更好地适应模型的几何形状。在网格参数设置中，根据模型的尺寸和计算精度要求，调整网格的大小和密度。对于超材料中的关键部位，如谐振单元的连接处、应力集中区域等，适当减小网格尺寸，增加网格密度，以提高这些部位的计算精度；对于模型中对计算结果影响较小的区域，可适当增大网格尺寸，降低网格密度，以减少计算量和计算时间。还可以利用软件的自适应网格划分功能，根据计算过程中物理量的变化情况，自动调整网格的疏密程度，进一步提高计算精度和效率。边界条件的设置对于准确模拟弹性波在超材料中的传播至关重要。在COMSOL软件中，根据实际的物理问题和研究需求，设置合适的边界条件。对于周期性结构的超材料模型，在模型的边界上设置周期性边界条件，以模拟超材料在无限空间中的周期性特性。在设置周期性边界条件时，需确保边界上的位移、应力等物理量满足周期性条件，即边界两侧的物理量在相应的位置上相等。对于弹性波的入射边界，设置入射波边界条件，定义入射波的类型（如平面波、点源波等）、频率、振幅和相位等参数。若研究平面波垂直入射到超材料模型的情况，在入射边界上选择平面波入射边界条件，并设置平面波的频率范围（如100-1000Hz）、振幅为1Pa、相位为0。为了避免弹性波在模型边界上的反射对计算结果的影响，在模型的输出边界上设置完美匹配层（PML）边界条件，PML边界条件能够有效地吸收传播到边界的弹性波，模拟无限空间的传播环境。在设置PML边界条件时，需根据模型的尺寸和计算频率范围，合理调整PML的厚度和吸收系数，以确保其能够有效地吸收弹性波。在计算过程中，还需设置求解器参数。选择合适的求解器类型，对于频域分析问题，可选择频域求解器；对于时域分析问题，可选择时域求解器。设置求解器的收敛精度，一般将收敛精度设置为1e-6-1e-8，以确保计算结果的准确性。根据模型的规模和计算机的性能，设置求解器的最大迭代次数和内存分配等参数，以避免计算过程中出现不收敛或内存不足的问题。3.4计算结果分析与验证通过COMSOL软件的计算，得到了所设计单相固体双负弹性波超材料的能带结构、等频率曲线、传输谱以及等效参数等关键结果。对这些结果进行深入分析，并与理论结果或已有研究进行对比，以验证计算方法的准确性。从能带结构计算结果来看，在特定的频率范围内，出现了明显的禁带。以基于局域共振原理设计的超材料为例，在共振频率附近，能带结构呈现出特殊的变化，这与理论预期相符。理论上，基于局域共振的超材料在共振频率处，由于谐振单元的强烈共振，会对弹性波产生强烈的散射和吸收，从而形成禁带。将计算得到的能带结构与[具体文献11]中基于相同设计原理的超材料能带结构进行对比，发现两者在禁带位置和宽度上具有较高的一致性。在[具体文献11]中，通过理论计算和实验测量，确定了超材料在150-250Hz频率范围内存在禁带，而本次COMSOL计算结果显示，该超材料在145-255Hz频率范围内出现禁带，两者误差在可接受范围内，这表明COMSOL计算得到的能带结构准确地反映了超材料的弹性波传播特性，验证了计算方法在能带结构计算方面的准确性。等频率曲线的计算结果也为超材料的弹性波传播特性提供了重要信息。通过对不同频率下的等频率曲线进行分析，发现超材料在某些频率范围内表现出明显的各向异性特性。在低频段，等频率曲线呈现出近似圆形的形状，表明弹性波在超材料中的传播具有各向同性的特点；而在高频段，等频率曲线发生了明显的变形，呈现出椭圆形或不规则形状，说明弹性波的传播速度在不同方向上存在差异，具有各向异性。与[具体文献12]中的研究结果相比，该文献通过理论分析和实验测量，得出了类似的等频率曲线变化规律，进一步验证了COMSOL计算结果的正确性。在[具体文献12]中，对等频率曲线的各向异性变化进行了详细的理论分析，指出超材料的微观结构和晶格排列方式是导致等频率曲线各向异性的主要原因，本次计算结果与该理论分析结果一致，证明了计算方法在等频率曲线计算方面的可靠性。传输谱的计算结果直观地展示了弹性波在超材料中的传输特性。在带隙频率范围内，弹性波的透射系数极低，几乎为零，表明弹性波在该频率范围内被有效地抑制，无法在超材料中传播；而在通带频率范围内，弹性波的透射系数较高，能够顺利通过超材料。将计算得到的传输谱与[具体文献13]中的实验测量传输谱进行对比，发现两者在带隙和通带的位置、宽度以及透射系数的变化趋势上基本一致。在[具体文献13]中，通过实验测量得到超材料在200-300Hz频率范围内为带隙，透射系数小于0.1，在其他频率范围内为通带，透射系数大于0.5，本次COMSOL计算结果在相应频率范围内的带隙和通带特性与实验测量结果相符，验证了计算方法在传输谱计算方面的准确性。等效参数的计算结果对于将超材料等效为均匀介质进行分析具有重要意义。通过均匀化方法计算得到的超材料等效质量密度和等效弹性模量，在特定频率范围内呈现出负值，这与单相固体双负弹性波超材料的设计目标一致。将计算得到的等效参数与[具体文献14]中基于理论模型计算得到的等效参数进行对比，发现两者在数值大小和变化趋势上具有较好的一致性。在[具体文献14]中，通过理论模型推导出了超材料等效参数与结构参数之间的关系，并计算得到了在特定结构参数下的等效参数值，本次COMSOL计算结果与该理论计算值的相对误差在10%以内，证明了计算方法在等效参数计算方面的可靠性。通过对能带结构、等频率曲线、传输谱以及等效参数等计算结果的分析，并与理论结果和已有研究进行对比，充分验证了基于COMSOL有限元软件的弹性波超材料计算方法的准确性和可靠性，为进一步研究超材料的性能和应用提供了有力的工具和依据。四、双负弹性波超材料的优化设计及成像特性研究4.1几何参数对材料性能的影响双负弹性波超材料的性能与几何参数密切相关，深入研究这些参数对材料性能的影响，对于优化超材料设计、提升其性能具有重要意义。几何参数主要包括结构形状、尺寸等，它们的变化会显著改变超材料的能带结构、等频率曲线及等效参数。以基于局域共振原理设计的双负弹性波超材料为例，其基本结构单元通常包含质量块与弹簧。质量块的形状对超材料性能有显著影响。当质量块为正方体时，在特定频率下，其与弹性波的相互作用使得超材料的等效质量密度在共振频率附近出现负值。若将质量块形状改为球体，由于球体的对称性和质量分布均匀性，在相同频率下，等效质量密度的负值范围和大小发生改变。通过数值模拟计算，在共振频率为200Hz时，正方体质量块的超材料等效质量密度在180-220Hz频率范围内为负值，最小值达到-0.5kg/m³；而球体质量块的超材料等效质量密度在190-210Hz频率范围内为负值，最小值为-0.3kg/m³。这表明质量块形状的改变会影响超材料等效质量密度的负特性频率范围和幅值，进而影响其对弹性波的调控能力。弹簧的形状也会对超材料性能产生重要影响。传统的圆柱形弹簧在弹性波作用下，其振动模式较为单一。若采用变截面弹簧，在弹性波传播过程中，由于不同截面处的刚度差异，弹簧的振动模式变得更加复杂。这种复杂的振动模式会增强弹簧与质量块之间的相互作用，从而改变超材料的等效弹性模量。在COMSOL软件中模拟发现，当采用等截面弹簧时，超材料在300Hz左右出现负等效弹性模量，负模量范围为-10--5GPa；而采用变截面弹簧后，负等效弹性模量出现在280-320Hz频率范围内，负模量范围扩大到-15--8GPa，且在该频率范围内，能带结构中的禁带宽度也有所增加，从原来的30Hz增加到40Hz，这表明变截面弹簧能够在更宽的频率范围内实现双负特性，并且增强了超材料对弹性波的抑制能力。尺寸参数同样对超材料性能有着关键影响。质量块的尺寸变化会直接影响其惯性，进而影响超材料的等效质量密度。当质量块尺寸增大时，其惯性增大，共振频率降低。在一个质量-弹簧结构的超材料单元中，质量块边长从10mm增大到20mm时，共振频率从300Hz降低到200Hz，等效质量密度在更低的频率范围内出现负值。这说明通过调整质量块尺寸，可以灵活地调控超材料双负特性的频率范围，以满足不同应用场景对频率的需求。弹簧的尺寸，如弹簧丝直径和弹簧圈数等，也会对超材料性能产生影响。弹簧丝直径增大，弹簧的刚度增大，共振频率升高。当弹簧丝直径从1mm增大到2mm时，共振频率从250Hz升高到350Hz，等效弹性模量的负值范围和频率位置也发生相应变化。弹簧圈数的增加会使弹簧的柔性增强，在相同外力作用下，弹簧的变形量增大，从而改变超材料的等效弹性模量。当弹簧圈数从5圈增加到8圈时，超材料的等效弹性模量在低频段的负值更加明显，这表明通过调整弹簧的尺寸参数，可以优化超材料的双负性能，使其在不同频率范围内更好地发挥对弹性波的调控作用。结构形状和尺寸等几何参数对双负弹性波超材料的能带结构、等频率曲线及等效参数有着显著的影响。通过合理设计和调整这些几何参数，可以精确调控超材料的性能，为其在声学和力学领域的应用提供更优化的材料选择。4.2模型的本征模态分析本征模态分析是深入理解双负弹性波超材料内部振动特性的重要手段，通过该分析可揭示材料内部振动模式与双负特性之间的内在联系。以基于局域共振原理设计的超材料模型为例，在COMSOL软件中进行本征模态分析。在COMSOL软件中，选择“结构力学”模块下的“模态分析”功能。对已建立的超材料模型进行分析设置，定义材料的密度、弹性模量等参数，确保参数设置与实际设计一致。在分析选项中，设置求解的模态阶数，一般选择前几阶模态进行重点分析，因为前几阶模态通常对材料的宏观性能影响较大。设置好参数后，软件开始计算超材料模型的本征频率和振型。计算结果显示，在不同的本征模态下，超材料内部呈现出不同的振动模式。在某一阶模态下，超材料中的质量块和弹簧结构表现出协同振动的模式。质量块在弹簧的作用下，以特定的频率和振幅进行振动，弹簧则在质量块的带动下发生拉伸和压缩变形。这种振动模式与双负特性密切相关，当弹性波的频率接近该本征模态的频率时，谐振单元的共振效应增强，导致超材料的等效质量密度和等效弹性模量发生变化，进而呈现出双负特性。通过对不同频率下的本征模态进行分析，发现随着频率的变化，超材料内部的振动模式也会发生改变。在低频段，质量块的振动幅度相对较小，弹簧的变形也较为均匀，此时超材料的等效参数变化相对较小；而在高频段，质量块的振动幅度增大，弹簧的变形更加复杂，超材料的等效参数在某些频率范围内出现明显的负值，表现出双负特性。进一步分析本征模态与双负特性的关系，绘制等效质量密度和等效弹性模量随本征频率的变化曲线。从曲线中可以看出，在某些本征频率附近，等效质量密度和等效弹性模量同时出现负值，这些频率区域正是超材料呈现双负特性的频率范围。这表明超材料的本征模态与双负特性之间存在着紧密的联系，本征模态的振动模式决定了超材料在不同频率下的等效参数变化，从而影响其双负特性的表现。本征模态分析揭示了双负弹性波超材料内部振动模式与双负特性的关系，为深入理解超材料的物理机制提供了重要依据。通过对本征模态的研究，可以进一步优化超材料的设计，提高其双负性能和应用效果。4.3弹性超材料成像及其分辨率分析弹性超材料成像基于其独特的双负特性，为实现高分辨率成像提供了新的途径。当弹性波从点源发出，在双负弹性超材料中传播时，由于其负等效质量密度和负等效弹性模量，弹性波的传播特性发生改变，能够实现对弹性波的聚焦和亚波长成像。在传统成像中，由于衍射极限的存在，成像分辨率受到限制，而双负弹性超材料可以突破这一限制，实现更高分辨率的成像。为了深入研究超材料成像特性，建立超材料成像模型，在模型中设置点源，模拟弹性波的发射。通过改变超材料的宽度，分析其对像点分辨率的影响。当超材料宽度较小时，弹性波在超材料中传播时，能量损失较大，像点的分辨率较低，像点呈现出较大的弥散斑，细节信息模糊；随着超材料宽度的增加，弹性波在超材料中的传播路径更加稳定，能量损失减小，像点的分辨率逐渐提高，像点的弥散斑变小，能够更清晰地分辨出点源的位置和形状。当超材料宽度从10mm增加到30mm时，像点的半高宽从0.8mm减小到0.3mm，分辨率提高了约62.5%，这表明增加超材料宽度可以有效提高成像分辨率。点源位置的变化也会对像点分辨率产生显著影响。当点源位于超材料的中心位置时，弹性波在超材料中对称传播，像点的分辨率较高；而当点源偏离中心位置时，弹性波在超材料中的传播路径发生不对称变化，像点的分辨率下降，像点出现偏移和变形。当点源偏离中心位置5mm时，像点的偏移量达到0.2mm，分辨率降低了约30%，这说明点源位置的精确控制对于提高成像分辨率至关重要。点源频率同样是影响像点分辨率的重要因素。在不同频率下，超材料的双负特性表现不同，从而影响成像分辨率。在低频段，超材料的双负特性较弱，弹性波的传播损耗较大，像点分辨率较低；随着频率的增加，超材料的双负特性增强，弹性波的传播损耗减小，像点分辨率提高。当频率从100Hz增加到300Hz时，像点的分辨率提高了约40%，但当频率继续增加到一定程度后，由于超材料内部结构的共振等因素，像点分辨率又会出现下降趋势。这表明存在一个最佳频率范围，在该范围内超材料能够实现最佳的成像分辨率。五、双负弹性波超材料负折射的数值模拟及实验研究5.1二维体波的负折射模拟为深入探究双负弹性波超材料中负折射现象的规律和特性，利用COMSOLMultiphysics软件对二维体波在超材料中的传播进行数值模拟，重点研究不同位置尺寸和频率下纵波的负折射情况。建立基于局域共振原理设计的二维双负弹性波超材料数值模型。模型采用正方形晶格排列，每个晶格单元包含一个质量块和连接质量块与基体的弹簧。在模拟不同位置尺寸下纵波的负折射时，固定超材料的其他参数，如材料属性、晶格常数等，仅改变质量块的位置和尺寸。当质量块位于晶格单元的中心位置时，设置质量块的边长分别为0.05m、0.1m和0.15m，分析纵波在超材料中的传播情况。模拟结果显示，在相同频率下，随着质量块尺寸的增大，纵波的负折射角度发生变化。当质量块边长为0.05m时，在频率为300Hz的纵波入射下，负折射角度为-30^{\circ}；当边长增大到0.1m时，负折射角度变为-35^{\circ}；边长为0.15m时，负折射角度达到-40^{\circ}。这表明质量块尺寸的增大，会增强超材料对纵波的负折射效果，使得负折射角度增大。将质量块的位置从晶格单元中心偏移至靠近晶格边界的位置，偏移距离分别为晶格边长的0.1倍、0.2倍和0.3倍，观察纵波负折射特性的变化。结果发现，随着质量块向边界偏移，纵波的负折射现象变得不稳定，在某些频率下，负折射角度出现波动，且负折射效率有所降低。这说明质量块在晶格单元中的位置对纵波的负折射有显著影响，偏离中心位置会破坏超材料微观结构的对称性，进而影响负折射性能。在模拟不同频率下纵波的负折射时，固定超材料的结构参数，将入射纵波的频率从100Hz逐渐增加到500Hz，步长为50Hz。模拟结果表明，在双负特性频率范围内，随着频率的增加，纵波的负折射角度先增大后减小。在频率为200Hz时，负折射角度达到最大值-45^{\circ}；当频率超过350Hz后，负折射角度逐渐减小，在500Hz时，负折射角度减小至-25^{\circ}。这是因为超材料的双负特性在不同频率下表现不同，随着频率的变化，超材料内部谐振单元的共振状态发生改变，导致等效参数变化，从而影响纵波的负折射角度。在高频段，由于超材料内部结构的阻尼和散射等因素的影响，纵波的传播损耗增大，负折射效果减弱。通过对不同位置尺寸和频率下纵波在双负弹性波超材料中负折射现象的模拟，揭示了超材料结构参数和频率对负折射特性的影响规律，为进一步优化超材料设计、提高负折射性能提供了重要的理论依据。5.2板中的弹性波负折射模拟在研究双负弹性波超材料的负折射特性时，板中的弹性波负折射模拟为深入理解弹性波在不同维度和结构中的传播行为提供了独特视角。通过与二维体波负折射进行对比分析，能够更全面地揭示弹性波超材料的负折射规律和特性差异。在COMSOLMultiphysics软件中建立双负弹性波超材料板的数值模型。模型采用二维平面应变假设，模拟弹性波在薄板中的传播情况。板的材料属性设置为具有双负特性的参数，晶格结构同样采用基于局域共振原理设计的单元结构，如质量-弹簧结构。在模拟过程中，设置平面波垂直入射到超材料板上，分析弹性波在板中的传播路径和折射情况。模拟结果显示，板中的弹性波负折射现象与二维体波负折射存在一定的相似性和差异。在相似性方面，两者在双负特性频率范围内都能够实现负折射，即折射波的传播方向与传统正折射率材料中的折射方向相反。在频率为250Hz时，无论是二维体波还是板中的弹性波，在双负弹性波超材料中都表现出明显的负折射，折射角均为负值。这表明双负弹性波超材料的负折射特性在不同维度和结构中具有一定的普适性，主要源于其负等效质量密度和负等效弹性模量对弹性波传播的特殊影响。两者也存在显著差异。在二维体波中，弹性波的传播空间更为自由，其负折射特性主要受超材料的微观结构和频率影响；而在板中，由于板的几何限制，弹性波的传播受到边界条件的约束，其负折射特性还与板的厚度、边界条件等因素密切相关。当板的厚度较小时，弹性波在传播过程中与板的上下表面相互作用增强，导致负折射角度发生变化。在板厚为5mm时，频率为300Hz的弹性波负折射角度为-32°；当板厚增加到10mm时，负折射角度变为-36°，这说明板的厚度变化会对弹性波在板中的负折射特性产生显著影响。板的边界条件对弹性波负折射也有重要作用。在固定边界条件下，弹性波在板中的传播受到更大的限制，负折射效率相对较低；而在自由边界条件下，弹性波的传播相对更自由，负折射效率有所提高。在固定边界条件下，频率为350Hz的弹性波透射系数为0.3，而在自由边界条件下，透射系数提高到0.45，这表明边界条件的改变会影响弹性波在板中的能量传输和负折射效果。通过板中的弹性波负折射模拟以及与二维体波负折射的对比，深入揭示了弹性波在不同结构中的负折射特性，为进一步优化超材料板的设计、拓展其在声学和力学领域的应用提供了重要的理论依据。5.3实验研究方案设计实验试件制备是验证双负弹性波超材料性能的关键步骤。根据数值模拟优化后的结构参数，选择合适的材料和制备工艺。材料选择方面，考虑到实验对材料力学性能和加工性能的要求，选用环氧树脂作为基体材料，其具有良好的成型性和力学性能，密度约为1180kg/m³，弹性模量约为3GPa。对于谐振单元中的质量块，选用密度较大的铅材料，密度为11340kg/m³，以增强谐振效果；弹簧选用不锈钢丝制作，其弹性模量约为193GPa，保证在弹性波作用下能够产生稳定的谐振。采用3D打印技术制备超材料试件。3D打印技术具有能够实现复杂结构制造的优势，能够精确控制超材料微观结构的尺寸和形状。利用专业的3D建模软件，如SolidWorks，根据设计的超材料微观结构模型，进行三维模型构建。在建模过程中，精确设置各个结构参数，如质量块的边长为15mm，弹簧的直径为2mm，弹簧圈数为6圈等。将构建好的三维模型导入3D打印机中，选择合适的打印参数。打印过程中，控制打印温度在150-180℃之间，以确保环氧树脂材料能够均匀熔化和成型；打印速度设置为30-50mm/s，保证打印精度和效率。打印完成后，对试件进行后处理，去除试件表面的支撑结构和杂质，对试件进行打磨和抛光处理，使其表面光滑，尺寸精度达到±0.1mm，以满足实验测试的要求。搭建高精度的实验测试平台是准确测量超材料弹性波传播特性的基础。采用激光超声技术作为主要的测试手段，该技术具有非接触、高分辨率和宽频带等优点，能够精确测量弹性波在超材料中的传播速度和衰减特性。实验平台主要由脉冲激光器、激光干涉仪、数据采集系统和计算机等组成。脉冲激光器用于产生短脉冲激光，照射到超材料试件表面，激发弹性波。激光干涉仪用于测量弹性波引起的试件表面位移变化，将位移信号转换为电信号。数据采集系统采集激光干涉仪输出的电信号，并将其传输到计算机中进行处理和分析。在实验步骤方面，首先将制备好的超材料试件放置在实验平台的测试区域内，调整试件的位置和角度，使其表面与激光束垂直。开启脉冲激光器，设置激光脉冲的能量为5-10mJ，脉冲宽度为10-20ns，激发弹性波。激光干涉仪实时测量弹性波在试件中传播时引起的表面位移变化，数据采集系统以100MHz的采样频率采集位移信号，并将其传输到计算机中。利用计算机中的数据分析软件，对采集到的位移信号进行处理，计算弹性波在超材料中的传播速度和衰减系数。在不同频率下重复上述实验步骤，测量弹性波在不同频率下的传播特性。实验过程中需注意一些事项。为了减少外界干扰对实验结果的影响，实验平台应放置在具有良好隔振和屏蔽功能的实验室内，实验区域周围应避免放置大型机械设备和电子设备。激光束的对准和聚焦对实验结果的准确性至关重要，在实验前需要仔细调整激光束的位置和角度，确保激光束能够准确地照射到试件表面的预定位置，并实现良好的聚焦。在数据采集和处理过程中，需要对采集到的数据进行滤波和降噪处理，去除噪声和干扰信号，提高数据的质量和可靠性。5.4实验结果分析与讨论将实验测量得到的弹性波在双负弹性波超材料中的传播特性与数值模拟结果进行对比，验证数值模拟的准确性。从传播速度对比来看，在100-300Hz频率范围内，数值模拟计算得到的弹性波在超材料中的传播速度与实验测量值较为接近。在频率为200Hz时，数值模拟得到的纵波传播速度为1500m/s，实验测量值为1480m/s，相对误差约为1.3%，处于可接受的误差范围内，这表明数值模拟能够较为准确地预测弹性波在超材料中的传播速度。在带隙特性方面，数值模拟预测超材料在250-350Hz频率范围内存在带隙，实验测量结果显示，在该频率范围内，弹性波的透射系数极低，几乎为零，与数值模拟结果一致，验证了数值模拟对超材料带隙特性预测的准确性。在实验过程中，也出现了一些问题。试件制备过程中，由于3D打印工艺的精度限制，超材料微观结构的实际尺寸与设计尺寸存在一定偏差。质量块的实际边长比设计尺寸偏差了±0.05mm，弹簧的直径偏差了±0.03mm。这种尺寸偏差会影响超材料的等效参数和弹性波传播特性，导致实验结果与数值模拟结果存在一定差异。为解决这一问题，在后续的试件制备中，对3D打印设备进行了更精确的校准，优化打印参数，如提高打印温度的稳定性，将温度波动控制在±2℃以内，减小打印速度的波动，将速度波动控制在±5mm/s以内，以提高打印精度；对打印后的试件进行更严格的尺寸检测，采用高精度的三坐标测量仪对试件的尺寸进行测量，对于尺寸偏差超出允许范围的试件，进行重新打印或修正，确保试件尺寸与设计尺寸的偏差控制在±0.02mm以内，有效提高了实验结果与数值模拟结果的一致性。实验测量过程中，外界环境的干扰也对实验结果产生了一定影响。实验室周围存在的大型机械设备产生的振动和电磁干扰，会影响激光超声测量系统的稳定性，导致测量信号出现噪声和波动。为减少外界干扰，将实验平台放置在具有良好隔振和屏蔽功能的实验室内，在实验平台周围设置了隔振垫和电磁屏蔽罩，有效降低了外界振动和电磁干扰对实验的影响；对激光超声测量系统进行了优化，增加了信号滤波和降噪处理环节，采用数字滤波算法对采集到的信号进行处理，去除噪声和干扰信号，提高了测量信号的质量和可靠性，进一步提高了实验结果的准确性。六、结论与展望6.1研究成果总结本研究围绕单相固体双负弹性波超材料的设计与实验展开，取得了一系列具有重要理论和实践意义的成果。在设计方面，基于局域共振原理，创新性地设计出一种新型的单相固体双负弹性波超材料微观结构。通过深入研究结构参数对等效质量密度和等效弹性模量的影响规律，运用参数化建模方法，成功建立起超材料结构参数与等效参数之间的定量关系模型。在一个基于质量-弹簧结构的超材料设计中，明确了质量块质量、弹簧刚度、单元尺寸、形状及排列方式等参数对等效参数的影响。当质量块质量在0.1-0.5kg范围内变化，弹簧刚度在10-50N/m范围内调整时，超材料的共振频率及双负特性