CN107451376B - 一种偏航变桨轴承滚道校核工具 - Google Patents

一种偏航变桨轴承滚道校核工具 Download PDF

Info

Publication number
CN107451376B
CN107451376B CN201710753617.4A CN201710753617A CN107451376B CN 107451376 B CN107451376 B CN 107451376B CN 201710753617 A CN201710753617 A CN 201710753617A CN 107451376 B CN107451376 B CN 107451376B
Authority
CN
China
Prior art keywords
load
bearing
rolling element
row
contact
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.)
Active
Application number
CN201710753617.4A
Other languages
English (en)
Other versions
CN107451376A (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.)
Shanwei Mingyang New Energy Technology Co.,Ltd.
Original Assignee
MingYang Smart Energy Group Co Ltd
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 MingYang Smart Energy Group Co Ltd filed Critical MingYang Smart Energy Group Co Ltd
Priority to CN201710753617.4A priority Critical patent/CN107451376B/zh
Publication of CN107451376A publication Critical patent/CN107451376A/zh
Application granted granted Critical
Publication of CN107451376B publication Critical patent/CN107451376B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G06Q10/06393Score-carding, benchmarking or key performance indicator [KPI] analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/04Ageing analysis or optimisation against ageing

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Geometry (AREA)
  • Development Economics (AREA)
  • Marketing (AREA)
  • Health & Medical Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • Entrepreneurship & Innovation (AREA)
  • General Business, Economics & Management (AREA)
  • Educational Administration (AREA)
  • General Health & Medical Sciences (AREA)
  • Computational Mathematics (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • Game Theory and Decision Science (AREA)
  • Primary Health Care (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Rolling Contact Bearings (AREA)
  • Coloring Foods And Improving Nutritive Qualities (AREA)

Abstract

本发明公开了一种偏航变桨轴承滚道校核工具,直接利用LDD载荷谱进行疲劳寿命评估,对时序载荷进等效处理,形式与LDD载荷谱形式一致;工具可根据经验公式或静力学方程法求解轴承静强度,也可通过积分公式法或静力学方程法求解轴承疲劳寿命,同时可求解最大动态剪切应力及硬化层深度。工具可以根据单排轴承的运行转数、最大受载滚道的运行转数、单排轴承的运行时间或最大受载滚道的运行时间来评估轴承寿命;工具可以根据双排轴承的运行转数、单排轴承的运行转数、最大受载滚道的运行转数、双排轴承的运行时间、单排轴承的运行时间或最大受载滚道的运行时间评估寿命。本发明工具可非常有效缩短整机开始设计周期,提高风机可靠性。

Description

一种偏航变桨轴承滚道校核工具
技术领域
本发明涉及风力发电机组部件强度分析领域,尤其是指一种偏航变桨轴承滚道校核工具。
背景技术
业内习知,偏航变桨轴承作为风力发电机组的回转支承部件,其结构的可靠性决定风机运行的情况,偏航变桨轴承的结构形式一般为单排四点接触球轴承、双排四点接触球轴承,目前每一家轴承供应商的计算方法、评价指标不一,这就为后期的载荷评估带来麻烦,因此有必要针对偏航变桨轴承建立一个集成不同算法、不同评价标准的工具。
发明内容
本发明的目的在于克服现有偏航变桨轴承滚道校核计算过程中的不足,提供了一种高效、快捷、系统的偏航变桨轴承滚道校核工具,该工具集成了不同的计算方法及评价指标,能够满足不同供应商计算轴承的要求,可非常有效缩短整机开始设计周期,更好服务于风力发电机设计、认证和评估,提高风机发电机组设计的可靠性,节约成本。
为实现上述目的,本发明所提供的技术方案为:一种偏航变桨轴承滚道校核工具,所述工具为基于Matlab开发的偏航变桨轴承滚道校核软件,主要用于风力发电机偏航、变桨轴承初期选型、认证以及后期的结构强度评估;
所述工具直接利用LDD载荷谱进行疲劳寿命评估,同时能够对时序载荷进等效处理,将每一个工况等效为一组载荷,等效的载荷谱的形式与LDD载荷谱形式一致;
偏航变桨轴承通常采用单排四点接触球轴承或双排四点接触球轴承的结构形式,所述工具能够根据经验公式或静力学方程法求解轴承静强度,也能够通过积分公式法或静力学方程法求解轴承疲劳寿命,同时求解最大动态剪切应力及硬化层深度;
所述工具针对不同轴承供应商的轴承寿命提出以下寿命评估方法:
对于单排四点接触球轴承,根据单排轴承的运行转数、最大受载滚道的运行转数、单排轴承的运行时间或最大受载滚道的运行时间来评估轴承寿命;
对于双排四点接触球轴承,根据双排轴承的运行转数、单排轴承的运行转数、最大受载滚道的运行转数、双排轴承的运行时间、单排轴承的运行时间或最大受载滚道的运行时间评估寿命;
其中,按转数进行评估轴承寿命时应以载荷谱中实际的转数为评价指标;按时间进行评估轴承寿命时陆上风电应以20年的运行时间为评价指标,海上风电应以25年的运行时间为评价指标。
所述时序载荷的处理情况如下:
时序载荷:
对每一个时序载荷工况的力、弯矩、速度通过下面方程1.1-1.3处理得到等效载荷工况表;
Figure GDA0002424747580000021
Figure GDA0002424747580000022
Figure GDA0002424747580000023
式1.1-1.3中,F为轴向力、径向力等效载荷,Fi为时序载荷中第i步的力,M为合弯矩等效载荷,Mi为时序载荷中第i步的弯矩,n为等效转速,ni为时序载荷中第i步的转速;
将每一个疲劳工况的发生时间按下述方法求解:
对于发生次数的工况:
T20/25=L×T1×0.05×b (1.4)
对于发生时间的工况:
T20/25=6×L×T1×0.05×b (1.5)
式1.4-1.5中,L为风机的寿命要求,对于海上工况为25年,陆上工况20年,T1为每个工况发生的次数或小时数,b为每个工况发生的步数,0.05表示步长,T20/25为20年或25年每个疲劳工况发生的时间;
对于LDD载荷则可直接导入工具进行寿命计算;
所述静力学方程的情况如下:
1)四点接触球轴承的受载后几何形变如下:
坐标系采用GL规范中叶根坐标系,假设轴承周边结构具有足够的刚性,外滚道沟渠率中心位置坐标系为xo,yo,内滚道沟渠率中心位置坐标系为xi,yi,轴承外圈固定,内圈相对于外圈产生轴向位移δa,径向位移δr,扭角θ;坐标系位置也由xi,yi转化xi',yi'再到xi”,yi”最后到xi”',yi”',沟渠率中心距离也由最初的MN变为MN”';
对于双排滚子轴承,由排距d与转角θ对内外滚道产生0.5dθ变形;
对于单排轴承,不存在排距的影响;
2)双排四点接触球轴承的力学模型与负载分布计算
四个接触对在位置角ψ处的接触力分别表示为Q、Q、Q、Q,根据Hertz点接触理论,法向接触载荷Q和接触变形δ的关系如下:
Figure GDA0002424747580000031
式中,在位置角ψj处,内圈受到轴向载荷Fa、径向载荷Fr、倾覆力矩M以及钢球对内滚道的接触载荷Q的作用;Kn为滚动体与内外圈总的负荷变形常数。
以内圈与滚动体为研究对象,内圈在外部载荷和所有滚动体载荷的作用下处于平衡状态,内圈的力学平衡方程为:
Figure GDA0002424747580000032
Figure GDA0002424747580000033
Figure GDA0002424747580000034
其中,α、α、α、α分别为受载后四个接触对在位置角ψ处的接触角;
上述方程通过Newton Raphson方法进行迭代求解,得到滚道精确的载荷分布;
3)单排四点接触球轴承的力学模型与负载分布计算
单排四点接触球轴承的静力学方程求解原理与上述双排四点接触球轴承的计算原理相同;
两个接触对在位置角ψ处的接触力分别表示为Q、Q,受载后两个接触对在位置角ψ处的接触角为α、α
故平衡方程为:
Figure GDA0002424747580000041
Figure GDA0002424747580000042
Figure GDA0002424747580000043
所述静强度计算的情况具体如下:
1)积分公式法的静强度根据下列公式求解
单排轴承:
Figure GDA0002424747580000044
双排轴承:
Figure GDA0002424747580000045
根据赫兹接触理论得:
Figure GDA0002424747580000046
Figure GDA0002424747580000047
式3.1-3.4中,Qmax为滚动体最大载荷,Fr为径向载荷,Fa为轴向载荷,M为合弯矩,Z为钢球总数,α为初始接触角,dm为节圆直径,a为接触椭圆短半轴,b为接触椭圆的长半轴,Smin为最小安全系数,[σ]为许用接触应力,σmax为最大接触应力;
2)静力学方程的静强度求解
根据方程2.2-2.7,通过Newton Raphson方法即可求出最大滚动体的载荷Qmax,其余的计算过程与积分公式法一致;
所述疲劳寿命的计算情况如下:
1)滚道的基本额定滚动体负荷:
Figure GDA0002424747580000051
其中,内滚道取上符号,外滚道取下符号;
Figure GDA0002424747580000052
Figure GDA0002424747580000053
式4.1-4.3中,Qci为内滚道的基本额定动负荷,Qce为外滚道的基本额定动负荷,fi为内滚道沟曲率半径系数,fo为外滚道沟曲率半径系数,γ为无量纲几何参数,为Dw为滚子直径,dm为滚道中心径,Z为钢球总个数;
2)积分公式求解单排四点接触球轴承当量滚动体负荷及接触滚道寿命
2.1)当量滚动体负荷
相对负荷方向旋转套圈的当量滚动体负荷积分为:
Figure GDA0002424747580000054
相对负荷方向静止套圈的当量滚动体负荷积分为:
Figure GDA0002424747580000055
相对负荷方向旋转套圈的当量滚动体负荷为:
Qev1=Qev2=QmaxJ1(ε) (4.6)
相对负荷方向静止套圈的当量滚动体负荷为:
Qeu1=Qeu2=QmaxJ2(ε) (4.7)
式4.4-4.7中,ε为载荷系数,ψ为位置角,Qmax为滚动体最大载荷,根据式3.1计算;
2.2)接触滚道的寿命,按百万转计算
相对负荷旋转套圈主推力滚道的额定寿命L10ev1为:
Figure GDA0002424747580000056
相对负荷静止套圈主推力滚道的额定寿命L10eu1为:
Figure GDA0002424747580000061
相对负荷旋转套圈主推力滚道的额定寿命L10ev2为:
Figure GDA0002424747580000062
相对负荷静止套圈副推力滚道的额定寿命L10eu2为:
Figure GDA0002424747580000063
3)积分公式求解双排四点接触球轴承当量滚动体负荷及接触滚道寿命
3.1)当量滚动体负荷
对于双排四点接触球轴承当量滚动体负荷积分如下:
相对负荷方向旋转套圈的当量滚动体负荷积分为:
Figure GDA0002424747580000064
相对负荷方向静止套圈的当量滚动体负荷积分为:
Figure GDA0002424747580000065
主推力滚道最大滚动体负荷为:
Figure GDA0002424747580000066
副推力滚道最大滚动体负荷为:
Figure GDA0002424747580000067
相对负荷方向旋转套圈的当量滚动体负荷为:
Qev13=Qmax1J1 (4.16)
Qev24=Qmax2J1 (4.17)
相对负荷方向静止套圈的当量滚动体负荷为:
Qeu13=Qmax1J2 (4.18)
Qeu24=Qmax2J2 (4.19)
3.2)接触滚道的寿命,按百万转计算
相对负荷旋转套圈主推力滚道的额定寿命L10ev13为:
Figure GDA0002424747580000071
相对负荷静止套圈主推力滚道的额定寿命L10eu13为:
Figure GDA0002424747580000072
相对负荷旋转套圈主推力滚道的额定寿命L10ev24为:
Figure GDA0002424747580000073
相对负荷静止套圈副推力滚道的额定寿命L10eu24为:
Figure GDA0002424747580000074
4)静力学方程求解单排四点接触球轴承当量滚动体负荷及接触滚道寿命
4.1)当量滚动体负荷
相对旋转套圈主推力滚道的当量滚动体负荷Qev1为:
Figure GDA0002424747580000075
相对静止套圈主推力滚道的当量滚动体负荷Qeu1为:
Figure GDA0002424747580000076
相对旋转套圈反推力滚道的当量滚动体负荷Qev2为:
Figure GDA0002424747580000081
相对静止套圈反推力滚道的当量滚动体负荷Qeu2为:
Figure GDA0002424747580000082
其中,Q是作用于主推力滚道各个滚动体上的负荷,Q是作用于反推力滚道各个滚动体上的负荷;
4.2)接触滚道的寿命,按百万转计算
接触滚道的寿命由方程4.8-4.11求解;
5)静力学方程求解双排四点接触球轴承当量滚动体负荷及接触滚道寿命
5.1)当量滚动体负荷
主推力滚道(接触对1、接触对3)的当量滚动体负荷为:
Figure GDA0002424747580000083
Figure GDA0002424747580000084
Figure GDA0002424747580000085
Figure GDA0002424747580000086
反推力滚道的当量滚动体负荷为:
Figure GDA0002424747580000087
Figure GDA0002424747580000088
Figure GDA0002424747580000091
Figure GDA0002424747580000092
其中,Q、Q是作用于主推力滚道各个滚动体上的负荷,Q、Q是作用于反推力滚道各个滚动体上的负荷,Qev1、Qev2、Qev3、Qev4分别为旋转套圈对应的四个接触位置的当量滚动体负荷,Qeu1、Qeu2、Qeu3、Qeu4分别为静止套圈对应的四个接触位置的当量滚动体负荷。
5.2)接触滚道的寿命,按百万转计算
相对负荷旋转套圈主推力滚道(接触对1、接触对3)的额定寿命分别为:
Figure GDA0002424747580000093
Figure GDA0002424747580000094
相对负荷静止套圈主推力滚道(接触对1、接触对3)的额定寿命分别为:
Figure GDA0002424747580000095
Figure GDA0002424747580000096
相对负荷旋转套圈副推力滚道(接触对2、接触对4)的额定寿命分别为:
Figure GDA0002424747580000097
Figure GDA0002424747580000098
相对负荷静止套圈副推力滚道(接触对2、接触对4)的额定寿命分别为:
Figure GDA0002424747580000099
Figure GDA00024247475800000910
所述寿命评估标准的情况如下:
1)积分公式求解单排四点接触球轴承寿命
整个轴承的运行转数为:
L10=[L10ev1 (-10/9)+L10eu1 (-10/9)+L10ev2 (-10/9)+L10eu2 (-10/9)]-0.9×106 (6.1)
最大受载滚道的运行转数为:
L10=[L10ev1 (-10/9)+L10eu1 (-10/9)]-0.9×106 (6.2)
2)积分公式求解双排四点接触球轴承寿命
整个轴承的运行转数为:
L10=[L10ev13 (-10/9)+L10eu13 (-10/9)+L10ev24 (-10/9)+L10eu24 (-10/9)]-0.9×106 (6.3)
单个轴承最大受载下的运行转数为:
L10=[L10ev13 (-10/9)+L10eu13 (-10/9)]-0.9×106 (6.4)
最大受载滚道的运行转数
L10=L10ev13×106 (6.5)
3)静力学方程求解单排四点接触球轴承寿命,其方法与式子6.1-6.2相同;
4)静力学方程求解双排四点接触球轴承寿命
整个轴承的运行转数为:
Figure GDA0002424747580000101
单个轴承最大受载下的运行转数为:
L10=[L10ev1 (-10/9)+L10eu1 (-10/9)+L10ev2 (-10/9)+L10eu2 (-10/9)]-0.9×106 (6.7)
最大受载滚道的运行转数为:
L10=[L10ei1 (-10/9)+L10ee1 (-10/9)]-0.9×106 (6.8)
5)按时间
根据运行时间评估轴承寿命,需要在上述方程的基础上按下式进行求解
Figure GDA0002424747580000102
其中,n为等效转速;
所述最大动态剪切应力及硬化层深度的情况如下:
Figure GDA0002424747580000111
Figure GDA0002424747580000112
Figure GDA0002424747580000113
式7.1-7.3中,τ0为最大动态剪切应力,Q为接触力,t为滚道与钢球的接触尺寸的结构参数,z0为硬化层深度。
所述工具的参数输入与结果输出的情况如下:
工具界面分为五部分:轴承输入参数、计算方法、载荷导入、计算结果以及结果输出;
在轴承输入参数中需要输入钢球直径、轴承节圆直径、内滚道沟曲率半径系数、外滚道沟曲率半径系数、初始接触角、单排钢球数目、列数、双排滚子间距、轴承游隙、钢球弹性模量、内套圈弹性模量、外套圈弹性模量、内套圈泊松比、外套圈泊松比、滚道硬度、许用接触应力;
在计算方法中选择积分公式或Newton Raphson方法;
在载荷导入中导入包含极限载荷与LDD载荷谱的文件,在载荷输入方法中选择LDD或时序载荷TS,当选择LDD时,无需选择时序路径,当选择TS时需要在载荷导入中选择时序的文件夹;
在计算结果中显示最终计算结果,包括钢球受到最大接触力、钢球受到最大接触应力、最小静安全系数、最大动态剪切应力、硬化层深度、计算寿命;
在结果输出中显示详细的计算结果报告的存放路径。
本发明与现有技术相比,具有如下优点与有益效果:
1、从功能角度:本工具可对时序载荷等效处理,也可直接导入LDD载荷谱,可根据经验公式或静力学方程法求解静强度,可分别通过积分公式法或静力学方程法求解轴承疲劳寿命,同时可求解最大动态剪切应力及硬化层深度。
2、从用户角度:本工具根据目前各大主流轴承供应商内控标准,提出不同寿命评估方法,对于单排四点接触球轴承可分为单排轴承的运行转数,最大受载滚道的运行转数,单排轴承的运行时间,最大受载滚道的运行时间。对于双排四点接触球轴承:双排轴承的运行转数、单排轴承的运行转数,最大受载滚道的运行转数,双排轴承的运行时间,单排轴承的运行时间,最大受载滚道的运行时间。按转数进行评估时应以载荷谱中实际的转数为评价指标。按时间进行评估时应以20年的运行时间为评价指标(对于海上风电25年),全面满足用户的需求。
3、从评估角度:大大缩短了整机产品的设计周期,原先供应商校核结果需要1个月(含报告),采用本系统理想情况下可缩短至1天且可自动输出计算报告。
附图说明
图1为本发明的轴承受载后几何形变示意图。
图2为本发明的轴承套圈上的受力图。
图3为本发明的软件主界面图。
图4为本发明的偏航变桨轴承滚道校核模块界面图。
具体实施方式
下面结合具体实施例对本发明作进一步说明。
本实施例所提供的偏航变桨轴承滚道校核工具,为基于Matlab开发的偏航变桨轴承滚道校核软件,主要用于风力发电机偏航、变桨轴承初期选型、认证以及后期的结构强度评估,其主界面如图3所示。所述工具可直接利用LDD载荷谱进行疲劳寿命评估,同时可对时序载荷进等效处理,将每一个工况等效为一组载荷,等效的载荷谱的形式与LDD载荷谱形式一致。偏航变桨轴承一般采用单排四点接触球轴承或双排四点接触球轴承的结构形式,工具可根据经验公式或静力学方程法求解轴承静强度,也可分别通过积分公式法或静力学方程法求解轴承疲劳寿命,同时可求解最大动态剪切应力及硬化层深度。所述工具针对不同轴承供应商的轴承寿命提出以下寿命评估方法:
对于单排四点接触球轴承:可以根据单排轴承的运行转数、最大受载滚道的运行转数、单排轴承的运行时间或最大受载滚道的运行时间来评估轴承寿命。
对于双排四点接触球轴承:可以根据双排轴承的运行转数、单排轴承的运行转数、最大受载滚道的运行转数、双排轴承的运行时间、单排轴承的运行时间或最大受载滚道的运行时间评估寿命。
其中,按转数进行评估轴承寿命时应以载荷谱中实际的转数为评价指标;按时间进行评估轴承寿命时陆上风电应以20年的运行时间为评价指标,海上风电应以25年的运行时间为评价指标。
以下为本实施例上述偏航变桨轴承滚道校核工具的具体情况:
所述载荷谱的处理的情况如下:
时序载荷:
对每一个时序载荷工况的力、弯矩、速度通过方程1.1-1.3处理得到等效载荷工况表。
Figure GDA0002424747580000131
Figure GDA0002424747580000132
Figure GDA0002424747580000133
式1.1-1.3中,F为轴向力、径向力等效载荷,Fi为时序载荷中第i步的力,M为合弯矩等效载荷,Mi为时序载荷中第i步的弯矩,n为等效转速,ni为时序载荷中第i步的转速。
在20年(海上风电25年)运行时间内,将每一个疲劳工况的发生时间按下述方法求解:
对于发生次数的工况:
T20/25=L×T1×0.05×b (1.4)
对于发生时间的工况:
T20/25=6×L×T1×0.05×b (1.5)
式1.4-1.5中,L为风机的寿命要求,对于海上工况为25年,陆上工况为20年,T1为每个工况发生的次数或小时数,b为每个工况发生的步数,0.05表示步长,T20/25为20年或25年每个疲劳工况发生的时间。
对于LDD载荷可直接导入工具进行寿命计算。
所述静力学方程的情况如下:
静力学方程求解以双排四点接触球轴承为例进行分析。
1)四点接触球轴承的受载后几何形变如下
如图1所示,坐标系采用GL规范中叶根坐标系,假设轴承周边结构具有足够的刚性。外滚道沟渠率中心位置坐标系为xo,yo,内滚道沟渠率中心位置坐标系为xi,yi,轴承外圈固定,内圈相对于外圈产生轴向位移δa,径向位移δr,扭角θ。坐标系位置也由xi,yi转化xi',yi'再到xi”,yi”最后到xi”',yi”',沟渠率中心距离也由最初的MN变为MN”'。
对于双排滚子轴承,由排距d与转角θ对内外滚道产生0.5dθ变形。
对于单排轴承,不存在排距的影响;
2)双排四点接触球轴承的力学模型与负载分布计算
如图2所示,四个接触对在在位置角ψ处的接触力分别表示为Q、Q、Q、Q,根据Hertz点接触理论,法向接触载荷Q和接触变形δ的关系如下:
Figure GDA0002424747580000141
式中,在位置角ψj处,内圈受到轴向载荷Fa、径向载荷Fr、倾覆力矩M以及钢球对内滚道的接触载荷Q的作用;Kn为滚动体与内外圈总的负荷变形常数。
以内圈与滚动体为研究对象,内圈在外部载荷和所有滚动体载荷的作用下处于平衡状态,内圈的力学平衡方程为:
Figure GDA0002424747580000142
Figure GDA0002424747580000143
Figure GDA0002424747580000144
其中,α、α、α、α分别为受载后四个接触对在位置角ψ处的接触角。
上述方程通过Newton Raphson方法进行迭代求解,得到滚道精确的载荷分布。
3)单排四点接触球轴承的力学模型与负载分布计算
单排四点接触球轴承的静力学方程求解原理与上述双排四点接触球轴承的计算原理相同。
两个接触对在位置角ψ处的接触力分别表示为Q、Q,受载后两个接触对在位置角ψ处的接触角为α、α
故平衡方程为:
Figure GDA0002424747580000151
Figure GDA0002424747580000152
Figure GDA0002424747580000153
所述静强度计算的情况如下:
1)积分公式法的静强度根据下列公式求解
单排轴承:
Figure GDA0002424747580000154
双排轴承:
Figure GDA0002424747580000155
根据赫兹接触理论得:
Figure GDA0002424747580000156
Figure GDA0002424747580000157
式3.1-3.4中,Qmax为滚动体最大载荷,Fr为径向载荷,Fa为轴向载荷,M为合弯矩,Z为钢球总数,α为初始接触角,dm为节圆直径,a为接触椭圆短半轴,b接触椭圆的长半轴,Smin为最小安全系数,[σ]为许用接触应力,σmax最大接触应力。
2)静力学方程的静强度求解
根据方程2.2-2.7,通过Newton Raphson方法即可求出最大滚动体的载荷Qmax,其余的计算过程与积分公式法一致。
所述疲劳寿命计算情况如下:
1)滚道的基本额定滚动体负荷为:
Figure GDA0002424747580000161
其中,内滚道取上符号,外滚道取下符号。
Figure GDA0002424747580000162
Figure GDA0002424747580000163
式4.1-4.3中,Qci为内滚道的基本额定动负荷,Qce为外滚道的基本额定动负荷,fi为内滚道沟曲率半径系数,fo为外滚道沟曲率半径系数,γ为无量纲几何参数,为Dw为滚子直径,dm为滚道中心径,Z为钢球总个数。
2)积分公式求解单排四点接触球轴承当量滚动体负荷及接触滚道寿命:
2.1)当量滚动体负荷
相对负荷方向旋转套圈的当量滚动体负荷积分为:
Figure GDA0002424747580000164
相对负荷方向静止套圈的当量滚动体负荷积分为:
Figure GDA0002424747580000165
相对负荷方向旋转套圈(接触对1和接触对2的内滚道)的当量滚动体负荷为:
Qev1=Qev2=QmaxJ1(ε) (4.6)
相对负荷方向静止套圈(接触对1和接触对2的外滚道)的当量滚动体负荷为:
Qeu1=Qeu2=QmaxJ2(ε) (4.7)
式4.4-4.7中,ε为载荷系数,ψ为方位角,Qmax为滚动体最大载荷,可根据3.1式计算。
2.2)接触滚道的寿命(百万转)
相对负荷旋转套圈主推力滚道(接触对1的内圈)的额定寿命L10ev1为:
Figure GDA0002424747580000166
相对负荷静止套圈主推力滚道(接触对1的外圈)的额定寿命L10eu1为:
Figure GDA0002424747580000171
相对负荷旋转套圈主推力滚道(接触对2的内圈)的额定寿命L10ev2为:
Figure GDA0002424747580000172
相对负荷静止套圈副推力滚道(接触对2的外圈)的额定寿命L10eu2为:
Figure GDA0002424747580000173
3)积分公式求解双排四点接触球轴承当量滚动体负荷及接触滚道寿命
3.1)当量滚动体负荷
对于双排四点接触球轴承当量滚动体负荷积分如下。
相对负荷方向旋转套圈的当量滚动体负荷积分
Figure GDA0002424747580000174
相对负荷方向静止套圈的当量滚动体负荷积分
Figure GDA0002424747580000175
主推力滚道最大滚动体负荷为:
Figure GDA0002424747580000176
副推力滚道最大滚动体负荷为:
Figure GDA0002424747580000177
相对负荷方向旋转套圈的当量滚动体负荷为:
Qev13=Qmax1J1 (4.16)
Qev24=Qmax2J1 (4.17)
相对负荷方向静止套圈的当量滚动体负荷为:
Qeu13=Qmax1J2 (4.18)
Qeu24=Qmax2J2 (4.19)
3.2)接触滚道的寿命,按百万转计算
相对负荷旋转套圈主推力滚道的额定寿命L10ev13为:
Figure GDA0002424747580000181
相对负荷静止套圈主推力滚道的额定寿命L10eu13为:
Figure GDA0002424747580000182
相对负荷旋转套圈主推力滚道的额定寿命L10ev24为:
Figure GDA0002424747580000183
相对负荷静止套圈副推力滚道的额定寿命L10eu24为:
Figure GDA0002424747580000184
4)静力学方程求解单排四点接触球轴承当量滚动体负荷及接触滚道寿命
4.1)当量滚动体负荷
相对旋转套圈主推力滚道的当量滚动体负荷Qev1为:
Figure GDA0002424747580000185
相对静止套圈主推力滚道的当量滚动体负荷Qeu1为:
Figure GDA0002424747580000186
相对旋转套圈反推力滚道的当量滚动体负荷Qev2为:
Figure GDA0002424747580000187
相对静止套圈反推力滚道的当量滚动体负荷Qeu2为:
Figure GDA0002424747580000191
其中,Q是作用于主推力滚道各个滚动体上的负荷,Q是作用于反推力滚道各个滚动体上的负荷;
4.2)接触滚道的寿命,按百万转计算
接触滚道的寿命由方程4.8-4.11求解;
5)静力学方程求解双排四点接触球轴承当量滚动体负荷及接触滚道寿命
5.1)当量滚动体负荷
主推力滚道(接触对1、接触对3)的当量滚动体负荷为:
Figure GDA0002424747580000192
Figure GDA0002424747580000193
Figure GDA0002424747580000194
Figure GDA0002424747580000195
反推力滚道的当量滚动体负荷为:
Figure GDA0002424747580000196
Figure GDA0002424747580000197
Figure GDA0002424747580000198
Figure GDA0002424747580000201
其中,Q、Q是作用于主推力滚道各个滚动体上的负荷,Q、Q是作用于反推力滚道各个滚动体上的负荷,Qev1、Qev2、Qev3、Qev4分别为旋转套圈对应的四个接触位置的当量滚动体负荷,Qeu1、Qeu2、Qeu3、Qeu4分别为静止套圈对应的四个接触位置的当量滚动体负荷。
5.2)接触滚道的寿命,按百万转计算
相对负荷旋转套圈主推力滚道(接触对1、接触对3)的额定寿命分别为:
Figure GDA0002424747580000202
Figure GDA0002424747580000203
相对负荷静止套圈主推力滚道(接触对1、接触对3)的额定寿命分别为:
Figure GDA0002424747580000204
Figure GDA0002424747580000205
相对负荷旋转套圈副推力滚道(接触对2、接触对4)的额定寿命分别为:
Figure GDA0002424747580000206
Figure GDA0002424747580000207
相对负荷静止套圈副推力滚道(接触对2、接触对4)的额定寿命分别为:
Figure GDA0002424747580000208
Figure GDA0002424747580000209
所述寿命评估标准的情况如下:
根据不同的评估标准,寿命计算结果的最终处理也不一样。
1)积分公式求解单排四点接触球轴承寿命
整个轴承的运行转数为:
L10=[L10ev1 (-10/9)+L10eu1 (-10/9)+L10ev2 (-10/9)+L10eu2 (-10/9)]-0.9×106 (6.1)
最大受载滚道的运行转数为:
L10=[L10ev1 (-10/9)+L10eu1 (-10/9)]-0.9×106 (6.2)
2)积分公式求解双排四点接触球轴承寿命
整个轴承的运行转数为:
L10=[L10ev13 (-10/9)+L10eu13 (-10/9)+L10ev24 (-10/9)+L10eu24 (-10/9)]-0.9×106 (6.3)
单个轴承(受载最大)的运行转数为:
L10=[L10ev13 (-10/9)+L10eu13 (-10/9)]-0.9×106 (6.4)
最大受载滚道的运行转数为:
L10=L10ev13×106 (6.5)
3)静力学方程求解单排四点接触球轴承寿命
方法与式子6.1-6.2一致。
4)静力学方程求解双排四点接触球轴承寿命
整个轴承的运行转数为:
Figure GDA0002424747580000211
单个轴承(受载最大)的运行转数为:
L10=[L10ev1 (-10/9)+L10eu1 (-10/9)+L10ev2 (-10/9)+L10eu2 (-10/9)]-0.9×106 (6.7)
最大受载滚道的运行转数为:
L10=[L10ei1 (-10/9)+L10ee1 (-10/9)]-0.9×106 (6.8)
5)按时间
如根据运行时间评估轴承寿命,需要在上述方程的基础上按下式进行求解
Figure GDA0002424747580000212
其中,n为等效转速。
所述最大动态剪切应力及硬化层深度的情况如下:
Figure GDA0002424747580000221
Figure GDA0002424747580000222
Figure GDA0002424747580000223
式7.1-7.3中,τ0为最大动态剪切应力,Q为接触力,t为滚道与钢球的接触尺寸的结构参数,z0为硬化层深度。
如图4所示,所述工具的参数输入与结果输出的情况如下:
工具界面分为五部分:轴承输入参数、计算方法、载荷导入、计算结果以及结果输出。
在轴承输入参数中需要输入钢球直径、轴承节圆直径、内滚道沟曲率半径系数、外滚道沟曲率半径系数、初始接触角、单排钢球数目、列数、双排滚子间距、轴承游隙、钢球弹性模量、内套圈弹性模量、外套圈弹性模量、内套圈泊松比、外套圈泊松比、滚道硬度、许用接触应力。
在计算方法中选择积分公式或Newton Raphson方法;
在载荷导入中导入包含极限载荷与LDD载荷谱的文件,在载荷输入方法中可以选择LDD或TS(时序载荷),当选择LDD时,无需选择时序路径,当选择TS时需要在载荷导入中选择时序的文件夹。
在计算结果中显示最终计算结果,其中包括钢球受到最大接触力、钢球受到最大接触应力、最小静安全系数、最大动态剪切应力、硬化层深度、计算寿命。
在结果输出中显示详细的计算结果报告的存放路径。
以上所述实施例只为本发明之较佳实施例,并非以此限制本发明的实施范围,故凡依本发明之形状、原理所作的变化,均应涵盖在本发明的保护范围内。

