一种齿轮传动系统非线性动力学建模方法技术方案

本发明专利技术涉及一种齿轮传动系统非线性动力学建模方法。针对传统的齿轮传动系统动力学模型中的时变啮合刚度和阻尼系数难以在通用多体动力学软件中直接体现的问题，本发明专利技术对动力学模型进行了如下改进：将齿轮副之间的时变啮合刚度和阻尼对系统的激励效果，等效为作用于主被动齿轮质心处的时变啮合力和时变啮合扭矩，并给出了以时变啮合刚度和阻尼系数为自变量的时变啮合力和啮合扭矩的函数表达式。应用该方法建立的齿轮传动系统非线性动力学模型在保证了计算精确度的同时，又有高效的计算速度。

【技术实现步骤摘要】
一种齿轮传动系统非线性动力学建模方法
本专利技术属于齿轮传动系统动力学仿真
，具体涉及齿轮传动系统的非线性动力学建模方法。
技术介绍
齿轮是机械传动设备中应用最广泛而且是最容易出现故障的零件之一。时变的啮合刚度是齿轮传动系统振动响应的主要动态激励源之一。在通用多体动力学软件中，难以直接将时变啮合刚度加入到齿轮系统动力学模型中。为了体现时变啮合刚度，现有的解决方法中，如把齿轮副以柔性体建模，然后在啮合轮齿之间添加柔性-柔性接触，这种方法可以较准确地体现时变啮合刚度对系统的影响，但缺点是计算量大，效率低。还有以纯刚体建模，通过定义“齿轮运动副”或添加刚体-刚体接触力来体现齿轮副的啮合关系，这种方法虽然计算量小，但是啮合刚度值与理论值相差甚远。
技术实现思路
鉴于上述
技术介绍
，本专利技术通过等效转换的方法，巧妙的把时变啮合刚度对系统的影响以啮合力和啮合扭矩的形式体现在纯刚体动力学模型中，使模型达到既精确又高效的目的。一种齿轮传动系统非线性动力学建模方法，其特征在于系统由以下部件构成：电机、输入轴、输入端轴承、主动齿轮、被动齿轮、输出端轴承、输出轴、负载。建模方法包括如下步骤：(1)确定齿轮传动系统的物理参数和运行参数为：主动齿轮和被动齿轮的齿数Z1和Z2，质量m1和m2，转动惯量I1和I2，输入轴扭转刚度和扭转阻尼系数kp和cp，输出轴扭转刚度和扭转阻尼系数kg和cg，输入端轴承的水平方向支撑刚度和阻尼系数kx1和cx1，竖直方向支撑刚度和阻尼系数ky1和cy1，输出端轴承的水平方向支撑刚度和阻尼系数kx2和cx2，竖直方向支撑刚度和阻尼系数ky2和cy2，主动齿轮和被动齿轮的基圆半径Rb1和Rb2；电机输入扭矩Mm，负载扭矩Mb；电机转频f；齿轮副的时变啮合刚度和时变啮合阻尼系数kt和ct。(2)以用啮合力和啮合扭矩等价表示齿轮副非线性啮合关系为核心，构建系统的非线性动力学模型：电机表示为一个集中质量的刚体，在其上施加绕轴的逆时针方向扭矩Mm；输入轴表示为一个单自由度扭转弹簧阻尼器，扭转刚度和扭转阻尼系数分别为kp和cp；输入端轴承表示为一对互相垂直的弹簧阻尼器，水平方向的刚度和阻尼系数分别为kx1和cx1，竖直方向的刚度和阻尼系数分别为ky1和cy1；主动齿轮表示为一个集中质量的刚体，在其质心处施加一个竖直方向的啮合力F1，一个绕轴的逆时针方向啮合扭矩M1；被动齿轮表示为一个集中质量的刚体，在其质心处施加一个竖直方向的力F2，一个绕轴的逆时针方向啮合扭矩M2；输出轴表示为一个单自由度扭转弹簧阻尼器，扭转刚度和扭转阻尼系数分别为kg和cg；输出端轴承表示为一对互相垂直的弹簧阻尼器，水平方向的刚度和阻尼系数分别为kx2和cx2，竖直方向的刚度和阻尼系数分别为ky2和cy2；负载表示为一个集中质量的刚体，在其上施加绕轴的逆时针方向扭矩Mb。(3)计算啮合力F1、F2，啮合扭矩M1和M2，如下：F1＝kt(y2-y1-Rb2θ2-Rb1θ1)+ct(y2-y1-Rb2θ2-Rb1θ1)(1)F2＝F1(2)M1＝F1.Rb1(3)M2＝F2.Rb2(4)以上各式中，kt和ct分别为齿轮副的时变啮合刚度和时变啮合阻尼系数；y1和y2分别为主动齿轮和被动齿轮的竖直方向位移，是方程组的因变量；θ1、θ2分别为主动齿轮和被动齿轮的角位移。(4)将由步骤(1)和步骤(3)获得的各参数数值代入到步骤(2)中的系统动力学模型中，采用牛顿第二定律推导出系统的非线性动力学微分方程，即可对系统进行非线性动力学求解。本专利技术在传统的齿轮系统非线性动力学模型的基础上，提出以时变啮合力和时变啮合扭矩来等效表示齿轮副之间的非线性啮合关系，克服了在多体动力学求解器软件平台中难以在齿轮副之间添加非线性时变弹簧阻尼器的问题。附图说明图1.本专利技术的工作流程图；图2.传统的齿轮传动系统动力学模型；图3.本专利技术所建立的齿轮传动系统动力学模型；图4.用本专利技术方法仿真的齿轮传动系统振动加速度信号时域和频域图。具体实施方式如图1所示，是本专利技术的一种齿轮传动系统非线性动力学建模方法的工作流程图。具体实施过程如下：(1)确定齿轮传动系统的物理参数和运行参数为：主动齿轮和被动齿轮的齿数Z1和Z2，质量m1和m2，转动惯量I1和I2，输入轴扭转刚度和扭转阻尼系数kp和cp，输出轴扭转刚度和扭转阻尼系数kg和cg，输入端轴承的水平方向支撑刚度和阻尼系数kx1和cx1，竖直方向支撑刚度和阻尼系数ky1和cy1，输出端轴承的水平方向支撑刚度和阻尼系数kx2和cx2，竖直方向支撑刚度和阻尼系数ky2和cy2，主动齿轮和被动齿轮的基圆半径Rb1和Rb2；电机输入扭矩Mm，负载扭矩Mb；电机转频f；齿轮副的时变啮合刚度和时变啮合阻尼系数kt和ct。(2)以用啮合力和啮合扭矩等效表示齿轮副非线性啮合关系为核心，构建系统的非线性动力学模型：电机表示为一个集中质量的刚体，在其上施加绕轴的逆时针方向扭矩Mm；输入轴表示为一个单自由度扭转弹簧阻尼器，扭转刚度和扭转阻尼系数分别为kp和cp；输入端轴承表示为一对互相垂直的弹簧阻尼器，水平方向的刚度和阻尼系数分别为kx1和cx1，竖直方向的刚度和阻尼系数分别为ky1和cy1；主动齿轮表示为一个集中质量的刚体，在其质心处施加一个竖直方向的啮合力F1，一个绕轴的逆时针方向啮合扭矩M1；被动齿轮表示为一个集中质量的刚体，在其质心处施加一个竖直方向的力F2，一个绕轴的逆时针方向啮合扭矩M2；输出轴表示为一个单自由度扭转弹簧阻尼器，扭转刚度和扭转阻尼系数分别为kg和cg；输出端轴承表示为一对互相垂直的弹簧阻尼器，水平方向的刚度和阻尼系数分别为kx2和cx2，竖直方向的刚度和阻尼系数分别为ky2和cy2；负载表示为一个集中质量的刚体，在其上施加绕轴的逆时针方向扭矩Mb。(3)计算啮合力F1、F2，啮合扭矩M1和M2，如下：F1＝kt(y2-y1-Rb2θ2-Rb1θ1)+ct(y2-y1-Rb2θ2-Rb1θ1)(1)F2＝F1(2)M1＝F1.Rb1(3)M2＝F2.Rb2(4)以上各式中，y1和y2分别为主动齿轮和被动齿轮的竖直方向位移，是方程组的因变量；θ1、θ2分别为主动齿轮和被动齿轮的角位移。(4)将由步骤(1)和步骤(3)获得的各参数数值代入到系统动力学模型中，采用牛顿第二定律推导出系统的非线性动力学微分方程，即可对系统进行非线性动力学求解。本实施例中的系统动力学方程为：...

