楊偉濤

(惠而浦(中國)股份有限公司,合肥 230088)

0 引 言

隨著科學(xué)技術(shù)的發(fā)展和人們生活水平的提高,大家對家用電器的性能等要求越來越高。另外,隨著工業(yè)自動化、航空航天、武器裝備、艦船潛艇等不斷發(fā)展,人們對電機(jī)性能以及節(jié)能環(huán)保的要求也越來越高,因此為了進(jìn)一步提高電機(jī)性能,對電機(jī)的多變量多目標(biāo)優(yōu)化設(shè)計(jì)就變得尤為重要。

ANSYS Maxwell有限元仿真軟件中的參數(shù)化掃描優(yōu)化方法是近年來在電機(jī)設(shè)計(jì)中常用的一種方法。特別是近幾年比較火熱的OptiSlang。文獻(xiàn)[3]采用此方法對電機(jī)進(jìn)行了優(yōu)化,但需要建立龐大的參數(shù)化模型,并進(jìn)行大量且復(fù)雜的運(yùn)算才能得到所需的優(yōu)化方案。過程操作繁雜,對運(yùn)算設(shè)備要求較高,且需要大量的計(jì)算時(shí)間。

Taguchi法是一種較為新穎的局部優(yōu)化算法,對比其他局部優(yōu)化算法,它可實(shí)現(xiàn)多目標(biāo)優(yōu)化,并且在運(yùn)用正交實(shí)驗(yàn)的實(shí)驗(yàn)方法后,能在較少的實(shí)驗(yàn)次數(shù)下給出多目標(biāo)設(shè)計(jì)的最優(yōu)參數(shù)組合。文獻(xiàn)[4] 采用此方法對電機(jī)進(jìn)行了優(yōu)化,但沒有給出各變量與優(yōu)化目標(biāo)間的定量關(guān)系,以及在預(yù)測刻畫模型的幫助下推導(dǎo)出精準(zhǔn)的最優(yōu)變量組合,給出最佳方案。

本文以一臺9槽6極洗碗機(jī)用水泵電機(jī)為例,選取了繞組匝數(shù)、鐵心高度、槽口大小、氣隙直徑、氣隙大小、定子齒寬、定子齒高7個(gè)參數(shù)為變量,以交直軸電感、反電動勢為優(yōu)化目標(biāo),通過基于實(shí)驗(yàn)設(shè)計(jì)(DOE)思想的多目標(biāo)多變量的電機(jī)快速參數(shù)化優(yōu)化方法,實(shí)現(xiàn)了通過16次的部分析因?qū)嶒?yàn)設(shè)計(jì)完成了常規(guī)參數(shù)化掃描最少2 187個(gè)節(jié)點(diǎn)的設(shè)計(jì)目標(biāo),大大節(jié)約了項(xiàng)目資源,提高了設(shè)計(jì)人員的設(shè)計(jì)效率,并縮短了項(xiàng)目的交付周期。另外,還科學(xué)地給出了各變量與優(yōu)化目標(biāo)間的定量關(guān)系,提高了設(shè)計(jì)人員對產(chǎn)品的理解以及經(jīng)驗(yàn)設(shè)計(jì)和理論設(shè)計(jì)的水平。

1 DOE在電機(jī)多目標(biāo)多變量參數(shù)化分析上的應(yīng)用

基于實(shí)驗(yàn)設(shè)計(jì)(以下簡稱DOE)的電機(jī)多目標(biāo)多變量參數(shù)化分析步驟如下:

1)選定優(yōu)化目標(biāo),并基于理論分析、設(shè)計(jì)經(jīng)驗(yàn)、專家意見等確定與優(yōu)化目標(biāo)相對應(yīng)的參數(shù)變量。

2)根據(jù)確定的參數(shù)變量,制作實(shí)驗(yàn)因子設(shè)置表,表中需要給出各實(shí)驗(yàn)因子的上下水平,此上下水平的設(shè)置一般基于理論公式推導(dǎo)、設(shè)計(jì)經(jīng)驗(yàn)、專家意見等。同時(shí)表中還需要給出各因子與優(yōu)化目標(biāo)之間的理論邏輯關(guān)系,以及最終的目標(biāo)變化趨勢。

3)根據(jù)以上確定的參數(shù)變量,制作各參數(shù)變量的實(shí)驗(yàn)因子關(guān)系圖,用以指導(dǎo)后期的實(shí)驗(yàn)工作,并且便于后期相關(guān)人員了解當(dāng)時(shí)的實(shí)驗(yàn)策略。

4)根據(jù)實(shí)驗(yàn)因子關(guān)系圖進(jìn)行相應(yīng)實(shí)驗(yàn),并正確記錄相應(yīng)數(shù)據(jù)。

