主電動(dòng)機(jī)風(fēng)道優(yōu)化仿真

蒲業(yè)偉 林鵬 劉為亞

摘 要：利用試驗(yàn)設(shè)計(jì)，通過(guò)對(duì)影響風(fēng)量分配比、直（彎）腿的沿程壓力損失等考核指標(biāo)的多組設(shè)計(jì)變量進(jìn)行了靈敏度分析，指出了入口處導(dǎo)流板與風(fēng)道的夾角是影響風(fēng)量分配比和直腿壓力損失最敏感因素，為后續(xù)主風(fēng)道的優(yōu)化重點(diǎn)指明了方向。

關(guān)鍵詞：靈敏度分析；風(fēng)量分配比；壓力損失

1 前言

目前，已經(jīng)完成了主電動(dòng)機(jī)風(fēng)道（下文簡(jiǎn)稱主風(fēng)道）的流動(dòng)特性及兩分支風(fēng)量分配特性的CFD仿真分析。在此基礎(chǔ)上，為了進(jìn)一步提高主風(fēng)道的性能，可嘗試進(jìn)行主風(fēng)道結(jié)構(gòu)優(yōu)化，使兩分支風(fēng)量的分配比例更接近1：1，同時(shí)兩分支沿程壓力損失更小。解決此類工程結(jié)構(gòu)優(yōu)化問(wèn)題通常采用基于試驗(yàn)設(shè)計(jì)（DOE，Design of Experiment）的響應(yīng)面（Response Surface）優(yōu)化設(shè)計(jì)[1]，如果問(wèn)題涉及的設(shè)計(jì)變量（也稱為設(shè)計(jì)變量或因子）較多，還需先對(duì)各變量進(jìn)行靈敏度分析，以確定各變量對(duì)響應(yīng)影響程度的主次順序，為下一步細(xì)致擬合主要設(shè)計(jì)變量與考核指標(biāo)的近似響應(yīng)面奠定基礎(chǔ)。

文章基于主風(fēng)道的CFD模型，利用試驗(yàn)設(shè)計(jì)（DOE，Design of Experiment）進(jìn)行風(fēng)道內(nèi)導(dǎo)流板的結(jié)構(gòu)對(duì)風(fēng)量分配比和兩分支沿程壓力損失的靈敏度分析，各導(dǎo)流板編號(hào)如圖1所示。涉及的設(shè)計(jì)變量有各導(dǎo)流板的長(zhǎng)度（其中忽略導(dǎo)流板4的長(zhǎng)度變化）及其與風(fēng)道的夾角，共計(jì)11個(gè)設(shè)計(jì)變量。所有相關(guān)的CFD仿真在Star-CCM+

中完成，DOE分析在Hyperstudy 中完成。

圖1 主風(fēng)道及各導(dǎo)流板的原結(jié)構(gòu)示意圖

2 部分變量試驗(yàn)設(shè)計(jì)方法

2.1 靈敏度分析

靈敏度分析是最優(yōu)化設(shè)計(jì)的重要組成部分，是研究與分析一個(gè)系統(tǒng)（或模型）的狀態(tài)或輸出變化對(duì)系統(tǒng)參數(shù)或周圍條件變化敏感程度的方法。通過(guò)靈敏度分析，可以研究原始數(shù)據(jù)不準(zhǔn)確或發(fā)生變化時(shí)最優(yōu)解的穩(wěn)定性，還可以決定哪些參數(shù)對(duì)系統(tǒng)或模型有較大的影響。對(duì)文章而言，直接利用DOE獲取11個(gè)設(shè)計(jì)變量對(duì)考核指標(biāo)的近似響應(yīng)面需要進(jìn)行大量的試驗(yàn)仿真（除去用于校驗(yàn)和修正的輔助計(jì)算，擬合二次響應(yīng)面至少需要計(jì)算77次，擬合三次響應(yīng)面則至少需要計(jì)算198次）。而通過(guò)靈敏度分析找出影響考核指標(biāo)最大或較大的設(shè)計(jì)變量，可極大節(jié)省計(jì)算耗費(fèi)，提高計(jì)算效率和精度。