【技术保护点】
一种齿轮传动系统非线性动力学建模方法，其特征在于系统由以下部件构成：电机、输入轴、输入端轴承、主动齿轮、被动齿轮、输出端轴承、输出轴、负载；建模方法包括如下步骤：（1）确定齿轮传动系统的物理参数和运行参数为：主动齿轮和被动齿轮的齿数Z1和Z2，质量m1和m2，转动惯量I1和I2，输入轴扭转刚度和扭转阻尼系数kp和cp，输出轴扭转刚度和扭转阻尼系数kg和cg，输入端轴承的水平方向支撑刚度和阻尼系数kx1和cx1，竖直方向支撑刚度和阻尼系数ky1和cy1，输出端轴承的水平方向支撑刚度和阻尼系数kx2和cx2，竖直方向支撑刚度和阻尼系数ky2和cy2，主动齿轮和被动齿轮的基圆半径Rb1和Rb2；电机输入扭矩Mm，负载扭矩Mb；电机转频f；齿轮副的时变啮合刚度和时变啮合阻尼系数kt和ct；（2）以用啮合力和啮合扭矩等效表示齿轮副非线性啮合关系为核心，构建系统的非线性动力学模型：电机表示为一个集中质量的刚体，在其上施加绕轴的逆时针方向扭矩Mm；输入轴表示为一个单自由度扭转弹簧阻尼器，扭转刚度和扭转阻尼系数分别为kp和cp；输入端轴承表示为一对互相垂直的弹簧阻尼器，水平方向的刚度和阻尼系数分别为kx1和cx1，竖直方向的刚度和阻尼系数分别为ky1和cy1；主动齿轮表示为一个集中质量的刚体，在其质心处施加一个竖直方向的啮合力F1，一个绕轴的逆时针方向啮合扭矩M1；被动齿轮表示为一个集中质量的刚体，在其质心处施加一个竖直方向的力F2，一个绕轴的逆时针方向啮合扭矩M2；输出轴表示为一个单自由度扭转弹簧阻 尼器，扭转刚度和扭转阻尼系数分别为kg和cg；输出端轴承表示为一对互相垂直的弹簧阻尼器，水平方向的刚度和阻尼系数分别为kx2和cx2，竖直方向的刚度和阻尼系数分别为ky2和cy2；负载表示为一个集中质量的刚体，在其上施加绕轴的逆时针方向扭矩Mb；（3）计算啮合力F1、F2，啮合扭矩M1和M2，如下：F1＝kt(y2?y1?Rb2θ2?Rb1θ1)+ct(y2?y1?Rb2θ2?Rb1θ1)??(1)F2＝?F1??(2)M1＝F1·Rb1??(3)M2＝?F2·Rb2??(4)以上各式中，kt和ct分别为齿轮副的时变啮合刚度和时变啮合阻尼系数；y1和y2分别为主动齿轮和被动齿轮的竖直方向位移，是方程组的因变量；θ1、θ2分别为主动齿轮和被动齿轮的角位移；（4）将由步骤（1）和步骤（3）获得的各参数数值代入到步骤（2）中的系统动力学模型中，采用牛顿第二定律推导出系统的非线性动力学微分方程，即可对系统进行非线性动力学求解。...