Claims (3)

1.一种偏航变桨轴承滚道校核工具,其特征在于:所述工具为基于Matlab开发的偏航变桨轴承滚道校核软件,主要用于风力发电机偏航、变桨轴承初期选型、认证以及后期的结构强度评估;
所述工具直接利用LDD载荷谱进行疲劳寿命评估,同时能够对时序载荷进等效处理,将每一个工况等效为一组载荷,等效的载荷谱的形式与LDD载荷谱形式一致;
偏航变桨轴承通常采用单排四点接触球轴承或双排四点接触球轴承的结构形式,所述工具能够根据经验公式或静力学方程法求解轴承静强度,也能够通过积分公式法或静力学方程法求解轴承疲劳寿命,同时求解最大动态剪切应力及硬化层深度;
所述工具针对不同轴承供应商的轴承寿命提出以下寿命评估方法:
对于单排四点接触球轴承,根据单排轴承的运行转数、最大受载滚道的运行转数、单排轴承的运行时间或最大受载滚道的运行时间来评估轴承寿命;
对于双排四点接触球轴承,根据双排轴承的运行转数、单排轴承的运行转数、最大受载滚道的运行转数、双排轴承的运行时间、单排轴承的运行时间或最大受载滚道的运行时间评估寿命;
其中,按转数进行评估轴承寿命时应以载荷谱中实际的转数为评价指标;按时间进行评估轴承寿命时陆上风电应以20年的运行时间为评价指标,海上风电应以25年的运行时间为评价指标。
2.根据权利要求1所述的一种偏航变桨轴承滚道校核工具,其特征在于:
所述时序载荷的处理情况如下:
时序载荷:
对每一个时序载荷工况的力、弯矩、速度通过下面方程1.1-1.3处理得到等效载荷工况表;
Figure FDA0002424747570000011
Figure FDA0002424747570000012
Figure FDA0002424747570000013
式1.1-1.3中,F为轴向力、径向力等效载荷,Fi为时序载荷中第i步的力,M为合弯矩等效载荷,Mi为时序载荷中第i步的弯矩,n为等效转速,ni为时序载荷中第i步的转速;
将每一个疲劳工况的发生时间按下述方法求解:
对于发生次数的工况:
T20/25=L×T1×0.05×b (1.4)
对于发生时间的工况:
T20/25=6×L×T1×0.05×b (1.5)
式1.4-1.5中,L为风机的寿命要求,对于海上工况为25年,陆上工况20年,T1为每个工况发生的次数或小时数,b为每个工况发生的步数,0.05表示步长,T20/25为20年或25年每个疲劳工况发生的时间;
对于LDD载荷则可直接导入工具进行寿命计算;
所述静力学方程的情况如下:
1)四点接触球轴承的受载后几何形变如下:
坐标系采用GL规范中叶根坐标系,假设轴承周边结构具有足够的刚性,外滚道沟渠率中心位置坐标系为xo,yo,内滚道沟渠率中心位置坐标系为xi,yi,轴承外圈固定,内圈相对于外圈产生轴向位移δa,径向位移δr,扭角θ;坐标系位置也由xi,yi转化xi',yi'再到xi”,yi”最后到xi”',yi”',沟渠率中心距离也由最初的MN变为MN”';
对于双排滚子轴承,由排距d与转角θ对内外滚道产生0.5dθ变形;
对于单排轴承,不存在排距的影响;
2)双排四点接触球轴承的力学模型与负载分布计算
四个接触对在位置角ψ处的接触力分别表示为
Figure FDA0002424747570000024
根据Hertz点接触理论,法向接触载荷
Figure FDA0002424747570000025
和接触变形
Figure FDA0002424747570000026
的关系如下:
Figure FDA0002424747570000021
式中,在角位置ψj处,内圈受到轴向载荷Fa、径向载荷Fr、倾覆力矩M以及钢球对内滚道的接触载荷
Figure FDA0002424747570000027
的作用;Kn为滚动体与内外圈总的负荷变形常数;
以内圈与滚动体为研究对象,内圈在外部载荷和所有滚动体载荷的作用下处于平衡状态,内圈的力学平衡方程为:
Figure FDA0002424747570000022
Figure FDA0002424747570000023
Figure FDA0002424747570000031
其中,
Figure FDA0002424747570000038
分别为受载后四个接触对在位置角ψ处的接触角;
上述方程通过Newton Raphson方法进行迭代求解,得到滚道精确的载荷分布;
3)单排四点接触球轴承的力学模型与负载分布计算
单排四点接触球轴承的静力学方程求解原理与上述双排四点接触球轴承的计算原理相同;
两个接触对在位置角ψ处的接触力分别表示为
Figure FDA0002424747570000039
受载后两个接触对在位置角ψ处的接触角为
Figure FDA00024247475700000310
故平衡方程为:
Figure FDA0002424747570000032
Figure FDA0002424747570000033
Figure FDA0002424747570000034
所述静强度计算的情况具体如下:
1)积分公式法的静强度根据下列公式求解
单排轴承:
Figure FDA0002424747570000035
双排轴承:
Figure FDA0002424747570000036
根据赫兹接触理论得:
Figure FDA0002424747570000037
Figure FDA0002424747570000041
式3.1-3.4中,Qmax为滚动体最大载荷,Fr为径向载荷,Fa为轴向载荷,M为合弯矩,Z为钢球总数,α为初始接触角,dm为节圆直径,a为接触椭圆短半轴,b为接触椭圆的长半轴,Smin为最小安全系数,[σ]为许用接触应力,σmax为最大接触应力;
2)静力学方程的静强度求解
根据方程2.2-2.7,通过Newton Raphson方法即可求出最大滚动体的载荷Qmax,其余的计算过程与积分公式法一致;
所述疲劳寿命的计算情况如下:
1)滚道的基本额定滚动体负荷:
Figure FDA0002424747570000042
其中,内滚道取上符号,外滚道取下符号;
Figure FDA0002424747570000043
Figure FDA0002424747570000044
式4.1-4.3中,Qci为内滚道的基本额定动负荷,Qce为外滚道的基本额定动负荷,fi为内滚道沟曲率半径系数,fo为外滚道沟曲率半径系数,γ为无量纲几何参数,Dw为滚子直径,dm为节圆直径,Z为钢球总个数;
2)积分公式求解单排四点接触球轴承当量滚动体负荷及接触滚道寿命
2.1)当量滚动体负荷
相对负荷方向旋转套圈的当量滚动体负荷积分为:
Figure FDA0002424747570000045
相对负荷方向静止套圈的当量滚动体负荷积分为:
Figure FDA0002424747570000046
相对负荷方向旋转套圈的当量滚动体负荷为:
Qev1=Qev2=QmaxJ1(ε) (4.6)
相对负荷方向静止套圈的当量滚动体负荷为:
Qeu1=Qeu2=QmaxJ2(ε) (4.7)
式4.4-4.7中,ε为载荷系数,ψ为位置角,Qmax为滚动体最大载荷,根据式3.1计算;
2.2)接触滚道的寿命,按百万转计算
相对负荷旋转套圈主推力滚道的额定寿命L10ev1为:
Figure FDA0002424747570000051
相对负荷静止套圈主推力滚道的额定寿命L10eu1为:
Figure FDA0002424747570000052
相对负荷旋转套圈主推力滚道的额定寿命L10ev2为:
Figure FDA0002424747570000053
相对负荷静止套圈副推力滚道的额定寿命L10eu2为:
Figure FDA0002424747570000054
3)积分公式求解双排四点接触球轴承当量滚动体负荷及接触滚道寿命
3.1)当量滚动体负荷
对于双排四点接触球轴承当量滚动体负荷积分如下:
相对负荷方向旋转套圈的当量滚动体负荷积分为:
Figure FDA0002424747570000055
相对负荷方向静止套圈的当量滚动体负荷积分为:
Figure FDA0002424747570000056
主推力滚道最大滚动体负荷为:
Figure FDA0002424747570000061
副推力滚道最大滚动体负荷为:
Figure FDA0002424747570000062
相对负荷方向旋转套圈的当量滚动体负荷为:
Qev13=Qmax1J1 (4.16)
Qev24=Qmax2J1 (4.17)
相对负荷方向静止套圈的当量滚动体负荷为:
Qeu13=Qmax1J2 (4.18)
Qeu24=Qmax2J2 (4.19)
3.2)接触滚道的寿命,按百万转计算
相对负荷旋转套圈主推力滚道的额定寿命L10ev13为:
Figure FDA0002424747570000063
相对负荷静止套圈主推力滚道的额定寿命L10eu13为:
Figure FDA0002424747570000064
相对负荷旋转套圈主推力滚道的额定寿命L10ev24为:
Figure FDA0002424747570000065
相对负荷静止套圈副推力滚道的额定寿命L10eu24为:
Figure FDA0002424747570000066
4)静力学方程求解单排四点接触球轴承当量滚动体负荷及接触滚道寿命
4.1)当量滚动体负荷
相对旋转套圈主推力滚道的当量滚动体负荷Qev1为:
Figure FDA0002424747570000071
相对静止套圈主推力滚道的当量滚动体负荷Qeu1为:
Figure FDA0002424747570000072
相对旋转套圈反推力滚道的当量滚动体负荷Qev2为:
Figure FDA0002424747570000073
相对静止套圈反推力滚道的当量滚动体负荷Qeu2为:
Figure FDA0002424747570000074
其中,
Figure FDA0002424747570000078
是作用于主推力滚道各个滚动体上的负荷,
Figure FDA0002424747570000079
是作用于反推力滚道各个滚动体上的负荷;
4.2)接触滚道的寿命,按百万转计算
接触滚道的寿命由方程4.8-4.11求解;
5)静力学方程求解双排四点接触球轴承当量滚动体负荷及接触滚道寿命
5.1)当量滚动体负荷
主推力滚道的当量滚动体负荷为:
Figure FDA0002424747570000075
Figure FDA0002424747570000076
Figure FDA0002424747570000077
Figure FDA0002424747570000081
反推力滚道的当量滚动体负荷为:
Figure FDA0002424747570000082
Figure FDA0002424747570000083
Figure FDA0002424747570000084
Figure FDA0002424747570000085
其中,
Figure FDA00024247475700000810
是作用于主推力滚道各个滚动体上的负荷,
Figure FDA00024247475700000811
是作用于反推力滚道各个滚动体上的负荷,Qev1、Qev2、Qev3、Qev4分别为旋转套圈对应的四个接触位置的当量滚动体负荷,Qeu1、Qeu2、Qeu3、Qeu4分别为静止套圈对应的四个接触位置的当量滚动体负荷;
5.2)接触滚道的寿命,按百万转计算
相对负荷旋转套圈主推力滚道的额定寿命分别为:
Figure FDA0002424747570000086
Figure FDA0002424747570000087
相对负荷静止套圈主推力滚道的额定寿命分别为:
Figure FDA0002424747570000088
Figure FDA0002424747570000089
相对负荷旋转套圈副推力滚道的额定寿命分别为:
Figure FDA0002424747570000091
Figure FDA0002424747570000092
相对负荷静止套圈副推力滚道的额定寿命分别为:
Figure FDA0002424747570000093
Figure FDA0002424747570000094
所述寿命评估标准的情况如下:
1)积分公式求解单排四点接触球轴承寿命
整个轴承的运行转数为:
L10=[L10ev1 (-10/9)+L10eu1 (-10/9)+L10ev2 (-10/9)+L10eu2 (-10/9)]-0.9×106 (6.1)
最大受载滚道的运行转数为:
L10=[L10ev1 (-10/9)+L10eu1 (-10/9)]-0.9×106 (6.2)
2)积分公式求解双排四点接触球轴承寿命
整个轴承的运行转数为:
L10=[L10ev13 (-10/9)+L10eu13 (-10/9)+L10ev24 (-10/9)+L10eu24 (-10/9)]-0.9×106 (6.3)
单个轴承最大受载下的运行转数为:
L10=[L10ev13 (-10/9)+L10eu13 (-10/9)]-0.9×106 (6.4)
最大受载滚道的运行转数为:
L10=L10ev13×106 (6.5)
3)静力学方程求解单排四点接触球轴承寿命,其方法与方程6.1-6.2相同;
4)静力学方程求解双排四点接触球轴承寿命
整个轴承的运行转数为:
Figure FDA0002424747570000095
单个轴承最大受载下的运行转数为:
L10=[L10ev1 (-10/9)+L10eu1 (-10/9)+L10ev2 (-10/9)+L10eu2 (-10/9)]-0.9×106 (6.7)
最大受载滚道的运行转数为:
L10=[L10ei1 (-10/9)+L10ee1 (-10/9)]-0.9×106 (6.8)
5)按时间
根据运行时间评估轴承寿命,需要在上述方程的基础上按下式进行求解
Figure FDA0002424747570000101
其中,n为等效转速;
所述最大动态剪切应力及硬化层深度的情况如下:
Figure FDA0002424747570000102
Figure FDA0002424747570000103
Figure FDA0002424747570000104
式7.1-7.3中,τ0为最大动态剪切应力,Q为接触力,t为滚道与钢球的接触尺寸的结构参数,z0为硬化层深度。
3.根据权利要求1所述的一种偏航变桨轴承滚道校核工具,其特征在于:所述工具的参数输入与结果输出的情况如下:
工具界面分为五部分:轴承输入参数、计算方法、载荷导入、计算结果以及结果输出;
在轴承输入参数中需要输入钢球直径、轴承节圆直径、内滚道沟曲率半径系数、外滚道沟曲率半径系数、初始接触角、单排钢球数目、列数、双排滚子间距、轴承游隙、钢球弹性模量、内套圈弹性模量、外套圈弹性模量、内套圈泊松比、外套圈泊松比、滚道硬度、许用接触应力;
在计算方法中选择积分公式或Newton Raphson方法;
在载荷导入中导入包含极限载荷与LDD载荷谱的文件,在载荷输入方法中选择LDD或时序载荷TS,当选择LDD时,无需选择时序路径,当选择TS时需要在载荷导入中选择时序的文件夹;
在计算结果中显示最终计算结果,包括钢球受到最大接触力、钢球受到最大接触应力、最小静安全系数、最大动态剪切应力、硬化层深度、计算寿命;
在结果输出中显示详细的计算结果报告的存放路径。
CN201710753617.4A 2017-08-29 2017-08-29 一种偏航变桨轴承滚道校核工具 Active CN107451376B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710753617.4A CN107451376B (zh) 2017-08-29 2017-08-29 一种偏航变桨轴承滚道校核工具

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710753617.4A CN107451376B (zh) 2017-08-29 2017-08-29 一种偏航变桨轴承滚道校核工具