5)通過專業(yè)的數(shù)據(jù)分析軟件對實(shí)驗(yàn)數(shù)據(jù)進(jìn)行分析處理,找出影響優(yōu)化目標(biāo)的關(guān)鍵因子,并確定其相互間的定量關(guān)系,最終通過預(yù)測刻畫器給出最優(yōu)方案。

6)對最優(yōu)方案進(jìn)行仿真或樣機(jī)驗(yàn)證,確認(rèn)是否滿足優(yōu)化目標(biāo)。

本文案例是基于一款洗碗機(jī)用水泵電機(jī),原始方案如圖1所示,客戶要求在保證電機(jī)槽極比、定子鐵心外徑、反電動勢(128 V)等參數(shù)不變的前提下將電感從98 mH降低到72 mH,以達(dá)到提高電機(jī)性能和最大輸出功率的目的。

圖1 原始方案

首先,通過以往設(shè)計(jì)經(jīng)驗(yàn)以及理論公式推導(dǎo),找出影響優(yōu)化目標(biāo)的相關(guān)因子。表述反電動勢與各參數(shù)間相互關(guān)系的公式推導(dǎo)過程如下:

E=4.44fkwNphΦ

(1)

(2)

將式(2)代入式(1),得到式(3):

(3)

式中:μ0為氣隙磁導(dǎo)率;f為頻率;i為電流;L為氣隙高度;r為氣隙半徑;kw為繞組系數(shù);Nph為每相串聯(lián)總匝數(shù);g為氣隙大小;p為極對數(shù)。

由式(3)可知,在不考慮飽和以及在極對數(shù)、頻率、繞組系數(shù)、電流等參數(shù)一定的情況下,反電動勢的大小與氣隙高度、氣隙半徑成正比,與每項(xiàng)串聯(lián)總匝數(shù)的平方成正比,與氣隙大小成反比。

表述交直軸電感與各參數(shù)間相互關(guān)系的公式推導(dǎo)過程如下:

(4)

(5)

將式(5)代入式(4),得到式(6):

(6)

式中:Lq為交軸電感;Ld為直軸電感;La1為漏電感分量;Laa0為定子自感的氣隙分量。

由式(6)可知,在不考慮飽以及在極對數(shù)、繞組系數(shù)等參數(shù)一定的情況下,交直軸電感的大小與氣隙高度、氣隙半徑成正比,與每項(xiàng)串聯(lián)總匝數(shù)的平方成正比,與氣隙大小成反比,同時(shí)受漏電感分量大小的影響。

考慮到受磁路飽和、槽滿率、漏磁通等的影響,以及由式(3)和式(6)推導(dǎo)分析得出的結(jié)論,本實(shí)驗(yàn)設(shè)計(jì)將繞組匝數(shù)、定子鐵心高度、氣隙直徑、氣隙大小、定子齒寬、定子齒高、槽口大小作為實(shí)驗(yàn)因子,具體實(shí)驗(yàn)因子的高低水平設(shè)置、理論趨勢預(yù)測等如表1所示。

基于以上分析,并考慮到實(shí)驗(yàn)資源和項(xiàng)目周期的影響,本實(shí)驗(yàn)設(shè)計(jì)采用7因子2水平16次的部分析因?qū)嶒?yàn)策略,另外,根據(jù)正交實(shí)驗(yàn)設(shè)計(jì)思想,在保證各因子上下水平設(shè)置平衡的情況下設(shè)計(jì)出了具體的實(shí)驗(yàn)因子關(guān)系圖(FRD),如圖2所示。

根據(jù)以上實(shí)驗(yàn)策略,通過專業(yè)的電機(jī)設(shè)計(jì)軟件,分別進(jìn)行上述16次的電機(jī)設(shè)計(jì)仿真,得出的實(shí)驗(yàn)數(shù)據(jù),具體如表2所示。

圖2 實(shí)驗(yàn)因子關(guān)系圖

表2 實(shí)驗(yàn)數(shù)據(jù)

基于以上實(shí)驗(yàn)數(shù)據(jù),通過專業(yè)的數(shù)據(jù)分析軟件,如Jmp、Minitab等,對交軸電感Lq和線線間反電動勢E_ll的實(shí)驗(yàn)數(shù)據(jù)分別進(jìn)行擬合模型分析,得出如圖3、圖4所示的參數(shù)估計(jì)值。

圖3 部分析因?qū)嶒?yàn)電感排序后的參數(shù)估計(jì)值

圖4 部分析因?qū)嶒?yàn)反電動勢排序后的參數(shù)估計(jì)值