2.2 部分變量試驗(yàn)設(shè)計(jì)

進(jìn)行DOE最直觀的方法是全變量試驗(yàn)設(shè)計(jì)，即考慮所有可能的試驗(yàn)條件組合，分析全部設(shè)計(jì)變量以及各變量水平的相互作用情況。但該方法的缺點(diǎn)也很明顯，因?yàn)殡S著變量個(gè)數(shù)或者變量變化水平的增加，所需的試驗(yàn)次數(shù)將成倍增多。不過(guò)，實(shí)際經(jīng)驗(yàn)表明，僅進(jìn)行全變量方法所要求的一部分試驗(yàn)便可以得到主效應(yīng)和低階（通常為二階，即兩變量之間）相互作用，這類試驗(yàn)稱為部分變量試驗(yàn)設(shè)計(jì)。很明顯，應(yīng)用部分變量設(shè)計(jì)可以在滿足主要求的同時(shí)，顯著降低成本和節(jié)省時(shí)間，因而這種方法在工業(yè)中得到了較廣泛的應(yīng)用。

2.3 Hyperstudy 介紹

Hyperstudy 是Hyperworks 軟件包中的一款主要產(chǎn)品，主要用于 CAE 環(huán)境下的DOE、曲線擬合、優(yōu)化及隨機(jī)分析的研究，可與多款外部求解器（如Abaqus、Radioss、Ansys等）合并使用。Hyperstudy中提供的試驗(yàn)設(shè)計(jì)方法，常用的包括全變量設(shè)計(jì)（Full Factorial Design）、部分變量設(shè)計(jì)（Fractional Factorial Design）、中心組合設(shè)計(jì)（Central Composite Design）以及拉丁超立方設(shè)計(jì)（Latin Hypercube Design）等[2]。

3 靈敏度分析試驗(yàn)設(shè)計(jì)

3.1 主風(fēng)道的CFD模型

為方便實(shí)現(xiàn)DOE中反復(fù)的模型修改，需簡(jiǎn)化用于校核方案設(shè)計(jì)的原CFD模型。簡(jiǎn)化模型主要忽略了導(dǎo)流板的厚度以及與風(fēng)道連接部分的材料，同時(shí)采用baffles屬性的交界面（interface）模擬導(dǎo)流板，如圖2所示。經(jīng)過(guò)驗(yàn)算發(fā)現(xiàn)，兩模型所得結(jié)果之間差別較小，各考核指標(biāo)的相對(duì)誤差在5%以內(nèi)，說(shuō)明對(duì)原模型的簡(jiǎn)化比較合理。

圖2 主風(fēng)道的CFD模型

3.2 試驗(yàn)設(shè)計(jì)參數(shù)

設(shè)計(jì)變量：各導(dǎo)流板的長(zhǎng)度比（與原導(dǎo)流板長(zhǎng)度的比值，0.9～1.1內(nèi)變化）以及各導(dǎo)流板與風(fēng)道的夾角（當(dāng)?shù)乜諝庵髁鲃?dòng)方向?yàn)檗D(zhuǎn)動(dòng)軸正方向，遵循右手定則，且順時(shí)針轉(zhuǎn)動(dòng)為正，±45°內(nèi)變化）。文章進(jìn)行靈敏度的部分變量試驗(yàn)設(shè)計(jì)時(shí)，兩設(shè)計(jì)變量的水平均為3，如表1所示，各設(shè)計(jì)變量名稱分別用英文字母A～K代替。

考核指標(biāo)：共有三個(gè)，兩出口風(fēng)量分配比（直腿/彎腿）、直腿沿程壓力損失（Pa）、彎腿沿程壓力損失（Pa），其中風(fēng)量分配比是首要考核指標(biāo)。

