CN113487180B - 一种基于云平台的齿轮齿面评价方法 - Google Patents

一种基于云平台的齿轮齿面评价方法 Download PDF

Info

Publication number
CN113487180B
CN113487180B CN202110758977.XA CN202110758977A CN113487180B CN 113487180 B CN113487180 B CN 113487180B CN 202110758977 A CN202110758977 A CN 202110758977A CN 113487180 B CN113487180 B CN 113487180B
Authority
CN
China
Prior art keywords
gear
tooth surface
tooth
deviation
point
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202110758977.XA
Other languages
English (en)
Other versions
CN113487180A (zh
Inventor
张登攀
李键
岳佳佳
闫勇刚
郑义博
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Henan University of Technology
Original Assignee
Henan University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Henan University of Technology filed Critical Henan University of Technology
Priority to CN202110758977.XA priority Critical patent/CN113487180B/zh
Publication of CN113487180A publication Critical patent/CN113487180A/zh
Application granted granted Critical
Publication of CN113487180B publication Critical patent/CN113487180B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393Score-carding, benchmarking or key performance indicator [KPI] analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Data Mining & Analysis (AREA)
  • Operations Research (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Educational Administration (AREA)
  • Mathematical Optimization (AREA)
  • Development Economics (AREA)
  • Computational Mathematics (AREA)
  • Game Theory and Decision Science (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Quality & Reliability (AREA)
  • Algebra (AREA)
  • Marketing (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明公开了一种基于云平台的齿轮齿面评价方法,包括以下具体步骤:A:用户通过齿轮测量云终端,向齿轮测量云平台提供待检测齿轮的齿面坐标点云数据;B:齿轮测量云平台计算各个齿面点对应的齿面偏差值得到齿面偏差点云数据;C:齿轮测量云平台将齿面点的坐标值数据和齿面偏差值进行处理,得到归一化处理后的齿面坐标系下的齿面偏差点云数据;D:齿轮测量云平台使用二维勒让德多项式和/或二维切比雪夫多项式的前6项作为基函数拟合得出的待检测齿轮的偏差曲面方程;E:齿轮测量云平台对待检测齿轮进行检测和评价;F:齿轮测量云平台通过齿轮测量云终端反馈结果。本发明能够对齿轮齿面实现更全面更快捷的测量及评价。

Description

一种基于云平台的齿轮齿面评价方法
技术领域
本发明涉及齿轮齿面测量领域,尤其涉及一种基于云平台的齿轮齿面评价方法。
背景技术
齿轮是标准化的典型传动零件,广泛应用于机器设备、仪器仪表的运动传递、负载传递和精密分度。齿轮的精度和动态特性直接决定了机器设备和仪器仪表的工作性能,齿轮的质量直接决定了装备的运行性能、服役寿命、安全性和可靠性。
目前,对一个齿轮进行测量需要使用到多台仪器,存在测量效率低、测量精度不一致等诸多问题。例如,现有的齿轮齿廓偏差测量是对整个齿面中的一条渐开线齿廓迹线进行测量与评价,作为整个齿面的齿廓偏差;齿轮齿距偏差测量是选取齿轮每个齿面上一点替代整个齿面,评价轮齿分布均匀状况。齿轮精度的评定皆基于轮齿表面上局部点的几何误差,由于其反映的是齿面复杂型面的部分信息,因此用其控制齿轮质量是不严格的,用来认识齿轮的误差特性是不完整的、用来分析误差来源是不准确的。为了全面反映齿轮加工质量,需要新的方法对齿轮齿面精度进行测量,进而为齿轮精度评价、齿轮加工工艺参数提供可靠数据。
随着我国互联网技术的发展,基于硬件资源和软件资源服务且提供计算、网络和存储能力的云平台得到了越来越广泛的应用。基于云平台的齿轮齿面评价作为新一代信息技术与齿轮测量技术相结合的产物,对实现数字化、信息化和网络大数据的齿轮测量水平的提高具有重要意义。若将基于云平台的数据互连统一技术引入测量评价***,在云端对齿轮测量评价各个环节中的数据进行采集、处理、分析以及可视化展示,将会使齿轮测量评价技术手段得到了质的飞越。
发明内容
本发明的目的是提供一种基于云平台的齿轮齿面评价方法,能够对齿轮齿面实现更全面更快捷的测量及评价,解决现有齿轮测量效率低、测量精度不一致等缺陷。
本发明采用下述技术方案:
一种基于云平台的齿轮齿面评价方法,依次包括以下具体步骤:
A:用户通过齿轮测量云终端,向齿轮测量云平台提供待检测齿轮的齿面坐标点云数据,齿面坐标点云数据包括待检测齿轮上若干个齿面点的坐标数据,每个齿面点的坐标数据中均包含该齿面点的三维坐标值;
其中,齿轮测量云终端即为用户所使用的本地客户端,用于与齿轮测量云平台实现数据交互,为用户提供齿轮齿面评价服务和齿轮齿面评价结果展示服务;齿轮测量云平台即为云上服务器端,用于对用户通过齿轮测量云终端提供的待检测齿轮进行计算及评价,并将评价结果反馈至齿轮测量云终端;
B:对于齿面坐标点云数据中的每一个齿面点的坐标值数据,齿轮测量云平台分别通过下述公式计算各个齿面点对应的齿面偏差值,然后将齿面坐标点云数据中所有的齿面点所对应的齿面偏差值组合构成齿面偏差点云数据;
Figure BDA0003148498570000021
其中,dlot表示齿面点对应的齿面偏差值,下角标lot表示根据铅垂线方向,K表示常数系数,ρ表示极径,ηb表示基本空间宽度半角,β表示螺旋角,rb表示基圆半径,X表示变位系数,αn表示齿轮法向压力角,Z表示齿数,齿数中内齿轮为负值,外齿轮为正值,inv表示渐开线函数,αt表示齿轮端面压力角,x、y和z分别表示齿面点在三维笛卡尔坐标系下在x轴、y轴和z轴上的坐标;
C:齿轮测量云平台对齿面坐标点云数据中的每一个齿面点的坐标值数据和齿面偏差值进行坐标转化处理,在从三维笛卡尔坐标系下转化到u-v-d齿面坐标系下后进行归一化处理,最终得到归一化处理后的u′-v′-d齿面坐标系下的齿面偏差点云数据;
D:齿轮测量云平台使用二维勒让德多项式和/或二维切比雪夫多项式的前6项作为基函数,利用归一化处理后的u′-v′-d齿面坐标系下的齿面偏差点云数据,通过最小二乘法拟合得出的待检测齿轮的偏差曲面方程;
其中,利用二维切比雪夫多项式的前6项作为基函数得到的待检测齿轮的偏差曲面方程为:
dC=AC0*1+AC1*u'+Ac2*v'+AC3*(2u'2-1)+AC4*u'v'+AC5*(2v'2-1);
方程系数AC0、AC1、AC2、AC3、AC4和AC5分别用于评价齿距偏差、齿廓斜率、螺旋线斜率、齿廓凸度、轮齿鼓度和齿面扭曲,计算得到的dc即为待检测齿轮的齿面上对应的齿面点或齿面线处的偏差值d;
利用二维勒让德多项式的前6项作为基函数得到的待检测齿轮的偏差曲面方程为:
Figure BDA0003148498570000031
系数AL0、AL1、AL2、AL3、AL4和AL5分别用于评价齿距偏差、齿廓斜率、螺旋线斜率、齿廓凸度、轮齿鼓度和齿面扭曲;计算得到的dL即为待检测齿轮的齿面上对应的齿面点或齿面线处的偏差值d;
E:齿轮测量云终端将用户输入的待评测齿轮上任意一个齿面点(u,v)的坐标或任意一条齿面线u=a的坐标传输至齿轮测量云平台,齿轮测量云平台首先按照步骤C中的方法,对用户输入的齿面点或齿面线的坐标进行坐标转化处理及归一化处理为对应(u′,v′)或u′;然后利用坐标转化处理后的坐标数据和步骤D中得到的待检测齿轮的偏差曲面方程,求出该齿面点或齿面线所对应位置处的偏差值d;齿轮测量云平台还依据待检测齿轮的偏差曲面方程中的系数求得待评测齿轮的齿距偏差、齿廓斜率、螺旋线斜率、齿廓凸度、轮齿鼓度和齿面扭曲,然后齿轮测量云平台结合用户选定的齿轮精度参数和齿轮测量评价标准,对待检测齿轮进行合格或不合格的评价;
F:齿轮测量云平台将计算得到的待评测齿轮的齿距偏差、齿廓斜率、螺旋线斜率、齿廓凸度、轮齿鼓度和齿面扭曲,以及用户输入的齿面点或齿面线所对应位置处的偏差值d和评价结果,传输至齿轮测量云终端并反馈给用户。
所述的步骤C中:
在进行坐标转化处理时,齿轮测量云平台首先将齿面坐标点云数据中的每一个齿面点的坐标值数据和齿面偏差值,从三维笛卡尔坐标系下转化到u-v-d齿面坐标系下,然后再对u-v-d齿面坐标系下的齿面点的坐标值(u,v,d)进行归一化处理得到坐标值(u′,v′,d),其中平面u′v′的取值范围在[-1,1]*[-1,1]内;
坐标转化公式如下:
Figure BDA0003148498570000041
Figure BDA0003148498570000042
d=dlot
其中,u表示沿齿廓方向的坐标值,φm表示极角,Λnom表示渐开线起始角,rb表示基圆半径,ξnom表示齿廓滚动角,θk表示展角,v表示沿齿宽方向的坐标值,Z表示齿数,β表示螺旋角,d表示垂直于uv平面方向的坐标值,dlot表示齿面点对应的齿面偏差值;渐开线起始角Λnom和齿廓滚动角ξnom均为待检测齿轮的基本参数。
所述的步骤C中:
在进行归一化处理时,对于齿面坐标系中任意一点(u,v,d),按照如下方式转化为(u′,v′,d),最终转化为在[-1,1]*[-1,1]矩形区域上的齿面点的坐标值(u′,v′,d);
Figure BDA0003148498570000043
Figure BDA0003148498570000044
d=dlot
其中,Lα表示u的取值范围,Lβ表示v的取值范围。
所述的步骤D中:
单独采用二维切比雪夫多项式作为基函数计算待检测齿轮的偏差曲面方程时,设齿面偏差点云数据中共有s个齿面偏差点数据;
用最小二乘法建立方程组,由该方程组求解出偏差曲面方程dC
Figure BDA0003148498570000051
通过二维切比雪夫多项式求出定义在u′-v′-d齿面坐标系下的偏差曲面方程为:
dC=AC0*1+AC1*u'+Ac2*v'+AC3*(2u'2-1)+AC4*u'v'+AC5*(2v'2-1);
齿距偏差等于AC0,齿廓斜率等于2AC1,螺旋线斜率等于-2AC2,齿廓凸度等于2AC3,轮齿鼓度等于-4AC4,齿面扭曲等于-2AC5,计算得到的偏差值dc即为待检测齿轮的齿面上点(u,v)归一化为(u′,v′)所对应的偏差值d。
所述的步骤D中:
单独采用二维勒让德多项式作为基函数计算待检测齿轮的偏差曲面方程时,设齿面偏差点云数据中共有s个齿面偏差点数据;
用最小二乘法建立方程组,由该方程组求解出偏差曲面方程dL
Figure BDA0003148498570000052
通过二维勒让德多项式所求dL为定义在u′-v′-d齿面坐标系下的偏差曲面方程为:
Figure BDA0003148498570000053
齿距偏差等于AL0,齿廓斜率等于AL1,螺旋线斜率等于AL2,齿廓凸度等于1.5AL3,轮齿鼓度等于2AL4,齿面扭曲等于1.5AL5;计算得到的偏差值dL即为待检测齿轮的齿面上点(u,v)归一化为(u′,v′)所对应的偏差值d。
所述的步骤D中:
分别采用二维勒让德多项式和二维切比雪夫多项式作为基函数得到对应的偏差曲面方程,然后根据对于不同种类的齿面偏差的评价需求,选择对应的偏差曲面方程的系数得到更为准确的不同种类的齿面偏差评价结果;
设齿面偏差点云数据中共有s个齿面偏差点数据;
采用二维切比雪夫多项式作为基函数得到的偏差曲面方程分别为:
dC=AC0*1+AC1*u'+AC2*v'+AC3*(2u'2-1)+AC4*u'v'+AC5*(2v'2-1);
采用二维切勒让德多项式作为基函数得到的偏差曲面方程为:
Figure BDA0003148498570000061
在需要准确评价齿廓凸度和齿面扭曲时,利用二维切比雪夫多项式拟合出的偏差曲面方程dc中第四项和第六项的系数进行求解,齿廓凸度等于2AC3,齿面扭曲等于-2AC5;在需要准确评价齿距偏差、齿廓斜率、螺旋线斜率和轮齿鼓度时,利用二维勒让德多项式拟合出的偏差曲面方程dL中第一项、第二项、第三项和第五项的系数进行求解,齿距偏差等于AL0,齿廓斜率等于AL1,螺旋线斜率等于AL2,轮齿鼓度等于2AL4
在需要更准确的拟合曲面峰值及细节时选用通过二维切比雪夫多项式拟合的偏差曲面方程dc,作为待检测齿轮的最终偏差曲面方程;在需要更准确的拟合曲面整体趋势及轮廓时选用通过二维勒让德多项式拟合的偏差曲面方程dL,作为待检测齿轮的最终偏差曲面方程。
所述的齿轮精度参数采用6级精度齿轮参数。
本发明采用下述技术方案:首先通过齿轮测量云终端向齿轮测量云平台提供待检测齿轮的齿面坐标点云数据,将通过齿面坐标点云数据计算获得的齿面偏差点云利用二维勒让德多项式和二维切比雪夫多项式分别拟合出偏差曲面方程。然后其拟合的偏差曲面方程的各项系数可以对齿距、齿廓斜率、螺旋线斜率、齿廓凸度、轮齿鼓度、齿面扭曲等进行评价。在u′-v′-d齿面坐标系下结合拟合的偏差曲面方程可以得出该齿面任意点线的偏差。最终将前述算法生成评价指标及数据在云平台通过前端界面可视化表达。
本发明使用包含众多点的齿面坐标点云数据计算相比于传统的以少数几个点计算的二维检测能更全面、更准确的评价整个齿面偏差。另外发明中使用两种拟合方法互补实现齿面偏差评价相比于单种拟合方法也更准确客观。而且本发明方法基于云平台可以更快捷的实现测量评价结果资源共享,进而为后续指导生产提供数据支撑。
附图说明
图1为本发明的流程示意图。
具体实施方式
以下结合附图和实施例对本发明作以详细的描述:
如图1所示,本发明所述的基于云平台的齿轮齿面评价方法,依次包括以下具体步骤:
A:用户通过齿轮测量云终端,向齿轮测量云平台提供待检测齿轮的齿面坐标点云数据,齿面坐标点云数据包括待检测齿轮上若干个齿面点的坐标数据,每个齿面点的坐标数据中均包含该齿面点的三维坐标值;
齿轮测量云终端即为用户所使用的本地客户端,用于与齿轮测量云平台实现数据交互,为用户提供齿轮齿面评价服务和齿轮齿面评价结果展示服务;
齿轮测量云平台即为云上服务器端,用于对用户通过齿轮测量云终端提供的待检测齿轮进行计算及评价,并将评价结果反馈至齿轮测量云终端。
本发明中,齿轮测量云终端和齿轮测量云平台,与现有众多软件程序中的本地客户端和云上服务器构建原理及信息交互原理类似,仅实现功能不同,在此不再赘述。
本发明中,待检测齿轮的齿面点为待检测齿轮的齿面上的一个点,设该齿面点在三维笛卡尔坐标系下的坐标为(x,y,z),齿面坐标点云数据为由待检测齿轮的齿面上的若干个齿面点的坐标值数据所组成的数据组;
B:对于齿面坐标点云数据中的每一个齿面点的坐标值数据,齿轮测量云平台分别通过下述公式计算各个齿面点对应的齿面偏差值,然后将齿面坐标点云数据中所有的齿面点所对应的齿面偏差值组合构成齿面偏差点云数据;
Figure BDA0003148498570000081
/>
其中,dlot表示齿面点对应的齿面偏差值,下角标lot表示根据铅垂线方向,K表示常数系数,ρ表示极径,ηb表示基本空间宽度半角,β表示螺旋角,rb表示基圆半径,X表示变位系数,αn表示齿轮法向压力角,Z表示齿数,齿数中内齿轮为负值,外齿轮为正值,inv表示渐开线函数,αt表示齿轮端面压力角,x、y和z分别表示齿面点在三维笛卡尔坐标系下在x轴、y轴和z轴上的坐标;
本发明中,齿数Z、变位系数X、螺旋角β、法向压力角αn、端面压力角αt和基圆半径rb均为待检测齿轮的基本参数,可通过查询待检测齿轮出厂规格得到。
C:齿轮测量云平台对齿面坐标点云数据中的每一个齿面点的坐标值数据和齿面偏差值进行坐标转化处理,在从三维笛卡尔坐标系下转化到u-v-d齿面坐标系下后进行归一化处理,最终得到归一化处理后的u′-v′-d齿面坐标系下的齿面偏差点云数据;
在进行坐标转化处理时,齿轮测量云平台首先将齿面坐标点云数据中的每一个齿面点的坐标值数据和齿面偏差值,从三维笛卡尔坐标系下转化到u-v-d齿面坐标系下,然后再对u-v-d齿面坐标系下的齿面点的坐标值(u,v,d)进行归一化处理得到坐标值(u′,v′,d),其中平面u′v′的取值范围在[-1,1]*[-1,1]内;
坐标转化公式如下:
Figure BDA0003148498570000082
Figure BDA0003148498570000083
d=dlot
其中,u表示沿齿廓方向的坐标值,φm表示极角,Λnom表示渐开线起始角,rb表示基圆半径,ξnom表示齿廓滚动角,θk表示展角,v表示沿齿宽方向的坐标值,Z表示齿数,β表示螺旋角,d表示垂直于uv平面方向的坐标值,dlot表示齿面点对应的齿面偏差值;渐开线起始角Λnom和齿廓滚动角ξnom均为待检测齿轮的基本参数,可通过查询待检测齿轮出厂规格得到。
由于三维笛卡尔坐标系下无法同时表示出齿面点的三维坐标(x,y,z)及齿面偏差值dlot的对应关系,因此需要将四个参数需要转化到u-v-d齿面坐标系下,以达到多项式拟合的使用条件,在后续用于将二维勒让德多项式和/或二维切比雪夫多项式的前6项作为基函数对齿面偏差点云数据进行最小二乘曲面拟合。
在完成齿面点的坐标值数据和齿面偏差值从三维笛卡尔坐标系到u-v-d齿面坐标系的转化后,此时坐标值(u,v)在[0,Lα]*[0,Lβ]的矩形区域,而二维勒让德多项式和/或二维切比雪夫多项式为定义在[-1,1]*[-1,1]矩形区域内的多项式,因此需要将齿面点在齿面坐标系下的坐标值(u,v)进行归一化处理得到坐标值(u′,v′),(u′,v′)的取值范围在[-1,1]*[-1,1]内;
在进行归一化处理时,对于齿面坐标系中任意一点(u,v,d),按照如下方式转化为(u′,v′,d),最终转化为在[-1,1]*[-1,1]矩形区域上的齿面点的坐标值(u′,v′,d)。
Figure BDA0003148498570000091
Figure BDA0003148498570000092
/>
d=dlot
其中,Lα表示u的取值范围,Lβ表示v的取值范围;
D:齿轮测量云平台使用二维勒让德多项式和/或二维切比雪夫多项式的前6项作为基函数,利用归一化处理后的u′-v′-d齿面坐标系下的齿面偏差点云数据,通过最小二乘法拟合得出的待检测齿轮的偏差曲面方程;
其中,利用二维切比雪夫多项式的前6项作为基函数得到的待检测齿轮的偏差曲面方程为:
dC=AC0*1+AC1*u'+Ac2*v'+AC3*(2u'2-1)+AC4*u'v'+AC5*(2v'2-1);
二维切比雪夫多项式的第一项为C0(u′,v′)=1,其方程系数AC0能够表示齿距偏差,第二项为C1(u′,v′)=u′,其方程系数AC1能够表示齿廓斜率,第三项为C2(u′,v′)=v′,其方程系数AC2能够于表示螺旋线斜率,第四项为C3(u′,v′)=2u′2-1,其方程系数AC3能够表示齿廓凸度,第五项为C4(u′,v′)=u′v′,其方程系数AC4能够表示轮齿鼓度,第六项为C5(u′,v′)=2v′2-1,其方程系数AC5能够齿面扭曲;即系数AC0、AC1、AC2、AC3、AC4和AC5可分别用于评价齿距偏差、齿廓斜率、螺旋线斜率、齿廓凸度、轮齿鼓度和齿面扭曲,计算得到的dc即为待检测齿轮的齿面上对应的齿面点或齿面线处的偏差值d。
上述C0(u′,v′)=1为数学领域中多项式的标准表达方式,含义为多项式第一项的值为1,后续表达式均与此表达含义相同,在此不再赘述。
利用二维勒让德多项式的前6项作为基函数得到的待检测齿轮的偏差曲面方程为:
Figure BDA0003148498570000101
二维勒让德多项式多项式的第一项为L0(u′,v′)=1,其方程系数AL0能够表示齿距偏差,第二项为L1(u′,v′)=u′,其方程系数AL1能够表示齿廓斜率,第三项为L2(u′,v′)=v′,其方程系数AL2能够于表示螺旋线斜率,第四项为
Figure BDA0003148498570000102
其方程系数AL3能够表示齿廓凸度,第五项为L4(u′,v′)=u′v′,其方程系数AL4能够表示轮齿鼓度,第六项为/>
Figure BDA0003148498570000103
其方程系数AL5能够表示齿面扭曲;即系数AL0、AL1、AL2、AL3、AL4和AL5可分别用于评价齿距偏差、齿廓斜率、螺旋线斜率、齿廓凸度、轮齿鼓度和齿面扭曲;计算得到的dL即为待检测齿轮的齿面上对应的齿面点或齿面线处的偏差值d。
上述L0(u′,v′)=1为数学领域中多项式的标准表达方式,含义为多项式第一项的值为1,后续表达式均与此表达含义相同,在此不再赘述。
本发明中,考虑到齿面可等效为齿面坐标系下的矩形表面,二维勒让德多项式和二维切比雪夫多项式本身为定义在[-1,1]*[-1,1]矩形区域上的多项式,故以二维勒让德多项式或二维切比雪夫多项式逼近齿面偏差方程。由于逼近时计算的是齿面偏差点云数据,相比于单个齿面点或齿面线的评价而言,所涉及的齿面数据更多更全面,能够更准确地表征齿面偏差。同时,二维勒让德多项式和二维切比雪夫多项式的前六项与齿面偏差具有相似特征,且多项式各项之间的正交特性使各项系数相互独立无串扰,有利于消除偶然因素造成的干扰。因此,本发明中选用二维勒让德多项式和/或二维切比雪夫多项式的前6项作为基函数,对齿面偏差点云数据进行最小二乘曲面拟合,通过拟合曲面的前六项各项系数来评价齿距偏差、齿廓斜率、螺旋线斜率、齿廓凸度、轮齿鼓度和齿面扭曲;同时根据得到偏差曲面方程dc和dL,即可求得待检测齿轮的齿面上任意齿面点或齿面线处的偏差值d。
本发明中:
1.可直接单独采用二维切比雪夫多项式拟合出偏差曲面方程dc,并利用偏差曲面方程dc中各项的对应系数AC0、AC1、AC2、AC3、AC4和AC5来评价齿距偏差、齿廓斜率、螺旋线斜率、齿廓凸度、轮齿鼓度和齿面扭曲。
2.可直接单独采用二维勒让德多项式拟合出偏差曲面方程dL,并利用偏差曲面方程dL中各项的对应系数AL0、AL1、AL2、AL3、AL4和AL5来评价齿距偏差、齿廓斜率、螺旋线斜率、齿廓凸度、轮齿鼓度和齿面扭曲。
3.由于二维切比雪夫多项式在表现偏差曲面的峰值以及细节处更优异,即系数AC3和AC5比系数AL3和AL5能够更为准确的评价齿廓凸度和齿面扭曲;而二维勒让德多项式则在表现曲面的整体趋势及轮廓上更优异,即AL0、AL1、AL2和AL4比AC0、AC1、AC2和AC4能够更为准确地评价齿距偏差、齿廓斜率、螺旋线斜率和轮齿鼓度。因此在需要准确评价齿廓凸度和齿面扭曲时,利用二维切比雪夫多项式拟合出偏差曲面方程dc更为准确;而在需要准确评价齿距偏差、齿廓斜率、螺旋线斜率和轮齿鼓度时,利用二维勒让德多项式拟合出偏差曲面方程dl更为准确。
因此,可以根据对于不同种类的齿面偏差(齿距偏差、齿廓斜率、螺旋线斜率、齿廓凸度、轮齿鼓度和齿面扭曲)的评价需求,同时使用二维勒让德多项式和/或二维切比雪夫多项式拟合出偏差曲面方程,在需要更准确的拟合曲面峰值及细节时选用通过二维切比雪夫多项式拟合的偏差曲面方程dc,作为待检测齿轮的偏差曲面方程;在需要更准确的拟合曲面整体趋势及轮廓时选用通过二维勒让德多项式拟合的偏差曲面方程dL,作为待检测齿轮的偏差曲面方程;并选择对应的系数用于得到更为准确的不同种类的齿面偏差评价结果。
实施例1
单独采用二维切比雪夫多项式作为基函数计算待检测齿轮的偏差曲面方程时,设齿面偏差点云数据中共有s个齿面偏差点数据,用最小二乘法求解偏差曲面方程时相当于对偏差曲面方程校正了s次,即得到了经s次校正后的最佳近似曲面。
下式为使用最小二乘法建立的方程组,由该方程组求解出偏差曲面方程dC
Figure BDA0003148498570000121
通过二维切比雪夫多项式求出定义在u′-v′-d齿面坐标系下的偏差曲面方程为:
dC=AC0*1+AC1*u'+Ac2*v'+AC3*(2u'2-1)+AC4*u'v'+AC5*(2v'2-1);
具体的偏差项和系数间的比例关系为:齿距偏差等于AC0,齿廓斜率等于2AC1,螺旋线斜率等于-2AC2,齿廓凸度等于2AC3,轮齿鼓度等于-4AC4,齿面扭曲等于-2AC5,计算得到的偏差值dc即为待检测齿轮的齿面上点(u,v)归一化为(u′,v′)所对应的偏差值d。
实施例2
单独采用二维勒让德多项式作为基函数计算待检测齿轮的偏差曲面方程时,设齿面偏差点云数据中共有s个齿面偏差点数据,用最小二乘法求解偏差曲面方程时相当于对偏差曲面方程校正了s次,即得到了经s次校正后的最佳近似曲面。
下式为使用最小二乘法建立的方程组,由该方程组求解出偏差曲面方程dL
Figure BDA0003148498570000131
通过二维勒让德多项式所求dL为定义在u′-v′-d齿面坐标系下的偏差曲面方程为:
Figure BDA0003148498570000132
具体的偏差项于系数间的比例关系为:齿距偏差等于AL0,齿廓斜率等于AL1,螺旋线斜率等于AL2,齿廓凸度等于1.5AL3,轮齿鼓度等于2AL4,齿面扭曲等于1.5AL5;计算得到的偏差值dL即为待检测齿轮的齿面上点(u,v)归一化为(u′,v′)所对应的偏差值d。
实施例3
分别采用二维勒让德多项式和二维切比雪夫多项式作为基函数得到对应的偏差曲面方程,然后根据对于不同种类的齿面偏差的评价需求,选择对应的偏差曲面方程的系数得到更为准确的不同种类的齿面偏差评价结果。
设齿面偏差点云数据中共有s个齿面偏差点数据,用最小二乘法求解偏差曲面方程时相当于对偏差曲面方程校正了s次,即得到了经s次校正后的最佳近似曲面。
采用二维切比雪夫多项式作为基函数得到的偏差曲面方程分别为:
dC=AC0*1+AC1*u'+AC2*v'+AC3*(2u'2-1)+AC4*u'v'+AC5*(2v'2-1);
采用二维切勒让德多项式作为基函数得到的偏差曲面方程为:
Figure BDA0003148498570000133
在需要准确评价齿廓凸度和齿面扭曲时,利用二维切比雪夫多项式拟合出的偏差曲面方程dc中第四项和第六项的系数进行求解,齿廓凸度等于2AC3,齿面扭曲等于-2AC5;在需要准确评价齿距偏差、齿廓斜率、螺旋线斜率和轮齿鼓度时,利用二维勒让德多项式拟合出的偏差曲面方程dL中第一项、第二项、第三项和第五项的系数进行求解,齿距偏差等于AL0,齿廓斜率等于AL1,螺旋线斜率等于AL2,轮齿鼓度等于2AL4
在需要更准确的拟合曲面峰值及细节时选用通过二维切比雪夫多项式拟合的偏差曲面方程dc,作为待检测齿轮的最终偏差曲面方程;在需要更准确的拟合曲面整体趋势及轮廓时选用通过二维勒让德多项式拟合的偏差曲面方程dL,作为待检测齿轮的最终偏差曲面方程。
E:齿轮测量云终端将用户输入的待评测齿轮上任意一个齿面点(u,v)的坐标或任意一条齿面线u=a的坐标传输至齿轮测量云平台,齿轮测量云平台首先按照步骤C中的方法,对用户输入的齿面点或齿面线的坐标进行坐标转化处理及归一化处理为对应(u′,v′)或u′;然后利用坐标转化处理后的坐标数据和步骤D中得到的待检测齿轮的偏差曲面方程,求出该齿面点或齿面线所对应位置处的偏差值d;齿轮测量云平台还依据待检测齿轮的偏差曲面方程中的系数求得待评测齿轮的齿距偏差、齿廓斜率、螺旋线斜率、齿廓凸度、轮齿鼓度和齿面扭曲,然后齿轮测量云平台结合用户选定的齿轮精度参数和齿轮测量评价标准,对待检测齿轮进行合格或不合格的评价;
齿轮精度参数为用户所选择的现有的齿轮评价精度标准,现有齿轮评测中,齿轮精度等级从0级至12级共分为13个精度等级,齿轮精度等级越高允许的偏差值越大,具体不同精度等级所允许的偏差值可通过查阅《机械设计手册》得知。本发明中使用的齿轮为6级精度齿轮,对于计算得出的偏差值只需要与《机械设计手册》中给出的6级精度齿轮允许的公差带相比较,偏差值在允许的公差带内即为合格。
F:齿轮测量云平台将计算得到的待评测齿轮的齿距偏差、齿廓斜率、螺旋线斜率、齿廓凸度、轮齿鼓度和齿面扭曲,以及用户输入的齿面点或齿面线所对应位置处的偏差值d和评价结果,传输至齿轮测量云终端并反馈给用户。

Claims (7)

1.一种基于云平台的齿轮齿面评价方法,其特征在于,依次包括以下具体步骤:
A:用户通过齿轮测量云终端,向齿轮测量云平台提供待检测齿轮的齿面坐标点云数据,齿面坐标点云数据包括待检测齿轮上若干个齿面点的坐标数据,每个齿面点的坐标数据中均包含该齿面点的三维坐标值;
其中,齿轮测量云终端即为用户所使用的本地客户端,用于与齿轮测量云平台实现数据交互,为用户提供齿轮齿面评价服务和齿轮齿面评价结果展示服务;齿轮测量云平台即为云上服务器端,用于对用户通过齿轮测量云终端提供的待检测齿轮进行计算及评价,并将评价结果反馈至齿轮测量云终端;
B:对于齿面坐标点云数据中的每一个齿面点的坐标值数据,齿轮测量云平台分别通过下述公式计算各个齿面点对应的齿面偏差值,然后将齿面坐标点云数据中所有的齿面点所对应的齿面偏差值组合构成齿面偏差点云数据;
Figure FDA0003922713390000011
其中,dlot表示齿面点对应的齿面偏差值,下角标lot表示根据铅垂线方向,K表示常数系数,ρ表示极径,ηb表示基本空间宽度半角,β表示螺旋角,rb表示基圆半径,X表示变位系数,αn表示齿轮法向压力角,Z表示齿数,齿数中内齿轮为负值,外齿轮为正值,inv表示渐开线函数,αt表示齿轮端面压力角,x、y和z分别表示齿面点在三维笛卡尔坐标系下在x轴、y轴和z轴上的坐标;
C:齿轮测量云平台对齿面坐标点云数据中的每一个齿面点的坐标值数据和齿面偏差值进行坐标转化处理,在从三维笛卡尔坐标系下转化到u-v-d齿面坐标系下后进行归一化处理,最终得到归一化处理后的u′-v′-d齿面坐标系下的齿面偏差点云数据;
D:齿轮测量云平台使用二维勒让德多项式和/或二维切比雪夫多项式的前6项作为基函数,利用归一化处理后的u′-v′-d齿面坐标系下的齿面偏差点云数据,通过最小二乘法拟合得出的待检测齿轮的偏差曲面方程;
其中,利用二维切比雪夫多项式的前6项作为基函数得到的待检测齿轮的偏差曲面方程为:
dC=AC0*1+AC1*u'+Ac2*v'+AC3*(2u'2-1)+AC4*u'v'+AC5*(2v'2-1);
方程系数AC0、AC1、AC2、AC3、AC4和AC5分别用于评价齿距偏差、齿廓斜率、螺旋线斜率、齿廓凸度、轮齿鼓度和齿面扭曲,计算得到的dC即为待检测齿轮的齿面上对应的齿面点或齿面线处的偏差值d;
利用二维勒让德多项式的前6项作为基函数得到的待检测齿轮的偏差曲面方程为:
Figure FDA0003922713390000021
系数AL0、AL1、AL2、AL3、AL4和AL5分别用于评价齿距偏差、齿廓斜率、螺旋线斜率、齿廓凸度、轮齿鼓度和齿面扭曲;计算得到的dL即为待检测齿轮的齿面上对应的齿面点或齿面线处的偏差值d;
E:齿轮测量云终端将用户输入的待评测齿轮上任意一个齿面点(u,v)的坐标或任意一条齿面线u=a的坐标传输至齿轮测量云平台,齿轮测量云平台首先按照步骤C中的方法,对用户输入的齿面点或齿面线的坐标进行坐标转化处理及归一化处理为对应(u′,v′)或u′;然后利用坐标转化处理后的坐标数据和步骤D中得到的待检测齿轮的偏差曲面方程,求出该齿面点或齿面线所对应位置处的偏差值d;齿轮测量云平台还依据待检测齿轮的偏差曲面方程中的系数求得待评测齿轮的齿距偏差、齿廓斜率、螺旋线斜率、齿廓凸度、轮齿鼓度和齿面扭曲,然后齿轮测量云平台结合用户选定的齿轮精度参数和齿轮测量评价标准,对待检测齿轮进行合格或不合格的评价;
F:齿轮测量云平台将计算得到的待评测齿轮的齿距偏差、齿廓斜率、螺旋线斜率、齿廓凸度、轮齿鼓度和齿面扭曲,以及用户输入的齿面点或齿面线所对应位置处的偏差值d和评价结果,传输至齿轮测量云终端并反馈给用户。
2.根据权利要求1所述的基于云平台的齿轮齿面评价方法,其特征在于,所述的步骤C中:
在进行坐标转化处理时,齿轮测量云平台首先将齿面坐标点云数据中的每一个齿面点的坐标值数据和齿面偏差值,从三维笛卡尔坐标系下转化到u-v-d齿面坐标系下,然后再对u-v-d齿面坐标系下的齿面点的坐标值(u,v,d)进行归一化处理得到坐标值(u′,v′,d),其中平面u′v′的取值范围在[-1,1]*[-1,1]内;
坐标转化公式如下:
Figure FDA0003922713390000031
Figure FDA0003922713390000032
d=dlot
其中,u表示沿齿廓方向的坐标值,φm表示极角,Λnom表示渐开线起始角,rb表示基圆半径,ξnom表示齿廓滚动角,θk表示展角,v表示沿齿宽方向的坐标值,Z表示齿数,β表示螺旋角,齿面点或齿面线所对应位置处的偏差值d在坐标转化过程中也表示垂直于uv平面方向的坐标值,dlot表示齿面点对应的齿面偏差值;渐开线起始角Λnom和齿廓滚动角ξnom均为待检测齿轮的基本参数。
3.根据权利要求2所述的基于云平台的齿轮齿面评价方法,其特征在于,所述的步骤C中:
在进行归一化处理时,对于齿面坐标系中任意一点(u,v,d),按照如下方式转化为(u′,v′,d),最终转化为在[-1,1]*[-1,1]矩形区域上的齿面点的坐标值(u′,v′,d);
Figure FDA0003922713390000033
Figure FDA0003922713390000034
d=dlot
其中,Lα表示u的取值范围,Lβ表示v的取值范围。
4.根据权利要求1所述的基于云平台的齿轮齿面评价方法,其特征在于,所述的步骤D中:
单独采用二维切比雪夫多项式作为基函数计算待检测齿轮的偏差曲面方程时,设齿面偏差点云数据中共有s个齿面偏差点数据;
用最小二乘法建立方程组,由该方程组求解出偏差曲面方程dC
Figure FDA0003922713390000041
通过二维切比雪夫多项式求出定义在u′-v′-d齿面坐标系下的偏差曲面方程为:
dC=AC0*1+AC1*u'+Ac2*v'+AC3*(2u'2-1)+AC4*u'v'+AC5*(2v'2-1);
齿距偏差等于AC0,齿廓斜率等于2AC1,螺旋线斜率等于-2AC2,齿廓凸度等于2AC3,轮齿鼓度等于-4AC4,齿面扭曲等于-2AC5,计算得到的偏差值dC即为待检测齿轮的齿面上点(u,v)归一化为(u′,v′)所对应的偏差值d。
5.根据权利要求1所述的基于云平台的齿轮齿面评价方法,其特征在于,所述的步骤D中:
单独采用二维勒让德多项式作为基函数计算待检测齿轮的偏差曲面方程时,设齿面偏差点云数据中共有s个齿面偏差点数据;
用最小二乘法建立方程组,由该方程组求解出偏差曲面方程dL
Figure FDA0003922713390000042
通过二维勒让德多项式所求dL为定义在u′-v′-d齿面坐标系下的偏差曲面方程为:
Figure FDA0003922713390000043
齿距偏差等于AL0,齿廓斜率等于AL1,螺旋线斜率等于AL2,齿廓凸度等于1.5AL3,轮齿鼓度等于2AL4,齿面扭曲等于1.5AL5;计算得到的偏差值dL即为待检测齿轮的齿面上点(u,v)归一化为(u′,v′)所对应的偏差值d。
6.根据权利要求1所述的基于云平台的齿轮齿面评价方法,其特征在于,所述的步骤D中:
分别采用二维勒让德多项式和二维切比雪夫多项式作为基函数得到对应的偏差曲面方程,然后根据对于不同种类的齿面偏差的评价需求,选择对应的偏差曲面方程的系数得到更为准确的不同种类的齿面偏差评价结果;
设齿面偏差点云数据中共有s个齿面偏差点数据;
采用二维切比雪夫多项式作为基函数得到的偏差曲面方程分别为:
dC=AC0*1+AC1*u'+AC2*v'+AC3*(2u'2-1)+AC4*u'v'+AC5*(2v'2-1);
采用二维切勒让德多项式作为基函数得到的偏差曲面方程为:
Figure FDA0003922713390000051
在需要准确评价齿廓凸度和齿面扭曲时,利用二维切比雪夫多项式拟合出的偏差曲面方程dC中第四项和第六项的系数进行求解,齿廓凸度等于2AC3,齿面扭曲等于-2AC5;在需要准确评价齿距偏差、齿廓斜率、螺旋线斜率和轮齿鼓度时,利用二维勒让德多项式拟合出的偏差曲面方程dL中第一项、第二项、第三项和第五项的系数进行求解,齿距偏差等于AL0,齿廓斜率等于AL1,螺旋线斜率等于AL2,轮齿鼓度等于2AL4
在需要更准确的拟合曲面峰值及细节时选用通过二维切比雪夫多项式拟合的偏差曲面方程dC,作为待检测齿轮的最终偏差曲面方程;在需要更准确的拟合曲面整体趋势及轮廓时选用通过二维勒让德多项式拟合的偏差曲面方程dL,作为待检测齿轮的最终偏差曲面方程。
7.根据权利要求1所述的基于云平台的齿轮齿面评价方法,其特征在于:所述的齿轮精度参数采用6级精度齿轮参数。
CN202110758977.XA 2021-07-05 2021-07-05 一种基于云平台的齿轮齿面评价方法 Active CN113487180B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110758977.XA CN113487180B (zh) 2021-07-05 2021-07-05 一种基于云平台的齿轮齿面评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110758977.XA CN113487180B (zh) 2021-07-05 2021-07-05 一种基于云平台的齿轮齿面评价方法

Publications (2)

Publication Number Publication Date
CN113487180A CN113487180A (zh) 2021-10-08
CN113487180B true CN113487180B (zh) 2023-03-24

Family

ID=77940296

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110758977.XA Active CN113487180B (zh) 2021-07-05 2021-07-05 一种基于云平台的齿轮齿面评价方法

Country Status (1)

Country Link
CN (1) CN113487180B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115222300B (zh) * 2022-09-20 2023-02-28 珠海翔翼航空技术有限公司 基于智能评价算法的航空模拟器复训时长分配***和方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010151657A (ja) * 2008-12-25 2010-07-08 Fuji Heavy Ind Ltd 歯車対の評価装置及びこの評価装置を用いて最適化された歯車対
CN103344210A (zh) * 2013-07-22 2013-10-09 北京工业大学 一种齿轮误差多自由度评定方法
CN105512442A (zh) * 2016-01-27 2016-04-20 北京工业大学 基于统计分析的齿轮精度评价方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011070927A1 (ja) * 2009-12-11 2011-06-16 株式会社トプコン 点群データ処理装置、点群データ処理方法、および点群データ処理プログラム
KR101399792B1 (ko) * 2013-02-25 2014-05-27 주식회사 만도 유성 기어 장치
CN107392351A (zh) * 2017-06-16 2017-11-24 山东科技大学 一种用于云制造环境下的齿轮箱设计方案的优选方法
CN109632293B (zh) * 2018-12-29 2020-07-17 河南理工大学 一种基于云平台的齿轮云测量***及齿轮云测量方法
CN110146033B (zh) * 2019-06-04 2020-08-07 西安工业大学 基于点云数据的接触线-展开线齿轮齿面误差表述方法
CN112580160B (zh) * 2020-12-03 2022-10-14 南京工业大学 一种用于成形磨齿机在机测量***标定方法
CN112903288B (zh) * 2021-01-25 2022-06-21 北京工业大学 一种齿轮三维误差的特征线统一表征方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010151657A (ja) * 2008-12-25 2010-07-08 Fuji Heavy Ind Ltd 歯車対の評価装置及びこの評価装置を用いて最適化された歯車対
CN103344210A (zh) * 2013-07-22 2013-10-09 北京工业大学 一种齿轮误差多自由度评定方法
CN105512442A (zh) * 2016-01-27 2016-04-20 北京工业大学 基于统计分析的齿轮精度评价方法

Also Published As

Publication number Publication date
CN113487180A (zh) 2021-10-08

Similar Documents

Publication Publication Date Title
CN113487180B (zh) 一种基于云平台的齿轮齿面评价方法
CN113496356B (zh) 一种齿轮数字化测量评价云***及测量评价方法
CN110849291B (zh) 大型弯管弯曲半径检测方法
CN111666643B (zh) 一种用于确定复杂齿面接触性能的方法
CN112903288A (zh) 一种齿轮三维误差的特征线统一表征方法
CN112539721A (zh) 机器人用减速机的三曲柄摆线轮关键加工误差测量方法
CN106980724A (zh) 一种基于复杂齿面网格规划原则的齿面精确模型截面放样构建方法
CN108645301A (zh) 一种直齿轮齿面偏差在机测量方法
CN113434817B (zh) 一种齿轮单项拓扑误差图谱的分析方法
CN115221655B (zh) 一种斜齿非圆齿轮的传动精度评价方法
CN109887079B (zh) 螺旋锥齿轮三维建模方法
JPH02649B2 (zh)
CN108731616B (zh) 一种基于云模型的螺旋锥齿轮齿面测量点自适应分布方法
CN112539722B (zh) 机器人用减速机的双曲柄摆线轮关键加工误差测量方法
CN104125873A (zh) 错口计算***
CN109886949B (zh) 一种基于机器视觉的直齿圆柱齿轮多参数评定方法
CN115270407A (zh) 抗剪强度参数计算方法及其评估方法
KR20030050469A (ko) 스파이럴 베벨기어의 형상오차 측정방법
CN109977530B (zh) 直齿面齿轮三维建模方法
CN110222382B (zh) 一种飞机坐标系优化拟合方法
CN112729206A (zh) 非渐开线车齿刀齿形检测方法
CN113566772A (zh) 一种基于坐标测量的局部齿面定位方法
CN112923890B (zh) 一种花键齿形误差测量评价方法
CN110457820B (zh) 一种确定齿轮副接触椭圆的方法
CN107992688A (zh) 机车牵引齿轮自动设计设备及自动设计方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant