劉榮梅,孫 梁,周克印
(1.南京航空航天大學(xué) 航空宇航學(xué)院, 南京 210016;2.南京機(jī)電液壓工程研究中心, 南京 210061)
?
敲擊錘參數(shù)對風(fēng)力機(jī)葉片聲振法檢測深度的影響
劉榮梅1,孫梁2,周克印1
(1.南京航空航天大學(xué) 航空宇航學(xué)院, 南京 210016;2.南京機(jī)電液壓工程研究中心, 南京 210061)
摘要:為提高風(fēng)力機(jī)葉片聲振法的檢測深度,利用有限元軟件,對風(fēng)力機(jī)葉片聲振法檢測過程進(jìn)行動力學(xué)模擬,選取沖擊過程中應(yīng)力波能量輻射系數(shù)為參數(shù),分析敲擊錘參數(shù)對應(yīng)力波能量輻射系數(shù)的影響,以期得到更高的應(yīng)力波能量輻射系數(shù),從而達(dá)到提高聲振法檢測缺陷深度的目的。有限元模擬結(jié)果表明,除了敲擊錘材料的泊松比外,其直徑、錘身長度、錘頭材料彈性模量以及沖擊速度均對應(yīng)力波能量輻射系數(shù)產(chǎn)生影響。在模擬分析結(jié)果的基礎(chǔ)上,通過比較不同敲擊錘尺寸下的應(yīng)力持續(xù)時間,得出適當(dāng)選取敲擊錘參數(shù)可提高結(jié)構(gòu)聲振法缺陷的檢測深度的結(jié)論。
關(guān)鍵詞:風(fēng)力葉片;聲振法;無損檢測;敲擊錘參數(shù)中圖分類號: TM315;TB302.5;TG115.28
文獻(xiàn)標(biāo)志碼:A
文章編號:1000-6656(2016)05-0005-05
葉片是風(fēng)力機(jī)的關(guān)鍵部件之一,其大型化使得人們對其結(jié)構(gòu)完整性和安全性的要求越來越高[1]。由于葉片長期工作在野外環(huán)境條件中,所承受的風(fēng)沙載荷、雷擊損傷和環(huán)境溫差等因素極易引起其內(nèi)部損傷的發(fā)生和擴(kuò)展導(dǎo)致的失效[2]。從國內(nèi)外風(fēng)機(jī)的實(shí)際應(yīng)用來看,葉片斷裂和損毀的情況屢見不鮮。為確保風(fēng)力機(jī)的使用安全和質(zhì)量,對其脫層或者內(nèi)部裂紋的檢測變得尤為重要[3]。
常見的葉片無損檢測技術(shù)主要有超聲波檢測[4]、射線檢測[5-6]和紅外熱成像檢測[7]等。射線和紅外熱成像檢測成本高[8],而超聲檢測雖然簡單,但由于需要耦合劑,且如果試件表面粗糙或有晶粒狀物質(zhì)存在易產(chǎn)生雜亂反射波,在現(xiàn)場損傷檢測中并不實(shí)用。聲振法相對簡單,且可在風(fēng)力機(jī)葉片較低速運(yùn)行時對其進(jìn)行實(shí)測,因此被廣泛采用。其原理是通過撞擊構(gòu)件表面,分析撞擊產(chǎn)生的機(jī)械振動(或稱為聲波)數(shù)據(jù),來判斷物體的固有頻率和阻尼是否發(fā)生了變化,進(jìn)而判斷物體是否存在損傷[9]。
早在19世紀(jì)80年代,英國的研究工作者選擇敲擊所產(chǎn)生的力的持續(xù)信號作為判斷特征量[10],建立了敲擊聲振法的理論模型。采用風(fēng)力葉片聲振法檢測時,在沖擊錘頭受力處內(nèi)置壓電薄膜傳感器,通過分析傳感器采集的碰撞持續(xù)信號,從而實(shí)現(xiàn)對葉片的損傷檢測[11]。
為得到較好的檢測結(jié)果,筆者采用有限元軟件,模擬分析了敲擊錘參數(shù)對葉片檢測深度的影響,通過選用葉片沖擊振動應(yīng)力波能量輻射系數(shù)為參數(shù),研究了沖擊速度及敲擊錘尺寸對該參數(shù)的影響;并對敲擊錘尺寸及沖擊速度等參數(shù)進(jìn)行了設(shè)計,提高了傳感器葉片損傷的探測深度。
1基本原理
1.1聲振法檢測原理
聲振法檢測原理如圖1所示,敲擊錘以一定速度沖擊被測物,假設(shè)被測構(gòu)件為瞬態(tài)激勵,忽略材料阻尼等影響因素,構(gòu)件的運(yùn)動可簡化為無阻尼簡諧振動。
圖1 聲振法檢測示意
當(dāng)被測對象受到敲擊錘的瞬態(tài)激勵時,可由式(1)得到試件振動周期:
(1)
式中:m為被測對象的質(zhì)量;k為構(gòu)件的剛度,k由構(gòu)件形狀、材料及約束條件確定。
敲擊錘敲擊葉片后反彈,其振動的第一個正弦波前半個周期即為傳感器檢測到的敲擊信號持續(xù)時間t,可由下式表示:
(2)
沖擊作用下的被檢測對象可簡化為彈簧模型,當(dāng)其存在缺陷時,相當(dāng)于在完好的結(jié)構(gòu)中串聯(lián)了一個因缺陷引起的彈簧,引起結(jié)構(gòu)的剛度改變,由式(1)和(2)可知,被測對象振動周期(頻率的倒數(shù))和敲擊持續(xù)時間t將發(fā)生變化。
敲擊錘端部的壓電傳感器采集被測物的沖擊運(yùn)動加速度,壓電元件表面產(chǎn)生的電荷q與沖擊加速度a之間滿足下式[12]:
(3)
式中:d11,KY及wn分別為壓電常數(shù)、壓電元件剛度系數(shù)及傳感器固有頻率。
由式(3)可知,傳感器采集到的壓電元件表面的電壓信號隨沖擊加速度的變化而變化。
1.2參數(shù)的選擇
為研究敲擊錘參數(shù)對應(yīng)力波輻射能量的影響,采用能量輻射系數(shù)λ來表征[13],其定義如下:
(4)
應(yīng)力波能量輻射系數(shù)λ的意義在于有多少初始的碰撞能量轉(zhuǎn)化為可用于檢測的應(yīng)力波能量,這些能量不但對試件產(chǎn)生激勵,還通過應(yīng)力波的傳遞不斷與內(nèi)部缺陷發(fā)生相互作用,輻射的能量越高,傳感器探測的缺陷距離也將越遠(yuǎn)。
通過模擬小錘撞擊無損傷葉片,研究沖擊速度、敲擊錘直徑及錘身長度等參數(shù)對聲振法檢測過程的影響,以獲得較高的能量輻射系數(shù)。
1.3有限元模型
為研究不同敲擊錘參數(shù)對傳感器檢測效果的影響,采用有限元軟件,模擬敲擊錘敲擊風(fēng)力葉片表面的檢測過程。
由于風(fēng)機(jī)葉片尺寸巨大[14],且敲擊過程為局部振動,為簡化問題,截取長寬高分別為1 m,75 cm及3.5 cm的板為研究對象,其邊界條件設(shè)為四邊固支。葉片材料為玻璃纖維增強(qiáng)復(fù)合材料層合板,其彈性模量、泊松比及密度分別為17 GPa、0.3和4 000 kg·m-3。為避免邊界效應(yīng)對模擬結(jié)果的影響,沖擊位置設(shè)在板中心。
模擬的葉片沖擊模型如圖2所示,在沖擊中心5.0 cm的方形范圍內(nèi)劃分為較細(xì)的網(wǎng)格,其余網(wǎng)格逐漸變稀。
在建立模型的基礎(chǔ)上,可以分析沖擊過程中板的振動情況。結(jié)合傳感器工作原理,選取錘頭位置在板上所對應(yīng)的點(diǎn)為對象,研究該點(diǎn)在不同工作狀況下的速度,圖3所示為有限元模擬分析得到的某工況下錘頭沖擊速度。根據(jù)沖擊前后的模擬速度,由式(4)可以得到不同工況下的應(yīng)力波能量輻射系數(shù)。
圖2 葉片沖擊模型
圖3 錘頭沖擊速度模擬結(jié)果
2敲擊錘參數(shù)對輻射應(yīng)力波能量的影響
不考慮沖擊參數(shù)對葉片結(jié)構(gòu)破壞的前提下,討論敲擊錘各設(shè)計參數(shù)對檢測結(jié)果的影響。
2.1敲擊錘尺寸的影響假定在錘頭材料為尼龍,初始沖擊速度為1 m·s-1的情形下,分別改變敲擊錘直徑及錘身長度,討論敲擊錘尺寸對沖擊應(yīng)力波輻射能量的影響。
2.1.1敲擊錘直徑的影響
在錘身長度為8 cm的條件下,改變敲擊錘直徑,研究無缺陷試件板在進(jìn)行垂直敲擊下的能量輻射系數(shù)λ,λ隨敲擊錘直徑變化曲線如圖4所示。
圖4 應(yīng)力波輻射能量系數(shù)隨敲擊錘直徑變化曲線
由圖4 可以看出,應(yīng)力波能量輻射系數(shù)λ先是隨敲擊錘直徑急劇變化,當(dāng)直徑大于3 cm時,λ隨直徑的增大呈緩慢增加趨勢;當(dāng)直徑大于10 cm時,λ幾乎不隨直徑變化。
2.1.2敲擊錘錘身長度的影響
假設(shè)錘頭直徑為5 cm,無缺陷試件板在垂直敲擊條件下,能量輻射系數(shù)隨錘身長度變化模擬結(jié)果如圖5所示。
圖5 應(yīng)力波輻射能量系數(shù)隨錘身長度變化曲線
由圖5可以看出,應(yīng)力波能量輻射系數(shù)λ隨錘身長度而變化,λ先隨錘身長度的增加而增大,當(dāng)錘身長度為7 cm左右時,λ數(shù)值最大,之后,λ隨著錘身長度的增加緩慢減小。
2.2碰撞速度的影響
在敲擊錘直徑為5 cm、錘身長度為8 cm、沖擊角度為垂直沖擊、錘頭材料為尼龍的條件下,以不同初速度對無缺陷試件板進(jìn)行垂直沖擊模擬,能量輻射系數(shù)λ隨碰撞速度的變化如圖6所示。
圖6 不同碰撞速度下的能量輻射系數(shù)
由圖6可知,隨著碰撞速度的增加,λ增大,并可分為三個區(qū)域。當(dāng)碰撞速度小于0.8 m·s-1時,能量輻射系數(shù)基本維持在一個恒定值 0.32;當(dāng)碰撞速度大于3 m·s-1時,能量輻射系數(shù)上升速率緩慢。
當(dāng)碰撞速度在 0.8~3 m·s-1范圍內(nèi)時,λ與初速度近似呈線性關(guān)系,這種關(guān)系與HUNTER的研究結(jié)果相類似[15],其關(guān)系如下:
(5)
式中:C0為縱波在板中的波速。
2.3錘頭材料參數(shù)的影響
設(shè)定沖擊速度和錘頭直徑分別為1 m·s-1和5 cm,錘身長度為8 cm,改變敲擊錘錘頭的材料參數(shù),對垂直沖擊下,板的應(yīng)力波能量輻射系數(shù)進(jìn)行分析。
2.3.1錘頭材料彈性模量的影響
首先假定錘頭材料泊松比為0.3,密度為7 800 kg·m-3,能量輻射系數(shù)隨著彈性模量的變化曲線如圖7所示。
圖7 能量輻射系數(shù)隨彈性模量的變化曲線
由圖7可見,錘頭的泊松比、密度不變,而彈性模量改變時,能量輻射系數(shù)隨錘頭的彈性模量E增加而增加,在E接近200 GPa時,λ逐漸趨于穩(wěn)定。
2.3.2錘頭材料泊松比的影響
假定錘頭材料的彈性模量為200 GPa,密度為7 800 kg·m-3,泊松比對能量輻射系數(shù)的影響見圖8。
圖8 泊松比對能量輻射系數(shù)的影響
由圖8可見,應(yīng)力波能量輻射系數(shù)隨錘頭材料泊松比增加而增大,但λ基本在0.7~0.71范圍內(nèi),可見泊松比對λ影響不大。
3參數(shù)優(yōu)化算例
根據(jù)以上研究結(jié)果來設(shè)計敲擊錘,研究了兩個算例:① 選用鋼材(假設(shè)泊松比為0.3,彈性模量為200 GPa,密度為7 800 kg·m-3)作為錘頭材料,錘頭直徑為5 cm,錘身長度為8 cm,從10 cm高度垂直下落,即落錘沖擊速度為1.4 m·s-1。② 錘頭材料選用同樣的鋼材,而錘頭直徑和錘身長度分別取10 cm和7 cm,敲擊錘以1 m·s-1沖擊速度垂直敲擊含缺陷的層合板。
為對敲擊錘設(shè)計結(jié)果檢測,利用有限元軟件分析敲擊錘垂直敲擊層合板,缺陷設(shè)計為直徑8 cm、厚4 mm、水平位置在板中心的圓柱體,檢測模擬結(jié)果如圖9所示。
圖9(a)為算例①下,敲擊錘與板接觸時間隨缺陷垂直位置變化。由圖9(a)可見,隨著缺陷在板內(nèi)層間位置的變化,敲擊錘與板接觸時間(即傳感元件檢測到的信號持續(xù)時間)發(fā)生變化,然而當(dāng)缺陷距表面距離超過10 mm之后,沖擊信號持續(xù)時間無明顯變化,所以采用該方案可檢測缺陷的最大深度為10 mm。
圖9 不同算例下時間-深度關(guān)系曲線
根據(jù)有限元計算分析結(jié)果,進(jìn)行錘頭優(yōu)化設(shè)計,選取參數(shù),即算例②,模擬結(jié)果見圖9(b)??梢?,相比算例①,最大檢測深度為15 mm。因此通過改變敲擊錘參數(shù),檢測效果可以得到較大幅度的提高。
4結(jié)語
以應(yīng)力波能量輻射系數(shù)λ為評價參數(shù),研究了聲振法檢測葉片缺陷中的敲擊錘參數(shù)的設(shè)計。有限元模擬結(jié)果表明,除了錘頭材料泊松比外,其他參數(shù)對λ也產(chǎn)生明顯影響。
在給定條件下,當(dāng)敲擊錘直徑小于10 cm時,λ隨直徑的增大而增加;當(dāng)錘身長度為7 cm左右時,λ最大;錘頭材料彈性模量E不大于200 GPa時,λ隨E的增大而增加;垂直沖擊速度v滿足1~3 m·s-1范圍時,λ隨v線性增加。
由于未曾驗算各沖擊條件下接觸點(diǎn)的沖擊應(yīng)力,應(yīng)用該方法設(shè)計敲擊錘時,應(yīng)考慮實(shí)際葉片結(jié)構(gòu)的強(qiáng)度??傊鶕?jù)參數(shù)對應(yīng)力波能量輻射系數(shù)的影響研究結(jié)果,可設(shè)計敲擊錘尺寸及沖擊初速度。有限元模擬結(jié)果表明,改進(jìn)參數(shù)后,采用聲振法檢測風(fēng)力葉片缺陷,其檢測深度可得到提高。
參考文獻(xiàn):
[1]FAUSTO P G M, ANDREW M T, JESUS M P P, et al. Condition monitoring of wind turbines: Techniques and methods[J].Renewable Energy, 2012,46:169-178.
[2]羅慧敏,張錦南,陳余岳,等.我國大型風(fēng)力發(fā)電葉片產(chǎn)業(yè)的現(xiàn)狀分析[J].玻璃鋼,2008(4):1-5.
[3]LIU P, GROVES R M, BENEDICTUS R. 3D monitoring of delamination growth in a wind turbine blade composite using optical coherence tomography[J].NDT&E International, 2014, 64: 52-58.
[4]JASIUNIENE E, RAISUTIS R, SLITERIS R, et al. Ultrasonic NDT of wind turbine blades using contact pulse-echo immersion testing with moving water container[J].Ultrasound, 2008, 63(1): 28-32.
[5]GHOSHAL A, SUNDARESAN M J, SCHULZ M J, et al. Structural health monitoring techniques for wind turbine blades[J].Journal of Wind Engineering and Industrial Aerodynamics, 2000, 85, 309-324.
[6]SCHILLING P J, KAREDLA B R, TATIPARTHI, et al. X-ray computed microtomography of internal damage in fiber reinforced polymer matrix composites[J]. Composite Science & Technology, 2005,65:2071-2078.
[7]岳大皓,李曉麗,張浩軍,等.風(fēng)電葉片紅外熱波無損檢測的實(shí)驗探究[J].紅外技術(shù),2011,33(10):614-617.
[8]YANG Wen-xian, TAVNER P J, CRABTREE C J, et al. Wind turbine condition monitoring: technical and commercial challenges[J].Wind Energy,2014,17(5): 673-693.
[9]冷勁松,杜善義, 顧震隆. 配橡膠內(nèi)襯復(fù)合材料板殼的敲擊法無損檢測[J].材料科學(xué)與工藝,1993(2):51-55.
[10]CAWLEY P, ADAMS R D.The mechanics of the coin-tap method of non-destructive testing[J]. Journal of Sound and Vibration, 1988, 122(2): 299-316.
[11]陶鵬,趙一中,姚恩濤,等.基于敲擊聲振法的風(fēng)機(jī)葉片脫層檢測系統(tǒng)設(shè)計[J].測控技術(shù),2014, 33(4):12-15.
[12]陳至坤,逄鵬,王福斌,等.基于壓電加速度傳感器的振動測量系統(tǒng)研究[J].計算機(jī)測量與控制,2014,22(12):3923-3925.
[13]雷剛.應(yīng)力波在橋梁預(yù)應(yīng)力結(jié)構(gòu)中的傳播特性及應(yīng)用[D].湘潭:湘潭大學(xué),2013.
[14]汪泉,陳進(jìn),郭小峰,等.基于性能及輕量化的新型風(fēng)力機(jī)葉片優(yōu)化設(shè)計研究[J].空氣動力學(xué)學(xué)報,2014,32(1):77-84.
[15]HUNTER S C. Energy absorbed by elastic waves during impact [J]. Journal of the Mechanics and Physics of Solids, 1957, 5(3): 162-171.
Influence of Impact Hammer on Parameter Acoustic Vibration Detecting Depth inside Wind Turbine Blade
LIU Rong-mei1, SUN Liang2, ZHOU Ke-yin1
(1.College of Aerospace Engineering, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China;2.Nanjing Electrical and Hydraulic Engineering Research Center, Nanjing 210061, China)
Abstract:Acoustic vibration has been applied in nondestructive examination of wind turbine blade because of its simplicity and applicability. In order to detect the deep defect, the dynamic procedure during the acoustic vibration of the wind turbine blade was simulated by finite element method (FEM). The influence of the hammer size and impact velocity on the stress wave energy radiation coefficient (SWERC) during the vibration was studied. The larger the radiation coefficient was, the deeper delamination could be recognized experimentally. The FEM analysis indicated that the parameters, such as the diameter and the length of the hammer, elastic modulus of the hammerhead and original impact velocity, had obvious influence on the SWERC. Based on the simulated results, two examples considering different hammering parameters were studied. The study indicates that an appropriate hammer could be designed to expand the detection depth.
Key words:Wind blade; Acoustic vibration; Nondestructive testing; Impact hammer parameter
收稿日期:2015-09-03
基金項目:國家科技支撐項目計劃資助項目(2012BAA01B02);江蘇省博士后基金資助項目(1302071B)
作者簡介:劉榮梅(1975-),女,博士,副教授,主要從事測試技術(shù)研究工作。通信作者: 劉榮梅, E-mail: romme@nuaa.edu.cn。
DOI:10.11973/wsjc201605002