基于果蝇优化算法的摆线轮工作齿廓计算方法
【专利摘要】基于果蝇优化算法的摆线轮工作齿廓计算方法，该方法包括如下步骤：根据摆线轮的三种修形方式，分别确定等距加移距以及转角修行的数学方程；阐述果蝇优化算法的基本思路和算法；基于果蝇优化算法，结合摆线针轮修形参数的特定，确定出优化算法中果蝇的初始位置、种群大小和迭代步长等参数，不同的参数设定将会直接影响到算法搜寻目标的能力，分别确定两条不同的修形曲线分别在0到θ0之间的交点θc和θ0到π之间的交点θd，进而确定出修形后摆线轮的工作齿廓为θ＝θd?θc；结合摆线轮修形参数的特点，多次比较果蝇优化算法中果蝇的初始位置、种群数量以及迭代步长等具体参数，最终确定本优化算法中所选定的优化计算参数，极大的提高了计算效率。
【专利说明】
基于果蜗优化算法的摆线轮工作齿廓计算方法
技术领域
[0001] 本发明设及摆线轮的设计与制造技术领域，特别是设及一种果蛹优化算法的摆线 轮工作齿廓计算方法。
【背景技术】
[0002] 摆线针轮晒合副是摆线液压马达、摆线全液压转向器的核屯、元件。由于其多方面 的显著优点而获得广泛的应用。标准齿形的摆线轮齿廓和针齿齿廓之间是没有间隙的，但 在摆线针轮晒合副中，为了补偿制造误差，保证合理的侧隙W利于润滑和保证有效的装拆 方便，实际上摆线轮必须进行修形。
[0003] 生产实践中，采用合适的正等距加负移距组合的方法可使磨削出的齿廓与用转角 修形得到齿廓相接近，等距加移距修形齿廓一方面应保证摆线针齿晒合副齿顶和齿根处有 合适间隙存在，同时又能保证摆线轮齿廓曲线在晒合区间内尽可能地近似于共扼齿廓，达 到多对齿同时晒合。采用运样的修形方式磨削摆线轮时要比转角加等距或转角加移距修形 方便得多。在摆线轮进一步受力分析的过程中，往往需要确认修行后摆线轮的工作齿廓和 非工作齿廓。
[0004] 文献[1]提出通过采用迭代法来得到工作齿廓的起始和终止的晒合点，但是并没 有具体的推导迭代的过程和思路。文献[2]提出通过设定同时晒合的齿数W及摆线齿廓与 针轮齿廓进入和退出晒合是的压力角相等来计算工作齿廓，也并没有给出具体的算法。
[0005] 目前的群智能优化算法主要有遗传算法、蚁群算法、鱼群算法、免疫算法、粒子群 算法，虽然运些算法各有优点，但它们共同的不足是实时性不好，果蛹算法也是一种群智能 算法，学者对其算法寻优特性与其他群智能算法进行了比较，发现果蛹算法具有算法简单、 参数少、易调节、计算量小、寻优精度高的特点，因此本文选用果蛹算法实现对修形后的摆 线轮工作齿廓范围的确定。
[0006] 上述文献[1]和[2]分别为：
[0007] [1]赵岩，关天民.摆线针摆传动转角组合修形下的受力分析新方法[J].机械设计 与制造，1998,(3):45-46.
[000引[2]陈振宇.RV减速器的误差建模与摆线齿廓修形[D].天津大学，2013.
【发明内容】
[0009] 本发明目的是：为了能够准确快速的确定摆线针轮修形之后工作齿廓的范围，通 过改进的果蛹优化算法迅速确定等距加移距修形曲线和转角修行曲线在制定区域的交点。 从而为进一步优化修形曲线提供理论依据。因此本发明提供一种计算周期短、执行效率高 的基于果蛹优化算法的摆线轮工作齿廓计算方法。
[0010] 本发明所采取的技术方案是：
[0011] -种果蛹优化算法的摆线轮工作齿廓计算方法，包括如下步骤：
[0012] 步骤一、根据摆线轮的Ξ种修形方式，分别确定等距加移距W及转角修行的数学 方程；
[0013] 步骤二、阐述果蛹优化算法的基本思路和算法；
[0014] 步骤Ξ、基于果蛹优化算法，结合摆线针轮修形参数的特定，确定出优化算法中果 蛹的初始位置、种群大小和迭代步长等参数，不同的参数设定将会直接影响到算法捜寻目 标的能力，分别确定两条不同的修形曲线分别在0到θ〇之间的交点Θ。和θ〇至I片之间的交点0d， 进而确定出修形后摆线轮的工作齿廓为Θ = Θ<1-Θ。；
[0015] 所述步骤一具体为：
[0016] 摆线轮未修形齿廓方程Fo为：
[0019]摆线轮的加入转角修形A δ之后的数学方程F功：
[0020]
[0021 ]其中Si = (l+ki2-化 1COS目）-〇'5，目 e (0, JI)，ki = azp/rp
[0022] 摆线轮引入等距修形量Δ rp和移距修形量Δ rrp之后的数学方程F2为：
[0023]
[00巧]所述步骤二具体为：
[0026] 果蛹优化算法又称为F0A算法，可W归纳为W下几个重要步骤：
[0027] 1)随机初始果蛹群体位置Χ，Υ。
[0028] 2)分配果蛹个体利用嗅觉捜寻食物随机方向与距离，rand()表示随机生成一个0 到1之间的数：
[0029] Xi = X+rand()，Yi = Y甘and()
[0030] 3)估计果蛹与原点的距离D，再计算味道浓度判定值S，S为距离D的倒数：
[003。 D=(Xi2 巧 i2)0'5，S=l/D
[0032] 4)将味道浓度判定值带入味道浓度判定函数W求出该果蛹个体位置味道浓度 smell:
[0033] 8111611=化]1。1:;[0]1(5)
[0034] 5)找出果蛹群体中味道浓度最大的果蛹：
[0035] [bestsmell bestindex]=max(smell)
[0036] 6)保留最佳味道浓度值与X、Y坐标，此时果蛹群体利用视觉向该位置飞去：
[0037] smellbest = bestsmell
[0038] X=X(bestindex)
[0039] Y = Y(bestindex)
[0040] 7)进入迭代寻优，重复步骤2)~5)，并判断味道浓度是否优于前一迭代味道浓度， 若是则进行步骤6)。
[0041] 所述步骤Ξ具体为：
[0042] 图1所示的0、1、2Ξ条曲线分别是方程F〇FiF2曲线在区间Θ e (0，31)上对应的齿廓方 程，现在要确定修形后的工作齿廓就是确定θ = θ<ι-θ。。由于确定交点C和交点D算法思路相 同，此算法W点C的计算为例说明基于果蛹优化算法求解的基本思路：
[0043] 1)在随机给定果蛹的初始位置时，将θ = 〇代入方程Fi中得到果蛹初始位置（Χ0， yo);在确定果蛹的数量时，经过比较3，10，20，30只果蛹的计算效率和能力，筛选出最优果 蛹数目20只作为果蛹的数量。
[0044] 2)在确定随机方向和距离时结合修形曲线的区间(0，θ〇)，得：
[0045]
[0046] 3)在估算果蛹的位置距离时方法不变：
[0047] D=(Xi2 巧 i2)〇'5，S=l/D
[0048] 4)将味道浓度判定值带入味道浓度判定函数W求出该果蛹个体位置味道浓度 smell:
[0049] 结合摆线轮两个修形曲线方程，此处的浓度判定函数为：
[0050] 8111611=化]1(31:;[0]1(5) = |Fi(S)-F2(S)
[0051] 5)找出果蛹群体中味道浓度最大的果蛹：
[0052] [bestsmell bestindex]=max(smell)
[0053] 6)保留最佳味道浓度值与Χ、Υ坐标，此时果蛹群体利用视觉向该位置飞去：
[0054] smellbest = bestsmell
[0055] X=X(bestindex)
[0056] Y = Y(bestindex)
[0057] 7)进入迭代寻优，重复步骤2)~5)，并判断味道浓度是否优于前一迭代味道浓度， 若是则进行步骤6)。
[0058] 本发明具有的优点和积极效果是：
[0059] 本发明通过分析摆线轮修形的数学方程和果蛹优化算法提出了一种基于果蛹优 化算法的摆线轮工作齿廓计算方法;结合摆线轮修形参数的特点，多次比较果蛹优化算法 中果蛹的初始位置、种群数量W及迭代步长等具体参数，最终确定本优化算法中所选定的 优化计算参数，极大的提高了计算效率，迅速而准确的确定了工作齿廓的范围。本发明采用 Matlab中的编写果蛹优化算法，取得良好的效果。
【附图说明】
[0060] 图1 F〇FiF2曲线在区间0e(〇,3i)上对应的齿廓曲线。
【具体实施方式】
[0061] 为能进一步了解本发明的
【发明内容】
、特点及功效，兹例举W下实施例，并配合附图 详细说明如下：
[0062] 请参阅图1，一种基于果蛹优化算法的摆线轮工作齿廓计算方法，包括如下步骤：
[0063] 步骤一根据摆线轮的Ξ种修形方式，分别确定等距加移距W及转角修行的数学方 程：
[0064] 摆线轮未修形齿廓方程Fo为：
[00 化]
[0066] 其中Si = (l+ki2-化 1COS目）-〇'5，目 e (0, JI)，ki = azp/rp
[0067] 摆线轮的加入转角修形Δ δ之后的数学方程Fi为：
[006引
[0069] 其中 Si=(l+ki2-2kicos 目厂 °'5,目 e(〇，ji)，ki = azp/rp
[0070] 摆线轮引入等距修形量Δ。和移距修形量Δ rrp之后的数学方程F2为：
[0071]
[0072] 其中S2 = (l+k22-2k2C〇s目）-〇'5,目 e (0, JI)，k2 = azp/(rp- Δ
[0073] 步骤二阐述果蛹优化算法的基本思路和算法：
[0074] 果蛹优化算法又称为F0A算法，可W归纳为W下几个重要步骤：
[0075] 1)随机初始果蛹群体位置Χ，Υ。
[0076] 2)分配果蛹个体利用嗅觉捜寻食物随机方向与距离，rand ()表示随机生成一个0 到1之间的数：
[0077] Xi = X+rand〇，Yi = Yhand〇
[0078] 3)估计果蛹与原点的距离D，再计算味道浓度判定值S，S为距离D的倒数：
[0079] D=(Xi2 巧 i2)〇'5，S=l/D
[0080] 4)将味道浓度判定值带入味道浓度判定函数W求出该果蛹个体位置味道浓度 smell:
[0081 ] smell 二 F'unction(S)
[0082] 5)找出果蛹群体中味道浓度最大的果蛹：
[0083] [bestsmell bestindex]=max(smell)
[0084] 6)保留最佳味道浓度值与Χ、Υ坐标，此时果蛹群体利用视觉向该位置飞去：
[00 化]smenbest = bestsmell
[00 化]X=X(bestindex)
[0087] Y = Y(bestindex)
[0088] 7)进入迭代寻优，重复步骤2)~5)，并判断味道浓度是否优于前一迭代味道浓度， 若是则进行步骤6)。
[0089] 步骤Ξ、基于果蛹优化算法，结合摆线针轮修形参数的特定，确定出优化算法中果 蛹的初始位置、种群大小和迭代步长等参数，不同的参数设定将会直接影响到算法捜寻目 标的能力。确定两条不同的修形曲线分别在0到θ〇之间的交点Θ。和θ〇至I片之间的交点0d，进而 确定出修形后摆线轮的工作齿廓为θ = θ<ι-θ。：
[0090] 图1所示的0，1，2Ξ条曲线分别是方程F〇FiF2曲线在区间0e(〇,3i)上对应的齿廓方 程，现在要确定修形后的工作齿廓就是确定θ = θ<ι-θ。。由于确定交点C和交点D算法思路相 同，此次算法W点C的计算为例说明基于果蛹优化算法求解的基本思路：
[0091] 1)在随机给定果蛹的初始位置时，将θ = 0代入方程Fi中得到果蛹初始位置（Χ0， yo);在确定果蛹的数量时，经过比较3，10，20，30只果蛹的计算效率和能力，筛选出最优果 蛹数目20只作为果蛹的数量。
[0092] 2)在确定随机方向和距离时结合修形曲线的区间(0，θ〇)，得：
[0093]
[0094] 3)在估算果蛹的位置距离时方法不变：
[00 巧]D=(Xi2 巧 i2)〇'5，S=l/D
[0096] 4)将味道浓度判定值带入味道浓度判定函数W求出该果蛹个体位置味道浓度 smell:
[0097] 结合摆线轮两个修形曲线方程，此处的浓度判定函数为：
[009引 8111611=化]1。1:;[0]1(5) = |Fi(S)-F2(S)
[0099] 5)找出果蛹群体中味道浓度最大的果蛹：
[0100] [bestsmell bestindex]=max(smell)
[0101] 6)保留最佳味道浓度值与Χ、Υ坐标，此时果蛹群体利用视觉向该位置飞去：
[0102] smellbest = bestsmell
[0103] X=X(bestindex)
[0104] Y = Y(bestindex)
[0105] 7)进入迭代寻优，重复步骤2)~5)，并判断味道浓度是否优于前一迭代味道浓度，
[0106] 若是则进行步骤6)。
[0107] 本发明具有的优点和积极效果是：
[0108] 本发明通过分析摆线轮修形的数学方程和果蛹优化算法提出了一种基于果蛹优 化算法的摆线轮工作齿廓计算方法;结合摆线轮修形参数的特点，多次比较果蛹优化算法 中果蛹的初始位置、种群数量W及迭代步长等具体参数，最终确定本优化算法中所选定的 优化计算参数，极大的提高了计算效率，迅速而准确的确定了工作齿廓的范围。本发明采用 Matlab中的编写果蛹优化算法，取得良好的效果。
【主权项】
1.基于果蛹优化算法的摆线轮工作齿廓计算方法，其特征在于:该计算方法包括如下 步骤， 步骤一、根据摆线轮的Ξ种修形方式，分别确定等距加移距W及转角修行的数学方程； 步骤二、阐述果蛹优化算法的基本思路和算法； 步骤Ξ、基于果蛹优化算法，结合摆线针轮修形参数的特定，确定出优化算法中果蛹的 初始位置、种群大小和迭代步长等参数，不同的参数设定将会直接影响到算法捜寻目标的 能力，分别确定两条不同的修形曲线分别在0到θ〇之间的交点Θ。和θ〇至I片之间的交点0d，进而 确定出修形后摆线轮的工作齿廓为θ = θ<ι-θ。； 所述步骤一具体为： 摆线轮未修形齿廓方程Fo为：摆线轮的加入转角修形A δ之后的数学方程Fi为：摆线轮引入等距修形量A。和移距修形量Δ rrp之后的数学方程F2为：所述步骤二具体为： 果蛹优化算法又称为F0A算法，可W归纳为W下几个重要步骤： 1) 随机初始果蛹群体位置Χ，Υ; 2) 分配果蛹个体利用嗅觉捜寻食物随机方向与距离，randO表示随机生成一个0到1之 间的数： Xi = X甘and ()，Yi = Y+rand () 3) 估计果蛹与原点的距离D，再计算味道浓度判定值S，S为距离D的倒数： D=(Xi2 巧 i2)〇'5，s=l/D 4) 将味道浓度判定值带入味道浓度判定函数W求出该果蛹个体位置味道浓度smell: smell =Funct ion (S) 5) 找出果蛹群体中味道浓度最大的果蛹： [bestsmell bestindex]=max(smell) 6) 保留最佳味道浓度值与X、Y坐标，此时果蛹群体利用视觉向该位置飞去： smellbest = bestsmell X = X(bestindex) Y = Y(bestindex) 7) 进入迭代寻优，重复步骤2)~5)，并判断味道浓度是否优于前一迭代味道浓度，若是 则进行步骤6); 所述步骤Ξ具体为： 设Ξ条曲线分别是方程F〇FiF2曲线在区间目e (0,31)上对应的齿廓方程，现在要确定修 形后的工作齿廓就是确定Θ = θ4-θ。；由于确定交点C和交点D算法思路相同，此算法W点C的 计算为例说明基于果蛹优化算法求解的基本思路： 1) 在随机给定果蛹的初始位置时，将θ = 0代入方程Fi中得到果蛹初始位置(xo，yo);在 确定果蛹的数量时，经过比较3,10,20,30只果蛹的计算效率和能力，筛选出最优果蛹数目 20只作为果蛹的数量； 2) 在确定随机方向和距离时结合修形曲线的区间(0，θ〇)，得：3) 在估算果蛹的位置距离时方法不变： D=(Xi2 巧 i2)〇'5，s=l/D 4) 将味道浓度判定值带入味道浓度判定函数W求出该果蛹个体位置味道浓度smell: 结合摆线轮两个修形曲线方程，此处的浓度判定函数为： 8111611=化]1(31：；[〇]1(5)= |Fi(S)-F2(S) 5) 找出果蛹群体中味道浓度最大的果蛹： [bestsmell bestindex]=max(smell) 6) 保留最佳味道浓度值与Χ、Υ坐标，此时果蛹群体利用视觉向该位置飞去： smellbest = best smell X = X(bestindex) Y = Y(bestindex) 7) 进入迭代寻优，重复步骤2)~5)，并判断味道浓度是否优于前一迭代味道浓度，若是 则进行步骤6)。
【文档编号】F16H55/08GK105972185SQ201610539273
【公开日】2016年9月28日
【申请日】2016年7月10日
【发明人】刘志峰, 张涛, 王冰
【申请人】北京工业大学