CN101770538A - 含损伤性单齿故障圆柱直齿轮啮合刚度仿真分析方法 - Google Patents
含损伤性单齿故障圆柱直齿轮啮合刚度仿真分析方法 Download PDFInfo
- Publication number
- CN101770538A CN101770538A CN201010034173A CN201010034173A CN101770538A CN 101770538 A CN101770538 A CN 101770538A CN 201010034173 A CN201010034173 A CN 201010034173A CN 201010034173 A CN201010034173 A CN 201010034173A CN 101770538 A CN101770538 A CN 101770538A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- alpha
- mfrac
- msubsup
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 87
- 238000004088 simulation Methods 0.000 title claims abstract description 40
- 238000004458 analytical method Methods 0.000 title claims abstract description 17
- 238000004364 calculation method Methods 0.000 claims abstract description 54
- 238000012937 correction Methods 0.000 claims abstract description 20
- 235000000621 Bidens tripartita Nutrition 0.000 claims abstract description 16
- 240000004082 Bidens tripartita Species 0.000 claims abstract description 16
- 208000006637 fused teeth Diseases 0.000 claims abstract description 16
- 238000010561 standard procedure Methods 0.000 claims abstract description 11
- 238000006073 displacement reaction Methods 0.000 claims description 18
- 239000011159 matrix material Substances 0.000 claims description 12
- 239000000463 material Substances 0.000 claims description 7
- 238000005452 bending Methods 0.000 claims description 3
- 230000006835 compression Effects 0.000 claims description 3
- 238000007906 compression Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 2
- 238000011160 research Methods 0.000 abstract description 5
- 230000007246 mechanism Effects 0.000 abstract description 4
- 230000004044 response Effects 0.000 abstract description 2
- 230000008685 targeting Effects 0.000 abstract 1
- 230000008859 change Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000009347 mechanical transmission Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Landscapes
- Gears, Cams (AREA)
Abstract
本发明涉及一种含损伤性单齿故障圆柱直齿轮啮合刚度仿真分析方法。该方法首先基于有限元法和国家标准方法的平均刚度计算结果,提出单、双齿啮合区间啮合刚度修正系数,用于改进能量法计算正常齿轮啮合刚度的精度。其次,针对齿轮故障部位,结合三维建模软件和有限元分析软件建立损伤性故障齿轮的有限元模型,采用计算机语言编制仿真计算程序,计算其时变啮合刚度;最后整合两部分的计算结果求解出故障齿轮的完整啮合刚度。该方法充分综合了修正能量法和有限元法的优势,既保证了计算精度,又提高了计算效率。应用本方法仿真计算的损伤性单齿故障齿轮啮合刚度可有效地用于齿轮系统的振动响应机理研究。
Description
技术领域
本发明属于齿轮测量技术领域,具体涉及一种故障齿轮啮合刚度仿真分析方法,特别是一种基于修正能量法和有限元法的含损伤性单齿故障的圆柱直齿轮时变啮合刚度的仿真分析方法。
背景技术
齿轮是机械传动设备中应用最广泛而且是最容易出现故障的零件之一。齿轮的啮合刚度对于齿轮传动的动力学性能有着显著影响。时变的啮合刚度是齿轮传动系统振动响应的主要动态激励源之一。轮齿变形与啮合刚度随啮合位置变化规律的研究是轮齿修形、动态特性、故障诊断以及寿命预测等研究的基础。所以,有必要深入地研究探讨故障齿轮的啮合刚度变化规律性及其快速有效地仿真分析方法。
能量法和有限元法是目前常用的齿轮啮合刚度计算方法。能量法基于弹性力学知识推导出刚度的积分公式,利用数值计算软件计算出时变啮合刚度的数值解,该方法计算效率较高但是精度和对于故障齿轮啮合刚度的计算不太理想。有限元法一般是利用有限元分析软件建立齿轮副接触有限元模型,通过非线性接触分析功能计算齿轮的时变啮合刚度,计算精度较高,但是计算量很大,效率不高。鉴于此,本发明利用国家标准方法的齿轮平均啮合刚度计算结果,先检验有限元法和能量法的计算误差,提出基于刚度修正系数的修正能量法,改进能量法的计算精度;再结合有限元法精确高效地计算损伤性故障齿轮的时变啮合刚度,本方法对研究损伤性故障圆柱直齿轮系统振动产生与扩展机理及有效的故障诊断技术有着非常重要的意义。
发明内容
本发明为了精确高效地求解含损伤性单齿故障齿轮副时变啮合刚度,提出了一种结合有限元法和修正能量法的齿轮啮合刚度仿真分析方法,采用此方法仿真计算损伤性单齿故障齿轮的时变啮合刚度达到了即保证精度,又提高计算效率的有益效果。
为实现上述目的,本发明的技术方案如下:
一种含损伤性故障圆柱直齿轮啮合刚度仿真分析方法,包括以下具体步骤:
步骤1正常齿轮啮合刚度仿真计算
步骤1.1用国家标准方法计算齿轮副平均啮合刚度
选定标准渐开线圆柱直齿轮的参数和材料特性:模数m,齿数z1、z2,齿轮宽度H,齿轮轴孔半径r,弹性模量E,泊松比ν,材料密度ρ;使用GB-T3480提供的计算公式:
c′=cthCMCRCB (1)
计算单对齿刚度,其中CM=0.8为理论修正系数,CR=1为轮柸结构系数,CB=1为基本齿廓系数,单对齿刚度的理论值:
再由公式:cγ=(0.75εα+0.25)H·c′×106 (3)
即可仿真计算出齿轮平均啮合刚度cγ,(3)式中εα为齿轮端面重合度,有:
步骤1.2用有限元法计算齿轮副时变啮合刚度
通过在齿轮副有限元模型的接触面上设定接触单元,将这些接触单元与齿轮副进行连接和组装,建立接触系统整体平衡方程为:
其中
K为总体刚度矩阵,u为节点位移矩阵,R为对角矩阵,N为单元形函数,L为算子,λ为节点处接触内力矩阵。
采用迭代法求解平衡方程,提取出轮齿轴孔节点的切向位移量uy,则齿轮副的啮合刚度:
式中N1为节点数,rb为小齿轮基圆半径,r为小齿轮轴孔半径,f为施加在每个轴孔节点上的切向力。求解(6)式即可求得一个啮合位置齿轮副的啮合刚度值。以小齿轮为主动轮,且顺时针转动,初始啮合位置定义为齿轮副刚进入双齿啮合区的临界位置。在齿轮副一个啮合周期内,对于双齿啮合区间和单齿啮合区间分别计算nys和nyd个啮合位置的刚度,第i个啮合位置对应的小齿轮转角位移为θi,啮合刚度为kθi 1,则一个啮合周期总的平均啮合刚度为:
为了验证有限元法的计算精度,做如下比较:
根据λ的取值不同,以小齿轮基于初始啮合位置的角位移θ为变量的有限元法最终啮合刚度仿真结果可表示为:
则双、单齿啮合区间的啮合刚度平均值cys和cyd可分别表示为:
步骤1.3用能量法计算齿轮副时变啮合刚度
齿轮的啮合刚度由赫兹刚度kh、弯曲刚度kb、径向压缩刚度ka和剪切刚度ks组成。分别可由下列公式求得:
其中,i=1,2分别对应于齿轮啮合时的第一对轮齿和第二对轮齿,α2、α′2分别是小、大齿轮的齿基半角;
invα0=tanα0-α (23)
以上各式中,θ为小齿轮的基于初始啮合位置的转角位移;则在小齿轮角位移为θ时齿轮副的总啮合刚度:
求解(24)式即可求得齿轮一个啮合位置的啮合刚度。在齿轮副一个啮合周期内,双齿啮合区间和单齿啮合区间分别计算nns和nnd个啮合位置的刚度,第i个啮合位置对应的小齿轮转角为θi,啮合刚度为kθi 2,则双、单齿啮合区间的啮合刚度平均值cns和cnd可分别表示为:
一个啮合周期总的平均啮合刚度为:
步骤2修正能量法的计算结果
基于步骤1中所得国家标准方法和有限元法的计算结果,定义并计算出用于修正能量法的单齿啮合区间啮合刚度修正系数μd和双齿啮合区间啮合刚度修正系数μs如下:
修正后的正常齿轮任意啮合位置啮合刚度为:
步骤3齿轮副故障部位的啮合刚度仿真计算
选定齿轮的损伤性单齿故障的类型和故障特性,建立故障齿轮副的有限元模型,之后采用步骤1中的有限元法仿真计算故障齿轮副在啮合位置θ的啮合刚度kθ 3;根据步骤1中λ的取值不同,故障部位的啮合刚度最终仿真结果可表示为:
因为损伤性单齿故障会影响齿轮副的两个啮合周期的啮合刚度,所以这个步骤中需要计算故障轮齿处于啮合时的两个连续地啮合周期的啮合刚度。
步骤4损伤性单齿故障齿轮完整啮合刚度仿真计算
以初始啮合的两对轮齿中左边的轮齿为基准,若故障轮齿是逆时针第a个轮齿;小齿轮旋转一周时,齿轮副有z1个啮合周期;其中在[1,a-1]和[a+2,z1]啮合周期的啮合刚度与无故障齿轮的啮合刚度相同,即采用步骤2中修正能量法所仿真计算的结果;在第a个和第a+1个啮合周期的啮合刚度则因故障轮齿处于啮合状态,所以齿轮副的啮合刚度采用步骤3中有限元法所仿真计算的结果;至此,即可得到含损伤性单齿故障齿轮的一个完整周期的啮合刚度。
本发明的有益效果是:利用国家标准方法的正常齿轮平均啮合刚度计算结果,检验有限元法和能量法的计算精度;并提出了针对能量法的刚度修正系数,改进了能量法的计算结果精度;再结合有限元法计算的故障部分齿轮啮合刚度,进而求解出损伤性故障齿轮完整的时变啮合刚度。该方法综合了修正能量法的高效性和有限元法的精确性,仿真分析的损伤性故障齿轮时变啮合刚度为齿轮系统故障的机理研究提供了准确、可靠的理论基础。
下面具体结合附图与实例对本发明作进一步的说明。
附图说明
图1是本发明的工作流程图;
图2是本发明正常齿轮副有限元模型;
图3是本发明仿真计算啮合刚度的APDL程序编制流程图;
图4是本发明有限元法仿真计算正常齿轮啮合刚度结果;
图5是本发明能量法仿真计算正常齿轮啮合刚度结果;
图6是本发明修正能量法和有限元法仿真计算正常齿轮啮合刚度结果;
图7是本发明有限元法仿真计算齿轮齿根裂纹故障部位啮合刚度结果;
图8是本发明齿根裂纹故障齿轮时变啮合刚度仿真结果;
具体实施方式
如图1所示,是本发明的一种含损伤性单齿故障圆柱直齿轮啮合刚度仿真分析方法的工作流程图。具体实施过程如下:
步骤1正常齿轮啮合刚度仿真计算
步骤1.1用国家标准方法计算齿轮副平均啮合刚度
选定标准渐开线圆柱直齿轮的参数和材料特性:模数m=3.175,齿数z1/z2=19/48,齿轮厚度H=16mm,齿轮轴孔半径r=20mm,弹性模量E=6.028Gpa,泊松比v=0.3,材料密度ρ=7800kg/m3。利用GB-T 3480提供的计算公式:
c′=cthCMCRCBcosβ
计算单对齿刚度,其中CM=0.8为理论修正系数,CR=1为轮柸结构系数,CB=1为基本齿廓系数,单对齿刚度的理论值:
再由公式:cγ=(0.75εα+0.25)H·c′×106 (3)
即可仿真计算出齿轮平均啮合刚度cγ,(3)式中εα为齿轮端面重合度,有:
本例中计算结果cγ=3.12×108N/m。
步骤1.2用有限元法计算齿轮副时变啮合刚度
利用三维建模软件SolidWorks建立上述步骤中选定齿轮的三维实体模型。将其保存为parasolid格式文件。在有限元分析软件Ansys中通过导入命令导入parasolid格式的齿轮副模型,经过材料属性定义,网格划分,设置接触对,定义约束,添加载荷后建立起齿轮副有限元仿真分析模型,最后设定分析类型,完成后的齿轮副有限元模型如图2。接触系统整体平衡方程为:
其中
K为总体刚度矩阵,u为节点位移矩阵,R为对角矩阵,N为单元形函数,L为算子,λ为节点处接触内力矩阵;
采用迭代法求解平衡方程,提取出轮齿轴孔节点的切向位移量uy,则齿轮副的啮合刚度:
式中Nl为节点数,rb为小齿轮基圆半径,r为小齿轮轴孔半径,f为施加在每个轴孔节点上的切向力;求解(6)式即可求得一个啮合位置齿轮副的啮合刚度值;以小齿轮为主动轮,且顺时针转动,初始啮合位置定义为齿轮副刚进入双齿啮合区的临界位置;在齿轮副一个啮合周期内,对于双齿啮合区间和单齿啮合区间分别计算nys和nyd个啮合位置的刚度,第i个啮合位置对应的小齿轮转角位移为θi,啮合刚度为kθi 1,则一个啮合周期总的平均啮合刚度为:
根据编制上述分析过程编制ansys的APDL程序,并加入循环结构和后处理计算分析命令,形成APDL计算分析程序,程序流程如图3所示。运行该程序可计算出cy=3.15×108N/m。
为了验证有限元法的计算精度,做如下比较:
所以以小齿轮基于初始啮合位置的角位移θ为变量的有限元法最终啮合刚度仿真结果可表示为:
有限元法仿真求解的正常齿轮时变啮合刚度,如图4。双、单齿啮合区间的啮合刚度平均值cys和cyd可分别表示为:
步骤1.3用能量法计算齿轮副时变啮合刚度
齿轮的啮合刚度由赫兹刚度kh、弯曲刚度kb、径向压缩刚度ka和剪切刚度ks组成;分别可由下列公式求得:
其中,i=1,2分别对应于齿轮啮合时的第一对轮齿和第二对轮齿,α2、α’2分别是小、大齿轮的齿基半角;
invα0=tanα0-α (23)
以上各式中,θ为小齿轮的基于初始啮合位置的转角位移;则在小齿轮角位移为θ时齿轮副的总啮合刚度:
求解(24)式即可求得齿轮一个啮合位置的啮合刚度;在齿轮副一个啮合周期内,双齿啮合区间和单齿啮合区间分别计算nns和nnd个啮合位置的刚度,第i个啮合位置对应的小齿轮转角为θi,啮合刚度为kθi 2,则双、单齿啮合区间的啮合刚度平均值cns和cnd可分别表示为:
一个啮合周期总的平均啮合刚度为:
本实例中计算结果cn=8.84×108N/m。
本实例中三种方法计算的正常齿轮平均啮合刚度结果比较如表1所示。从表中可以明显知道,有限元法的计算结果较能量法精确,其与标准方法的误差为0.9%。而能量法的计算结果偏大,但因能量法的结果也能反映时变啮合刚度的主要特征,而且在实际实验中发现其计算速度是有限元法的十倍左右。
(将下表从说明书附图移到了说明书中)好的。
表1正常齿轮啮合刚度c(108N/m)
方法 | 最大值 | 最小值 | 平均值 |
有限元法 | 3.66 | 2.38 | 3.15 |
能量法 | 10.55 | 5.78 | 8.84 |
国家标准方法 | ------ | ------ | 3.12 |
步骤2修正能量法的计算结果
基于步骤1中所得国家标准方法和有限元法的计算结果,定义并计算出用于修正能量法的单齿啮合区间啮合刚度修正系数μd和双齿啮合区间啮合刚度修正系数μs如下:
修正后的正常齿轮任意啮合位置啮合刚度为:
利用此修正能量法计算的正常齿轮啮合刚度和有限元法的计算结果如图6所示。从图中可以看出,两种方法的计算结果已经非常的接近,这充分地说明了修正能量法的有益效果:改进计算精度的同时,又保持了其计算效率。
步骤3齿轮副故障部位的啮合刚度仿真计算
本实例采用了实际情况中常见的齿根裂纹故障(其它损伤性故障亦适用),裂纹深度为3mm。利用SolidWorks建立故障齿轮高精度实体模型,导入ansys,步骤1中的有限元法,划分网格时采用混合分网的方式划分有限元网格:用切割命令将齿轮的正常部分和故障部分分离,正常部位用体扫掠划分网格,而有故障部位则用自由网格划分方法。这样即能快速简易地建立故障齿轮副的有限元模型,计算出故障齿轮副在啮合位置θ的啮合刚度kθ 3;因为步骤1中λ<1%,所以故障部位的啮合刚度最终仿真结果为:
因为损伤性单齿故障会影响齿轮副的两个啮合周期的啮合刚度,所以这个步骤中需要计算故障轮齿处于啮合时的两个连续地啮合周期的啮合刚度。本实例有限元法计算齿轮故障部位啮合刚度结果如图7所示。
步骤4损伤性单齿故障齿轮完整啮合刚度仿真计算
以初始啮合的两对轮齿中左边的轮齿为基准,故障轮齿是逆时针第2个轮齿;小齿轮旋转一周时,齿轮副有19个啮合周期;其中在1和[4,19]啮合周期的啮合刚度与无故障齿轮的啮合刚度相同,即采用步骤2中修正能量法所仿真计算的结果;在第2个和第3个啮合周期的啮合刚度则因故障轮齿处于啮合状态,所以齿轮副的啮合刚度采用步骤3中有限元法所仿真计算的结果;至此,即可得到含损伤性单齿故障齿轮的一个完整周期的啮合刚度,如图8。
本发明修正了能量法计算齿轮正常部位啮合刚度的结果精度,用有限元法计算齿轮故障部位啮合刚度的结果,将两者整合成故障齿轮的完整啮合刚度,达到了高效精确的目的。仿真分析的损伤性单齿故障齿轮时变啮合刚度为齿轮系统故障的机理研究提供了准确、可靠的理论基础。
Claims (1)
1.一种含损伤性单齿故障圆柱直齿轮啮合刚度仿真分析方法,其特征在于,包括以下步骤:
步骤1正常齿轮啮合刚度计算
步骤1.1用国家标准方法计算齿轮副平均啮合刚度
选定标准渐开线圆柱直齿轮的参数和材料特性:模数m,齿数z1、z2,齿轮宽度H,齿轮轴孔半径r,弹性模量E,泊松比v,材料密度ρ;使用GB-T3480提供的计算公式:
c′=cthCMCRCB (1)
计算单对齿刚度,其中CM=0.8为理论修正系数,CR=1为轮柸结构系数,CB=1为基本齿廓系数,单对齿刚度的理论值:
再由公式:cγ=(0.75εα+0.25)H·c′×106 (3)
即可计算出齿轮平均啮合刚度cγ,(3)式中εα为齿轮端面重合度,有:
步骤1.2用有限元法计算齿轮副时变啮合刚度
通过在齿轮副有限元模型的接触面上设定接触单元,将这些接触单元与齿轮副进行连接和组装,建立接触系统整体平衡方程为:
其中
K为总体刚度矩阵,u为节点位移矩阵,R为对角矩阵,N为单元形函数,L为算子,λ为节点处接触内力矩阵;
采用迭代法求解平衡方程,提取出轮齿轴孔节点的切向位移量uy,则齿轮副的啮合刚度:
式中N1为节点数,rb为小齿轮基圆半径,f为施加在每个轴孔节点上的切向力;
以小齿轮为主动轮,且顺时针转动,初始啮合位置定义为齿轮副刚进入双齿啮合区的临界位置;在齿轮副一个啮合周期内,对于双齿啮合区间和单齿啮合区间分别计算nys和nyd个啮合位置的刚度,第i个啮合位置对应的小齿轮转角位移为θi,啮合刚度为kθi 1;一个啮合周期总的平均啮合刚度为:
为了校准有限元法的计算精度,做如下比较:
根据λ的取值不同,以小齿轮基于初始啮合位置的角位移θ为变量的有限元法最终啮合刚度仿真结果可表示为:
则双、单齿啮合区间的啮合刚度平均值cys和cyd可分别表示为:
步骤1.3用能量法计算齿轮副时变啮合刚度
齿轮的啮合刚度由赫兹刚度kh、弯曲刚度kb、径向压缩刚度ka和剪切刚度ks组成;分别可由下列公式求得:
其中,i=1,2分别对应于齿轮啮合时的第一对轮齿和第二对轮齿,α2、α′2分别是小、大齿轮的齿基半角;
invα0=tanα0-α (23)
以上各式中,θ为小齿轮的基于初始啮合位置的转角位移;则在小齿轮角位移为θ时齿轮副的总啮合刚度:
在齿轮副一个啮合周期内,双齿啮合区间和单齿啮合区间分别计算nns和nnd个啮合位置的刚度,第i个啮合位置对应的小齿轮转角为θi,啮合刚度为kθi 2,则双、单齿啮合区间的啮合刚度平均值cns和cnd可分别表示为:
一个啮合周期总的平均啮合刚度为:
步骤2修正能量法的计算结果
基于步骤1中所得国家标准方法和有限元法的计算结果,定义并计算出用于修正能量法的单齿啮合区间啮合刚度修正系数μd和双齿啮合区间啮合刚度修正系数μs如下:
修正后的正常齿轮任意啮合位置啮合刚度为:
步骤3齿轮副故障部位的啮合刚度仿真计算
选定齿轮的损伤性单齿故障的类型和故障特性,建立故障齿轮副的有限元模型,之后采用步骤1中的有限元法仿真计算故障齿轮副在啮合位置θ的啮合刚度kθ 3;根据步骤1中λ的取值不同,故障部位的啮合刚度最终仿真结果可表示为:
因为损伤性单齿故障会影响齿轮副的两个啮合周期的啮合刚度,所以这个步骤中需要计算故障轮齿处于啮合时的两个连续地啮合周期的啮合刚度;
步骤4损伤性单齿故障齿轮完整啮合刚度仿真计算
以初始啮合的两对轮齿中左边的轮齿为基准,若故障轮齿是逆时针第a个轮齿;小齿轮旋转一周时,齿轮副有z1个啮合周期;其中在[1,a-1]和[a+2,z1]啮合周期的啮合刚度与无故障齿轮的啮合刚度相同,即采用步骤2中修正能量法所仿真计算的结果;在第a个和第a+1个啮合周期的啮合刚度则因故障轮齿处于啮合状态,所以齿轮副的啮合刚度采用步骤3中有限元法所仿真计算的结果;至此,即可得到含损伤性单齿故障齿轮的一个完整周期的啮合刚度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010100341737A CN101770538B (zh) | 2010-01-15 | 2010-01-15 | 含损伤性单齿故障圆柱直齿轮啮合刚度仿真分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010100341737A CN101770538B (zh) | 2010-01-15 | 2010-01-15 | 含损伤性单齿故障圆柱直齿轮啮合刚度仿真分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101770538A true CN101770538A (zh) | 2010-07-07 |
CN101770538B CN101770538B (zh) | 2011-11-09 |
Family
ID=42503394
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010100341737A Expired - Fee Related CN101770538B (zh) | 2010-01-15 | 2010-01-15 | 含损伤性单齿故障圆柱直齿轮啮合刚度仿真分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101770538B (zh) |
Cited By (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102235465A (zh) * | 2011-06-30 | 2011-11-09 | 重庆大学 | 行星齿轮机构运动学分析方法及分析系统 |
CN102262697A (zh) * | 2011-07-20 | 2011-11-30 | 上海师范大学 | 一种用于斜齿圆锥齿轮的建模方法 |
WO2013000261A1 (zh) * | 2011-06-30 | 2013-01-03 | 重庆大学 | 行星齿轮机构机械效率分析方法及系统 |
CN103234021A (zh) * | 2013-04-23 | 2013-08-07 | 北京工业大学 | 一种确定克林跟贝尔格锥齿轮时变啮合刚度的方法 |
CN103577687A (zh) * | 2013-09-23 | 2014-02-12 | 北京工业大学 | 一种局部缺陷齿轮啮合刚度的时变特性定量计算方法 |
CN104091007A (zh) * | 2014-07-01 | 2014-10-08 | 重庆工商职业学院 | 一种少齿差行星减速器动力学仿真分析方法 |
CN104502095A (zh) * | 2015-01-05 | 2015-04-08 | 盐城工学院 | 测量直齿轮啮合阻尼及其阻尼成分的方法 |
CN104535318A (zh) * | 2014-12-29 | 2015-04-22 | 盐城工学院 | 测量齿轮啮合时变刚度的方法 |
CN105181327A (zh) * | 2015-08-26 | 2015-12-23 | 北京工业大学 | 一种裂纹轮齿啮合刚度计算的方法 |
CN105784360A (zh) * | 2016-05-12 | 2016-07-20 | 重庆长安汽车股份有限公司 | 一种基于啮合接触线长度变化确定齿轮啮合动刚度的方法 |
CN107420523A (zh) * | 2017-09-14 | 2017-12-01 | 东北大学 | 一种具有齿面裂纹缺陷的斜齿轮副啮合刚度计算方法 |
CN107436982A (zh) * | 2017-07-27 | 2017-12-05 | 东北大学 | 考虑基体刚度修正的剥落斜齿轮副的啮合特性分析方法 |
CN108534966A (zh) * | 2017-03-02 | 2018-09-14 | 武汉理工大学 | 一种齿轮时变啮合刚度测量计算方法 |
CN109063300A (zh) * | 2018-07-24 | 2018-12-21 | 北京工业大学 | 一种基于改进能量法的行星齿轮时变啮合刚度求解方法 |
CN110657986A (zh) * | 2019-10-11 | 2020-01-07 | 北京工业大学 | 齿轮双面啮合测量中测量力引入轮齿变形计算方法 |
CN114354187A (zh) * | 2022-01-05 | 2022-04-15 | 上海交通大学 | 基于辨识啮合刚度的齿轮故障分类检测方法及系统 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103971006B (zh) * | 2014-05-16 | 2017-06-13 | 清华大学 | 一种考虑主减速器壳的驱动桥齿轮动力学特性确定方法 |
-
2010
- 2010-01-15 CN CN2010100341737A patent/CN101770538B/zh not_active Expired - Fee Related
Cited By (22)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013000261A1 (zh) * | 2011-06-30 | 2013-01-03 | 重庆大学 | 行星齿轮机构机械效率分析方法及系统 |
CN102235465B (zh) * | 2011-06-30 | 2014-10-08 | 重庆大学 | 行星齿轮机构运动学分析方法及分析系统 |
CN102235465A (zh) * | 2011-06-30 | 2011-11-09 | 重庆大学 | 行星齿轮机构运动学分析方法及分析系统 |
CN102262697A (zh) * | 2011-07-20 | 2011-11-30 | 上海师范大学 | 一种用于斜齿圆锥齿轮的建模方法 |
CN103234021A (zh) * | 2013-04-23 | 2013-08-07 | 北京工业大学 | 一种确定克林跟贝尔格锥齿轮时变啮合刚度的方法 |
CN103234021B (zh) * | 2013-04-23 | 2015-09-02 | 北京工业大学 | 一种确定克林跟贝尔格锥齿轮时变啮合刚度的方法 |
CN103577687A (zh) * | 2013-09-23 | 2014-02-12 | 北京工业大学 | 一种局部缺陷齿轮啮合刚度的时变特性定量计算方法 |
CN103577687B (zh) * | 2013-09-23 | 2017-05-03 | 北京工业大学 | 一种局部缺陷齿轮啮合刚度的时变特性定量计算方法 |
CN104091007A (zh) * | 2014-07-01 | 2014-10-08 | 重庆工商职业学院 | 一种少齿差行星减速器动力学仿真分析方法 |
CN104535318B (zh) * | 2014-12-29 | 2017-02-22 | 盐城工学院 | 测量齿轮啮合时变刚度的方法 |
CN104535318A (zh) * | 2014-12-29 | 2015-04-22 | 盐城工学院 | 测量齿轮啮合时变刚度的方法 |
CN104502095A (zh) * | 2015-01-05 | 2015-04-08 | 盐城工学院 | 测量直齿轮啮合阻尼及其阻尼成分的方法 |
CN105181327A (zh) * | 2015-08-26 | 2015-12-23 | 北京工业大学 | 一种裂纹轮齿啮合刚度计算的方法 |
CN105181327B (zh) * | 2015-08-26 | 2018-05-08 | 北京工业大学 | 一种裂纹轮齿啮合刚度计算的方法 |
CN105784360A (zh) * | 2016-05-12 | 2016-07-20 | 重庆长安汽车股份有限公司 | 一种基于啮合接触线长度变化确定齿轮啮合动刚度的方法 |
CN108534966A (zh) * | 2017-03-02 | 2018-09-14 | 武汉理工大学 | 一种齿轮时变啮合刚度测量计算方法 |
CN107436982A (zh) * | 2017-07-27 | 2017-12-05 | 东北大学 | 考虑基体刚度修正的剥落斜齿轮副的啮合特性分析方法 |
CN107436982B (zh) * | 2017-07-27 | 2020-04-14 | 东北大学 | 考虑基体刚度修正的剥落斜齿轮副的啮合特性分析方法 |
CN107420523A (zh) * | 2017-09-14 | 2017-12-01 | 东北大学 | 一种具有齿面裂纹缺陷的斜齿轮副啮合刚度计算方法 |
CN109063300A (zh) * | 2018-07-24 | 2018-12-21 | 北京工业大学 | 一种基于改进能量法的行星齿轮时变啮合刚度求解方法 |
CN110657986A (zh) * | 2019-10-11 | 2020-01-07 | 北京工业大学 | 齿轮双面啮合测量中测量力引入轮齿变形计算方法 |
CN114354187A (zh) * | 2022-01-05 | 2022-04-15 | 上海交通大学 | 基于辨识啮合刚度的齿轮故障分类检测方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN101770538B (zh) | 2011-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101770538B (zh) | 含损伤性单齿故障圆柱直齿轮啮合刚度仿真分析方法 | |
CN104573196B (zh) | 一种斜齿圆柱齿轮时变啮合刚度解析计算方法 | |
Wan et al. | Mesh stiffness calculation using an accumulated integral potential energy method and dynamic analysis of helical gears | |
CN104820756B (zh) | 一种考虑延长啮合的裂纹齿轮转子系统动力参数确定方法 | |
Ma et al. | Time-varying mesh stiffness calculation of spur gears with spalling defect | |
CN109726520B (zh) | 考虑复杂基体与裂纹扩展路径的直齿轮啮合刚度计算方法 | |
CN106844818A (zh) | 基于粗糙表面的直齿轮三维接触刚度计算方法 | |
CN107944174B (zh) | 一种圆柱齿轮齿向载荷分布系数获取方法 | |
Wang et al. | Simulating coupling behavior of spur gear meshing and fatigue crack propagation in tooth root | |
CN103577687A (zh) | 一种局部缺陷齿轮啮合刚度的时变特性定量计算方法 | |
CN104408220A (zh) | 一种改进的轮齿加载接触分析方法 | |
Wang et al. | A model to determine mesh characteristics in a gear pair with tooth profile error | |
CN109190227A (zh) | 基于解析-有限元的复杂齿基直齿圆柱齿轮啮合刚度计算方法 | |
CN102706306A (zh) | 渐开线齿轮波度样板 | |
CN105181327A (zh) | 一种裂纹轮齿啮合刚度计算的方法 | |
Sun et al. | Mesh stiffness and dynamic response analysis of modified gear system with thin web and weight reduction holes | |
CN105787149A (zh) | 一种由弧齿锥齿轮传动系统轴上功率谱向齿面应力谱精确转换的方法 | |
CN106295015B (zh) | 一种渐开线直齿圆柱齿轮副的齿廓修形方法及与其配套的专用参数化cad系统 | |
CN104535318A (zh) | 测量齿轮啮合时变刚度的方法 | |
Li et al. | Tooth profile modification and simulation analysis of involute spur gear | |
Li et al. | A Loaded Analysis Method for RV Cycloidal-pin Transmission Based on the Minimum Energy Principle. | |
CN113868755B (zh) | 基于基节误差的斜齿轮啮合刚度计算方法 | |
CN112214845A (zh) | 一种基于新的点蚀模型的齿轮时变啮合刚度计算方法 | |
Li et al. | STATICS ANALYSIS OF CYLINDRICAL GEAR DRIVE AT ANY MESHING POSITION BASED ON ANSYS. | |
Wu et al. | A revised method to calculate time-varying mesh stiffness of helical gear |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20111109 Termination date: 20130115 |
|
CF01 | Termination of patent right due to non-payment of annual fee |