蘇慧雯 強鵬 王乙 余俊豪

摘 要：航天器在軌運行時受到高能帶電粒子的作用而引起深層介質(zhì)充放電。當(dāng)高能粒子進入材料深部時，引起充放電過程，對器件造成直接的損傷，同時充放電過程中產(chǎn)生的放電電磁脈沖會干擾航天器電子系統(tǒng)的正常工作。本文根據(jù)麥克斯韋方程組，通過有限差分法對放電脈沖進行了數(shù)值計算，得到由放電脈沖產(chǎn)生的電磁輻射脈沖及放電電場很強，場強最高可達1.8×104V/m。對放電脈沖進行傅里葉變換后獲得電磁脈沖頻率譜，可以得出放電脈沖頻譜較寬，高頻部分頻率最高可達2GHz。由于高頻電磁脈沖穿透力很強，且難以抑制，高頻航天器件空間輻射損傷及防護問題需要進一步研究。

關(guān)鍵詞：深層放電 放電脈沖 分析方法

中圖分類號：V52 文獻標(biāo)識碼：A 文章編號：1674-098X（2017）12（c）-0012-04

航天器在空間運行過程中，會遭受空間帶電粒子輻射、空間微流星體及碎片碰撞和空間氧原子侵蝕等。而對航天器影響最嚴(yán)重的則是空間帶電粒子輻射，其會引起航天器材料損傷、器件性能退化等影響?？臻g帶電粒子輻射會引起單粒子效應(yīng)、總劑量效應(yīng)、充放電效應(yīng)和位移損傷效應(yīng)這4種主要效應(yīng)。空間等離子體與航天器材料表面或內(nèi)部發(fā)生相互作用，使電荷在材料表面或內(nèi)部不斷積累，直至不同部件間或同一部件不同部位形成的電位差超過擊穿閾值時，就會發(fā)生放電現(xiàn)象，放電過程中會釋放出電荷、熱量、電磁脈沖以及輝光等，這就是充放電效應(yīng)。

航天器充放電效應(yīng)又可分為表面充放電效應(yīng)和深層充放電效應(yīng)[1]。其中，高能粒子輻射引起的介質(zhì)材料深層充放電效應(yīng)對于航天器損傷最為嚴(yán)重，其又被稱為衛(wèi)星內(nèi)帶電。深層充電是由高能電子（0.1～10MeV）引起的，主要發(fā)生在材料內(nèi)部。高能帶電粒子入射到材料表面，穿透表面屏蔽層，在材料中沉積，形成深層充電，會導(dǎo)致材料擊穿，進而引起性能的改變；材料中沉積的大量高能粒子后隨即發(fā)生靜電放電現(xiàn)象形成干擾更加強烈的電子系統(tǒng)，同時充放電過程中產(chǎn)生的放電電磁脈沖也會干擾衛(wèi)星的電子系統(tǒng)正常工作，對在軌航天器造成了巨大的傷害。

1 放電電流脈沖

航天器在經(jīng)過深層充電后，介質(zhì)內(nèi)部由于高能電子形成的強電場會隨充電的持續(xù)不斷升高，當(dāng)其超過介質(zhì)可承受的放電閾值時，便會發(fā)生局部放電，而電子在這個過程中受強電場影響，引起雪崩效應(yīng)，產(chǎn)生巨大的放電電流，此即為深層介質(zhì)放電過程。

由于靜電放電對航天器會造成產(chǎn)生的巨大的影響，國際上為分析靜電放電干擾實驗提出了技術(shù)標(biāo)準(zhǔn)IEC61000-4-2[2]。如圖1所示為放電電流脈沖的技術(shù)標(biāo)準(zhǔn)時域圖。

放電脈沖中的快成分在納秒數(shù)量級，其對應(yīng)的頻譜則延伸到GHz頻率范圍，如圖2所示。

由此可以歸納出靜電放電產(chǎn)生的瞬變電磁場主要有以下特性：

（1）脈沖上升很快，并會形成陡峭的上升沿，根據(jù)靜電放電情況的不同而有所不同，通常在ns級。

（2）頻帶比較寬，放電脈沖有從直流分量到GHz的所有頻率分量，也正是因為陡峭的上升沿，使它具有很寬的頻譜分量。

（3）近場幅值很高，在靜電放電發(fā)生的附近，電場可達幾千V/m，放電電流可達10A數(shù)量級，近場磁場可達毫特斯拉數(shù)量級。

（4）脈沖波形具有較長的衰減尾沿，說明放電頻譜中存在極低頻率的成分，脈沖持續(xù)時間長，對電子器件的充電時間長，對電子器件帶來的損傷會增大。

