ACMP 凝聚态物理进展 1687-8124 1687-8108 Hindawi 10.1155 / 2021/3914210 3914210 研究文章 一种新型的多谐振子弹性超材料 诺拉尔丁万德 萨曼艾哈迈迪 1 https://orcid.org/0000-0002-3852-5473 Sedighi 哈米德M。 1 2 Yaghootian 阿明 1 夏 Shengxuan 1 机械工程系 工学院 阿瓦士大学的Shahid Chamran 阿瓦兹 伊朗 scu.ac.ir 2 卓越钻井中心和研究中心 阿瓦士大学的Shahid Chamran 阿瓦兹 伊朗 scu.ac.ir 2021 2 6 2021 2021 6 4 2021 19 5 2021 2 6 2021 2021 版权所有©2021 Saman Ahmadi Nooraldinvand等人。 这是一篇根据知识共享署名许可证发布的开放获取文章，允许在任何媒体中不受限制地使用、分发和复制，前提是原创作品被正确引用。

本文提出了包含一个谐振腔和两个谐振腔的弹性超材料的两种模型，以获得带隙，以提供宽带振动抑制。该模型由多个单元组成，每个单元由一个矩形框架和一个齿轮齿条机构组成，齿轮齿条机构在框架两侧用线性弹簧连接。在二自由度模型中，增加一个小质量，其中心一端附在小齿轮的中心，另一端通过线性弹簧与矩形框架相连。在第一种机构中，小齿轮被视为单个谐振器，而在二自由度模型中，小齿轮和小质量被视为多个谐振器。通过获得各模型中单细胞的运动控制方程，深入研究了两种元结构的动力学行为。因此，将两种模型的运动方程写成矩阵形式，并在此基础上建立色散关系，分析系统参数对带隙起止频率的影响。最后对两个模型进行了比较，并利用MATLAB-SIMULINK和MSC-ADAMS软件进行了数值模拟。借助于开始/结束频率的闭合表达式，可以很容易地识别系统参数与带隙间隔之间的相关性。

 阿瓦士大学的Shahid Chamran SCU.EM99.98 1.介绍

