CN107729626A - 一种基于概率分布的齿轮点蚀模型建模方法 - Google Patents

一种基于概率分布的齿轮点蚀模型建模方法 Download PDF

Info

Publication number
CN107729626A
CN107729626A CN201710881432.1A CN201710881432A CN107729626A CN 107729626 A CN107729626 A CN 107729626A CN 201710881432 A CN201710881432 A CN 201710881432A CN 107729626 A CN107729626 A CN 107729626A
Authority
CN
China
Prior art keywords
mrow
msub
alpha
mfrac
msup
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.)
Granted
Application number
CN201710881432.1A
Other languages
English (en)
Other versions
CN107729626B (zh
Inventor
雷亚国
王德龙
刘宗尧
杨筱
刘欢
林京
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Jiaotong University
Original Assignee
Xian Jiaotong University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xian Jiaotong University filed Critical Xian Jiaotong University
Priority to CN201710881432.1A priority Critical patent/CN107729626B/zh
Publication of CN107729626A publication Critical patent/CN107729626A/zh
Application granted granted Critical
Publication of CN107729626B publication Critical patent/CN107729626B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/17—Mechanical parametric or variational design
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F16—ENGINEERING ELEMENTS AND UNITS; GENERAL MEASURES FOR PRODUCING AND MAINTAINING EFFECTIVE FUNCTIONING OF MACHINES OR INSTALLATIONS; THERMAL INSULATION IN GENERAL
    • F16H—GEARING
    • F16H57/00—General details of gearing
    • F16H2057/0087—Computer aided design [CAD] specially adapted for gearing features ; Analysis of gear systems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/04—Ageing analysis or optimisation against ageing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Gears, Cams (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

一种基于概率分布的齿轮点蚀模型建模方法，首先建立了齿面概率点蚀平面分布模型，点蚀产生按照沿齿宽方向随机分布、沿齿高方向正态分布；其次，采用图像处理方法来获取啮合线长度变化；然后，建立啮合线长度变化值与压力角之间的一一对应关系；最后，综合考虑赫兹接触刚度、弯曲刚度、剪切刚度和轴向压缩刚度，精确计算齿轮时变啮合刚度，求解齿轮系统内部激励；本发明克服了传统点蚀模型与实际相差大、无法处理点蚀复杂分布等缺点，同时提高了模型计算精度，实现点蚀齿轮系统内部激励的精确计算，完成基于概率分布的齿轮点蚀模型建模。

Description

一种基于概率分布的齿轮点蚀模型建模方法
技术领域
本发明属于齿轮系统动态建模技术领域，具体涉及一种基于概率分布的齿轮点蚀模型建模方法。
背景技术
齿轮传动系统以其适用范围广、传动比准确、传动效率高、结构紧凑、可靠寿命长等诸多优点，被广泛应用于仪器、仪表、冶金、矿山等行业。然而，由于齿轮系统存在制造误差、安装误差，以及常常工作于过载、欠润滑等恶劣的工况下，导致齿面点蚀、剥落、磨损、齿根裂纹等机械故障频繁出现。作为一种常见的齿轮故障类型，齿面点蚀会导致齿轮啮合刚度减小；任其发展，将导致齿面剥落，甚至严重磨损，进而给整个传动链带来严重的安全隐患，造成严重的经济损失，甚至人员伤亡。
在齿轮系统的众多激励形式中，内部刚度激励是齿轮传动的重要激励形式，准确计算时变啮合刚度是齿轮系统动力学建模的首要任务，齿轮系统动力学模型的准确性直接影响时变啮合刚度计算的精度。点蚀发生位置的不确定性给齿轮时变啮合刚度的研究带来诸多不便，传统模型中，研究人员将点蚀进行简化，主要存在以下缺点：(1)点蚀位置与实际不符，点蚀通常集中在节线偏下的位置上，但在简化时，多将点蚀坑中心定位于节线上；(2)点蚀形状与实际不符，点蚀坑的形状是不规则的，传统模型多将点蚀坑形状简化为圆形、椭圆形和矩形；(3)点蚀数量与实际不符，为简化计算，采用面积大，数量少的点蚀坑来模拟点蚀。综上，已有模型与实际相差较远，必然导致齿轮时变啮合刚度计算结果不精确，造成计算所得结果与实际结果存在较大差异。
发明内容
为了克服上述现有技术的缺点，本发明的目的在于提供一种基于概率分布的齿轮点蚀模型建模方法，采用圆柱凹坑模拟点蚀，点蚀位置基于概率分布，还原点蚀分布对时变啮合刚度的影响，减小传统方法带来的计算误差，实现对点蚀齿轮时变啮合刚度的精确计算。
为了达到上述目的，本发明采取的技术方案为：
一种基于概率分布的齿轮点蚀模型建模方法，包括以下步骤：
1)建立齿面概率点蚀平面分布模型，点蚀故障在啮合齿面上的分布具有规律性：点蚀沿轮齿齿宽方向服从随机分布，沿轮齿齿高方向服从正态分布，正态分布的具体参数如下式所示，
σ＝h/6
μ＝r-σ/2
式中，μ为数学期望，即齿高方向点蚀分布的中心位置到齿轮中心的距离，σ为齿高方向正态分布标准差，h为发生点蚀的轮齿齿高，r为点蚀齿轮节圆半径，点蚀故障的严重程度由点蚀面积占齿面面积的比例决定，划分为正常、一级点蚀、二级点蚀、三级点蚀共4个等级，一级点蚀的点蚀面积占齿面面积的比例不超过3％，二级点蚀的点蚀面积占齿面面积的比例为3％-15％，三级点蚀的点蚀面积占齿面面积的比例为15％-50％，第i(i＝1,2,3)等级新产生的点蚀坑的数量ni为：
式中，L为齿宽，Pi为第i等级点蚀面积占齿面面积的比例，ri为第i等级点蚀圆柱凹坑的半径，综合齿面点蚀概率和各等级点蚀坑数量，得到各等级齿面概率点蚀平面分布模型，保存为不失真的位图；
2)依次读取步骤1)中各等级齿面概率点蚀平面分布模型图像各像素点灰度值，灰度值grey≥128，此像素点不在点蚀区域内；grey＜128，认为此像素点在点蚀区域内，分别读取各等级齿面概率点蚀平面分布模型图像，获取点蚀引起的接触线长度减小值ΔL；
3)建立ΔL与啮合压力角α之间的联系，计算齿廓上与α相对应的啮合点对中间面投影与节圆投影之间的距离Δx1，
式中，rbp为主动轮基圆半径，α2为半基圆齿厚角，m为齿轮模数，Np为主动轮齿数，同样地，计算ΔL对应位置到节圆投影之间的距离Δx2，
Δx2＝x-r
式中，x为ΔL的位置到圆心的距离，r为节圆半径，寻找与Δx1最接近的Δx2，则此Δx2对应的ΔL与此Δx1对应的α之间就建立了一一对应关系；
4)建立刚度求解方法，求解齿轮系统内部激励，齿轮时变啮合刚度包括：赫兹接触刚度kh、弯曲刚度kb、剪切刚度ks、轴向压缩刚度ka，齿轮时变啮合刚度计算公式如下：
式中，下标的1和2分别代表主动轮和从动轮，j＝1和j＝2分别代表参与啮合的第一对齿轮和第二对齿轮，其中，kh，kb，ks，ka的计算公式如下，
式中，rb为基圆半径，rr为齿根圆半径，α1为啮合压力角，E为齿轮材料杨氏模量，RR为齿根过渡圆弧半径，x1为过渡部分对中见面投影到圆心的距离，H为点蚀坑深度，υ为齿轮材料泊松比，由此，便得到了齿轮系统内部激励，求得齿轮时变啮合刚度，完成基于概率分布的齿轮点蚀模型建模。
本发明的有益效果为：
本发明考虑了实际中点蚀产生按照沿齿宽方向随机分布、沿齿高方向正态分布的特点，并采用概率函数进行模拟；提出了采用图像处理的方式来获取啮合线长度变化的方法，建立了啮合线长度变化值与积分变量之间的联系，摆脱了传统方法中对模拟点蚀区域必须进行分段的思想，大大简化计算过程；此外，综合考虑赫兹接触刚度、弯曲刚度、剪切刚度和轴向压缩刚度，实现了齿轮时变啮合刚度的精确计算，求解齿轮系统内部激励。因此，一种基于概率分布的齿轮点蚀模型建模方法克服了传统点蚀模型与实际相差大、无法处理点蚀复杂分布情况、计算方法复杂等缺点，提高了模型计算精度，实现点蚀齿轮系统内部激励的精确计算，完成基于概率分布的齿轮点蚀模型建模。
附图说明
图1为本发明方法的流程图。
图2为概率点蚀分布模型示意图。
图3为点蚀四个等级点蚀坑分布图，图(a)为正常齿面，图(b)为一级点蚀齿面，图(c)为二级点蚀齿面，图(d)为三级点蚀齿面。
图4为模型时变啮合刚度图，图(a)为四个等级对比图，图(b)为图(a)中单齿啮合区放大图，图(c)为图(a)中双齿啮合区放大图。
图5为点蚀四个等级的有限元模型图，图(a)为正常模型，图(b)为一级点蚀模型，图(c)为二级点蚀模型，图(d)为三级点蚀模型。
图6为有限元模型计算结果与计算结果对比图，图(a)为正常齿面有限元法与势能法计算所得刚度对比图，图(b)为一级点蚀齿面有限元法与势能法计算所得刚度对比图，图(c)为二级点蚀齿面有限元法与势能法计算所得刚度对比图，图(d)为三级点蚀齿面有限元法与势能法计算所得刚度对比图。
具体实施方式
下面结合附图对本发明做进一步详细描述：
参照图1，一种基于概率分布的齿轮点蚀模型建模方法，包括以下步骤：
1)建立齿面概率点蚀平面分布模型，忽略齿轮加工和安装误差，如图2所示，点蚀故障在啮合齿面上的分布具有规律性：点蚀沿轮齿齿宽方向服从随机分布，沿轮齿齿高方向服从正态分布，正态分布的具体参数如下式所示，
σ＝h/6
μ＝r-σ/2
式中，μ为数学期望，即齿高方向点蚀分布的中心位置到齿轮中心的距离，σ为齿高方向正态分布标准差，h为发生点蚀的轮齿齿高，r为点蚀齿轮节圆半径，点蚀故障的严重程度由点蚀面积占齿面面积的比例决定，划分为正常、一级点蚀、二级点蚀、三级点蚀共4个等级，一级点蚀的点蚀面积占齿面面积的比例不超过3％，二级点蚀的点蚀面积占齿面面积的比例为3％-15％，三级点蚀的点蚀面积占齿面面积的比例为15％-50％，第i(i＝1,2,3)等级新产生的点蚀坑的数量ni为：
式中，L为齿宽，Pi为第i等级点蚀面积占齿面面积的比例，ri为第i等级点蚀圆柱凹坑的半径，综合齿面点蚀概率和各等级点蚀坑数量，得到各等级齿面概率点蚀平面分布模型，如图3所示，保存为不失真的位图；
2)依次读取步骤1)中各等级齿面概率点蚀平面分布模型图像各像素点灰度值，灰度值grey≥128，此像素点不在点蚀区域内；grey＜128，认为此像素点在点蚀区域内，分别读取各等级齿面概率点蚀平面分布模型图像，获取点蚀引起的接触线长度减小值ΔL；
3)建建立ΔL与啮合压力角α之间的联系，计算齿廓上与α相对应的啮合点对中间面投影与节圆投影之间的距离Δx1，
式中，rbp为主动轮基圆半径，α2为半基圆齿厚角，m为齿轮模数，Np为主动轮齿数，同样地，计算ΔL对应位置到节圆投影之间的距离Δx2，
Δx2＝x-r
式中，x为ΔL的位置到圆心的距离，r为节圆半径，寻找与Δx1最接近的Δx2，则此Δx2对应的ΔL与此Δx1对应的α之间就建立了一一对应关系；
4)建立刚度求解方法，求解齿轮系统内部激励，齿轮时变啮合刚度包括：赫兹接触刚度kh，弯曲刚度kb，剪切刚度ks，轴向压缩刚度ka，齿轮时变啮合刚度计算公式如下：
式中，下标的1和2分别代表主动轮和从动轮，j＝1和j＝2分别代表参与啮合的第一对齿轮和第二对齿轮，其中，kh，kb，ks，ka的计算公式如下，
式中，rb为基圆半径，rr为齿根圆半径，α1为啮合压力角，E为齿轮材料杨氏模量，RR为齿根过渡圆弧半径，x1为过渡部分对中见面投影到圆心的距离，H为点蚀坑深度，υ为齿轮材料泊松比，由此，便得到了齿轮系统内部激励，求得齿轮时变啮合刚度，如图4所示，完成基于概率分布的齿轮点蚀模型建模。
为了进一步证明该基于概率分布的齿轮点蚀模型建模方法的准确性，建立齿轮点蚀有限元模型求解齿轮内部激励，齿轮结构参数如表1所示。
表1：齿轮结构参数
结合各等级点蚀分布情况，利用Solidworks三维建模软件建立基于概率分布的齿轮点蚀模型，将基于概率分布的齿轮点蚀模型导入ANSYS软件进行有限元分析。采用ANSYS软件分析时，单元类型为Solid 185，为减少有限元分析计算量，建立基于概率分布的齿轮点蚀模型时，只建立基础部分和包括点蚀轮齿在内的三个轮齿，建立的基于概率分布的齿轮点蚀有限元模型如图5所示。
基于概率分布的齿轮点蚀有限元模型中，采用静力学分析。将每个啮合周期均分为20份，在每个位置上分别进行静力学分析计算，最终将多个点分别计算得到的啮合刚度进行拟合，得到齿轮时变啮合刚度。进行静力学分析时，主要工作有：网格划分，添加约束和载荷，计算求解。划分网格时，正常轮齿划分六面体扫略单元，由于故障轮齿点蚀坑的不规则性，将故障轮齿划分为四面体单元进行求解。约束条件有三：从动轮轮轴固定；主动轮轮轴施加圆柱约束，保留切向自由度；两齿面之间为接触约束。载荷施加在主动轮轮轴上，求解模型应力、应变以及位移，根据静变形分别计算每个角度对应的啮合刚度。四个等级下基于概率分布的齿轮点蚀有限元模型数值解与能量法解析解的对比如图6所示，从图中可以看出，能量法解析解与有限元模型数值解相吻合，验证了一种基于概率分布的齿轮点蚀模型建模方法的准确性。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施方式仅限于此，对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单的推演或替换，都应当视为属于本发明由所提交的权利要求书确定专利保护范围。

