ypp207 发表于 2008-11-21 12:33

~关于强非线性时变参数动力学系统的求解~

这是我在一篇论文上看到的方程，原文用的是变步长5~6阶r-k法求解。简便起见，用MATLAB里的ode45先算了µ=0时相图，但得到的结果是发散的，而不是原文中的周期解。后用c编的4阶变步长r-k也是得到类似的结果。在此有几个问题向师兄们请教：
1.r-k法步长的设置对结果的影响，是否会导致最终结果的发散？
2.初值的设置是否会导致结果的发散？又该怎样选择初值呢？
3.如果算法失败又该怎样选择其他算法呢？（如原文中的变步长5~6阶r-k法）

health 发表于 2008-11-21 18:03

你仿真的是齿轮非线性啮合力吧？考虑了齿轮间隙和时变啮合刚度？

ypp207 发表于 2008-11-24 11:31

回复二楼

对，这是考虑齿面摩擦，齿侧间隙和时变啮合刚度，现在就是考虑摩擦这一因素，但是感觉很难出结果，就请各位讨论一下关于齿轮动力学的问题

huaijuliu 发表于 2011-12-10 03:05

请问楼主 您这些公式出自哪篇文献？能否给出题目？
页: [1]
查看完整版本: ~关于强非线性时变参数动力学系统的求解~