李 躍,肖小兵,蔡永翔,張紫嫣,魯彩江

(1. 貴州電網(wǎng)有限責(zé)任公司電力科學(xué)研究院， 貴陽(yáng) 550005；2. 西南交通大學(xué) 機(jī)械工程學(xué)院， 成都 610031)

0 引 言

雷電測(cè)量在雷電防護(hù)技術(shù)中至關(guān)重要，是雷電研究的基礎(chǔ)。因此，在新型電力系統(tǒng)建設(shè)背景下，要防止和減少電網(wǎng)事故，就要準(zhǔn)確獲知雷電流大小以及波形，在此基礎(chǔ)上結(jié)合相關(guān)風(fēng)險(xiǎn)評(píng)估系統(tǒng)，構(gòu)建電網(wǎng)全方位預(yù)警評(píng)估系統(tǒng)，其中蘊(yùn)含著顯著的學(xué)術(shù)價(jià)值、巨大的社會(huì)和經(jīng)濟(jì)效益，實(shí)現(xiàn)這一切不可或缺的工具就是瞬態(tài)磁場(chǎng)傳感器。

在瞬態(tài)電流測(cè)量領(lǐng)域，最常用的是羅氏線圈[1-3]和光學(xué)傳感器[4]。Rogowski線圈克服了磁性材料非線性飽和引起的信號(hào)失真問(wèn)題。然而，當(dāng)檢測(cè)磁場(chǎng)時(shí)，它需要一個(gè)積分器。光纖電流傳感器比電磁感應(yīng)變壓器和羅氏線圈具有更寬的頻帶響應(yīng)范圍和更好的電氣絕緣特性。然而，它容易受到溫度和振動(dòng)的影響，已成為制約其推廣的瓶頸。它需要激光器和調(diào)制器，因此系統(tǒng)復(fù)雜、成本高、體積大。

近年來(lái)，由于磁致伸縮/壓電復(fù)合材料具有很強(qiáng)的磁電耦合且易于制備，受到了廣泛研究[5-6]。磁電復(fù)合材料是制作高靈敏度、高分辨率磁傳感器件的理想功能材料，其靈敏度高于霍爾元件與超磁阻傳感器，頻率范圍比磁通門(mén)大；雖然靈敏度比超導(dǎo)量子干涉儀小，但是其結(jié)構(gòu)簡(jiǎn)單、制備容易、成本低、可在室溫下使用[7]。由于磁致伸縮材料的偏置磁場(chǎng)依賴特性，其可作為直流磁場(chǎng)傳感器[8-10]；在外部交流磁場(chǎng)下，磁致伸縮相施加機(jī)械應(yīng)力，然后轉(zhuǎn)移到壓電相，產(chǎn)生交流電荷輸出，其可作為交流磁場(chǎng)傳感器[11-12]、位移傳感器[13-14]、電流傳感器[15-21]、轉(zhuǎn)速傳感器[22]。然而，迄今為止，對(duì)磁電復(fù)合傳感器進(jìn)行的大多數(shù)理論和實(shí)驗(yàn)研究都是在靜態(tài)磁場(chǎng)和標(biāo)準(zhǔn)正弦波磁場(chǎng)激勵(lì)的前提下進(jìn)行的[1-22]。因此，為了擴(kuò)展磁電復(fù)合材料在電力系統(tǒng)雷擊、故障等瞬態(tài)電流測(cè)量領(lǐng)域中的潛在應(yīng)用，有必要對(duì)磁電復(fù)合材料在瞬態(tài)磁場(chǎng)/激勵(lì)激勵(lì)下的磁電響應(yīng)開(kāi)展研究。

本文利用有限元仿真元件，分析FeGa/PZT圓柱體磁電復(fù)合材料在8/20 μs雷電激勵(lì)下磁電響應(yīng)，為磁電復(fù)合材料在雷電測(cè)量領(lǐng)域提供理論分析基礎(chǔ)。

1 有限元仿真模型建立

為充分研究FeGa/PZT磁電復(fù)合材料在8/20 μs雷電激勵(lì)下磁電響應(yīng)，本文將利用有限元分析軟件Comsol對(duì)其進(jìn)行模擬分析。本文分析的FeGa/PZT圓柱體磁電復(fù)合材料結(jié)構(gòu)如圖1(a)所示。利用Comsol軟件直接對(duì)FeGa/PZT復(fù)合結(jié)構(gòu)進(jìn)行建模，采用二位軸對(duì)稱進(jìn)行模型建立，結(jié)構(gòu)如圖2(b)所示。其參數(shù)為12 mm×16 mm(Lp×Lm)，其中tm和tp分別為2和0.5 mm，結(jié)構(gòu)兩端均為自由狀態(tài)。圖中隱藏了所建立的空氣域與線圈幾何模型，空氣域?yàn)榘霃綖?0 mm、高度為180 mm的圓柱體。線圈為內(nèi)徑10 mm，外徑10.5 mm，高度為50 mm的圓管形，磁電復(fù)合結(jié)構(gòu)置于其中心處，線圈中的施加的激勵(lì)電流采用外部電流密度。磁致伸縮材料根據(jù)FeGa材料參數(shù)自行設(shè)置，壓電材料選取庫(kù)中的PZT-5H。