如圖3所示,根據(jù)P小于等于0.05的原則,我們可以得出影響電感的關(guān)鍵因子從大到小分別為定子鐵心高度、繞組匝數(shù)和槽口大小,其中交軸電感Lq與定子鐵心高度、繞組匝數(shù)成正相關(guān),預(yù)測公式系數(shù)分別為23.56和20.17,與槽口大小成負(fù)相關(guān),預(yù)測公式系數(shù)為-18.81。

如圖4所示,根據(jù)P小于等于0.05的原則,我們可以得出影響反電動勢的關(guān)鍵因子從大到小分別為定子鐵心高度、繞組匝數(shù)、氣隙直徑和氣隙大小,其中線線間反電動勢E_ll與定子鐵心高度、繞組匝數(shù)、氣隙直徑成正相關(guān),預(yù)測公式系數(shù)分別為26.51、14.28和8.71,與氣隙大小成負(fù)相關(guān),預(yù)測公式系數(shù)為-7.87。

通過數(shù)據(jù)分析,我們還得到了交軸電感Lq和線線間反電動勢E_ll與各實(shí)驗(yàn)因子的交互作用刻畫圖,如圖5和圖6所示。由此可知,各實(shí)驗(yàn)因子間無交互作用,可以不用考慮各因子間交互作用的影響。

圖5 電感各因子間交互作用刻畫圖

圖6 反電動勢各因子間交互作用刻畫圖

根據(jù)以上分析,考慮到各因子的預(yù)測系數(shù)以及成本因素,通過如圖7所示的部分析因?qū)嶒?yàn)預(yù)測刻畫器可知,要得到較大的反電動勢、較小的電感值,即E_ll-Lq取較大值時(shí),其最關(guān)鍵因子為槽口大小,且與其成正相關(guān),因此槽口大小B0的系數(shù)取1,由表1實(shí)驗(yàn)因子設(shè)置可知,槽口大小B0系數(shù)取1時(shí)即B0=6 mm。

根據(jù)圖7的預(yù)測刻畫器,我們可以計(jì)算出,在繞組匝數(shù)取470匝、定子鐵心高度14 mm、氣隙直徑40 mm、氣隙大小1.2 mm、定子齒寬3.3 mm、定子齒高12.5 mm時(shí),所得方案滿足客戶要求,此時(shí)Lq=69.8 mH,E_ll=128.6 V。

圖7 部分析因?qū)嶒?yàn)預(yù)測刻畫器

為驗(yàn)證7因子2水平16次的部分析因?qū)嶒?yàn)策略的可信度,本文還進(jìn)行了7因子2水平128次的完全析因?qū)嶒?yàn),并對實(shí)驗(yàn)數(shù)據(jù)進(jìn)行分析,得到了完全析因?qū)嶒?yàn)排序后的參數(shù)估計(jì)值,如圖8、圖9所示,我們發(fā)現(xiàn),影響電感的關(guān)鍵因子從大到小分別為定子鐵心高度、繞組匝數(shù)和槽口大小,其中與定子鐵心高度、繞組匝數(shù)成正相關(guān),預(yù)測公式系數(shù)分別為23.87和20.25,與槽口大小成負(fù)相關(guān),預(yù)測公式系數(shù)為-18.97。

圖8 完全析因?qū)嶒?yàn)電感排序后的參數(shù)估計(jì)值

影響反電動勢的關(guān)鍵因子從大到小分別為定子鐵心高度、繞組匝數(shù)、氣隙直徑和氣隙大小,其中與定子鐵心高度、繞組匝數(shù)、氣隙直徑成正相關(guān),預(yù)測公式系數(shù)分別為26.7、14.49和8.64,與氣隙大小成負(fù)相關(guān),預(yù)測公式系數(shù)為-8.01。以上結(jié)果與部分析因?qū)嶒?yàn)策略下得到的結(jié)果基本相同。

圖9 完全析因?qū)嶒?yàn)反電動勢排序后的參數(shù)估計(jì)值

另外,根據(jù)完全析因?qū)嶒?yàn)的預(yù)測刻畫圖,如圖10所示,同樣在槽口大小6 mm繞組匝數(shù)470匝、定子鐵心高度14 mm、氣隙直徑40 mm、氣隙大小1.2 mm、定子齒寬3.3 mm、定子齒高12.5 mm時(shí),Lq=69.1 mH,E_ll=128 V,其結(jié)果與部分析因?qū)嶒?yàn)策略下得到的結(jié)果基本相同。

圖10 完全析因?qū)嶒?yàn)預(yù)測刻畫器

