风电偏航轴承寿命计算与承载能力分析

本资料包含pdf文件1个，下载需要1积分

Life Calculation and Carrying Capacity Analysis of Yaw Bearingsin W ind TurbinesFan Zeng-hui，Li Xiu-zhen，Gu Xiao-hui，Ouyang Hua，Zhu Wei-liang(CSR Zhuzhou Institute Co.，Ltd.，Zhuzhou 412001，China)Abstract：Based on comprehensive consideration of the basic design concept for bal beatings and rotary supporting，combining with GL design specification and ISO standard of wind turbine industry，the related calculation methods isreasonably used based on real loading conditions under yaw bearings coordinate system of wind turbines.The concreteprocess is accurately combed for calculation of life and carrying capacity of yaw bearings.The diference is compared be-tween the different standards and specifications.The results show that despite the calculation results of various methodshave differences，but does not affect the basic judgment for yaw bearings，SO the simple algorithm is able to be used torapidly check and analyze yaw bearings。

Key words：wind turbine；yaw beating；life calculation；carying capacity；check风力发电机偏航轴承作为风力发电机(以下简称风机)迎风捕获能量的回转支承，安装于风机机舱底部，是偏航系统的关键部件。由于风机吊装维护费用极高，且风机工况复杂，偏航轴承-般承受较大的轴向力、径向力和倾覆力矩的联合作用，因此对偏航轴承的承载能力、寿命及可靠性提出了很高的要求。

目前主流的兆瓦级以上风机的偏航轴承主要采用四点接触球转盘轴承，因此对风机偏航轴承的设计计算大多运用经典的球轴承设计理论或者回转支承的设计方法 。但由于风机工况的特收稿 日期：2012-12-21；修回日期：2013-02-16作者简介：凡增辉(1986-)，男，工程师，主要从事风机机械传动设计工作，E-mail：fanzh###teg.en。

殊性与复杂性，风机偏航轴承的设计计算实际上有其独特的地方。下文综合运用球轴承的基本原理与风机设计的相关标准与规范，在基于风机真实工况与载荷的前提下，对风机偏航轴承的疲劳寿命与承载能力进行分析计算。

1 载荷工况偏航轴承外圈与偏航制动盘-起固定在塔简上，偏航轴承内圈与主框架通过螺栓连接(图 1)。

当偏航电动机工作时，带动主框架随着偏航轴承内圈旋转，从而使整个机舱和风轮处于最佳迎风状态。

在风机设计时，参照 GL规范与 IEC标准 ]，运用行业通用的风机设计软件 Bladed，对风机的各工况进行模拟计算，经过后处理，可以提取出针对《轴承)2013.No.7偏航轴承坐标系下的载荷工况，如图2所示。

l-偏航轴承；2-塔筒上凸缘；3 滞I动盘，4 制动器 机架；6-偏航驱动图1 偏航轴承安装示意图图2 偏航轴承坐标系示意图图2中偏航轴承坐标系下的6个载荷分量依据偏航轴承的载荷特点可以简化合成为 3个载荷：与 F 相等的沿z轴的轴向力 F ；由F ，F 合成的径向力 F ；由 ， 合成的倾覆力矩 ，如图1所示(偏航轴承不承受 ，因此，在此不予考虑)。

风机偏航轴承的承载能力与寿命主要考虑对象是沟道，且其承受的主要载荷为倾覆力矩，因此需要在 Bladed软件中利用Channel Combination”提取不同倾覆力矩 的累积转数(或时间)分布统计，以计算偏航轴承全寿命周期内的等效疲劳载荷，如图3所示。为了使统计数据更保守，轴向力F 和径向力 F 固定为最大倾覆力矩时对应的值。

0 1 000 2 000 3 000 4 000 5 000倾覆力矩/(kN-m)图 3 偏航 轴承倾覆 力矩值 的累积 时间分布而相对于疲劳载荷的统计处理，极限载荷的提取相对来说比较简单，即从 Bladed后处理程序中提取各种工况下全寿命内最大 对应的载荷即可，见表 1。其中载荷工况中各工况具体描述见文献[8]。

