摘 要：针对风力发电机组偏航轴承结构与受力的复杂性，提出了一种新的偏航轴承简化建模方法。简要阐述了某MW级风电机组偏航轴承内圈与主机架连接螺栓、外圈与塔筒连接螺栓的强度，分析模型利用GAP单元模拟偏航轴承滚珠与内外圈之间的接触来实现力的传递。结果表明，GAP单元可以有效模拟轴承实际受力和力的传递，为偏航轴承建模提供一种有效的简化方法。
关键词：风力发电机组；GAP单元；偏航轴承；有限元
中图分类号：TM315；TH133.3 文献标识码：A DOI：10.15913/j.cnki.kjycx.2017.05.016
轴承是风电机组的重要组成部件，主要包括主轴轴承、变桨轴承、偏航轴承、发电机轴承和齿轮箱轴承。其中，偏航轴承通过螺栓分别与主机架和塔筒连接，维持偏航系统的稳定运转，对整个机组的受力特性、稳定性、震动和噪声等都起着至关重要的作用。当风向改变时，偏航轴承旋转使风轮正对风向，以捕获更多的风能。
偏航轴承是四点接触球轴承，尺寸庞大且结构复杂，当计算与偏航轴承相连零部件的强度时，如果对偏航轴承采用实体建模，计算量比较大，且计算过程不易收敛。这就需要寻求一种简化偏航轴承的有效方法，既能提高计算效率，同时，又不影响与之相连零部件的强度计算结果。
文章利用有限元分析软件MSC.Marc/Mentat中的GAP单元代替滚动体，利用模拟滚动体与内外圈之间的接触特性来实现偏航轴承的简化和力的传递。同时，要对模型进行加载，计算轴承内外圈连接螺栓的强度。计算出的强度结果比较理想，可以为后期偏航轴承建模提供一种新的思路。
1 GAP单元的属性特征
1.1 GAP单元的模式
GAP单元有3种模式：真实距离模式，用于几何非线性分析；方向固定模式，用于接触体不穿透给定平面的线性或非线性几何分析中；第三种模式通过用户子程序GAPU中指定閉合距离和间隙方向，间隙方向和距离可在分析过程中被更新，以模拟沿曲面的滑动。
1.2 GAP单元特性描述
假设GAP单元由A，B两点组成，A，B为内外圈与滚动体接触的对应节点对，A，B两点之间的初始间隙设为滚动体的直径Dw。在外力作用下，GAP单元的载荷-位移关系可用接触面的法向分量和切向分量来描述。在未接触区，它不影响分析对象的运动状态；在接触区域，GAP单元的法向刚度将变得足够大，以阻止接触体相互侵入。
当点对的相对位移小于初始间隙时，表明接触表面处于接触状态，GAP单元的法线方向就会存在法向力Fn，且法向力为负值。此时，GAP单元犹如一个线性弹簧，其法向接触刚度为KA，内外圈就会通过GAP单元在A，B两点之间传递载荷；否则相对位移大于初始间隙时，接触点对分离，接触表面为未接触状态，法向力Fn值为零，GAP单元不受力也不传递载荷，不会影响分析对象的运动状态。
在切向方向，当法向力Fn<0且切向力Ft≤μFn（μ为接触面间摩擦系数）时，接触面保持接触且无相对滑动；当Fn<0且Ft>μFn时，接触面间将产生相对滑动。
2 偏航轴承连接部件的有限元建模
2.1 偏航系统计算各部套连接三维模型
在Pro/E中建立偏航系统计算各部套的三维模型，主要包含定子、主机架、偏航轴承、刹车盘和塔筒，省略了偏航轴承中保持架的作用。螺栓将在后面的有限元模型中建立，如图1所示。主机架通过螺栓84-M30与偏航轴承内圈连接，螺纹啮合长度为45 mm；塔筒与刹车盘也通过螺栓84-M30与偏航轴承外圈连接，螺纹啮合长度为55 mm。
2.2 有限元模型
风力机偏航轴承是四点接触球轴承，在受载过程中，载荷传递路线经过3个关键点，即滚珠与外圈接触点、滚珠中心点以及滚珠与内圈接触点，载荷是以直线的形式传递的，如图2所示。
根据四点接触球轴承的传力特点，采用GAP单元代替滚动体模拟轴承内外圈的接触和力的传递。GAP单元采用真实距离模式，在接触体之间建立接触点对，不考虑滚子变形，将GAP单元设置为刚体。初始接触点处的GAP单元设置为初始接触，其余GAP单元初始设置为非接触，GAP单元长度为滚子直径。模型中，GAP单元过多就不易收敛，此例中滚动体与内外圈接触的点对有14对，所以，采用14个GAP单元代替一个滚动体进行接触应力分析，所得模型收敛性好。
在有限元分析软件MSC.Marc/Mentat中建立有限元模型，模型采用8节点六面体7号单元划分，单元大小为10～30 mm，单元数总共112 902个，如图3和图4所示。
3 计算模型的接触设置
模型中只需设置7个可变形接触体，包括主机架（main_con）、偏航轴承内圈（ph_in）、内圈上螺栓（ph_in_bolt）、偏航轴承外圈（ph_out）、刹车盘（brake）、塔筒（tower）和外圈上螺栓（ph_out_bolt）。内圈与主机架的接触面、内圈上螺栓与内圈的接触面、外圈上螺栓与外圈的接触面、刹车器与外圈的接触面、刹车器与塔筒的接触面均设置为touching，摩擦系数取0.15.
这里主要分析偏航轴承及其连接螺栓的强度，因为GAP单元对其影响最为显著。不考虑定子与主机架连接螺栓，将定子与主机架接触面通过粘接（Glue）关系来耦合。
内圈螺栓与偏航轴承内圈在螺纹联接处有45 mm的螺纹啮合，外圈螺栓与偏航轴承外圈在螺纹联接处有55 mm的螺纹啮合，螺纹啮合部位通过粘接（Glue）关系来耦合两部分的自由度，以模拟螺纹联接。
4 模型的边界条件及载荷施加
边界条件设置为：约束塔筒底部6个方向的自由度，约束轴承内圈绕其轴向的转动自由度，约束螺栓接触体的径向位移自由度（消除计算中多余的刚体位移），锁定螺栓轴向位移（预紧过后锁定螺栓长度）。在轮毂中心点处建立一MPC点与定子外表面相连，将轮毂中心处载荷添加到MPC点上，如图5所示。
分析过程按2个工况进行，第一个工况中给所有螺栓施加预紧力，预紧力取350 kN，建立各部件间的接触关系；第二个工况中施加轮毂中心处的极限载荷（由BLADED软件计算得到，共16种工况载荷），求得各螺栓的强度。
5 有限元分析结果
5.1 螺栓与偏航轴承应力结果
采用MSC.Marc2005r2求解器计算各极限工况下螺栓与偏航轴承的应力结果。其中，最大工况下的应力云图如图6、图7、图8、图9、图10所示。
由应力结果可知，偏航轴承内圈连接螺栓在轮毂极限载荷下的最大应力值为574.5 MPa，偏航轴承外圈连接螺栓在轮毂极限载荷下的最大应力值为575.4 MPa，偏航轴承在轮毂极限载荷下的最大应力值为597.6 MPa。
5.2 螺栓与偏航轴承安全裕度计算
因螺栓与偏航轴承的屈服极限分别为900 MPa和930 MPa，所以，可计算出内圈螺栓的安全裕度为0.424，外圈螺栓的安全裕度为0.422，偏航轴承的安全裕度为0.415.
螺栓和偏航轴承的安全裕度均大于零，说明零件强度满足要求。
6 结论
文章考虑偏航轴承的结构特点和接触方式，采用GAP单元模拟滚珠来简化偏航轴承。根据风电机组的实际受力情况对模型进行加载计算，得出螺栓与偏航轴承的应力结果。结果显示，螺栓强度满足要求，与实际情况相符。GAP单元对偏航轴承的简化建模可以有效地模拟轴承的实际受力和力的传递，为风电机组偏航轴承建模提供一种新的简化方法。
参考文献
[1]杨静，廖辉.风力机组偏航变桨轴承载荷分布分析及数值求解[J].东方汽轮机，2012（3）.
[2]陈火红，杨剑，薛小香，等.新编Marc有限元实例教程[M].北京：机械工业出版社，2007.
[3]陈达亮，郝志勇.GAP单元在气门弹簧动态特性有限元分析中的应用[J].拖拉机与农用运输车，2004（2）.
[4]皮亚南，董懿琼，戴莉莉，等.滚动轴承保持架不稳定性分析[J].江西科学，2008，10（5）.
[5]闻邦椿.机械设计手册[M].北京：机械工业出版社，2010.
〔编辑：白洁〕
