在某型动力伺服刀塔的传动系统设计中, 合理调整系统中传动齿轮的变位系数和根切量, 能够显著改善齿轮啮合对的传动性能, 进而降低换刀过程中刀盘的残余振动; 因而准确建立根切量和根切点位置与齿轮设计参数的对应关系, 并计算变位齿轮的根切量, 成为动力刀塔传动系统设计研究中的重要问题.本文给出一种较为简便的动力刀塔传动齿轮根切量的计算方法.
众多学者对齿轮根切量的计算问题进行过深入研究.陈奇等[1]在建立变位齿轮的齿根过渡曲线方程的基础上, 通过程序完成变位齿轮的全齿廓数据的精确计算, 并利用宏命令接口将数据输入到三维设计软件中, 但文中并未具体说明齿面渐开线与齿根过渡曲线交点的确定方法.齿轮根切量及根切点位置的确定方式主要有间接法和直接法.间接方法是先计算根切量, 再根据参数方程确定根切点位置.间接方法计算量较小, 但涉及变量较多, 求解精度不易控制.直接方法是计算参数曲线的交点得到根切点的位置坐标, 进而计算根切量.直接方法计算量较大, 便于模块化计算, 容易控制求解精度.本文采用直接方式计算齿轮根切量及根切点位置.
本文给出一种基于多项式变异粒子群优化的齿轮根切量计算方法.给出多项式变异粒子群优化(PMOPSO)算法的基本公式及算法参数的选取原则; 建立根切点位置和根切量的计算模型; 讨论计算模型解的情况与根切点位置的关系, 以及变位系数和模数对根切量的影响情况.
1 根切量的常用计算方法齿轮的根切量hu是齿轮齿廓根切点到基圆的径向距离.由图 1所示的几何关系, 给出关于θ和φ的方程, 其中M是根切点.
(1) |
如图 1所示, 计算根切点的初始切入量u:
式中:刀具直线切削刃齿顶高hs=m(hα*-χ), 这里一般有变位系数χ<0;m是模数; hα*是齿顶高系数; rb是基圆半径; α是分度圆压力角; λ是中间变量.
求解式(1) 中的方程可得到φ*, 但式(1) 中的超越方程求解困难.叶松林[2]给出了一种处理方法, 由式(1) 得
于是将式(1) 的方程转化为关于θ的一元非线性方程:
(2) |
其中
于是计算齿轮的根切量hu:
会田俊夫[3]还给出了hu的近似计算公式:
(3) |
其中λ由式(1) 确定.
由hu可得根切点到基圆圆心距离rM和根切点的压力角αM:
(4) |
于是由式(6) 得到根切点M的位置坐标.
叶松林方法(式(2))和会田俊夫公式(式(3))是传统的齿轮根切量计算方法, 其仅能确定齿轮根切量的值, 并不能给出根切点的位置坐标, 这显然无法满足动力刀塔传动齿轮设计中需要精确分析根切情况的需求.
2 PMOPSO算法多项式变异粒子群优化(PMOPSO)方法的基本思想是, 在粒子群演化过程中, 利用多项式变异算法对当前粒子进行扰动, 并以一定的概率接受劣化解.随着时间的推移, 劣化解的接受概率逐渐趋于零, PMOPSO算法也逐渐收敛于最优解.在搜索后期, 多项式变异算子产生的新点逐渐逼近当前点, 因而可提升PMOPSO算法的局部搜索精度.PMOPSO算法是对标准粒子群优化[4-5]方法的改进, 适用于含有复杂非线性约束的连续型变量的数学优化问题.
本文采用含有固定边界约束的多项式变异算子[6].已知当前点p和变异新点x满足x, p∈[xl, xu], 则定义扰动因子
于是得到扰动因子δ的计算公式[7]:
其中随机数u~U(0, 1), v=0.5(1-β)η+1, 常数β的计算公式是
粒子群优化算法设计中, 这里采用带有惯性权重ωt和收缩因子χ0的粒子群进化方程:
其中, 惯性权重采用线性递减的计算公式:
式中:t是进化时间; ωt是t时刻惯性权重的值; ωmin和ωmax是惯性权重的取值范围; tmax是进化时间的最大值.
PMOPSO方法中的算法参数和种群拓扑结构对优化设计问题的求解精度和计算复杂度的影响很大.在满足计算精度要求的条件下, 尽量降低种群规模和最大进化次数, 以降低计算复杂度.这里采用环形邻域拓扑结构, 个体的邻域半径约为种群规模的1/5.
3 数学模型如图 3所示, 以齿轮圆心为原点O, 齿轮轮齿中心线为y轴建立直角坐标系xOy.齿廓包含4段曲线:齿根圆弧曲线AB, 齿根过渡曲线BM, 齿面渐开线MD和齿顶圆弧曲线DE.齿根过渡曲线与齿面渐开线的交点是根切点, 并记作M.M点是齿根过渡曲线与齿面渐开线的过渡点.当齿根过渡曲线与齿面渐开线相交于M点时, 范成产生的渐开线齿面的一部分被刀具切去, 产生根切现象.当齿根过渡曲线与齿面渐开线相切时, M是切点, 而且无根切.M点又可称为过渡点.
文献[8]给出了图 3中齿根过渡曲线和齿面渐开线的参数方程.关于参数α′的齿根过渡曲线的参数方程是
(5) |
式中:x1和y1是齿根过渡曲线上点的坐标;φ是中间变量并与α′相关.齿根过渡曲线参数的取值范围是
关于参数αk的齿面渐开线的参数方程是
(6) |
其中x2和y2是齿面渐开线上点的坐标.齿面渐开线参数的取值范围是αk∈[αM, αa].αM是根切点M处的齿面渐开线压力角, αa是齿顶圆压力角.
式中:δ, rb是常数; γ=tanαk.
如图 3所示, M点是过渡曲线与齿面渐开线的交点.计算根切点M位置的问题可以转化为求齿根过渡曲线与齿面渐开线上最小距离点的问题, 这里采用PMOPSO算法求解.建立确定根切点M位置的数学模型:
(7) |
其中
以下讨论式(7) 中优化模型的解的数量与根切点位置的关系.定义根切临界齿数z0*:
基圆临界齿数z1*是
一般地, 对于ha*=1.0和c*=0.25可以得到z0*, z1*和变位系数χ的函数关系, 如图 4所示, 其中满足z≥z1*>z0*.如图 5所示, 根切量hu是根切点M到基圆的距离.根据齿轮齿数和临界齿数的关系, 可以分三种情况讨论.
当z<z0*时, 有rb>rf.如图 5a所示, 过渡曲线与齿面渐开线只有一个交点, 存在根切.此时式(7) 存在唯一解, 过渡点M是根切点.
当z0*≤z<z1*时, 有rb>rf.如图 5b所示, 过渡曲线与齿面渐开线相切, 无根切.过渡点M是切点, 其对应的过渡曲线参数是α′M=α.
当z≥z1*时, 有rb<rf.如图 5c所示, 过渡曲线与齿面渐开线相切, 无根切.过渡点M是切点, 其对应的过渡曲线参数是α′M=α.
式(7) 中, 若z<z0*则存在根切点, 且存在唯一解; 若z≥z0*,则不存在根切, 过渡点M是切点, 且有αM′=α.
综上所述, 当齿数小于根切临界齿数时, 式(7) 的解是唯一的, 过渡曲线与齿面渐开线有唯一交点; 否则, 两曲线相切, 且切点处的压力角等于齿轮分度圆压力角.因而, 式(7) 中的优化模型总是存在唯一解, 这给问题的求解提供了很大的便利.
4 计算实例本节以动力伺服刀塔传动系统中的某关键齿轮为研究对象, 其基本设计参数和PMOPSO算法参数如表 1所示.PMOPSO算法采用了较小的种群规模.根切点位置与齿面渐开线的起始位置较为接近, 选取计算初值时可以认为根切点与齿面渐开线的起始位置重合.在本文的计算方法基础上开发了动力刀塔传动齿轮的辅助设计程序[9].
图 6a是粒子种群中最佳个体在解空间中的运动轨迹.图 6b是粒子种群中最佳个体所对应的目标函数值的演化曲线.在搜索初期, 种群最佳个体与理论最优解的距离较远, 种群中个体的目标函数值较大.随着搜索的推进, 最佳个体的运动轨迹迅速逼近理论最优解, 种群中个体的目标函数值迅速下降, 种群整体性能逐步提升.当时间t>20时, 种群中最佳个体对应的目标函数值已逐渐收敛于零, 种群演化接近尾声时, 种群性能也已逐渐稳定.
图 7是m=3.0、变位系数不同时, 齿面渐开线轮廓与根切点位置的演化情况.当χ=-0.5时, 齿根过渡曲线已深入齿廓内部, 齿面渐开线长度缩短, 齿根过渡曲线逐渐伸长.相比于χ=0的情形, 根切点的位置略有升高, 根切量增大.同时, 齿廓顶部的厚度增加.这里可以根据设计要求灵活地调整根切点位置的计算精度.
如图 8所示, 分别采用叶松林方法[2]、会田俊夫公式[3]和本文提出的方法计算齿轮的根切量, 得到根切量hu随变位系数χ的变化曲线.图 8中, χ=0附近, 齿轮的根切量较小, 不同方法的根切量计算结果一致性最佳; 在χ=-0.5附近, 根切量较大, 不同方法的计算值存在微小差异, 其中本文方法得到的值略高于其他方法.会田俊夫公式是一个近似计算公式, 当变位量增大时, 根切量的计算精度有所下降.相比之下叶松林方法和本文方法的计算结果一致性较好, 稳定性较高.
图 9是在模数m和变位系数χ取不同的值时, 采用本文提出的方法得到的根切量变化曲线.模数m一定时, 根切量hu与变位系数χ的变化关系接近于线性.随着变位系数χ的增大, 根切量hu逐渐减小.变位系数χ一定时, 随着模数m的增大, 根切量hu逐渐增大.变位系数χ值较小时, 模数m的变化对根切量hu的影响较大.
图 10是动力刀塔传动齿轮辅助设计程序的图形界面.用户通过图形界面输入设计参数, 本程序输出相应的齿轮几何模型.程序支持含有根切和几何偏心的齿轮模型.本文中提出的计算方法应用于该程序中的根切点位置计算模块.
本文给出了一种基于PMOPSO技术的齿轮根切量的计算方法.首先, 以渐开线齿轮齿廓参数方程为基础, 建立根切点的计算模型并确定齿轮根切量.然后, 以动力刀塔中关键齿轮为研究对象并计算其根切量, 同时分析变位系数和模数对根切量的影响规律.最后由分析结果可知, 相比于已有的计算方法, 本文提出的方法在应用灵活性和计算精度上具有显著优势.本文方法可以准确分析关键齿轮的根切情况, 这有助于提升动力刀塔中齿轮副的传动性能.
[1] | ChenQi, ZhaoHan, HuangKang. Accurate calculation of dedendum transition curve of profile shifted gear and its fast realization of engineering drawing in CATIA[J]. Journal of Engineering Graphics, 2010, 31(5): 163–168. |
[2] |
叶松林.
渐开线齿廓根切点位置的精确计算[J]. 仪表技术与传感器, 1988(3): 22–24.
( Ye Song-lin. Accurate calculation of undercutting point location of involute tooth profile[J]. Instrument Technology and Sensor, 1988(3): 22–24. ) |
[3] |
会田俊夫. 齿轮的设计和制造: 第一卷圆柱齿轮的设计[M]. 金公望, 译. 北京: 中国农业机械出版社, 1983.
( Toshio A.Gear design and manufacturing:vol.1 design of cylindrical gear[M].Translated by Jin Gong-wang.Beijing:China Agricultural Machinery Press, 1983. ) |
[4] | Kennedy J.The particle swarm:social adaptation of knowledge[C]//Proceedings of IEEE International Conference on Evolutionary Computation.Indianapolis:IEEE, 1997:303-308. |
[5] | Shi Y, Eberhart R.A modified particle swarm optimizer[C]//Proceedings of IEEE International Conference on Evolutionary Computation.IEEE World Congress on Computational Intelligence.Washington D C:IEEE, 1998:69-73. |
[6] | Deb K, Beyer H G. Self-adaptive genetic algorithms with simulated binary crossover[J]. Evolutionary Computation, 2001, 9(2): 197–221. DOI:10.1162/106365601750190406 |
[7] | Deb K, Agrawal S. A niched-penalty approach for constraint handling in genetic algorithms[M]. Vienna: Springer-Vienna, 1999: 235-243. |
[8] |
梁松, 李海波, 张义民.
基于C和gnuplot的渐开线齿轮辅助几何设计程序开发[J]. 工程设计学报, 2014, 21(1): 80–86.
( Liang Song, Li Hai-bo, Zhang Yi-min. Development of auxiliary geometric design software for involute gear based on the C and gnuplot[J]. Chinese Journal of Engineering Design, 2014, 21(1): 80–86. ) |
[9] |
东北大学. 动力伺服刀塔直齿圆柱齿轮辅助参数化设计软件1. 0: 2013SR127443[CP]. 2013.
( Northeastern University.Aided parametric design software for spur cylindrical gear in power tool turrets 1.0:2013SR127433[CP].2013. ) |