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

/风力发电机塔筒顶部法兰的有限元分析李杰摘要：使用非线性有限元软件MSC．Ｍarc／Mｅｎtat建立风力发电机塔筒顶部法兰联接的有限元模型,施加合理的边界条件和载荷后，通过非线性接触分析得到在预紧工况和极限工况下各组件的应力分布和变化情况，并对塔筒顶部法兰接触面进行安全性校核和塔筒顶部法兰与塔筒筒体焊缝处的疲劳寿命分析。分析结果表明,该塔筒顶部法兰的强度、安全性和焊缝处的疲劳寿命均满足设计要求，且结果为大型风力发电机法兰合理设计和性能强化提供了科学依据。关键词:有限元法;塔筒顶部法兰;非线性接触;风力发电机FiｎiteｅlemeｎtanalyｓｉsofwindtｕrbinetowertopflaｎｇeAbstract:UsｉngnonliｎｅaｒｆiniteｅｌementｓoftwareMSC。Mａrc/MｅntatestabliｓhaＦEAmoｄｅlofwinｄtuｒbinetopｆlangeconｎｅctｉoｎ。Ｂaseｄoｎaｐpｌyinｇaｐpｒopｒiａtｅboundaryｃonditioｎanｄlｏads,thestresｓdｉstrｉbutｉonanｄtｈｅstrｅssvariatiｏｎｏfeachｃompoｎｅntｗereobｔａinedwiｔhtheｍethodofnonlｉｎearcｏntactaｎalｙsisundertｈｅcaｓｅｏfprｅ-tighｔeｎwｏｒkcondiｔiｏnａndthelimiｔedｃondition。Theanalysisresuｌtsｉndiｃａtethａtthestreｎgｔhａnｄsｅcuriｔyoｆfｌangeaｎｄｔhefatiｇueliｆefｉlｌetｗeldmeeｔdesｉgｎrequiｒemenｔs。Tｈeresｕltｈasａguidingroｌeforrationalｄesｉgnandpｅrfｏrmａnceimproｖｅｍentｏflarｇｅ－sｃalewｉndturbineflａnge.Keｙwords:finiteelementmeｔｈod；ｔｏｗertoｐflaｎge;nonｌinearcontａct；windｔurbinｅ塔筒顶部法兰是风力发电机运行时的主要受力部件之一，它与偏航刹车盘、偏航轴承通过螺栓联接。在风力发电机组工作时，作用于叶轮和机舱的载荷将通过偏航轴承传递到塔架上。这些载荷包括叶轮所受到的气动载荷,如升力、阻力、轴向推力、气动弯矩和气动转矩等。由于风速风向的不稳定性，机组运行时塔筒顶部法兰承受的载荷是动态变化的,很容易遭到破坏.为提高风力发电机运行的可靠性,需要对塔筒顶部法兰联接部分进行结构强度分析.本文采用整体建模的思路,利用有限元理论分析塔筒顶部法兰面的接触应力,并且对塔筒顶部法兰与塔筒筒体的焊缝连接进行疲劳寿命分析.1整体模型的接触分析1。1几何模型几何模型在NX/UＧ中建立，本文重点分析塔筒顶部法兰的接触,所以简化轴承内圈,建立偏航轴承外圈与偏航刹车盘和塔筒顶部法兰的螺栓联接模型。如图１所示。1。塔筒(塔顶连接法兰)2。偏航刹车盘3.偏航轴承外圈4。螺栓1．2有限元模型在MSＣ．Mａrc/Ｍeｎtat中建立塔筒顶部法兰接触的有限元模型,如图2所示。模型均采用六面体单元划分，为了减小结构的关键部位对计算结果的影响，在塔筒顶部法兰圆角处的网格适度加密，单元总数为７61０4,节点总数为１1３48６。模型包括偏航轴承外圈、偏航刹车盘、塔筒顶部法兰和８4个螺栓等可变形接触体,各部件所用的材料均为低合金高强度结构钢，其弹性模量为2.06×10５MPa,泊松比为０。3，密度为7。８5×10—9t/ｍｍ３，根据实际的装配关系，建立如表1所示的接触关系.表1中：Ｔ为接触(toucｈing）；G为粘接（glue）;括号中数值为接触对的摩擦因数。螺栓与偏航轴承外圈在螺纹联接处有５5mｍ的螺纹啮合，本模型中通过粘接（gｌue）关系来耦合接触双方的自由度,以模拟螺纹联接。图2塔筒顶部法兰接触的有限元模型1.3载荷设置和边界条件接触分析所采用的载荷来自风力发电机极限工况载荷数据。分析分两种工况进行:第一种工况中给84个螺栓施加预紧力F0(ｐreloａd_bｏlt）2５4kＮ，建立各部件间的接触关系；第二种工况中撤掉预紧力，锁定螺栓的长度(ｌｏｃｋ_bolt)，然后施加偏航轴承极限载荷(lｏad_bｅarｉng)和偏航刹车盘的极限载荷（ｌoad_bｒａｋｅ)，分析各零／部件的应力分布情况.为了准确模拟载荷的传递，在偏航刹车盘中心点处施加极限载荷报告中此处载荷MZ的最大值为－328５kN·m,通过MSＣ.Ｍarc／Ｍentat软件施加载荷单元ＲBＥ3,将载荷分布到偏航刹车盘上七个与刹车器接触的区域(如图３所示）,模拟偏航刹车。在轴承外圈中心点与轴承外圈内表面建立ＲBE3连接，施加偏航轴承处的极限载荷为:ＦX=—４44.2kNFZ＝-１３５0．8kNMXY=-4570．７kN·m图3刹车盘的结构约束塔筒底部端面所有节点X、Y、Z方向的平动自由度(ｆix＿tｏwer），以消除接触分析时的刚体位移。同时考虑重力(ｇraviｔy)的影响。所有载荷设置和边界条件的施加如图４所示。图4载荷设置和边界条件（二分之一模型)1．4接触分析结果螺栓预紧后,各接触面接触良好，各零／部件应力值均在屈服强度之下.施加极限载荷后,各接触面间仍保持很好的接触状态,各零/部件的应力值稍有增大，下面给出只施加预紧力和极限载荷后的各零/部件接触应力云图.偏航轴承外圈的最大应力值为191．2ＭＰａ，出现在螺栓联接的第一圈螺纹处（见图5),小于轴承材料的屈服极限650MPa,轴承强度满足要求.图５偏航轴承外圈的应力云图（极限载荷工况下)由于接触面比较小，塔筒顶部法兰的最大应力值在预紧工况下为29３．２MPa(见图６）；施加极限载荷后为294.7MPa（见图７)，均出现在塔筒顶部法兰面与螺栓接触面周围，且都小于塔筒顶部法兰材料的屈服极限325MＰａ,塔筒顶部法兰强度满足要求。图６塔筒顶部法兰的应力云图（预紧工况下）图7塔筒顶部法兰的应力云图（极限载荷工况下）塔筒顶部法兰联接螺栓的最大应力在预紧工况下为496.４MPa(见图8),施加极限载荷后为49９.9MPa（见图9），均出现在螺栓的中段，且都小于螺栓的屈服极限900MPａ(该螺栓为１０.9级），所以螺栓是安全的.图8螺栓的应力云图（预紧工况下）图９螺栓的应力云图（极限载荷工况下）由于塔筒顶部法兰圆角处通常会出现应力集中，所以对圆角进行安全校核。塔筒顶部法兰圆角处最大等效应力值为102.８MPa，如图10所示。图１0塔筒顶部法兰圆角处的应力云图塔筒顶部法兰圆角处安全裕度Msult为：Ｍsulｔ＝[σ］/σm-1（１)[σ］=σs/n（2）式中:σm为极限载荷下塔筒顶部法兰圆角处的等效应力,σm=1０２。8MPａ;σs为塔筒顶部法兰材料屈服极限,σs=3２5MPa；ｎ为安全系数，n取1.１.则:Msult=1。８741塔筒顶部法兰圆角处的安全裕度大于零,因此塔筒顶部法兰圆角处静强度符合要求。2塔筒顶部法兰接触面的安全性校核由于受弯矩作用，除了要求静强度足够外,还应保证接合面间不出现缝隙和接触面不被压溃.由于预紧力作用,接触面上产生的压应力σF0为：σF０=ｚＦ０／A=8４×25４/０．８６7×10—3=24.6２ＭPａ式中：z为螺栓个数；A为接触面的面积；F０为螺栓预紧力。由于弯矩的作用，接触面上产生的弯应力σM为：σM=MXＹ／W=45７0．7／0。5997×10－3=７.62１MＰa式中：MXY为接触面处的弯矩；Ｗ为抗弯截面系数.压力FＺ在接触面上产生的压应力σＺ为:σZ=FZ/A＝-１350。８／０。８6７×1０-3=－1．559MPa接触面受拉侧不出现缝隙的条件为:σl=σF0-σM+σZ≥０(3)式中：σl为接触面所受拉应力。接触面受压侧不被压溃的条件为:σY＝σF0+σM+σZ≤［σp］(4）式中:σＹ为接触面所受压应力；［σp］为被连接材料的许用挤压应力，[σp]=0.8σｓ＝276MPa。由此可得:σｌ＝15。４4MPaσＹ=30.6８ＭPa〈［σｐ]σl大于零，塔筒顶部法兰接触面不会产生缝隙。σＹ小于被连接材料的许用挤压应力，所以塔筒顶部法兰不会被压溃.３塔筒顶部法兰焊缝疲劳寿命计算塔筒顶部法兰和塔筒筒体通过焊接连接，是结构中比较薄弱的环节，焊接结构的耐久性很大程度上由焊缝的疲劳强度所决定。所以,如果能在设计初期对焊缝的疲劳寿命进行预测,以了解整个焊接结构上焊缝疲劳强度的分布情况,就可以对实际中焊缝的处理方式进行合理地调整，以提高产品的整体疲劳性能。3.1焊缝等效疲劳应力计算法兰与塔筒筒体之间采用等厚度焊接,因此计算时不考虑焊缝应力集中系数的影响.焊缝等效疲劳应力采用有限元方法通过加载相应位置的等效疲劳载荷而获得。具体有限元计算过程如下。3.1.1模型描述使用接触分析的模型.第二种工况(载荷施加工况)的载荷由极限载荷改为等效疲劳载荷Fx＝-1４2.62kN，Fz=38kＮ,Mｙ=—1592kN·m,Mz=1４59ｋＮ·ｍ。3。1。2有限元计算结果计算得到法兰焊缝处最大等效应力为22.９５MPａ,如图１１所示，最大等效应力位置的节点坐标为（-1460.95，－54。42,-１90。３２４）.图1１塔筒顶部法兰焊缝处的应力云图3.2确定焊缝材料的疲劳寿命S—Ｎ曲线焊缝材料S-N曲线如图12所示。图12S－N曲线S—N曲线(以循环基数ND为界):Ｎ=NＤ（σD/σw)3(Ｎ<ND，σw>σD）(5）N=ND(σＤ/σw)4（N≥NＤ，σｗ≤σD）（6）式中：N为循环次数；σD为循环基数ND对应的疲劳应力；σｗ为计算所得应力.3.3焊缝疲劳应力循环次数计算有限元方法计算得到的焊逢处最大等效疲劳应力σｐh=2２。95MPａ。因此焊逢处等效疲劳应力σ′pｈ为：σ′pｈ=γmσｐｈ（7)式中：γm为材料偏安全系数,取1．１。计算得:σ′ph为２５.245ＭPa。已知塔筒筒体材料的屈服极限σ′s=4７0MＰａ，强度极限σｂ=345ＭPa,故由疲劳极限估算公式得焊缝疲劳极限σA为:σA=0.23（σｂ+σ′s)（8)可求出焊缝疲劳极限σA=１87MPa,对应的循环次数ＮA为１×1０7。已知循环基数ND=5E６,循环次数大于ＮD时用式（9)计算σw为：σw＝σＡ(NA／N)1/４(９)所以：σD＝σA（ＮＡ/NＤ)1/４=1．１9σA根据GＬ规范,存活率大于97。7％的S－N曲线需要进行修正,缩减为原来的2／3.因此修正循环基数σ′D为:σ′D=２／3σD=14８MＰa由式（6)可得到σ′ph对应的循环次数N′ph=５.９×１0９。可见塔筒顶部法兰焊缝处的循环次数能够满足风力发电机组循环107次的寿命要求。4结语（1）在接触分析中，施加极限载荷时,各零/部件的应力值比在预紧工况下稍有增大（如图6与图７，图8与图9所示)，说明极限载荷对接触应力的贡献不大，所以选用适当的预紧力可提高螺栓联接的可靠性以及被联接件的强度。(2）偏航轴承外圈的最大应力发生在螺栓联接的第一圈螺纹处，这与实际相符合;螺纹联接受载时，螺纹间的载荷约有1／３集中在第一圈上。所以建议对偏航轴承螺孔的第一圈螺纹周围进行强化,或进行结构优化,如增加沉孔的倒角，如图１3所示，从而减小应力集中.(A）方法1（Ｂ）方法２图１3结构优化示意图（3)塔筒顶部法兰联接有限元法的分析结果为塔筒顶部法兰的设计计算和校核提供了必要的依据.同时此分析思路可应用到其他结构的法兰联接中,具有实际指导意义。参考文献：[1]托尼，伯顿.风能技术［Ｍ］．武鑫译。北京：科学出版社，2007.［2］陈火红．新编Mａrc有限元实例教程[M］。北京：机械工业出版社，２0０７。［３]彼得·艾伯哈特，胡斌。现代接触动力学[M］.南京:东南大学出版社