做開(kāi)關(guān)電器、等離子體焊槍或者高壓斷路器設(shè)計(jì)的朋友應(yīng)該都有同一個(gè)感受電弧是工程問(wèn)題里最復(fù)雜、最棘手、也最迷人的物理現(xiàn)象之一。一個(gè)小小的放電通道溫度輕輕松松上萬(wàn)K電流密度、氣流速度、熱輻射和材料燒蝕全部擠在一個(gè)毫米甚至亞毫米級(jí)的空間里互相糾纏。早年我在學(xué)校做高壓開(kāi)關(guān)的仿真課題時(shí)被導(dǎo)師要求用COMSOL搭建一個(gè)“多物理場(chǎng)電弧模型”一開(kāi)始我以為只是把電流場(chǎng)、流場(chǎng)、溫度場(chǎng)同時(shí)算一遍后來(lái)才明白真正的難點(diǎn)在于磁流體動(dòng)力學(xué)MHD耦合、瞬態(tài)電弧演化以及電極燒蝕深度的準(zhǔn)確定量。這幾年做過(guò)的電弧仿真項(xiàng)目中COMSOL是少數(shù)能把“電-磁-熱-流-固”這幾塊密耦合放在一個(gè)環(huán)境里統(tǒng)一求解的工具配合MATLAB做聯(lián)合仿真還能反向驅(qū)動(dòng)參數(shù)掃描和結(jié)果分析。這篇文章就把我從模型搭建、物理場(chǎng)選擇、邊界條件設(shè)計(jì)、瞬態(tài)求解到燒蝕深度提取的完整經(jīng)驗(yàn)和踩坑記錄整理出來(lái)希望能給正在搞電弧仿真、等離子體仿真或者開(kāi)關(guān)設(shè)備研發(fā)的同行一些參考。1. 電弧模型的整體設(shè)計(jì)與物理場(chǎng)耦合思路1.1 多物理場(chǎng)耦合的基本邏輯很多人第一次接觸電弧仿真習(xí)慣先把它當(dāng)成“一個(gè)電場(chǎng)問(wèn)題”去思考這個(gè)起點(diǎn)就偏了。電弧并不是單純的電流傳導(dǎo)它是一個(gè)閉環(huán)的自持放電過(guò)程外加電場(chǎng)讓氣體電流導(dǎo)通電流產(chǎn)生焦耳熱焦耳熱把氣體溫度快速抬升到上萬(wàn)K高溫氣體的電導(dǎo)率隨之飆升電導(dǎo)率升高又讓電流密度進(jìn)一步集中電流通道變窄之后焦耳熱更加集中這種正反饋會(huì)在極短時(shí)間內(nèi)形成一條高溫、高導(dǎo)電率的等離子體通道也就是我們看到的電弧。與此同時(shí)電流在自身電磁場(chǎng)中會(huì)受到洛倫茲力洛倫茲力驅(qū)動(dòng)電弧通道內(nèi)等離子體流動(dòng)形成強(qiáng)烈的射流這個(gè)“電流-磁場(chǎng)-流動(dòng)”的相互作用才是磁流體動(dòng)力學(xué)的本質(zhì)。把這條因果鏈翻譯成COMSOL的物理場(chǎng)節(jié)點(diǎn)標(biāo)準(zhǔn)的映射是這樣的電場(chǎng)模塊求解電流守恒方程得到電勢(shì)和電流密度分布焦耳熱項(xiàng)作為能量方程的熱源驅(qū)動(dòng)溫度場(chǎng)升高溫度場(chǎng)通過(guò)溫度相關(guān)的電導(dǎo)率、導(dǎo)熱系數(shù)重新影響電場(chǎng)分布電流密度和磁矢勢(shì)耦合出洛倫茲力洛倫茲力作為體積力進(jìn)入流體動(dòng)量方程流體流動(dòng)改變熱量輸運(yùn)電弧通道被吹偏、收縮或者拉伸高溫金屬表面發(fā)生熔化、蒸發(fā)、飛濺形成燒蝕坑燒蝕形貌再反過(guò)頭來(lái)改變電場(chǎng)邊界和電流密度分布。這個(gè)閉環(huán)耦合看起來(lái)復(fù)雜但COMSOL的模塊化結(jié)構(gòu)恰好適合這種多物理場(chǎng)聯(lián)立求解。我習(xí)慣把模型拆成幾條“物理總線”每一條都用單獨(dú)的物理場(chǎng)接口去定義最后通過(guò)多物理場(chǎng)耦合節(jié)點(diǎn)把它們綁在一起。這樣做的最大好處是任何一個(gè)環(huán)節(jié)出問(wèn)題都能單獨(dú)調(diào)試不會(huì)出現(xiàn)一個(gè)報(bào)錯(cuò)導(dǎo)致整個(gè)系統(tǒng)無(wú)從下手。1.2 為什么選COMSOL而不是Fluent或完全自編程其實(shí)同行之間討論電弧仿真經(jīng)常繞不開(kāi)另一個(gè)工業(yè)流體軟件Fluent還有很多課題組選擇自己寫(xiě)MHD求解器。我的實(shí)際體感是Fluent在純流場(chǎng)和燃燒計(jì)算上有優(yōu)勢(shì)但要處理電弧里的“電路驅(qū)動(dòng)”需要自己寫(xiě)好幾份UDF去實(shí)現(xiàn)電流連續(xù)方程、電勢(shì)求解和洛倫茲力源項(xiàng)而且在迭代過(guò)程中UDF與流場(chǎng)的耦合穩(wěn)定性全看個(gè)人功底。完全自編程就更不必說(shuō)了光是搞定非結(jié)構(gòu)網(wǎng)格上的電磁場(chǎng)求解器、DAT文件交換和前處理已經(jīng)足夠消磨掉大半年時(shí)間。COMSOL的優(yōu)勢(shì)在于“物理場(chǎng)接口之間的耦合是原生的”。比如我可以直接在材料屬性里定義電導(dǎo)率為溫度的函數(shù)sigma(T)電場(chǎng)接口自動(dòng)讀取當(dāng)前溫度場(chǎng)去更新導(dǎo)電特性洛倫茲力可以作為一個(gè)體積力表達(dá)式直接寫(xiě)入Navier-Stokes方程的源項(xiàng)不需要額外寫(xiě)任何接口程序。“Flux pinning”式的多層物理場(chǎng)管理讓復(fù)雜耦合的建模過(guò)程保持清晰調(diào)試迭代時(shí)可以隨時(shí)打開(kāi)某個(gè)物理場(chǎng)節(jié)點(diǎn)查看源項(xiàng)是否正確這對(duì)快速驗(yàn)證物理假設(shè)至關(guān)重要。當(dāng)然COMSOL也有自己的短板。電弧等離子體實(shí)際上處于局部熱力學(xué)平衡與非平衡的過(guò)渡區(qū)尤其靠近電極邊界層需要更細(xì)致的碰撞-輻射模型COMSOL默認(rèn)接口不會(huì)覆蓋這么深。所以我會(huì)把COMSOL定位成“工程仿真驗(yàn)證平臺(tái)”而不是“等離子體物理研究全套工具”需要精確電子溫度和化學(xué)反應(yīng)機(jī)理的研究人員通常會(huì)把COMSOL計(jì)算出的宏觀熱流場(chǎng)和流場(chǎng)作為輸入再拿專(zhuān)門(mén)等離子體程序或?qū)嶒?yàn)數(shù)據(jù)去校準(zhǔn)近壁面區(qū)域。這個(gè)定位要提前想清楚否則你在細(xì)節(jié)上糾結(jié)太久項(xiàng)目進(jìn)度就完全失控了。1.3 模型維度與計(jì)算代價(jià)的權(quán)衡電弧模型到底應(yīng)該建2D還是建3D這個(gè)問(wèn)題幾乎每三個(gè)同行就會(huì)問(wèn)一次。我的經(jīng)驗(yàn)非常明確如果電極結(jié)構(gòu)是軸對(duì)稱(chēng)的比如常見(jiàn)的觸頭結(jié)構(gòu)、TIG焊槍噴嘴、均勻介質(zhì)間隙那就老老實(shí)實(shí)先用2D軸對(duì)稱(chēng)做一切前期研究。2D軸對(duì)稱(chēng)模型的網(wǎng)格量只有3D模型的幾十分之一求解速度大幅提升特別適合做“參數(shù)掃描初篩”和“物理模型校準(zhǔn)”。只有在明確存在強(qiáng)非對(duì)稱(chēng)因素的情況下才建議升級(jí)到3D模型。典型場(chǎng)景包括電弧受到外部橫向磁場(chǎng)驅(qū)動(dòng)發(fā)生彎曲、分閘過(guò)程中電弧被氣流橫吹、電極表面存在隨機(jī)粗糙度導(dǎo)致燒蝕點(diǎn)偏移。這種三維仿真不是簡(jiǎn)單把幾何拉伸一下就完事你需要額外關(guān)注幾個(gè)問(wèn)題3D磁矢勢(shì)方程在COMSOL中的收斂性會(huì)差很多建議使用“Coil”接口或直接耦合磁場(chǎng)邊界條件3D網(wǎng)格規(guī)模容易突破內(nèi)存瓶頸電弧核心區(qū)域網(wǎng)格尺寸需要控制在毫米以下一整片放電空間剖分下來(lái)幾百萬(wàn)單元是常態(tài)瞬態(tài)求解的壓力基OCE/BDF組合在3D狀況下時(shí)間步進(jìn)極慢最好在2D模型上確認(rèn)好物理參數(shù)和時(shí)間步長(zhǎng)再遷移到3D做最終驗(yàn)證。我踩過(guò)的坑就是貪圖“一步到位”第一次直接建3D斷路器燃弧室模型結(jié)果網(wǎng)格剖分花了三天第一次求解在第二十分鐘就發(fā)散了。后來(lái)老老實(shí)實(shí)回到2D軸對(duì)稱(chēng)把動(dòng)網(wǎng)格、邊界條件、材料參數(shù)逐一驗(yàn)證完畢再去擴(kuò)展到3D才順利跑通。給沒(méi)經(jīng)驗(yàn)的同行一個(gè)建議先2D后3D不要讓前期的求知欲在上機(jī)時(shí)變成折磨。2. 核心物理場(chǎng)搭建與實(shí)操細(xì)節(jié)2.1 電磁場(chǎng)模塊的方程與控制在COMSOL中搭建電弧模型的電磁場(chǎng)環(huán)境最基礎(chǔ)的是電流守恒方程。電弧等離子體被近似為導(dǎo)體電勢(shì)分布滿(mǎn)足?·(σ?φ)0其中σ是電導(dǎo)率。表征氣體從常溫絕緣體到上萬(wàn)K導(dǎo)電體的轉(zhuǎn)變時(shí)電導(dǎo)率的值會(huì)相差十幾個(gè)數(shù)量級(jí)。空氣在常溫下電導(dǎo)率可以低到10?1?量級(jí)而在20000K等離子體狀態(tài)下電導(dǎo)率高達(dá)幾千S/m。這么大的跨數(shù)量級(jí)跨度如果直接采用物理值求解方程矩陣的條件數(shù)會(huì)惡化到讓求解器徹底失效。所以我建議在初始階段為電導(dǎo)率設(shè)置一個(gè)邊界下限例如取1e-3 S/m作為數(shù)值截?cái)嚅撝当苊獗尘皻怏w區(qū)域被電場(chǎng)方程強(qiáng)行當(dāng)成絕緣體而引發(fā)數(shù)值奇異問(wèn)題。磁場(chǎng)的引入上COMSOL中“磁場(chǎng)”接口默認(rèn)求解磁矢勢(shì)A。在二維軸對(duì)稱(chēng)和三維情況下通過(guò)如下關(guān)系得到磁通密度B ?×A實(shí)際計(jì)算中電磁場(chǎng)和流場(chǎng)通過(guò)洛倫茲力F J×B相互作用。這里的J由電流守恒方程直接給出。需要注意在高頻或強(qiáng)瞬變條件下趨膚效應(yīng)會(huì)影響電流分布但如果電弧頻率不太高準(zhǔn)靜態(tài)近似已經(jīng)足夠這也是大多數(shù)開(kāi)關(guān)電器電弧仿真采納的簡(jiǎn)化假設(shè)。2.2 磁流體動(dòng)力學(xué)的實(shí)現(xiàn)方法分離求解還是全耦合這里應(yīng)該是全文最重要、也最容易被卡住的一個(gè)技術(shù)分歧點(diǎn)。COMSOL里的多物理場(chǎng)求解器有兩種策略全耦合求解和分離式求解。初學(xué)者往往默認(rèn)選擇全耦合因?yàn)榭雌饋?lái)最“嚴(yán)謹(jǐn)”它把所有的待求變量放在一個(gè)巨大的雅可比矩陣?yán)镆淮吻蠼?。但電弧的物理?chǎng)之間非線性極強(qiáng)全耦合矩陣的規(guī)模會(huì)迅速爆炸在時(shí)間步進(jìn)過(guò)程中還極易出現(xiàn)“殘差振蕩-發(fā)散”的惡性循環(huán)。我自己經(jīng)過(guò)大量對(duì)比后穩(wěn)定使用的策略是分離式求解。典型做法是把變量分成三組電場(chǎng)與磁場(chǎng)變量電勢(shì)φ、磁矢勢(shì)A為一組流體流動(dòng)變量速度u、壓力p為一組能量變量溫度T為一組。每個(gè)時(shí)間步內(nèi)按順序求解這三組并在每一組之間執(zhí)行一定次數(shù)的迭代循環(huán)直到整體殘差滿(mǎn)足容差。這種分離式策略的好處是內(nèi)存消耗低而且當(dāng)某個(gè)物理場(chǎng)發(fā)散時(shí)可以從對(duì)應(yīng)的子步驟中快速定位問(wèn)題而不是面對(duì)一個(gè)全耦合矩陣無(wú)從下手。代價(jià)是收斂速度在強(qiáng)耦合場(chǎng)景下會(huì)比全耦合多耗一些迭代次數(shù)但由于每一步的規(guī)模小整體時(shí)間成本仍然可接受。實(shí)操時(shí)我還會(huì)在分離求解器中手動(dòng)調(diào)整迭代次數(shù)上限。默認(rèn)設(shè)置往往只顧“最大迭代次數(shù)”卻沒(méi)有考慮單元質(zhì)量較差的區(qū)域會(huì)在某次迭代中持續(xù)震蕩從而浪費(fèi)大量時(shí)間。通過(guò)實(shí)驗(yàn)對(duì)比把“Maximum Number of Iterations”降到合理區(qū)間并配合“Termination Criterion”可以讓整體求解時(shí)間縮短20%以上。2.3 邊界條件與電弧初始種子電弧仿真最難處理的兩類(lèi)邊界電極-氣體界面和開(kāi)放邊界。電極表面既是電流進(jìn)出邊界又是強(qiáng)熱流作用的界面還要考慮材料熔化蒸發(fā)。常見(jiàn)做法是在電極邊界施加電勢(shì)邊界條件而氣體區(qū)域的上端接地或接地為0電位。在COMSOL流體接口中開(kāi)放邊界應(yīng)謹(jǐn)慎使用“Open Boundary”在電弧射流強(qiáng)烈向外擴(kuò)展的情況下建議使用具有一定長(zhǎng)度的“壓力遠(yuǎn)場(chǎng)”區(qū)并設(shè)置壓力遠(yuǎn)場(chǎng)邊界條件配合弱約束以抑制回流在出口處產(chǎn)生的數(shù)值振蕩。初始種子問(wèn)題也值得展開(kāi)講。電弧無(wú)法在完全低溫絕緣的氣體中憑空建立因?yàn)槌跏茧妼?dǎo)率太低的區(qū)域電流根本無(wú)法注入。仿真中常見(jiàn)做法是在兩個(gè)電極之間的某個(gè)微小區(qū)域預(yù)先設(shè)置一個(gè)“高溫高導(dǎo)電種子區(qū)”通常做法是在初始值或初始階段直接寫(xiě)入一個(gè)局部溫度分布例如20000K的高斯分布區(qū)域。這個(gè)種子區(qū)域模擬了真實(shí)的擊穿起始過(guò)程比如觸頭分離瞬間的場(chǎng)致發(fā)射形成的初始導(dǎo)電通道讓電流有了流通路徑焦耳熱才能順勢(shì)把放電通道維持和擴(kuò)展起來(lái)。如果偷懶不設(shè)種子電流無(wú)法導(dǎo)通整個(gè)模型就永遠(yuǎn)是“空氣泡里的絕緣體”怎么迭代都點(diǎn)不燃。對(duì)于觸頭分離的瞬態(tài)過(guò)程有時(shí)還需要用移動(dòng)邊界來(lái)模擬觸頭之間距離隨時(shí)間增加。這種情況下“移動(dòng)網(wǎng)格”接口與“變形幾何”的功能就派上了用場(chǎng)這是后面第3節(jié)要展開(kāi)的內(nèi)容。2.4 溫度相關(guān)的材料屬性管理電弧模型的材料屬性不是固定常量而是強(qiáng)溫度非線性函數(shù)?？諝饣騍F6氣體在不同溫度下的電導(dǎo)率、導(dǎo)熱系數(shù)和比熱容往往跨越好幾個(gè)數(shù)量級(jí)。COMSOL允許設(shè)置“插值函數(shù)”讀取外部實(shí)驗(yàn)數(shù)據(jù)或求解出的熱力學(xué)表格數(shù)據(jù)。我最常用的做法是把來(lái)自文獻(xiàn)分段的“等離子體熱物理屬性表”整理成文本文件通過(guò)“Interpolation”節(jié)點(diǎn)導(dǎo)入再在材料屬性中引用這些函數(shù)。這里必須提醒一個(gè)問(wèn)題當(dāng)你把溫度相關(guān)物理屬性施加到整個(gè)計(jì)算域時(shí),會(huì)因?yàn)楸尘皻怏w區(qū)域的溫度從常溫到高溫變化導(dǎo)致物性表的外推不準(zhǔn)確。常溫區(qū)電導(dǎo)率近似為0導(dǎo)熱系數(shù)顯著較低如果函數(shù)插值不到位會(huì)引起低溫區(qū)的虛假加熱或散熱。建議在函數(shù)定義時(shí)設(shè)置好上邊界和下邊界讓低溫區(qū)穩(wěn)定趨向飽和值而不是無(wú)限外推導(dǎo)致數(shù)值溢出。我通常會(huì)把300K以下和30000K以上的數(shù)值用固定的截?cái)嘀当Ｗo(hù)起來(lái)確保求解器面對(duì)的永遠(yuǎn)是一個(gè)連續(xù)、可微的屬性函數(shù)。3. 瞬態(tài)電弧與移動(dòng)網(wǎng)格技術(shù)3.1 瞬態(tài)過(guò)程的物理必要性不少人圖省事直接把電弧仿真做成穩(wěn)態(tài)或者只做簡(jiǎn)單的參數(shù)掃描。但電弧的最大特點(diǎn)就是“瞬態(tài)性”——擊穿的發(fā)生只在微秒量級(jí)電弧通道的擴(kuò)張和收縮在毫秒量級(jí)電極燒蝕材料的蒸發(fā)和熔滴飛濺則持續(xù)貫穿整個(gè)過(guò)程。如果我們把模型簡(jiǎn)化成穩(wěn)態(tài)等于默認(rèn)電弧從一開(kāi)始就處于一種“理想的平衡狀態(tài)”這會(huì)完全忽略電壓恢復(fù)過(guò)程、電流過(guò)零重燃、電弧熄滅-重燃的復(fù)雜動(dòng)態(tài)。開(kāi)關(guān)電器里的實(shí)際應(yīng)用場(chǎng)景更是如此。斷路器開(kāi)斷故障電流時(shí)需要知道電弧在弧觸頭間是如何從擊穿、發(fā)展到大電流電弧再被氣流吹熄的每一步的時(shí)間尺度都不同。瞬態(tài)仿真的意義正是在這些動(dòng)態(tài)過(guò)程中捕捉電流密度和能量密度在空間上的實(shí)時(shí)分布為燒蝕深度和滅弧室設(shè)計(jì)的評(píng)估提供依據(jù)。3.2 COMSOL移動(dòng)網(wǎng)格ALE與單元活化電弧仿真中移動(dòng)網(wǎng)格至少有兩個(gè)典型應(yīng)用。第一個(gè)是電極間距變化真實(shí)斷路器中觸頭在分閘過(guò)程一直在拉開(kāi)電弧長(zhǎng)度不斷增加。第二個(gè)是電極表面的熔化燒蝕高溫作用下電極表面材料發(fā)生燒蝕而形成凹坑形貌隨過(guò)程演化。這兩個(gè)場(chǎng)景都需要用“變形幾何”接口處理幾何形狀隨時(shí)間的變化。COMSOL的移動(dòng)網(wǎng)格Moving Mesh / Deformed Geometry通過(guò)求解Laplace平滑方程控制節(jié)點(diǎn)位移。實(shí)現(xiàn)時(shí)要注意邊界區(qū)域的網(wǎng)格畸變當(dāng)電極形貌變化幅度很大時(shí)網(wǎng)格會(huì)出現(xiàn)折疊和穿透導(dǎo)致單元質(zhì)量降到負(fù)值求解崩潰。遇到這種情況有幾種處理方案限制形變上限每次求解的時(shí)間步不長(zhǎng)讓形變?cè)隽勘３州^小采用“Automatic Remeshing”選項(xiàng)當(dāng)網(wǎng)格質(zhì)量低于設(shè)定閾值時(shí)自動(dòng)重新剖分并映射解用單元活化/去活化技術(shù)在電弧增材制造或電極燒蝕建模中把發(fā)生燒蝕的單元做“生死單元”處理使其退出計(jì)算域同時(shí)讓下方與之接觸的新單元被激活參與后續(xù)計(jì)算。移動(dòng)網(wǎng)格和單元活化這兩個(gè)技術(shù)搭配使用幾乎可以應(yīng)對(duì)所有電弧燒蝕形貌演化問(wèn)題。但要記住COMSOL中的自動(dòng)重剖分不是魔法它在大應(yīng)變時(shí)會(huì)引入插值誤差重剖分前后的守恒量也可能有微小損失。我的經(jīng)驗(yàn)是如果燒蝕深度達(dá)到毫米級(jí)別盡量在變形幾何中把最大單元位移限制在單元尺寸的1/3以?xún)?nèi)寧可多加幾個(gè)重剖分步也不要貪大引起的精度震蕩。3.3 時(shí)間步進(jìn)與CFL條件的工程把握瞬態(tài)求解器的時(shí)間步長(zhǎng)設(shè)定直接決定仿真能不能跑完。COMSOL默認(rèn)的自適應(yīng)時(shí)間步進(jìn)一般比較保守但電弧模型中強(qiáng)對(duì)流傳熱和快速電磁響應(yīng)組合會(huì)讓默認(rèn)時(shí)間步長(zhǎng)積極“亂跳”。我的建議是給時(shí)間步長(zhǎng)設(shè)置一個(gè)明確上限通常取網(wǎng)格單元尺寸與典型流體速度比值的若干分之一也就是CFL條件。具體計(jì)算方式可以這樣估算電弧核心區(qū)域網(wǎng)格尺寸假設(shè)為0.1mm等離子體射流速度可能達(dá)到每秒幾百到幾千米對(duì)應(yīng)的真實(shí)時(shí)間步長(zhǎng)上限大概在0.1微秒到1微秒級(jí)別?？此茦O小的步長(zhǎng)實(shí)際電流建立過(guò)程的特征時(shí)間也確實(shí)在這個(gè)量級(jí)。如果把時(shí)間步長(zhǎng)設(shè)置太大流場(chǎng)中就會(huì)出現(xiàn)明顯的偽振蕩溫度場(chǎng)也會(huì)出現(xiàn)“跳變式”的非物理分布。在COMSOL求解器設(shè)置中我會(huì)使用BDFBackward Differentiation Formula方法階數(shù)設(shè)為2這能在穩(wěn)定性和精度之間取得較好的平衡。更高階數(shù)的BDF在非線性極強(qiáng)的電弧場(chǎng)景中反而容易出現(xiàn)過(guò)沖讓溫度場(chǎng)局部失真。另一點(diǎn)實(shí)際經(jīng)驗(yàn)是如果發(fā)現(xiàn)某段時(shí)間殘差曲線振蕩得很厲害第一反應(yīng)不是調(diào)低默認(rèn)容差而是試著把時(shí)間步長(zhǎng)限制得更小通常效果立竿見(jiàn)影。4. 電場(chǎng)燒蝕深度的計(jì)算與后處理4.1 燒蝕的物理模型與簡(jiǎn)化電弧對(duì)電極表面造成的“電場(chǎng)燒蝕”其實(shí)是多個(gè)物理過(guò)程疊加的結(jié)果高溫引起表層材料的固態(tài)加熱、熔化過(guò)程吸收潛熱、熔融金屬在電磁力作用下被吹離電極表面、蒸發(fā)氣體形成金屬蒸氣進(jìn)入等離子體弧柱、濺射與飛濺進(jìn)一步帶走材料。把所有物理機(jī)制都顯式建模是不現(xiàn)實(shí)的工程上通常用兩個(gè)層次的簡(jiǎn)化模型。第一層次是純熱傳導(dǎo)燒蝕模型忽略熔化和飛濺細(xì)節(jié)把電極看作固體材料表面接收到電弧輸入的熱流密度表面溫度達(dá)到熔點(diǎn)或沸點(diǎn)的深度隨時(shí)間演化即得到燒蝕深度。這個(gè)模型簡(jiǎn)單適合工程估算。第二層次是考慮相變和移動(dòng)邊界使用COMSOL的“固體-流體相變”接口把熔化潛熱作為等效熱容加入能量方程同時(shí)配合變形幾何讓熔化后的材料邊界隨熱流繼續(xù)后退或熔池振蕩。這個(gè)模型更接近真實(shí)但計(jì)算量大得多。針對(duì)項(xiàng)目目標(biāo)——探究電場(chǎng)燒蝕深度——我更推薦“熱穿透深度估算相變修正”的組合方案。熱穿透深度可以粗略通過(guò)以下關(guān)系估算d ≈ √(4αt)其中α是熱擴(kuò)散系數(shù)t是電弧加熱持續(xù)時(shí)間。通過(guò)溫度云圖和邊界熱流密度準(zhǔn)確確定加熱持續(xù)時(shí)間和熱量輸入再用熱穿透深度估算燒蝕量級(jí)再結(jié)合實(shí)驗(yàn)結(jié)果修正系數(shù)。這種方式比起完全大動(dòng)干戈做相變仿真不僅在項(xiàng)目前中期更省時(shí)而且更容易在標(biāo)定后達(dá)到工程精度。4.2 燒蝕深度后處理的三類(lèi)指標(biāo)在COMSOL中提取燒蝕深度時(shí)不要只盯著“表面最高溫度”這一項(xiàng)指標(biāo)?？煽壳腋咝У奶崛》绞接腥齻€(gè)溫度場(chǎng)穿透深度在后處理中沿電極材料厚度方向創(chuàng)建一條截線觀察溫度分布達(dá)到熔點(diǎn)的深度這個(gè)深度直接給出了潛在燒蝕區(qū)域的最大范圍熱流累積量對(duì)邊界熱流密度做時(shí)間積分得到能量密度單位為J/mm2再用能量密度除以材料體積熔化/蒸發(fā)所需單位體積能量估算燒蝕質(zhì)量表面形貌位移在采用移動(dòng)網(wǎng)格的模型中通過(guò)提取邊界節(jié)點(diǎn)的總位移量直接得到燒蝕坑的瞬態(tài)深度曲線。這里給一個(gè)核心經(jīng)驗(yàn)燒蝕深度不要只看最終時(shí)刻的數(shù)值一定要繪制“時(shí)間-最大深度”曲線。因?yàn)殡娀崃魇撬矐B(tài)的燒蝕深度往往存在明顯的“飽和趨勢(shì)”到后期材料表面達(dá)到熱平衡后深度增長(zhǎng)變緩。只有理解了這條曲線的變化規(guī)律才能正確評(píng)估電極材料耐受性和滅弧室壽命。4.3 MHD耦合對(duì)燒蝕深度的重要性燒蝕深度為什么會(huì)和磁流體動(dòng)力學(xué)扯上關(guān)系關(guān)鍵在于電弧收縮效應(yīng)。電弧電流通過(guò)自身磁場(chǎng)產(chǎn)生的洛倫茲力會(huì)強(qiáng)烈擠壓電弧通道讓電流密度集中在中心位置形成極高的局部熱流密度。洛倫茲力同時(shí)驅(qū)動(dòng)等離子體射流沿軸向高速?zèng)_刷電極表面這種沖擊力會(huì)加強(qiáng)熔融金屬的飛濺和機(jī)械蝕除。因此一個(gè)只考慮溫度場(chǎng)和電場(chǎng)的簡(jiǎn)化模型很可能?chē)?yán)重低估電極燒蝕深度。加入MHD耦合后電弧通道收縮、射流沖擊以及熔池飛濺被一并納入熱流密度集中在更小的區(qū)域燒蝕深度自然大幅提升。實(shí)際工程中電弧仿真模型如果不把MHD的收縮效應(yīng)考慮進(jìn)去算出的燒蝕深度能比實(shí)驗(yàn)結(jié)果低一個(gè)數(shù)量級(jí)。這也是“磁流體動(dòng)力學(xué)”在電弧仿真中處于不可替代位置的最直接原因。我在自己項(xiàng)目里的做法是先用無(wú)MHD的簡(jiǎn)化模型跑一批估算結(jié)果再打開(kāi)MHD耦合重新仿真兩組結(jié)果做對(duì)比用來(lái)衡量MHD效應(yīng)在特定工況下的貢獻(xiàn)比例。如果貢獻(xiàn)低比如電弧電流比較小、電弧通道較寬簡(jiǎn)化模型還能用如果貢獻(xiàn)高大電流、細(xì)電弧通道就必須以MHD結(jié)果為準(zhǔn)。4.4 結(jié)果驗(yàn)證與誤差控制仿真結(jié)果拿到手千萬(wàn)別直接信。電弧的物理過(guò)程極其復(fù)雜數(shù)值模型必然存在誤差必須做“對(duì)標(biāo)-修正”閉環(huán)。至少要做三件事網(wǎng)格收斂性研究對(duì)關(guān)鍵工況分別用粗、中、細(xì)三套網(wǎng)格計(jì)算燒蝕深度如果三套結(jié)果差距太大說(shuō)明當(dāng)前網(wǎng)格不夠差距小說(shuō)明網(wǎng)格對(duì)結(jié)果的依賴(lài)已經(jīng)降到可接受范圍熱流密度與總電流校驗(yàn)把仿真結(jié)果中的總電流與設(shè)定的外部電路電流對(duì)比若偏差超過(guò)幾個(gè)百分點(diǎn)往往是邊界條件設(shè)置或材料電導(dǎo)率取值的問(wèn)題實(shí)驗(yàn)對(duì)標(biāo)如果條件允許把預(yù)測(cè)的燒蝕深度和實(shí)測(cè)的觸頭質(zhì)量損失或坑深做對(duì)比建立修正系數(shù)這個(gè)修正系數(shù)在后續(xù)同類(lèi)型仿真中可以直接復(fù)用。最后再提一點(diǎn)后處理層面的事COMSOL導(dǎo)出燒蝕深度數(shù)據(jù)時(shí)建議使用“截線-投影”的方式保存為CSV文件然后用MATLAB或Python做后續(xù)曲線擬合而不是直接在COMSOL界面里讀數(shù)值。這樣既能保留完整的過(guò)程數(shù)據(jù)也便于在不同項(xiàng)目之間橫豎對(duì)比。5. 聯(lián)合仿真工作流搭建MATLAB/Simulink COMSOL5.1 為什么要做聯(lián)合仿真COMSOL很少是產(chǎn)品開(kāi)發(fā)中唯一使用的軟件工具。在真實(shí)的開(kāi)關(guān)設(shè)備、弧焊電源或等離子體系統(tǒng)研發(fā)中外部往往有完整的控制回路比如電弧電流由電力電子拓?fù)涮峁┤蓟∵^(guò)程與機(jī)械操動(dòng)機(jī)構(gòu)聯(lián)動(dòng)。把COMSOL的多物理場(chǎng)模型嵌入Simulink的系統(tǒng)級(jí)仿真環(huán)境里才能全面考察電弧暫態(tài)行為與外部電路、控制算法之間的相互影響。聯(lián)合仿真帶來(lái)的直接好處有三個(gè)其一通過(guò)MATLAB腳本批量修改COMSOL模型參數(shù)并運(yùn)行求解這對(duì)數(shù)十上百組工況的掃描極其有效其二COMSOL模型可以作為“高保真被控對(duì)象”嵌入Simulink用于控制器設(shè)計(jì)中的硬件在環(huán)或軟件在環(huán)測(cè)試其三后期可以對(duì)COMSOL仿真結(jié)果做降階處理提煉出等效電路模型或響應(yīng)面模型再用這個(gè)降階模型嵌入系統(tǒng)級(jí)仿真大幅降低大系統(tǒng)仿真成本。5.2 基于Livelink for MATLAB的聯(lián)合仿真設(shè)置COMSOL官方提供了與MATLAB的Livelink接口配合“COMSOL for MATLAB”的使用方式非常直接。常見(jiàn)的調(diào)用套路是% 載入已建好的電弧模型 model mphload(arc_model.mph); % 修改模型參數(shù)比如外部施加電壓 model.param.set(U0, 50000); % 修改物理場(chǎng)設(shè)置比如電極表面導(dǎo)入熱流 model.physics(ht).feature(heatflux1).set(q0, 500e6); % 運(yùn)行瞬態(tài)求解 model.study(std1).run(); % 提取特定邊界或截線上的數(shù)據(jù) [t, data] mphinterp(model, T, coord, [x_array; y_array; z_array], dataset, dset1);這段腳本最大的威力在于“參數(shù)驅(qū)動(dòng)”。你可以很自然地寫(xiě)一個(gè)for循環(huán)對(duì)一組不同的電極距離或者氣壓值逐一運(yùn)行模型最后把所有結(jié)果匯總到一個(gè)數(shù)組里做優(yōu)化分析。我在做燒蝕深度與電流幅值關(guān)系研究時(shí)就是用這種循環(huán)方式一次性跑了二十多個(gè)工況節(jié)省了數(shù)周人工操作時(shí)間。在Simulink側(cè)的集成可以使用“COMSOL with MATLAB”接口或者將COMSOL模型編譯為共享庫(kù)/應(yīng)用程序在仿真步進(jìn)中調(diào)用。需要注意的是不要讓Simulink在每一個(gè)微小步長(zhǎng)里都去調(diào)用COMSOL求解這種跨度懸殊的模型交換會(huì)造成嚴(yán)重的性能瓶頸。正確做法是先在COMSOL端把典型的暫態(tài)電弧求解結(jié)果生成數(shù)據(jù)表或者輸入輸出映射再用這個(gè)映射搭建Simulink子系統(tǒng)。5.3 COMSOL 6.4新特性與Linux環(huán)境下的批處理COMSOL 6.4版本在等離子體模塊和多物理場(chǎng)耦合體驗(yàn)上有不小提升。比如對(duì)移動(dòng)網(wǎng)格自動(dòng)重剖分的控制更加精細(xì)在求解器方面新的多核并行和內(nèi)存管理對(duì)大規(guī)模模型的支撐明顯增強(qiáng)參數(shù)估計(jì)功能也更好用可以直接把仿真結(jié)果對(duì)實(shí)測(cè)曲線做校準(zhǔn)來(lái)確定熱流分配系數(shù)這類(lèi)未知參數(shù)。如果項(xiàng)目服務(wù)器跑在Linux上那就要善用命令行批處理模式。典型的命令是comsol batch -inputfile arc_model.mph -outputfile arc_result.mph -study std1 -param U050000批量提交多個(gè)仿真作業(yè)到服務(wù)器上并行運(yùn)行能大幅壓縮參數(shù)掃描周期。需要注意License類(lèi)型是網(wǎng)絡(luò)浮動(dòng)許可證多節(jié)點(diǎn)并行要合理設(shè)置MPI進(jìn)程數(shù)過(guò)高的進(jìn)程數(shù)在小模型上反而會(huì)因通信開(kāi)銷(xiāo)變得效率更低這個(gè)在集群上務(wù)必實(shí)測(cè)對(duì)比。Linux環(huán)境還有一個(gè)經(jīng)驗(yàn)之談COMSOL GUI占用的窗口性能可能不如Windows流暢但求解器性能相當(dāng)。在服務(wù)器上建議全部用batch模式跑模型本地用COMSOL桌面版做后處理和分析這樣效率最好。5.4 聯(lián)合仿真的典型落地場(chǎng)景我目前做過(guò)的聯(lián)合仿真落地方案包括斷路器開(kāi)斷過(guò)程Simulink負(fù)責(zé)電弧電流電路方程和操動(dòng)機(jī)構(gòu)運(yùn)動(dòng)COMSOL負(fù)責(zé)電弧流場(chǎng)-電磁場(chǎng)耦合瞬態(tài)仿真兩個(gè)模型之間通過(guò)每1ms交換一次“電弧電壓-電流-弧長(zhǎng)”參數(shù)電弧增材制造COMSOL模擬焊縫熔池和電弧熱流MATLAB負(fù)責(zé)路徑規(guī)劃和工藝參數(shù)優(yōu)化通過(guò)優(yōu)化算法自動(dòng)迭代仿真參數(shù)找到最佳焊接速度和電流匹配。標(biāo)準(zhǔn)化參數(shù)標(biāo)定先跑COMSOL得到不同電流、氣壓組合下的電弧溫度場(chǎng)和燒蝕深度再用MATLAB擬合代理模型生成可以用作產(chǎn)品手冊(cè)的性能曲線。每一步都談不上神秘但把這些模塊串起來(lái)的工程意義很大COMSOL不再是孤島而是真正融入了產(chǎn)品研發(fā)的數(shù)字化流程。6. 常見(jiàn)問(wèn)題與排查技巧實(shí)錄6.1 收斂失敗原因速查表電弧模型的收斂問(wèn)題幾乎每個(gè)新項(xiàng)目都會(huì)遇到。為了減少排查時(shí)間我把常見(jiàn)的導(dǎo)致發(fā)散或長(zhǎng)時(shí)間不收斂的原因整理成一個(gè)速查表在實(shí)際協(xié)作中這條路幫助很大。癥狀主要原因解決思路計(jì)算直接發(fā)散電導(dǎo)率跨數(shù)量級(jí)導(dǎo)致矩陣奇異給電導(dǎo)率設(shè)置下限截?cái)嘀当热?e-3 S/m電弧無(wú)法建立初始種子區(qū)缺失或溫度過(guò)低在電極間隙設(shè)置高溫高導(dǎo)電初始種子區(qū)溫度場(chǎng)出現(xiàn)負(fù)值或超過(guò)物理上限材料屬性外推超界給材料插值函數(shù)增加上下邊界保護(hù)截?cái)鄽埐钋€振蕩不收斂時(shí)間步長(zhǎng)過(guò)大或CFL條件破壞降低時(shí)間步長(zhǎng)上限檢查網(wǎng)格尺寸移動(dòng)網(wǎng)格負(fù)體積位移增量過(guò)大限制單步位移開(kāi)啟自動(dòng)重剖分流體壓力振蕩出口邊界設(shè)置不當(dāng)改用壓力遠(yuǎn)場(chǎng)或吸收邊界條件電磁和流場(chǎng)交替發(fā)散分離求解器迭代次數(shù)不足增加子步驟內(nèi)部迭代次數(shù)放寬容差瞬態(tài)結(jié)果出現(xiàn)尖刺BDF階數(shù)過(guò)高導(dǎo)致過(guò)沖將BDF階數(shù)調(diào)整為2排查時(shí)不要一次改多個(gè)參數(shù)一定要“單因子調(diào)試”。如果發(fā)散問(wèn)題仍不能定位可以把物理場(chǎng)一個(gè)個(gè)地關(guān)閉測(cè)試只保留電場(chǎng)-熱場(chǎng)確認(rèn)結(jié)果合理后再加入流體場(chǎng)最后再加入磁場(chǎng)。這種二分的定位方法在復(fù)雜的電弧模型里是最高效的。6.2 網(wǎng)格設(shè)計(jì)與流場(chǎng)振蕩的細(xì)節(jié)網(wǎng)格質(zhì)量往往是隱藏的最大瓶頸。電弧核心區(qū)域的網(wǎng)格必須足夠細(xì)但整個(gè)計(jì)算域的網(wǎng)格又不能浪費(fèi)得太多。實(shí)際操作中我會(huì)在網(wǎng)格序列里添加一個(gè)“尺寸”控件把電弧核心通道和電極表面附近單獨(dú)細(xì)化其他區(qū)域使用較稀疏的網(wǎng)格。典型策略是電弧直徑方向布置至少20個(gè)單元網(wǎng)格增長(zhǎng)率控制在1.2以?xún)?nèi)電極表面邊界層至少布置5~8層以保證近壁面熱流和黏性效應(yīng)解析足夠。對(duì)于流場(chǎng)振蕩還有一個(gè)不太常被提及的原因電弧模型中的氣體密度在高溫下變化很大如果采用不可壓縮流或者馬赫數(shù)較低的假設(shè)會(huì)在強(qiáng)熱源區(qū)產(chǎn)生虛假壓力波動(dòng)。正確的做法是使用可壓縮流動(dòng)選項(xiàng)或者至少開(kāi)啟“弱可壓縮”計(jì)算讓密度隨溫度瞬態(tài)變化。這個(gè)細(xì)節(jié)直接決定了射流速度和電弧收縮形態(tài)的準(zhǔn)確性千萬(wàn)不能忽略。6.3 參數(shù)掃描與批量計(jì)算的高效管理當(dāng)模型需要大批量參數(shù)掃描時(shí)單靠GUI逐個(gè)修改參數(shù)是不可行的。推薦做法是把參數(shù)文件獨(dú)立出來(lái)在COMSOL中設(shè)定好“全局參數(shù)”然后用MATLAB腳本循環(huán)修改并批量求解。我在項(xiàng)目中通常按如下思路組織把所有可調(diào)參數(shù)U0、氣壓、電極間距、材料屬性等集中到COMSOL的Global Definitions參數(shù)表中用MATLAB腳本批量修改通過(guò)mphsave保存為不同名稱(chēng)的模型文件避免覆蓋每個(gè)模型求解后自動(dòng)提取關(guān)鍵結(jié)果例如最大燒蝕深度、最大溫度、輸出電流把結(jié)果匯總到Excel或MATLAB表格進(jìn)行可視化后處理。這套流程保證了幾百組參數(shù)掃描也能在幾天內(nèi)完成而不是靠人工盯屏幕盯到崩潰。再往后結(jié)合降階模型算法你甚至可以把這些仿真數(shù)據(jù)訓(xùn)練成響應(yīng)面模型實(shí)現(xiàn)毫秒級(jí)快速預(yù)測(cè)那就完全是另一種應(yīng)用層次了。6.4 經(jīng)驗(yàn)沉淀如何搭建可復(fù)用的電弧仿真模板經(jīng)過(guò)多個(gè)項(xiàng)目的沉淀與打磨我的建議是把你調(diào)通的第一套電弧模型立刻升級(jí)為模板而不是每次從零開(kāi)始。這個(gè)模板應(yīng)該包括下面幾部分固定單位的定義電極尺寸、氣體區(qū)域尺寸、網(wǎng)格尺寸、物性參數(shù)全部使用統(tǒng)一單位制標(biāo)準(zhǔn)物理場(chǎng)設(shè)置電磁場(chǎng)、流體場(chǎng)、熱場(chǎng)的全部邊界條件都設(shè)置成可復(fù)用的命名供后續(xù)改參數(shù)標(biāo)準(zhǔn)結(jié)果處理把燒蝕深度、總電流、電弧電壓等關(guān)鍵結(jié)果用“派生值”和“報(bào)告”方式固化一鍵導(dǎo)出調(diào)試文檔記錄每次收斂失敗的成因和解決辦法形成團(tuán)隊(duì)內(nèi)部的知識(shí)庫(kù)。這樣做的價(jià)值在項(xiàng)目周期緊張時(shí)尤為明顯。你視覺(jué)上看到的可能是“我只是復(fù)制了上次的模型”實(shí)際上模板的價(jià)值是把完整的經(jīng)驗(yàn)知識(shí)沉淀下來(lái)后續(xù)每個(gè)新項(xiàng)目都不需要去交一次“智商稅”。我個(gè)人在實(shí)際操作中的體會(huì)是電弧仿真的成敗不在某一個(gè)孤立的技術(shù)點(diǎn)上而是在從物理建模到求解策略再到結(jié)果認(rèn)知的整個(gè)鏈條上。COMSOL這類(lèi)多物理場(chǎng)軟件給了我們一個(gè)優(yōu)秀的平臺(tái)但真正決定仿真深度的是建模者對(duì)物理機(jī)制的理解和對(duì)數(shù)值方法的駕馭。燒蝕深度也好瞬態(tài)電弧也好它們從來(lái)不是“點(diǎn)擊求解”就能水到渠成的結(jié)果而是需要在反復(fù)調(diào)試、反復(fù)對(duì)比中逐漸逼近真實(shí)狀態(tài)的過(guò)程。如果這篇文章能幫你少踩幾個(gè)坑自然好得很。更我希望的是你在真正跑通自己的第一個(gè)電弧模型之后能體會(huì)到那種“物理規(guī)律在屏幕上呈現(xiàn)”的興奮感然后帶著這個(gè)感覺(jué)去處理更復(fù)雜的多物理場(chǎng)問(wèn)題。仿真的下一步也許就是實(shí)驗(yàn)把模型和真實(shí)工況的差距再拉近一點(diǎn)——這比糾結(jié)于軟件按鈕要有意思得多。