MgB2高溫超導(dǎo)發(fā)電機(jī)設(shè)計(jì)及轉(zhuǎn)子結(jié)構(gòu)優(yōu)化

洪天琪，余海濤，胡敏強(qiáng)

(東南大學(xué)電氣工程學(xué)院，江蘇南京210096)

1 引言

隨著人們對(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ī)。

2 50kW高溫超導(dǎo)電機(jī)設(shè)計(jì)

2.1 高溫超導(dǎo)電機(jī)設(shè)計(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

2.2 電機(jī)結(jié)構(gòu)設(shè)計(jì)

在超導(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)題。

3 電機(jī)結(jié)構(gòu)優(yōu)化設(shè)計(jì)方法

影響磁場(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ì)目的。

4 結(jié)論

通過(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.