发明名称 基于弧齿锥齿轮的四阶传动误差曲线设计方法
摘要 本发明公开了一种基于弧齿锥齿轮的四阶传动误差曲线设计方法,用于解决现有四阶传动误差曲线设计方法存在换齿冲击的技术问题。技术方案是将传动误差方程在齿面参考点处按Taylor级数展开成四阶形式,从四阶Taylor级数展开式方程中分离出小轮相对于参考点啮合时转角的表达式,建立拟合二阶曲线的目标函数。建立控制峰峰值之间差距的目标函数,使得当前齿传动误差曲线向左或向右平移一个小轮齿距的距离时,刚好得到先导齿和后继齿的传动误差曲线,即此时相邻两传动误差曲线峰峰之间的距离为一个小轮齿距的距离,交点刚好落在极大值点,使交点处的切线夹角趋近于180°,消除了换齿时的振动和冲击。
申请公布号 CN103438186A 申请公布日期 2013.12.11
申请号 CN201310376786.2 申请日期 2013.08.27
申请人 西北工业大学 发明人 刘光磊;李德庚;王乐赟
分类号 F16H55/17(2006.01)I 主分类号 F16H55/17(2006.01)I
代理机构 西北工业大学专利中心 61204 代理人 王鲜凯
主权项 1.一种基于弧齿锥齿轮的四阶传动误差曲线设计方法,其特征在于包括以下步骤:步骤一、弧齿锥齿轮为点接触局部共轭传动,传动误差定义为<img file="FDA0000372398850000011.GIF" wi="1312" he="127" />式中,<img file="FDA0000372398850000012.GIF" wi="56" he="59" />是小轮的实际转角,<img file="FDA0000372398850000013.GIF" wi="65" he="59" />是小大轮的实际转角;<img file="FDA0000372398850000014.GIF" wi="202" he="76" />表示齿面参考点啮合时,小轮和大轮的实际转角;z<sub>1</sub>、z<sub>2</sub>分别小轮、大轮的齿数;<img file="FDA0000372398850000015.GIF" wi="216" he="85" />为大轮相对于参考点啮合时的转角,<img file="FDA0000372398850000016.GIF" wi="206" he="83" />为小轮相对于参考点啮合时的转角,<img file="FDA0000372398850000017.GIF" wi="251" he="121" />为按名义传动比确定的大轮名义转角;其中,大轮相对转角<img file="FDA0000372398850000018.GIF" wi="205" he="84" />是小轮相对转角<img file="FDA0000372398850000019.GIF" wi="199" he="86" />的函数,写为<img file="FDA00003723988500000110.GIF" wi="1365" he="84" />步骤二、将大轮实际转角式(2)在参考点处展开成Taylor级数<img file="FDA00003723988500000111.GIF" wi="1154" he="170" />(3)<img file="FDA00003723988500000112.GIF" wi="1330" he="167" />取参考点处的瞬时传动比等于名义传动比,则级数形式的传动误差为<img file="FDA00003723988500000113.GIF" wi="764" he="173" /><img file="FDA00003723988500000114.GIF" wi="537" he="130" />(4)<img file="FDA00003723988500000115.GIF" wi="992" he="164" /><img file="FDA00003723988500000116.GIF" wi="840" he="159" />或简记为<img file="FDA00003723988500000117.GIF" wi="1557" he="113" />式中,A是传动误差曲线的二阶导数或瞬时传动比的一阶导数;B是传动误差曲线的三阶导数或瞬时传动比的二阶导数;C是传动误差曲线的四阶导数或瞬时传动比的三阶导数;<img file="FDA00003723988500000118.GIF" wi="309" he="109" />是高于五阶各项之和;当仅取式(5)的右端第一项时,即得到Litvin提出的二阶传动误差曲线<img file="FDA0000372398850000021.GIF" wi="1318" he="89" />当取式(5)至右端第三项时,即得到Stadtfeld四阶传动误差曲线<img file="FDA0000372398850000022.GIF" wi="1358" he="95" />步骤三、采用笛卡尔直角坐标系描述四阶传动误差曲线,笛卡尔直角坐标系的横坐标为小轮转角,纵坐标为大轮传动误差;当前齿对啮合的齿面参考点为坐标原点,其他齿对的齿面参考点依据齿距角沿横坐标轴向两侧平移延拓而得;定义相邻传动误差曲线交点处的传动误差为传动误差曲线幅值,其值根据大轮转动角加速度和跃度在2"~20"范围内选取;步骤四、四阶传动误差曲线设计方法;1)四阶传动误差曲线生成;选择若干预控点A<sub>1</sub>,A<sub>2</sub>,…A<sub>n</sub>,采用最小二乘法拟合四阶传动误差曲线即获得四阶传动误差曲线的三个系数A、B和C;2)预控点的选取;选择10个预控点A<sub>i</sub>(φ<sub>1i</sub>,Δφ<sub>2i</sub>),i=1,2,...,10;其中,预控点A<sub>2</sub>、A<sub>5</sub>、A<sub>6</sub>和A<sub>9</sub>位于横坐标轴上,即Δφ<sub>2j</sub>=0,j=2,5,6,9,φ<sub>15</sub>=-(0.05~0.1)φ<sub>1P</sub>,φ<sub>1P</sub>为小齿轮齿距角,φ<sub>16</sub>=-φ<sub>15</sub>;φ<sub>12</sub>=-(0.6~0.7)φ<sub>1P</sub>,φ<sub>19</sub>=-φ<sub>12</sub>;A<sub>3</sub>和A<sub>8</sub>用于调节传动误差的极大值,初步计算时,A<sub>3</sub>和A<sub>8</sub>对称分布在纵轴两侧,φ<sub>13</sub>=-(0.4~0.5)φ<sub>1P</sub>,φ<sub>18</sub>=-φ<sub>13</sub>,纵坐标Δφ<sub>23</sub>=Δφ<sub>28</sub>=2′′~20′′;A<sub>4</sub>位于<img file="FDA0000372398850000023.GIF" wi="109" he="89" />的中部,A<sub>7</sub>位于<img file="FDA0000372398850000024.GIF" wi="110" he="89" />的中部,A<sub>1</sub>和A<sub>10</sub>用于控制传动误差曲线波浪形顶部以下的曲线走向,φ<sub>11</sub>=-(0.9~0.98)φ<sub>1P</sub>,φ<sub>110</sub>=-φ<sub>11</sub>,Δφ<sub>21</sub>=Δφ<sub>210</sub>=(1.5~3)Δφ<sub>23</sub>;3)系数A、B和C的确定;从式(7)右端提出<img file="FDA0000372398850000025.GIF" wi="237" he="92" />得到<img file="FDA0000372398850000026.GIF" wi="1387" he="111" />从中抽取出二阶曲线<img file="FDA0000372398850000027.GIF" wi="99" he="81" />方程,即<img file="FDA0000372398850000028.GIF" wi="1239" he="93" />用前面的10个预控点,建立拟合二阶曲线的目标函数<img file="FDA0000372398850000031.GIF" wi="1150" he="143" />式中,Δφ<sub>2'</sub>(φ<sub>1i</sub>)是二阶曲线<img file="FDA0000372398850000032.GIF" wi="96" he="78" />在φ<sub>1i</sub>处的函数值;Δφ<sub>2i</sub>是预控点A<sub>i</sub>处的误差控制值;式(10)中有三个未知参数A、B和C,按照最小二乘法,计算目标函数F对A、B和C的偏导数<img file="FDA0000372398850000039.GIF" wi="282" he="118" />令其等于零<maths num="0001"><![CDATA[<math><mrow><mfenced open='{' close=''><mtable><mtr><mtd><mfrac><mrow><mo>&PartialD;</mo><mi>F</mi></mrow><mrow><mo>&PartialD;</mo><mi>A</mi></mrow></mfrac><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><mfrac><mrow><mo>&PartialD;</mo><mi>F</mi></mrow><mrow><mo>&PartialD;</mo><mi>B</mi></mrow></mfrac><mo>=</mo><mn>0</mn></mtd></mtr><mtr><mtd><mfrac><mrow><mo>&PartialD;</mo><mi>F</mi></mrow><mrow><mo>&PartialD;</mo><mi>C</mi></mrow></mfrac><mo>=</mo><mn>0</mn></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mi></mi><mn>11</mn><mo>)</mo></mrow></mrow></math>]]></maths>求解这个方程组,即得二阶传动误差曲线中的三个系数A、B和C;4)Stadtfeld四阶传动误差曲线;计算式(7)的极大值<img file="FDA0000372398850000034.GIF" wi="212" he="86" />将传动误差曲线向下平移这个极值的距离,即<img file="FDA0000372398850000035.GIF" wi="1604" he="94" />则形成单齿对啮合时的传动误差曲线;当预控点A<sub>3</sub>和A<sub>8</sub>的纵坐标值近乎相等时,就得到Stadtfeld四阶传动误差曲线;Stadtfeld当前齿对的啮入点和啮出点基本以坐标原点为对称,所得接触印痕基本上位于齿高和齿宽的中部;5)改进的传动误差曲线;继续保证预控点A<sub>3</sub>和A<sub>8</sub>的纵坐标值相等,以便得到改进的四阶传动误差曲线;对峰峰点之间的距离进行控制;建立控制峰峰之间差距的目标函数;拆分10个点的横纵坐标,将四阶传动误差降为二阶,得到ax<sup>2</sup>+bx+c=0,求得<img file="FDA0000372398850000036.GIF" wi="399" he="132" />对应的分别是A<sub>3</sub>和A<sub>8</sub>的横坐标;令<img file="FDA0000372398850000037.GIF" wi="619" he="79" />并且使右峰向左侧平移一个齿距,此时,两峰间距离为:<img file="FDA0000372398850000038.GIF" wi="588" he="149" />由于预控点的横坐标是由a、b、c三个参数所表示的,对预控点的横坐标中的参数进行优化,不断调整预控点的横坐标,使得相对应的两个右峰或两个左峰之间的距离刚好为一个小轮齿距的距离,即Δx<sub>peak</sub>=0,此时x=0,相邻传动误差曲线的两凸峰重合,得到改进后的传动误差曲线。
地址 710072 陕西省西安市友谊西路127号