CN110457820B - 一种确定齿轮副接触椭圆的方法 - Google Patents
一种确定齿轮副接触椭圆的方法 Download PDFInfo
- Publication number
- CN110457820B CN110457820B CN201910742874.7A CN201910742874A CN110457820B CN 110457820 B CN110457820 B CN 110457820B CN 201910742874 A CN201910742874 A CN 201910742874A CN 110457820 B CN110457820 B CN 110457820B
- Authority
- CN
- China
- Prior art keywords
- tooth surface
- axis
- contact ellipse
- reference circle
- points
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 25
- 239000013598 vector Substances 0.000 claims abstract description 35
- 238000004364 calculation method Methods 0.000 claims abstract description 25
- 238000013507 mapping Methods 0.000 claims abstract description 4
- 230000013011 mating Effects 0.000 claims description 12
- 238000004422 calculation algorithm Methods 0.000 description 6
- 238000013461 design Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- 238000004458 analytical method Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Landscapes
- Image Analysis (AREA)
- Gears, Cams (AREA)
Abstract
本发明公开了一种确定齿轮副接触椭圆的方法,首先,根据配合齿面啮合点的坐标计算配合齿面啮合点处的公共法向量;然后,将配合齿面上的点映射到以啮合点处两齿面公共法向量为坐标轴啮合点为坐标原点建立的计算坐标系中;接着,在过啮合点并且垂直于啮合点公共法向量的平面内选定参考圆;最后,计算圆周上各点处两配合齿面相应点间的距离,将参考圆圆周上对应的齿面距离最大及最小的4个点确定为接触椭圆的长轴和短轴点,进一步确定接触椭圆长轴和短轴的长度以及方向向量。该方法在没有配合齿面啮合点处两曲面的主曲率、主方向等信息的情况下就能够确定齿面接触椭圆;克服了传统方法不能准确确定复杂修形齿面接触椭圆和鲁棒性差的难题。
Description
技术领域
本发明属于齿轮设计领域,特别是齿轮复杂修形设计领域,涉及一种确定齿轮副接触椭圆的方法。
技术背景
齿轮轮齿接触分析在齿轮专业领域内具有重要作用,广泛应用于高度复杂的弧齿锥齿轮、复杂修形齿轮的设计、加工等方面。齿面接触椭圆是评价齿面接触状况的重要指标,对于齿面设计参数选取具有重要指导意义。
目前,在齿轮齿面设计过程中,齿面接触椭圆的计算主要是基于曲面理论,其要点是通过计算配合齿面啮合点处齿面变量的1阶、2阶以及2阶混合偏导数确定配合齿面的啮合点处主曲率和主方向,最终确定接触椭圆的长半轴和短半轴的方向以及长度。
因此,现有的计算接触椭圆的方法基于曲面理论,尤其依赖配合齿面的解析表达式,通过复杂的偏导计算,确定接触椭圆。这种方法一方面建模复杂,整个计算过程涉及复杂的公式推导和复杂的偏导计算,对一般轮齿设计人员来说,掌握这一套理论和计算方法有很大难度;另一方面,计算过程中计算参数相差数个数量级,算法鲁棒性差。
此外,由于现有方法无法确定复杂修形齿轮或者没有解析式的齿面的接触椭圆,导致在设计复杂修形齿面时,需要通过齿面三维接触仿真甚至是齿轮副对滚实验来确定不同设计参数下齿面的接触状态(接触椭圆),严重降低了复杂修形齿面等的设计效率。
发明内容
针对现有方法存在的缺陷,本发明提供了一种能够确定复杂修形齿面和无解析方程齿面啮合传动过程中的接触椭圆,并且具有高鲁棒性。
为了解决上述技术问题,本发明采用了如下的技术方案:
一种确定齿轮副接触椭圆的方法,首先,根据配合齿面啮合点的坐标,计算配合齿面啮合点处的公共法向量;然后,将配合齿面的点映射到以啮合点处公共法向量为坐标轴啮合点为坐标原点建立的计算坐标系中;接着,在过啮合点并且垂直于啮合点公共法向量的平面内选定半径一定的参考圆,参考圆的半径可以任意选取;最后,计算参考圆圆周上各点处两配合齿面相应点间的距离,通过逼近算法寻找圆周上两齿面距离最大和最小的4个点,通过所得的4个点确定接触椭圆长轴、短轴的长度以及方向向量。
作为本发明的一种优选方案,通过确定参考圆圆周上两齿面间距离最大和最小的4个极值点,根据所得的4个极值点与计算坐标系原点确定接触椭圆的长轴和短轴的方向;
作为本发明的另一种优选方案,通过确定参考圆圆周上两齿面间距离最大和最小的4个极值点,根据所得的4个极值点与计算坐标系原点确定接触椭圆的长轴和短轴的长度;
式中:LenLA为接触椭圆长轴的长度;
LenSA为接触椭圆短轴的长度;
δ为接触深度;
r为参考圆半径;
Dismin,Dismax为参考圆上配合齿面间最小和最大距离。
作为本发明的一种改进方案,所述的通过啮合点确定啮合平面公共法向量,对于有齿面解析式的,对齿面解析式进行偏微分得到;对于复杂修形齿轮或没有明确解析式的齿面,通过从啮合点确定齿面在啮合点处的任意两切向量得到。
作为本发明的另一种改进方案,所述的通过逼近算法确定圆周上两齿面距离最大和最小的2个点。
与现有技术相比,本发明具有如下技术效果:克服了传统解析法建模复杂,计算困难,算法鲁棒性低的缺点;并且打破了传统解析法无法准确计算复杂修形齿面和没有解析表达式的齿面的接触椭圆的技术瓶颈;另外,具有过程简单、鲁棒性高等优点。
附图说明
图1为计算复杂修形齿面/无解析式齿面啮合点处公共法向量的流程图;
图2为配合齿面啮合模型图;
图3为计算平面参考圆示意图;
图4为参考圆圆周上任意点对配合齿面距离示意图;
图5为基于参考圆的,搜寻圆周上的对应接触椭圆长半轴/短半轴端点算法图。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细地描述。
一种确定齿轮副接触椭圆的方法,该方法包括如下步骤:
1、根据配合齿面啮合点,计算啮合点处两齿面的公共法向量,以确定配合齿面的公切面;对于有明确的解析式的齿面可以通过齿面解析式进行偏微分求得:
对于复杂修形齿面或没有明确解析式的齿面,配合齿面公共法向量,可按图1求得。
式中:f(u,θ)为配合齿面中任意齿面的方程,
u、θ为齿面的方程的自变量。
2、在计算得到配合齿面啮合点处的公共法向量后,以该公共法向量为一坐标轴、以啮合点为原点建立计算坐标系Sc;并将两齿面方程或者齿面点位置矢量从参考坐标系Sf映射到计算坐标系中,如图2所示。计算坐标系Sc与参考坐标系Sf间的转换关系如下:
式中:Mcf为参考坐标系Sf转换到Sc的转换矩阵;
θx,θy,θz分别是参考坐标系Sf转换到Sc过程中绕各坐标轴的旋转角度;
3、过坐标系Sc原点建立垂直于公共法向量的计算平面,在计算平面上确定一半径任意的参考圆,如图3所示。
4、如图4、图5所示,在参考圆圆周上任意选定一起始点;选定一合适步长;计算圆周上的点对应的两齿面上的点到计算平面的距离之和;在搜寻过程中,一旦发现圆周上点对应的齿面间的距离的变化趋势发生改变,则以当前点与前两圆周点为区间(接触椭圆的长半轴/短半轴存在于该区间),然后,在区间内用逼近算法寻找到满足精度要求的接触椭圆的长半轴/短半轴点。以此类推,直到找到长半轴/短半轴对应的所有点。
5、在找到参考圆轴上的对应的配合齿面距离最大和最小的4个极值点后,根据下式计算接触椭圆的长半轴、短半轴的方向以及长度。
其中:为计算坐标系中接触椭圆短轴的方向向量;为为计算坐标系中接触椭圆长轴的方向向量;为参考圆圆周上配合齿面距离最大和最小点的位置向量;为相应向量的模;LenLA为接触椭圆长轴的长度;LenSA为接触椭圆短轴的长度;δ为接触深度;r为参考圆半径;Dismin,Dismax为参考圆上配合齿面间最小和最大距离。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (3)
1.一种确定齿轮副接触椭圆的方法,其特征在于:首先,根据配合齿面啮合点的坐标确定配合齿面啮合点处的公共法向量;然后,以接触点处两啮合齿面的法向量为一坐标轴,啮合点为原点建立计算笛卡尔坐标系,并将配合齿面的点映射到新建的计算坐标系中;接着,在过啮合点并且垂直于啮合点公共法向量的平面内选定参考圆;最后,计算参考圆圆上的点对应的配合齿面间的距离,以参考圆圆周上对应的齿面间距离最大和最小的4个极值点作为接触椭圆长半轴和短半轴的点,利用所确定的4个极值点确定接触椭圆的长半轴、短半轴方向向量以及长度;
通过确定参考圆圆周上两齿面间距离最大和最小的4个极值点,根据所得的4个极值点与计算坐标系原点确定接触椭圆的长轴和短轴的方向;
通过确定参考圆圆周上两齿面间距离最大和最小的4个极值点,根据所得的4个极值点与计算坐标系原点确定接触椭圆的长轴和短轴的长度;
式中:LenLA为接触椭圆长轴的长度;
LenSA为接触椭圆短轴的长度;
δ为接触深度;
r为参考圆半径;
Dismin,Dismax为参考圆上配合齿面间最小和最大距离。
2.根据权利要求1所述的一种确定齿轮副接触椭圆的方法,其特征在于,仅需要计算参考圆圆周上各点所对应的配合齿面相应点间的距离,参考圆的半径可以任意选取。
3.根据权利要求1所述的一种确定齿轮副接触椭圆的方法,其特征在于,啮合面上的点需要映射到计算坐标系中,计算坐标系的原点是配合齿面的啮合点,计算坐标系一坐标轴的方向向量与配合齿面的啮合点处公共法向量共线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910742874.7A CN110457820B (zh) | 2019-08-12 | 2019-08-12 | 一种确定齿轮副接触椭圆的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910742874.7A CN110457820B (zh) | 2019-08-12 | 2019-08-12 | 一种确定齿轮副接触椭圆的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110457820A CN110457820A (zh) | 2019-11-15 |
CN110457820B true CN110457820B (zh) | 2022-11-18 |
Family
ID=68486196
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910742874.7A Active CN110457820B (zh) | 2019-08-12 | 2019-08-12 | 一种确定齿轮副接触椭圆的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110457820B (zh) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3824873A (en) * | 1970-12-18 | 1974-07-23 | W Rouverol | Rolling contact gearing |
GB1500793A (en) * | 1974-03-22 | 1978-02-08 | Rouverol W | Toothed gearing |
CN101476614A (zh) * | 2009-01-20 | 2009-07-08 | 河南科技大学 | 直线-渐开线齿廓内啮合圆柱齿轮副 |
CN101876589A (zh) * | 2009-04-28 | 2010-11-03 | 上海工程技术大学 | 一种点啮合齿面性能的检测方法 |
CN101937482A (zh) * | 2010-09-08 | 2011-01-05 | 北京航空航天大学 | 面齿轮插齿刀齿形设计方法 |
CN104408220A (zh) * | 2014-10-08 | 2015-03-11 | 西北工业大学 | 一种改进的轮齿加载接触分析方法 |
CN106570278A (zh) * | 2016-11-08 | 2017-04-19 | 江苏大学 | 融合齿距偏差的圆柱齿轮啮合刚度计算方法 |
CN107917176A (zh) * | 2017-10-20 | 2018-04-17 | 淮阴工学院 | 一种球面渐开线弧齿锥齿轮的修形方法 |
CN107992698A (zh) * | 2017-12-13 | 2018-05-04 | 长安大学 | 一种螺旋锥齿轮齿面接触分析方法 |
CN109376455A (zh) * | 2018-11-12 | 2019-02-22 | 中南大学 | 一种弧齿锥齿轮载荷弹性接触变形数值计算方法 |
CN109376456A (zh) * | 2018-11-12 | 2019-02-22 | 中南大学 | 弧齿锥齿轮含安装误差的齿面载荷接触性能数值计算方法 |
CN109992877A (zh) * | 2019-03-28 | 2019-07-09 | 长安大学 | 一种齿轮副几何接触分析的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20170011142A1 (en) * | 2015-07-07 | 2017-01-12 | Gert Heirman | Gear Contact Modeling |
-
2019
- 2019-08-12 CN CN201910742874.7A patent/CN110457820B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3824873A (en) * | 1970-12-18 | 1974-07-23 | W Rouverol | Rolling contact gearing |
GB1500793A (en) * | 1974-03-22 | 1978-02-08 | Rouverol W | Toothed gearing |
CN101476614A (zh) * | 2009-01-20 | 2009-07-08 | 河南科技大学 | 直线-渐开线齿廓内啮合圆柱齿轮副 |
CN101876589A (zh) * | 2009-04-28 | 2010-11-03 | 上海工程技术大学 | 一种点啮合齿面性能的检测方法 |
CN101937482A (zh) * | 2010-09-08 | 2011-01-05 | 北京航空航天大学 | 面齿轮插齿刀齿形设计方法 |
CN104408220A (zh) * | 2014-10-08 | 2015-03-11 | 西北工业大学 | 一种改进的轮齿加载接触分析方法 |
CN106570278A (zh) * | 2016-11-08 | 2017-04-19 | 江苏大学 | 融合齿距偏差的圆柱齿轮啮合刚度计算方法 |
CN107917176A (zh) * | 2017-10-20 | 2018-04-17 | 淮阴工学院 | 一种球面渐开线弧齿锥齿轮的修形方法 |
CN107992698A (zh) * | 2017-12-13 | 2018-05-04 | 长安大学 | 一种螺旋锥齿轮齿面接触分析方法 |
CN109376455A (zh) * | 2018-11-12 | 2019-02-22 | 中南大学 | 一种弧齿锥齿轮载荷弹性接触变形数值计算方法 |
CN109376456A (zh) * | 2018-11-12 | 2019-02-22 | 中南大学 | 弧齿锥齿轮含安装误差的齿面载荷接触性能数值计算方法 |
CN109992877A (zh) * | 2019-03-28 | 2019-07-09 | 长安大学 | 一种齿轮副几何接触分析的方法 |
Non-Patent Citations (2)
Title |
---|
Contact characteristics of spherical gears;Li-chi chao et al.;《Mechanism and Machine Theory》;20081031;1317-1331 * |
螺旋锥齿轮多体多自由度非线性动力学研究;冯治恒;《中国博士学位论文全文数据库 工程科技II辑》;20120115;C029-3 * |
Also Published As
Publication number | Publication date |
---|---|
CN110457820A (zh) | 2019-11-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108564605B (zh) | 一种三维测量点云优化配准方法 | |
CN111080684B (zh) | 一种点邻域尺度差异描述的点云配准方法 | |
CN102393678B (zh) | 适用于五轴数控装置的轨迹平滑处理方法 | |
CN108062073B (zh) | 一种用于高质量加工的圆弧平滑压缩插补方法 | |
CN108446524B (zh) | 一种用于数控加工的弧齿锥齿轮的建模方法 | |
CN110849291B (zh) | 大型弯管弯曲半径检测方法 | |
CN108073138B (zh) | 适用于高速高精加工的椭圆弧平滑压缩插补算法 | |
CN104898556B (zh) | 由成形刀具确定被加工螺旋曲面廓形的点矢量包络法 | |
CN106958643B (zh) | 齿轮轮齿接触性能分析方法 | |
CN111238386B (zh) | 发动机管道最小间距测量方法 | |
CN106874580A (zh) | 一种基于点云数据的弯管模型重建方法 | |
CN111489432A (zh) | 一种基于点云数据的弯管重建及余量计算方法 | |
CN101196394A (zh) | 小段圆弧圆度的优化最小二乘评价方法 | |
CN111666643B (zh) | 一种用于确定复杂齿面接触性能的方法 | |
CN106774153A (zh) | 一种s形球头铣刀的轮廓度补偿方法 | |
CN112017293A (zh) | 一种关于圆钢管几何初始缺陷的测量方法 | |
CN109855587A (zh) | 基于轮廓度的二维凸轮类零件型面数字化测量方法 | |
CN103272901A (zh) | 一种导管实样取制方法 | |
CN110457820B (zh) | 一种确定齿轮副接触椭圆的方法 | |
CN106774154A (zh) | 一种基于密切面理论的空间曲线插补方法 | |
CN106292531B (zh) | 一种计算加工zn1蜗杆盘状成形刀具廓形边界的算法 | |
CN109389626B (zh) | 一种基于采样球扩散的复杂异形曲面点云配准方法 | |
CN110765621B (zh) | 一种nx渐开线花键建模方法 | |
Ye et al. | A Method of Binocular Laser 3-D Scanning Imaging for Reflective Workpieces | |
CN113487180B (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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20240426 Address after: 317699 No.1 Shengyuan Road, Mechanical and Electrical Industry Functional Zone, Yuhuan City, Wenzhou City, Zhejiang Province Patentee after: ZHEJIANG SHUANGHUAN DRIVELINE Co.,Ltd. Country or region after: China Address before: 400044 No. 174 Sha Jie street, Shapingba District, Chongqing Patentee before: Chongqing University Country or region before: China |