基于解析-有限元的复杂齿基直齿圆柱齿轮啮合刚度计算方法技术

本发明专利技术提出了基于解析‑有限元的复杂齿基直齿圆柱齿轮啮合刚度计算方法，旨在提高含有复杂齿基的直齿圆柱齿轮啮合刚度的计算效率。实现步骤为：基于解析法计算直齿圆柱齿轮轮齿的弯曲、剪切和径向压缩刚度；计算接触刚度；计算单齿啮合刚度；建立含有复杂齿基的直齿圆柱齿轮副的有限元模型；在啮合点处用Mass21单元建立主质量点，并将主质量点与对应的齿基部分进行刚性耦合；将处于啮合状态的轮齿对应的主质量点用Combin14单元连接；计算时啮合刚度。本发明专利技术考虑了含有腹板和减重孔的复杂齿基对直齿圆柱齿轮啮合刚度的影响，采用了解析‑有限元法计算了直齿圆柱齿轮的时变啮合刚度，提高了计算效率，可用于直齿圆柱齿轮的动态性能分析及优化设计。

Meshing stiffness calculation method of complex spur gear based on analytic-finite element method

The method for calculating meshing stiffness of complex spur gear based on analytical and finite element method is proposed in order to improve the efficiency of calculating meshing stiffness of spur gear with complex tooth base. The realization steps are as follows: calculating bending, shearing and radial compression stiffness of spur gear teeth based on analytical method; calculating contact stiffness; calculating meshing stiffness of single tooth; establishing finite element model of spur gear pair with complex tooth base; establishing main mass point with Mass21 element at meshing point, and rigidly coupling the main mass point with corresponding tooth base part; and being in rodent position. The main mass points corresponding to the closed gear teeth are connected by Combin14 element, and meshing stiffness is calculated. The method considers the influence of complex tooth base with web and weight-reducing hole on the meshing stiffness of spur gear. The time-varying meshing stiffness of spur gear is calculated by analytic finite element method, which improves the calculation efficiency and can be used for dynamic performance analysis and optimal design of spur gear.

【技术实现步骤摘要】
基于解析-有限元的复杂齿基直齿圆柱齿轮啮合刚度计算方法
本专利技术属于齿轮精密传动技术和力学分析领域，涉及基于解析-有限元的复杂齿基直齿圆柱齿轮啮合刚度计算方法，可用于直齿圆柱齿轮的动态性能分析及优化设计。
技术介绍
在齿轮传统系统中，主要包括直齿轮传动、斜齿轮传动和蜗轮蜗杆传动等多种传动方式，其中的直齿轮传动又可进一步分为直齿圆柱齿轮传动、直齿圆锥齿轮传动和直齿端面齿轮传动等形式，而直齿圆柱齿轮传动是在传动领域中运用最为广泛的传动形式之一，这种传动结构不仅具有传动平稳，冲击、振动和噪声较小等特点，还具备体积小、重量轻、传递转矩大、起动平稳和传动比分级精细的优点，因此广泛应用于高端数控装备、自动化精密机械和航空航天等领域。近几年来，高端数控装备正不断向着高速、高精度和长寿命的方向发展，并对直齿圆柱齿轮这一关键部件的动力学特性提出了更高的要求。而啮合刚度周期性的变化是传动系统的主要激励形式之一，这种啮合刚度周期性的变化称之为时变啮合刚度，其直接影响着传动过程中的动力学特性。时变啮合刚度主要是由于直齿圆柱齿轮的重合度一般不为整数，即直齿圆柱齿轮同时参与啮合的齿数随时间呈周期性变化造成的。在直齿圆柱齿轮的实际加工中，为了获得良好的机械性能，直齿圆柱齿轮的齿基往往被加工成腹板带孔的结构。因此能否精确求解这种含有复杂齿基的直齿圆柱齿轮时变啮合刚度是研究直齿圆柱齿轮传动系统动力学特性，进行动态性能分析及优化设计的前提，所以研究含有复杂齿基的直齿圆柱齿轮时变啮合刚度的计算方法很有必要的。从目前公开的资料来看，在计算含有复杂齿基的直齿圆柱齿轮的时变啮合刚度方面主要采用有限元法进行模拟计算。例如，2010年Abbes在其发表的名称为“Dynamicbehaviourmodellingofaflexiblegearsystembytheelasticfoundationtheoryinpresenceofdefects”的期刊论文中公开了一种考虑复杂齿基的直齿圆柱齿轮的动力学分析计算方法。该方法主要是采用有限元仿真的方法获得含有复杂齿基的直齿圆柱齿轮的动力学响应。但是这种方法需要建立完整的有限元模型来模拟真实的情况，需要耗费大量的计算资源，因此存在着计算效率低的不足。
技术实现思路
本专利技术的目的在于克服上述现有技术存在的缺陷，提出了基于解析-有限元的复杂齿基直齿圆柱齿轮啮合刚度计算方法，旨在提高含有复杂齿基的直齿圆柱齿轮啮合刚度的计算效率和精度。为了实现上述目的，本专利技术采用的技术方案包括如下步骤：(1)基于解析法，分别计算直齿圆柱齿轮轮齿刚度，包括轮齿弯曲刚度值kb、剪切刚度值ks和径向压缩刚度值ka；(2)基于赫兹接触理论，计算轮齿的接触刚度值kh；(3)基于刚度串并联理论，计算单齿对啮合刚度值kt；(4)根据含复杂齿基的直齿圆柱齿轮几何参数建立含复杂齿基直齿圆柱齿轮的几何模型，并基于有限单元法，采用Solid185单元进行网格映射划分，形成含复杂齿基的直齿圆柱齿轮的有限元模型；(5)根据齿轮啮合理论，建立刚性耦合区域；(5.1)基于直齿圆柱齿轮的啮合理论，分别在主动齿轮和从动齿轮的啮合点位置处建立关键点，并采用Mass21单元划分关键点，形成主质量点；(5.2)基于主质量点与复杂齿基上轮齿的对应关系，采用Mpc184单元将主质量点与对应的齿基部分进行刚性耦合，形成刚性耦合区域；(6)根据齿轮啮合理论，判断主质量点对应的轮齿是否处于啮合状态，若处于啮合状态，则将处于啮合状态的轮齿对应的主质量点用Combin14单元连接，并将Combin14单元的刚度值设置为步骤(3)中得到的kt值；(7)基于静态分析，计算含复杂齿基的直齿圆柱齿轮的啮合刚度值k；(8)判断是否完成一个啮合周期的刚度计算，若没有完成一个啮合周期的刚度计算，则直齿圆柱齿轮旋转一定角度，重复步骤(1)至步骤(7)，若完成了一个啮合周期的刚度计算，则停止计算，得到含复杂齿基直齿圆柱齿轮的时变啮合刚度。作为本专利技术的一种优选，步骤(1)中所述的分别计算直齿圆柱齿轮轮齿刚度，包括轮齿弯曲刚度值kb、剪切刚度值ks和径向压缩刚度值ka，其计算表达式分别为：其中，kb为轮齿的弯曲刚度，ks为轮齿的剪切刚度，ka为轮齿的径向压缩刚度，E为弹性模量，G为剪切模量，为分度圆压力角，rb为基圆半径，τc＝αc-(θb-invαc)，为接触点处的压力角，θb为基齿角的一半，为齿顶高系数，m为模数，γ和τ为积分范围，y1为过渡曲线上任一点处的水平坐标，y2为渐开线上任一点处的水平坐标，Iy1为过渡曲线上任意位置处的截面惯性矩，Iy2为渐开线上任意位置处的截面惯性矩，Ay1为过渡曲线上任意位置处的截面面积，Ay2为渐开线上任意位置处的截面面积。作为本专利技术的一种优选，步骤(2)中所述的计算轮齿的接触刚度值kh，其计算表达式为：其中，kh为轮齿之间的接触刚度，E为弹性模量，B为切片的宽度，v为切片的泊松比。作为本专利技术的一种优选，步骤(3)中所述的计算单齿对啮合刚度kt，其计算表达式为：其中，kt为单齿对的啮合刚度，kb1为主动轮轮齿的弯曲刚度，ks1为主动轮轮齿的剪切刚度，ka1为主动轮轮齿的径向压缩刚度，kb2为从动轮轮齿的弯曲刚度，ks2为从动轮轮齿的剪切刚度，ka2为从动轮轮齿的径向压缩刚度，kh为轮齿之间的接触刚度。作为本专利技术的一种优选，步骤(4)中所述的根据含复杂齿基的直齿圆柱齿轮几何参数建立含复杂齿基直齿圆柱齿轮的几何模型，其具体建立过程为：首先，将直齿圆柱齿轮简化为一个圆柱体，圆柱体的高度等于直齿圆柱齿轮的齿宽，圆柱体的半径等于直齿圆柱齿轮的齿根圆半径；其次，根据直齿圆柱齿轮的参数，在圆柱体上完成腹板部分的建模；最后，根据直齿圆柱齿轮的参数，在腹板上完成减重孔的建模，从而形成含复杂齿基的直齿圆柱齿轮的几何模型。作为本专利技术的一种优选，步骤(7)中所述的计算含复杂齿基的直齿圆柱齿轮的啮合刚度值k，其计算表达式为：其中，T为施加在主动轮中心孔处的扭矩，rb1为主动轮中心孔半径，Δθ为主动轮中心孔的变形。本专利技术与现有技术相比，具有如下优点：本专利技术采用解析法与有限元法相结合的方式来计算含有复杂齿基的直齿圆柱齿轮的时变啮合刚度，并且在主质量点的连接方式上采用了Combin14单元进行连接，相比现有的只通过有限元法来计算含有复杂齿基的直齿圆柱齿轮的时变啮合刚度，本专利技术在计算效率上得到了很大的提高，从原来的10个小时左右缩短至2-3分钟。附图说明图1为本专利技术的实现流程图；图2为本专利技术实施例的直齿圆柱齿轮受力的示意图；图3为本专利技术实施例的直齿圆柱齿实体齿基结构的示意图；图4为本专利技术实施例的直齿圆柱齿只含有腹板的齿基结构的示意图；图5为本专利技术实施例的直齿圆柱齿含有腹板和腹板孔的复杂齿基的示意图；图6为本专利技术实施例的含复杂齿基的直齿圆柱齿轮有限元模型的示意图；图7为本专利技术实施例的齿对活动线的示意图；图8为本专利技术实施例的时变啮合刚度结果对比图。具体实施方式下面结合附图，并以表1所示的直齿圆柱齿轮结构参数为例，对本专利技术作进一步详细说明。表1参考图1，基于解析-有限元的复杂齿基直齿圆柱齿轮啮合刚度计算方法，包括如下步骤：步骤1)基于解析法，分别计算直齿圆柱齿轮轮齿刚度，包括轮齿弯曲刚度值kb、剪本文档来自技高网...

