1. 引言
近年来研究人员发现非线性能量阱(NES)具有质量轻、振动抑制频带宽、可以实现靶能传递 [1] 等优点,对其进行了大量研究。Wierschem等人 [2] 研究了不同NES质量对主系统能量传递的影响,得到NES质量与主系统质量比远远小于1的结论。张也弛和孙斌等人 [3] [4] 分析了NES系统在宽频激励下主系统结构的幅频响应,结果表明NES系统具有优良的宽频减振效果。在此基础上张也弛等人 [5] 进一步研究了两自由度非线性耦合振子间的靶能量传递现象,并通过数值仿真验证了NES系统最优靶能量传递所需的初始能量。此外,NES的耦合系统相关性能也成为研究热点。Fang和Zhang等人 [6] [7] 将NES与超磁致伸缩材料相耦合实现了系统的振动控制和能量采集。Xiong等人 [8] 将NES与压电振动能量采集器(PVEH)相耦合,实现了系统宽频能量采集。Pennisi等人 [9] 将NES与磁电能量转换器相耦合,通过实验证明了耦合系统具有优良的振动抑制和能量转移特性。由于NES系统复杂的动力学特性,常需要运用一些近似解析的方法对系统进行数值仿真。Sun等人 [10] 结合复化平均法研究了NES系统对飞轮的微振动影响。Gendelman等人 [11] 运用平均–多尺度法研究了NES系统在简谐激励下的振动抑制特性。Zang和国忠金等人 [12] [13] 运用谐波平衡法研究了NES系统的动力学特性。这些研究方法都很好地反映了NES在系统减振方面的优良性能,为以后的科学实验提供理论基础。
2. 系统模型
如图1给出了NES-电磁能量采集耦合系统连接在一个线性主系统上的力学模型,以及电磁能量采集的电学示意图,由牛顿第二定律列出如下非线性控制方程:
(1)
其中,
、
和
分别为单自由度系统的质量、刚度和阻尼;
、
和
分别为NES的质量、非线性刚度和阻尼,
为单自由度系统质量的位移,
为NES质量的位移,
表示电磁能量采集器提供的力(电磁能量采集器的反电动势),
、
和
分别表示电机常数、电动势系数和内阻,位移激励u作用在基座上,可以表示为:
(2)
其中
是基础位移幅值,f是激励频率。
根据文献 [14] [15] [16] 可以将磁能量采集器的系统参数等效为一个电磁阻尼
,由相对位移关系和电磁能量采集的工作原理可以求得耦合系统平均电功率p,分别表示为:
(3)
(4)
其中
是主系统质量
位移幅值,
是NES系统质量
位移幅值。
3. 分析方法
由于系统是强非线性激励系统,通过谐波平衡法对方程(1)进行近似求解,假设方程(1)的谐波解是满足如下形式:
(5)
其中
和
表示用谐波平衡法时解的傅里叶系数,将方程(5)代入方程(1)然后对代数方程组进行推导和简化得到
、
、
和
系数分别相等的8个代数方程式。
(6)
在给定的激励频率下,通过牛顿迭代法和弧长延拓法进一步求解方程(6),最终可以求得控制方程的解。其中方程(6)的Jacobi矩阵如下:
(7)
根据文献 [7] [14] [15] 选取如表1所示的系统参数,将
代入控制方程研究一体化耦合系统在振动抑制和能量采集方面的性能。
图2为通过谐波平衡–弧长延拓法与4阶龙格库塔法分别求解系统幅频响应得到的结果对比图,分析结果表明解析解和数值解一致重合,证明了研究方法的正确性。
图3给出了当系统参数
,
时,结合谐波平衡–弧长延拓法和4阶龙格库塔法分析系统幅频响应的结果对比图,其解析解和数值解仍然重合,再次证明了研究方法的正确性。
4. 对比分析
4.1. 电磁能量采集器
图4给出了电磁能量采集器对主系统稳态位移幅频响应的影响结果对比图,通过对比NES-电磁能量采集耦合系统和单一NES系统的最值,结果表明耦合系统的位移幅值远小于NES系统的位移幅值。证明一体化耦合系统比NES具有更优良的减振性能。
![](//html.hanspub.org/file/2-2650141x46_hanspub.png)
Figure 4. Comparative analysis of the presence or absence of an electromagnetic energy harvester
图4. 有无电磁能量采集器的对比分析
4.2. 系统参数分析
图5至图8给出了各参数对系统位移幅频响应的影响结果对比图。图5 NES阻尼
和图6等效电磁阻尼
的参数变化对系统位移幅频响应的影响较大,而图7非线性刚度
和图8 NES质量
参数的改变对系统位移幅频响应无明显影响。结果表明阻尼的改变是影响系统位移幅值的主要原因,因此通过调节电磁能量采集器的电阻
可以改变等效电磁阻尼的大小,进而影响系统的振动,以达到所需要的振动效果。
![](//html.hanspub.org/file/2-2650141x52_hanspub.png)
Figure 5. NES damping effects on the frequency domain
图5. NES阻尼对频域影响
![](//html.hanspub.org/file/2-2650141x53_hanspub.png)
Figure 6. Effect of electromagnetic damping on the frequency domain
图6. 电磁阻尼对频域影响
![](//html.hanspub.org/file/2-2650141x54_hanspub.png)
Figure 7. Effect of nonlinear stiffness on the frequency domain
图7. 非线性刚度对频域影响
![](//html.hanspub.org/file/2-2650141x55_hanspub.png)
Figure 8. NES quality impact on the frequency domain
图8. NES质量对频域影响
4.3. 时域分析
图9对比了主系统质量
和NES质量
的位移关系,但这些位移关系与系统参数的选取有关,因此图10对比了电磁阻尼
分别为50、200、350 Ns/m时对主系统位移的影响,结果表明电磁阻尼参数越增对主系统振动抑制效果越好,因此通过增大电磁阻尼参数可以更好地抑制主系统的振动。
![](//html.hanspub.org/file/2-2650141x59_hanspub.png)
Figure 9. Displacement of the main system and the NES system
图9. 主系统和NES系统的位移
![](//html.hanspub.org/file/2-2650141x60_hanspub.png)
Figure 10. Effect of equivalent electromagnetic damping on the time domain
图10. 等效电磁阻尼对时域影响
4.4. 平均电磁功率和电压
图11为电磁阻尼
分别为50、200、350 Ns/m时,不同频率下平均电功率的输出关系对比图,结果表明当频率增大到20 Hz开始有功率的输出,30 Hz附近上升到最大功率,随后迅速下降;局部放大图表明功率随着电磁阻尼增大的变化,会呈现出先增大后减小的变化趋势,因此通过设定电磁阻尼参数可以在特定频率下输出更多功率。
由功率和耦合系统电阻的关系可以进一步分析系统的输出电压U,表示为:
(8)
图12给出了一体化耦合系统电磁阻尼
时输出电压和频率的对比关系,在15 Hz时有电压产生,在30 Hz处电压达到最大值然后快速下降并趋于稳定。局部放大图表明输出电压随着电磁阻尼的增大,会呈现出先增大后减小的变化趋势,因此通过设定电磁阻尼参数可以在特定频率下输出更高的电压。进一步研究表明在较大频率范围内,随着电磁阻尼的增大可以输出更高的电压,因此在较高频率下可以实现系统电压的稳定输出。
5. 能量关系
在系统受基础简谐位移激励时,采用4阶龙格库塔法对系统的能量吸收情况进行相关计算。耦合系统吸收的能量占输入总能量的百分比
表示如下 [7] :
(9)
NES系统吸收的能量占输入总能量的百分比
表示为:
(10)
基础系统吸收的能量占系统总能量的百分比
可以表示为:
(11)
一体化耦合系统的能量采集率
可以表示为:
(12)
其中
和
分别为NES-电磁能量采集系统的动能和势能,
为NES阻尼的能量耗散,
为电磁能量采集器从系统中吸收的能量,
为电磁能量采集器从系统中采集到的能量,
为系统输入的总能量。根据能量守恒原理系统的总能量
表示为:
(13)
其中
和
分别为基础结构响应的动能和势能,
为基础阻尼的耗散能。各能量计算如下 [7] :
(14)
(15)
(16)
(17)
(18)
(19)
(20)
(21)
对比分析图13和图14可知在激励开始时,基础系统中吸收的能量百分比急剧上升,而耦合系统中吸收的能量百分比急剧下降;随后主系统的能量吸收减少,而耦合系统的能量吸收增加。在30 s的时间内96.1%的能量分布在耦合系统,3.9%的能量分布在基础系统中,结果表明受到外激励时输入的能量大部分被耦合系统稳定吸收。通过基础系统的能量快速转移到耦合系统,可以更好地实现耦合系统的能量采集。
图15为耦合系统的能量采集百分比随时间的变化图。在激励开始时,耦合系统的能量采集百分比快速上升到0.78%,并且随后稳定在这一数值附近。结果表明电磁能量采集器可以将系统的振动能转化为电能输出,
![](//html.hanspub.org/file/2-2650141x93_hanspub.png)
Figure 13. Energy absorption rate of the main system
图13. 主系统能量吸收率
![](//html.hanspub.org/file/2-2650141x94_hanspub.png)
Figure 14. Energy absorption rate of the NES-electromagnetic energy acquisition coupling system
图14. NES-电磁能量采集耦合系统能量吸收率
![](//html.hanspub.org/file/2-2650141x95_hanspub.png)
Figure 15. Energy harvesting rate of coupled systems
图15. 耦合系统能量采集率
因此当耦合系统和主系统有相对运动时可以为外系统某些不易更换电池的小功率元件提供源源不断的电能。
6. 结论
本文研究NES-电磁能量采集一体化耦合系统的非线性动力学响应和能量采集。通过谐波平衡–弧长延拓法和四阶龙格库塔法验证分析结果的准确性。讨论有无电磁能量采集器对主系统位移幅值的影响;并对比分析系统参数对频域、时域和功率等的影响;进一步研究耦合系统能量吸收和采集的变化情况。
根据本文的研究,得到了以下结论:
1) 谐波平衡–弧长延拓法的近似解与控制方程的数值解高度一致,证明了本文采用的分析方法是正确的。
2) 对比分析主系统的幅频响应,结果表明NES-电磁能量采集耦合系统比单一NES系统,具有更优良的振动抑制性能。
3) 耦合系统平均电功率和电压随电磁阻尼的增大,在特定频率下会呈现先增大后减小的变化趋势,因此通过设定电磁阻尼参数可以输出更高功率和电压。进一步研究表明,耦合系统可以将系统振动的能量转化为电能输出,实现了对系统的能量采集。