刘恩彩 方鑫 温激鸿 郁殿龙
(国防科技大学, 装备综合保障技术重点实验室, 长沙410073)
(2019 年 7 月 15日收到; 2020 年 1 月 4日收到修改稿)
以典型的双稳态系统—屈曲梁结构为例, 基于等效模型, 结合解析、数值和实验手段, 研究了双稳态结构中的1/2次谐波共振特性、演化过程、参数调节规律及其对隔振特性的影响. 研究发现, 当非线性刚度系数或激励幅值增加到一定程度时, 系统会在一定带宽下产生显著的1/2次谐波共振; 随着激励幅值增加, 阻尼系统的1/2次谐波遵循“产生-增强-衰退-消失”的过程, 该过程对峰值频率和峰值传递率有重要影响; 适当提高非线性强度能有效改善双稳态结构隔振特性. 针对双稳态屈曲梁结构开展的实验验证了1/2次谐波特性和隔振特性变化规律.
非线性系统的独特性质, 如跳跃、分叉、混沌等, 为新型功能结构设计提供了新的思路与机理[1−5].双稳态系统是一种典型的非线性系统, 其具有两个稳定平衡状态和一个不稳定平衡状态, 稳定平衡点之间的刚度为负值. 当能量输入较大时, 系统产生在两个稳定平衡状态之间的跳变跃迁运动, 导致强非线性现象[6]. 应用双稳态结构设计的隔振系统具有高静态、低动态刚度特征, 可在保证高承载性能的情况下实现低频隔振[7], 是双稳态系统的重要应用场景; 双稳态振子也可以实现高效的能量捕获和振动能量采集[8−13]. 双稳态结构中存在次谐波共振, 次谐波能用于质量检测[14,15]和拓宽双稳态能量采集器的工作频带[16,17]. 目前, 国内外学者对次谐波行为开展了一系列的研究工作[18−20], 然而, 双稳态结构中1/2次谐波的产生机理、演化过程及其对隔振特性的影响尚待深入分析.
双稳态可利用磁铁、栓锁、铰链、弹簧、屈曲梁等结构产生[21], 其中, 屈曲梁结构简单、适应性强,是最具代表性的双稳态结构[22−26], 大量学者对其隔振特性开展了分析. 研究表明, 双稳态屈曲梁结构[27]的突弹跳变[28]、混沌响应[29]、驱动力的位置[25]和量值[30]等都会影响其隔振特性. 然而, 其中的次谐波共振及其对隔振特性的影响的研究尚不完善,参数影响规律有待揭示.
本文以典型的双稳态屈曲梁结构为对象, 基于解析、数值和实验方法系统研究双稳态系统中的1/2次谐波共振的产生机理、分岔过程、参数影响规律及其对隔振特性的影响. 研究结果为双稳态系统的设计与应用提供有益参考.
本文研究的对象如图1(a)所示, 主要部分为两个并列的双稳态圆弧. 双稳态圆弧可保证实际实验中加载配重后样件的稳定性. 该结构可视作两个非线性弹簧并联. 在结构振动的过程中, 两弹簧位移相同, 共同承受附加在结构上的质量. 所研究对象可简化为图1(b)所示的非线性弹簧振子模型,为两组相同的弹簧振子结构并联. 屈曲梁几何参数为: 跨距l= 35 mm, 厚度t= 1 mm, 高度h=5.25 mm, 宽度b= 56 mm. 双稳态屈曲梁结构试样使用3D打印的方法制备, 材料为TPU (热塑性聚氨酯弹性体橡胶), 其材料参数通过压缩实验测得:密度r= 916.7 kg/m3, 弹性模量E= 17.9396 MPa,泊松比n= 0.385, 总配重M为 0.93 Kg.
图1 (a) 双稳态屈曲梁试样; (b)弹簧振子模型Fig. 1. (a) Prototype of bistable buckling beam; (b) spring oscillator structure.
如图1所示, 外力p作用在屈曲梁结构中间位置, 该位置的位移为x. 每根屈曲梁都可等效为一个包含立方刚度非线性的双稳态系统, 但在理论解析和数值仿真分析过程中仅需取其中一个分析其动力学特性. 系统微分方程为
其中质量m=M/2,c为阻尼系数,k0为结构负刚度,kn为立方刚度系数,u为激励位移, 则系统运动的绝对位移为x+u. 根据 Jin 等[31]的研究, 屈曲梁所受外力与所发生位移的关系为
求解得到的参数值分别为kn= 0.28 N/mm3,k0=–7.5 N/mm. 在静态时, 系统处于其中的一个稳态平衡点. 当响应幅值非常微小时, 质点仅在该平衡点附近运动, 响应特性近似为线性的. 该平衡点附近的线性刚度为
本文采用谐波平衡法求解该系统在给定激励下响应的解析解. 因阻尼对系统动力学特性的定性规律影响较弱, 为简化分析, 在解析分析中未考虑其影响. 实际试样中存在阻尼, 后文通过实验与仿真对比, 通过曲线拟合确定阻尼参数c= 30 kg/s.设系统的激励为u=Ucos(wt), 其中U为激励幅值,w= 2πf为激励频率. 由于双稳态系统包含偏离0的响应项, 将系统的解设为x=A0+A1cos(wt).代入方程(1)得
令方程(5)的常数项和coswt项的系数为零, 得
给定U可通过方程组(6)求得A0,A1, 即可求得系统在激励u下的相对位移x. 因此, 结构的振动传递率为
图2 Simulink 数值仿真模型Fig. 2. Numerical simulation model in Simulink.
本文采用MATLAB软件的Simulink模块对系统进行数值仿真分析. 针对系统的微分方程(1)建立系统的数值仿真模型, 如图2所示. 系统的激励u为正弦波,y=x+u为系统响应的绝对位移. 利用该模型, 可求解系统在给定激励下响应的数值解.
给定U= 0.01 mm, 该系统的数值仿真结果与解析结果对比如图3所示. 实数解析结果具有3个分支, 双稳态系统的非线性使传递率曲线弯向左侧, 与软刚度非线性系统类似; 0—29 Hz 频段,解析解有两个根, 系统稳定状态下的响应对应的一般为幅值较低的解. 仿真结果与解析解的分支1和分支2一致性较好, 然而仿真结果中传递率在58 Hz处出现奇异峰值, 此频率恰好为系统主共振频率f0的2倍, 因此可能产生了1/2次谐波共振现象. 但是根据(5)式得到的解析结果并没有预测出该奇异峰值, 原因在于以上理论没有考虑次谐波响应, 所以预测精度较低.
图3 解析与数值分析结果对比Fig. 3. Comparison between analytical and numerical results.
为了更准确地分析双稳态系统的响应特性, 在解设定时考虑1/2次谐波成分, 将系统的周期解表示为x=A0+A12cos(wt/2) +A1cos(wt) 的形式, 其中A12cos(wt/2)为1/2次谐波成分. 代入方程(1)得
令方程(8)的常数项、cos(wt/2)项系数、cos(wt)项系数为 0, 得
给定U可通过方程组 (9)求得A0,A12,A1,即可求得系统在激励位移u下的相对位移x. 因此,考虑1/2次谐波共振时系统传递率的表达式变为
此时, 基波与1/2次谐波的解相互耦合.
无阻尼情况下, 依然取U= 0.01 mm, 对比结果如图4所示. 解析结果除了具有与图3相似的3个分支外, 还具有另外一条向左弯曲的分支4,且该分支实际上是由两条曲线构成. 分支3与分支 4在时源于相同的分岔, 当时, 分支4的两条曲线与分支2重合. 数值分析表明, 在0—57 Hz频率范围内, 真实响应对应于幅值较低的解, 与未考虑 1/2次谐波相同; 在 2f0= 58 Hz附近, 系统的响应发生突变对应的是分支2到分支4的跳变; 当激励频率远离2f0时, 解再次从分支4跳变到分支2. 解析分析的主共振峰和1/2次谐波共振导致的奇异峰恰好与数值仿真结果十分符合. 系统在58 Hz正弦激励下的响应频谱如图5所示, 可以看出, 系统频域响应曲线除在 58 Hz处有尖峰外, 在主共振29 Hz处也有一个明显的峰值, 甚至主共振的响应大于基波的响应. 这些分析都证明, 图3中2f0的峰值是1/2次谐波共振引起的, 1/2次谐波共振导致高频能量向低频转移(而不是低频向高频转移), 无阻尼时导致主共振能量急剧增加从而使系统在2f0激励下的总响应幅值传递率增加80 dB以上, 为系统隔振带来不利影响.因此, 在实际使用双稳态隔振系统时, 需要设计恰当的结构和参数, 避免系统产生次谐波共振现象.当系统存在强阻尼时, 高频向低频转移的能量能够被有效耗散.
图4 考虑 1/2次谐波时的解析-数值结果对比Fig. 4. Comparison between analytical and numerical results with considering 1/2 sub-harmonics.
图5 次谐波共振发生时的频谱Fig. 5. Spectra for sub-harmonic resonance.
上文所分析的传递率为总幅值的传递率. 为更好地分析基波与1/2次谐波的耦合行为, 定义基波传递率为激励频率处的响应幅值与激励幅值之比,即T1=A1/U; 定义 1/2次谐波传递率为 1/2激励频率处的响应幅值与激励幅值之比, 即T12=A12/U.
利用(9)式求得幅值变化引起的传递率分岔见图6(a), 图6(a)中的基波传递率曲线1和1/2次谐波传递率曲线突出呈现在图6(b)中, 图6(a)和图6(b)的Y轴分别为对数坐标和线性坐标. 图中粉色虚线所代表的T1值对应的T12= 0, 因实际系统存在次谐波响应, 此虚线值可忽略. 基波传递率有两条连续曲线(曲线1和曲线 2), 但是1/2次谐波的解仅有单个分支.U→ 0时, 曲线 1上的T1→ –50 dB, 但曲线 2上的T1起始于很大的值; 由于激励幅值很小时系统呈线性规律, 因此曲线1代表了小幅值激励下系统的真实解. 0
图6 解析求解的幅值变化对次谐波共振的影响(a)对数坐标; (b) 线性坐标Fig. 6. Analytical results of the influence of amplitude on sub-harmonic resonance: (a) The Y coordinate of panel is logarithmic; (b) the Y coordinate of panel is linear.
图7 (a)呈现了无阻尼情况下的数值仿真结果,可见基波传递率呈现的规律与分岔分析中图6(b)(即图6(a)中曲线 1)的规律恰好符合; 当U<0.3 mm时, 图6(b)解析解与图7(a)数值解呈现的T12变化规律一致; 但是当U> 0.3 mm 时,1/2次谐波共振的数值解出现波动, 单频激励下的响应频谱为连续谱, 主要的能量并不是集中在单频的1/2次谐波上, 而是集中在一个较宽的频带10—30 Hz内 (图7(a)内插图), 因此系统的响应不再是周期或拟周期的, 这是典型的混沌响应特征.
图7 幅值变化对次谐波共振影响的数值解(a)无阻尼结果; (b) 有阻尼结果Fig. 7. Numerical results of the impact of amplitude on subharmonic resonance: (a) Results without damping; (b) results with damping.
图7 (b)呈现了阻尼系数c= 30 kg/s时的数值仿真结果. 可见当系统中出现阻尼时, 基波传递率在 0 4 参数变化对系统隔振特性的影响
为了研究参数变化对双稳态系统隔振特性的影响规律, 本节分别阐明: 1)激励幅值U对频率偏移量的影响; 2)固定立方刚度系数kn, 研究负刚度系数k0变化对非线性系统隔振特性的影响; 3)固定负刚度系数k0, 研究立方刚度系数kn变化对系统隔振特性的影响. 得到参数变化对系统1/2次谐波共振现象和隔振特性的影响规律, 以此作为系统动力学特性优化的依据.
在以下分析中, 解析分析不考虑阻尼的影响(c=0),在数值分析中,c= 30 kg/s, 分析系统参数变化的影响时, 激励幅值为0.3 mm.
共振频率和共振峰峰值是双稳态系统振动特性的两个重要指标, 也是隔振应用中十分关注的两个变量. 前文分析表明, 在分析主共振时解析计算中可以不考虑次谐波的影响, 不同幅值下的数值和解析结果分别如图8中实线和点线所示. 在共振区域, 数值与解析结果趋势一致: 随激励幅值增加,共振峰峰值降低, 共振点的位置向低频移动, 峰值传递率下降. 数值仿真结果还表明, 在一定带宽的激励下 (图中为 55—58 Hz), 均有1/2次谐波共振产生.
当考虑阻尼效应时, 上文所述共振频率偏移量会显著增加. 数值分析表明(图9), 随着激励幅值增加, 其共振频率从 28 Hz降低到 23 Hz, 频率前移了17.8%, 峰值传递率也随之下降. 但是激励幅值增加到一定程度 (此系统为U> 1.05 mm之后), 系统的峰值频率(不一定是主共振)急剧跳变增大, 峰值传递率也跳变升高. 由图7(b)可知, 这一跳变点恰好对应于T12=T1, 峰值频率和传递率升高的原因在于高幅值的1/2次谐波改变了系统的本质特性.
图8 无阻尼条件下幅值变化对隔振特性的影响Fig. 8. Influence of amplitude on vibration isolation characteristics without damping.
图9 有阻尼条件下幅值变化对频率偏移(a)和共振峰峰值(b)的影响Fig. 9. Influences of amplitude on frequency shifting (a) and the peaks of harmonic resonance (b) with damping.
固定立方刚度系数kn= 0.2 N/mm3, 通过改变双稳态屈曲梁结构的跨距l、厚度t、宽度b, 改变系统的负刚度参数, 变化k0= –0.2—–80 N/mm,结果如图10所示. 解析和数值仿真结果均表明,随|k0|的增加, 系统共振点的位置向高频移动, 系统共振峰升高. 共振频率的变化规律与公式相近, 但值略大于
图10 固定 kn = 0.2 N/mm3, 负刚度系数 k0变化对系统隔振特性的影响(a)解析结果; (b)数值仿真结果Fig. 10. When kn = 0.2 N/mm3, influence of k0 on vibration isolation: (a) Analytical results; (b) numerical simulations.
固定负刚度参数k0= –7.5 N/mm, 通过改变双稳态屈曲梁结构的高度h、跨距l、厚度t、宽度b,改变系统的立方刚度系数, 变化kn= 0.01—4 N/mm3,结果如图11所示. 解析和数值结果均表明, 随着kn的增加, 系统共振峰值频率略向低频偏移, 共振峰峰值降低. 与参数k0理论分析结果相比, 参数kn变化引起的效果远小于参数k0变化引起的效果.从图11(b)还可以看出,kn的大小与1/2次谐波共振是否出现有关, 当kn增大到一定值(此系统参数下约为 2.6 N/mm3)时, 系统出现次谐波共振. 且传递率被1/2次谐波共振增强的频带在kn=2.6 N/mm3达到 52—56 Hz, 进一步增加kn还会拓展该频带, 表明强非线性下产生1/2次谐波共振现象.
综上可得, 激励幅值影响系统的主共振频率;负刚度系数k0的变化主要引起系统共振点位置和共振峰大小的变化; 立方刚度系数kn主要与系统次谐波共振的出现有关.
图11 固定 k0 = –7.5 N/mm, 立方刚度系数 kn 变化对系统隔振特性的影响(a)解析结果; (b)数值仿真结果Fig. 11. When k0 = –7.5 N/mm, influence of kn on vibration isolation: (a) Analytical results; (b) numerical simulations.
实验装置如图12所示, 通过功率放大器利用电磁激振器在屈曲梁底部施加激励, 利用激光扫描测振系统拾取结构的激励和响应, 利用电脑实现控制和数据采集等功能.
为研究并验证该结构的1/2次谐波共振现象,本文在实验中施加频率为55 Hz的正弦激励并研究系统的响应特性, 结果如图13所示. 图13(a)表明, 系统频域响应曲线除在 55 Hz处有尖峰外, 在接近27.5 Hz处也有一个明显的峰值. 由此可知,本文设计的双稳态结构确实可以产生显著的1/2次谐波共振响应, 时域响应曲线每2倍周期响应幅值被加强. 此外, 本文实验还表明, 1/2次谐波共振现象能在基波为40—60 Hz的宽带范围内观察到, 证明了1/2次谐波共振在较宽的带宽下均可产生.
图12 (a)实验示意图; (b)实验装置Fig. 12. (a) Experimental schematic diagram; (b) experimental setups.
图13 频率为 55 Hz的正弦激励信号下系统的响应(a) 频域响应; (b)时域响应Fig. 13. Response under sinusoidal excitation signal with frequency of 55 Hz: (a) Frequency domain; (b) time domain.
为从实验中验证激励幅值对系统次谐波共振的影响规律, 本文以 40 Hz为例开展实验研究, 结果如图14和图15所示.
图14呈现了激励幅值逐渐增加时基波和1/2次谐波的传递率. 随激励幅值的增加, 系统的基波传递率整体呈下降的趋势, 1/2次谐波传递率在激励幅值较小时为0, 在激励幅值为1.251 mm时向上跳变, 然后逐渐降低直至消失. 这一规律与图7呈现的特性完全一致, 证明了1/2次谐波共振的耦合特性以及次谐波的“产生-增强-衰退-消失”过程.
图15呈现了跳变前后的响应特性. 在跳变发生前 (见图15(a)激励幅值U= 1.251 mm), 响应信号中仅包含基波成分. 但只需将激励幅值进一步增加 (图15(b)U= 1.351 mm), 响应频谱 (图15(c))除在 40 Hz处有峰值外, 在 20和 60 Hz处也有明显的共振峰, 且 20 Hz的频谱值更高, 1/2次谐波成分对系统响应的影响最大, 此时, 其时域波形产生2倍周期的显著畸变, 如图15(d)所示.
图14 激励幅值变化对 1/2次谐波共振的影响Fig. 14. Influence of excitation amplitude on the 1/2 subharmonic resonance.
为验证激励幅值对双稳态屈曲梁结构隔振特性的影响, 实验中施加扫频激励信号得到结构在不同激励幅值下的隔振特性, 实验与数值仿真结果对照如图16所示, 实验结果与数值仿真结果符合得较好. 实验证明, 随激励幅值增加,双稳态屈曲梁结构的共振频率向低频移动且共振峰峰值降低. 因此, 适当提高非线性系统的激励幅值, 能有效改善系统的隔振特性.
图15 系统 40 Hz 处的次谐波共振现象(a) U = 1.251 mm 时的响应和激励频谱; (b), (c)U = 1.351 mm 时的响应和激励频谱;(d)U = 1.351 mm 时的时域波形Fig. 15. Sub-harmonic resonance phenomena at 40 Hz: (a) Response and excitation spectrum with U = 1.251 mm; (b), (c) response and excitation spectrum with U = 1.351 mm; (d) time-domain waveform with U = 1.351 mm.
图16 实验-数值仿真结果对比Fig. 16. Comparison between numerical and experimental results.
本文利用谐波平衡、数值积分和实验研究了典型双稳态系统—屈曲梁结构中的1/2次谐波共振、演化过程、参数调节规律及其对隔振的影响,得到如下结论.
非线性强度增加到一定程度时会使系统产生显著的1/2次谐波共振, 使二倍固有频率附近的幅值传递率增加, 降低隔振效果, 且1/2次谐波共振能在一定带宽的激励下产生. 激励幅值较小时, 无次谐波共振现象; 激励幅值增加到一定程度时次谐波共振传递率急速增加, 但继续增加激励幅值会使次谐波传递率降低. 在激励幅值变化的过程中, 次谐波共振条件下基波和次谐波传递率产生分岔: 无阻尼时1/2次谐波的产生伴随着基波传递率增强过程, 但1/2次谐波传递率达到峰值并不会引起基波传递率的突变; 有阻尼时, 随着幅值增加, 1/2次谐波并不会在产生后一直存在, 而是遵循着“产生-增强-衰退-消失”的过程. 在这一过程中, 系统的峰值频率和峰值传递率先降低, 但是1/2次谐波传递率大于基波传递率时就引起峰值频率和峰值传递率的跳变增加.
负刚度参数k0对系统主共振特性影响较大,在给定激励幅值下, 随着|k0|的增大, 系统共振峰向高频偏移且传递率升高; 非线性系数kn主要影响系统次谐波共振特性, 对系统共振峰的峰值和共振点频率影响较小. 实验验证了理论分析结果.
非线性系统的次谐波共振会影响系统隔振区的隔振效果, 使系统在该频段的隔振效果变差, 在实际使用中应避开次谐波共振传递率增强的幅值区间; 而当需要应用次谐波的幅值增强效应时, 需调节参数使系统运行在“增强”幅值区间. 本文研究所揭示的规律和特性为双稳态结构的设计和应用提供指导.
附录A系统(1)的两个重要参数
利用(2)式求解双稳态屈曲梁机械结构理论模型的参数, 将其化为普通形式
将零点位置化到原点
取屈曲梁双稳态结构的几何、材料参数, 代入(A2)式得x3项系数
x项系数