代码主要为四阶龙格库塔求解四自由度动力学模型，可无缝衔接时变刚度的导入以及后续振动加速度、位移的提取，可出相图等非线性结果。稍加修改方程即可完成简单的六自由度动力学模型求解。主要适用于刚学习，齿轮动力学同学。

4.4 内啮合齿轮副计算 一对内啮合齿轮副，模数=0.6，小齿轮齿数 25 齿、内齿圈 75 齿，齿宽 b=7，中心距 15。 齿轮的转速 n=3000 r/min，小齿轮输入功率 P=1KW，无冲击 KA=1，每天工作８小时，使用 期限不低于 3 年。 小齿轮选用 20MnCrTi 材料制造并经渗碳淬火磨齿，齿圈采用 40Cr 钢制 造并经高频淬火，A 型齿面齿根均淬火，表面硬度均应在 50HRC 以上，请验证其安全系数 及 佳的变位系数， 大过载载荷。 解： KISSSOFT 内啮合齿轮幅与外啮合部分参数用负号区别，如：中心距以负值表示 a= -15，内齿圈齿数= -75，内齿圈变位系数以负表示齿变薄（齿顶圆增大），X= -0.3 不同与 我国大部分软件。 普通硬齿面齿轮损坏一般是疲劳损坏，因此变位系数一般考虑尺寸的强度与磨损，而 等滑动率刚好可以满足要求，对小齿轮大的变位。 1.打开KISSSFOFT界面，进入圆柱齿轮幅单元。输入指定的参数如图4.8所示。 图 4.8 2.单击变位系数自动调整按钮 ，选择第一个选项 佳优化滑动率，如图4.9所示。

采用齿形啮合式的方法连接破碎机锤头的柄部和头部,其柄部和头部分别采用高锰钢和合金高铬铸铁经消失模铸造而成,该工艺能够取代镶嵌铸造的方法生产中、小型的破碎机锤头。制作工艺实现了锤头的柄部长时间使用,锤头的头部快速更换。节约金属材料,提高生产效率。

求齿轮时变啮合刚度的程序，包含子程序，能跑出图，大家多多下载

该模型考虑时变啮合刚度、传递误差、齿侧间隙、齿轮偏心及齿面摩擦等非线性因素，应用Runge-Kutta算法对系统的微分方程进行求解，结合系统时域图、FFT频谱图、相图、Poincaré截面图和三维频谱图，分析齿侧间隙与偏心量对系统响应的影响。

斜齿轮时变啮合刚度计算

计算齿轮时变啮合刚度的函数，里面包括了斜齿轮的，可以使用

用于计算无故障的斜齿轮对啮合时接触线长度随时间的变化

采用能量法计算时变啮合刚度，通过计算各部分变形及各部分刚度，最终求得综合啮合刚度

比较全的齿轮啮合理论！.....................................