Mathematica软件在非线性振动力学中的一些应用_2.pdf

第二章 Duffing方程 参看资料：《Mathematica 全书》（第四版）中文版，（美）沃尔夫雷姆著，赫孝良，周义仓译，西安：西安交通 大学出版社，2002.7 《摄动方法导引》中文版，（美）Ali H. Nayfeh 著 《Perturbation Method with Mathematica》 Ali H. Nayfeh and Char-Ming Chin,1999. Off @General::spell1D 为了使某些记号符号化，我们需要加载程序包 Needs @Notation`D 2.1 Duffing方程 许多单自由度保守系统的自由振动的控制方程的形式为 eq21a x @tD+f @x @tDD 0; f x @tD x≥@tD f @x @tDD x x 式中 是 的非线性函数。这里 是系统的加速度，而 是恢复力。设 0是系统的平衡位置，那么 x 是一个常数，因此 f @x D 0。为了使Mathematica 把 x 解释为常数，我们需要把它符号化。否则，如果我们 0 0 0 用表达式g @tD代替x ，Mathematica将用下标为0的 g[t]代替x0 ，而不是把它作为一个常数处理。例如，在表达式 x2 中用 代替 得到 +x0 Sin @tD x x2 +x0 ê. x −Sin @tD 2 Sin @t D +Sin @t D0 显然地，x 被处理为下标为0的变量x ，这是错误的。因此，我们需要符号化 x (enter and select x followed by 0 0 0 clicking Symbolize@†D from NotationPalette) Symbolize @x0 D; 然后， Mathematica 就把符号化的表达式x0 处理为一个常数。在表达式x2 +x0 中用 Sin @tD代替x 得到正确的表 达式 x2 + x0 ê. x −Sin @tD 2 x0 +Sin @t D 接下来，我们假设f 是一个关于x 在x x0 处的解析函数，并在x0 附近展开为Taylor级数 2 Mathematica软件在非线性振动力学中的一些应用_2.nb fexp Series @f @x @tDD, 8x @tD, x , 3D 0 1 1 f 2 f H3L 3 4 f @x0 D+f @x0 DHx @t D−x0 L+ @x0 DHx @t D−x0 L + @x0 DHx @t D−x0 L +O @x @t D−x0 D 2 6 式中保留到立方项。使用Normal 从fexp 中截断高阶项，我们得到 fpoly fexp êêNormal 1