申请/专利权人：南京航空航天大学

申请日：2023-10-18

公开（公告）日：2024-03-01

公开（公告）号：CN117634057A

主分类号：G06F30/17

分类号：G06F30/17;G06F30/23;G06F17/11;G06F17/16

优先权：

专利状态码：在审-实质审查的生效

法律状态：2024.03.19#实质审查的生效;2024.03.01#公开

摘要：本发明公开了一种含剥落故障的弧齿锥齿轮时变啮合刚度计算方法，包括采集弧齿锥齿轮的属性参数和剥落故障参数、轮齿接触分析、承载轮齿接触分析的步骤，并输出时变啮合刚度的计算结果，本发明所提供的方法仅需提供弧齿锥齿轮的属性参数和剥落故障参数即可运行，不受限于操作人员的技术基础，且耗时短，在保证准确性的前提下，大幅提高了计算效率。

主权项：1.一种含剥落故障的弧齿锥齿轮时变啮合刚度计算方法，其特征在于，包括以下步骤：1采集弧齿锥齿轮的属性参数和剥落故障参数；2同时联立求解啮合方程和旋转投影面方程得到弧齿锥齿轮大齿轮齿面上任意一点的坐标以及弧齿锥齿轮小齿轮齿面上任意一点的坐标；3将弧齿锥齿轮小齿轮齿面方程和法向量从轮坯坐标系转化到装配坐标系；将弧齿锥齿轮大齿轮齿面方程和法向量从轮坯坐标系转化到装配坐标系；4根据弧齿锥齿轮大、小齿轮连续相切方程组计算出齿轮副在装配坐标系下第一次啮合时各自所需要的旋转的角度φg0、φp0；其中φg0为弧齿锥齿轮大齿轮第一次啮合时所需要的旋转的角度；φp0为弧齿锥齿轮小齿轮第一次啮合时所需要的旋转的角度；5以所述的φp0为中间值，上下波动一定量φrange假定小齿轮转角的区间[φp0-φrange,φp0+φrange]，在该区间内等间距均分得到不同啮合时刻下弧齿锥齿轮小齿轮的转角φp；通过转角φp解得初始接触点的未知变量值，代入弧齿锥齿轮大、小齿轮齿面方程得到初始接触点的坐标，每个初始接触点代表一个啮合时刻和齿轮啮合位置；6弧齿锥齿轮大、小齿轮齿面上任意一点的坐标为离散点坐标，根据该齿面离散点坐标通过线性插值生成齿轮内部网格；采用三维双线性等参单元建立有限元模型，在齿基底部添加除绕轴旋转外的约束；在不同啮合时刻下大小齿轮接触面齿面点分别施加接触点法向的单位载荷，提取在加载点齿面以下位置的变形，组集成弧齿锥齿轮大齿轮3维整体柔度矩阵[λg]n*i*j及弧齿锥齿轮小齿轮3维整体柔度矩阵[λp]n*i*j和，n表示啮合时刻的个数，i表示齿面施力点的个数，j表示齿面提取位移点的个数；7当受载时，弧齿锥齿轮的啮合状态由接触点变为接触椭圆，将接触椭圆表示为沿接触椭圆长轴的接触曲线；在接触点附近建立极坐标系Lp，θ，将初始接触点法向量沿着切平面移动Lp，分别与大小齿轮相交于两点，两点之间距离即为齿面间隙dp，计算不同θ下的齿面间隙，寻找最小间隙的θmin值，此方向即为接触椭圆的长轴方向，最大间隙的θmax值为接触椭圆的短轴方向；沿着长轴方向，计算不同距离Lp下的齿面间隙ε；8根据[λg]n*i*j、[λp]n*i*j将网格点以样条插值的方法加密，然后用距离潜在接触点最近的插值点或网格点的柔度作为潜在接触点的柔度，得到弧齿锥齿轮小齿轮潜在接触点柔度矩阵λp和弧齿锥齿轮大齿轮潜在接触点柔度矩阵λg，然后整合为潜在接触点柔度矩阵λb；9根据步骤1中采集的剥落顶点，在投影平面内寻找剥落区域穿过的潜在接触点，其为剥落影响的范围点，在剥落影响范围内的潜在接触点间隙增加剥落深度值，得到新的潜在接触点间隙εn×1；10计算时变啮合刚度Kmesh： 式中，εi表示第i个潜在接触点的间隙，Fn_all表示接触力的合力；Ste表示承载传递误差。

全文数据：

权利要求：

百度查询： 南京航空航天大学 一种含剥落故障的弧齿锥齿轮时变啮合刚度计算方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD