质量变化与曲柄滑块机构几何对参数激励和DC激励下摆结构稳定性的耦合效应

《Mechanics Research Communications》：A note on the Coupled Effect of Mass Variation and the Slider-Crank Mechanism's Geometry on the Structural Stability of a Pendulum under Parametric and DC Excitation

　　

研究人员研究了质量变化与曲柄滑块机构几何对参数激励和DC激励下摆结构稳定性的耦合效应。分析了参数激励下的摆动力学，通过Floquet乘子（μ）构建了稳定性图，揭示了Flip分岔与Neimark-Sacker分岔，数值模拟显示了亚谐波周期倍化，并表征了通往混沌的转变路径。

论文解读：质量变化与曲柄滑块机构几何对参数激励摆结构稳定性的耦合影响

参数激励摆是非线性动力学研究中的经典模型，其垂直平衡位置的稳定性传统上由Mathieu方程和Hill方程描述。当摆的支撑点与曲柄滑块机构耦合时，运动学约束引入了多频参数调制，从根本上改变了控制方程的时变系数。这种机构在非线性能量采集系统和调谐质量阻尼器设计中有重要应用，例如利用环境激励驱动的参数摆式能量采集器，或用于被动抑制结构振动。已有研究基于纯无量纲参数模型，在归一化驱动频率与几何耦合比（ω, ε）参数空间内进行了广泛映射，但将物理项合并为抽象参数会掩盖工程变量之间的直接相互作用。实际机械系统中，改变摆端质量m或连杆长度比ε会同时影响多个无量纲系数。为此，研究人员采用物理参数化模型，显式保留质量m在方程中的作用，在固定几何ε的条件下考察(ω, m)稳定平面，以揭示质量变化与曲柄滑块几何的耦合效应对结构稳定性的影响。

研究人员通过欧拉-拉格朗日方程推导了单自由度运动方程，引入无量纲时间τ=ω₀t，得到无量纲阻尼比ζ=c/(mω₀)和参数激励幅值p=aω²/?。在状态空间中线性化后，采用基于修正第一类Chebyshev多项式和Picard迭代法的高精度数值格式，计算了单值矩阵及Floquet乘子（μ）。稳定性图在四个几何场景下扫掠质量0.010=5.7174rad/s，黏性阻尼系数c=0.0571N·s/m，连杆长度b=0.25m，摆长?=0.3m，对应ε=a/b分别为0.04、0.4、0.6和0.9。

4.1 场景一：准谐波机制（a=0.01m）
该场景下ε=0.04，参数激励函数接近于纯谐波。Floquet乘子分析显示，在(ω, m)平面上仅存在一个位于ω≈2的主不稳定舌区。沿m=1kg扫描时，边界处实乘子在μ=-1处穿出单位圆，对应Flip（周期倍化）分岔。舌区最低质量阈值为m=0.15kg，位于ω=2处。

4.2 场景二：非线性过渡机制（a=0.1m）
ε=0.4时，激励曲线出现两个局部最小值和中心局部最大值。稳定性图显示两个不稳定舌区：第一个出现在ω≈1附近，边界乘子穿过μ=+1，构成Fold（鞍结）分岔，带宽Δω=0.0375；第二个舌区覆盖主共振ω≈2，边界同样为μ=-1的Flip分岔，带宽大幅扩展。第一舌区最低质量阈值为m=0.26kg，第二舌区延伸至m=0.01kg。

4.3 场景三：新不稳定区域（a=0.15m）
ε=0.6时，激励函数呈明显多谐波“W”形。稳定性图新增一个位于ω≈0.66的中间舌区，边界乘子穿过μ=-1，属Flip分岔；第二个舌区（ω≈1）仍由Fold分岔支配，带宽增大至Δω=0.0868。新增舌区最低质量阈值为m=0.57kg，第二舌区最低质量阈值降至m=0.10kg。

4.4 场景四：接近几何极限（a=0.225m）
ε=0.9时，激励函数畸变加剧，所有不稳定舌区显著扩张。第一个中间舌区阈值带宽增至Δω≈0.05，仍为Flip分岔；第二个舌区（ω≈1）带宽大幅增至Δω≈0.22，为Fold分岔；且从ω=1.31056开始，乘子在μ=-1处穿出后不再回到单位圆内部，系统对高频激励呈现永久不稳定。两个中间舌区的最低质量阈值分别降至m=0.13kg和m=0.04kg。

4.5 动态机制综合与物理讨论
随着ε增大，激励函数的Fourier级数不再纯余弦，而是引入高次谐波。依据参数共振准则，不稳定区域出现在归一化频率接近2ω₀/n处（n=1,2,3），这与数值结果吻合。Flip分岔（μ=-1）导致亚谐波周期倍化；Fold分岔（μ=+1）则造成对称破缺和跳跃现象。因此，忽略高次谐波的经典Mathieu模型会严重低估结构失稳风险，尤其在较大ε时，即使低质量系统中因ζ=c/(mω₀)增大而具有高相对阻尼，多谐波能量注入仍足以破坏垂直平衡。

结论部分可翻译为：本工作研究了质量变化（m）与曲柄滑块几何（ε=a/b）对参数激励摆结构稳定性的耦合效应。采用基于欧拉-拉格朗日方程的物理参数化模型，在理想运动学激励θ(t)=Ωt下，质量明确决定无量纲阻尼比ζ=c/(mω₀)，几何参数则控制激励函数的幅值及多谐波轮廓。结合修正Chebyshev多项式与Picard迭代法的解析-数值方案，构造单值矩阵并评估Floquet乘子。在(ω, m)平面上的稳定性图揭示了Hill不稳定舌区的渐进拓扑演化：场景一（ε=0.04）呈准谐波性，只有主共振ω≈2处的单舌区；场景二（ε=0.4）激活几何非线性，在ω≈1处出现第二舌区并伴随主舌区持续加宽；场景三（ε=0.6）在ω≈0.66处出现第三亚谐波舌区；极限场景四（ε=0.9）中，严重非线性耦合使各不稳定区域显著放大。物理上，边界处乘子穿过μ=-1（Flip/周期倍化分岔）引起稳定的亚谐波振荡，而ω≈1处的边界乘子穿过μ=+1（Fold/循环鞍结分岔），导致对称破缺和向大幅振荡或持续旋转的跳跃。结果表明，经典Mathieu类公式（ε→0）忽视有限曲柄滑块比激活的高次谐波，会严重低估结构风险；在ε=0.9时临界质量阈值剧烈降低，即便轻质高阻尼系统也因耦合机构的强能量注入而易于失稳。上述现象与Hill方程的一般理论一致，为理想激励下曲柄滑块摆系统提供了安全设计图，并为考虑驱动电机功率约束和摆反作用力矩的非理想系统研究奠定了基础。