梁濤,謝彥召,席志豪
(1.西安交通大學(xué)電力設(shè)備電氣絕緣國家重點實驗室,710049,西安;2.西安交通大學(xué)電氣工程學(xué)院,710049,西安)
有意電磁干擾(intentional electromagnetic interference,IEMI)是以干擾、毀傷電子設(shè)備為目的,人為惡意制造的電磁環(huán)境,對設(shè)備的電磁安全構(gòu)成了巨大威脅[1-4]。隨著脈沖功率、輻射天線等技術(shù)的迅速發(fā)展,IEMI干擾源的設(shè)計門檻不斷降低,相關(guān)犯罪案件、恐怖活動、軍事應(yīng)用時有報導(dǎo),引發(fā)了學(xué)界對關(guān)鍵設(shè)備IEMI抗擾度的廣泛關(guān)注[5-8]。
為定量化評價受擾設(shè)備在輻射電磁環(huán)境下的損傷性,前人基于對大量效應(yīng)試驗的總結(jié)歸納,發(fā)現(xiàn)了電子設(shè)備電磁損傷程度同其響應(yīng)波形范數(shù)間的強(qiáng)相關(guān)性,由此發(fā)展了基于范數(shù)的設(shè)備抗擾度量化表征評價體系[9-12]。具體而言,針對響應(yīng)波形f的p范數(shù)定義為
(1)
表1 響應(yīng)波形f的p范數(shù)、對應(yīng)的物理含義及效應(yīng)現(xiàn)象Table 1 p-norm of measurement f,the associated physical quantity and possible susceptibility effect
上述針對HEMP、UWB、HPM等電磁環(huán)境的典型設(shè)備抗擾度研究建立在明確輻射場時域波形的基礎(chǔ)之上,然而,IEMI人為制造的特點決定了其波形更為隨機(jī)、智能,設(shè)計者可利用多種技術(shù),產(chǎn)生不同波形的IEMI(例如,不同頻段窄帶脈沖、雙指數(shù)脈沖、高斯脈沖、正弦衰減波等)[18]。不同輻射IEMI波形作用于電子設(shè)備的耦合效率差異巨大,難以通過窮舉所有可能波形的方法評定系統(tǒng)的響應(yīng)范數(shù)。
為克服輻射波形隨機(jī)多變給設(shè)備IEMI抗擾度評定帶來的挑戰(zhàn)。作者在前序研究中提出了在約束IEMI輻射能量密度的條件下,視其波形為連續(xù)隨機(jī)變量,優(yōu)化研究任意可能波形IEMI作用下受擾設(shè)備響應(yīng)范數(shù)的上確界(損傷程度最大的情況),以此表征設(shè)備IEMI抗擾度。該思路體現(xiàn)了對設(shè)備受IEMI輻照最極端情況的分析,適用于對電磁安全性具有較高要求的電子設(shè)備。依照該思路,前期對整流沖擊(1范數(shù))、能量(2范數(shù))及峰值(∞范數(shù))的范數(shù)定界問題進(jìn)行了研究[19-22],發(fā)現(xiàn)IEMI輻照下系統(tǒng)響應(yīng)的1范數(shù)不存在上確界,窄帶脈沖作用下會發(fā)散至無窮;2范數(shù)上確界在IEMI為調(diào)諧窄帶波形時可獲得,其大小由系統(tǒng)耦合效率函數(shù)諧振峰值決定;波形峰值范數(shù)(∞范數(shù))在IEMI頻譜同系統(tǒng)耦合效率函數(shù)共軛匹配時達(dá)到最大,表現(xiàn)為復(fù)雜的時域?qū)拵Рㄐ?。對比?可知,尚缺失針對斜率峰值、沖擊峰值的范數(shù)定界分析,對∞范數(shù)敏感型目標(biāo)設(shè)備的IEMI抗擾度評定有待深入研究。
由此,為完善IEMI輻照下的設(shè)備響應(yīng)范數(shù)定界方法理論體系,本文進(jìn)一步研究了針對斜率峰值、沖擊峰值的上確界評定問題,通過構(gòu)建泛函優(yōu)化問題,研究對設(shè)備最具毀傷能力的入射場波形特征,補(bǔ)全了面向∞范數(shù)敏感型設(shè)備的IEMI抗擾度評估方法。建立了IEMI輻照設(shè)備模型,通過構(gòu)建并求解等效最優(yōu)化問題,獲得了∞范數(shù)上確界解析解;將方法應(yīng)用于Vivaldi天線及對數(shù)周期天線(LPDA)兩類超寬帶天線,對比分析了最大化斜率峰值、沖擊峰值的IEMI波形特征,研究了范數(shù)上確界評定方法的準(zhǔn)確性。
(2)
(3)
以波形峰值∞范數(shù)為最大化目標(biāo),可將IEMI最優(yōu)波形辨識等價為以下優(yōu)化問題
(4)
式中tp為波形達(dá)到峰值的時刻。文獻(xiàn)[5]利用拉格朗日乘子法獲得了該問題的最優(yōu)解析解,即
(5)
(6)
(7)
(8)
(9)
(10)
式(9)中等號成立的條件是存在常數(shù)C1,使得
(11)
聯(lián)立式(10)和(11),可求得常數(shù)為
(12)
由此,可求得問題式(8)的最優(yōu)解析解為
(13)
VL(t)=
(14)
(15)
(16)
由傅里葉變換性質(zhì)可知
(17)
(18)
Hilbert空間中,由Cauchy-Schwarz不等式可知,電壓波形的沖擊峰值滿足
(19)
式(19)中等號成立的條件是存在常數(shù)C2,使得
(20)
式(20)同約束條件式(10)聯(lián)立,可得常數(shù)C2應(yīng)滿足
(21)
代入式(20)可得所辨識的IEMI最優(yōu)波形頻譜及其負(fù)載響應(yīng)電壓波形分別為
(22)
VL(t)=
(23)
(24)
本節(jié)選取了Vivaldi及LPDA兩類超寬帶天線作為受擾電子設(shè)備[25-26],利用所提范數(shù)定界方法,分析了IEMI輻照受擾天線的極端情況,對比討論了不同范數(shù)目標(biāo)下IEMI最優(yōu)波形的時頻域特征。
Vivaldi、LPDA天線均為典型的超寬帶天線,基于PCB技術(shù)設(shè)計的兩類受擾超寬帶天線結(jié)構(gòu)如圖2所示。選用FR-4環(huán)氧玻璃布層壓板(相對介電常數(shù)為4.3,介質(zhì)損耗角正切為0.025)作為PCB基板材料,基材厚度ts=2 mm,天線貼片及微帶饋線分別位于基板兩側(cè),材質(zhì)均為厚度tt=38 μm、電導(dǎo)率為58 MS/m的退火純銅。
如圖2(a)所示,Vivaldi天線總長度la=41 mm,寬度2lb+s=30 mm。天線開口按照指數(shù)函數(shù)變化規(guī)律延伸至基板末端,貼片層關(guān)于y軸對稱,其中l(wèi)b=14.5 mm,le=31 mm。其他微帶饋線結(jié)構(gòu)參數(shù)為lf=5 mm,ld=3 mm,lp=20 mm。LPDA天線由10根上下兩層交錯布放的陣子組成,結(jié)構(gòu)如圖2(b)所示,由饋源端至遠(yuǎn)端,將陣子由1~10編號,則陣子長度lo1~lo10分別為[3.34,4.08,4.98,6.08,7.41,9.04,11.02,13.45,16.40,20.00] mm;寬度wo1~wo10分別為[0.57,0.69,0.84,1.03,1.26,1.54,1.87,2.29,2.79,3.40] mm;間距分別為[11.97,2.05,2.50,3.05,3.71,4.53,5.53,6.74,8.19,10.03] mm。天線長度lc=70 mm,總寬度ld=50 mm;wf=1 mm。
假設(shè)IEMI是頻譜能量分布在[3 GHz,6 GHz]內(nèi)的任意波形平面電磁場干擾,輻射場的能量密度不超過WE=1 mJ/m2??紤]最大化負(fù)載波形的斜率峰值、沖擊峰值兩種情況,假設(shè)相應(yīng)范數(shù)均在t=td=ti=10 ns時達(dá)到最大。由式(13)、式(22)可分別計算兩種極端耦合情況下所需施加的IEMI最優(yōu)波形頻譜,結(jié)果如圖4所示。
利用逆傅里葉變換,可進(jìn)一步獲得最優(yōu)IEMI電場分量的時域波形,結(jié)果如圖5所示,為突出波形的細(xì)節(jié)差異,僅繪制了[8,12] ns內(nèi)的電場波形,可見,IEMI最優(yōu)時域波形并非雙指數(shù)、高斯等常見簡單波形,而是不斷振蕩衰減的雙極性復(fù)雜波形,針對Vivaldi及LPDA天線的最優(yōu)IEMI波形十分相似,但是細(xì)節(jié)上存在差異,大多能量集中于2~3個振蕩周期內(nèi),波形幅值約為45 kV/m。若考慮不同能量密度約束WE,僅需等比例增減輻射場幅值,無需改變具體時域波形的細(xì)節(jié)特征。
將圖4、圖5中給出的波形作為IEMI激勵場,輻照受擾Vivaldi天線,可得到天線負(fù)載上感應(yīng)電壓時域波形如圖6所示。針對同一天線,斜率峰值/沖擊峰值最優(yōu)波形所感應(yīng)的響應(yīng)電壓均相對于極值時刻t=10 ns奇對稱,且該時刻對應(yīng)電壓波形的過零點。由于LPDA耦合效率高于Vivaldi天線,可發(fā)現(xiàn)響應(yīng)波形的幅值及對應(yīng)∞范數(shù)也明顯增加。
為驗證上述IEMI最優(yōu)波形具備對受擾天線的最強(qiáng)損傷能力,可用于極端條件下的設(shè)備IEMI輻照抗擾度的量化表征。本節(jié)選取了均勻頻譜波形、寬帶高斯波形、窄帶高斯波形3種典型IEMI環(huán)境作為輻射激勵場,通過對比典型IEMI波形及最優(yōu)波形作用下的負(fù)載電壓響應(yīng)波形及其范數(shù),驗證方法的最優(yōu)性。選取的3種典型IEMI電場頻譜分布如圖7所示,其能量密度均歸一化至約束條件,即WE=1 mJ/m2。
具體而言,均勻頻譜波形在[3 GHz,6 GHz]頻段內(nèi)的幅值完全相同,相位均為0;寬帶高斯波形服從文獻(xiàn)[24]中定義,為減小頻譜泄漏的影響,將90%的能量強(qiáng)制約束在所考慮頻段內(nèi)[24];窄帶高斯波形中心頻譜為系統(tǒng)諧振點fc=5 GHz,帶寬比為1.01,其對應(yīng)的時域波形為快速上升并衰減的正弦波形。
表2 不同IEMI激勵下的天線負(fù)載電壓波形∞范數(shù)Table 2 ∞-norm of voltage waveform across antenna load under illumination of different IEMI environment
本文將IEMI輻照電子設(shè)備的范數(shù)定界問題拓展到斜率峰值、沖擊峰值兩類∞范數(shù)敏感型設(shè)備的IEMI輻射抗擾度評估中,通過嚴(yán)格求解,獲得了能最大化設(shè)備相應(yīng)范數(shù)的波形解析解,進(jìn)一步分析發(fā)現(xiàn)以下幾點:
(1)峰值、斜率峰值、沖擊峰值3類∞范數(shù)對應(yīng)的IEMI最優(yōu)波形均由系統(tǒng)電磁耦合效率函數(shù)決定,呈現(xiàn)復(fù)雜寬帶脈沖的波形特征;
(3)針對Vivaldi、LPDA天線等典型寬帶受擾設(shè)備,設(shè)計的寬帶IEMI波形較其他常規(guī)寬帶、窄帶波形更具對設(shè)備的毀傷能力,范數(shù)上確界適合作為表征設(shè)備IEMI損傷性的量化指標(biāo)。
綜上所述,從底線思維角度出發(fā),解決了IEMI波形復(fù)雜多變條件下,無窮范數(shù)敏感性設(shè)備的輻射抗擾度量化評價問題。研究結(jié)論可進(jìn)一步被抗擾度測試、防護(hù)效能評估等研究所借鑒,對于提升關(guān)鍵設(shè)備的IEMI生存能力具有重要基礎(chǔ)理論意義。