亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        電弧仿真收斂性研究

        2017-10-13 05:18:02吳國林張懷亮
        船電技術(shù) 2017年6期
        關(guān)鍵詞:主程序子程序電弧

        吳國林,張懷亮

        ?

        電弧仿真收斂性研究

        吳國林1,張懷亮2

        (1. 中國人民解放軍91003部隊(duì),海南三亞572000; 2. 海軍裝備部裝備采購中心,北京170001)

        電弧仿真的研究,對(duì)于電器開關(guān)的發(fā)展改進(jìn)有著重要意義。本文以穩(wěn)態(tài)空氣電弧為研究對(duì)象,對(duì)其中心軸附近近極區(qū)域建立了二維軸對(duì)稱模型,利用MATLAB軟件編程,采用局部細(xì)化網(wǎng)格的電弧仿真方法,對(duì)電弧近極的電位場進(jìn)行了仿真,給出了一種解決電弧仿真過程中發(fā)散問題的新思路。

        電弧仿真 收斂性 MATLAB

        0 引言

        隨著電網(wǎng)中故障短路電流水平的不斷提高[1-3],大容量分?jǐn)嘣O(shè)備的應(yīng)用也越來越廣泛。分?jǐn)嘣O(shè)備在分?jǐn)喽搪冯娏鲿r(shí),會(huì)在分?jǐn)嗟耐髦w之間產(chǎn)生電弧[4],大容量分段設(shè)備的廣泛應(yīng)用也使得電弧的研究越來越重要。電弧近極區(qū)溫度的劇變使得近極區(qū)電弧物性參數(shù)變化劇烈,對(duì)該區(qū)域進(jìn)行仿真時(shí),若采用常規(guī)的網(wǎng)格劃分方法,則極易使得計(jì)算發(fā)散。目前為避免發(fā)散問題的出現(xiàn),大都將整個(gè)模型或該區(qū)域的網(wǎng)格劃分得十分密集[5-7],使得電弧的仿真計(jì)算動(dòng)輒數(shù)十小時(shí),造成很大的資源浪費(fèi)。

        本文以穩(wěn)態(tài)空氣電弧為研究對(duì)象,對(duì)其中心軸附近近極區(qū)域建立了二維軸對(duì)稱模型,利用MATLAB軟件編程,給出了一種局部細(xì)化網(wǎng)格的方法,使得仿真計(jì)算大幅簡化。

        1 模型建立

        1.1 物理模型

        如圖1所示,給出了電弧模型示意圖,電弧為陽極與陰極之間的虛線部分,目前普遍將電弧作為圓柱體仿真。本文為簡化模型,將電弧轉(zhuǎn)化為二維模型,又由于本章的目的主要是解決電弧近極部分由于物性參數(shù)劇變帶來的發(fā)散問題,從而僅取電弧中心軸附近區(qū)域進(jìn)行仿真,即圖1中斜線區(qū)域,并且仿真區(qū)域軸向長度為電弧軸向長度1/2。如圖2所示,給出近極電弧模型的尺寸及邊界條件,由于模型選取自電弧中心軸處,可認(rèn)為電弧模型與電弧其他部分沒有熱量交換,模型熱量只能由電極傳出。

        1.2 數(shù)學(xué)模型

        1.2.1 電場方程

        式(1)中,及分別為圓柱坐標(biāo)徑向及軸向分量;表示電位,單位;表示電導(dǎo)率,單位1/s。

        圖1 電弧模型示意圖

        圖2 模型邊界

        1.2.2 導(dǎo)熱微分方程

        式(2)中,表示溫度,單位為K;表示單位體積單位時(shí)間內(nèi)產(chǎn)生的熱量,單位J;表示熱導(dǎo)率,單位W/m?K[8]。

        1.3 離散方法

        本文采用有限差分法對(duì)所建立的數(shù)學(xué)模型進(jìn)行離散。首先對(duì)仿真區(qū)域進(jìn)行網(wǎng)格劃分。用(i,j)表示點(diǎn)(i,j)的溫度,用(i,j)表示點(diǎn)(i,j)的電位,如圖3所示。

        圖3 網(wǎng)格剖分

        進(jìn)行網(wǎng)格劃分后,任一內(nèi)部節(jié)點(diǎn)(i,j)在方向上的偏導(dǎo)可以由其四周的點(diǎn)之間的關(guān)系進(jìn)行描述,從而得到內(nèi)部節(jié)點(diǎn)(i,j)的離散公式(以電位的離散公式為例):

        2 收斂條件

        由文獻(xiàn)[9-11]可知,對(duì)于式(3)所示線性方程組,若系數(shù)矩陣中的每一行中不在主對(duì)角線上的所有元素的絕對(duì)值之和小于同一行中主對(duì)角線元素之絕對(duì)值,即:

        則稱(3)為系數(shù)矩陣按行嚴(yán)格對(duì)角占優(yōu),并且對(duì)于該方程組,簡單迭代法及Seidel迭代法都收斂。

        (3)

        3 MATLAB程序設(shè)計(jì)

        3.1 MATLAB程序設(shè)計(jì)

        如圖4及圖5所示,給出了近極電弧仿真計(jì)算的MATLAB程序框圖,包括主程序及子程序兩部分,主程序?qū)崿F(xiàn)大網(wǎng)格劃分模型的仿真計(jì)算,并實(shí)現(xiàn)對(duì)各網(wǎng)格節(jié)點(diǎn)收斂性的實(shí)時(shí)監(jiān)測。一旦發(fā)現(xiàn)不滿足嚴(yán)格對(duì)角占優(yōu)的點(diǎn),即調(diào)用子程序,將該點(diǎn)坐標(biāo)輸入給子程序,子程序?qū)⒉环€(wěn)定點(diǎn)處以密集網(wǎng)格劃分,并以迭代法求解不穩(wěn)定點(diǎn)的溫度值及電位值,并將計(jì)算結(jié)果返回給主程序,主程序再進(jìn)行后續(xù)的迭代計(jì)算。

        圖4 主程序

        3.2 局部密集網(wǎng)格劃分方法

        如圖6所示,為局部密集網(wǎng)格細(xì)化方法示意圖,下面以電位的處理為例,介紹實(shí)現(xiàn)方法。當(dāng)主程序以迭代法解方程組時(shí),檢測到點(diǎn)(0,0)不滿足按行嚴(yán)格對(duì)角占優(yōu)時(shí),即調(diào)用子程序,并將1,,3,12,3輸入給子程序。若子程序中將網(wǎng)格尺寸大小變?yōu)樵瓉淼?/n,則a, b, c, d四點(diǎn)電位分別為(以n=3為例):

        圖6 局部網(wǎng)格細(xì)化

        同理,紅色虛線框上其他點(diǎn)的電位均以此方法給出,紅色虛線框上點(diǎn)的電位值實(shí)際上相當(dāng)于局部小網(wǎng)格區(qū)域的第一類邊界條件。紅色虛線框內(nèi)部點(diǎn)的電位值也以這種方法給出,這就相當(dāng)于給所有點(diǎn)賦初值。局部小網(wǎng)格區(qū)域?qū)嶋H上相當(dāng)于一個(gè)新的模型,滿足電位微分方程及導(dǎo)熱微分方程,其邊界條件均為第一類邊界條件,同樣用迭代法求解該方程組,可得點(diǎn)()的電位值及溫度值,子程序?qū)Ⅻc(diǎn)()的電位值及溫度值返回給主程序,主程序即得到了滿足按行嚴(yán)格對(duì)角占優(yōu)求解得到的收斂值。

        4 近極電弧仿真結(jié)果

        如圖7所示,給出了模型電位的軸向分布,即圖2中電弧模型軸AB上的電位分布。仿真中將電弧在軸向上劃分為101個(gè)大網(wǎng)格,網(wǎng)格0對(duì)應(yīng)電弧的中心點(diǎn)A,網(wǎng)格101對(duì)應(yīng)于點(diǎn)B。圖7中,縱軸表示電位值,橫軸表示網(wǎng)格。由電位分布可以看出,近極電弧遠(yuǎn)離觸頭表面的大部分區(qū)域的電位梯度都很小,而在靠近觸頭表面的小距離內(nèi)發(fā)生了巨大的躍升。

        圖7 電位分布

        文獻(xiàn)[12]中給出了普遍認(rèn)為的電弧電位的軸向分布,如圖8所示。電弧由三部分組成:陰極區(qū)、陽極區(qū)以及等離子區(qū)。陰極區(qū)緊鄰陰極表面,并且厚度很小,但是大部分電弧電壓都施加在這一個(gè)狹小的區(qū)域,這使得陰極區(qū)具有很高的電場強(qiáng)度。當(dāng)電流足夠大時(shí),會(huì)在靠近陽極表面的位置出現(xiàn)陽極區(qū),陽極區(qū)也分擔(dān)一部分電弧電壓。等離子區(qū)是位于陰極區(qū)和陽極區(qū)之間的廣大區(qū)域,等離子區(qū)中充滿著等離子體,等離子區(qū)中的正、負(fù)粒子都具有很高的密度,使得等離子區(qū)的導(dǎo)電性非常好,從而等離子區(qū)電位梯度較小。

        圖8 電弧電位軸向分布

        1-陰極區(qū) 2-等離子區(qū) 3-陽極區(qū)

        5 結(jié)論

        本文以MATLAB軟件自編程,提出了一種仿真模型局部采用密集網(wǎng)格劃分的方法,可以有效避免電弧仿真易發(fā)散問題,給出了一種電弧仿真的新思路。

        [1] 袁志方, 莊勁武, 王晨, 等. 石英砂對(duì)于電磁斥力高速開斷器介質(zhì)恢復(fù)特性的影響[J].高電壓技術(shù), 2014, 40(1): 294-299.

        [2] 陳博, 莊勁武, 肖翼洋, 等. 10KV/2KA混合型限流熔斷器用電弧觸發(fā)器的分析與設(shè)計(jì)[J]. 高電壓技術(shù), 2012, 38(8): 1948-1955.

        [3] 袁志方, 莊勁武, 王晨, 等. 窄縫滅弧法提升電磁斥力高速開斷器電弧電壓的分析與試驗(yàn)[J]. 中國電機(jī)工程學(xué)報(bào), 2013, 33(33): 139-144.

        [4] 顧雪晨, 王晨. 直流電力系統(tǒng)短路限流裝置研究[J]. 船舶, 2012, 23(3): 60-62.

        [5] Alferov F, Nevrovskii A, Sidorov A. Anode mode of vacuum arc in multirod electrode system[J]. High Temperature, 2002, 40(1): 15-20.

        [6] Wang L, Jia S, Yang D, et al. Modeling and simulation of anode activity in high-current vacuum Arc[J]. Phys. D: Appl. Phys., 2009, 42(5): 62-67.

        [7] 王立軍, 賈申利, 劉宇, 等. 縱磁下真空電弧陽極熱過程的仿真[J]. 中國電機(jī)工程學(xué)報(bào), 2011, 26(3): 66-72.

        [8] 戴鍋生.傳熱學(xué)(第二版)[M]. 北京: 高等教育出版社, 2011: 22-25.

        [9] 何漢林, 梅家斌. 數(shù)值分析[M]. 北京: 科學(xué)出版社, 2007: 79-82.

        [10] 何漢林, 魏汝祥, 李衛(wèi)軍. 數(shù)值分析[M]. 武漢: 湖北科學(xué)技術(shù)出版社, 1999.

        [11] 羅家洪, 方衛(wèi)東. 矩陣分析引論[M]. 廣州: 華南理工大學(xué)出版社, 2006.

        [12] 王其平. 電器電弧理論[M]. 北京: 機(jī)械工業(yè)出版社, 1992.

        Research on the Convergence of Arc Simulation

        Wu Guolin1, Zhang Huailiang2

        (1. Unit 91003 of PLA, Shanya 572000, Hainan, China; 2. Armaments Procurement Agency, Naval Equipment Ministry, Beijing 100071, China)

        TM5 TP391.9

        A

        1003-4862(2017)06-0006-03

        2017-03-13

        吳國林(1990-),男,碩士。專業(yè)方向:電氣工程。E-mail: wuguolin1990@163.com

        猜你喜歡
        主程序子程序電弧
        故障電弧探測器與故障電弧保護(hù)裝置在工程中的應(yīng)用分析
        淺談數(shù)控銑削技術(shù)代碼程序的嵌套方式研究
        2219鋁合金激光電弧復(fù)合焊接及其溫度場的模擬
        電控冰箱軟件模塊化設(shè)計(jì)
        日用電器(2019年12期)2019-02-26 01:53:06
        時(shí)光倒流 換回PotPlayer老圖標(biāo)
        航空電氣系統(tǒng)中故障電弧的分析
        電子制作(2017年22期)2017-02-02 07:10:15
        淺談子程序在數(shù)控車編程中的應(yīng)用
        子程序在數(shù)控車加工槽中的應(yīng)用探索
        西門子840D系統(tǒng)JOG模式下PLC調(diào)用并執(zhí)行NC程序
        弓網(wǎng)離線電弧電磁干擾機(jī)理及防護(hù)
        国产在线视频一区二区三区不卡| 久久久久无码国产精品不卡| 午夜短视频日韩免费| 国产亚洲青春草在线视频| 亚洲av成人一区二区| 国产精品久久久久高潮| 国产剧情av麻豆香蕉精品| 精品三级久久久久久久| 强迫人妻hd中文字幕| 国产区精品一区二区不卡中文| 国产av人人夜夜澡人人爽| 午夜福利影院不卡影院| 最新日本久久中文字幕| 噜噜综合亚洲av中文无码| 丰满少妇被猛烈进入| 亚洲国产一区二区三区在线视频 | 日韩AVAV天堂AV在线| 日韩av一区在线播放| 亚洲国产国语在线对白观看| 曰韩人妻无码一区二区三区综合部| 亚洲区在线播放| 一区二区三区手机看片日本韩国 | 丁香综合网| 国产一区二区三区经典| 99e99精选视频在线观看| 欧美精品一区二区蜜臀亚洲| 无码一区久久久久久久绯色AV| 国产精品一品二区三区| 国产成人av一区二区三区在线观看| 欧美成a人片在线观看久| 久久久99久久久国产自输拍| 青青草成人免费在线视频| 午夜福利一区二区三区在线观看| 丰满岳乱妇在线观看中字无码 | 亚洲国产一区二区av| 亚洲av永久无码一区二区三区| 国产精品一区二区久久| 91久久精品一区二区喷水喷白浆 | 青青视频在线播放免费的| 狠狠躁18三区二区一区| 免费人成视频在线观看视频|