风力发电机塔筒顶部法兰的有限元分析.docx

1、风力发电机塔筒顶部法兰的有限元分析风力发电机塔筒顶部法兰的有限元分析李杰摘要:使用非线性有限元软件 MSC. Marc /Mentat建立风力发电机塔筒顶部法兰联接的有限元模型 ,施加合理的边界条件和载荷后 ,通过非线性接触分析得到在预紧工况和极限工况下各组件的应力分布和变化情况 ,并对塔筒顶部法兰接触面进行安全性校核和塔筒顶部法兰与塔筒筒体焊缝处的疲劳寿命分析。分析结果表明 ,该塔筒顶部法兰的强度、安全性和焊缝处的疲劳寿命均满足设计要求 ,且结果为大型风力发电机法兰合理设计和性能强化提供了科学依据。关键词:有限元法;塔筒顶部法兰;非线性接触;风力发电机Finite element analy

2、sis of wind turbine tower top flangeAbstract: Using nonlinear finite element software MSC. Marc/Mentat establish a FEA model of wind turbine top flange connection. Based on applying appropriate boundary condition and loads, the stress distribution and the stress variation of each component were obta

3、ined with the method of nonlinear contact analysis under the case of pre-tighten work condition and the limited condition. The analysis results indicate that the strength and security of flange and the fatigue life fillet weld meet design requirements. The result has a guiding role for rational desi

4、gn and performance improvement of large-scale wind turbine flange.Key words: finite element method; tower top flange; nonlinear contact; wind turbine塔筒顶部法兰是风力发电机运行时的主要受力部件之一 ,它与偏航刹车盘、偏航轴承通过螺栓联接。在风力发电机组工作时 ,作用于叶轮和机舱的载荷将通过偏航轴承传递到塔架上。这些载荷包括叶轮所受到的气动载荷 ,如升力、阻力、轴向推力、气动弯矩和气动转矩等。由于风速风向的不稳定性 ,机组运行时塔筒顶部法兰承受的载

5、荷是动态变化的 ,很容易遭到破坏。为提高风力发电机运行的可靠性 ,需要对塔筒顶部法兰联接部分进行结构强度分析。本文采用整体建模的思路 ,利用有限元理论分析塔筒顶部法兰面的接触应力 ,并且对塔筒顶部法兰与塔筒筒体的焊缝连接进行疲劳寿命分析。1整体模型的接触分析1.1几何模型几何模型在NX/UG中建立 ,本文重点分析塔筒顶部法兰的接触 ,所以简化轴承内圈 ,建立偏航轴承外圈与偏航刹车盘和塔筒顶部法兰的螺栓联接模型。如图 1所示。1.塔筒 (塔顶连接法兰 ) 2.偏航刹车盘 3.偏航轴承外圈 4.螺栓1.2有限元模型在 MSC. Marc/Mentat中建立塔筒顶部法兰接触的有限元模型 ,如图 2所

6、示。模型均采用六面体单元划分 ,为了减小结构的关键部位对计算结果的影响 ,在塔筒顶部法兰圆角处的网格适度加密 ,单元总数为76104,节点总数为 113486。模型包括偏航轴承外圈、偏航刹车盘、塔筒顶部法兰和 84个螺栓等可变形接触体 ,各部件所用的材料均为低合金高强度结构钢,其弹性模量为 2.06105MPa,泊松比为 0.3,密度为 7.85 10-9t/mm3,根据实际的装配关系 ,建立如表 1所示的接触关系。表1中: T为接触 (touching)； G为粘接 (glue) ；括号中数值为接触对的摩擦因数。螺栓与偏航轴承外圈在螺纹联接处有 55mm的螺纹啮合 ,本模型中通过粘接 (gl

7、ue)关系来耦合接触双方的自由度 ,以模拟螺纹联接。 图2 塔筒顶部法兰接触的有限元模型1.3载荷设置和边界条件接触分析所采用的载荷来自风力发电机极限工况载荷数据。分析分两种工况进行:第一种工况中给 84个螺栓施加预紧力 F0 ( preload _ bolt)254kN ,建立各部件间的接触关系；第二种工况中撤掉预紧力 ,锁定螺栓的长度 (lock_bolt) ,然后施加偏航轴承极限载荷 (load_bearing)和 偏航刹车盘的极限载荷(load_brake) ,分析各零 /部件的应力分布情况。为了准确模拟载荷的传递 ,在偏航刹车盘中心点处施加极限载荷报告中此处载荷 MZ的最大值为-32

8、85kNm ,通过MSC. Marc/Mentat软件施加载荷单元RBE3,将载荷分布到偏航刹车盘上七个与刹车器接触的区域 (如图 3所示),模拟偏航刹车。在轴承外圈中心点与轴承外圈内表面建立 RBE3连接 ,施加偏航轴承处的极限载荷为: FX = - 444.2kNFZ = - 1350.8kNMXY =- 4570.7kN m图3 刹车盘的结构约束塔筒底部端面所有节点 X、Y、Z 方向的平动自由度 (fix _ tower) , 以消除接触分析时的刚体位移。同时考虑重力 (gravity)的影响。所有载荷设置和边界条件的施加如图 4所示。图4 载荷设置和边界条件(二分之一模型)1.4接触分

