~关于强非线性时变参数动力学系统的求解~
这是我在一篇论文上看到的方程,原文用的是变步长5~6阶r-k法求解。简便起见,用MATLAB里的ode45先算了µ=0时相图,但得到的结果是发散的,而不是原文中的周期解。后用c编的4阶变步长r-k也是得到类似的结果。在此有几个问题向师兄们请教:1.r-k法步长的设置对结果的影响,是否会导致最终结果的发散?
2.初值的设置是否会导致结果的发散?又该怎样选择初值呢?
3.如果算法失败又该怎样选择其他算法呢?(如原文中的变步长5~6阶r-k法) 你仿真的是齿轮非线性啮合力吧?考虑了齿轮间隙和时变啮合刚度?
回复二楼
对,这是考虑齿面摩擦,齿侧间隙和时变啮合刚度,现在就是考虑摩擦这一因素,但是感觉很难出结果,就请各位讨论一下关于齿轮动力学的问题 请问楼主 您这些公式出自哪篇文献?能否给出题目?
页:
[1]