超材料是指综合制成的先进材料类型，包括一般表现为连续材料的小型子结构。声波和弹性波不能传播的频带称为带隙，这是超材料最显著的特征。不同波长的波的传播由超材料产生的低/高频带隙控制。鉴于这种独特的特性，超材料可以用于隔音、波过滤、减少振荡和声波传播领域[ 1- - - - - - 3.］.声学超材料的胞体由一个或多个局部振动部件所在的基本结构组成。基底环境和局部振动子之间的相互作用使有趣的物理现象成为可能，如带宽、负有效密度、负弹性模量、定向滤波和波传导、增加的耗散和振动声衰减[ 4， 5］.基本采用两种局部谐振器:(1)平移谐振器和(2)旋转谐振器，它们分别产生负的有效密度和负的有效弹性模量[ 6］.Huang et al. [ 7证明了使用不同的等效模型来表示由质量中质量单位组成的格系统的结果。然后他们研究了色散波的传播，并与各种等效模型进行了比较。他们发现，如果用经典的弹性连续体来表示原始质量中的质量结构，有效质量密度将与频率有关，并且在内部质量共振频率附近的频率可能变为负的。同时，如果采用多位移微结构连续体模型，可以充分描述波传播的色散行为和带隙结构。

主动周期结构和被动周期结构由相似的子结构或单元细胞组成，这些子结构或单元细胞以相似的方式相互连接。这些结构的周期性表现出独特的动态特性，因为它们可以阻止波的传播，并起到机械过滤器的作用[ 8］.波在无限周期结构中以特定的频率间隔传播是不可能的。光子带隙是声波/弹性波不能传播的频段。与无限周期结构不同的是，如果在带隙内激发有限周期结构，则波的透射或振动很少发生，可以通过检查带隙的频响函数图来获得波或振动的透射[ 9］.传播或通带允许波沿这些结构在一定的频率范围内传播而不会衰减;而阻带、衰减带或带隙是波可以衰减的频带[ 10］.在周期结构中，产生声子带隙的两种常用方法是布拉格散射和局域共振[ 9］.上面提到的带结构和波折减是由于这两种机制明显不同[ 11］.在高频模式下，可以使用布拉格带隙对波进行滤波，但在低频模式下，这是不可能的[ 12].另一方面，在比布拉格散射低得多的频率下，可以通过局部共振产生带隙[ 9］.由于局部共振(LR)机制，谐振型带隙的控制是很容易实现的，不像布拉格型带隙。为了调整带隙特性，可以改变局域谐振结构的特性。陈及谦[ 12]提出了一种质量网格系统，在这一新思想应用于弹性领域之前，它是在声学领域提出的。

Lazarov和Jensen [ 10]以线性和非线性两种模式进行研究，以在非线性振子连接的一维链中传播波。他们发现，在非线性状态下，带隙可以根据局部非线性的幅度和程度发生位移。此外，他们的研究结果表明，在线性模式中，间隙的位置是围绕谐振频率，这使得间隙产生在一个较低的频率域。Casalotti等人[ 13]研究了一个带有非线性质量-弹簧子结构的欧拉-伯努利梁，以研究非线性超材料梁吸收多态振动的能力。作为局部谐振器或吸振器的梁连接子结构产生了带隙。首先，他们识别了单细胞的频率停止带，然后为了研究多频率停止带，提取了非线性超材料光束的频率响应。他们指出，对于超材料光束，适当调整局部吸收器的组成参数可以显著降低与最低三阶振动模态相关的振荡。Zhou et al. [ 14提出了带有柔性带隙和高静刚度的局部谐振子的声超材料束的两种模型。首先考虑了变截面声超材料梁的质量-弹性梁等效模型，并进行了分析。提取了声超材料束的色散关系，研究了不同控制参数对带隙的影响。然后，为了在COMSOL Multiphysics中验证，他们模拟和检验了声学超材料梁的二维有限元模型。结果表明，解析模型与二维模型具有良好的相容性。最后，提出了声学超材料束的三维模型，并讨论了等厚度和不等厚度两种形式测量波的衰减。黄和孙[ 15研究了多腔质量网格系统，得到了它的色散曲线和带隙。结果表明，通过改变内部质量的大小和弹簧常数可以使带隙发生位移。他们提出了一个与原始系统等效的单原子模型，发现带隙频率处的有效质量为负。最后，他们引入了一个微观结构连续模型，可以获得原始系统的色散行为和带隙结构。根据Sun等人的研究[ 16研究表明，基于超材料的吸声器可以被认为是一个均匀的各向同性梁，在不同的点上有许多非常小的质量-弹簧子系统。他们还讨论了如何通过子系统产生带隙、负有效质量和刚度。结果表明，一般的机械吸振器是超材料梁实际工作机理的基础。然后产生剪力和弯矩，使梁保持直，防止波传播。Xiao等人在活性超材料和超器件领域进行了重要和深入的专题综述。[ 17］.他们报告了从微波到可见波长的有源超器件和超材料的进展，包括里程碑和最先进的技术，并最后提出了未来的前景，以及几种新兴的调谐策略和材料。

大多数超材料模型的研究都是基于阻尼器-弹簧-阻尼器组合。然而，本文使用齿轮齿条机构和质量的新组合来实现更宽的带隙，并研究新组合的动态行为。在当前的研究中，提出了两种弹性超材料模型，以获得宽带bandgaps用于振动抑制。模型由连接到集中质量的齿轮齿条机构组合而成。使用MATLAB-SIMULINK求解一个和两个自由度模型的控制方程，并绘制频散曲线和频率响应图。然后，分析控制参数的影响最后，在MSC-ADAMS软件中对这两个模型进行了建模和仿真，并对所需的图形进行了数值提取，仿真结果与数值算法得到的结果吻合良好。

 2.型号:1DOF和2DOF

本节将介绍和讨论由一自由度和二自由度弹性超材料(EMM)组成的两种模型。模型包括包含齿条-齿轮机构和集中质量的矩形单元。齿条的存在防止在小齿轮和电池之间滑动，从而将传递运动转换为纯滚动运动。本研究的目的是分析获得所提出的模型的带隙，以减少或消除机械振动。针对一自由度和二自由度模型，建立了相应的数学模型，并以闭合表达式对带隙起始频率和结束频率进行了数学提取。

2.1.EMM的单自由度模型 2．1．1．一自由度模型:数学公式

建议的弹性超材料单自由度模型如图所示 1. 如图所示，该模型的主体由内部带有齿条和齿轮的矩形框架组成。然后，矩形框架通过线性弹簧连接。另一方面，小齿轮通过两个线性弹簧连接到矩形框架的壁上。

一自由度模型的超材料排列示意图。

为该问题定义的参数如图所示 2，表示单自由度模型的单元格。因此，利用牛顿第二定律，推导出单元格的控制方程: （1） ∑ F x ＝ 米 1 x ¨ n ⇒ − k 1 x n − x n − 1 − k 1 x n − x n + 1 − 米 2 x ¨ n + r θ ¨ ＝ 米 1 x ¨ n ， (2） ∑ 米 O ＝ 我 O θ ¨ ⇒ − 2 k 2 r θ r − 米 2 x ¨ n r ＝ 我 O θ ¨ ， 哪里 米 1 和 米 2 分别为框架质量和小齿轮质量， r 和 θ 分别为小齿轮的半径和转动位移， x n 的平移位移 nth单元，以及 我 O 表示小齿轮相对于该点的质量惯性矩 o．

一自由度模型的单细胞示意图。

现在,方程( 1)和( 2)以矩阵形式表示如下: （3) 米 1 + 米 2 米 2 r 米 2 r 我 O x ¨ n θ ¨ + 2 k 1 0 0 2 k 2 r 2 x n θ + − k 1 x n − 1 + x n + 1 0 ＝ 0 0 ， 模型中使用的无量纲参数定义为 （4) 米 ＝ 米 1 米 2 ， 我 ⌢ ＝ 我 O 米 2 r 2 ， K ⌢ ＝ k 1 k 2 ， ω ⌢ n ＝ k 2 米 2 ， r x ¯ n ＝ x n ， τ ＝ ω t ， ω ˜ ＝ ω ω ⌢ n ， 哪里 ω ⌢ n 为本地固有频率， 米 ， 我 ⌢ ， K ⌢ ， x ¯ n ， τ ， 和 ω ˜ 分别为无量纲质量、转动惯量、刚度、位移、时间、频率。代入式中无量纲参数( 3.)时，新的控制方程形式可表示为: (5） 米 + 1 ω ˜ 2 ω ˜ 2 ω ˜ 2 我 ⌢ ω ˜ 2 d 2 x ¯ n d τ 2 d 2 θ d τ 2 + 2 K ⌢ 0 0 2 x ¯ n θ + − K ⌢ x ¯ n − 1 + x ¯ n + 1 0 ＝ 0 0 ．

2.1.2。一自由度模型:色散特性

本节定义1DOF模型谐波对应解，提取色散关系和带隙频率: （6） x ¯ n θ ＝ x ¯ n ∘ θ ∘ e 我 质量保证 ． n + τ + c ． c ， x ¯ n ± 1 ＝ x ¯ n ± 1 ∘ e 我 质量保证 ． n + τ + c ． c ， 哪里 x ¯ n ∘ 和 θ ∘ 为解的稳态振幅， 质量保证 是相位因数，和 n 表示周期数。此外, 一个 ． n 表示距离 n th 细胞起源(支持)[ 18- - - - - - 22］.用方程( 6) ( 5)，其矩阵形式的控制方程的新形式为: （7） 米 + 1 ω ˜ 2 ω ˜ 2 ω ˜ 2 我 ⌢ ω ˜ 2 x ¯ n ∘ − 1 e 我 质量保证 ． n + τ + c ． c θ ∘ − 1 e 我 质量保证 ． n + τ + c ． c + 2 K ⌢ 0 0 2 x ¯ n ∘ e 我 质量保证 ． n + τ + c ． c θ ∘ e 我 质量保证 ． n + τ + c ． c + − K ⌢ x ¯ n ∘ e 我 质量保证 n − 1 e 我 τ + x ¯ n ∘ e 我 质量保证 n + 1 e 我 τ + c ． c 0 ＝ 0 0 ．

现在，通过分解和整理公式( 7)时，所考虑系统的运动方程的最终形式为 （8） − 米 + 1 ω ˜ 2 + 2 K ⌢ − K ⌢ 2 因为 质量保证 − ω ˜ 2 − ω ˜ 2 − 我 ⌢ ω ˜ 2 + 2 x ¯ n ∘ θ ∘ e 我 质量保证 ． n + τ + c ． c ＝ 0 0 ．

单自由度模型的色散关系可由式( 8）.对于一个非平凡解，系数矩阵的行列式必须等于零。因此，将色散关系提取为 (9） 我 ⌢ 1 + 米 − 1 ω ˜ 4 + 2 − 米 − 1 − 我 ⌢ K ⌢ 1 + 因为 质量保证 ω ˜ 2 + 4 K ⌢ 1 − 因为 质量保证 ＝ 0

通过解方程( 9)为了 ω ˜ 时，EMM的声频率与光频率的关系可由式( 10)和( 11）.然后，利用得到的封闭表达式，得到1DOF模型的色散曲线如图所示 3.． (10） ω ˜ 声 ＝ Δ − Φ Ψ ， （11） ω ˜ 光学 ＝ Δ + Φ Ψ ， 哪里 （12） Δ ＝ 米 + 我 ⌢ K ⌢ − 我 ⌢ K ⌢ 因为 质量保证 + 1 ， Ψ ＝ 我 ⌢ + 我 ⌢ 米 − 1 ， Φ ＝ 我 ⌢ 2 K ⌢ 2 因为 2 质量保证 + − 2 我 ⌢ 2 K ⌢ 2 + 4 我 ⌢ K ⌢ − 4 K ⌢ 因为 质量保证 + 1 + 4 K ⌢ + 2 米 − 2 我 ⌢ K ⌢ − 我 ⌢ K ⌢ 米 ．

一自由度模型的色散曲线。

方程( 10)和( 11)分别表示声学和光学频率分支。低频和高频色散曲线由 ω ˜ 声 和 ω ˜ 光学 分别为( 1].在上述方程式中，在 质量保证 ＝ ± π ，色散曲线的斜率，即行波的速度，等于零。这意味着波不能传播，就像在带隙处发生的情况一样。这个区域是第一个布里渊区边界。因此,用 质量保证 ＝ π 和 质量保证 ＝ 0 方程( 10)和( 11)然后求解这些表达式，提取出1DOF模型带隙的起始频率和终止频率[ 1， 20］.

 2．2．EMM的二自由度模型 2.2.1. 二自由度模型：数学公式

本节给出了EMM的二自由度模型(见图) 4)。此模型是上一节讨论的第一个模型的扩展模型。该模型的单元由矩形框架构成，框架内有一个齿条和齿轮，与一个集中质量相连。如图所示 5，小齿轮和集中质量通过刚度为的线性弹簧连接到矩形框架上 k 2 . 此外，矩形框架通过刚度为的线性弹簧相互连接 k 1 并创建一个完整的元结构。因此，单元的控制方程可以描述如下： （13) ∑ F x ＝ 米 1 x ¨ n ⇒ − k 1 x n − x n − 1 − k 1 x n − x n + 1 − 米 2 x ¨ n + r θ ¨ − 米 3. x ¨ n + x ¨ ＝ 米 1 x ¨ n ， （14） ∑ 米 O ＝ 我 O θ ¨ ⇒ − 2 k 2 r θ r + k 2 x r − 米 2 x ¨ n r ＝ 我 O θ ¨ ， （15） ∑ F x ＝ 米 3. x ¨ ⇒ − 米 3. x ¨ n − 2 k 2 x + k 2 r θ ＝ 米 3. x ¨ ， 哪里 米 3. 是集中的质量吗 x 表示集中质量的平移位移。

二自由度模型的超材料排列示意图。

二维自由度模型的单元示意图。

现在,方程( 13) ( 15)以矩阵形式编写，如下所示： （16） 米 1 + 米 2 + 米 3. 米 2 r 米 3. 米 2 r 我 O 0 米 3. 0 米 3. x ¨ n θ ¨ x ¨ + 2 k 1 0 0 0 2 k 2 r 2 − k 2 r 0 − k 2 r 2 k 2 x n θ x + − k 1 x n − 1 + x n + 1 0 0 ＝ 0 0 0 ．

为了以无量纲形式表示问题，为2DOF模型定义了以下无量纲参数： (17） 米 ¯ 1 ＝ 米 2 + 米 3. 米 1 ， 米 ¯ 2 ＝ 米 2 米 1 ， 我 ¯ ＝ 我 O 米 2 r 2 ， K ¯ ＝ k 2 k 1 ， ω n 1 ＝ k 1 米 1 ω n 2 ＝ 2 k 2 米 3. ， r x ¯ n ＝ x n ， r x ¯ ＝ x ， τ ＝ ω t ， ω ¯ 1 ＝ ω ω n 1 ， ω ¯ 2 ＝ ω ω n 2 ， 哪里 ω n 1 和 ω n 2 分别为第一模态和第二模态的固有频率， 米 ¯ 1 和 米 ¯ 2 分别定义为小齿轮总质量和集中质量与矩形框架质量之比和小齿轮与框架质量之比， 我 ¯ ， K ¯ ， x ¯ n ， x ¯ , τ 是否无量纲惯性矩、刚度、单元位移、集中质量位移和无量纲时间 ω ¯ 1 和 ω ¯ 2 分别表示第一个和第二个无量纲频率。使用引入的无量纲参数，方程( 16)可以以无因次形式重写为: （18） 1 + 米 ¯ 1 ω ¯ 1 2 米 ¯ 2 ω ¯ 1 2 米 ¯ 1 − 米 ¯ 2 ω ¯ 1 2 ω ¯ 1 2 我 ¯ ω ¯ 1 2 0 ω ¯ 2 2 0 ω ¯ 2 2 d 2 x ¯ n d τ 2 d 2 θ d τ 2 d 2 x ¯ d τ 2 + 2 0 0 0 2 K ¯ 米 ¯ 2 − K ¯ 米 ¯ 2 0 − 1 2 1 x ¯ n θ x ¯ + − x ¯ n − 1 + x ¯ n + 1 0 0 ＝ 0 0 0 ．

2.2.2.两自由度模型：色散特性

在本节中，为提取色散关系和带隙频率，给出了2DOF模型谐波对应的解如下: (19) x ¯ n θ x ¯ ＝ x ¯ n ∘ θ ∘ x ¯ ∘ e 我 质量保证 ． n + τ + c ． c ， x ¯ n ± 1 ＝ x ¯ n ± 1 ∘ e 我 质量保证 ． n + τ + c ． c ．

与上一节类似， x ¯ n ∘ ， θ ∘ , x ¯ ∘ 表示解的稳态振幅。代入方程( 19) ( 18)，给出了用矩阵符号表示的控制方程的新形式

通过简化方程( 1)时，控制方程的最终形式得到 (20） 1 + 米 ¯ 1 ω ¯ 1 2 米 ¯ 2 ω ¯ 1 2 米 ¯ 1 − 米 ¯ 2 ω ¯ 1 2 ω ¯ 1 2 我 ¯ ω ¯ 1 2 0 ω ¯ 2 2 0 ω ¯ 2 2 x ¯ n ∘ − 1 e 我 质量保证 ． n + τ + c ． c θ ∘ − 1 e 我 质量保证 ． n + τ + c ． c x ¯ ∘ − 1 e 我 质量保证 ． n + τ + c ． c + 2 0 0 0 2 K ¯ 米 ¯ 2 − K ¯ 米 ¯ 2 0 − 1 2 1 x ¯ n ∘ e 我 质量保证 ． n + τ + c ． c θ ∘ e 我 质量保证 ． n + τ + c ． c x ¯ ∘ e 我 质量保证 ． n + τ + c ． c + − x ¯ n ∘ e 我 质量保证 n − 1 e 我 τ + x ¯ n ∘ e 我 质量保证 n + 1 e 我 τ + c ． c 0 0 ＝ 0 0 0 ． (21） − 1 + 米 ¯ 1 ω ¯ 1 2 + 2 − 2 因为 质量保证 − 米 ¯ 2 ω ¯ 1 2 − 米 ¯ 1 − 米 ¯ 2 ω ¯ 1 2 − ω ¯ 1 2 − 我 ¯ ω ¯ 1 2 + 2 K ¯ 米 ¯ 2 − K ¯ 米 ¯ 2 − ω ¯ 2 2 − 1 2 1 − ω ¯ 2 2 x ¯ n ∘ θ ∘ x ¯ ∘ e 我 质量保证 ． n + τ + c ． c ＝ 0 0 0 ．

二自由度模型的色散关系可由( 21）.对于非零解，系数矩阵的行列式应设为零。因此，我们将色散关系提取为: (22） 一个 ω 6 + B ω n 1 2 + C ω n 2 2 ω 4 + D ω n 1 2 ω n 2 2 + E ω n 1 4 ω 2 + F ＝ 0 ， 哪里 （23) 一个 ＝ 2 米 ¯ 2 米 ¯ 2 我 ¯ − 1 − 我 ¯ ， B ＝ − 2 K ¯ 2 + 米 ¯ 2 + 4 米 ¯ 2 我 ¯ 1 − 因为 质量保证 ， C ＝ 米 ¯ 2 米 ¯ 2 + 2 我 ¯ + 米 ¯ 1 米 ¯ 2 1 + 2 我 ¯ ， D ＝ 3. K ¯ 1 + 米 ¯ 1 − 4 米 ¯ 2 我 ¯ 1 − 因为 质量保证 ， E ＝ 8 K ¯ 1 − 因为 质量保证 ， F ＝ − 6 K ¯ 1 − 因为 质量保证 ω n 1 4 ω n 2 2 ．

与单自由度模型的解释类似，式( 22)显示了两者之间的关系 ω 和qa。通过绘制这个方程，就可以得到色散曲线。通过替换 质量保证 ＝ π 在方程( 22)，则起动频率表达式为: (24） F 年代 1 ＝ H 3. 一个 − 1 2 C B + B ， F 年代 2 ＝ H 3. 一个 + C B + B ．

将qa = 0代入式( 22)时，系统参数与结束频率的相关性为 （25) F e 1 ＝ K − l 2 一个 ， F e 2 ＝ K + l 2 一个 ， 哪里 （26) 一个 ＝ 2 米 ¯ 2 2 我 ¯ + 2 米 ¯ 2 2 + 2 米 ¯ 2 我 ¯ ， B ＝ E + E + D − F 2 − C 3. 1 / 2 + D − F 1 / 3. ， C ＝ 1 3. 一个 H 2 3. 一个 + G ， D ＝ G H 6 一个 2 ， E ＝ H 3. 27 一个 3. ， F ＝ 6 K ω n 1 4 ω n 2 2 一个 ， G ＝ 16 K ω n 1 4 + 3. K + 3. K 米 ¯ 1 − 8 米 ¯ 2 我 ¯ ω n 1 2 ω n 2 2 ， H ＝ 米 ¯ 2 ω n 2 2 米 ¯ 2 + 米 ¯ 1 + 2 我 ¯ + 2 米 ¯ 1 我 ¯ + ω n 1 2 − 4 K ¯ − 2 K ¯ 米 ¯ 2 + 8 米 ¯ 2 我 ¯ ， 我 ＝ 米 ¯ 1 2 米 ¯ 2 2 2 我 ¯ + 1 2 + 2 米 ¯ 1 米 ¯ 2 3. + 4 米 ¯ 1 米 ¯ 2 2 我 ¯ 2 我 ¯ + 1 + 2 我 ¯ 米 ¯ 2 + 米 ¯ 2 2 2 ， J ＝ 4 K ¯ 米 ¯ 1 米 ¯ 2 2 4 我 ¯ + 5 + 8 K ¯ 米 ¯ 1 米 ¯ 2 我 ¯ − 1 + 16 K ¯ 米 ¯ 2 2 我 ¯ + 1 + 8 K ¯ 米 ¯ 2 我 ¯ − 4 K ¯ 米 ¯ 2 3. ， K ＝ H − 8 米 ¯ 2 我 ¯ ω n 1 2 ， l ＝ 16 K ¯ 米 ¯ 2 2 4 + 米 ¯ 2 + 1 ω n 1 4 + 我 ω n 2 4 + J ω n 1 2 ω n 2 2 ．

3.一自由度和二自由度模型的完整性

为了验证数学建模提取的研究结果，在MSC-ADAMS软件中，在刚体动力学假设下对弹性超材料进行建模和仿真，如图所示 6- - - - - - 9．此外，对模型进行了频响分析，并绘制了频响图。在下一节中，系统的瞬态参数用参数表示 Π 以及第5单元位移与第1单元位移的比值，然后与数值解进行比较。

一自由度模型的EMM布置示意图- msc - adams软件。

单自由度单元模型的原理图- msc - adams软件。

二自由度模型的EMM布置示意图- msc - adams软件。

二自由度模型的单元示意图- msc - adams软件。

采用MATLAB/Simulink软件进行数值仿真，并对所建模型的控制方程进行数值求解。用MATLAB/Simulink模型绘制的框图如图所示 10．值得一提的是，图中的每个部分都包含控制单个单元格动态的方程。

用MATLAB/Simulink绘制控制方程的数值解框图。(一)1自由度模型。(b) 2DOF模型，单元阵列(上)和单元(下)。

对于这两个模型，数值解和仿真得到的频响图如图所示 11和 12．如上所述，两种分析提供了相同的系统带隙间隔，证明了两种方法的结果令人满意的一致性，这表明我们的数学建模和数值计算的完整性。

FRA为一自由度模型，数值求解(MATLAB/Simulink)与仿真(MSC-ADAMS)比较。

二维自由度模型的FRA，数值解（MATLAB/Simulink）与仿真（MSC-ADAMS）的比较。

 4.结果与讨论

在验证数值解的可靠性后，全面研究了色散特性以及不同系统参数对带隙间隔的影响。在本节中，它的目的是检查引入的模型的带隙，并演示在哪些条件下它们可以扩展或传输。

4．1.EMM的一自由度模型 以下4.4.1。一自由度模型的色散曲线

身材 2显示了单自由度模型的色散曲线，其中显示了声波和光波模式以及带隙间隔。所选参数的值为 米 ＝ 1.56 ， 我 ⌢ ＝ 1．5 ， K ⌢ ＝ 0．35 ， 和 ω ⌢ n ＝ 8.9。 考虑这些参数，单自由度模型的带隙位于频率6.15 ~ 11.94 (Hz)之间，与图中得到的频率范围相同 11．

 4.1.2。一自由度模型中系统参数的影响

如图所示的表面 13展示了不同的无因次参数对其相应带隙的起始和结束频率的影响。顶面和底面分别表示开始频率和结束频率，它们之间的距离表示1DOF模型的带隙间隔。如图所示 13，带隙与无量纲刚度有直接关系，即随着无量纲刚度的增大，带隙也增大。另一方面，随着无量纲质量的增加，带隙减小。推导出刚度参数对带隙的扩展有积极的影响，质量参数的影响是削弱带隙的间隔。此外，从图中可以看出，一般情况下，无量纲惯性矩直接影响带隙范围，即随着无量纲惯性矩的增大，模型的带隙也会增大。虽然带隙间隔在较低值处略有减小，但当无量纲惯性矩取较高值时，带隙间隔得到了令人满意的增强。

改变的影响(a) K ⌢ 和 米 和(b) K ⌢ 和 我 ⌢ 对1DOF模型的带隙的起始和结束频率进行了分析。

 4．2．EMM的二自由度模型 4.2.1。准备二自由度模型的色散曲线

为了显示带有多个谐振器的修正模型的色散曲线，可解式( 22)为了 ω 并提取所需的曲线。所选参数的值为 米 ¯ 1 ＝ 1.04 ， 米 ¯ 2 ＝ 0.64 ， 我 ¯ ＝ 1．5 ， 和 K ¯ ＝ 2.86 ， ω n 1 ＝ 4．2 rad / 年代 ， ω n 2 ＝ 15.9 rad / 年代 ． 考虑以上参数，带隙启动频率( F 年代 1 和 F 年代 2 )通过代入计算 质量保证 ＝ ± π 分别在声波和光波模式下。结果图如图所示 14分别为Ds(1)和Ds(2)。此外，带隙结束频率( F e 1 和 F e 2 )也是通过替换来提取的 质量保证 ＝ 0 式( 22),分别。绘制的结果如图所示 14分别为Ds(2)和Ds(3)。如图所示 14，所述第一起始频率与第一结束频率之差表示第一带隙，所述第二起始频率与第二结束频率之差表示第二带隙。根据图中所示结果 12，可以发现在第一和第二带隙之间有一个小区域。这一结果也在图中得到了证明 14Ds(2)的最高点和最低点略有差异。微小的差异导致Ds(2)曲线变平。

二自由度模型的色散曲线。

 4.2.2。两种模型色散曲线的比较

如图所示 15，通过增加集中质量并将单自由度模型修改为双自由度模型，在新模型的动力学中产生了两个带隙，这使我们能够在更多的频率范围内减少或抑制改进系统的振动。对比两种模型的色散曲线，可以得出二自由度模型的第一个带隙的起始和结束频率都比单自由度模型的低。

两种模型中色散曲线的比较。

 4.2.3。二自由度模型中系统参数的影响

身材 16当参数 米 ¯ 1 向上移动，第一个起始频率也增加，第一个结束频率降低，这使修正模型的第一个带隙得到了满意的扩展。另一方面，通过增加第一个带隙，第二个带隙减小，起始和结束频率最终达到一个恒定值。

的影响 米 ¯ 1 参数对二自由度模型的带隙。

随着参数的增加 K ¯ ，第一个带隙可以扩展到更大的频率范围，如图所示 17，尽管第一个启动频率最终将接近一个恒定值。但是，模型的第二个带隙的变化并没有表现出相同的趋势，最终收敛到一个几乎恒定的值，变化不再明显。

的影响 K ¯ 参数对二自由度模型的带隙。

为了扩展修正模型的第一和第二带隙，对惯性矩进行了无量纲化 我 ¯ 必须增加。如图所示 18，继续增加无量纲转动惯量的过程，第一次结束频率和第二次开始频率不变，取恒定值。

的影响 我 ¯ 参数对二自由度模型的带隙。

随着第一和第二局域频率的增加，该模型的第一带隙并没有出现同样的行为，但其值较小时减小，且在较高时呈增加趋势(见图) 19和 20）.一般情况下，可以推断第一个带隙是展开的。最后，发现随着第一局域频率的增加，第二带隙的扩展并不均匀，但总的来说会变大。需要指出的是，当第二局域频率增大时，修正模型的第二带隙减小。

的影响 ω n 1 参数对二自由度模型的带隙。

的影响 ω n 2 参数对二自由度模型的带隙。

 5.结束语

本文首先提出了一种新的涉及齿轮齿条机构的弹性超材料模型(1DOF模型)。然后，通过增加集中质量，将模型转换为具有多个谐振腔的新元结构，其中产生更多的宽带带隙。不同的系统参数允许我们通过改变它们来改变或扩大系统振动减少或消除的频率范围。不同系统参数对所考虑的超材料带隙间隔的影响总结如下:

通过增加集中质量并将模型由一自由度转换为二自由度，系统的动力学过程中产生了两个满意的带隙

在1DOF模型中，为了实现更宽的带隙，一般采用无量纲参数 我 ⌢ 和 K ⌢ 是增加还是无量纲质量参数 米 应该减少

在二自由度模型中，随着参数的增加，第一个带隙的宽度得到了提高 K ¯ ， 我 ¯ ， 米 ¯ 1 , ω n 1 和减少 ω n 2

在二自由度模型中，第二带隙的宽度通过增加来提高 我 ¯ ， ω n 1 , ω n 2 或者减小参数 米 ¯ 1 . 通过改变参数，未显示出显著效果 K ¯

数据可用性

这些数据将在通讯作者要求时提供。

 的利益冲突

作者声明他们没有利益冲突。

 致谢

H. M. Sedighi感谢阿瓦士沙希德·查姆兰大学研究委员会的财政支持。SCU.EM99.98)。

 高 M。 吴 Z。 温 Z。 含纯立方振荡器的有效负质量非线性声学超材料 土木工程进展 2018 2018 15 3081783 10.1155 / 2018/3081783 2 - s2.0 - 85055593995 Mendhe s E。 科斯塔 Y.P。 超材料的特性和应用 国际信息技术与管理杂志 2011 4 1 85 89 Valipour 一个。 Kargozarfard m . H。 Rakhshi M。 Yaghootian 一个。 Sedighi h . M。 超材料及其应用:综述 机械工程师学会学报，第L部分:材料学报:设计与应用 2021 10.1177 / 1464420721995858 周 X。 王 J。 王 R。 林 J。 相关参数对多腔多腔声学超材料带隙的影响 应用物理 2016 122 4 1 8 10.1007/s00339-016-9978-x 2-s2.0-84961575126 阿尔巴 H。 迪堡 D。 辛格 T。 Nouh M。 有限的惯性放大声学超材料的色散跃迁和极零特性 应用物理学杂志 2018 123 1 12 唱 美国。 为声学超材料的应用提供了一种产生局部谐振子的新方法 美国声学学会杂志 2017 142 黄 H.H。 宣儿 c . T。 黄 g . L。 声学超材料中的负有效质量密度 国际工程科学杂志 2009 47 4 610 617 10.1016 / j.ijengsci.2008.12.007 2 - s2.0 - 60449109154 Nouh M。 Aldraihem O。 巴兹 一个。 具有周期性局部共振的超材料结构 结构和生物系统的健康监测 2014 9064 1 11 Orta A.H。 Yilmaz C。 惯性放大诱导的声子带隙产生的顺应轴向旋转运动转换机构 声音与振动杂志 2019 439 329 343 10.1016 / j.jsv.2018.10.014 2 - s2.0 - 85054808025 Lazarov b S。 詹森 j·S。 附加非线性振荡器链中的低频带隙 非线性力学学报 2007 42 10 1186 1193 10.1016 / j.ijnonlinmec.2007.09.007 2 - s2.0 - 36148988057 王 K。 周 J。 徐 D。 欧阳 H。 利用几何非线性原理研究一维周期杆中纵波的低带隙 机械系统和信号处理 2019 124 664 678 10.1016 / j.ymssp.2019.02.008 2-s2.0-85061388064 陈 j·S。 简 我。 嵌入膜质结构的超材料梁的动力行为 应用力学学报 2017 84 12 1 7 10.1115/1.4038146 2 - s2.0 - 85032002700 Casalotti 一个。 El-Borgi 美国。 Lacarbonara W。 嵌入非线性减振器的超材料梁 非线性力学学报 2018 98 32 42 10.1016/j.ijnonlinmec.2017.10.002 2-s2.0-85033676200 周 X。 小君 W。 王 R。 带有周期性局部谐振子的网格结构的带隙 现代物理字母B 2017 31 25 黄 g . L。 太阳 c . T。 多谐振器中的带隙 振动与声学学报 2010 132 1 6 10.1115/1.4000784 2 - s2.0 - 77955215753 太阳 H。 杜 X。 拜 P.F。 超材料光束宽带吸振的理论与实验研究 智能材料系统与结构学报 2011 21 11 1085 1101 肖 美国。 王 T。 刘 T。 周 C。 姜 X。 张 J。 有源超材料与超器件:综述 物理学报D:应用物理 2020 53 50 503002 10.1088 / 1361 - 6463 /列线图 陈 Y Y。 巴恩哈特 m V。 陈 j·K。 胡 g·K。 太阳 c . T。 黄 g . L。 亚波长尺度宽带波衰减耗散弹性超材料 复合结构 2016 136 358 371 10.1016 / j.compstruct.2015.09.048 2 - s2.0 - 84946223202 李 Z。 王 C。 王 X。 基于平动共振的负质量和模量弹性超材料模型 国际固体与结构学报 2019 162 271 284 10.1016 / j.ijsolstr.2018.12.015 2 - s2.0 - 85058564231 设拉子 A.H。 Sedighi h . M。 带旋转谐振器的齿条和小齿轮激励的超材料宽带抑制振动的物理学 欧洲物理杂志附加 2020 135 3. 1 23 Comi C。 马里戈 J-J。 二维局域共振超材料带隙预测的均匀化方法和bloch-floquet理论 杂志的弹性 2020 139 1 61 90 10.1007 / s10659 - 019 - 09743 - x 2 - s2.0 - 85069688685 O 'Faolain l 用于光学互连的光子晶体腔 用于数据中心的光互连 2016 英国剑桥 瑞斯出版 121 156