【技术保护点】
1.基于解析‑有限元的复杂齿基直齿圆柱齿轮啮合刚度计算方法，其特征在于，包括如下步骤：(1)基于解析法，分别计算直齿圆柱齿轮轮齿刚度，包括轮齿弯曲刚度值kb、剪切刚度值ks和径向压缩刚度值ka；(2)基于赫兹接触理论，计算轮齿的接触刚度值kh；(3)基于刚度串并联理论，计算单齿对啮合刚度值kt；(4)根据含复杂齿基的直齿圆柱齿轮几何参数建立含复杂齿基直齿圆柱齿轮的几何模型，并基于有限单元法，采用Solid185单元进行网格映射划分，形成含复杂齿基的直齿圆柱齿轮的有限元模型；(5)根据齿轮啮合理论，建立刚性耦合区域：基于直齿圆柱齿轮的啮合理论，分别在主动齿轮和从动齿轮的啮合点位置处建立关键点，并采用Mass21单元划分关键点，形成主质量点；基于主质量点与复杂齿基上轮齿的对应关系，采用Mpc184单元将主质量点与对应的齿基部分进行刚性耦合，形成刚性耦合区域；(6)根据齿轮啮合理论，判断主质量点对应的轮齿是否处于啮合状态，若处于啮合状态，则将处于啮合状态的轮齿对应的主质量点用Combin14单元连接，并将Combin14单元的刚度值设置为步骤(3)中得到的kt值；(7)基于静态分析，计算含复杂齿基的直齿圆柱齿轮的啮合刚度值k；(8)判断是否完成一个啮合周期的刚度计算，若没有完成一个啮合周期的刚度计算，则直齿圆柱齿轮旋转一定角度，重复步骤(1)至步骤(7)，若完成了一个啮合周期的刚度计算，则停止计算，得到含复杂齿基直齿圆柱齿轮的时变啮合刚度。...

