《Journal of Low Frequency Noise, Vibration and Active Control》：Nonlinear response analysis of the Duffing-Van der Pol oscillator with external excitation using a modified Lindstedt-Poincaré method
编辑推荐：
Van der Pol振子以在无外力状态下产生自持极限环而著称；然而当施加外部周期力时，其动力学响应可能严重偏离这一自然节律。现有Lindstedt-Poincaré方法在强非线性、阻尼及外力并存时变得繁琐。为克服该局限，研究人员引入一种替代性改进Lindst
Van der Pol振子以在无外力状态下产生自持极限环而著称；然而当施加外部周期力时，其动力学响应可能严重偏离这一自然节律。现有Lindstedt-Poincaré方法在强非线性、阻尼及外力并存时变得繁琐。为克服该局限，研究人员引入一种替代性改进Lindstedt-Poincaré方法，以获取Duffing-Van der Pol振子的稳态响应。通过将近似结果与Runge-Kutta四阶法(RK4)和谐波平衡法(HBM)所得结果比较，获得了优异精度。所提出的改进Lindstedt-Poincaré方法提供了简洁的求解流程，并对应用科学与工程中广泛的非线性振动问题表现出普适适用性。
研究背景方面，Van der Pol振子是研究自激振荡、共振及电学、机械、生物等系统中非线性现象的基准模型。传统数值技术如Runge-Kutta四阶法(RK4)虽能计算近似解，但不提供显式通解；经典扰动法（如Lindstedt-Poincaré法、Krylov-Bogoliubov-Mitropolskii法、多尺度法）在强非线性下失效。Cheung等提出的改进Lindstedt-Poincaré法虽有效，但在处理带阻尼与外激励的Duffing-Van der Pol振子时，扰动计算冗长且伴随耦合非线性代数方程，分析代价大。为此，研究人员在《Journal of Low Frequency Noise, Vibration and Active Control》发表本研究，引入相展开时间变换，构建替代性改进Lindstedt-Poincaré列式，使各阶扰动系数以显式解析形式顺序确定，在宽阻尼与外激励参数范围内以低计算量维持高精度。
关键技术方法上，研究人员以含线性阻尼μ、非线性阻尼λ、三次刚度项及外激励E cos(νt)的Duffing-Van der Pol方程（式13）为对象；先回顾Cheung等以α=εω?/(ω?2+εω?)为重标度参数的改进Lindstedt-Poincaré框架（式6–12），指出其一阶、二阶近似已显繁琐；继而引入新变量τ=νt?φ（φ为相位参数）作相展开时间变换，将x(τ)、频率ν2、相位φ按ε幂级数展开（式26–28），代入重整化方程（式25）后逐阶配平ε幂次，得到x?″+x?=0及后续线性递推方程（式31–33），以显式解析形式顺序求解各阶系数，避免耦合代数方程组。样本来源为数学方程本身，无实验队列。
研究结果部分，按原文小节叙述如下：
- 2.
Modified Lindstedt-Poincaré method for Van der Pol oscillator一节中，研究人员重述Cheung等保守系统列式（式1–12），展示α参数化下ω2=(ω?2/(1?α))(1+δ?α2+δ?α3+…)的频率展开，以及x(τ)=x?+αx?+α2x?+…的位移展开；将其推广至受迫阻尼Duffing-Van der Pol方程（式13–16），导出α幂次配平后的x?、x?、x?线性方程（式17–19），由初值x?(0)=a?、x?(0)=0求解，获得一阶近似幅值R=√(A?2+B?2)（式22–24）。研究人员指出该路径下高阶近似急剧繁复，需替代流程。
在引入τ=νt?φ后（式25），研究人员将x、ν2、φ展为ε级数（式26–30），代回得零阶x?=a cos τ（式35），一阶方程（式32）右端含?x?3?μ?x?′+μ?x?′?μ?x?2x?′+E cos(τ+φ?)，通过消除长期项条件可解析定出ω?、φ?；二阶方程（式33）同理顺推ω?、φ?。该相展开改进法使所有扰动系数显式顺序确定，免除Cheung路径中的耦合非线性代数方程。
讨论与结论部分，研究人员将原改进法与新相展开改进法对照：前者在α参数化下虽收敛但每阶需解耦合代数系统，后者借相位自由度将频率修正与相位修正解耦，各ε阶仅解常系数线性方程。通过与RK4及HBM在宽μ、λ、E、ν范围内比对，稳态幅频与相频曲线吻合优良，证明该替代列式在强非线性受迫Duffing-Van der Pol振子上以更简流程获高精度的结论。研究人员强调，所提改进Lindstedt-Poincaré方法流程简洁、计算量低，对应用科学与工程中的广类非线性振荡问题具备可迁移的解析近似价值。
