[发明专利]一种齿轮传动系统非线性动力学建模方法

说明书

技术领域

本发明属于齿轮传动系统动力学仿真技术领域，具体涉及齿轮传动系统的非线性动力学建模方法。

背景技术

齿轮是机械传动设备中应用最广泛而且是最容易出现故障的零件之一。时变的啮合刚度是齿轮传动系统振动响应的主要动态激励源之一。在通用多体动力学软件中，难以直接将时变啮合刚度加入到齿轮系统动力学模型中。为了体现时变啮合刚度，现有的解决方法中，如把齿轮副以柔性体建模，然后在啮合轮齿之间添加柔性-柔性接触，这种方法可以较准确地体现时变啮合刚度对系统的影响，但缺点是计算量大，效率低。还有以纯刚体建模，通过定义“齿轮运动副”或添加刚体-刚体接触力来体现齿轮副的啮合关系，这种方法虽然计算量小，但是啮合刚度值与理论值相差甚远。

发明内容

鉴于上述背景技术，本发明通过等效转换的方法，巧妙的把时变啮合刚度对系统的影响以啮合力和啮合扭矩的形式体现在纯刚体动力学模型中，使模型达到既精确又高效的目的。

一种齿轮传动系统非线性动力学建模方法，其特征在于系统由以下部件构成：电机、输入轴、输入端轴承、主动齿轮、被动齿轮、输出端轴承、输出轴、负载。

建模方法包括如下步骤：

（1）确定齿轮传动系统的物理参数和运行参数为：主动齿轮和被动齿轮的齿数Z₁和Z₂，质量m₁和m₂，转动惯量I₁和I₂，输入轴扭转刚度和扭转阻尼系数k_p和c_p，输出轴扭转刚度和扭转阻尼系数k_g和c_g，输入端轴承的水平方向支撑刚度和阻尼系数k_x1和c_x1，竖直方向支撑刚度和阻尼系数k_y1和c_y1，输出端轴承的水平方向支撑刚度和阻尼系数k_x2和c_x2，竖直方向支撑刚度和阻尼系数k_y2和c_y2，主动齿轮和被动齿轮的基圆半径R_b1和R_b2；电机输入扭矩M_m，负载扭矩M_b；电机转频f；齿轮副的时变啮合刚度和时变啮合阻尼系数k_t和c_t。

（2）以用啮合力和啮合扭矩等价表示齿轮副非线性啮合关系为核心，构建系统的非线性动力学模型：电机表示为一个集中质量的刚体，在其上施加绕轴的逆时针方向扭矩M_m；输入轴表示为一个单自由度扭转弹簧阻尼器，扭转刚度和扭转阻尼系数分别为k_p和c_p；输入端轴承表示为一对互相垂直的弹簧阻尼器，水平方向的刚度和阻尼系数分别为k_x1和c_x1，竖直方向的刚度和阻尼系数分别为k_y1和c_y1；主动齿轮表示为一个集中质量的刚体，在其质心处施加一个竖直方向的啮合力F₁，一个绕轴的逆时针方向啮合扭矩M₁；被动齿轮表示为一个集中质量的刚体，在其质心处施加一个竖直方向的力F₂，一个绕轴的逆时针方向啮合扭矩M₂；输出轴表示为一个单自由度扭转弹簧阻尼器，扭转刚度和扭转阻尼系数分别为k_g和c_g；输出端轴承表示为一对互相垂直的弹簧阻尼器，水平方向的刚度和阻尼系数分别为k_x2和c_x2，竖直方向的刚度和阻尼系数分别为k_y2和c_y2；负载表示为一个集中质量的刚体，在其上施加绕轴的逆时针方向扭矩M_b。

（3）计算啮合力F₁、F₂，啮合扭矩M₁和M₂，如下：

F₁＝k_t(y₂-y₁-R_b2θ₂-R_b1θ₁)+c_t(y₂-y₁-R_b2θ₂-R_b1θ₁)  (1)