2 放電電磁脈沖數(shù)值的計算

在充電過程中，電子在介質(zhì)中積累并在其中形成特殊的空間分布，在這種空間電荷分布下會形成特定的電場分布。當(dāng)入射電子通量不變時，充電形成的電場、空間電荷分布最終會達到一個動態(tài)平衡，即達到穩(wěn)態(tài)。當(dāng)存在擾動或材料缺陷等因素時，充電形成的穩(wěn)態(tài)會被打破，進而發(fā)生放電過程，形成放電電流。在電荷分布，電流的共同作用下，形成電磁脈沖。

對電流脈沖產(chǎn)生電磁脈沖的可根據(jù)麥克斯韋電磁場理論進行求解。在電荷分布、電流分布及介質(zhì)參數(shù)和邊界條件確定的情況下，就能夠準(zhǔn)確獲得放電產(chǎn)生的電磁脈沖。

麥克斯韋方程組的積分形式和微分形式如式（1）和式（2）所示。

（1）

（2）

麥克斯韋方程的數(shù)值求解有多種方法，常用方法有矩量法、有限元法和時域有限差分法。本文采用時域有限差分法對的麥克斯韋方程進行差分離散化，利用蛙跳算法對空間領(lǐng)域內(nèi)的電場和磁場進行交替計算，進而迭代求解出電磁場見圖3。

利用時域有限差分法推導(dǎo)出直角坐標(biāo)系中麥克斯韋方程組的的差分方程[3]如式（3）所示。

（3）

其中Δx，Δy，Δz為空間步長，Δt為時間步長。

平板介質(zhì)模型如圖4所示，介質(zhì)版厚度1.5mm。高能電子垂直入射進入介質(zhì)平板時，將入射背面接地。空間步長取0.15mm。

經(jīng)過一段時間的充電后，介質(zhì)深層充電電場逐漸增強。當(dāng)充電電場強度超過放電閾值時，會產(chǎn)生快電流激勵，電流激勵方程如式（4）所示。

（4）

放電時域波形如圖5所示。

從圖6中可以看出，當(dāng)放電開始發(fā)生時，介質(zhì)電導(dǎo)率快速變大；當(dāng)放電電流達到峰值后，可近似認(rèn)為電導(dǎo)率不再改變。由時域有效差分可以求得電磁脈沖分布：

該放電脈沖是包含多種頻率的電磁脈沖，脈沖上升時間在10ns之間，脈沖衰減時間則相對較長。通過數(shù)值分析，可以得到深層充電引發(fā)的放電過程包含多種頻率的電磁輻射，其高頻成分可以到達GHz頻率范圍。

對放電脈沖進行傅里葉變換后得到放電電磁脈沖頻譜如圖7所示。

從圖7中可以看出，放電過程中釋放出來的電磁脈沖包含多種成分，空間電子器件在電磁脈沖的作用下會出現(xiàn)擾動、軟錯誤及損傷，特別是近場區(qū)域的電子器件影響更大。

3 結(jié)語

本文以Maxwell方程為基礎(chǔ)，利用有限元差分方法編程計算放電過程產(chǎn)生的電磁脈沖，通過差分計算獲得了放電產(chǎn)生的電磁脈沖波形，在此基礎(chǔ)上通過傅里葉變換得出放電電磁脈沖的頻譜，分析了電磁脈沖的頻譜特征。分析得出：傅里葉變換后得到的放電電磁脈沖頻譜中，高頻部分頻率很高，最高可達2GHz；放電過程中形成的電場強度極強，最高可達1.8×104V/m。由于高頻成份穿透力極強，且難以抑制，所以隨著航天技術(shù)的發(fā)展，許多航天器件工作的頻率向高頻方向延伸，開展充放電過程中電磁脈沖的分析研究很有意義。

參考文獻

[1] 秦曉剛.介質(zhì)深層帶電數(shù)值模擬與應(yīng)用研究[D].蘭州大學(xué)，2010.

[2] Hsing-Yi Chen，Yu-Ching Chu，ESD currents induced in a charged human body approaching a vehicle[J].Journal of Electrostatics，2011（69）：54-59.

[3] 王依剛，陳彬，石立華，等.任意形狀金屬腔體內(nèi)電場脈沖的三維計算[J].強激光與粒子束，2012，24（5）：1141-1145.

[4] 全榮輝.航天器介質(zhì)深層充放電特征及其影響[D].中國科學(xué)院，2009.

[5] 方美華.深空輻射粒子在介質(zhì)材料中的輸運及損傷研究[D].南京航空航天大學(xué)，2011.

[6] 胡紅波，韓金林.宇宙線的百年研究[J].物理，2012（2）：111-113.