試驗(yàn)方案：利用Hyperstudy 生成的方案如表2所示，根據(jù)各方案給出的變量數(shù)值組合，在Hypermesh 中調(diào)整原模型的網(wǎng)格結(jié)構(gòu)，然后倒入Star-ccm+ 中進(jìn)行流體仿真和結(jié)果提取。文章需要進(jìn)行27次試驗(yàn)仿真才可得到全部變量對(duì)指標(biāo)的主效應(yīng)，模型修改和計(jì)算的總耗時(shí)為10天。

表1 設(shè)計(jì)變量與水平

3.3 試驗(yàn)結(jié)果及分析

將得到的27次試驗(yàn)仿真結(jié)果填入表2，然后采用極差法（又稱直觀分析法或R法，所謂極差，就是平均效果中最大值和最小值的差）對(duì)表2進(jìn)行統(tǒng)計(jì)分析，首先得到如圖3所示的各設(shè)計(jì)變量和可靠指標(biāo)的數(shù)值分布情況，然后基于圖3便可得到各設(shè)計(jì)變量對(duì)風(fēng)量分配比等考核指標(biāo)的主效應(yīng)曲線，如圖4所示。以圖4中的變量A（導(dǎo)流板1與風(fēng)道的夾角）為例，在（-45°～0°）之間，風(fēng)量分配比隨變量A的增大而減小，而在（0°～45°）之間，風(fēng)量分配比則隨變量A的增加而增加。另外，宏觀整個(gè)圖4，不難看出變量J（導(dǎo)流板6與風(fēng)道的夾角）對(duì)風(fēng)量分配比的影響最明顯，而變量B（導(dǎo)流板1的長(zhǎng)度比）對(duì)該指標(biāo)的影響最不明顯。同樣，對(duì)影響直腿壓力損失及彎腿壓力損失的設(shè)計(jì)變量可按同樣規(guī)律分析，文章不再詳述。

最后可對(duì)27組實(shí)驗(yàn)數(shù)據(jù)進(jìn)行多項(xiàng)式擬合（文章各指標(biāo)的擬合誤差均在±3%以內(nèi)），得到靈敏度分析結(jié)果如圖5所示，可知影響風(fēng)量分配比和直腿壓力損失最敏感的兩變量分別是J和A，影響彎腿壓力損失最敏感的兩變量分別是F和E。鑒于風(fēng)量分配比是首要考核指標(biāo)，因此接下來(lái)對(duì)主風(fēng)道的結(jié)構(gòu)優(yōu)化可著重優(yōu)化導(dǎo)流板6與風(fēng)道的夾角。

4 結(jié)束語(yǔ)

文章將部分變量試驗(yàn)設(shè)計(jì)方法用于主風(fēng)道結(jié)構(gòu)優(yōu)化設(shè)計(jì)的靈敏度分析，綜合考慮了影響風(fēng)量分配比、直腿及彎腿的沿程壓力損失等考核指標(biāo)的多個(gè)設(shè)計(jì)變量。通過(guò)試驗(yàn)仿真和統(tǒng)計(jì)分析，確定了影響各考核指標(biāo)的最敏感或較敏感設(shè)計(jì)變量，指出入口處導(dǎo)流板與風(fēng)道的夾角是影響風(fēng)量分配比這一首要指標(biāo)最大的設(shè)計(jì)變量，從而使后續(xù)主風(fēng)道結(jié)構(gòu)的優(yōu)化更有針對(duì)性，并可為今后類似結(jié)構(gòu)的風(fēng)道設(shè)計(jì)提供依據(jù)。

參考文獻(xiàn)

[1]郝琪，楊林松.基于HYPERSTUDY的車身剛度響應(yīng)面優(yōu)化設(shè)計(jì)[J].2009中國(guó)汽車工程學(xué)會(huì)年會(huì)論文集，2009.

[2]洪清泉.OptiStruct&HyperStudy理論基礎(chǔ)與工程應(yīng)用[M].北京：機(jī)械工業(yè)出版社，2013.