Claims (1)

1.一种基于概率分布的齿轮点蚀模型建模方法，其特征在于，包括以下步骤：
1)建立齿面概率点蚀平面分布模型，点蚀故障在啮合齿面上的分布具有规律性：点蚀沿轮齿齿宽方向服从随机分布，沿轮齿齿高方向服从正态分布，正态分布的具体参数如下式所示，
σ＝h/6
μ＝r-σ/2
式中，μ为数学期望，即齿高方向点蚀分布的中心位置到齿轮中心的距离，σ为齿高方向正态分布标准差，h为发生点蚀的轮齿齿高，r为点蚀齿轮节圆半径，点蚀故障的严重程度由点蚀面积占齿面面积的比例决定，划分为正常、一级点蚀、二级点蚀、三级点蚀共4个等级，一级点蚀的点蚀面积占齿面面积的比例不超过3％，二级点蚀的点蚀面积占齿面面积的比例为3％-15％，三级点蚀的点蚀面积占齿面面积的比例为15％-50％，第i(i＝1,2,3)等级新产生的点蚀坑的数量ni为：
<mrow> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mi>h</mi> <mo>&amp;CenterDot;</mo> <mi>L</mi> <mo>&amp;CenterDot;</mo> <msub> <mi>P</mi> <mn>1</mn> </msub> </mrow> <mrow> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <msubsup> <mi>r</mi> <mn>1</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> </mrow>
<mrow> <msub> <mi>n</mi> <mn>2</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mi>h</mi> <mo>&amp;CenterDot;</mo> <mi>L</mi> <mo>&amp;CenterDot;</mo> <msub> <mi>P</mi> <mn>2</mn> </msub> <mo>-</mo> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>&amp;CenterDot;</mo> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <msubsup> <mi>r</mi> <mn>2</mn> <mn>2</mn> </msubsup> </mrow> <mrow> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <msubsup> <mi>r</mi> <mn>1</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> </mrow>
<mrow> <msub> <mi>n</mi> <mn>3</mn> </msub> <mo>=</mo> <mfrac> <mrow> <mi>h</mi> <mo>&amp;CenterDot;</mo> <mi>L</mi> <mo>&amp;CenterDot;</mo> <msub> <mi>P</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>n</mi> <mn>1</mn> </msub> <mo>&amp;CenterDot;</mo> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <msubsup> <mi>r</mi> <mn>3</mn> <mn>2</mn> </msubsup> <mo>-</mo> <msub> <mi>n</mi> <mn>2</mn> </msub> <mo>&amp;CenterDot;</mo> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <msubsup> <mi>r</mi> <mn>2</mn> <mn>2</mn> </msubsup> </mrow> <mrow> <mi>&amp;pi;</mi> <mo>&amp;CenterDot;</mo> <msubsup> <mi>r</mi> <mn>1</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> </mrow>
式中，L为齿宽，Pi为第i等级点蚀面积占齿面面积的比例，ri为第i等级点蚀圆柱凹坑的半径，综合齿面点蚀概率和各等级点蚀坑数量，得到各等级齿面概率点蚀平面分布模型，保存为不失真的位图；
2)依次读取步骤1)中各等级齿面概率点蚀平面分布模型图像各像素点灰度值，灰度值grey≥128，此像素点不在点蚀区域内；grey＜128，认为此像素点在点蚀区域内，分别读取各等级齿面概率点蚀平面分布模型图像，获取点蚀引起的接触线长度减小值ΔL；
3)建立ΔL与啮合压力角α之间的联系，计算齿廓上与α相对应的啮合点对中间面投影与节圆投影之间的距离Δx1，
<mrow> <msub> <mi>&amp;Delta;x</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>r</mi> <mrow> <mi>b</mi> <mi>p</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mrow> <mo>(</mo> <mi>&amp;alpha;</mi> <mo>+</mo> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;alpha;</mi> <mo>+</mo> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mi>&amp;alpha;</mi> <mo>&amp;rsqb;</mo> <mo>-</mo> <mfrac> <mrow> <mi>m</mi> <mo>&amp;CenterDot;</mo> <msub> <mi>N</mi> <mi>p</mi> </msub> </mrow> <mn>2</mn> </mfrac> <mo>&amp;CenterDot;</mo> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mfrac> <mi>&amp;pi;</mi> <mrow> <mn>2</mn> <mo>&amp;CenterDot;</mo> <msub> <mi>N</mi> <mi>p</mi> </msub> </mrow> </mfrac> <mo>)</mo> </mrow> </mrow>
式中，rbp为主动轮基圆半径，α2为半基圆齿厚角，m为齿轮模数，Np为主动轮齿数，同样地，计算ΔL对应位置到节圆投影之间的距离Δx2，
Δx2＝x-r
式中，x为ΔL的位置到圆心的距离，r为节圆半径，寻找与Δx1最接近的Δx2，则此Δx2对应的ΔL与此Δx1对应的α之间就建立了一一对应关系；
4)建立刚度求解方法，求解齿轮系统内部激励，齿轮时变啮合刚度包括：赫兹接触刚度kh、弯曲刚度kb、剪切刚度ks、轴向压缩刚度ka，齿轮时变啮合刚度计算公式如下：
<mrow> <mi>k</mi> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mfrac> <mn>1</mn> <mrow> <mfrac> <mn>1</mn> <msub> <mi>k</mi> <mi>h</mi> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>k</mi> <mrow> <mi>b</mi> <mn>1</mn> </mrow> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>k</mi> <mrow> <mi>a</mi> <mn>1</mn> </mrow> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>k</mi> <mrow> <mi>s</mi> <mn>1</mn> </mrow> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>k</mi> <mrow> <mi>b</mi> <mn>2</mn> </mrow> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>k</mi> <mrow> <mi>a</mi> <mn>2</mn> </mrow> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>k</mi> <mrow> <mi>s</mi> <mn>2</mn> </mrow> </msub> </mfrac> </mrow> </mfrac> </mtd> </mtr> <mtr> <mtd> <mrow> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>j</mi> <mo>=</mo> <mn>1</mn> </mrow> <mn>2</mn> </munderover> <mfrac> <mn>1</mn> <mrow> <mfrac> <mn>1</mn> <msub> <mi>k</mi> <mrow> <mi>h</mi> <mo>.</mo> <mi>j</mi> </mrow> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>k</mi> <mrow> <mi>b</mi> <mn>1.</mn> <mi>j</mi> </mrow> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>k</mi> <mrow> <mi>a</mi> <mn>1.</mn> <mi>j</mi> </mrow> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>k</mi> <mrow> <mi>s</mi> <mn>1.</mn> <mi>j</mi> </mrow> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>k</mi> <mrow> <mi>b</mi> <mn>2.</mn> <mi>j</mi> </mrow> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>k</mi> <mrow> <mi>a</mi> <mn>2.</mn> <mi>j</mi> </mrow> </msub> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msub> <mi>k</mi> <mrow> <mi>s</mi> <mn>2.</mn> <mi>j</mi> </mrow> </msub> </mfrac> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced> </mrow>
式中，下标的1和2分别代表主动轮和从动轮，j＝1和j＝2分别代表参与啮合的第一对齿轮和第二对齿轮，其中，kh，kb，ks，ka的计算公式如下，
<mrow> <msub> <mi>k</mi> <mi>h</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mi>&amp;pi;</mi> <mi>E</mi> <mi>L</mi> </mrow> <mrow> <mn>4</mn> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msup> <mi>v</mi> <mn>2</mn> </msup> <mo>)</mo> </mrow> </mrow> </mfrac> </mrow>
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <mfrac> <mn>1</mn> <msub> <mi>k</mi> <mi>b</mi> </msub> </mfrac> <mo>=</mo> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mrow> <msub> <mi>r</mi> <mi>b</mi> </msub> <mo>-</mo> <msub> <mi>r</mi> <mi>r</mi> </msub> </mrow> </msubsup> <mfrac> <mrow> <mn>3</mn> <msup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <msub> <mi>r</mi> <mi>b</mi> </msub> <mo>-</mo> <msub> <mi>r</mi> <mi>b</mi> </msub> <msub> <mi>cos&amp;alpha;</mi> <mn>1</mn> </msub> <msub> <mi>cos&amp;alpha;</mi> <mn>2</mn> </msub> <mo>+</mo> <msub> <mi>x</mi> <mn>1</mn> </msub> <msub> <mi>cos&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mn>2</mn> </msup> </mrow> <mrow> <mn>2</mn> <mi>E</mi> <mi>L</mi> <msup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <msub> <mi>r</mi> <mi>b</mi> </msub> <msub> <mi>sin&amp;alpha;</mi> <mn>2</mn> </msub> <mo>+</mo> <mi>R</mi> <mi>R</mi> <mo>-</mo> <msqrt> <mrow> <msup> <mi>RR</mi> <mn>2</mn> </msup> <mo>-</mo> <msubsup> <mi>x</mi> <mn>1</mn> <mn>2</mn> </msubsup> </mrow> </msqrt> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mn>3</mn> </msup> </mrow> </mfrac> <msub> <mi>dx</mi> <mn>1</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <msubsup> <mo>&amp;Integral;</mo> <mrow> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> </msubsup> <mfrac> <mrow> <mn>3</mn> <msup> <mrow> <mo>{</mo> <mrow> <mn>1</mn> <mo>+</mo> <msub> <mi>cos&amp;alpha;</mi> <mn>1</mn> </msub> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <mi>sin</mi> <mi>&amp;alpha;</mi> <mo>-</mo> <mi>cos</mi> <mi>&amp;alpha;</mi> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mrow> <mo>}</mo> </mrow> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <mi>cos</mi> <mi>&amp;alpha;</mi> </mrow> <mrow> <mn>2</mn> <mi>E</mi> <mi>L</mi> <msup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mi>sin</mi> <mi>&amp;alpha;</mi> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <mi>cos</mi> <mi>&amp;alpha;</mi> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mn>3</mn> </msup> </mrow> </mfrac> <mi>d</mi> <mi>&amp;alpha;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msubsup> <mo>&amp;Integral;</mo> <mrow> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> </msubsup> <mfrac> <mrow> <mn>3</mn> <msup> <mrow> <mo>{</mo> <mrow> <mn>1</mn> <mo>+</mo> <msub> <mi>cos&amp;alpha;</mi> <mn>1</mn> </msub> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <mi>sin</mi> <mi>&amp;alpha;</mi> <mo>-</mo> <mi>cos</mi> <mi>&amp;alpha;</mi> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mrow> <mo>}</mo> </mrow> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <mi>cos</mi> <mi>&amp;alpha;</mi> </mrow> <mrow> <mn>2</mn> <mi>E</mi> <mi>L</mi> <msup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mi>sin</mi> <mi>&amp;alpha;</mi> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <mi>cos</mi> <mi>&amp;alpha;</mi> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mn>3</mn> </msup> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>L</mi> </mrow> <mi>L</mi> </mfrac> <mi>d</mi> <mi>&amp;alpha;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <msubsup> <mo>&amp;Integral;</mo> <mrow> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> </msubsup> <mfrac> <mrow> <mn>3</mn> <msup> <mrow> <mo>{</mo> <mrow> <mn>1</mn> <mo>+</mo> <msub> <mi>cos&amp;alpha;</mi> <mn>1</mn> </msub> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <mi>sin</mi> <mi>&amp;alpha;</mi> <mo>-</mo> <mi>cos</mi> <mi>&amp;alpha;</mi> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mrow> <mo>}</mo> </mrow> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <mi>cos</mi> <mi>&amp;alpha;</mi> </mrow> <mrow> <mn>2</mn> <mi>E</mi> <mi>L</mi> <msup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mi>sin</mi> <mi>&amp;alpha;</mi> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <mi>cos</mi> <mi>&amp;alpha;</mi> <mo>-</mo> <mfrac> <mi>H</mi> <mrow> <mn>2</mn> <msub> <mi>r</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mn>3</mn> </msup> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>L</mi> </mrow> <mi>L</mi> </mfrac> <mi>d</mi> <mi>&amp;alpha;</mi> </mrow> </mtd> </mtr> </mtable> </mfenced>
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <mfrac> <mn>1</mn> <msub> <mi>k</mi> <mi>s</mi> </msub> </mfrac> <mo>=</mo> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mrow> <msub> <mi>r</mi> <mi>b</mi> </msub> <mo>-</mo> <msub> <mi>r</mi> <mi>r</mi> </msub> </mrow> </msubsup> <mfrac> <mrow> <mn>1.2</mn> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>&amp;upsi;</mi> </mrow> <mo>)</mo> </mrow> <msup> <mi>cos</mi> <mn>2</mn> </msup> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> <mrow> <mi>E</mi> <mi>L</mi> <msup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <msub> <mi>r</mi> <mi>b</mi> </msub> <msub> <mi>sin&amp;alpha;</mi> <mn>2</mn> </msub> <mo>+</mo> <mi>R</mi> <mi>R</mi> <mo>-</mo> <msqrt> <mrow> <msup> <mi>RR</mi> <mn>2</mn> </msup> <mo>-</mo> <msubsup> <mi>x</mi> <mn>1</mn> <mn>2</mn> </msubsup> </mrow> </msqrt> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mn>3</mn> </msup> </mrow> </mfrac> <msub> <mi>dx</mi> <mn>1</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <msubsup> <mo>&amp;Integral;</mo> <mrow> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> </msubsup> <mfrac> <mrow> <mn>1.2</mn> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>&amp;upsi;</mi> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <msup> <mi>cos&amp;alpha;cos</mi> <mn>2</mn> </msup> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> <mrow> <mi>E</mi> <mi>L</mi> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mi>sin</mi> <mi>&amp;alpha;</mi> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <mi>cos</mi> <mi>&amp;alpha;</mi> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mrow> </mfrac> <mi>d</mi> <mi>&amp;alpha;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msubsup> <mo>&amp;Integral;</mo> <mrow> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> </msubsup> <mfrac> <mrow> <mn>1.2</mn> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>&amp;upsi;</mi> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <msup> <mi>cos&amp;alpha;cos</mi> <mn>2</mn> </msup> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> <mrow> <mi>E</mi> <mi>L</mi> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mi>sin</mi> <mi>&amp;alpha;</mi> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <mi>cos</mi> <mi>&amp;alpha;</mi> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>L</mi> </mrow> <mi>L</mi> </mfrac> <mi>d</mi> <mi>&amp;alpha;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <msubsup> <mo>&amp;Integral;</mo> <mrow> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> </msubsup> <mfrac> <mrow> <mn>1.2</mn> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>+</mo> <mi>&amp;upsi;</mi> </mrow> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <msup> <mi>cos&amp;alpha;cos</mi> <mn>2</mn> </msup> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> <mrow> <mi>E</mi> <mi>L</mi> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mi>sin</mi> <mi>&amp;alpha;</mi> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <mi>cos</mi> <mi>&amp;alpha;</mi> <mo>-</mo> <mfrac> <mi>H</mi> <mrow> <mn>2</mn> <msub> <mi>r</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>L</mi> </mrow> <mi>L</mi> </mfrac> <mi>d</mi> <mi>&amp;alpha;</mi> </mrow> </mtd> </mtr> </mtable> </mfenced>
<mfenced open = "" close = ""> <mtable> <mtr> <mtd> <mrow> <mfrac> <mn>1</mn> <msub> <mi>k</mi> <mi>a</mi> </msub> </mfrac> <mo>=</mo> <msubsup> <mo>&amp;Integral;</mo> <mn>0</mn> <mrow> <msub> <mi>r</mi> <mi>b</mi> </msub> <mo>-</mo> <msub> <mi>r</mi> <mi>r</mi> </msub> </mrow> </msubsup> <mfrac> <mrow> <msup> <mi>sin</mi> <mn>2</mn> </msup> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> <mrow> <mn>2</mn> <mi>E</mi> <mi>L</mi> <msup> <mrow> <mo>&amp;lsqb;</mo> <mrow> <msub> <mi>r</mi> <mi>b</mi> </msub> <msub> <mi>sin&amp;alpha;</mi> <mn>2</mn> </msub> <mo>+</mo> <mi>R</mi> <mi>R</mi> <mo>-</mo> <msqrt> <mrow> <msup> <mi>RR</mi> <mn>2</mn> </msup> <mo>-</mo> <msubsup> <mi>x</mi> <mn>1</mn> <mn>2</mn> </msubsup> </mrow> </msqrt> </mrow> <mo>&amp;rsqb;</mo> </mrow> <mn>3</mn> </msup> </mrow> </mfrac> <msub> <mi>dx</mi> <mn>1</mn> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <msubsup> <mo>&amp;Integral;</mo> <mrow> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> </msubsup> <mfrac> <mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <msup> <mi>cos&amp;alpha;sin</mi> <mn>2</mn> </msup> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> <mrow> <mn>2</mn> <mi>E</mi> <mi>L</mi> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mi>sin</mi> <mi>&amp;alpha;</mi> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <mi>cos</mi> <mi>&amp;alpha;</mi> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mrow> </mfrac> <mi>d</mi> <mi>&amp;alpha;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>-</mo> <msubsup> <mo>&amp;Integral;</mo> <mrow> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> </msubsup> <mfrac> <mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <msup> <mi>cos&amp;alpha;sin</mi> <mn>2</mn> </msup> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> <mrow> <mn>2</mn> <mi>E</mi> <mi>L</mi> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mi>sin</mi> <mi>&amp;alpha;</mi> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <mi>cos</mi> <mi>&amp;alpha;</mi> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>L</mi> </mrow> <mi>L</mi> </mfrac> <mi>d</mi> <mi>&amp;alpha;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mo>+</mo> <msubsup> <mo>&amp;Integral;</mo> <mrow> <mo>-</mo> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> </msubsup> <mfrac> <mrow> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <msup> <mi>cos&amp;alpha;sin</mi> <mn>2</mn> </msup> <msub> <mi>&amp;alpha;</mi> <mn>1</mn> </msub> </mrow> <mrow> <mn>2</mn> <mi>E</mi> <mi>L</mi> <mrow> <mo>&amp;lsqb;</mo> <mrow> <mi>sin</mi> <mi>&amp;alpha;</mi> <mo>+</mo> <mrow> <mo>(</mo> <mrow> <msub> <mi>&amp;alpha;</mi> <mn>2</mn> </msub> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> <mi>cos</mi> <mi>&amp;alpha;</mi> <mo>-</mo> <mfrac> <mi>H</mi> <mrow> <mn>2</mn> <msub> <mi>r</mi> <mi>b</mi> </msub> </mrow> </mfrac> </mrow> <mo>&amp;rsqb;</mo> </mrow> </mrow> </mfrac> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <mi>&amp;Delta;</mi> <mi>L</mi> </mrow> <mi>L</mi> </mfrac> <mi>d</mi> <mi>&amp;alpha;</mi> </mrow> </mtd> </mtr> </mtable> </mfenced>
式中，rb为基圆半径，rr为齿根圆半径，α1为啮合压力角，E为齿轮材料杨氏模量，RR为齿根过渡圆弧半径，x1为过渡部分对中见面投影到圆心的距离，H为点蚀坑深度，υ为齿轮材料泊松比，由此，便得到了齿轮系统内部激励，求得齿轮时变啮合刚度，完成基于概率分布的齿轮点蚀模型建模。
CN201710881432.1A 2017-09-26 2017-09-26 一种基于概率分布的齿轮点蚀模型建模方法 Active CN107729626B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710881432.1A CN107729626B (zh) 2017-09-26 2017-09-26 一种基于概率分布的齿轮点蚀模型建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710881432.1A CN107729626B (zh) 2017-09-26 2017-09-26 一种基于概率分布的齿轮点蚀模型建模方法

Publications (2)

Publication Number Publication Date
CN107729626A true CN107729626A (zh) 2018-02-23
CN107729626B CN107729626B (zh) 2020-03-17

Family

ID=61208049

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710881432.1A Active CN107729626B (zh) 2017-09-26 2017-09-26 一种基于概率分布的齿轮点蚀模型建模方法

Country Status (1)

Country Link
CN (1) CN107729626B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109002592A (zh) * 2018-06-26 2018-12-14 安徽理工大学 一种随机点蚀坑在板实体有限元模型中全自动建模方法
CN110567985A (zh) * 2019-10-14 2019-12-13 重庆大学 一种基于深度学习的自适应齿轮点蚀定量评估与检测装置
CN110657986A (zh) * 2019-10-11 2020-01-07 北京工业大学 齿轮双面啮合测量中测量力引入轮齿变形计算方法
CN112214845A (zh) * 2020-09-04 2021-01-12 燕山大学 一种基于新的点蚀模型的齿轮时变啮合刚度计算方法
CN112417738A (zh) * 2020-12-08 2021-02-26 江苏科技大学 一种含随机坑点蚀缺陷的球形耐压壳的数值计算方法
CN114997000A (zh) * 2022-05-24 2022-09-02 湖南大学 不同类型裂纹下的多级齿轮传动系统动态响应分析方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090301244A1 (en) * 2008-06-04 2009-12-10 Anderson George E gearbox apparatus and method of manufacture
CN104484489A (zh) * 2014-07-24 2015-04-01 江苏科技大学 一种点蚀损伤圆柱壳的四边形有限元网格自动生成方法
CN105224744A (zh) * 2015-09-29 2016-01-06 西安交通大学 一种基于啮合刚度的剥落齿轮啮合模型
CN105736676A (zh) * 2016-04-21 2016-07-06 合肥工业大学 一种基于蝴蝶鳞片表面形貌的仿生齿面齿轮
CN105736675A (zh) * 2016-04-21 2016-07-06 合肥工业大学 一种基于龟壳表面形貌的新型耐磨齿轮
CN107092750A (zh) * 2017-04-24 2017-08-25 北京航空航天大学 一种内啮合齿轮轴减速器轮齿损伤故障的非线性动力学建模方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090301244A1 (en) * 2008-06-04 2009-12-10 Anderson George E gearbox apparatus and method of manufacture
CN104484489A (zh) * 2014-07-24 2015-04-01 江苏科技大学 一种点蚀损伤圆柱壳的四边形有限元网格自动生成方法
CN105224744A (zh) * 2015-09-29 2016-01-06 西安交通大学 一种基于啮合刚度的剥落齿轮啮合模型
CN105736676A (zh) * 2016-04-21 2016-07-06 合肥工业大学 一种基于蝴蝶鳞片表面形貌的仿生齿面齿轮
CN105736675A (zh) * 2016-04-21 2016-07-06 合肥工业大学 一种基于龟壳表面形貌的新型耐磨齿轮
CN107092750A (zh) * 2017-04-24 2017-08-25 北京航空航天大学 一种内啮合齿轮轴减速器轮齿损伤故障的非线性动力学建模方法

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109002592A (zh) * 2018-06-26 2018-12-14 安徽理工大学 一种随机点蚀坑在板实体有限元模型中全自动建模方法
CN109002592B (zh) * 2018-06-26 2022-09-16 安徽理工大学 一种随机点蚀坑在板实体有限元模型中全自动建模方法
CN110657986A (zh) * 2019-10-11 2020-01-07 北京工业大学 齿轮双面啮合测量中测量力引入轮齿变形计算方法
CN110567985A (zh) * 2019-10-14 2019-12-13 重庆大学 一种基于深度学习的自适应齿轮点蚀定量评估与检测装置
CN112214845A (zh) * 2020-09-04 2021-01-12 燕山大学 一种基于新的点蚀模型的齿轮时变啮合刚度计算方法
CN112417738A (zh) * 2020-12-08 2021-02-26 江苏科技大学 一种含随机坑点蚀缺陷的球形耐压壳的数值计算方法
WO2022121200A1 (zh) * 2020-12-08 2022-06-16 江苏科技大学 一种含随机坑点蚀缺陷的球形耐压壳的数值计算方法
CN112417738B (zh) * 2020-12-08 2024-03-19 江苏科技大学 一种含随机坑点蚀缺陷的球形耐压壳的数值计算方法
CN114997000A (zh) * 2022-05-24 2022-09-02 湖南大学 不同类型裂纹下的多级齿轮传动系统动态响应分析方法
CN114997000B (zh) * 2022-05-24 2024-04-16 湖南大学 不同类型裂纹下的多级齿轮传动系统动态响应分析方法

Also Published As

Publication number Publication date
CN107729626B (zh) 2020-03-17

Similar Documents

Publication Publication Date Title
CN107729626A (zh) 一种基于概率分布的齿轮点蚀模型建模方法
CN101937211A (zh) 渐开线直齿圆锥齿轮修形方法
Tiwari et al. Stress analysis of mating involute spur gear teeth
Wei et al. Effects of flank deviation on load distributions for helical gear
CN103727204B (zh) 摆线马达摆线轮的修形方法
CN207830505U (zh) 一种rv减速器摆线轮
CN100414464C (zh) 渐开线直齿圆锥齿轮修形方法
CN110362858A (zh) 一种高压内啮合齿轮泵齿轮副的可靠性评估方法
Liu et al. Concave and convex modifications analysis for skewed beveloid gears considering misalignments
CN108052763A (zh) 一种基于遗传算法的摆线轮齿廓修形量优化方法
CN104408241A (zh) 一种修形圆柱齿轮的有限元网格自动生成方法
Dong et al. Rapid hob tip corner optimization of spur gears for increasing bending strength
CN113832267A (zh) 一种基于spc控制理论的高炉仪表数据智能监测系统和方法
CN111079300B (zh) 一种考虑齿向误差的直齿轮啮合刚度计算方法
Oh et al. A study on modeling and optimization of tooth microgeometry for a helical gear pair
CN206036178U (zh) 双点接触型斜齿轮
CN104565289B (zh) 大齿向修形量齿轮的设计方法及大齿向修形量齿轮
Zheng et al. Research on the mold release motion for spiral bevel gear forging
CN105156634B (zh) 渐开线直齿齿轮修形的优化方法
CN204784574U (zh) 一种驱动桥轮边行星减速器齿轮
CN105537479A (zh) 一种齿条冷摆辗精密成形方法
CN103234020A (zh) 一种消除渐开线齿轮根切的方法
CN208634285U (zh) 一种直齿锥齿轮副
CN109766562B (zh) 基于遗传算法和粒子群组合算法的摆线轮齿廓修形方法
CN1167881C (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