各解析逼近周期解与.pptVIP

吉林大学博士学位论文答辩会 答 辩 人：李鹏松 指导教师：吴柏生 教授 Mickens et al. 应用有理谐波平衡法建立了奇非线性一维保守系统振动问题的解析逼近周期及逼近周期解。Sarma et al.将Mickens的方法进一步扩展，建立了恢复力为一般函数的一维保守系统非线性振动问题的解析逼近周期及周期解。上述有理谐波平衡法在构造非线性振动问题的高阶解析逼近解时，由于涉及到求复杂的非线性代数方程组解析解，求解是很困难的。 二.求解大振幅非线性振动问题的修正谐波平衡法 例1. Duffing方程 例2. Duffing-harmonic振子 三. 求解大振幅非线性振动问题的迭代法 图2.6各解析逼近周期解与“精确”周期解的比较 图2.7 各解析逼近周期解与“精确”周期解的比较 3.非自然系统 通过以上的讨论可以看到，本文提出的修正谐波平衡法对于恢复力为位移奇函数、恢复力为位移一般函数、Duffing-harmonic 振子及非自然系统等强非线性振动问题都能给出高精度解析逼近周期与周期解。 考虑如下非线性振动方程 (2.13) 的情形，此时修正谐波平衡法依然能给 出高精度解析逼近周期与周期解。 本文我们只考虑 考虑如下的一维保守系统非线性振动方程 (3.1) 引入新变量 ，将方程(3.1)写为 (3.2) 本文提出如下的迭代过程 (3.3a) (3.3b) 1. 非线性恢复力为位移奇函数的情形 此时初始输入函数可取为 (3.4) 对于给定的 及 ，式 (3.3a) 是关于 的二阶非齐次微分方程，它的右端能被展开成如下的 Fourier级数形式： (3.5) 为正整数或无穷，取值依赖于函数 式中 。 微分方程 (3.3a) 的周期解中不含长期项，这要求 ，解这个方程能够给出 方程(3.5)右端系数 非线性振动的解析逼近频率 与逼近周期 。 方程(3.3a)的通解可写为 (3.6) 式中未知常量 可通过引入方程(3.3b) 中初始条件来 确定，它是 的函数。 以上过程的 可以取任意值，但是， 对于许多非 线性问题只需计算到 即可。 * 0.072547 0.074262 0.074158 100.0 * 0.145065 0.148490 0.148283 50.0 * 0.720731 0.737273 0.736289 10.0 * 1.413827 1.443611 1.441937 5.00 5.291104 3.141593 3.181337 3.179723 2.00 4.859503 4.749642 4.768407 4.768022 1.00 5.771682 5.765846 5.768867 5.768846 0.50 6.259763 6.259755 6.259762 6.259762 0.10 6.282950 6.282950 6.282950 6.282950 0.01 表3.1 各解析逼近周期与“精确”周期的比较 例子：考虑如下Duffing方程 (3.7) 此外，当振幅趋于正无穷时还有以下结果 (3.10) (3.11) 从表3.1及式(3.10,11)知本文提出的迭代法得到的两个 解析逼近周期公式无论对小振幅还是大振幅都具有非 常好的逼近精度。而 Mickens 提出的迭代方法给出的 解析逼近公式仅适用于小振幅的情况。本文新方法给 出的解析逼近周期解的精度同样很高。 2. 非线性恢复力为位移一般函数的情形 初始输入函数可取为 (3.12) 此时式(3.3a)的右端能被展成如下的Fourier级数形式 (3.13) 为正整数或无穷，取值依赖于函数 式中 。 利用与上节相同的迭代方法能构造出这类问题的 解析逼近周期与周期解。 6.935163 7.619846 8.397068 8.336019 -0.7951 0.4750 6.862182 7.277470 7.725633 7.706477 -0.7046 0.4500 6.731984 6.907137 7.129134 7.124567 -0.5725 0.4000 6.621139 6.693309 6.821991 6.820465 -0.4688 0.3500 6.527985 6.552044 6.629916 6.629357 -0.3804 0.3000 6.389680 6.384659 6.411459 6.411392 -0.2315 0.2000 6.309475 6.305863 6.31