具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

在本发明的描述中，需要明确的是，术语“垂直”、“横向”、“纵向”、“前”、“后”、“左”、“右”、“上”、“下”、“水平”等指示方位或位置关系为基于附图所示的方位或位置关系，仅仅是为了便于描述本发明，而不是意味着所指的装置或元件必须具有特有的方位或位置，因此不能理解为对本发明的限制。

下述介绍提供了本发明的多个实施例，不同实施例之间可以替换或者合并组合，因此本发明也可认为包含所记载的相同和/或不同实施例的所有可能组合。因而，如果一个实施例包含特征A、B、C，另一个实施例包含特征B、D，那么本申请也应视为包含有A、B、C、D的一个或多个所有其他可能的组合的实施例，尽管该实施例可能并未在以下内容中有明确的文字记载。

本发明实施例提供一种压电声学超材料的主动调控方法，如图1所示，具体按照如下步骤实施：

S1，建立压电声学超材料的机电耦合模型；

S2，根据所述S1的机电耦合模型，在压电声学超材料产生负动态刚度的前提下，计算压电声学超材料的带隙范围和形成条件；

S3，根据所述S1的机电耦合模型，推导传递函数并对谐振频率进行补偿，在满足所述S2的带隙形成条件的前提下，可形成范围及深度可调控的带隙；

S4，将所述S3的传递函数数字化，实现压电声学超材料的主动调控。

这样，采用上述方法，本实施例针对有限长度和有限数量局域谐振单元的一维压电声学超材料，首先建立机电耦合模型，之后在压电声学超材料产生负动态刚度的前提下，计算压电声学超材料的带隙范围和形成条件，接着推导传递函数并对谐振频率进行补偿，最后将传递函数数字化，实现了数字可编程压电声学超材料。

本发明所提出的主动调控方法可以实时地、方便地根据目标频率进行参数配置，在低频范围内实现了明显的带隙和超过40dB的非传统超强振动衰减；并且，压电声学超材料是一个多输入多输出(MIMO)系统，本实施例提出的主动调控方法极大降低了MIMO系统的设计复杂度，为主动控制算法在压电声学超材料中的大规模应用搭建了桥梁。

在具体实施例中：

所述S1中建立压电声学超材料的机电耦合模型，具体为：

根据控制单元，得到外部输入电压作用下的复合梁结构的动力学方程：

其中，EI是复合梁结构在压电元件短路条件下的挠曲刚度，w(x,t)是梁在位置x与时间t处的横向挠度，m是复合梁结构单位长度的质量，θ为物理坐标下的机电耦合项，v_j(t)是第j对压电元件电极间的电势差，压电元件的编号为j＝1,2,...,S，每对压电元件从x＝x_j ^L开始，到x_j ^R结束，总长度Δx_j＝x_j ^R-x_j ^L，H(s)是Heaviside函数，f(x,t)是分布作用在复合梁上的单位长度横向力。

将复合梁结构的横向挠度在模态空间展开，并截取前N阶振型，则有：

其中，模态的编号为r＝1,2,...,N，η_r(t)是第r阶模态对应的模态坐标表达式，φ_r(x)是第r阶模态对应的模态振型的质量归一化本征函数。

将分布作用在复合梁上的单位长度横向力f(x,t)用集中力进行等效替换，从而简化振动问题，得到复合梁结构在x＝x_f处的模态激振力q_r(t)：

建立压电声学超材料的机电耦合模型，其中超材料结构的机电耦合方程经过解耦处理，第r阶模态响应可以表示为：

式中，H_r(s)是第r阶模态对应的模态坐标的拉普拉斯变换，Q_r(s)是复合梁结构模态激振力的拉普拉斯变换，s是拉普拉斯算子，ζ_r是复合梁结构的第r阶模态对应的机械阻尼比，ω_r是复合梁结构的第r阶固有频率，α是与机电耦合作用相关的无量纲参数，C_p是压电片的内部等效电容，β是电压放大电路的电压放大倍数；

Z₀(s)＝1/(Cs)，Z(s)＝Ls+R，其中C是电荷放大器的反馈电容，L、R分别是增强型分流电路中与压电片串联的电感和电阻。

更具体：

所述控制单元包括一对压电元件和增强型分流电路，一对所述压电元件和增强型分流电路连接。

所述增强型分流电路1包括串联设置的电荷放大电路2、电压放大电路3、电阻以及电感，一对所述压电元件中的一片压电元件连接电荷放大电路，另一片压电元件利用逆压电效应形成压电作动器。

一对压电元件和增强型分流电路共同构成局域谐振单元，如图2所示为增强型分流电路局域谐振单元电路原理图，所述局域谐振单元包含一对压电元件和增强型分流电路，一对所述压电元件中的一片压电元件连接电荷放大电路，作用相当于压电传感器，另一片利用逆压电效应形成压电作动器，每个压电元件可以表示为与其内部电容并联的电流源。

