非线性二自由度范德波尔振动：基于双时间域法的El-Dib优化频率公式

《Journal of Low Frequency Noise, Vibration and Active Control》：Nonlinear 2DOF van der Pol vibrations: El-Dib’s optimized frequency formula via the two-time-domains method

　　

本文提出一种针对强耦合二自由度（2-DOF）非线性振动系统的系统性解耦框架。此处耦合被归类为强耦合而非弱耦合，因为控制方程中的跨模态相互作用项（如x2y型、xy2型和xy型项）与主导的立方刚度非线性项同阶，无小参数对其进行缩放，而弱耦合（摄动）表述中会假设存在此类小参数。所提方法通过双时间域Lindstedt-Poincaré型变换，为每个模态分配各自独立伸缩的时间域ξ = Ω1t与ζ = Ω2t。结合模态特异性加权平均，将控制方程解耦为两个解析可处理的子系统，每个子系统携带各自的闭式、幅值依赖频率。该方法将El-Dib频率公式拓展至双项耦合构型，产生准确刻画非线性系统动态的闭式频率-幅值关系。解析解在幅值、相位与振荡包络方面与直接数值模拟高度吻合；残差差异仅归因于高阶截断效应。分析表明，模态间耦合与模态阻尼均显著影响振动模态间能量交换及系统整体稳定性。所提方法可直接应用于流固耦合、旋转机械、微机电系统（MEMS）谐振器与磁流体界面中出现的若干非线性问题，从而提供一个连接非线性动力学理论与定量工程预测的严谨解析平台。

研究背景方面，耦合非线性多自由度系统的解析处理长期依赖多尺度法等摄动技术，但在强耦合情形下跨模态项（如x²y、xy²）与立方刚度同阶，无小参数可缩放，经典一阶摄动难以捕捉内共振与模态分岔；同时单时间轴表述在ω₁₀≠ω₂₀时会引入伪共振与人工耦合。因此研究人员开展本研究，旨在建立一种非摄动、系统性解耦强耦合2-DOF非线性振动（含van der Pol型阻尼与Duffing型刚度）的解析框架，并拓展El-Dib频率公式至耦合构型。论文发表于《Journal of Low Frequency Noise, Vibration and Active Control》。

关键技术方法上，研究人员采用双时间域Lindstedt-Poincaré变换，为x、y两模态分别引入独立伸缩时间ξ=Ω₁t、ζ=Ω₂t；将耦合非线性刚度写为奇乘偶形式F₁=x f₁(x,y)、F₂=y f₂(x,y)；构造试解x₀=A cos Ω₁t、y₀=B cos Ω₂t，定义模态权函数W₁(ξ)、W₂(ζ)；通过均方误差最小化导出扩展El-Dib保守频率公式（式14、15），以加权平均等效线性化替代立方项，获得闭式幅值依赖频率Ω₁、Ω₂，全程无小参数展开。

研究结果部分，按原文小节叙述如下。第2节“Significance of the problem”中，研究人员由Newton第二定律与Euler-Lagrange方程推导含线性粘性阻尼、立方自刚度及x²y、xy²、?²x等交叉项的2-DOF van der Pol–Duffing方程（式1、2），明确质量归一化系数μ₁、μ₂、ω₁₀²、ω₂₀²及α_i、β_j的物理来源，指出强耦合源于势能V_c(x,y)四次多项式微分所得交叉力。第3节“Decoupling strategy for interconnected nonlinear two-degree-of-freedom vibrations”中，研究人员论证单时间变量仅在ω₁₀=ω₂₀或简并内共振时合法，否则须赋各模态独立时间域；经变换得式12、13，代数耦合项仅作独立变量代换不引入交叉导数，两方程干净分离。第3.1节“Two-time domains (Lindstedt–Poincaré) transformation”中，确认Ω₁、Ω₂为非摄动闭式幅值函数，由等效线性化直接给出。第4节“Theorem”中，研究人员给出扩展El-Dib公式（式14、15），以W₁²W₂²加权双积分比定义ω₁₁²、ω₂₂²，经均方残差（式18–20）置零导出解耦线性方程（式16、17），每模态携自身频率关系F(Ω_i)=0。

讨论与结论部分，研究人员指出：双时间域法消除了单时间轴对异频模态的人为频率共享，避免伪共振；扩展El-Dib公式在强耦合下仍给出闭式幅值依赖频率，解析解与数值模拟在幅值、相位、包络上吻合，残差仅来自高阶截断；跨模态耦合系数与模态阻尼μ₁、μ₂共同调控能量转移与稳定性。该框架可直推至流固耦合、旋转机械、MEMS谐振器及磁流体界面等工程场景，为非线性动力学理论至定量预测提供解析桥梁。结论原文要义为：所提双时间域解耦与加权等效线性化将El-Dib频率公式成功推广至强耦合2-DOF van der Pol–Duffing系统，获得闭式、幅值依赖、非摄动频率，验证精度由数值对照确认，方法具备跨工程领域的直接适用性。