本发明提供一种基于自然齿面活动标架的点啮合齿面的齿轮设计方法，所述设计方法包括步骤，S10：齿面沿接触点迹线的局部结构的设计；S20：基于自然活动标架的齿面接触点迹线以外的齿面整体设计。本发明提供的设计方法引入曲面自然活动标架作为齿面设计的坐标系，将点啮合齿面的设计问题，归结为以第一齿面上的包络特征曲线作为样条曲线的曲面拟合逼近问题，通过第一齿面的运动不变量描述和设计点啮合齿面的参数，使点啮合齿面的整体设计摆脱对于特定机床结构的依赖，因而具有通用性和精确性。
1.一种基于自然齿面活动标架的点啮合齿面的齿轮设计方法，所述齿轮包括互相点啮合的第一齿轮和第二齿轮，所述第一齿轮上具有第一齿面，所述第二齿轮上具有第二齿面，所述设计方法包括以下步骤：S10：齿面沿接触点迹线的局部结构的设计
S11：给定第一齿面∑1：r(1)＝r(1)(u1，v1)并在所述第一齿面上指定一条接触点迹线L1：由齿面共轭原理设计第二齿面∑2上与接触点迹线L1共轭的接触点迹线L2：及所述接触点迹线L2的单位法线矢量参数方程 使所述第一
齿面∑1与所述第二齿面∑2沿着接触点迹线L1啮合传动，并满足式(1)给定的相对运动A式中 为使啮合的起始位置与第二齿轮转动的起始位置重合，即 时，第一齿面∑1相对固定坐标系S1[O1；x1y1z1]转过的角度；ci(i＝2，…，n)为多项式系数，按预定的两齿轮的相对运动规律给出； 为第二齿轮与第一齿轮的齿数比±Z2/Z1； 是第一齿面∑1自转的角位移； 是第二齿面∑2自转的角位移；
S12：给定啮合点的位置对安装误差的敏感度ε和齿面接触许用应力值[σ]，再设计第二齿面∑2沿接触点迹线L2的二阶微分结构参数：接触点迹线L2上各点沿其切线方向的法曲率及其短程挠率 接触点迹线L2的切线在第二齿面∑2切平面上的垂直方向的法曲率S20：基于自然活动标架的齿面接触点迹线以外的齿面整体设计
S21：第一齿面∑1与第二齿面∑2的相对运动速度的曲面运动不变量参数描述根据第一齿面∑1：r(1)＝r(1)(u1，υ1)，则参数u1，υ1构成正交参数网线，在第一齿面∑1上任意一接触点M(u1M，υ1M)取定一个活动标架场当接触点M沿正交参数网线u1，υ1以速度 移动时，活动标架Sft也随之移动并绕接触点M转动，由微分几何的曲面论的基本定理，活动标架Sft的角速度运动方程如式(2)ω1＝ω23e1-ω13e2+ω12n1M (2)
设接触点迹线L2是第一齿面∑1在相对运动下的共轭齿面∑p：r(p)＝r(p)(up，υp)上的一条曲线，其方程为 sp为接触点迹线L2的弧长参数；接触点迹线L2上任意一接触点M的单位切向量为α，曲面在M点的单位法向量为npM，在接触点迹线L2上接触点M处取定一个正交的活动标架Sfp[M； α ν＝npM×α npM]；
当接触点M沿接触点迹线L2以速度 移动时，活动标架Sfp[M； α ν＝npM×α npM]也随之移动并绕接触点M转动，由微分几何的曲面论的基本定理，活动标架Sfp[M； α ν＝npM×α npM]的角速度运动方程如式(3)式中， 分别是第二齿面∑2在接触点M沿α方向的短程挠率和法曲率， 是第二齿面∑2在接触点M沿α方向的短程曲率；接触点迹线L2
与固定坐标系S(1)[O(1)； x(1) y(1) z(1)]固结的第一齿面∑1通过相对于与固定坐标系S(2)[O(2) ；x(2) y(2) z(2)]固结的共轭齿面∑p的相对运动Ψ包络出第二齿面∑2；设定Lg是形成包络共轭齿面∑p的特征线，则在第一齿面∑1沿接触点迹线L2相对共轭齿面∑p按相对运动Ψ的任意时刻，第一齿面∑1与第二齿面∑2相切于特征线Lg，特征线Lg和接触点迹线L2交于接触点M，而在接触点M，对于第一齿面∑1与第二齿面∑2分别有活动标架Sft和活动标架Sfp，共轭齿面∑p与第二齿面∑2同固结于固定坐标系S(2)[O(2)； x(2) y(2) z(2)]，共轭齿面∑p是第二齿面∑2的拟合逼近目标，由式(2)和式(3)知第一齿面∑1相对第二齿面∑2的角速度为式(4)考虑S(1)、Sft、Sfp和S(2)之间的相对运动，则在接触点M处第一齿面∑1相对第二齿面∑2的速度为式(5)按式(5)的方式，将式(4)中的 因子提出来，并将 因子指定为基准速率，不失一般性，令 而由于接触点迹线L2的方程已经由与接触点迹线L1相对运动的共轭条件确定，故 也是已经确定的函数量，则在第一齿面∑1上任意点P(u1，υ1)处，第一齿面∑1相对第二齿面∑2的速度函数表示为式(6)式(6)中的r1(ft)表示点P(ut，υt)在活动标架Sft下的位置矢量，Δ是活动标架Sft与活动标架Sfp之间的夹角；
S22：从第一齿面∑1到第二齿面∑2的运动变换及第二齿面∑2方程的曲面运动不变几何量描述从第一齿面∑1出发，分别由相对运动Ψ展成第二齿面∑2和由相对运动A展成共轭齿面∑p的坐标变换过程以及第二齿面∑2的设计思路；其中相对运动Ψ是在第一齿面∑1和共轭齿面∑p的活动标架下由第一齿面∑1和共轭齿面∑p的运动不变量描述；相对运动A是在点啮合齿轮副固定坐标系下由齿轮副结构参数描述的；沿相对运动Ψ和相对运动A，从固定坐标系S(1)到固定坐标系S(2)的变换矩阵为式(7)式中
上任意接触点M在固定坐标系S(1)[O(1)； x(1) y(1) z(1)]中的位置矢量函数，此接触点M与共轭齿 面 ∑ p上 的 接触 点 迹 线L 2 上的 接 触 点M 关 于 相对 运 动Ψ 共 轭 ；
表示接触点迹线L2在固定坐标系S2[O(2)； x(2)
第一齿面∑1在相对运动Ψ下的第二齿面∑2在固定坐标系S2[O(2)； x(2) y(2) z(2)]下的方程为式(8)式中， 表示第一齿面∑1在固定坐标系S(1)
[O(1)； x(1) y(1) z(1)]下的矢量参数方程，N1(u1，υ1)是第一齿面∑1上任意一点P处单位法矢的矢量参数方程，此点P是第一齿面∑1上不同于接触点M的另一任意点；
设第一齿面∑1和第二齿面∑2间的，与最小相对曲率对应的相对主方向为αm，沿αm的第一齿面∑1、第二齿面∑2的法曲率为 短程挠率为 则沿相对主方向αm的第一齿面∑1、第二齿面∑2的相对短程挠率为零，即满足式(9)由欧拉-贝特朗公式，对于共轭齿面∑p沿相对主方向αm满足式(10)、(11)其中，φ2和φ3满足式(12)，
在式(10)至式(13)中， 分别是第一齿面∑1在接触点M沿ei(i＝1，2)方向的短程挠率和法曲率； 和 分别是共轭齿面∑p在接触点M沿接触点迹线L2的切线方向即单位切矢α方向的短程挠率和法曲率及α方向在公切平面上的垂直方向ν的法曲率；
和 分别是共轭齿面∑p在接触点M沿接触点迹线L1的切线方向即单位切矢α1方向的短程挠率和法曲率及α1方向在公切平面上的垂直方向v1的法曲率；由于在步骤S10中已经确定了接触点迹线L1和接触点迹线L2，那么接触点迹线L1的单位切矢α1和接触点迹线L2的单位切矢α就已经确定，故φ2-φ3是已经确定的已知量；
基于最终形成的第二齿面∑2需满足步骤S10的设计要求，故对于第二齿面∑2，沿相对主方向αm满足式(14)活动标架Sft和活动标架Sfp之间的夹角Δ满足式(15)
Δ(u1M，υ1M，sp)＝φ2(sp)-φ1(u1M，υ1M) (15)对式(9)～(14)联立求解，用接触点M处第一齿面∑1和第二齿面∑2的各个方向的法曲率和短程挠率表示φ1、φ2和φ3，具体表示为φ1(utM，υtM)，φ2(sp)和φ3(s1)，从而确定相对主方向αm，将φ1(utM，υtM)，φ2(sp)代入式(15)，即可将夹角Δ表达为运动不变量utM，υtM，sp的函数；
将上述方法中求得的Δ分别代入式(6)和(7)，分别求得在任意点P(u1，υ1)处第一齿面∑1相对第二齿面∑2的速度函数和从第一齿面∑1到第二齿面∑2的变换矩阵函数，再将式(6)和式(7)代入式(8)并消去参数u1，就得到第一齿面∑1在相对运动Ψ下的第二齿面∑2在与共轭齿面∑p固结的固定坐标系S(2)[O(2)； x(2) y(2) z(2)]下的矢量函数方程(16)S24：第二齿面∑2在接触点迹线L2以外的区域与第一齿面∑1的间隙在共轭齿面∑p上任意一点Qg点处第二齿面∑2对共轭齿面∑p的间隙满足式(17)式中 是共轭齿面∑p上任意一点Qg处单位法矢的矢量参数方程，如果选定了点Qp，即确定了矢量 则式(17)所表达的三个标量方程正好可以求解参数u1、 和δk，从而确定第二齿面∑2与共轭齿面∑p在任意位置的间隙δk以及点Qg的位置矢量S25：第二齿面∑2整体优化设计模型
在过任意接触点M的第一齿面∑1在相对运动A下的特征线Lg上，分别位于接触点M的两侧的取两点 和 其对应的曲坐标分别设为(sp，υ11)、(sp，υ12)，特征线Lg上点 和 之间的各点到第二齿面∑2的距离δk与该点到接触点M的距离成正比，设曲线Lgp、 点 和分别是特征线Lg、点 和 在第二齿面∑2上的投影，当第一齿面∑1对第二齿面∑2的间隙小于一个给定值时，第一齿面∑1上点 和 之间的狭长齿面与第二齿面∑2上 和 点之间的狭长曲面在承受满载的情况下将贴合在一起，并在第二齿面∑2上形成一个瞬时接触区，从啮合的起始点M1到啮合的终点M2的整个啮合过程将在第二齿面∑2上连续形成一系列这样的瞬时接触区，并拼接成齿面接触区，在齿面承受满载的情况下，控制齿面接触区的面积大小占整个第二齿面面积满足式(18)式中 是用参数(sp，υ1)表示的共轭齿面∑p方程，并由
式中δ0是啮合的第一齿面∑1和第二齿面∑2在承受满载的情况下的齿面接触区的边界各点在齿面无载接触时的间隙；
求解由式(18)和(19)组成的不等式方程式组，得到两个函数闭区间u1M(sp)∈[u1M1(sp) u1M2(sp)]和υ1M(sp)∈[υ1M1(sp) υ1M2(sp)]，将其代入式(16)即可求得待求的第二齿面∑2的方程的设计区间
基于自然齿面活动标架的点啮合齿面设计方法
技术领域
[0001] 本发明涉及齿轮技术领域，尤其涉及到相互啮合齿轮的啮合齿面，具体为一种基于自然齿面活动标架的点啮合齿面的设计方法。
背景技术
[0002] 点啮合齿面是为了满足高精高速重载的机械动力传递要求，通过修形将共轭齿面失配为点接触齿面而得到的，它的设计必须考虑齿轮承载弹性变形、制造和安装误差对齿面啮合特性的影响，它的类型包括圆柱鼓形齿面齿轮、蜗轮蜗杆和螺旋锥齿轮等等。当前点啮合齿面的主动设计仅止于上述沿齿面接触点迹线的局部设计，而接触点迹线以外的大范围齿面结构则通过齿面接触分析(Tooth Contact Analysis，TCA)来确定；但是，这种方法通过所谓“齿面修形”将点啮合齿面的描述与特定机床上齿面的加工调整参数联系在一起，以模拟“试切”的方式最终确定齿面的整体结构。这种齿面设计方法以满足特定摇台式齿轮加工机床的加工可行性为优先目标，因而，不具通用性，更无法据此对齿面啮合特性进行全局优化设计，造成这种情况的原因在于：当前的点啮合齿面的设计是在特定机床坐标系下以特定机床的加工调整参数描述和设计齿面。
发明内容
[0003] 本发明旨在解决现有技术中存在的技术问题。
[0004] 本发明通过引入第一齿面上的自然活动标架作为描述和设计点啮合齿面的坐标系，将点啮合齿面的设计问题，归结为以第一齿面上的包络特征曲线作为样条曲线的曲面拟合逼近问题，通过第一齿面的运动不变量描述和设计点啮合齿面的参数，使点啮合齿面的整体设计摆脱对于特定机床结构的依赖，因而具有通用性和精确性，并能够据此方法对齿面啮合特性进行整体优化设计。
[0005] 本发明提出的一种基于自然齿面活动标架的点啮合齿面的齿轮设计方法，所述齿轮包括互相点啮合的第一齿轮和第二齿轮，所述第一齿轮上具有第一齿面，所述第二齿轮上具有第二齿面，所述设计方法包括以下步骤：
[0006] S10：齿面沿接触点迹线的局部结构的设计
[0007] S11：给定第一齿面∑1：r(1)＝r(1)(u1，v1)并在所述第一齿面上指定一条接触点迹线L1： 由齿面共轭原理设计第二齿面∑2上与接触点迹线L1共轭的接触点迹线L2： 及所述接触点迹线L2的单位法线矢量参数方程 使所述
第一齿面∑1与所述第二齿面∑2沿着接触点迹线L1啮合传动，并满足式(1)给定的相对运动A
[0008]
[0009] 式中 为使啮合的起始位置与第二齿轮转动的起始位置重合，即 时，第一齿面∑1相对固定坐标系S1[O1；x1y1z1]转过的角度；ci(i＝2，…，n)为多项式系数，按预定的两齿轮的相对运动规律给出；为第二齿轮与第一齿轮的齿数比±Z2/Z1；是第一齿面∑1自转的角位移； 是第二齿面∑2自转的角位移；
[0010] S12：给定啮合点的位置对安装误差的敏感度ε和齿面接触许用应力值[σ]，再设计第二齿面∑2沿接触点迹线L2的二阶微分结构参数：接触点迹线L2上各点沿其切线方向的法曲率 及其短程挠率 接触点迹线L2的切线在第二齿面∑2切平面上的垂直方向的法曲率
[0011] S20：基于自然活动标架的齿面接触点迹线以外的齿面整体设计
[0012] S21：第一齿面∑1与第二齿面∑2的相对运动速度的曲面运动不变量参数描述[0013] 根据第一齿面∑1：r(1)＝r(1)(u1，υ1)，则参数u1，υ1构成正交参数网线，在第一齿面∑1上任意一接触点M(u1M，υ1M)取定一个活动标架场
[0014] 当接触点M沿正交参数网线u1，υ1以速度 移动时，活动标架Sft也随之移动并绕接触点M转动，由微分几何的曲面论的基本定理，活动标架Sft的角速度运动方程如式(2)
[0015] ω1＝ω23e1-ω13e2+ω12n1M (2)
[0016] 式中，
[0017] 设接触点迹线L2是第一齿面∑1在相对运动下的共轭齿面∑p：r(p)＝r(p)(up，υp)上的一条曲线，其方程为 sp为接触点迹线L2的弧长参数；接触点迹线L2上任意一接触点M的单位切向量为α，曲面在M点的单位法向量为npM，在接触点迹线L2上接触点M处取定一个正交的活动标架Sfp[M；α ν＝npM×α npM]；
[0018] 当接触点M沿接触点迹线L2以速度 移动时，活动标架Sfp[M；α ν＝npM×α npM]也随之移动并绕接触点M转动，由微分几何的曲面论的基本定理，活动标架Sfp[M；α ν＝npM×α npM]的角速度运动方程如式(3)
[0019]
[0020] 式中， 分别是第二齿面∑2在接触点M沿α方向的短程挠率和法曲率， 是第二齿面∑2在接触点M沿α方向的短程曲率；
[0021] 与固定坐标系S(1)[O(1)；x(1) y(1) z(1)]固结的第一齿面∑1通过相对于与固定坐标系S(2)[O(2)；x(2) y(2) z(2)]固结的共轭齿面∑p的相对运动Ψ包络出第二齿面∑2；设定Lg是形成包络共轭齿面∑p的特征线，则在第一齿面∑1沿接触点迹线L2相对共轭齿面∑p按相对运动Ψ的任意时刻，第一齿面∑1与第二齿面∑2相切于特征线Lg，特征线Lg和接触点迹线L2交于接触点M，而在接触点M，对于第一齿面∑1与第二齿面∑2分别有活动标架Sft和活动标架Sfp，共轭齿面∑p与第二齿面∑2同固结于坐标系S(2)[O(2)；x(2) y(2) z(2)]，共轭齿面∑p是第二齿面∑2的拟合逼近目标，由式(2)和式(3)知第一齿面∑1相对第二齿面∑2的角速度为式(4)
[0022]
[0023] 考虑S(1)、Sft、Sfp和S(2)之间的相对运动，则在接触点M处第一齿面∑1相对第二齿面∑2的速度为式(5)
[0024]
[0025] 按式(5)的方式，将式(4)中的 因子提出来，并将 因子指定为基准速率，不失一般性，令 而由于接触点迹线L2的方程已经由与接触点迹线L1相对运动的共轭条件确定，故 也是已经确定的函数量，则在第一齿面∑1上任意点P(ut，υt)处，第一齿面∑1相对第二齿面∑2的速度函数表示为式(6)
[0026]
[0027] 式(6)中的r1(ft)表示点P(ut，υt)活动标架Sft下的位置矢量，Δ是活动标架Sft与活动标架Sfp之间的夹角；
[0028] S22：从第一齿面∑1到第二齿面∑2的运动变换及第二齿面∑2方程的曲面运动不变几何量描述
[0029] 从第一齿面∑1出发，分别由相对运动Ψ展成第二齿面∑2和由相对运动A展成共轭齿面∑p的坐标变换过程以及第二齿面∑2的设计思路；其中相对运动Ψ是在第一齿面∑1和共轭齿面∑p的活动标架下由第一齿面∑1和共轭齿面∑p的运动不变量描述；相对运动A是在点啮合齿轮副固定坐标系下由齿轮副结构参数描述的；沿相对运动Ψ和相对运动A，从固定坐标系S(1)到固定坐标系S(2)的变换矩阵为式(7)
[0030]
[0031] 式中
[0032]
[0033]
[0034] 其中， 表示第一齿面∑1上任意接触点M在固定坐标系S(1)[O(1)；x(1) y(1) z(1)]中的位置矢量函数，此接触点M与共轭齿面∑p上的接触点迹线L2上的接触点M关于相对运动Ψ共轭；
[0035] 表示接触点迹线L2在固定坐标系
[0036] S2[O(2)；x(2) y(2) z(2)]中的矢量参数方程；
[0037] 第一齿面∑1在相对运动Ψ下的第二齿面∑2在固定坐标系S2[O(2)；x(2) y(2) z(2)]下的方程为式(8)
[0038]
[0039] 式中， 表示第一齿面∑1在固定坐标系S(1)[O(1)；x(1) y(1) z(1)]下的矢量参数方程，N1(u1，υ1)是第一齿面∑1上任意一点P处单位法矢的矢量参数方程，此点P是第一齿面∑1上不同于接触点M的另一任意点；
[0040] S23：运动参数Δ的确定
[0041] 设第一齿面∑1和第二齿面∑2间的，与最小相对曲率对应的相对主方向为αm，沿αm的第一齿面∑1、第二齿面∑2的法曲率为 短程挠率为 则沿相对主方向αm的第一齿面∑1、第二齿面∑2的相对短程挠率为零，即满足式(9)
[0042]
[0043] 由欧拉-贝特朗公式，对于共轭齿面∑p沿相对主方向αm满足式(10)、(11)[0044]
[0045]
[0046] 其中，φ2和φ3满足式(12)，
[0047] φ2-φ3＝arc cos(α1·α) 12)
[0048] 对于第一齿面∑1沿相对主方向αm有式(13)
[0049]
[0050] 在式(10)至式(13)中， 分别是第一齿面∑1在接触点M沿ei(i＝1，2)方向的短程挠率和法曲率； 和 分别是共轭齿面∑p在接触点M沿接触点迹线L2的切线方向即单位切矢α方向的短程挠率和法曲率及α方向在公切平面上的垂直方向ν的法曲率；
和 分别是共轭齿面∑p在接触点M沿接触点迹线L1的切线方向即单位切矢α1方
向的短程挠率和法曲率及α1方向在公切平面上的垂直方向v1的法曲率；由于在步骤S10中已经确定了接触点迹线L1和接触点迹线L2，那么接触点迹线L1的单位切矢α1和接触点迹线L2的单位切矢α就已经确定，故φ2-φ3是已经确定的已知量；
[0051] 基于最终形成的第二齿面∑2需满足步骤S10的设计要求，故对于第二齿面∑2，沿相对主方向αm满足式(14)
[0052]
[0053] 活动标架Sft与活动标架Sfp之间的夹角Δ满足式(15)
[0054] Δ(u1M，υ1M，sp)＝φ2(sp)-φ1(u1M，υ1M) (15)
[0055] 对式(9)～(14)联立求解，用接触点M处第一齿面∑1和第二齿面∑2的各个方向的法曲率和短程挠率表示φ1、φ2和φ3，具体表示为φ1(utM，υtM)，φ2(sp)和φ3(s1)，从而确定相对主方向αm，将φ1(utM，υtM)，φ2(sp)代入式(15)，即可将夹角Δ表达为运动不变量utM，υtM，sp的函数；
[0056] 将上述方法中求得的Δ分别代入式(6)和(7)，分别求得在任意点P(u1，υ1)处第一齿面∑1相对第二齿面∑2的速度函数和从第一齿面∑1到第二齿面∑2的变换矩阵函数，再将式(6)和式(7)代入式(8)并消去参数u1，就得到第一齿面∑1在相对运动Ψ下的第二齿面∑2在与共轭齿面∑p固结的固定坐标系S(2)[O(2)；x(2) y(2) z(2)]下的矢量函数方程(16)[0057]
[0058] S24：第二齿面∑2在接触点迹线L2以外的区域与第一齿面∑1的间隙
[0059] 在共轭齿面∑p上任意一点Qg点处第二齿面∑2对共轭齿面∑p的间隙满足式(17)[0060]
[0061] 式中 是共轭齿面∑p上任意一点Qg处单位法矢的矢量参数方程，如果选定了点Qp，即确定了矢量 则式(17)所表达的三个标量方程正好可以求解参数u1、和δk，从而确定第二齿面∑2与共轭齿面∑p在任意位置的间隙δk以及点Qg的位置矢量[0062] S25：第二齿面∑2整体优化设计模型
[0063] 在过任意接触点M的第一齿面∑1在相对运动A下的特征线Lg上，分别位于接触点M的两侧的取两点 和 其对应的曲坐标分别设为(sp，υ11)、(sp，υ12)，特征线Lg上点 和之间的各点到第二齿面∑2的距离δk与该点到接触点M的距离成正比，设曲线Lgp、点 和分别是特征线Lg、点 和 在第二齿面∑2上的投影，当第一齿面∑1对第二齿面∑2的间隙 小于一个给定值时，第一齿面∑1上点 和 之间的狭长齿面与第二齿面∑2上 和 点之间的狭长曲面在承受满载的情况下将贴合在一起，并在第二齿面∑2上形成一个瞬时接触区，从啮合的起始点M1到啮合的终点M2的整个啮合过程将在第二齿面∑2上连续形成一系列这样的瞬时接触区，并拼接成齿面接触区，在齿面承受满载的情况下，控制齿面接触区的面积大小占整个第二齿面面积满足式(18)
[0064]
[0065] 式中 是用参数(sp，υ1)表示的共轭齿面∑p方程，并由式(17)求得的 并代入
方程 得到；而式中
[0066] 再按式(19)进行优化设计约束：
[0067]
[0068] 式中δ0是啮合的第一齿面∑1和第二齿面∑2在承受满载的情况下的齿面接触区的边界各点在齿面无载接触时的间隙；
[0069] 求解由式(18)和(19)组成的不等式方程式组，得到两个函数闭区间u1M(sp)∈[u1M1(sp)u1M2(sp)]和υ1M(sp)∈[υ1M1(sp)υ1M2(sp)]，将其代入式(16)即可求得待求的第二齿面∑2的方程的设计区间
[0070] 本发明基于齿面自然活动标架，由第一齿面运动不变量参数描述的点啮合齿面整体设计方法，通过引入曲面自然活动标架作为齿面设计的坐标系，将点啮合齿面的设计问题，转换归结为以第一齿面上的包络特征曲线作为样条曲线的曲面拟合逼近问题，通过第一齿面的运动不变量描述和设计点啮合齿面的参数，使点啮合齿面的整体设计摆脱对于特定机床结构的依赖，因而具有通用性和精确性，并能够据此方法对齿面啮合特性进行全局优化设计。
附图说明
[0071] 本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
[0072] 图1本发明方案中的点啮合齿轮传动运动描述所涉及的坐标系及其相对位置关系图；
[0073] 图2本发明方案中的第一齿面上沿接触点迹线的自然活动标架相对齿面的运动关系示意图；
[0074] 图3是本发明方案中的第二齿面上沿接触点迹线的自然活动标架相对齿面的运动关系示意图；
[0075] 图4是本发明方案中的点啮合齿面沿接触点迹线相互啮合时的相对运动关系示意图；
[0076] 图5是本发明方案中的齿面坐标系运动变换及齿面设计流程图；
[0077] 图6是本发明方案中的第一齿面与第二齿面在啮合接触点M处的公切平面上活动标架的相互关系示意图；
[0078] 图7是本发明方案中的第二齿面∑2沿接触线Lg与共轭齿面∑p的间隙关系示意图；
[0079] 图8是本发明方案中的点啮合齿面接触区的控制关系示意图。
具体实施方式
[0080] 为了能够更清楚地理解本发明的上述目的、特征和优点，下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是，在不冲突的情况下，本申请的实施例及实施例中的特征可以相互组合。
[0081] 在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发明还可以采用其他不同于在此描述的方式来实施，因此，本发明的保护范围并不受下面公开的具体实施例的限制。
[0082] 下面参照图1至图7对本发明实施例作进一步的描述。
[0083] 一种基于自然齿面活动标架的点啮合齿面的齿轮设计方法，所述齿轮包括互相点啮合的第一齿轮和第二齿轮，所述第一齿轮上具有第一齿面，所述第二齿轮上具有第二齿面，所述第一齿面与第二齿面实现点啮合，如图1所示，引入了一个自然齿面活动标架体系，包括固定标架S1[O1；x1y1z1]，固定标架S(2)[O(2)；x(2) y(2) z(2)]，所述设计方法包括以下步骤：
[0084] S10：齿面沿接触点迹线的局部结构的设计
[0085] S11：给定第一齿面∑1：r(1)＝r(1)(u1，ν1)并在所述第一齿面上指定一条接触点迹线L1： 由齿面共轭原理设计第二齿面∑2上与接触点迹线L1共轭的接触点迹线L2： 及所述接触点迹线L2的单位法线矢量参数方程 使所述
第一齿面∑1与所述第二齿面∑2沿着接触点迹线L1啮合传动，并满足式(1)给定的相对运动A
[0086]
[0087] 式中 为使啮合的起始位置与第二齿轮转动的起始位置重合，即 时，第一齿面∑1相对固定坐标系S1[O1；x1y1z1]转过的角度，正如图1所示；ci(i＝2，…，n)为多项式系数，按预定的两齿轮的相对运动规律给出；为第二齿轮与第一齿轮的齿数比±Z2/Z1；是第一齿面∑1自转的角位移； 是第二齿面∑2自转的角位移；
[0088] S12：给定啮合点的位置对安装误差的敏感度ε和齿面接触许用应力值[σ]，再设计第二齿面∑2沿接触点迹线L2的二阶微分结构参数：接触点迹线L2上各点沿其切线方向的法曲率 及其短程挠率 接触点迹线L2的切线在第二齿面∑2切平面上的垂直方向的法曲率 其中设计二阶微分结构参数可以参照现有技术进行，可以参阅吴训成的文章，《基于功能需求的弧齿锥齿轮齿面主动设计与先进制造技术研究》([D].西安：西安交通大学，2000)。
[0089] S20：基于自然活动标架的齿面接触点迹线以外的齿面整体设计
[0090] S21：第一齿面∑1与第二齿面∑2的相对运动速度的曲面运动不变量参数描述[0091] 如图2所示，根据第一齿面∑1：r(1)＝r(1)(u1，υ1)，则参数u1，υ1构成正交参数网线，在第一齿面∑1上任意一接触点M(u1M，υ1M)取定一个活动标架场
[0092]
[0093] 当接触点M沿正交参数网线u1，υ1以速度 移动时，活动标架Sft也随之移动并绕接触点M转动，由微分几何的曲面论的基本定理，活动标架Sft的角速度运动方程如式(2)
[0094] ω1＝ω23e1-ω13e2+ω12n1M (2)
[0095] 式中，
[0096] 如图3所示，设接触点迹线L2是第一齿面∑1在相对运动A下的共轭齿面∑p：r(p)＝r(p)(up，υp)上的一条曲线，其方程为 sp为接触点迹线L2的弧长参数；接触点迹线L2上任意一接触点M的单位切向量为α，曲面在M点的单位法向量为npM，在接触点迹线L2上接触点M处取定一个正交的活动标架Sfp[M；α ν＝npM×α npM]；
[0097] 当接触点M沿接触点迹线L2，以速度 移动时，活动标架Sfp[M；α ν＝npM×α npM]也随之移动并绕接触点M转动，由微分几何的曲面论的基本定理，活动标架Sfp[M；α ν＝npM×α npM]的角速度运动方程如式(3)
[0098]
[0099] 式中， 分别是第二齿面∑2在接触点M沿α方向的短程挠率和法曲率， 是第二齿面∑2在接触点M沿α方向的短程曲率；
[0100] 如图4所示，与固定坐标系S(1)[O(1)；x(1) y(1) z(1)]固结的第一齿面∑1通过相对于与固定坐标系S(2)[O(2)；x(2) y(2) z(2)]固结的共轭齿面∑p的相对运动Ψ包络出第二齿面∑2；设定Lg是形成包络共轭齿面∑p的特征线，则在第一齿面∑1沿接触点迹线L2相对共轭齿面∑p按相对运动Ψ的任意时刻，第一齿面∑1与第二齿面∑2相切于特征线Lg，特征线Lg和接触点迹线L2交于接触点M，而在接触点M，对于第一齿面∑1与第二齿面∑2分别有活动标架Sft和活动标架Sfp，共轭齿面∑p与第二齿面∑2同固结于坐标系S(2)[O(2)；x(2) y(2) z(2)]，共轭齿面∑p是第二齿面∑2的拟合逼近目标，由式(2)和式(3)知第一齿而∑1相对第二齿面∑2的角速度为式(4)
[0101]
[0102] 考虑S(1)、Sft、Sfp和S(2)之间的相对运动，则在接触点M处第一齿面∑1相对第二齿面∑2的速度为式(5)
[0103]
[0104] 按式(5)的方式，将式(4)中的 因子提出来，并将 因子指定为基准速率，不失一般性，令 而由于接触点迹线L2的方程已经由与接触点迹线L1相对运动的共轭条件确定，故 也是已经确定的函数量，则在第一齿面∑1上任意点P(u1，υ1)处，第一齿面∑1相对第二齿面∑2的速度函数表示为式(6)
[0105]
[0106] 式(6)中的r1(ft)表示点P(ut，υt)在活动标架Sft下的位置矢量，Δ是活动标架Sft与活动标架Sfp之间的夹角；
[0107] S22：从第一齿面∑1到第二齿面∑2的运动变换及第二齿面∑2方程的曲面运动不变几何量描述
[0108] 如图5所示，从第一齿面∑1出发，分别由相对运动Ψ展成(拟合)第二齿面∑2和由相对运动A展成共轭齿面∑p的坐标变换过程以及第二齿面∑2的设计思路；图5中两个虚线方框内的内容分别描述的是相对运动Ψ和相对运动A的坐标变换过程，其中相对运动Ψ是如图4所示的在第一齿面∑1和共轭齿面∑p的活动标架下由第一齿面∑1和共轭齿面∑p的运动不变量描述；相对运动A是如图1所示的在点啮合齿轮副固定坐标系下由齿轮副结构参数描述的；沿相对运动Ψ和相对运动A，由图3可知，沿相对运动Ψ和相对运动A，从固定坐标系S(1)到固定坐标系S(2)的变换矩阵为式(7)
[0109]
[0110] 式中
[0111]
[0112]
[0113] 其中， 表示第一齿面∑x上任意接触点M在固定坐标系S(1)[O(1)；x(1) y(1) z(1)]中的位置矢量函数，此接触点M与共轭齿面∑p上的接触点迹线L2上的接触点M关于相对运动Ψ共轭；
[0114] 表示接触点迹线L2在固定坐标系S2[O(2)；x(2) y(2) z(2)]中的矢量参数方程；
[0115] 第一齿面∑1在相对运动Ψ下的第二齿面∑2在固定坐标系S(2)[O(2)；x(2) y(2) z(2)]下的方程为式(8)
[0116] S(2)[O(2)；x(2) y(2) z(2)] (8)
[0117] 式中， 表示第一齿面∑1在固定坐标系S(1)[O(1)；x(1) y(1) z(1)]下的矢量参数方程，N1(u1，υ1)是第一齿面∑1上任意一点P处单位法矢的矢量参数方程，如图2所示，此点P是第一齿面∑1上不同于接触点M的另一任意点；
[0118] S23：运动参数Δ的确定
[0119] 如图6所示，设第一齿面∑1和第二齿面∑2间的，与最小相对曲率对应的相对主方向为αm，沿αm的第一齿面∑1、第二齿面∑2的法曲率为 短程挠率为则沿相对主方向αm的第一齿面∑1、第二齿面∑2的相对短程挠率为零，即满足式(9)
[0120]
[0121] 由欧拉-贝特朗公式，对于共轭齿面∑p沿相对主方向αm满足式(10)、(11)[0122]
[0123]
[0124] 其中，φ2和φ3满足式(12)，
[0125] φ2-φ3＝arc cos(α1·α) (12)
[0126] 对于第一齿面∑1沿相对主方向αm有式(13)
[0127]
[0128] 在式(10)至式(13)中， 分别是第一齿面∑1在接触点M沿ei(i＝1，2)方向的短程挠率和法曲率； 和 分别是共轭齿面∑p在接触点M沿接触点迹线L2的切线方向即单位切矢α方向的短程挠率和法曲率及α方向在公切平面上的垂直方向ν的法曲率；
和 分别是共轭齿面∑p在接触点M沿接触点迹线L1的切线方向即单位切矢α1方
向的短程挠率和法曲率及α1方向在公切平面上的垂直方向ν1的法曲率；由于在步骤S10中已经确定了接触点迹线L1和接触点迹线L2，那么接触点迹线L1的单位切矢α1和接触点迹线L2的单位切矢α就已经确定，故φ2-φ3是已经确定的已知量；
[0129] 基于最终形成的第二齿面∑2需满足步骤S10的设计要求，故对于第二齿面∑2，沿相对主方向αm满足式(14)
[0130]
[0131] 活动标架Sft与活动标架Sfp之间的夹角Δ满足式(15)
[0132] Δ(u1M，υ1M，sp)＝φ2(sp)-φ1(u1M，υ1M) (15)
[0133] 对式(9)～(14)联立求解，用接触点M处第一齿面∑1和第二齿面∑2的各个方向的法曲率和短程挠率表示φ1、φ2和φ3，具体表示为φ1(utM，υtM)，φ2(sp)和φ3(s1)，从而确定相对主方向αm，将φ1(utM，υtM)，φ2(sp)代入式(15)，即可将夹角Δ表达为运动不变量φ1(utM，υtM)，φ2(sp)的函数；
[0134] 将上述方法中求得的Δ分别代入式(6)和(7)，分别求得在任意点P(u1，υ1)处第一齿面∑1相对第二齿面∑2的速度函数和从第一齿面∑1到第二齿面∑2的变换矩阵函数，再将式(6)和式(7)代入式(8)并消去参数u1，就得到第一齿面∑1在相对运动Ψ下的第二齿面∑2在与共轭齿面∑p固结的固定坐标系S(2)[O(2)；x(2) y(2) z(2)]下的矢量函数方程(16)[0135]
[0136] S24：第二齿面∑2在接触点迹线L2以外的区域与第一齿面∑1的间隙
[0137] 如图7所示，在共轭齿面∑p上任意一点Qg点处第二齿面∑2对共轭齿面∑p的间隙满足式(17)
[0138]
[0139] 式中 是共轭齿面∑p上任意一点Qg处单位法矢的矢量参数方程，如果选定了点Qp，即确定了矢量 则式(17)所表达的三个标量方程正好可以求解参数u1、和δk，从而确定第二齿面∑2与共轭齿面∑p在任意位置的间隙δk以及点Qg的位置矢量[0140] S25：第二齿面∑2整体优化设计模型
[0141] 如图8所示，在过任意接触点M的第一齿面∑1在相对运动A下的特征线Lg上，分别位于接触点M的两侧的取两点 和 其对应的曲坐标分别设为(sp，υ11)、(sp，υ12)，特征线Lg上点 和 之间的各点到第二齿面∑2的距离δk与该点到接触点M的距离成正比，设曲线Lgp、点 和 分别是特征线Lg、点 和 在第二齿面∑2上的投影，根据齿面弹性力学实验可知，当第一齿面∑1对第二齿面∑2的间隙 小于一个给定值时，第一齿面∑1上点 和 之间的狭长齿面与第二齿面∑2上 和 点之间的狭长曲面在承受满载的情况下将贴合在一起，并在第二齿面∑2上形成一个瞬时接触区，逐点控制每一个接触点所形成的瞬时接触区的大小、形状和方向显然是不可行的，因为这样设计的第二齿面将很可能不是连续可微的。由图8分析可知：从啮合的起始点M1到啮合的终点M2的整个啮合过程将在第二齿面∑2上连续形成一系列这样的瞬时接触区，并拼接成如图8阴影区域所示的齿面接触区，控制整个齿面啮合过程中在第二齿面∑2上所形成的接触区，并获得连续可微的第二齿面是可能的。因此，齿面优化设计目标是：在齿面承受满载的情况下，接触区的面积应在实际齿面的范围内尽量大，但不能超出实际齿面。因此，控制齿面接触区的面积大小占整个第二齿面面积满足式(18)
[0142]
[0143] 式中 是用参数(sp，υ1)表示的共轭齿面∑p方程，并由式(17)求得的 并代入
方程 得到；而式中 齿面优化
设计目标是：在齿面承受满载的情况下，接触区的面积应在实际齿面的范围内尽量大，但不能超出实际齿面。因此，控制接触区的面积大小占整个齿面面积的2/3至8/9为宜；
[0144] 再按式(19)进行优化设计约束：
[0145]
[0146] 式中δ0是啮合的第一齿面∑1和第二齿面∑2在承受满载的情况下的齿面接触区的边界各点在齿面无载接触时的间隙；
[0147] 求解由式(18)和(19)组成的不等式方程式组，得到两个函数闭区间u1M(sp)∈[u1M1(sp)u1M2(sp)]和υ1M(sp)∈[υ1M1(sp)υ1M2(sp)]，将其代入式(16)即可求得待求的第二齿面∑2的方程的设计区间
[0148] 本发明的有益效果在于：
[0149] 本发明是一种基于齿面自然活动标架的由第一齿面运动不变量参数描述的点啮合齿面整体设计方法，与现有技术相比，本发明引入曲面自然活动标架作为齿面设计的坐标系，将点啮合齿面的设计问题，转换归结为以第一齿面上的包络特征曲线作为样条曲线的曲面拟合逼近问题，通过第一齿面的运动不变量描述和设计点啮合齿面的参数，使点啮合齿面的整体设计摆脱对于特定机床结构的依赖，因而具有通用性和精确性，并能够据此方法对齿面啮合特性进行全局优化设计。
[0150] 以上显示和描述了本发明的基本原理和主要特征及本发明的优点，本行业的技术人员应该了解，以上仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的创造性精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