Publications (2)

Publication Number Publication Date
CN107451376A CN107451376A (zh) 2017-12-08
CN107451376B true CN107451376B (zh) 2020-07-14

Family

ID=60494219

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710753617.4A Active CN107451376B (zh) 2017-08-29 2017-08-29 一种偏航变桨轴承滚道校核工具

Country Status (1)

Country Link
CN (1) CN107451376B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103256908A (zh) * 2013-05-22 2013-08-21 河南科技大学 一种变桨轴承径向游隙的确定方法
CN103324805A (zh) * 2013-07-04 2013-09-25 沈阳华创风能有限公司 风机底架极限强度校核的计算方法
CN105550403A (zh) * 2015-12-04 2016-05-04 广东明阳风电产业集团有限公司 一种风机外购件强度校核工具

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103256908A (zh) * 2013-05-22 2013-08-21 河南科技大学 一种变桨轴承径向游隙的确定方法
CN103324805A (zh) * 2013-07-04 2013-09-25 沈阳华创风能有限公司 风机底架极限强度校核的计算方法
CN105550403A (zh) * 2015-12-04 2016-05-04 广东明阳风电产业集团有限公司 一种风机外购件强度校核工具

Also Published As

Publication number Publication date
CN107451376A (zh) 2017-12-08

Similar Documents

