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]
查看完整版本: ~关于强非线性时变参数动力学系统的求解~