羅柏文, 楊華科, 周知進(jìn), 王朝輝
(1.湖南科技大學(xué) 海洋礦產(chǎn)資源探采裝備技術(shù)湖南省工程實(shí)驗(yàn)室,湖南 湘潭 411201; 2.貴州理工學(xué)院 機(jī)械工程學(xué)院, 貴陽 550000)
?
基于改進(jìn)型NSGA-II算法的深海高頻變壓器優(yōu)化研究
羅柏文1, 楊華科1, 周知進(jìn)2, 王朝輝1
(1.湖南科技大學(xué) 海洋礦產(chǎn)資源探采裝備技術(shù)湖南省工程實(shí)驗(yàn)室,湖南 湘潭 411201; 2.貴州理工學(xué)院 機(jī)械工程學(xué)院, 貴陽 550000)
針對(duì)深海高頻變壓器建立其體積和損耗的多目標(biāo)優(yōu)化數(shù)學(xué)模型。利用遺傳算法在求解非線性多目標(biāo)優(yōu)化問題中的優(yōu)勢(shì),解決深海高頻變壓器損耗目標(biāo)和體積目標(biāo)相互矛盾的問題。改進(jìn)NSGA-II算法中的非支配排序策略和選擇截?cái)嗖呗裕苟咛幱谕环侵鋵?,從而增加了周圍密度小的個(gè)體的遺傳概率,提高了算法的全局搜索能力和運(yùn)行速度。運(yùn)用改進(jìn)型NSGA-II算法對(duì)6 kW高頻變壓器進(jìn)行優(yōu)化設(shè)計(jì)并根據(jù)優(yōu)化結(jié)果制造出高頻變壓器實(shí)體。試驗(yàn)測(cè)試結(jié)果表明,采用改進(jìn)型NSGA-II算法能有效地降低深海高頻變壓器的總損耗和體積。
深海;高頻變壓器;多目標(biāo)優(yōu)化;NSGA-II算法
NSGA-II algorithm
深海高頻變壓器要求體積小、重量輕、耐高壓和高效率,其優(yōu)化設(shè)計(jì)是在給定的設(shè)計(jì)要求下,選擇磁芯和繞組參數(shù),以降低變壓器磁芯損耗和繞組損耗所引起的溫升,同時(shí)減小變壓器的體積。由于磁芯材料的幾何結(jié)構(gòu)、繞組導(dǎo)線的大小、層間距、高度及初、次級(jí)的布局等設(shè)計(jì)參數(shù)與其損耗、體積等設(shè)計(jì)目標(biāo)之間相互影響,同時(shí),體積和損耗等目標(biāo)之間存在相互矛盾的關(guān)系[1-4]。因此,深海高頻變壓器優(yōu)化設(shè)計(jì)是一個(gè)多變量、多約束、多目標(biāo)的非線性優(yōu)化問題。而目前針對(duì)多目標(biāo)優(yōu)化問題常采用的智能算法有NSGA-II算法、粒子群算法、蟻群算法等[5-8]。NSGA-II算法采用非支配排序策略,選擇截?cái)嗖呗?,基于分層擁擠距離選擇算子策略,其SBX交叉算子策略使解集的全局搜索能力得到很大的提高,是解決多目標(biāo)優(yōu)化問題最廣泛的方法。由于傳統(tǒng)的高頻變壓器設(shè)計(jì)時(shí)參數(shù)選擇主觀性較大,重復(fù)設(shè)計(jì)計(jì)算工作量大,而遺傳算法在非線性多目標(biāo)優(yōu)化問題方面的具有一定的優(yōu)勢(shì)。因此,該文在對(duì)NSGA-II進(jìn)行改進(jìn)的基礎(chǔ)上,將其運(yùn)用到深海高頻變壓器的優(yōu)化設(shè)計(jì)之中。
1.1 優(yōu)化數(shù)學(xué)模型
深海高頻變壓器的優(yōu)化數(shù)學(xué)模型如下:
(1)
(2)
(3)
式中:x為決策向量,x=[x1,x2];y為目標(biāo)向量;gi(x)為不等式約束因子;gi為約束量界定常數(shù);hj(x)為等式約束因子。
1.2 目標(biāo)函數(shù)、決策變量和約束條件
考慮到深海高頻變壓器的性能、結(jié)構(gòu)及應(yīng)用場(chǎng)合的多樣性,對(duì)設(shè)計(jì)對(duì)象進(jìn)行如下設(shè)定:
(1) 高頻變壓器磁芯結(jié)構(gòu)采用圓環(huán)形磁芯,選擇內(nèi)外徑差值較小、厚度適中的磁芯,磁芯材料選用鐵基非晶合金。
(2) 針對(duì)特定的工作頻率,選擇繞組圓導(dǎo)線直徑為4個(gè)~6個(gè)集膚深度[8]。
(3) 繞組布局采用初級(jí)與次級(jí)交叉布局。初級(jí)、次級(jí)繞組均勻密繞在圓環(huán)形磁芯的外表面。
深海高頻變壓器的主要指標(biāo)分為變壓器負(fù)載損耗和體積兩個(gè)部分。在滿足設(shè)計(jì)要求的前提下使其負(fù)載損耗最小,以提高產(chǎn)品的可靠性,與此同時(shí)盡量減小變壓器的體積。因此,針對(duì)上述兩個(gè)指標(biāo),基于深海高頻變壓器設(shè)計(jì)相關(guān)公式分別推導(dǎo)出其負(fù)載損耗和體積的兩個(gè)目標(biāo)函數(shù)如下:
體積目標(biāo)函數(shù):
(4)
負(fù)載損耗目標(biāo)函數(shù):
(5)
式中:AP為磁芯幾何常數(shù),其值與變壓器體積成正比;K為體積與AP值的相關(guān)系數(shù);Fcu為變壓器銅損;Ffe為變壓器鐵損;Fo為附加損耗;k為交流電阻與直流電阻之比; Pout為變壓器最大輸出功率;Lo為初、次級(jí)繞組平均匝數(shù)長度;K0為繞組導(dǎo)線直流電阻率;f為導(dǎo)線電流密度;Kf為波形系數(shù);f為工作頻率;B為工作磁感應(yīng)強(qiáng)度;Ae為磁芯窗口利用系數(shù);t為AP值系數(shù);D為最大占空比;η為變壓器轉(zhuǎn)換效率; F為磁芯損耗波形因素,方波取π/4;K1為磁芯損耗系數(shù);a、b為指數(shù)系數(shù);V為磁芯體積; ρ為磁芯材料的密度;決策變量x為工作磁感應(yīng)強(qiáng)度B和導(dǎo)線電流密度J組成的向量空間,其值與變壓器動(dòng)態(tài)損耗及體積都成非線性相關(guān)。
約束條件為高頻變壓器效率下限:
(6)
若滿足條件,則結(jié)果保留;若不滿足條件,該結(jié)果則舍去。
圖1 NSGA-II進(jìn)化圖
2.1 NSGA-II算法分析
圖2 個(gè)體周圍密度信息
NSGA-II算法種群進(jìn)化如圖1所示。從圖1中可以看出:父代種群經(jīng)過基于非支配和擁擠度的選擇、交叉和變異算子后得到與父代種群個(gè)體數(shù)目N相等的子代種群。父代和子代種群合并后得到一個(gè)種群數(shù)目為2N的新種群,對(duì)該種群進(jìn)行非支配排序策略,得到非支配個(gè)體集Q1,Q2,Q3,Q4并計(jì)算非支配個(gè)體集的擁擠距離,然后通過選擇截?cái)嗖呗缘玫椒侵鋫€(gè)體集F1,F(xiàn)2,F(xiàn)3。
在NSGA-II算法非支配排序策略中,種群中的個(gè)體是將非支配排序中所處的層數(shù)值賦予為個(gè)體的適應(yīng)度評(píng)價(jià)函數(shù)來進(jìn)行遺傳的[9,10]。這將有利于克服在早期多目標(biāo)優(yōu)化的過程中,通過線性加權(quán)法形成適應(yīng)度評(píng)價(jià)函數(shù)時(shí),造成的客觀性較差、計(jì)算量較大等缺陷,但同時(shí)也帶來了以下兩點(diǎn)局限性:(1) 每一級(jí)非支配排序都要考慮該個(gè)體集中所有的個(gè)體,計(jì)算量較大;(2) 對(duì)處于同一非支配層級(jí)的個(gè)體,周圍密度不同的個(gè)體進(jìn)化到下一代種群的概率相同,使下一代種群個(gè)體分布可能會(huì)出現(xiàn)局部集中現(xiàn)象,降低了種群分布的多樣性,從而影響了算法的全局搜索能力。
個(gè)體周圍密度信息如圖2所示。從圖2中可以看出:個(gè)體a,b分別所受支配個(gè)體數(shù)分別為3和2,且個(gè)體a周圍個(gè)體的密度大于個(gè)體b,然而由于a,b處于同一非支配層級(jí),其適應(yīng)度評(píng)價(jià)函數(shù)相等,遺傳到下一代的概率相同,從而影響了算法的全局搜索能力。
在NSGA-II算法選擇截?cái)嗖呗灾?,將父代種群與同規(guī)模的子代種群進(jìn)行合并后,再進(jìn)行選擇截?cái)嘁赃x擇占優(yōu)個(gè)體。這將增加備選個(gè)體的數(shù)量,有利于解集的多樣性。但在合并的種群中,處于較高非支配層、擁擠度較大且重復(fù)的個(gè)體比處于低非支配層、擁擠度較小的個(gè)體擁有更高的遺傳概率,這將導(dǎo)致種群收斂分布不均。
2.2 改進(jìn)策略
針對(duì)2.1節(jié)中所述非支配排序策略的局限性,在個(gè)體適應(yīng)度評(píng)價(jià)函數(shù)計(jì)算時(shí)加入種群個(gè)體周圍密度信息,個(gè)體適應(yīng)度評(píng)價(jià)函數(shù)為個(gè)體處的非支配層級(jí)與該個(gè)體所受支配個(gè)體數(shù)之和,其定義如下:
(7)
圖3 算法流程圖
根據(jù)公式6,首先將種群中n為零的個(gè)體存入第一級(jí)非支配解集,對(duì)于同一Pareto值中的個(gè)體p以及p所支配的個(gè)體集sp,將sp-1的值賦給sp。個(gè)體q的改進(jìn)排序值為Pareto值與其所支配個(gè)體累計(jì)排序之和,若n為零,則存入下一級(jí)非支配解集,直至算法終止。
針對(duì)選擇截?cái)嗖呗缘木窒扌?,在NSGA-II改進(jìn)型算法中采取下列措施:在對(duì)所有種群進(jìn)行非支配排序之前先進(jìn)行擁擠度計(jì)算并排序,選取擁擠度較大的個(gè)體組成一個(gè)新的種群。在非支配排序時(shí)設(shè)置一個(gè)終止信號(hào),即當(dāng)非支配個(gè)體集的數(shù)目等于N時(shí),跳出循環(huán),運(yùn)算終止。在圖1所示過程中,非支配子集Q4在選擇截?cái)噙^程應(yīng)被修剪,所以設(shè)置非支配排序到Q3子集時(shí),將跳出排序過程,有利于簡化非支配排序運(yùn)算量。NSGA-II改進(jìn)型算法簡化流程如圖3所示。
對(duì)開關(guān)電源用的6 kW高頻變壓器進(jìn)行實(shí)例分析,以高頻變壓器小體積、低損耗為目標(biāo)函數(shù),運(yùn)用改進(jìn)型NSGA-II算法對(duì)其參數(shù)B、J進(jìn)行優(yōu)化設(shè)計(jì)。
設(shè)計(jì)中,高頻變壓器磁芯選擇環(huán)形鐵基非晶合金,最大占空比為0.9,輸出功率為6 kW,最大工作頻率為20 kHz,初級(jí)出入電壓240 V,匝數(shù)比為1:1.5。算法設(shè)置如下:決策變量B、J的上下限分別為0,0.5;0,500。約束條件變壓器效率下限為0.95,溫升上限為50 K,父代群體設(shè)置為200,子代粒子群大小設(shè)置為30,最大迭代次數(shù)設(shè)置為200。通過計(jì)算得到高頻變壓器改進(jìn)NSGA-II算法和基本型算法的優(yōu)化解空間分布分別如圖4、圖5所示。
圖4 改進(jìn)型NSGA-II算法優(yōu)化解空間分布
圖5 NSGA-II算法優(yōu)化解空間分布
利用傳統(tǒng)設(shè)計(jì)方法(AP法)得出優(yōu)化前的計(jì)算結(jié)果,并對(duì)比分析改進(jìn)型算法、基本型算法和傳統(tǒng)設(shè)計(jì)的計(jì)算結(jié)果見表1。
表1 6 kW高頻變壓器設(shè)計(jì)結(jié)果對(duì)比
從表1中可以看出:基于NSGA-II改進(jìn)型算法和基本型算法的優(yōu)化設(shè)計(jì)結(jié)果相比于基于傳統(tǒng)設(shè)計(jì)的結(jié)果在高頻變壓器損耗方面分別下降了33.62 W、33.41 W;高頻變壓器的效率分別提高了0.561%、0.557%;高頻變壓器的AP值較傳統(tǒng)設(shè)計(jì)分別降低了18.9%、18.7%,對(duì)應(yīng)的高頻變壓器的體積也相應(yīng)地降低;改進(jìn)型NSGA-II算法的優(yōu)化速度和優(yōu)化結(jié)果明顯優(yōu)于基本型NSGA-II算法。
為了驗(yàn)證NSGA-II改進(jìn)型算法在高頻變壓器優(yōu)化設(shè)計(jì)方面的正確性,根據(jù)優(yōu)化設(shè)計(jì)后的參數(shù)設(shè)計(jì)出優(yōu)化后的6 kW高頻變壓器實(shí)物并對(duì)其進(jìn)行動(dòng)態(tài)損耗測(cè)試。由于現(xiàn)階段對(duì)大功率高頻變壓器的動(dòng)態(tài)損耗無法實(shí)現(xiàn)直接測(cè)量,因此本次測(cè)試采用間接測(cè)量的方法,測(cè)量方法和測(cè)試現(xiàn)場(chǎng)分別如圖6、圖7所示。
圖6 大功率移相全橋開關(guān)電源動(dòng)態(tài)損耗測(cè)試原理圖
圖7 測(cè)試現(xiàn)場(chǎng)
根據(jù)記錄的輸入輸出電壓、電流數(shù)據(jù),整理出全橋移相拓?fù)潆娫纯傮w損耗的測(cè)量數(shù)據(jù),見表2。
表2 移相全橋開關(guān)電源測(cè)試總體損耗數(shù)據(jù)
由于開關(guān)電源中磁性元件分為高頻變壓器和電感兩個(gè)部分,磁性元件動(dòng)態(tài)損耗值占電路總損耗的5%。根據(jù)表2中測(cè)得的總損耗數(shù)據(jù)可以得出6 kW高頻變壓器和電感的總損耗在不同輸出功率時(shí)最大損耗估算值。如忽略電感損耗值,則6 kW高頻變壓器在不同功率輸出時(shí)的最大損耗如圖8所示,6 kW高頻變壓器實(shí)物如圖9所示。
圖8 6 kW高頻變壓器在不同功率輸出時(shí)的最大損耗
圖9 6 kW高頻變壓器實(shí)物
從圖8可以看出:測(cè)試電路輸出功率增加時(shí),6 kW高頻變壓器損耗呈上升趨勢(shì)。當(dāng)電路輸出功率為6 kW、其他測(cè)試條件滿足高頻變壓器的設(shè)計(jì)要求時(shí),此時(shí)6 kW變壓器的測(cè)試最大損耗值為20.655 W,測(cè)試最低效率為99.639%,為最終測(cè)試結(jié)果。將實(shí)驗(yàn)測(cè)得的損耗結(jié)果與傳統(tǒng)設(shè)計(jì)以及兩種優(yōu)化設(shè)計(jì)結(jié)果進(jìn)行對(duì)比,見表3。
表3 6 kW高頻變壓器損耗數(shù)據(jù)對(duì)比
從表3可以看出:實(shí)驗(yàn)測(cè)量結(jié)果比兩種優(yōu)化算法結(jié)果差,但比傳統(tǒng)設(shè)計(jì)結(jié)果好,這是因?yàn)樵摐y(cè)量值為高頻變壓器和電感的損耗值之和,高頻變壓器損耗實(shí)際值應(yīng)低于該測(cè)量值;其次,變壓器的繞組之間存在兩層亞胺膜,層間距大于理論設(shè)計(jì)值,變壓器繞制過程中也會(huì)增加一定的附加損耗,導(dǎo)致最終測(cè)試結(jié)果低于理論設(shè)計(jì)值。因此,實(shí)驗(yàn)測(cè)量結(jié)果在傳統(tǒng)設(shè)計(jì)值與優(yōu)化后計(jì)算值之間符合理論預(yù)期,實(shí)驗(yàn)測(cè)試結(jié)果進(jìn)一步驗(yàn)證了優(yōu)化設(shè)計(jì)方法和結(jié)果的正確性。
(1) 在分析NSGA-II算法的基礎(chǔ)上,提出在非支配排序策略前先進(jìn)行擁擠度排序。選擇擁擠度較大的個(gè)體,在非支配排序時(shí)考慮了周圍密度信息和設(shè)置終止信號(hào),提高了算法的運(yùn)行速度。
(2) 將NSGA-II改進(jìn)型算法運(yùn)用到高頻變壓器的優(yōu)化設(shè)計(jì)中,實(shí)驗(yàn)結(jié)果表明改進(jìn)型算法的優(yōu)化結(jié)果好于基本型算法,相對(duì)于傳統(tǒng)設(shè)計(jì)方法在降低損耗和減小體積方面有了明顯提高。
因此,基于NSGA-II改進(jìn)型算法能在深海高頻變壓器優(yōu)化設(shè)計(jì)中可針對(duì)特定的工作要求選擇合理的設(shè)計(jì)參數(shù),從而降低總損耗和體積,提高工作效率。
[1] Ouyang Z W,Thomson O C, Michael A E. Optimal Design and Tradeoff Analysis of Planar Transformer in High-Power DC-DC Converters[J]. Industrial Electronics,2012,59(7):2800-2810.
[2] MU M K,SU Y P, LI Q,et al. Magnetic Characterization of Low temperature Co-fired Ceramic (LTCC) Ferrite Materials for High Frequency Power Converters[C]. Phoenix, AZ: Energy Conversion Congress and Exposition,2011.
[3] Agheb E, Bahmani A, Hans K H,etal. Core Loss Behavior in High Frequency High Power Transformers—II: Arbitrary excitation[J]. Renewable Sustainable Energy,2012,4(3):31-44.
[4] 姚偉,張周勝,肖登明.開關(guān)頻率及線圈層數(shù)、高度對(duì)變壓器銅損的影響[J].上海電機(jī)學(xué)院學(xué)報(bào),2009,12(1):5-9.
[5] 史峰,王輝,胡斐,等.MATLAB智能算法30個(gè)案例分析[M].北京:北京航空航天大學(xué)出版社,2011.
[6] 付立,竇明罡,朱建凱.非支配排序遺傳算法的改進(jìn)[J]. 計(jì)算機(jī)與數(shù)字工程,2011,39(02):11-14.
[7] 曲敏,高岳林,江巧永.基于Pareto鄰域交叉算子的多目標(biāo)粒子群優(yōu)化算法[J].計(jì)算機(jī)應(yīng)用,2011,31 (7):1789-1792.
[8] 王朝輝,羅柏文,周知進(jìn).高頻變壓器繞組渦流損耗等效因子分析[J].電子器件,2014,37(4):754-756.
[9] 陳婕,熊盛武,林婉如.NSGA-II算法的改進(jìn)策略研究[J].計(jì)算機(jī)工程與應(yīng)用,2011,49(19):42-45.
[10] 凌海風(fēng),周獻(xiàn)中,江勛林,等.改進(jìn)的約束多目標(biāo)粒子群算法[J]. 計(jì)算機(jī)應(yīng)用,2012,32(5):1320-1324.
The Deep-sea High-frequency Transformer Optimization Research Based on the Improved NSGA- II Algorithm
LUO Bo-wen1, YANG Hua-ke1, ZHOU Zhi-jin2, WANG Zhao-hui1
(1. Marine Mineral Resources Exploration and Mining Equipment and Technology of Hunan Province Engineering Laboratory, Hunan University of Science and Technology, Hunan Xiangtan 411201, China; 2. Guizhou Institute of Technology College of Mechanical Engineering, Guiyang 550000, China)
Aimed at the deep-sea high-frequency transformer, the multi-objective optimization mathematical model of the volume and the depletion can be established. It deals with the optimization design of contradiction between loss and volume target of the deep-sea high-frequency transformer, taking the advantage of genetic algorithm into solving nonlinear multi-objective optimization problems. By changing and improving the non-dominated sorting and choosing truncation strategies in the NSGA-Ⅱalgorithm, individuals with low density in one non-dominated level can possess larger genetic probability, which improves the global search ability of the algorithm and increases its speed. The improved NSGA-Ⅱalgorithm was used to make optimization designs to 6 kW high-frequency transformer, and the high-frequency transformer entitycould be manufactured according to the results of optimization. The test proves that the improved NSGA - II algorithm can effectively reduce the total loss and the volume of deep-sea high-frequency transformer.
deep sea; high frequency transformer; multi-objective optimization;
2015-10-15
國家自然科學(xué)基金面上項(xiàng)目(51174087),湖南省教育廳資助項(xiàng)目(10C0681),湘潭市產(chǎn)業(yè)創(chuàng)新研究項(xiàng)目(潭財(cái)企發(fā)(2015)12號(hào)文)。
羅柏文(1975-),男,副教授。
1001-4500(2016)05-0051-06
TM412
A