圖1 (a)磁電復(fù)合結(jié)構(gòu)幾何模型 (b)網(wǎng)格劃分后幾何模型Fig.1 Geometric model of the magnetoelectric composite before and after lattice division

仿真中，對(duì)磁致伸縮材料、壓電材料以及線圈采用自由四邊形網(wǎng)格劃分，設(shè)置最大單元大小為0.75 mm，其余剩余空氣域采用默認(rèn)的自由三角形網(wǎng)格劃分。網(wǎng)格劃分完成后如圖1(c)所示，圖中隱藏了線圈以及空氣域。

2 正弦交流激勵(lì)下磁電響應(yīng)仿真分析

首先，利用COMSOL中頻域求解器對(duì)所建立磁電復(fù)合材料的諧振頻率進(jìn)行求解。設(shè)置線圈中外部電流密度為1 000 A/m2，交流電流頻率從1 kHz變化到400 kHz，在求解器中設(shè)置所選頻段內(nèi)選取400個(gè)點(diǎn)。求解結(jié)束后，為獲得磁電復(fù)合材料的輸出電壓，對(duì)壓電管內(nèi)壁建立二維數(shù)據(jù)集，其性質(zhì)選擇加和。

在COMSOL軟件的結(jié)果選項(xiàng)中選擇一維繪圖組，數(shù)據(jù)源選擇已經(jīng)建立的二維數(shù)據(jù)集，選擇點(diǎn)圖繪制之后得到壓電管正極電勢(shì)隨交流激勵(lì)頻率變化如圖得2所示。由于壓電管外壁為接地，電勢(shì)為0，因此壓電材料的輸出電壓即為其正極電勢(shì)。從圖2中可以看到，磁電復(fù)合材料在一階諧振頻率100.36 kHz時(shí)輸出電壓達(dá)到最大，約為1.15 V，在二階、三階諧振頻率處輸出電壓約為0.1 V，遠(yuǎn)遠(yuǎn)小于一階諧振頻率處的輸出電壓。

隨后，對(duì)磁電復(fù)合結(jié)構(gòu)諧振振動(dòng)模態(tài)進(jìn)行了分析。復(fù)合材料在一階諧振、二階諧振和三階諧振處的振型如圖3所示。從圖中可知，在一階諧振頻率100.36 kHz時(shí)，整個(gè)結(jié)構(gòu)明顯處于沿軸向的拉伸狀態(tài)，而在二階、三階諧振頻率下，有很明顯的彎曲和扭轉(zhuǎn)現(xiàn)象。綜上可知，磁電復(fù)合材料的主振模態(tài)為縱波振動(dòng)，頻率為100.36 kHz，或固有周期為9.964 μs。

圖2 磁電復(fù)合材料輸出電壓隨激勵(lì)頻率變化曲線Fig.2 Output voltage of the magnetoelectric composite as a function of the AC current frequency

圖3 磁電復(fù)合結(jié)構(gòu)(a)一階、(b)二階和(c)三階振動(dòng)模態(tài)Fig.3 First, second and third order vibration modes of the magnetoelectric composite

3 8/20 μs雷電流雙指數(shù)函數(shù)模型

描述8/20 μs雷電流的函數(shù)有多種，常用的函數(shù)模型有雙指數(shù)函數(shù)模型、Heidler函數(shù)模型和脈沖函數(shù)模型。本文仿真中，采用雙指數(shù)函數(shù)模型作為雷電流瞬態(tài)激勵(lì)信號(hào)，如公式(1)。

(1)

根據(jù)公式(1)，繪制雷電流時(shí)域波形如圖4所示，由于在仿真分析中只能使用函數(shù)作為電流激勵(lì)源，因此其波形與國(guó)家標(biāo)準(zhǔn)8/20 μs雷電流的波形不太相同。從圖中可知，電流從0 μs開(kāi)始迅速上升，在約12 μs時(shí)達(dá)到峰值，隨后逐漸下降，在約100 μs時(shí)已經(jīng)完全降為零。

圖4 8/20 μs雷電流雙指數(shù)函數(shù)模型Fig.4 Double exponential function model of 8/20 μs lightning current

3.1 8/20 μs雷電流激勵(lì)下，F(xiàn)eGa/PZT磁電復(fù)合材料時(shí)域波形