為了驗(yàn)證實(shí)驗(yàn)設(shè)計(jì)思想在設(shè)計(jì)領(lǐng)域應(yīng)用的可行性,作者還進(jìn)行了常規(guī)的7個(gè)變量,每個(gè)變量設(shè)置3段,共2 187個(gè)節(jié)點(diǎn)的參數(shù)化掃描仿真實(shí)驗(yàn),其結(jié)果如圖11和圖12所示。當(dāng)槽口大小6 mm、繞組匝數(shù)470匝、定子鐵心高度14 mm、氣隙直徑40 mm、氣隙長度1.2 mm、定子齒寬3.3 mm、定子齒高12.5 mm時(shí),Lq=68.9 mH,E_ll=129.8 V,與基于實(shí)驗(yàn)設(shè)計(jì)思想下的多參數(shù)化多目標(biāo)仿真結(jié)果基本相同,證明了實(shí)驗(yàn)設(shè)計(jì)思想在設(shè)計(jì)領(lǐng)域應(yīng)用的可行性。

圖11 線間反電動勢

圖12 交軸電感

2 實(shí)驗(yàn)設(shè)計(jì)分析結(jié)果與樣機(jī)實(shí)測結(jié)果對比分析

根據(jù)16次部分析因?qū)嶒?yàn)設(shè)計(jì)、128次完全析因?qū)嶒?yàn)設(shè)計(jì)以及常規(guī)2 187個(gè)節(jié)點(diǎn)的參數(shù)化掃描,其所得結(jié)果基本相同,最終我們確定最優(yōu)方案的相關(guān)參數(shù)取值如下:槽口大小6 mm、繞組匝數(shù)470匝、定子鐵心高度14 mm、氣隙直徑40 mm、氣隙長度1.2 mm、定子齒寬3.3 mm、定子齒高12.5 mm。根據(jù)以上參數(shù)取值,我們制作了相應(yīng)的手板樣機(jī),如圖13所示。并將實(shí)測結(jié)果和仿真結(jié)果整理如表3所示。

圖13 樣機(jī)

表3 DOE分析結(jié)果與樣機(jī)實(shí)測結(jié)果對比

由表3可以看出,16次部分析因?qū)嶒?yàn)設(shè)計(jì)結(jié)果、128次完全析因?qū)嶒?yàn)設(shè)計(jì)結(jié)果、常規(guī)2 187個(gè)節(jié)點(diǎn)的參數(shù)化掃描結(jié)果與樣機(jī)實(shí)測結(jié)果的誤差均小于>5%,說明基于實(shí)驗(yàn)設(shè)計(jì)思想下的電機(jī)多目標(biāo)多變量的快速參數(shù)化設(shè)計(jì)分析的方法是可靠的,并且在保證同等設(shè)計(jì)精度的同時(shí)可以節(jié)約大量的設(shè)計(jì)資源和設(shè)計(jì)時(shí)間,大大提高了設(shè)計(jì)人員的設(shè)計(jì)效率,縮短了項(xiàng)目的交付周期,同時(shí),本方法還科學(xué)的給出了各因子間的定量關(guān)系,提高了設(shè)計(jì)人員的理論和經(jīng)驗(yàn)設(shè)計(jì)水平。

3 結(jié) 語

通過以上理論分析和實(shí)際測量結(jié)果對比,我們可以得出以下結(jié)論。

(1)基于DOE的電機(jī)多目標(biāo)多變量的參數(shù)化優(yōu)化方法可以快速地完成參數(shù)化尋優(yōu)設(shè)計(jì),方法是可行的,結(jié)果是準(zhǔn)確的;

(2)科學(xué)地給出各因子與優(yōu)化目標(biāo)間的定量關(guān)系,提高設(shè)計(jì)人員對產(chǎn)品的理解以及經(jīng)驗(yàn)設(shè)計(jì)和理論設(shè)計(jì)的水平。

(3)本文案例用16次的部分析因?qū)嶒?yàn)設(shè)計(jì)策略完成了常規(guī)參數(shù)化掃描2 187個(gè)節(jié)點(diǎn)的設(shè)計(jì)目標(biāo),大大節(jié)約了開發(fā)資源,提高了設(shè)計(jì)人員的設(shè)計(jì)效率,縮短了項(xiàng)目的交付周期。

(4)由于只需要較少的實(shí)驗(yàn)次數(shù)就可以完成多變量的尋優(yōu)目標(biāo),在不建立龐大的二維或三維參數(shù)化模型的情況下,可以采用二維有限元,甚至是三維有限元的設(shè)計(jì)仿真方法進(jìn)行多目標(biāo)多變量的參數(shù)化分析,大大提高了設(shè)計(jì)精度。

(5)由于其只需要較少的實(shí)驗(yàn)次數(shù)就可以完成多變量的實(shí)驗(yàn)?zāi)繕?biāo),所以對設(shè)計(jì)仿真的硬件要求不高,可大大降低公司對設(shè)計(jì)仿真硬件的投資。