洪天琪,余海濤,胡敏強(qiáng)
(東南大學(xué)電氣工程學(xué)院,江蘇南京210096)
隨著人們對(duì)單電機(jī)功率和性價(jià)比等要求的提高,傳統(tǒng)的電機(jī)設(shè)計(jì)理念受到了嚴(yán)重的挑戰(zhàn)。為了克服單機(jī)功率增大導(dǎo)致體積和重量成幾何倍數(shù),以及運(yùn)輸、裝配等成本的增加[1-3],國(guó)內(nèi)外學(xué)者和技術(shù)人員已將目光投向維護(hù)成本較低、損耗較少的高溫超導(dǎo)電機(jī),而高溫超導(dǎo)材料的日益成熟,使在電機(jī)中采用高溫超導(dǎo)材料成為可能。
目前,在高溫超導(dǎo)電機(jī)設(shè)計(jì)中,國(guó)內(nèi)外學(xué)者和技術(shù)人員已在高溫超導(dǎo)材料的科學(xué)選擇和電機(jī)內(nèi)部結(jié)構(gòu)優(yōu)化設(shè)計(jì)等方面進(jìn)行了大量研究,取得了許多有價(jià)值的成果[4,7]。然而,由于現(xiàn)場(chǎng)采用的 YBCO與鉍系高溫超導(dǎo)材料,在性能特征、加工工藝和成本等方面的局限性,本文使用加工性能優(yōu)良MgB2材料作為高溫超導(dǎo)電機(jī)的勵(lì)磁繞組,克服了最小彎曲半徑等諸多限制的不足。
在電機(jī)結(jié)構(gòu)設(shè)計(jì)方面,借助商用電磁場(chǎng)軟件建立了電機(jī)結(jié)構(gòu)設(shè)計(jì)模型,并進(jìn)行了有限元仿真計(jì)算;通過(guò)引入灰色系統(tǒng)理論,對(duì)仿真結(jié)果進(jìn)行優(yōu)化建模,解決了勵(lì)磁繞組形狀和位置對(duì)繞組端部氣隙磁場(chǎng)畸變影響的優(yōu)化問(wèn)題?;谏鲜鲅芯砍晒晒υO(shè)計(jì)出功率為50kW高溫超導(dǎo)電機(jī)。
雖然超導(dǎo)電機(jī)與傳統(tǒng)電機(jī)在設(shè)計(jì)流程上有很多相似之處,但是由于超導(dǎo)電機(jī)需要在超導(dǎo)部分裝配保持低溫的制冷設(shè)備。同時(shí),無(wú)鐵芯結(jié)構(gòu)將導(dǎo)致傳統(tǒng)電機(jī)設(shè)計(jì)賴以支撐的磁路法失效,致使在設(shè)計(jì)方式上產(chǎn)生許多不同點(diǎn)[5]。
基于制冷保溫層要求和數(shù)值模擬計(jì)算,研制的50kW新型高溫超導(dǎo)電機(jī)氣隙寬度減小至2.5mm[1],并采用整體制冷的方案,以減少定、轉(zhuǎn)子間的輻射傳熱,降低轉(zhuǎn)子部分的制冷壓力。其中,定子利用液氮維持在77K的低溫,轉(zhuǎn)子利用制冷機(jī)維持在35K的超導(dǎo)工作溫度。定子尺寸選用現(xiàn)行的標(biāo)準(zhǔn)尺寸,材料為硅鋼片,因?yàn)樵诘蜏叵鹿桎撈愯F磁材料有助于減小其中渦流損耗[6]。轉(zhuǎn)子主要由傳熱部分與保溫部分兩部分構(gòu)成,起支撐作用的傳熱部分,采用輕型鋁鎂合金材料;保溫部分采用G10材料。新型高溫超導(dǎo)發(fā)電機(jī)的合理制冷、保溫設(shè)計(jì),使轉(zhuǎn)子重量大幅度減少,同時(shí),增加了電機(jī)鏤空部分,保證了制冷介質(zhì)的冷量傳遞[7]。
由于本次研制的功率為50kW的新型高溫超導(dǎo)發(fā)電機(jī),為了使定子制冷和絕緣成本之和降低,電機(jī)額定電壓選擇為1.05kV。另外,考慮到轉(zhuǎn)速過(guò)高會(huì)增加轉(zhuǎn)子及轉(zhuǎn)子制冷設(shè)備成本,故采用6極結(jié)構(gòu)。50kW新型高溫超導(dǎo)電機(jī)的主要電氣參數(shù)和結(jié)構(gòu)參數(shù)見(jiàn)表1。
表1 高溫超導(dǎo)電機(jī)主要電氣與結(jié)構(gòu)參數(shù)Tab.1 Main parameters of HTS generator
在超導(dǎo)電機(jī)結(jié)構(gòu)設(shè)計(jì)中,轉(zhuǎn)子結(jié)構(gòu)的設(shè)計(jì)是整個(gè)設(shè)計(jì)的重要環(huán)節(jié),轉(zhuǎn)子結(jié)構(gòu)的設(shè)計(jì)實(shí)質(zhì)就是超導(dǎo)勵(lì)磁繞組的設(shè)計(jì)。根據(jù)新型高溫超導(dǎo)發(fā)電機(jī)的功率要求,超導(dǎo)勵(lì)磁繞組采用的是一次成型長(zhǎng)度為600m的超導(dǎo)線材,并設(shè)定有效長(zhǎng)為400mm以節(jié)約超導(dǎo)線材,彎曲部分長(zhǎng)度為200mm,則超導(dǎo)勵(lì)磁繞組一極最大繞線匝數(shù)為600匝。超導(dǎo)線材的橫截面積為2mm2,因此,高溫超導(dǎo)繞組單極截面積為1200mm2。高溫超導(dǎo)勵(lì)磁繞組設(shè)計(jì)圖如圖1所示。
圖1 超導(dǎo)勵(lì)磁繞組設(shè)計(jì)圖Fig.1 Design graph of HTS excitation coils
由于高溫超導(dǎo)發(fā)電機(jī)的轉(zhuǎn)子未使用鐵磁材料,傳統(tǒng)磁路分析法無(wú)法在無(wú)鐵磁材料部分使用,因此,筆者們借助商用軟件建立了50kW新型高溫超導(dǎo)電機(jī)的結(jié)構(gòu)設(shè)計(jì)模型,并進(jìn)行了有限元仿真計(jì)算。電機(jī)結(jié)構(gòu)模型如圖2所示。
有限元仿真計(jì)算結(jié)果如圖3所示,從仿真計(jì)算結(jié)果上看,氣隙磁場(chǎng)磁感應(yīng)強(qiáng)度最大值為0.82T(忽略齒槽影響),且氣隙磁感應(yīng)強(qiáng)度曲線周期為256mm。由于轉(zhuǎn)子未使用鐵磁材料,磁感應(yīng)強(qiáng)度曲線在零值附近產(chǎn)生了畸變,當(dāng)x=144.2時(shí),根據(jù)標(biāo)準(zhǔn)磁感應(yīng)強(qiáng)度曲線公式得:
圖2 電機(jī)結(jié)構(gòu)模型圖Fig.2 Graph of generator's structure
式中,Bair為理想條件下氣隙磁感應(yīng)強(qiáng)度;Bmax為氣隙磁感應(yīng)強(qiáng)度最大值;T為周期;x為所求磁感應(yīng)強(qiáng)度在氣隙圓周上的對(duì)應(yīng)距離;x0為偏移距離??傻玫交冏畲笾?
圖3 電機(jī)初始設(shè)計(jì)氣隙磁感應(yīng)強(qiáng)度曲線圖Fig.3 Air-gap magnetic flux density graph in original design of generator
顯然,這種磁場(chǎng)畸變無(wú)疑會(huì)導(dǎo)致電機(jī)電壓波形產(chǎn)生畸變,是電機(jī)設(shè)計(jì)中亟待解決的問(wèn)題。
影響磁場(chǎng)畸變的因素很多,有制造工藝、超導(dǎo)材料、工作狀態(tài)等。而從電機(jī)內(nèi)電磁場(chǎng)角度考慮,這種磁場(chǎng)畸變主要影響因素是轉(zhuǎn)子勵(lì)磁繞組在電機(jī)內(nèi)的位置,即:勵(lì)磁繞組的縱向高度和彎曲半徑。實(shí)質(zhì)上優(yōu)化結(jié)構(gòu)設(shè)計(jì)減少磁場(chǎng)畸變是一個(gè)多參數(shù)選擇問(wèn)題,筆者們通過(guò)引入灰色系統(tǒng)理論,借助有限元仿真手段,給出一種電機(jī)結(jié)構(gòu)優(yōu)化設(shè)計(jì)方法,并對(duì)50kW新型高溫超導(dǎo)電機(jī)設(shè)計(jì)中,減少勵(lì)磁繞組形狀和位置對(duì)繞組端部氣隙磁場(chǎng)畸變影響的優(yōu)化設(shè)計(jì)參數(shù)進(jìn)行了具體計(jì)算。
針對(duì)減少磁場(chǎng)畸變目的建立事件集A={A1,…,Am},根據(jù)電機(jī)結(jié)構(gòu)設(shè)計(jì)中諸多影響因素建立對(duì)策集B={B1,…,Br};利用有限元仿真計(jì)算結(jié)果建立措施決策的目標(biāo)集C={C1,…,Cn},即評(píng)價(jià)方案優(yōu)劣的指標(biāo)。針對(duì)本文討論問(wèn)題,事件集A={勵(lì)磁繞組高度,彎曲半徑}。勵(lì)磁繞組高度為圖1中轉(zhuǎn)子繞組優(yōu)化移動(dòng)點(diǎn)與電機(jī)中心軸的距離。同時(shí)繞組彎曲半徑為勵(lì)磁繞組最小彎曲半徑,即內(nèi)彎曲半徑。
由于轉(zhuǎn)子制冷裝置限制,轉(zhuǎn)子勵(lì)磁繞組高度不能低于80mm,故轉(zhuǎn)子勵(lì)磁繞組高度允許范圍為80~89mm。彎曲半徑范圍可由勵(lì)磁繞組高度通過(guò)計(jì)算獲取。勵(lì)磁繞組的彎曲半徑變化范圍為3~16mm。由于在勵(lì)磁繞組高度變化過(guò)程中,其彎曲半徑不變。并認(rèn)為在勵(lì)磁繞組高度變化過(guò)程中,其彎曲半徑對(duì)氣隙磁場(chǎng)變化趨勢(shì)無(wú)影響,故在勵(lì)磁繞組高度從80~89mm變化時(shí),彎曲半徑為9mm。同理設(shè)定在彎曲半徑變化時(shí),勵(lì)磁繞組高度為85mm。考慮到無(wú)論是勵(lì)磁繞組高度,還是彎曲半徑的取值區(qū)間都是連續(xù)區(qū)間。為了便于分析,采用二分法將各對(duì)策區(qū)間離散成三種對(duì)策。對(duì)應(yīng)關(guān)系如表2所示。借助商用軟件對(duì)所有對(duì)策情況進(jìn)行有限元仿真計(jì)算,求得該對(duì)策下電機(jī)內(nèi)氣隙磁場(chǎng)參數(shù)。
表2 事件集與對(duì)策對(duì)應(yīng)關(guān)系表Tab.2 Relationship between events and solutions
勵(lì)磁繞組縱向移動(dòng)氣隙磁場(chǎng)變化情況如圖4所示。顯然,當(dāng)勵(lì)磁繞組不斷遠(yuǎn)離定子時(shí),氣隙磁場(chǎng)畸變問(wèn)題將被有效緩解,但氣隙磁感應(yīng)強(qiáng)度最大值也在不斷減小。根據(jù)圖3中氣隙磁場(chǎng)最大值和氣隙畸變最大值的定義,由圖4可知:當(dāng)勵(lì)磁繞組縱向高度90mm時(shí),磁感應(yīng)強(qiáng)度最大值為0.85T,畸變處最大畸變值為0.35T。當(dāng)勵(lì)磁繞組縱向高度85mm時(shí),磁感應(yīng)強(qiáng)度最大值為0.74T,畸變處最大畸變值為0.28T。當(dāng)繞組高度變?yōu)?0mm,磁感應(yīng)強(qiáng)度最大值為0.65T,畸變處最大畸變值為0.21T。
圖4 勵(lì)磁繞組高度變化氣隙磁場(chǎng)曲線圖Fig.4 Air-gap magnetic flux density graph under height changing
通過(guò)圖5中有限元仿真的計(jì)算結(jié)果可得,在彎曲半徑增加過(guò)程中,氣隙磁感應(yīng)強(qiáng)度曲線得到改善,同時(shí),氣隙磁密最大值在增大。根據(jù)圖3中氣隙磁場(chǎng)最大值和氣隙畸變最大值的定義,當(dāng)彎曲半徑為3mm時(shí),磁感應(yīng)強(qiáng)度最大值為0.72T,氣隙畸變最大值為0.3T。當(dāng)繞組彎曲半徑為9mm,磁感應(yīng)強(qiáng)度最大值為0.74T,氣隙畸變最大值為0.20T。當(dāng)繞組彎曲半徑為16mm,磁感應(yīng)強(qiáng)度最大值增加至0.76T,氣隙畸變最大值降低至0.15T。當(dāng)彎曲半徑增加到16mm時(shí),相鄰兩極勵(lì)磁繞組相碰。
圖5 勵(lì)磁繞組彎曲半徑變化氣隙磁場(chǎng)曲線圖Fig.5 Air-gap magnetic flux density graph under bend radius changing
為了考慮電機(jī)制造的經(jīng)濟(jì)性,在50kW新型高溫超導(dǎo)電機(jī)優(yōu)化設(shè)計(jì)中,引入了加工難度和成本費(fèi)用影響因素。研制過(guò)程顯示,彎曲半徑增加,實(shí)際生產(chǎn)過(guò)程中各極固定繞組的中心柱加工難度將減小。另外,若勵(lì)磁繞組幾何尺寸和超導(dǎo)線材長(zhǎng)度、斷面積確定,并忽略勵(lì)磁繞組在位置調(diào)整過(guò)程對(duì)磁場(chǎng)畸變的影響,制作成本費(fèi)用與加工難度成正比,即:當(dāng)勵(lì)磁繞組離定子鐵心越接近時(shí),加工難度越高,制作成本費(fèi)用也越高。同樣,當(dāng)勵(lì)磁繞組彎曲半徑越小,加工難度與制作成本費(fèi)用也越高。根據(jù)上述實(shí)踐經(jīng)驗(yàn),將制作成本費(fèi)用與加工難度用百分?jǐn)?shù)量化,即:100%表示最難加工、費(fèi)用最高,0%表示最容易加工、費(fèi)用最低。為了統(tǒng)一測(cè)量標(biāo)準(zhǔn),即在表3中的評(píng)價(jià)數(shù)值越小表示越符合設(shè)計(jì)期望,故將氣隙磁場(chǎng)最大值取倒數(shù)。量化后的費(fèi)用和加工難度與仿真計(jì)算結(jié)果綜合整理,可構(gòu)建參數(shù)統(tǒng)計(jì)見(jiàn)表3。
表3 參數(shù)統(tǒng)計(jì)表Tab.3 Statistics of parameters
表3中各指標(biāo)的實(shí)際措施中最小值和最大值分別為:umax(1)=100,umin(1)=20,umax(2)=100,=20,umax(3)=1/0.65,umin(3)=1/0.85,=0.35,umin(4)=0.15,則可列出矩陣:
矩陣中各參數(shù)使用一階線性規(guī)劃計(jì)算得出:
為了增加轉(zhuǎn)子勵(lì)磁繞組優(yōu)化結(jié)果的可靠性,采用特而菲方法確定各目標(biāo)的權(quán)重,即:ω1=0.2,ω2=0.15,ω3=0.3,ω4=0.35,通過(guò)計(jì)算可得到局勢(shì)決策綜合矩陣:
對(duì)D(Σ)按效果測(cè)度大小從左到右、從上到下排列成遞減序列矩陣D*,該矩陣即為評(píng)價(jià)決策的優(yōu)序化矩陣:
從電機(jī)內(nèi)電磁場(chǎng)角度考慮,主要影響因素是勵(lì)磁繞組位置指標(biāo),因此,采用行決策標(biāo)準(zhǔn)進(jìn)行判斷。在決策矩陣中,根據(jù)事件挑選最好的對(duì)策,稱為行決策,即在決策矩陣的行向量中選擇效果測(cè)度最大的決策元,表達(dá)式如下:
則最佳行決策為
式中,sij*為最佳決策的局勢(shì),表示事件Ai的最佳對(duì)策是Bj*。
因此,由式(4)和式(5)可知,50kW新型高溫超導(dǎo)電機(jī)中,轉(zhuǎn)子勵(lì)磁繞組位置的最優(yōu)局勢(shì)決策方案為:勵(lì)磁繞組內(nèi)切圓半徑為85mm,彎曲半徑為16mm。從圖4和圖5中的氣隙磁感應(yīng)強(qiáng)度曲線也可以直觀地得到該決策選擇對(duì)應(yīng)的勵(lì)磁繞組位置,氣隙畸變較小,磁感應(yīng)強(qiáng)度較大,同時(shí),加工難度和制作成本費(fèi)用最合理。最終電機(jī)氣隙磁場(chǎng)設(shè)計(jì)結(jié)果與初始結(jié)果對(duì)比圖如圖6所示。
圖6 超導(dǎo)發(fā)電機(jī)最終磁場(chǎng)優(yōu)化結(jié)果Fig.6 Final optimizing result of air-gap magnetic flux density graph
圖7 電機(jī)空載輸出電壓曲線圖Fig.7 Output voltage curve under no load
在確定勵(lì)磁繞組在電機(jī)內(nèi)的優(yōu)化位置后,筆者們對(duì)50kW新型高溫超導(dǎo)電機(jī)模型進(jìn)行了有限元瞬態(tài)仿真計(jì)算,得到空載輸出電壓曲線圖,其結(jié)果如圖7所示,從圖上可以看出:空載相電壓幅值為1.05kV,即輸出電壓為1.28kV。且輸出電壓已十分接近正弦曲線,達(dá)到了優(yōu)化設(shè)計(jì)目的。
通過(guò)本文的分析研究可知在高溫超導(dǎo)電機(jī)轉(zhuǎn)子不采用鐵磁材料,傳統(tǒng)磁路分析法無(wú)法應(yīng)用的情況下,借助灰色系統(tǒng)理論和有限元仿真手段,可以解決電機(jī)結(jié)構(gòu)的優(yōu)化設(shè)計(jì)問(wèn)題。本文給出了一種電機(jī)結(jié)構(gòu)優(yōu)化設(shè)計(jì)方法,并對(duì)50kW新型高溫超導(dǎo)電機(jī)設(shè)計(jì)中,減少勵(lì)磁繞組形狀和位置對(duì)繞組端部氣隙磁場(chǎng)畸變影響的優(yōu)化設(shè)計(jì)參數(shù)進(jìn)行了具體計(jì)算。
[1]Li Xiaohang,Zhou Yigang,Han Lei,et al.Design of a high temperature superconducting generator for wind power applications[J].IEEE Transactions on Applied Superconductivity,2011,21(3):1155-1158.
[2]Sharad K Singh,David W Scherbarth,Ernest S Ortoli,et al.Conceptual design of a high temperature superconducting generator[J].IEEE Transactions on Applied Superconductivity,1999,9(2):1237-1240.
[3]謝寶昌 (Xie Baochang).兆瓦級(jí)風(fēng)力發(fā)電機(jī)綜述(Overview of megawatt wind generator)[J].電機(jī)與控制應(yīng)用 (Electric Machines&Control Application),2007,34(2):1-4.
[4]謝寶昌,張涵 (Xie Baochang,Zhang Han).超導(dǎo)電機(jī)及其關(guān)鍵技術(shù)綜述 (Superconducting machine and its key technology)[A].2011-第十六屆中國(guó)小電機(jī)技術(shù)研討會(huì)論文集 (2011-The 16th China Small Motor Technology Symposium)[C].上海:上海交通大學(xué)風(fēng)電研究中心 (Shanghai:Wind Power Research Center of Shanghai Jiaotong University),2011.18-22.
[5]Naoki Maki,Tomoaki Takao,F(xiàn)uchino S,et al.Study of practical applications of HTS synchronous machines[J].IEEE Transactions on Applied Superconductivity,2005,15(2):2166-2169.
[6]陳敏,丘明 (Chen Min,Qiu Ming).鐵芯材料在低溫下的磁性能的研究 (Study on magnetic characteristics of the ferromagnetic materials at 77K)[J].電工電能新技術(shù) (Advanced Technology of Electrical Engineering and Energy),2003,22(1):35-38.
[7]Jean Leveque, Denis Netter, Bruno Douine, et al.Some considerations about the cooling of the rotor of a superconducting motor[J].IEEE Transactions on Applied Superconductivity,2007,17(1):44-51.