在8/20 μs雷電流激勵(lì)下，將瑞利阻尼中質(zhì)量阻尼參數(shù)設(shè)置為200，剛度阻尼參數(shù)設(shè)置為8e-8。計(jì)算求解之后，得到FeGa/PZT磁電復(fù)合材料輸出電壓曲線如圖5所示。從圖中可以看到，磁電復(fù)合材料的輸出電壓從零點(diǎn)開(kāi)始迅速增加，由于振蕩的影響，在9 μs時(shí)達(dá)到峰值，之后逐漸下降到零點(diǎn)，并在下降過(guò)程中伴隨很明顯的振蕩現(xiàn)象，振蕩周期為10 μs。輸出電壓到達(dá)峰值的時(shí)間比激勵(lì)電流源提前3 μs，這是由于在上升過(guò)程中也伴隨振蕩現(xiàn)象，在9 μs剛好處于振蕩中，因此造成了峰值提前。

圖5 8/20 μs雷電流激勵(lì)下，磁電復(fù)合材料的輸出電壓Fig.5 The output voltage of magnetoelectric composite driven by 8/20 μs lightning current

由此可知，在8/20 μs雷電流激勵(lì)下，F(xiàn)eGa/PZT磁電復(fù)合材料存在阻尼振蕩現(xiàn)象，振蕩頻率等于材料一階縱波振動(dòng)頻率(100.36 kHz)。這可以解釋為，F(xiàn)eGa/PZT磁電復(fù)合材料是機(jī)械慣性器件，在脈沖激勵(lì)下會(huì)發(fā)生阻尼振蕩。

3.2 8/20 μs雷電流激勵(lì)下，F(xiàn)eGa/PZT磁電復(fù)合材料不同時(shí)間的位移與電勢(shì)分布

圖5中(a)、(b)、(c)和(d) 4點(diǎn)分別為步點(diǎn)0、9、14、19 μs，經(jīng)過(guò)后處理后得到該四點(diǎn)的電勢(shì)分布如圖6所示。可以看到在0 μs時(shí)，壓電層電勢(shì)分布整體為0；在9 μs時(shí)電勢(shì)達(dá)到最大值，通過(guò)旁邊圖例可看出其電勢(shì)分布大于另3點(diǎn)；在14和19 μs時(shí)分別處于振蕩的波谷與波峰，因此19 μs時(shí)的電勢(shì)分布明顯大于14 μs。

圖6 分別在0、9、14和19 μs時(shí)壓電層電勢(shì)分布云圖Fig.6 Cloud diagram of piezoelectric layer potential distribution at 0, 9, 14 and 19 μs

圖7 分別在61、66、71和76 μs時(shí)位移圖Fig.7 Displacement diagram at 61, 66, 71 and 76 μs

圖5中(e)、(f)、(g)和(h)4點(diǎn)為步點(diǎn)的61、66、71、76 μs，均處于輸出電壓趨于穩(wěn)定過(guò)程中，通過(guò)后處理對(duì)其機(jī)械狀態(tài)進(jìn)行查看得到其位移圖如圖7所示?？梢钥吹酱烹姀?fù)合材料處于不斷伸縮狀態(tài)中，且由于在趨于穩(wěn)定的過(guò)程中，因此每一個(gè)峰值或波谷的位移均比上一個(gè)峰值或波谷的位移小。

3.3 8/20 μs雷電流激勵(lì)下，F(xiàn)eGa/PZT磁電復(fù)合材料幅值傳感特性分析

接下來(lái)，對(duì)8/20 μs雷電流激勵(lì)下磁電復(fù)合材料的輸入輸出特性進(jìn)行仿真分析，通過(guò)設(shè)置雷電流峰值由0.5 kA增加到10 kA，得到其輸出電壓峰值變化如圖8所示。在0.5 ～10 kA范圍內(nèi)，F(xiàn)eGa/PZT磁電復(fù)合材料輸出電壓峰值Vp與8/20 μs雷電流峰值Ip擬合關(guān)系為

Vp=0.0018×Ip-2.3406,R2=0.9855

(2)

由此可知，F(xiàn)eGa/PZT磁電復(fù)合材料可用于8/20 μs雷電流峰值測(cè)量。

圖8 8/20 μs雷電流激勵(lì)下，磁電復(fù)合材料的幅值傳感特性Fig.8 Amplitude sensing characteristics of magnetoelectric composites driven by 8/20 μs lightning current

4 結(jié) 論

通過(guò)時(shí)域分析研究了在8/20 μs雷電流激勵(lì)下FeGa/PZT磁電復(fù)合材料的瞬態(tài)響應(yīng)特性。在8/20 μs雷電流激勵(lì)下，磁電復(fù)合材料的輸出有較為明顯的振蕩現(xiàn)象，振蕩頻率等于材料一階縱波振動(dòng)頻率(100.36 kHz)。這可以解釋為，F(xiàn)eGa/PZT磁電復(fù)合材料是機(jī)械慣性器件，在脈沖激勵(lì)下，會(huì)發(fā)生阻尼振蕩。在0.5kA到10kA范圍內(nèi)，F(xiàn)eGa/PZT對(duì)8/20 μs雷電流峰值靈敏度達(dá)到1.8 mV/A，線性度達(dá)到98.55%。以上仿真分析結(jié)果證明磁電復(fù)合材料在瞬態(tài)電流測(cè)量領(lǐng)域有良好的應(yīng)用前景。