Publication Publication Date Title
CN107729597B (zh) 一种主轴轴承滚道校核工具
CN107704652B (zh) 一种风力发电机组轴承刚度计算工具
CN111896255B (zh) 基于多轴随机道路载荷下轮毂轴承服役寿命快速评估方法
EP2532904A2 (en) Structure of slewing ring bearing
CN102819635B (zh) 特大型双排四点接触球轴承静力学模型建立方法
CN105550403B (zh) 一种风机外购件强度校核工具
WO2007061439A3 (en) Wind turbine pitch bearing and method
CN109753723B (zh) 一种向心滚动轴承疲劳寿命计算方法
Stammler et al. Friction torque of wind-turbine pitch bearings–comparison of experimental results with available models
CN108843760A (zh) 变速装置
CN103256908B (zh) 一种变桨轴承径向游隙的确定方法
CN108692938B (zh) 一种获取滚动轴承寿命的方法
CN103674546B (zh) 风电转盘轴承加速疲劳寿命试验方法
CN112861271B (zh) 风电机组主轴承力学特性的计算方法、系统及介质
CN107451376B (zh) 一种偏航变桨轴承滚道校核工具
CN103174740A (zh) 一种四点接触球轴承沟曲率半径系数的设计方法
CN107563081B (zh) 高速空心圆柱滚子轴承的参数优化方法
RU2626800C1 (ru) Способ комплектования шарикоподшипников
CN106934165A (zh) 一种交叉圆锥滚子轴承静刚度的快速提取方法
Rabréau et al. Influence of bearing kinematics hypotheses on ball bearing heat generation
CN111859575A (zh) 一种单列四点接触球转盘轴承额定动静载荷计算方法
Sun et al. Rolling-element bearings in China: From ancient times to the 20th century
CN106640559A (zh) 非对称双排四点接触球轴承非对称角的设置方法以及装置
CN211288484U (zh) 一种s形定位回转径结构的偏航轴承
CN211174993U (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
TR01 Transfer of patent right

Effective date of registration: 20201216

Address after: No.2, Luohu West Road, Houpu village committee, Jieshi Town, Lufeng City, Shanwei City, Guangdong Province

Patentee after: Shanwei Mingyang New Energy Technology Co.,Ltd.

Address before: 528437 No. 22 Torch Road, Torch Development Zone, Zhongshan City, Guangdong Province

Patentee before: MING YANG SMART ENERGY GROUP Co.,Ltd.

TR01 Transfer of patent right