CN112861271A - 风电机组主轴承力学特性的计算方法、系统及介质 - Google Patents
风电机组主轴承力学特性的计算方法、系统及介质 Download PDFInfo
- Publication number
- CN112861271A CN112861271A CN201911182046.9A CN201911182046A CN112861271A CN 112861271 A CN112861271 A CN 112861271A CN 201911182046 A CN201911182046 A CN 201911182046A CN 112861271 A CN112861271 A CN 112861271A
- Authority
- CN
- China
- Prior art keywords
- main bearing
- contact
- load
- wind turbine
- calculating
- 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 44
- 238000005096 rolling process Methods 0.000 claims abstract description 51
- 238000013459 approach Methods 0.000 claims abstract description 33
- 230000014509 gene expression Effects 0.000 claims abstract description 21
- 238000006073 displacement reaction Methods 0.000 claims abstract description 8
- 238000004590 computer program Methods 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 16
- 238000011156 evaluation Methods 0.000 claims description 4
- 230000008569 process Effects 0.000 claims description 4
- 230000005489 elastic deformation Effects 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 230000002035 prolonged effect Effects 0.000 claims description 2
- 238000013461 design Methods 0.000 abstract description 13
- 230000005540 biological transmission Effects 0.000 abstract description 11
- 230000006872 improvement Effects 0.000 description 6
- 238000011161 development Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Landscapes
- Rolling Contact Bearings (AREA)
- Wind Motors (AREA)
Abstract
本发明公开了一种风电机组主轴承力学特性的计算方法及系统，包括：根据风轮载荷、传动链布局参数、仰角α，主轴重量Ms，齿轮箱重量Mg，求解作用于主轴承上的轴向载荷Fa和径向载荷Fr；以轴向载荷Fa和径向载荷Fr作为输入量，分析主轴承内圈、滚动体的位移，计算滚动体的弹性趋近量δ的表达式，建立内圈的平衡方程；模拟套圈与滚动体的接触，建立滚动体接触载荷Q与弹性趋近量δ、最大接触应力ph的方程；联立各方程，运用迭代法进行各参数求解。本发明可实现风电机组主轴承接触应力、接触载荷、接触角、寿命等关键力学特性参数的精确求解，实现主轴承的快速选型以及传动链的精益设计，同时为主轴承的设计提供一定的理论支撑。
Description
技术领域
本发明主要涉及风力发电技术领域，特指一种风电机组主轴承力学特性的计算方法、系统及介质。
背景技术
风电是一种无污染的绿色能源，开发利用风能，大力发展风力发电技术，对于解决全球性能源和环境污染具有重大的意义。在风电行业竞争日益激烈的今天，不同功率等级、叶轮直径、风速等级的整机产品层出不穷，要求整机厂商能快速完成风电机组迭代、优化设计，同时保证整机成本最优化。2MW+双馈机组中通过设置主轴承来承担风轮传递过来的大部分非扭矩的载荷(全部的推力载荷、绝大部分的风轮重力载荷)。主轴承作为风电机组传动系统的重要部件，其选型直接影响到自身的可靠性、低速轴系的重量及成本，主轴承的选型校核速度的快慢很大程度上影响了机组设计周期。
风电机组中常选用双列球面滚子轴承作为主轴承，具有自调心功能，可承受齿轮箱所要求的较大偏转角，并可承受较大的径向载荷及2个方向的重轴向载荷。但其校核方法一直由轴承供应商掌握，不对外公开，整机厂商掌握的数据很少，很大程度上制约了风电机组的设计开发周期，由于整机厂商对双列球面滚子轴承的了解不够深入，也导致传动链的降本效率低下，甚至无法入手。
发明内容
本发明要解决的技术问题就在于：针对现有技术存在的技术问题，本发明提供一种能够完成主轴承接触应力、寿命关键评价指标的计算，实现主轴承的快速选型以及风电机组传动链精益设计的风电机组主轴承力学特性的计算方法、系统、介质及设备。
为解决上述技术问题，本发明提出的技术方案为：
一种风电机组主轴承力学特性的计算方法，包括步骤：
S01、根据风轮载荷、传动链布局参数、仰角α，主轴重量Ms，齿轮箱重量Mg，求解作用于主轴承上的轴向载荷Fa和径向载荷Fr；
S02、以轴向载荷Fa和径向载荷Fr作为输入量，分析主轴承内圈、滚动体的位移，计算滚动体的弹性趋近量δ的表达式，建立内圈的平衡方程；
S03、模拟套圈与滚动体的接触，建立滚动体接触载荷Q与弹性趋近量δ、最大接触应力ph的方程；
S04、联立步骤S02和步骤S03中的方程，运用迭代法进行各参数求解。
作为上述技术方案的进一步改进，在步骤S01中，主轴承上的轴向载荷Fa和径向载荷Fr的表达式为：
式中：[Mx，My，Mz，Fx，Fy，Fz]为风轮载荷；坐标系采用GL坐标系；A、B、C、D为传动链布局参数。
作为上述技术方案的进一步改进，在步骤S02中，滚动体的弹性趋近量δ的表达式为：
式中：δ1j，δ2j分别为第1，2列滚子弹性趋近量；A0为轴承内、外圈滚道曲率中心初始距离；α0为初始接触角；
为第j个滚子方位角；Dw为滚子直径，cd为双列球面滚子轴承的径向工作游隙。
作为上述技术方案的进一步改进，在步骤S02中，所述内圈平衡方程表达式如下：
式中：Z为滚子数量；α1j、α2j分别为第j个滚子受载后的接触角，其中
作为上述技术方案的进一步改进，在步骤S03中，滚动体接触载荷Q与弹性趋近量δ、最大接触应力p的关系式为：
其中：
式中：bk为两切片间距长度的一半，N为切片的总数量；z(yl)为滚子第l个切片处与滚道表面的初始间距，Dl,k为柔度系数。
作为上述技术方案的进一步改进，步骤S04的具体过程为：
根据式(1)，求解出轴承的轴向载荷Fa和径向载荷Fr，作为后续方程组的输入；
先赋初值δa、δr，根据式(2)、(3)、(7)、(8)、(9)求解两列滚动体的弹性趋近量δ，再对接触区域上各切片的接触半宽ak赋初值，求解柔度系数Dl,k，运用式(12)的变形式p＝πE’×(D-1*W)得到滚动体上的接触应力，其中W为每个切片中心的弹性变形矩阵，可根据各方位角滚子与内、外圈之间的弹性趋近量及滚子与滚道之间的凸度来计算；
然后根据有限长一维接触问题的前提条件及Hertz线接触理论推导出各切片的接触半宽计算公式ak＝2Rpkmax/E’，其中R为滚子-滚道接触体在切片处的综合曲率半径；通过与接触半宽ak的初值或上次的计算值进行反复迭代，完成给定条件下滚子-滚道之间的接触应力求解p；
根据式(11)求解Q，再根据已求出的弹性趋近量δ求解轴向载荷Fa和径向载荷Fr，通过对比利用式(1)求解出来的Fa、Fr，进行反复迭代，直到满足精度要求；
求解完成后可输出主轴承接触应力关键评价指标及寿命，对比预设标准应力和预设标准寿命，可完成主轴承的校核。
作为上述技术方案的进一步改进，在步骤S03中，采用一维有限长接触理论模拟套圈与滚动体的接触。
本发明还公开了一种风电机组主轴承力学特性的计算系统，包括
第一模块，用于根据风轮载荷、传动链布局参数、仰角α，主轴重量Ms，齿轮箱重量Mg，求解作用于主轴承上的轴向载荷Fa和径向载荷Fr；
第二模块，用于以轴向载荷Fa和径向载荷Fr作为输入量，分析主轴承内圈、滚动体的位移，计算滚动体的弹性趋近量δ的表达式，建立内圈的平衡方程；
第三模块，用于模拟套圈与滚动体的接触，建立滚动体接触载荷Q与弹性趋近量δ、最大接触应力ph的方程；
第四模块，用于联立各方程，运用迭代法进行各参数求解。
本发明进一步公开了一种计算机可读存储介质，其上储存有计算机程序，所述计算机程序在被处理器运行时执行如上所述的风电机组主轴承力学特性的计算方法的步骤。
本发明还公开一种计算机设备，包括存储器和处理器，所述存储器上存储有计算机程序，所述计算机程序在被处理器运行时执行如上所述的风电机组主轴承力学特性的计算方法的步骤。
与现有技术相比，本发明的优点在于：
本发明的风电机组主轴承力学特性的计算方法，针对目前主轴承力学特性计算缺乏使用理论模型，主轴承校核依赖供应商，导致校核周期长，影响整机开发进度等问题，基于有限长接触理论，建立了风电机组双列球面滚子轴承力学特性计算模型，准确对主轴承的接触应力、接触载荷、接触角、寿命等关键力学特性参数的精确求解，能够实现主轴承的自主快速校核，为实现风电机组的精益设计、主轴承的设计提供一定理论支撑；使设计人员更加深入了解主轴承的承载特性，帮助设计人员进行风电机组主轴承的快速选型，实现风电机组传动链的精益设计，最终实现高可靠性、低成本风电机组的优化设计。
(2)本发明采用一维有限长接触理论来模拟分析，相对比球轴承的赫兹接触理论，可模拟滚动体的边缘接触应力，计算更加精确；相对比二维有限长接触理论，计算速度更快，更适用于工程计算。
(3)本发明在实现主轴承力学特性求解的基础上，加入了机组风轮载荷、传动链布局参数等，可用于研究整机载荷、布局参数对主轴承力学特性的影响，实现风电机组主轴承的快速选型、传动链的正向设计。
附图说明
图1为本发明的方法在实施例的流程图。
图2为本发明的主轴承的结构示意图。
图3为本发明的风电机组传动链布局参数示意图。
图4为风电机组主轴承滚动体接触载荷示意图。
图5为风电机组主轴承接触应力示意图。
具体实施方式
以下结合说明书附图和具体实施例对本发明作进一步描述。
如图1至图5所示，本发明的风电机组主轴承力学特性的计算方法，包括步骤：
S01、根据风轮载荷、传动链布局参数、仰角α，主轴重量Ms，齿轮箱重量Mg，求解作用于主轴承上的轴向载荷Fa和径向载荷Fr；
S02、以轴向载荷Fa和径向载荷Fr作为输入量，分析主轴承内圈、滚动体的位移，计算滚动体的弹性趋近量δ的表达式，建立内圈的平衡方程；
S03、模拟套圈与滚动体的接触，建立滚动体接触载荷Q与弹性趋近量δ、最大接触应力ph的方程；
S04、联立步骤S02和步骤S03中的方程，运用迭代法进行各参数求解。
本发明的风电机组主轴承力学特性的计算方法，针对目前主轴承力学特性计算缺乏使用理论模型，主轴承校核依赖供应商，导致校核周期长，影响整机开发进度等问题，基于有限长接触理论，建立了风电机组双列球面滚子轴承力学特性计算模型，准确对主轴承的接触应力、接触载荷、接触角、寿命等关键力学特性参数的精确求解，能够实现主轴承的自主快速校核，为实现风电机组的精益设计、主轴承的设计提供一定理论支撑；使设计人员更加深入了解主轴承的承载特性，帮助设计人员进行风电机组主轴承的快速选型，实现风电机组传动链的精益设计，最终实现高可靠性、低成本风电机组的优化设计。
本实施例中，在步骤S01中，根据风轮载荷、传动链布局参数(包括A、B、C、D，如图3所示)、仰角α、主轴重量Ms、齿轮箱重量Mg，利用理论力学模型求解作用于主轴承上的轴向载荷Fa和径向载荷Fr，表达式为：
式中：[Mx，My，Mz，Fx，Fy，Fz]为风轮载荷，坐标系采用GL坐标系。
本实施例中，在步骤S02中，以轴向载荷Fa和径向载荷Fr作为输入量，分析主轴承内圈与滚动体的位移，计算滚动体的弹性趋近量δ的表达式，建立内圈的平衡方程，如图2所示，主轴承包括外圈、内圈以及滑动设置于外圈与内圈之间的滚动体。其中滚动体的弹性趋近量δ的表达式为：
式中：δ1j，δ2j分别为第1,2列滚子弹性趋近量；A0为轴承内、外圈滚道曲率中心初始距离；α0为初始接触角；
为第j个滚子方位角；Dw为滚子直径，cd为双列球面滚子轴承的径向工作游隙。
其中内圈平衡方程表达式如下：
式中：Z为滚子数量；α1j、α2j分别为第j个滚子受载后的接触角，其中
如图4所示，本实施例中，在步骤S03中，采用一维有限长接触理论模拟套圈与滚动体的接触，采用切片法，建立滚动体接触载荷Q与弹性趋近量δ、最大接触应力p的关系式，如下：
其中：
式中：bk为两切片间距长度的一半，N为切片的总数量；z(yl)为滚子第l个切片处与滚道表面的初始间距，Dl,k为柔度系数。
本发明采用一维有限长接触理论来模拟分析，相对比球轴承的赫兹接触理论，可模拟滚动体的边缘接触应力，计算更加精确；相对比二维有限长接触理论，计算速度更快，更适用于工程计算。
本实施例中，在步骤S04中，对各方程进行求解，具体过程为：
根据式(1)，求解出轴承的轴向载荷Fa和径向载荷Fr，作为后续方程组的输入；
由于δ、Q、p、Dl,k、α等变量的相互耦合，重难点是进行式(2)-式(13)联立的求解，具体方法为：
先赋初值δa、δr，根据式(2)、(3)、(7)、(8)、(9)求解两列滚动体的弹性趋近量δ，再对接触区域上各切片的接触半宽ak赋初值，求解柔度系数Dl,k，运用(12)式的变形式p＝πE’×(D-1*W)得到滚动体上的接触应力，其中W为每个切片中心的弹性变形矩阵，可根据各方位角滚子与内、外圈之间的弹性趋近量及滚子与滚道之间的凸度来计算；
然后根据有限长一维接触问题的前提条件及Hertz线接触理论可推导出各切片的接触半宽计算公式ak＝2Rpkmax/E’(R为滚子-滚道接触体在切片处的综合曲率半径)；
通过与接触半宽ak的初值(或上次的计算值)进行反复迭代，完成给定条件下滚子-滚道之间的接触应力求解p，如图5所示；
根据式(11)求解Q，再根据已求出的弹性趋近量δ求解轴向载荷Fa、径向载荷Fr，通过对比利用方程1求解出来的Fa、Fr，进行反复迭代，直到满足精度要求，求解完成后可输出主轴承接触应力等关键评价指标，寿命可根据ISO16281标准求解，对比GL标准中轴承的应力、寿命要求，即可完成主轴承的校核。
以某双馈风电机组为研究对象，其中布局参数A＝2240mm，B＝470mm，C＝2740mm，D＝505mm，仰角α＝5°，主轴重量Ms＝12t，齿轮箱重量Mg＝22t，其主轴承结构参数初始接触角αo＝10.4°，滚子直径Dw＝82mm，节圆直径dm＝928mm，滚子长度L＝134mm，内圈沟道曲率半径Ri＝516mm，外圈沟道曲率半径Ro＝513mm，滚子曲率半径r＝506mm，滚子数量2×30个，材料泊松比ν＝0.3，弹性模量E＝207Gpa，径向工作游隙cd＝224μm。输入载荷如下：
表1：风轮载荷
| 工况 | Fx(kN) | Fy(kN) | Fz(kN) | My(kNm) | Mz(kNm) | 持续时间(h) |
| 1 | 306.32 | -13.55 | -670.99 | -692.87 | -444.32 | 15790.67 |
| 2 | 306.32 | -13.55 | -670.99 | -692.87 | 1017.44 | 87353.37 |
| 3 | 306.32 | -13.55 | -670.99 | 735.02 | -444.32 | 4112.94 |
| 4 | 306.32 | -13.55 | -670.99 | 735.02 | 1017.44 | 22752.6 |
| 5 | 306.32 | 8.45 | -670.99 | -692.87 | -444.32 | 5644.43 |
| 6 | 306.32 | 8.45 | -670.99 | -692.87 | 1017.44 | 31224.76 |
| 7 | 306.32 | 8.45 | -670.99 | 735.02 | -444.32 | 1470.18 |
| 8 | 306.32 | 8.45 | -670.99 | 735.02 | 1017.44 | 8132.99 |
根据公式1可以计算出该载荷条件下，主轴承所承受的轴向载荷Fa、径向载荷Fr计算结果如表2所示：
表2：轴向载荷Fa、径向载荷Fr
| 载荷 | 工况1 | 工况2 | 工况3 | 工况4 | 工况5 | 工况6 | 工况7 | 工况8 |
| F<sub>a</sub>(kN) | 335.36 | 335.36 | 335.36 | 335.36 | 335.36 | 335.36 | 335.36 | 335.36 |
| F<sub>r</sub>(kN) | 1535.84 | 1580 | 1018.43 | 1083.88 | 1539.93 | 1570.47 | 1024.59 | 1069.94 |
最后可以得到载荷作用下，滚动体的最大接触应力，如表3所示：
表3：各工况载荷下滚动体最大接触应力
根据ISO281:2007标准，轴承的基本额定寿命L10＝304023h，根据ISO16281标准，轴承的修正额定参考寿命L10mr＝4507237h。
根据GL标准，疲劳工况下，主轴承的校核需满足以下条件：最大接触应力ph≤1650Mpa，L10≥130000h，L10mr≥1750000h，通过校核结果发现，该轴承的校核不通过。
本发明还公开了一种风电机组主轴承力学特性的计算系统，包括
第一模块，用于根据风轮载荷、传动链布局参数、仰角α，主轴重量Ms，齿轮箱重量Mg，求解作用于主轴承上的轴向载荷Fa和径向载荷Fr；
第二模块，用于以轴向载荷Fa和径向载荷Fr作为输入量，分析主轴承内圈、滚动体的位移，计算滚动体的弹性趋近量δ的表达式，建立内圈的平衡方程；
第三模块，用于模拟套圈与滚动体的接触，建立滚动体接触载荷Q与弹性趋近量δ、最大接触应力ph的方程；
第四模块，用于联立上述方程，运用迭代法进行各参数求解。
本发明进一步公开了一种计算机可读存储介质，其上储存有计算机程序，所述计算机程序在被处理器运行时执行如上所述的风电机组主轴承力学特性的计算方法的步骤。同时，本发明还公开一种计算机设备，包括存储器和处理器，所述存储器上存储有计算机程序，所述计算机程序在被处理器运行时执行如上所述的风电机组主轴承力学特性的计算方法的步骤。本发明实现上述实施例方法中的全部或部分流程，也可以通过计算机程序来指令相关的硬件来完成，计算机程序可存储于一个计算机可读存储介质中，该计算机程序在被处理器执行时，可实现上述各个方法实施例的步骤。其中，计算机程序包括计算机程序代码，计算机程序代码可以为源代码形式、对象代码形式、可执行文件或某些中间形式等。计算机可读介质可以包括：能够携带计算机程序代码的任何实体或装置、记录介质、U盘、移动硬盘、磁碟、光盘、计算机存储器、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、电载波信号、电信信号以及软件分发介质等。存储器可用于存储计算机程序和/或模块，处理器通过运行或执行存储在存储器内的计算机程序和/或模块，以及调用存储在存储器内的数据，实现各种功能。存储器可以包括高速随机存取存储器，还可以包括非易失性存储器，例如硬盘、内存、插接式硬盘，智能存储卡(Smart MediaCard,SMC)，安全数字(Secure Digital,SD)卡，闪存卡(Flash Card)、至少一个磁盘存储器件、闪存器件、或其它易失性固态存储器件等。
虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明。任何熟悉本领域的技术人员，在不脱离本发明技术方案范围的情况下，都可利用上述揭示的技术内容对本发明技术方案做出许多可能的变动和修饰，或修改为等同变化的等效实施例。因此，凡是未脱离本发明技术方案的内容，依据本发明技术实质对以上实施例所做的任何简单修改、等同变化及修饰，均应落在本发明技术方案保护的范围内。
Claims (10)
1.一种风电机组主轴承力学特性的计算方法，其特征在于，包括步骤：
S01、根据风轮载荷、传动链布局参数、仰角α，主轴重量Ms，齿轮箱重量Mg，求解作用于主轴承上的轴向载荷Fa和径向载荷Fr；
S02、以轴向载荷Fa和径向载荷Fr作为输入量，分析主轴承内圈、滚动体的位移，计算滚动体的弹性趋近量δ的表达式，建立内圈的平衡方程；
S03、模拟套圈与滚动体的接触，建立滚动体接触载荷Q与弹性趋近量δ、最大接触应力ph的方程；
S04、联立步骤S02和步骤S03中的方程，运用迭代法进行各参数求解。
2.根据权利要求1所述的风电机组主轴承力学特性的计算方法，其特征在于，在步骤S01中，主轴承上的轴向载荷Fa和径向载荷Fr的表达式为：
式中：[Mx，My，Mz，Fx，Fy，Fz]为风轮载荷；坐标系采用GL坐标系；A、B、C、D为传动链布局参数。
3.根据权利要求2所述的风电机组主轴承力学特性的计算方法，其特征在于，在步骤S02中，滚动体的弹性趋近量δ的表达式为：
式中：δ1j，δ2j分别为第1，2列滚子弹性趋近量；A0为轴承内、外圈滚道曲率中心初始距离；α0为初始接触角；
为第j个滚子方位角；Dw为滚子直径，cd为双列球面滚子轴承的径向工作游隙。
4.根据权利要求3所述的风电机组主轴承力学特性的计算方法，其特征在于，在步骤S02中，所述内圈平衡方程表达式如下：
式中：Z为滚子数量；α1j、α2j分别为第j个滚子受载后的接触角，其中
5.根据权利要求4所述的风电机组主轴承力学特性的计算方法，其特征在于，在步骤S03中，滚动体接触载荷Q与弹性趋近量δ、最大接触应力p的关系式为：
其中：
式中：bk为两切片间距长度的一半，N为切片的总数量；z(yl)为滚子第l个切片处与滚道表面的初始间距，Dl,k为柔度系数。
6.根据权利要求5所述的风电机组主轴承力学特性的计算方法，其特征在于，步骤S04的具体过程为：
根据式(1)，求解出轴承的轴向载荷Fa和径向载荷Fr，作为后续方程组的输入；
先赋初值δa、δr，根据式(2)、(3)、(7)、(8)、(9)求解两列滚动体的弹性趋近量δ，再对接触区域上各切片的接触半宽ak赋初值，求解柔度系数Dl,k，运用式(12)的变形式p＝πE’×(D-1*W)得到滚动体上的接触应力，其中W为每个切片中心的弹性变形矩阵，可根据各方位角滚子与内、外圈之间的弹性趋近量及滚子与滚道之间的凸度来计算；
然后根据有限长一维接触问题的前提条件及Hertz线接触理论推导出各切片的接触半宽计算公式ak＝2Rpkmax/E’，其中R为滚子-滚道接触体在切片处的综合曲率半径；通过与接触半宽ak的初值或上次的计算值进行反复迭代，完成给定条件下滚子-滚道之间的接触应力求解p；
根据式(11)求解Q，再根据已求出的弹性趋近量δ求解轴向载荷Fa和径向载荷Fr，通过对比利用式(1)求解出来的Fa、Fr，进行反复迭代，直到满足精度要求；
求解完成后可输出主轴承接触应力关键评价指标及寿命，对比预设标准应力和预设标准寿命，可完成主轴承的校核。
7.根据权利要求1至6中任意一项所述的风电机组主轴承力学特性的计算方法，其特征在于，在步骤S03中，采用一维有限长接触理论模拟套圈与滚动体的接触。
8.一种风电机组主轴承力学特性的计算系统，其特征在于，包括
第一模块，用于根据风轮载荷、传动链布局参数、仰角α，主轴重量Ms，齿轮箱重量Mg，求解作用于主轴承上的轴向载荷Fa和径向载荷Fr；
第二模块，用于以轴向载荷Fa和径向载荷Fr作为输入量，分析主轴承内圈、滚动体的位移，计算滚动体的弹性趋近量δ的表达式，建立内圈的平衡方程；
第三模块，用于模拟套圈与滚动体的接触，建立滚动体接触载荷Q与弹性趋近量δ、最大接触应力ph的方程；
第四模块，用于联立各方程，运用迭代法进行各参数求解。
9.一种计算机可读存储介质，其上储存有计算机程序，其特征在于，所述计算机程序在被处理器运行时执行如权利要求1-7中任意一项所述的风电机组主轴承力学特性的计算方法的步骤。
10.一种计算机设备，包括存储器和处理器，所述存储器上存储有计算机程序，其特征在于，所述计算机程序在被处理器运行时执行如权利要求1-7中任意一项所述的风电机组主轴承力学特性的计算方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911182046.9A CN112861271B (zh) | 2019-11-27 | 2019-11-27 | 风电机组主轴承力学特性的计算方法、系统及介质 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911182046.9A CN112861271B (zh) | 2019-11-27 | 2019-11-27 | 风电机组主轴承力学特性的计算方法、系统及介质 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112861271A true CN112861271A (zh) | 2021-05-28 |
| CN112861271B CN112861271B (zh) | 2024-04-16 |
Family
ID=75984667
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911182046.9A Active CN112861271B (zh) | 2019-11-27 | 2019-11-27 | 风电机组主轴承力学特性的计算方法、系统及介质 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112861271B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113496091A (zh) * | 2021-09-07 | 2021-10-12 | 西安航天动力研究所 | 液体火箭发动机高速重载球轴承接触动态特性仿真方法 |
| CN113883022A (zh) * | 2021-10-29 | 2022-01-04 | 新疆金风科技股份有限公司 | 风力发电机组轴系和风力发电机组 |
| CN114624021A (zh) * | 2022-02-22 | 2022-06-14 | 东风汽车集团股份有限公司 | 试验台、变速器锁止机构扭转疲劳试验方法及相关设备 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015133582A1 (ja) * | 2014-03-06 | 2015-09-11 | Ntn株式会社 | 計算方法、計算装置、計算システム、アンギュラ玉軸受および機械装置 |
| CN107729597A (zh) * | 2017-08-29 | 2018-02-23 | 明阳智慧能源集团股份公司 | 一种主轴轴承滚道校核工具 |
-
2019
- 2019-11-27 CN CN201911182046.9A patent/CN112861271B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2015133582A1 (ja) * | 2014-03-06 | 2015-09-11 | Ntn株式会社 | 計算方法、計算装置、計算システム、アンギュラ玉軸受および機械装置 |
| CN107729597A (zh) * | 2017-08-29 | 2018-02-23 | 明阳智慧能源集团股份公司 | 一种主轴轴承滚道校核工具 |
Non-Patent Citations (1)
| Title |
|---|
| 付大鹏;褚加瑞;: "基于Romax的风电机组主轴轴承接触应力分析", 东北电力大学学报, no. 06 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113496091A (zh) * | 2021-09-07 | 2021-10-12 | 西安航天动力研究所 | 液体火箭发动机高速重载球轴承接触动态特性仿真方法 |
| CN113496091B (zh) * | 2021-09-07 | 2021-12-24 | 西安航天动力研究所 | 液体火箭发动机高速重载球轴承接触动态特性仿真方法 |
| CN113883022A (zh) * | 2021-10-29 | 2022-01-04 | 新疆金风科技股份有限公司 | 风力发电机组轴系和风力发电机组 |
| CN114624021A (zh) * | 2022-02-22 | 2022-06-14 | 东风汽车集团股份有限公司 | 试验台、变速器锁止机构扭转疲劳试验方法及相关设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112861271B (zh) | 2024-04-16 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112861271A (zh) | 风电机组主轴承力学特性的计算方法、系统及介质 | |
| Chen et al. | Load performance of large-scale rolling bearings with supporting structure in wind turbines | |
| CN107704652A (zh) | 一种风力发电机组轴承刚度计算工具 | |
| CN105160067A (zh) | 风机变桨轴承与轮毂和叶片连接螺栓强度性能的校核方法 | |
| CN111444647B (zh) | 一种风机变桨轴承连接螺栓建模与强度校核方法 | |
| Waghole et al. | Optimization of needle roller bearing design using novel hybrid methods | |
| CN105550403B (zh) | 一种风机外购件强度校核工具 | |
| CN107590356B (zh) | 一种风力发电机组主轴轴承自动选型方法及存储设备 | |
| CN107729597A (zh) | 一种主轴轴承滚道校核工具 | |
| CN105022868A (zh) | 风机偏航轴承连接螺栓极限强度和疲劳强度的计算方法 | |
| CN106503367A (zh) | 一种风力发电机组主机架强度计算方法及装置 | |
| CN113515865A (zh) | 一种风电齿轮箱行星级均载性能优化方法 | |
| Graßmann et al. | Validation of a finite-element model of a wind turbine blade bearing | |
| CN112035979A (zh) | 一种风电机组整机一体化优化设计方法、平台及寻优方法 | |
| CN108256704A (zh) | 风力发电机的子系统动态特征的仿真方法及仿真设备 | |
| Behnke et al. | Exploring limiting factors of wear in pitch bearings of wind turbines with real scale tests | |
| Aguirrebeitia et al. | Wire twisting stiffness modelling with application in wire race ball bearings. Derivation of analytical formula and finite element validation | |
| CN108595765A (zh) | 风力发电机组双列圆锥滚子轴承载荷分布及寿命计算方法 | |
| Song et al. | Enhanced test strategy of pitch bearing based on detailed motion profile | |
| CN110598241B (zh) | 一种风机轮毂强度校核方法 | |
| CN112364551A (zh) | 一种行星传动齿轮箱疲劳分析方法 | |
| CN107563041A (zh) | 一种风电机组大部件静强度快速评估方法 | |
| CN114021288B (zh) | 一种风电机组偏航轴承寿命预测的方法 | |
| Jin et al. | Structural optimization of pitch bearings for offshore wind turbines | |
| Menck | The Finite Segment Method—A Numerical Rolling Contact Fatigue Life Model for Bearings Subjected to Stochastic Operating Conditions |
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 |