表 1 偏航轴承中心极限载荷2 寿命计算2.1 理论分析 型 。在该模型中，寿命的绝对值和离散度决定证明轴承拥有足够寿命的多数方法是运用数 Weibul1分布的曲线形状。从 Lundberg和Palmgren咖 咖 咖加 mu/厘害旺 酶凡增辉 ，等：风电偏航轴承寿命计算与承载能力分析的寿命试验结果可知，对于球轴承表达离散度的e值取 10/9，Lundberg和 Palmgren从理论与试验推导出如下关系。c(古) ， (1)式中：P为寿命指数，对于球轴承P3；P为当量动载荷。

对于(1)式，如将比例常数设为 c，则(1)式为 ( ) ， (2)式中： 。为基本额定寿命，10。r(内圈转数)；C为基本额定动载荷 ，在实际应用中， 。寿命用时间(h)表示时的关系式为L 石1o0 6( C) ， (3)式中： 为轴承的转速，r/min。

2.2 基本额定动载荷Lundberg和 Palmgren将 Hertz弹性接触理论与Weibul损伤概率理论相结合，并考虑了轴承的尺寸及应力循环次数的影响，在此基础上依据试验结果修正相关系数与指数，最终推导出了球轴承基本额定动载荷为C 1[1.044( 州·( m ×40.21( 甜·L 上 (icos 0 7方D ， (4)(1 )式中：A为形状误差参数；Ot为接触角；i为球的列数；z为单列球数； ，re分别为内、外圈沟曲率半径； 为轴承结构参数， ：D-w cos a；D 为钢球上，pw直径；D。 为球组节圆直径。

由于(4)式较为复杂，在ISO标准中将其进行了整理 。对于如偏航轴承等球径 D >22.5mln的轴承，C3.647b (/cos ) Z D7～ ， (5)式中：b 为额定系数，对于偏航轴承，b 1.3；厂c为 的函数，可通过文献[12]计算或查表插值得到。ANSI/ABMA标准 中，b 和 被合并为， 从而得到基本额定动载荷的计算式为C3.64 (icos ) z D 。 (6)2.3 当量动载荷Lundberg和Palmgren给出了将联合载荷换算成当量动载荷的计算式及其系数，其核心是通过轴承内部载荷分布与外载荷的静力平衡关系给出等效载荷。由于该关系式过于复杂，ISO等标准采用下列简化公式来计算当量动载荷PXF ， (7)式中： ，y分别为径向、轴向载荷系数，可以依据ISO 281提供的数据查表得到。但如上文所述，偏航轴承主要承受的是倾覆力矩，如果依据(7)式推导当量动载荷，就需要将倾覆力矩转化为偏心轴向载荷。

在将倾覆力矩转化为偏心轴向载荷的方法上，还有基于相对偏心距，利用 Simpson近似积分法通过积分参数获得当量动载荷的计算式P [ xy· ]， (8) J(81. ) 厶 u u式中：厶( ， ( 以)， ( )为载荷计算的积分系数l引。显然，以上计算公式较为繁琐 ，在 NREL提供的风机设计规范 中，提供了更为简便的计算公式，P：0.75F F 。 (9)上，pw3 承载能力在风机设计时，由于 20年设计寿命的考虑，常常需要考虑 10年甚至 50年-遇的极端风况的载荷影响，因此，在通过 Bladed输出的风机载荷中，极限载荷与疲劳载荷的数值差距较大，如图 3和表 1所示。因此，对偏航轴承进行设计校核时需计算其静承载能力。

3.1 静态安全系数轴承的静承载能力以滚动体允许的最大变形时的应力(对于球轴承，该值为4 200 MPa)与其实际承受的最大 Hertz接触应力 的比值来衡量，因此，转盘轴承的承载能力-般用静态安全系数来衡量，表示为：(4 200)。， (10)or ： ， (11) m ---- L ll式中：Q 为偏航轴承最大滚动体载荷；o，b分别为 Hertz接触椭圆的长、短半轴，可以通过计算曲率函数 F(p)获得 J。

依据IEC61400的要求，在进行风机结构设计《轴承)2013.No.7的最大极限状态分析时，还需要考虑材料的不确定性和易变性，分析方法的不确定性以及失效时结构部件的重要性。因此，在针对偏航轴承的承载能力分析时，在(10)式的基础上还需要引入材料局部安全系数 与重要失效系数 进行修正，则： ( )，。 (12)m norma3.2 最大滚动体载荷对于F ，F 及 联合作用的偏航轴承的最大滚动体载荷，可通过下面方法求得。

Jones建立的载荷 -位移关系为QKfi ， (13)式中： 为轴承的载荷 -位移系数； 为接触变形量；f为载荷 -位移指数。

建立以轴承相对轴向、径向与角位移为未知量的静力平衡方程为F ∑Q sin ，∞ 0(14)F ∑Q cos sin 0f， (15)Mx ∑Q COSpsin ， (16)式中：Q。为位置角为 的球的法向载荷。该方程的未知量可以用数值方法如 Newton-Raphson法进行求解。获得轴承相对轴向、径向与角位移后，在 0处即可利用(13)式求得 Q 。

但是上面的计算方法过于繁琐，在 NREL提供的风力设计规范里，提供了非常简单的最大滚动体载荷 Q 的计算公式Q 2Frzgia4 Mxy， (17)式中：F ，F 及 选取表1中 最大的-行数值。

4 工程 示例以某型号2.5 MW风机为例，其偏航轴承结构参数见表2。

表2 偏航轴承结构参数内沟曲率半径/ram外沟曲率半径/mm球径/mm球数球组节圆直径/mm初始接触角31.6731.676o.3251423 05045。

在进行当量动载荷计算时，需要用到图3所示的疲劳载荷(倾覆力矩 )。依据 GL规范的要求，在计算偏航轴承寿命时，需将疲劳载荷依据其作用的累积时间(或转数)分布进行等效，则等效疲劳弯矩 为(18)式中： 为倾覆力矩 的第 i个值；t 为 作用的时间。

利用(18)式计算出等效疲劳弯矩，分别基于(7)～(9)式计算当量动载荷，寿命计算结果见表3。

表3 寿命计算结果根据实际运行统计，偏航轴承运行时间约为风机寿命的 10%，结合偏航的速度，可以估算出在 20年的设计寿命内，偏航轴承需要运行的转数约为55 000 r。从表3可知，所有计算方法的结果均大于该值，判定该偏航轴承寿命符合要求。

该型轴承的承载能力见表4。根据表4可以看出，在考虑了材料与失效安全系数的基础上，该偏航轴承的静态安全系数大于GL规范的要求值1.1，因此该偏航轴承的承载能力满足要求。

表 4 承载能力计算结果参数 数值最大滚动体载荷 Q /kN接触椭圆长半轴a/nun接触椭圆短半轴b/mm最大接触应力 - /MPa静态安全系数，。

5 结束语在综合考虑球轴承以及回转支承基本设计理念的基础上，结合风机行业的GL设计规范与 ISO等标准，基于风机偏航轴承坐标系下的真实载荷工况，合理运用相关的计算方法，得到了完整的计算流程与较为准确的计算结果。

(1)偏航轴承疲劳载荷谱与其寿命计算强相关，需要依据偏航轴承坐标系下不同弯矩作用的累积转数的统计分布，按照 GL规范的要求，将疲劳载荷平均等效。

(2)通过不同计算方法在工程实例中的应用可以看出，各种方法的计算结果虽有差异，但不影(下转第8页)啷蕈i ～· 8· 《轴承)2013.No.7QlO- UlU -0.05 U UI.15 UlU位移/rma图8 轴向力为 150 N时的保持架质心运动轨迹钢球与套圈能够始终保持接触，保证较高的几何耦合，并且也能够减小轴承的轴向及径向振动。

3 结论(1)波纹度数影响保持架质心运动轨迹，随波纹度数增加，保持架质心涡动变得剧烈；(2)波纹度数为奇数时，保持架质心运动轨迹比波纹度数为偶数时的情况混乱；(3)波纹度数相同的情况下，内沟道波纹度对保持架质心运动轨迹的影响大于外沟道；(4)波纹度数-定的情况下，随着轴向力增大，保持架进入稳定涡动时的轨迹更为有序。