15自由度的齿轮系方程数值解法的问题

我不是小偷

给大家问个好先！
本人目前已经编好了一个行星齿轮系的15自由度的方程，利用状态方程表述后有30个方程组，初值设定都是0，看到很多文献中有使用ode45、ode15s的方法对方程进行求解从而求时间历程图，相图，庞加莱图等等。但是这两种方法在我运行程序后，程序一直处于busy状态，ode45我运行了一晚上还没有出来结果，这里边一定有问题，想同各位老师进行讨教，是我的方程有问题还是这个比较大的的自由度方程组应该换一种解法进行求解？我已经能简化的都简化了。
方程中考虑了时变啮合刚度（KLRP=Kmr+Kar*sin(w0*t)）、和齿侧间隙、啮合阻尼是常数。
方程已经作了无量纲处理，参数也是对照文献没有问题。还请各位大牛不惜赐教。在此谢过了！
自由度如图。