9、析结果螺栓预紧后 ,各接触面接触良好 ,各零/部件应力值均在屈服强度之下。施加极限载荷后 ,各接触面间仍保持很好的接触状态 ,各零/部件的应力值稍有增大 ,下面给出只施加预紧力和极限载荷后的各零/部件接触应力云图。偏航轴承外圈的最大应力值为 191.2MPa,出现在螺栓联接的第一圈螺纹处(见图 5) ,小于轴承材料的屈服极限 650MPa,轴承强度满足要求。图5偏航轴承外圈的应力云图(极限载荷工况下)由于接触面比较小 ,塔筒顶部法兰的最大应力值在预紧工况下为 293.2MPa(见图6)；施加极限载荷后为294.7MPa(见图7) ,均出现在塔筒顶部法兰面与螺栓接触面周围 ,且都小于塔筒顶部法兰

10、材料的屈服极限325MPa,塔筒顶部法兰强度满足要求。图6 塔筒顶部法兰的应力云图(预紧工况下)图7 塔筒顶部法兰的应力云图(极限载荷工况下)塔筒顶部法兰联接螺栓的最大应力在预紧工况下为496.4MPa (见图8) ,施加极限载荷后为499.9MPa(见图9) ,均出现在螺栓的中段 ,且都小于螺栓的屈服极限900MPa (该螺栓为10.9级) ,所以螺栓是安全的。图8 螺栓的应力云图 (预紧工况下 )图9 螺栓的应力云图 (极限载荷工况下 )由于塔筒顶部法兰圆角处通常会出现应力集中 ,所以对圆角进行安全校核。塔筒顶部法兰圆角处最大等效应力值为102.8MPa,如图 10所示。图10 塔筒顶部法

11、兰圆角处的应力云图塔筒顶部法兰圆角处安全裕度 Msult为:Msult =/m-1 (1) =s/n (2) 式中:m 为极限载荷下塔筒顶部法兰圆角处的等效应力,m = 102.8MPa；s为塔筒顶部法兰材料屈服极限,s = 325MPa；n为安全系数, n取 1.1。 则: Msult = 1.8741塔筒顶部法兰圆角处的安全裕度大于零, 因此塔筒顶部法兰圆角处静强度符合要求。2塔筒顶部法兰接触面的安全性校核由于受弯矩作用, 除了要求静强度足够外, 还应保证接合面间不出现缝隙和接触面不被压溃。由于预紧力作用,接触面上产生的压应力F0 为:F0 =zF0/A=84254/0.86710-3=

12、24.62MPa式中: z为螺栓个数；A为接触面的面积；F0为螺栓预紧力。由于弯矩的作用,接触面上产生的弯应力M 为:M =MXY/W=4570.7/0.599710-3= 7.621MPa式中:MXY为接触面处的弯矩；W 为抗弯截面系数。压力FZ在接触面上产生的压应力Z为:Z =FZ/A=- 1350.8/0 .86710-3= -1.559MPa接触面受拉侧不出现缝隙的条件为:l =F0-M +Z 0 (3)式中:l为接触面所受拉应力。接触面受压侧不被压溃的条件为:Y =F0 +M +Z p (4) 式中:Y为接触面所受压应力；p 为被连接材料的许用挤压应力, p = 0.8s = 276

13、MPa。由此可得:l = 15.44MPaY = 30.68MPa pl大于零,塔筒顶部法兰接触面不会产生缝隙。Y小于被连接材料的许用挤压应力,所以塔筒顶部法兰不会被压溃。3塔筒顶部法兰焊缝疲劳寿命计算塔筒顶部法兰和塔筒筒体通过焊接连接,是结构中比较薄弱的环节,焊接结构的耐久性很大程度上由焊缝的疲劳强度所决定。所以,如果能在设计初期对焊缝的疲劳寿命进行预测,以了解整个焊接结构上焊缝疲劳强度的分布情况,就可以对实际中焊缝的处理方式进行合理地调整,以提高产品的整体疲劳性能。3.1焊缝等效疲劳应力计算法兰与塔筒筒体之间采用等厚度焊接, 因此计算时不考虑焊缝应力集中系数的影响。焊缝等效疲劳应力采用有限

14、元方法通过加载相应位置的等效疲劳载荷而获得。具体有限元计算过程如下。3.1.1模型描述使用接触分析的模型。第二种工况(载荷施加工况)的载荷由极限载荷改为等效疲劳载荷Fx = -142.62kN , Fz = 38kN , My = - 1592kNm , Mz =1459kNm。3.1.2有限元计算结果计算得到法兰焊缝处最大等效应力为 22.95MPa,如图11所示,最大等效应力位置的节点坐标为(- 1460.95, - 54.42, - 190.324) 。图11 塔筒顶部法兰焊缝处的应力云图3.2确定焊缝材料的疲劳寿命S-N曲线焊缝材料S-N曲线如图 12所示。图12 S-N曲线S-N曲线

15、(以循环基数 ND为界):N=ND(D/w)3 (ND) (5) N=ND(D/w)4 (NND ,w D) (6) 式中: N为循环次数；D 为循环基数ND对应的疲劳应力；w为计算所得应力。3.3焊缝疲劳应力循环次数计算有限元方法计算得到的焊逢处最大等效疲劳应力ph = 22.95MPa。因此焊逢处等效疲劳应力 ph为:ph =mph (7) 式中:m 为材料偏安全系数,取1.1。计算得:ph 为 25.245MPa。已知塔筒筒体材料的屈服极限 s = 470MPa,强度极限b = 345MPa,故由疲劳极限估算公式得焊缝疲劳极限A 为:A = 0.23(b +s) (8) 可求出焊缝疲劳极限A = 187MPa,对应的循环次数NA为 1107 。已知循环基数 ND = 5E6,循环次数大于N D 时用式(9)计算w 为:w =A(NA/N)1/4 (9) 所以:D =A(NA/ND)1/4= 1.19A根据 GL规范 ,存活率大于 97.7%的 S-N 曲线需要进行修正 ,缩减为原来的 2/3 。因此修正循环基数