增强型分流电路1由电荷放大电路2、电压放大电路3和电感电阻串联电路组成，由增强型分流电路和压电元件组成的整个框架可视为典型的主动控制系统，谐振原理是主要的控制策略。

在具体实施例中：

压电元件按照如下表1参数进行选择：

表1

在具体实施例中：

复合梁结构基体材料选用不锈钢，八对压电贴片紧密贴附在基体梁的上下表面，压电贴片之间的距离相同，并且压电贴片对于复合梁结构长度方向上的对称轴对称分布。压电贴片选用PZT-5H，成对贴附的压电元件的极化方向相反，基体梁尺寸与结构参数如表2所示：

表2

所述S2中根据机电耦合模型和压电声学超材料产生负动态刚度，计算出压电声学超材料的带隙范围和形成条件，具体为：

带隙范围是：

带隙形成条件：

上式中，ω_t是电感电阻串联分流电路的谐振频率，满足ω_t ²＝1/(LC_p)。

由带隙范围、形成条件及参数对比实验可知，α和β的值越大，压电声学超材料形成的带隙宽度和深度增加；R越大，带隙的宽度和深度逐渐减小，当R值过大时，压电超材料不存在带隙。

具体地，令

将s＝jω带入上式，在频域范围内有

ω′_r ²＝a_r(ω)+b_r(ω)j

其中

如果a_r(ω)<0，则ω'_r ²的实部小于零，这可以看作是压电声学超材料在第r个模态下实现了等效负动态刚度，此时ω所在的频率范围就是带隙范围。同时b_r(ω)始终大于零，其作用可以理解成为压电声学超材料增加了等效阻尼。因此，增强型分流电路通过等效负动态刚度和等效阻尼的联合作用实现了振动抑制。

在其中一个实施例中：

所述S3中根据带隙范围和形成条件，推导传递函数并对谐振频率进行补偿，具体为：

S31，根据所述控制单元得到传递函数的一般形式为：

上式中，γ是控制输出的附加增益。

图3所示为数字可编程压电声学超材料控制单元原理图，所述数字可编程压电声学超材料由若干个周期性分布的压电超材料控制单元组成，每个控制单元都相同。通过对增强型分流电路的电路原理进行数字化等效，得到一种数字可编程压电声学超材料，其中所述控制单元的传递函数的一般形式如上所示；

S32，根据带隙范围和形成条件，对所述S31中的谐振频率ω_t进行频率补偿：

ω_t＝ω_c+Δω_c

上式中，Δω_c是补偿频率，ω_c是带隙的中心频率。

在利用传递函数进行主动调控时，给定目标频率，期望带隙的中心频率与目标频率重合，但是数字电路形成的局域谐振带隙的中心频率会低于目标频率，因此必须对传递函数中的谐振频率ω_t进行补偿。根据传递函数，C_p和C已知，调节控制输出的附加增益γ实质上等效于调节机电耦合强度α，调节γ、β和R，实现了对带隙的宽度和深度的调控。需要注意的是，控制系统达到动态平衡后，继续增大γ和β无法再加大数字控制器的控制输出电压，对性能不再有增强效果，而且会导致系统发散。

在其中一个实施例中：

所述S4中将传递函数数字化，实现压电声学超材料的主动调控，具体为：

首先对传递函数根据采样时间使用零阶保持器离散化，然后将离散化后的二阶传递函数转化为差分方程的形式并输入至数字控制器中，实现压电声学超材料的主动调控。

具体地：

通过合适的电荷放大电路就可以将电荷信号转化为适合ADC采集范围的电压信号，利用嵌入式程序完全实现控制器逻辑。程序逻辑计算过程是在微控制器内部实时完成的。将所述传递函数通过微控制器数字化，首先对连续时间传递函数根据采样时间使用零阶保持器离散化，然后将离散化后的二阶传递函数转化为差分方程的形式，最终通过相应的滤波器结构实现了数字化：

y(k)＝a₀x(k)+a₁x(k-1)+a₂x(k-2)-b₁y(k-1)-b₂y(k-2)

式中，x(k)是k时刻的输入信号，y(k)是k时刻的输出信号，a₀，a₁，a₂是前馈滤波器系数，b₁，b₂是反馈滤波器系数。

在本实施例中，采样频率选取10kHz。

在所有压电元件接地实现短路的条件下，对压电声学超材料进行扫频实验，由信号发生器提供扫频信号，频率从5Hz线性变化至200Hz，数据采集装置采集信号发生器输出的激振电压作为系统输入，采集加速度计输出的电压信号作为系统输出，根据输入输出数据，利用信号处理方法可以得到系统的幅频特性如图4所示。选取固有频率105.6Hz作为数字可编程压电声学超材料的目标频率，在其附近进行振动控制实验，可以看出超材料产生了明显的带隙，取得超过40dB的显著振动衰减。

数字控制器可以无线、实时地调整参数，从而轻易地控制带隙的频率、宽度和深度。本发明提出的主动调控方法极大降低了MIMO系统的设计复杂度，为主动控制算法在压电声学超材料中的大规模应用搭建了桥梁。

以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。