嚴 勇,劉楚達
(1.長沙航空職業(yè)技術(shù)學(xué)院 航空機械制造工程系,湖南 長沙 410124;2.長沙飚能信息科技有限公司,湖南長沙 410205)
近幾年,我國風(fēng)電行業(yè)發(fā)展迅猛,風(fēng)電機組的質(zhì)量安全也越來越受到重視[1]。一方面,隨著海上風(fēng)電的發(fā)展,風(fēng)力發(fā)電機的單機容量越來越大,2020年由三峽集團和東方電氣集團聯(lián)合研發(fā)的亞太地區(qū)最大、全球第二大10 MW海上風(fēng)電機組已經(jīng)成功并網(wǎng)發(fā)電;另一方面,關(guān)于風(fēng)力機一些重要部件的強度分析存在偏差,缺乏對風(fēng)力機的深入分析與研究,這些因素是導(dǎo)致目前機組運行過程中經(jīng)常發(fā)生事故的主要原因[2]。葉片是風(fēng)力機的關(guān)鍵核心部件,葉片承受的動態(tài)載荷和靜態(tài)載荷均通過葉根螺栓傳遞,隨著風(fēng)機功率的提高,葉根螺栓受力也越來越大,螺栓的安全關(guān)系到整個風(fēng)電機組的安全可靠運行,葉根高強度螺栓的斷裂問題已經(jīng)成為風(fēng)力機安全的研究重點。
葉片根部連接形式主要分為T型螺栓和螺栓套預(yù)埋式兩種[3],大多采用42CrMoA材質(zhì)的10.9級高強螺栓。目前,國內(nèi)針對葉根螺栓的強度分析、疲勞壽命分析以及斷裂失效研究多采用基于ANSYS等軟件的有限元仿真手段,其計算模型過于理想化,且缺乏葉根螺栓的實測載荷數(shù)據(jù)支撐。本文致力于螺栓載荷的在線監(jiān)測技術(shù)研究,開發(fā)了基于聲彈性原理的螺栓軸向力超聲測量及在線監(jiān)測設(shè)備,并在葉片實驗臺上針對某型號葉片極限載荷試驗中的葉根螺栓組載荷進行了實時監(jiān)測,通過實時監(jiān)測數(shù)據(jù)分析了葉根螺栓組的載荷分布狀況。
風(fēng)力機運轉(zhuǎn)時,葉根是葉片上載荷最大的部位,承受著復(fù)雜的剪切、擠壓和彎扭載荷的組合作用,這些作用力最終將通過葉根螺栓傳遞給輪轂,因此葉根螺栓必須具有足夠的強度、剛度和局部穩(wěn)定性以及疲勞斷裂強度,才能保證整機的安全運行。統(tǒng)計表明,因為連接螺栓疲勞強度不夠而導(dǎo)致螺栓斷裂的安全事故屢見不鮮[4],[5]。
本文針對某2 MW級風(fēng)力機葉片的極限載荷進行試驗,分別采用應(yīng)變片法和超聲法實時監(jiān)測葉根螺栓軸向力,結(jié)合ANSYS有限元仿真分析,研究了葉片在最小擺振、最大擺振、最小揮舞和最大揮舞4種不同工況時的葉根螺栓載荷變化。本次試驗葉片葉根節(jié)圓直徑φ為2 300 mm,共計92根螺栓,具體試驗內(nèi)容如下。
①葉片用吊車懸掛在實驗臺法蘭上,采用扭矩扳手對螺栓進行預(yù)緊,用超聲法和應(yīng)變片法測試空載狀態(tài)下的螺栓預(yù)緊力。
②依照設(shè)計數(shù)據(jù),對葉片施加最小擺振、最大擺振、最小揮舞和最大揮舞4種工況的載荷,分別采用超聲測試法、應(yīng)變片測試法以及ANSYS有限元仿真分析法研究葉根螺栓載荷變化。
本次試驗遵循的技術(shù)規(guī)范主要有《Guideline for the Certification of Wind Turbines Edition 2010》和《IEC 61400 Wind Turbines--Parts 23:Full-Scale Structural Testing of Rotors Blades》。
2.1.1螺栓軸向力超聲測試原理
根據(jù)聲彈性原理,超聲波的速度會因材料中緊固力的存在而產(chǎn)生微小的變化,通過研究被測螺栓軸向應(yīng)力與超聲波傳播時間變化率的關(guān)系,利用超聲波來測量被測螺栓中的軸向力[6],測量原理圖如圖1所示。該方法利用超聲波傳播速度或傳播時間等參量來表征螺栓應(yīng)力的大小。實際應(yīng)用中一般不直接測量超聲波波速,而是轉(zhuǎn)為測量渡越時間,即超聲波沿螺栓軸向傳播時所需的往返時間。
圖1 螺栓軸向力超聲測量原理圖Fig.1 Schematic diagram of ultrasonic measurement of bolt axial force
螺栓的軸向應(yīng)力σ可表示為[7]
式中:E為螺栓材料的彈性模量;V0為零應(yīng)力下的超聲聲速;L1為螺栓長度;t0為測量零應(yīng)力聲時的溫度;t1為測量應(yīng)力為σ時的溫度;S0為零應(yīng)力聲時;S1為測量應(yīng)力為σ時的聲時;Ks為應(yīng)力系數(shù),是機械伸長量和聲程增加量的比值,在誤差許可范圍內(nèi)該比值為材料系數(shù),與應(yīng)力無關(guān);Kt為溫度系數(shù),每10℃溫度變化引起的聲程變化率。
風(fēng)電葉片的葉根螺栓主要采用42CrMoA材料制造,可通過試驗標定Ks,Kt和V0的值,實測時只需測量S1,S0,t1和t0以及螺栓的有效受力長度L0,即可由式(1)計算出螺栓的σ和軸向力F。
2.1.2超聲測試點的布置
除了前述應(yīng)變測試中,葉根端面呈米字分布的8個測點外,擺振工況在葉片合??p前后緣分別多布置了5個測點,揮舞工況則在合??p垂直方向上多布置了5個測點。在全部92根螺栓中,布置了18個超聲螺栓軸向力測試點,其中1,12,24,35,47,58,70和81監(jiān)測點疊加了應(yīng)變片監(jiān)測,具體螺栓編號如圖2所示。
圖2 葉根螺栓軸向力測試點分布圖Fig.2 Distribution of test points for axial force of blade root bolt
探頭布置在螺栓端部,連接導(dǎo)線從葉根法蘭端面開槽,向葉片內(nèi)部引出,監(jiān)測設(shè)備置于葉根孔蓋內(nèi)部。圖3為螺栓軸向力監(jiān)測超聲傳感器的安裝現(xiàn)場。
圖3 螺栓軸向力監(jiān)測超聲傳感器的安裝Fig.3 Installation of ultrasonic sensor for bolt axial force monitoring
應(yīng)變測試方法測試螺栓軸力,是利用螺栓受力變形捕捉螺栓軸向應(yīng)變,再用平均應(yīng)變乘以螺栓模量和截面積,獲得螺栓的平均軸力。
式中:A為螺栓光桿截面面積;r為螺栓光桿截面半徑;ε為螺栓光桿截面平均測試應(yīng)變。
葉根螺栓為42CrMoA材質(zhì)的10.9級M36螺栓,光桿部分直徑為28 mm,由式(2)得出1μm的應(yīng)變對應(yīng)軸力為0.13 kN,經(jīng)過精密校準后的8根測試螺栓,應(yīng)變最小測試刻度為1μm,因此應(yīng)變片軸力測試的最小精度為0.13 kN。
應(yīng)變測試螺栓在螺栓光桿同一截面上貼4個應(yīng)變片,取4個應(yīng)變的平均值來換算軸力,標定螺栓如圖4所示。
圖4 應(yīng)變片測試標定螺栓Fig.4 Strain gauge test calibration bolt
2.3.1有限元模型的建立
對試驗葉片葉根連接載荷分析有限元模型進行簡化,將葉片部分的工況載荷等效作用到葉根,可大大簡化參數(shù)化模型,因此葉根模型主要包含3部分:葉根(含葉根法蘭)、連接螺栓和連接法蘭。有限元仿真網(wǎng)格采用Hypermesh 11.0軟件劃分,運用ABAQUS 6.11進行計算。連接法蘭、螺栓及葉根法蘭的材料參數(shù)如表1所示。
表1 葉根材料屬性Table 1 Root material properties
螺栓采用實體建模,其他部件均采用三維網(wǎng)格,單元類型采用C3D8R,總網(wǎng)格數(shù)目約為220萬個,整體網(wǎng)格和各個部件網(wǎng)格如圖5所示。
圖5 葉根有限元模型的構(gòu)建Fig.5 Construction of finite element model of blade root
2.3.2載荷及邊界
試驗葉片葉根連接主要承受的載荷,一方面來自螺栓預(yù)緊力,另一方面來自葉片外載荷彎矩。計算模型中,連接法蘭與實驗臺連接處采用剛性約束。根據(jù)前述試驗中采集的螺栓預(yù)緊力,施加預(yù)緊力平均載荷為175 kN。葉片外載荷簡化為4種工況,即最大揮舞、最小揮舞、最大擺振和最小擺振,在實驗臺上通過鋼絲繩模擬加載葉片靜態(tài)載荷。根據(jù)施加的外部載荷,可計算出上述工況下作用在葉根處的最大彎矩(表2),有限元模型載荷如圖6所示。
表2 不同工況下葉根處最大彎矩載荷Table 2 Maximum bending moment load at blade root under different working conditions
圖6 葉根有限元模型邊界及載荷定義Fig.6 Boundary and load definition of blade root finite element model
分別采用應(yīng)變片技術(shù)和超聲技術(shù)對本次風(fēng)力機葉片的靜載荷加載試驗中預(yù)緊工況以及外載荷工況下的葉根螺栓軸向力進行實時監(jiān)測。由于應(yīng)變片敷設(shè)在螺桿中部,需要在葉根部分打孔引出連接導(dǎo)線,施工條件較差,最終8根校準螺桿中只有6根存活。在安裝的24個螺栓軸向力超聲傳感器中,由于葉根法蘭與實驗臺接合面的碾壓作用,導(dǎo)致多個超聲傳感器失效。
預(yù)緊工況下,葉片由吊車吊起,92根葉根連接螺栓僅承受預(yù)緊力,分3次擰緊,實測預(yù)緊力為150~200 kN,平均值為175 kN。從應(yīng)變測試結(jié)果來看,每個測點同一根螺栓上的4個測點應(yīng)力均大小不一,相鄰的應(yīng)變片差達到最大應(yīng)變的一半,但對角應(yīng)變之和均相等,說明在預(yù)緊過程中,螺栓發(fā)生了較大的彎曲。在預(yù)緊過程中,螺栓彎曲的程度與試驗法蘭平整度、螺栓對中情況、螺紋潤滑、各個接觸表面的潤滑程度以及葉根端面平整度相關(guān)。從超聲監(jiān)測數(shù)據(jù)來看,與應(yīng)變片監(jiān)測數(shù)據(jù)吻合較好,相差在5%以內(nèi),完全滿足一般工程應(yīng)用需求,監(jiān)測數(shù)據(jù)如表3所示。
表3 預(yù)緊工況下螺栓軸向力監(jiān)測數(shù)據(jù)Table 3 Monitoring data of bolt axial force under preloading condition
分別用螺栓軸力應(yīng)變測試法、超聲測試法以及有限元計算對比外載荷作用下螺栓軸力的變化值,并以應(yīng)變法測試值作為比較基準計算了超聲測試法及有限元法的相對誤差(表4,5)。
由表4,5可知,超聲測試值與應(yīng)變測試值吻合相對較好,誤差基本在10%以內(nèi),有限元仿真分析的結(jié)果誤差較大。有限元計算誤差原因有以下幾點:①從預(yù)緊工況可知,此次測試螺栓預(yù)緊力存在一定的差異,有限元計算采用統(tǒng)一的平均值175 kN,預(yù)緊力誤差會引起變化值的誤差;②根據(jù)本次測試情況,試驗法蘭、螺紋以及各個接觸面潤滑情況與有限元計算有差異;③雖然有限元計算采用實體單元,考慮了外載荷對螺栓彎曲作用,但是實際安裝誤差、螺栓對中及端面垂直度、螺栓預(yù)緊等引起的螺栓彎曲無法考慮,也會對外載荷引起螺栓軸力變化帶來誤差。
表4 擺振工況下螺栓軸向力變化值應(yīng)變片測試、超聲測試及有限元仿真對比Table 4 Strain gauge test,ultrasonic test and finite element simulation of axial force variation of bolt under shimmy condition
表5 揮舞工況下螺栓軸向力變化值應(yīng)變片測試、超聲測試及有限元仿真對比Table 5 Comparison of strain gauge test,ultrasonic test and finite element simulation of bolt axial force variation under swing condition
從本次試驗數(shù)據(jù)來看,只需在關(guān)鍵位置合理布置傳感器,螺栓軸向力超聲監(jiān)測系統(tǒng)就能實時捕捉葉根彎曲的中性面,也能通過計算實時擬合出葉片機械載荷,為風(fēng)力發(fā)電機機械載荷監(jiān)測、運行安全保障及疲勞壽命預(yù)測提供支撐。
在風(fēng)力機運行過程中,其塔筒、輪轂、葉根等各個法蘭面承受的均為多變的動載荷,在這些載荷作用下,機械承載系統(tǒng)可能會發(fā)生過載或者疲勞,造成潛在的失效隱患。運用基于聲彈性原理的螺栓軸向力超聲在線監(jiān)測技術(shù),對風(fēng)力機各法蘭面螺栓進行實時監(jiān)測,通過監(jiān)測數(shù)據(jù)可實時分析螺栓受到的工作載荷及應(yīng)力幅值的變化,不但能夠直接判斷出當(dāng)前螺栓的工作狀態(tài),還能間接計算出螺栓所在法蘭面的工作狀態(tài),從而為設(shè)備的安全運行提供遠程數(shù)據(jù)支撐。