【技术特征摘要】
2018.06.12 CN 20181060094401.基于解析-有限元的复杂齿基直齿圆柱齿轮啮合刚度计算方法，其特征在于，包括如下步骤：(1)基于解析法，分别计算直齿圆柱齿轮轮齿刚度，包括轮齿弯曲刚度值kb、剪切刚度值ks和径向压缩刚度值ka；(2)基于赫兹接触理论，计算轮齿的接触刚度值kh；(3)基于刚度串并联理论，计算单齿对啮合刚度值kt；(4)根据含复杂齿基的直齿圆柱齿轮几何参数建立含复杂齿基直齿圆柱齿轮的几何模型，并基于有限单元法，采用Solid185单元进行网格映射划分，形成含复杂齿基的直齿圆柱齿轮的有限元模型；(5)根据齿轮啮合理论，建立刚性耦合区域：基于直齿圆柱齿轮的啮合理论，分别在主动齿轮和从动齿轮的啮合点位置处建立关键点，并采用Mass21单元划分关键点，形成主质量点；基于主质量点与复杂齿基上轮齿的对应关系，采用Mpc184单元将主质量点与对应的齿基部分进行刚性耦合，形成刚性耦合区域；(6)根据齿轮啮合理论，判断主质量点对应的轮齿是否处于啮合状态，若处于啮合状态，则将处于啮合状态的轮齿对应的主质量点用Combin14单元连接，并将Combin14单元的刚度值设置为步骤(3)中得到的kt值；(7)基于静态分析，计算含复杂齿基的直齿圆柱齿轮的啮合刚度值k；(8)判断是否完成一个啮合周期的刚度计算，若没有完成一个啮合周期的刚度计算，则直齿圆柱齿轮旋转一定角度，重复步骤(1)至步骤(7)，若完成了一个啮合周期的刚度计算，则停止计算，得到含复杂齿基直齿圆柱齿轮的时变啮合刚度。2.根据权利要求1所述的基于解析-有限元的复杂齿基直齿圆柱齿轮啮合刚度计算方法，其特征在于：步骤(1)中所述的分别计算直齿圆柱齿轮轮齿刚度，包括轮齿弯曲刚度值kb、剪切刚度值ks和径向压缩刚度值ka，其计算表达式分别为：其中，kb为轮齿的弯曲刚度，ks为轮齿的剪切刚度，ka为轮齿的径向压缩刚度，E为弹性模量，G为剪切模量，为分度圆压力角，...

【专利技术属性】
技术研发人员：付金涛，王奇斌，朱次节，闫海鹏，袁剑，
申请(专利权)人：南京聚能传动设备有限公司，
类型：发明
国别省市：江苏,32