【技术特征摘要】
1.一种齿轮传动系统非线性动力学建模方法，其特征在于系统由以下部件构成：电机、输入轴、输入端轴承、主动齿轮、被动齿轮、输出端轴承、输出轴、负载；建模方法包括如下步骤：(1)确定齿轮传动系统的物理参数和运行参数为：主动齿轮和被动齿轮的齿数Z1和Z2，质量m1和m2，转动惯量I1和I2，输入轴扭转刚度和扭转阻尼系数kp和cp，输出轴扭转刚度和扭转阻尼系数kg和cg，输入端轴承的水平方向支撑刚度和阻尼系数kx1和cx1，竖直方向支撑刚度和阻尼系数ky1和cy1，输出端轴承的水平方向支撑刚度和阻尼系数kx2和cx2，竖直方向支撑刚度和阻尼系数ky2和cy2，主动齿轮和被动齿轮的基圆半径Rb1和Rb2；电机输入扭矩Mm，负载扭矩Mb；电机转频f；齿轮副的时变啮合刚度和时变啮合阻尼系数kt和ct；(2)以用啮合力和啮合扭矩等效表示齿轮副非线性啮合关系为核心，构建系统的非线性动力学模型：电机表示为一个集中质量的刚体，在其上施加绕轴的逆时针方向扭矩Mm；输入轴表示为一个单自由度扭转弹簧阻尼器，扭转刚度和扭转阻尼系数分别为kp和cp；输入端轴承表示为一对互相垂直的弹簧阻尼器，水平方向的刚度和阻尼系数分别为kx1和cx1，竖直方向的刚度和阻尼系数分别为ky1和cy1；主动齿轮表示为一个集中质量的刚体，在其质心处施加一个竖直方向的啮合力F1，一个绕轴的逆时针方向啮合扭矩M1；被动齿轮表示为一个集中质量的刚体，在其质心处施加一个竖直方向的力F2，一个绕轴的逆时针方向啮合扭矩M2；输出轴表示为一个单自由度扭转弹簧阻尼器，扭转刚度和扭转阻尼系数分别为kg和cg；输出端轴承表示为一对互相垂直的弹簧阻尼器，水平方向的刚度和阻尼系数分别为kx2和cx2，竖直方向的刚度和阻尼系数分别为ky2和cy2；负载表示为一个集中质量的刚体，在其上施加绕轴的逆时针方向扭矩Mb；(3)计算啮合力F1、F2，啮合扭矩M1和M2，如下：F1＝kt(y2-y1-Rb2θ2-Rb1θ1)+ct(y2-y1-Rb2θ2-Rb1θ1)(1)F2＝F1(2)M1＝F1.Rb1(3)M2＝F2.Rb2(4)以上各式中，kt和ct分别为齿轮副的时变啮合刚度和时变啮合阻尼系数；y1和y2分别为主...

【专利技术属性】
技术研发人员：高立新，张飞斌，蔡力钢，王锋，刘傲，徐海，赵鑫，雷涛涛，
申请(专利权)人：北京工业大学，
类型：发明
国别省市：