雷飛 張?zhí)礻?/p>

摘 ? 要：針對(duì)常規(guī)模糊推理算法無法有效求解具有復(fù)雜或未知傳熱規(guī)律的傳熱學(xué)反問題，提出一種反饋-模糊推理全局算法. 基于模糊推理的基本流程，將變論域方法、反饋思想和模擬退火算法結(jié)合起來，通過反饋單元降低模糊規(guī)則對(duì)傳熱規(guī)律的依賴，并通過模擬退火單元防止陷入局部最優(yōu). 采用此方法對(duì)風(fēng)冷散熱器翅片的幾何結(jié)構(gòu)設(shè)計(jì)問題進(jìn)行求解，并與常規(guī)模糊推理算法和模擬退火算法結(jié)果進(jìn)行對(duì)比. 同時(shí)，對(duì)不同初始值、不同輸入誤差下反饋-模糊推理全局算法的計(jì)算結(jié)果進(jìn)行驗(yàn)證和對(duì)比. 結(jié)果表明，該算法可以解決傳熱學(xué)規(guī)律復(fù)雜或未知的傳熱學(xué)反問題，計(jì)算結(jié)果不受初始值影響. 該算法在解決此類問題時(shí)擁有良好的魯棒性和抗不適定性，可以為反問題、結(jié)構(gòu)設(shè)計(jì)和優(yōu)化提供參考.

關(guān)鍵詞：傳熱;反問題;模糊推理;純電動(dòng)汽車;全局優(yōu)化

中圖分類號(hào)：TK124 ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)志碼：A

A Feedback-fuzzy Inference Global Algorithm for Thermal Design

LEI Fei?，ZHANG Tianhao

（State Key Laboratory of Advanced Design and Manufacturing for Vehicle Body，Hunan University，Changsha 410082，China）

Abstract：Conventional fuzzy inference algorithms cannot effectively solve the inverse heat transfer problems with complex or unknown heat transfer laws. In this paper， a feedback-fuzzy inference global method is proposed. Based on the basic process of fuzzy inference， variable universe method， feedback thought and simulated annealing algorithm are combined. The dependence of fuzzy rules on heat transfer law is reduced by feedback unit， and the local optimum is prevented by simulated annealing unit. This method was employed to solve the geometric structural design problem of air-cooled radiator fins. The results were compared with those of conventional fuzzy inference and simulated annealing algorithm. And the results of feedback-fuzzy inference global method under different initial values and input errors were verified and compared. The results showed that the proposed method could solve the inverse heat transfer problem with complicated or unknown heat transfer law， and the results were not affected by initial values. The method demonstrates robustness and resistance to ill-posed problems， which can provide a reference for inverse problems， structural design and optimization.

Key words：heat transfer;inverse problem;fuzzy inference;battery electric vehicle;global optimization

電動(dòng)汽車電機(jī)控制器內(nèi)部的集成化程度越來越高，控制器內(nèi)元器件的發(fā)熱量不斷增加. 在有限的物理空間內(nèi)，多合一控制器中大功率元件的熱通量可達(dá)200 ～ 300 W/cm2. 過高的溫度會(huì)直接損壞元器件，熱設(shè)計(jì)和熱管理成為集成控制器設(shè)計(jì)的關(guān)鍵環(huán)節(jié)[1].

熱設(shè)計(jì)問題通常是已知目標(biāo)溫度來求解相關(guān)設(shè)計(jì)參數(shù)，本質(zhì)上屬于已知目標(biāo)溫度求解幾何邊界的傳熱學(xué)反問題[2]. 通?；趥鳠嵯到y(tǒng)的局部測(cè)量溫度信息，求解系統(tǒng)的熱流密度、熱源參數(shù)、熱物理參數(shù)和幾何參數(shù)等難以測(cè)量的未知參數(shù)的問題. 在控制器熱設(shè)計(jì)問題中，由于大量電子電器元件高度集成導(dǎo)致其傳熱規(guī)律復(fù)雜，難以得到傳熱學(xué)規(guī)律的定性描述. 復(fù)雜的傳熱規(guī)律是該問題的典型特點(diǎn)之一.

不適定性是傳熱學(xué)反問題的另一特點(diǎn). 這種不適定性是指該問題在Hadamard意義下解的唯一性、存在性和穩(wěn)定性不能同時(shí)得到滿足[3]. 一方面，由于輸入信息不全面，反問題的解可能不存在或不唯一;另一方面，由于輸入信息存在偏差，輸入信息的誤差在反問題的求解過程中會(huì)被放大，從而導(dǎo)致反問題的解往往存在不穩(wěn)定性. 反問題的不適定性為求解帶來了較大挑戰(zhàn).

傳熱學(xué)反問題的求解方法可分為非迭代法和迭代法. 非迭代法依據(jù)已知的測(cè)點(diǎn)信息和待求解量之間的數(shù)學(xué)關(guān)系，直接通過已知量計(jì)算待反演參數(shù)，無需反復(fù)計(jì)算正問題以修正計(jì)算結(jié)果，可實(shí)現(xiàn)較高的計(jì)算速度. 吳國鵬等[4]基于最小二乘法和誤差函數(shù)反演了功能梯度材料熱源發(fā)熱功率;Yu等[5]將最小二乘法和邊界元法相結(jié)合提高算法的抗不適定性并對(duì)鍋爐內(nèi)壁的幾何邊界進(jìn)行了反演;Yu等[6]利用邊界元法分別在穩(wěn)態(tài)下和瞬態(tài)下對(duì)爐內(nèi)壁的邊界條件進(jìn)行識(shí)別. 對(duì)于復(fù)雜傳熱學(xué)問題，由于已知條件和待求解量之間的數(shù)學(xué)關(guān)系建立存在困難，該類方法應(yīng)用有限，同時(shí)存在不適定性問題.

迭代法可分為梯度法和非梯度法，其中經(jīng)典的梯度法有最速下降法（Steepest Descent Method，SDM）、共軛梯度法（Conjugate Gradient Method，CGM）、Levenberg-Marquardt 法（L-M）等. Huang等[7]通過SDM成功反演了翅片換熱器的換熱系數(shù);Yang等[8]采用改進(jìn)的L-M法估算了傳熱學(xué)反問題中的導(dǎo)熱系數(shù);Xiong等[9]應(yīng)用CGM反求了二維傳熱系統(tǒng)中的熱流量;王堃等[10]用正則化CGM估計(jì)了二維穩(wěn)態(tài)傳熱邊界溫度分布.

梯度法可以解決部分傳熱學(xué)反問題，但當(dāng)傳熱規(guī)律較為復(fù)雜時(shí)會(huì)陷入局部最優(yōu)，且問題的可微性不能保證. 此外，如果測(cè)量信息不完整或存在測(cè)量誤差，梯度法的求解結(jié)果可能出現(xiàn)惡化[11].

非梯度法具有良好的抗不適定性能力和全局搜索能力，但計(jì)算量較大. 這類方法有遺傳算法、蟻群算法等，近年來隨著計(jì)算技術(shù)的發(fā)展而得到較廣泛應(yīng)用. Lee等人[12]利用粒子群優(yōu)化方法估計(jì)介質(zhì)的輻射特性;Liu[13]用一種改進(jìn)的遺傳算法估算瞬態(tài)傳熱系統(tǒng)的平板熱源功率;Parwani[14]等人利用差分進(jìn)化算法估算了傳熱系統(tǒng)中熱源的發(fā)熱功率和位置.

模糊推理是20世紀(jì)60年代Zadeh提出的一種啟發(fā)式算法[15]. 它是在模糊集合論基礎(chǔ)上對(duì)人腦的決策和判斷做定量研究的一種不確定推理方法. 該算法只需從單一起始點(diǎn)開始計(jì)算，相對(duì)于蟻群算法、遺傳算法等需要多個(gè)初始點(diǎn)的非梯度類算法而言計(jì)算成本較低，具有良好的抗不適定性，能有效抵抗輸入信息的干擾[16]. 此外，模糊推理算法可以有效利用不精確和不完整的信息[17]. 這些特點(diǎn)在解決包括傳熱學(xué)反問題在內(nèi)的不適定問題時(shí)具備優(yōu)勢(shì). Wang等[18]提出了分散模糊推理算法用來求解熱傳導(dǎo)反問題中的幾何參數(shù);Krzywanski等[19]利用模糊推理算法預(yù)測(cè)了大型循環(huán)流床的導(dǎo)熱系數(shù);Chen等[20]提出一種具有目標(biāo)跟蹤能力和降噪性能的智能加權(quán)模糊推理算法，并將其應(yīng)用于空心復(fù)合管的熱傳導(dǎo)幾何參數(shù)反演問題中;姜曙等[21]利用分散模糊推理對(duì)圓筒熱傳導(dǎo)問題的幾何邊界進(jìn)行了反演. 在以上傳熱規(guī)律明確的問題中，模糊推理方法均取得了較好的結(jié)果.

但是，常規(guī)模糊推理方法中模糊規(guī)則的制定對(duì)傳熱學(xué)規(guī)律的認(rèn)知依賴程度較高，在求解傳熱學(xué)規(guī)律未知或其定性描述為非單調(diào)的傳熱學(xué)反問題時(shí)，常規(guī)模糊推理算法會(huì)由于無法制定模糊規(guī)則而失效.

本文針對(duì)傳熱學(xué)規(guī)律未知或非單調(diào)的傳熱學(xué)反問題提出一種改進(jìn)的模糊推理算法，該算法以變論域模糊推理方法為基礎(chǔ)，結(jié)合反饋算法和模擬退火算法的優(yōu)勢(shì)，不僅降低對(duì)初始論域依賴性，同時(shí)降低模糊規(guī)則對(duì)傳熱規(guī)律認(rèn)知的依賴性. 使用該方法對(duì)電動(dòng)汽車集成控制器的風(fēng)冷散熱器進(jìn)行設(shè)計(jì)，并討論初始值和測(cè)量誤差對(duì)該方法的影響.

1 ? 算法描述

反饋-模糊推理全局算法由模糊推理單元，反饋單元和模擬退火單元組成.

該算法的輸入輸出為：輸入量e = Tmea0 - Tcal，其中Tmea0為測(cè)點(diǎn)溫度，Tcal為每次迭代過程的計(jì)算溫度. 對(duì)應(yīng)輸入、輸出論域的語言值有零（ZO）、正?。≒S）、正中（PM）、正大（PB）. 輸出量的大小由其對(duì)應(yīng)的論域和模糊規(guī)則決定，符號(hào)由反饋單元根據(jù)前一步迭代的計(jì)算結(jié)果決定. 這種輸出決定方式可以降低模糊推理輸出值對(duì)傳熱學(xué)規(guī)律的依賴.

若傳熱學(xué)反問題的傳熱規(guī)律較復(fù)雜或未知，則可以通過模擬退火單元防止結(jié)果陷入局部最優(yōu). 記uk（i）為模糊推理單元第i次迭代的輸出，n（i）為待求解量第i次迭代的計(jì)算結(jié)果，Tcal（i）為第i次迭代的計(jì)算溫度，迭代結(jié)束條件為Tmea0 - Tcal≤ε，ε為收斂條件. 算法的流程如圖1所示.

2 ? 問題描述

圖2為某電動(dòng)汽車集成控制器的簡(jiǎn)化幾何模型，該控制器箱體外形尺寸為388 mm×288 ?mm×124 mm，芯片為方形均勻體熱源，為模擬芯片封裝，在元件外部增加聚苯硫醚（Polyphenylene sulfide，PPS）材料外罩;用各向異性材料模擬印制電路板（Printed Circuit Board，PCB）;IGBT為導(dǎo)熱性良好的均勻方形體熱源;為模擬IGBT的水冷，將箱體底部設(shè)為333 K恒溫邊界，其余的面設(shè)為298 K恒溫邊界;箱體材料及翅片散熱器材料選用鋁材，各材料物性參數(shù)見表1;系統(tǒng)內(nèi)熱源及其發(fā)熱功率見表2;風(fēng)冷散熱器尺寸結(jié)構(gòu)如圖3所示.

未安裝風(fēng)冷散熱器時(shí)1、2、3號(hào)芯片溫度分別為392、411、389 K，其中2號(hào)芯片的溫度超過了設(shè)計(jì)要求的393 K，故對(duì)安裝在PCB下方風(fēng)冷散熱器的翅片數(shù)量進(jìn)行設(shè)計(jì)，使2號(hào)芯片溫度降至378 K. 該問題本質(zhì)上屬于已知測(cè)點(diǎn)溫度求解幾何邊界的傳熱學(xué)幾何反問題，即求解翅片數(shù)量n使得Tmea - Tcal≤ε，其中Tmea為2號(hào)芯片的目標(biāo)溫度，Tcal為每次迭代過程的計(jì)算溫度，ε為迭代停止閾值或收斂條件，由于該問題中的目標(biāo)溫度精度要求為1 K，故在此ε取值為1.

該傳熱系統(tǒng)中同時(shí)存在熱傳導(dǎo)、熱對(duì)流、熱輻射，假定在控制器內(nèi)的空氣為理想氣體，其流動(dòng)為三維湍流流動(dòng)，且忽略輻射換熱，則正問題的控制方程如下.

熱傳導(dǎo)方程：

[Δ]2T = 0 ? ? ?（1）

熱對(duì)流質(zhì)量守恒方程：

[Δ]V = 0 ? ? ?（2）

熱對(duì)流動(dòng)量守恒方程：

ρ = F - [Δ]P + η[Δ]2V ? ? ? ? ? （3）

熱對(duì)流能量守恒方程：

ρcp = λ[Δ]2T ? ? ? ? ? （4）

式中：V為流體速度;T為流體溫度;F為體積力;P為壓力;τ為時(shí)間;ρ、λ、η和cp分別為流體的密度、熱導(dǎo)率、動(dòng)力黏度和比熱容.

由式（1）～式（4）可知，由于該問題為熱傳導(dǎo)和熱對(duì)流共存的復(fù)合傳熱問題，傳熱規(guī)律較為復(fù)雜.

3 ? 反饋-模糊推理全局算法

3.1 ? 模糊推理單元

模糊推理是以模糊集合理論為基礎(chǔ)，結(jié)合計(jì)算機(jī)語言規(guī)則和先驗(yàn)知識(shí)，模仿人類思維進(jìn)行判決的不確定推理方法，其本質(zhì)是將一個(gè)給定輸入空間通過模糊邏輯映射到一個(gè)特定的輸出空間的計(jì)算過程.

模糊推理算法的具體模塊有模糊化、模糊規(guī)則的制定、模糊推理和解模糊化，如圖4所示.

3.1.1 ? 模糊化

模糊推理的輸入變量e = Tmea - Tcal，其取值范圍是[0，15]，單位為K. 翅片數(shù)量為n，其取值范圍為[5，25]，模糊推理的輸出變量uk 是n的補(bǔ)償量，取值范圍為[0，q]. q的取值方式如下：

q = 5e0.3，e > 2

1，e ≤ 2 ? ? ? ?（5）

上述可變輸出域可以保證輸出域的大小與輸入量e正相關(guān)，從而在e較大時(shí)提高迭代速度，在e較小時(shí)提高推理精度并且保證推理的穩(wěn)定性，減少初始域?qū)敵龅挠绊?

輸入變量e和輸出變量uk的論域分別劃分為4個(gè)模糊等級(jí)，對(duì)應(yīng)的語言值為零（ZO）、正?。≒S）、正中（PM）、正大（PB）. 采用三角形隸屬度函數(shù)，此時(shí)輸入論域、輸出論域及其對(duì)應(yīng)的隸屬度函數(shù)分別如圖5和圖6所示.

3.1.2 ? 模糊規(guī)則的制定

一般來說，帶翅片的風(fēng)冷散熱器的傳熱規(guī)律可定性描述為：當(dāng)?shù)装迕娣e一定，翅片數(shù)量在一定數(shù)量范圍內(nèi)時(shí)，散熱效果隨翅片數(shù)量增加而增強(qiáng);超過某一數(shù)量范圍，增加散熱翅片時(shí)，散熱效果會(huì)因風(fēng)阻而降低. 理論上，這是一個(gè)單極值問題，但由于本傳熱系統(tǒng)中的眾多其他熱源會(huì)對(duì)翅片散熱效果產(chǎn)生影響，因此傳熱規(guī)律不能簡(jiǎn)單地做定性描述. 這種情況下常規(guī)模糊推理會(huì)因模糊規(guī)則不明確而失效. 本文提出的反饋-模糊推理全局算法中模糊規(guī)則僅用來確定輸出值的大小，此時(shí)模糊規(guī)則的制定僅影響計(jì)算效率. 不考慮傳熱規(guī)律的模糊規(guī)則如下：

1）If e is PB then uk is PB，

2）If e is PM then uk is PM，

3）If e is PS then uk is PS，

4）If e is ZO then uk is ZO

3.1.3 ? 模糊推理和解模糊化

本文采用Mamdani法作為模糊推理方法，定義μA（e）、 μB（μk）分別為輸入變量和輸出變量對(duì)應(yīng)模糊集的隸屬度函數(shù)，μC（uk）為模糊推理結(jié)果，計(jì)算過程可表示為：

μC（uk） = max{min[ μA（e）μB（uk）]} ? ? ?（6）

解模糊化方法采用重心法，計(jì)算公式為：

uk =μC（uk）ukdu /μC（uk）du ? ? ?（7）

3.2 ? 反饋單元

根據(jù)反饋的基本思想，每次迭代的計(jì)算結(jié)果都將傳遞給該單元進(jìn)行判斷，若優(yōu)于當(dāng)前結(jié)果，則采納該結(jié)果;若劣于當(dāng)前結(jié)果，則重新選取補(bǔ)償方向. 利用反饋單元對(duì)補(bǔ)償量符號(hào)進(jìn)行判斷的代碼如下，其中r為隨機(jī)數(shù).

Start

if i > 1，

if Tcal（i） - Tcal（i - 1） < 0，

uk（i）/uk（i）= uk（i - 1）/uk（i - 1）;

else if Tcal（i） - Tcal（i - 1） > 0，

uk（i）/uk（i）= -uk（i - 1）/uk（i - 1）;

else if Tcal（i） - Tcal（i - 1） = 0，

uk（i） = 2（r - 0.5）;

n（i） = n（i - 1） + uk（i）;

if n（i） < -4，

n（i + 1） = -4 + 0.5r;

if n（i + 1） > 16，

n（i + 1） = 16 - r;

End

End

3.3 ? 模擬退火單元

模擬退火算法是一種基于蒙特卡羅迭代求解策略的隨機(jī)優(yōu)化方法. 該算法模擬了固體退火規(guī)律，結(jié)合溫度下降過程中的概率跳躍，隨機(jī)搜索目標(biāo)函數(shù)的全局最優(yōu)解. 當(dāng)問題的傳熱學(xué)規(guī)律定性描述為單極值問題時(shí)可不使用此單元. 由于本研究中反問題的傳熱學(xué)規(guī)律未知，為避免陷入局部最優(yōu)，故在反饋單元的基礎(chǔ)上加入模擬退火單元進(jìn)行全局求解.

設(shè)置當(dāng)前退火溫度為T，退火初始溫度 T2 = 90 ℃，退火結(jié)束溫度T3 = 89.9 ℃，降溫方式采用等比例降溫，比例系數(shù)j為0.999，結(jié)果接受準(zhǔn)則為：

當(dāng)Tcal（i） - Tcal（i-1） < 0時(shí)，接受當(dāng)前結(jié)果;當(dāng)Tcal（i） - Tcal（i-1） ≥ 0時(shí)，如果rand < e，接受結(jié)果，否則不接受.

3.4 ? 計(jì)算過程

Step1 ? 設(shè)置模糊推理的輸入量e與輸入論域，輸出量uk與輸出論域，初始化迭代次數(shù)i = 1;

Step2 ? 進(jìn)入模糊推理模塊求解uk（i），令n（i+1）= n（i）+ uk（i）;

Step3 ? 根據(jù)n（i+1）計(jì)算正問題，得到Tcal（i+1）和e（i+1）;

Step4 ? 判斷e（i + 1）≤1或T < T3是否滿足，如果滿足，終止迭代并輸出n（i + 1），否則i = i + 1，迭代繼續(xù);

Step5 ? 重復(fù)Step2—Step4，當(dāng)i ≥ 2時(shí)進(jìn)入反饋單元和模擬退火單元;

Step6 ? T2 = T2 × j;

Step7 ? 重復(fù)Step4;

Step8 ? 重復(fù)Step2—Step7.

4 ? 計(jì)算結(jié)果及討論

4.1 ? 計(jì)算結(jié)果及驗(yàn)證

當(dāng)初始翅片數(shù)為25時(shí)，使用反饋-模糊推理全局算法和模擬退火算法計(jì)算翅片數(shù)和目標(biāo)芯片溫度的迭代結(jié)果分別如圖7和圖8所示. 在相同的輸入、輸出論域條件下，使用常規(guī)模糊推理計(jì)算翅片數(shù)和目標(biāo)芯片溫度的迭代結(jié)果分別如圖9和圖10所示.

為了驗(yàn)證算法的準(zhǔn)確性，本文通過仿真計(jì)算了采用不同翅片數(shù)時(shí)芯片溫度的情況，如圖11所示.

由圖7和圖8可見，使用模擬退火算法得到翅片數(shù)量n=6.5，2號(hào)芯片溫度為378.6 K，取整數(shù)得n=7;使用反饋-模糊推理全局算法得到翅片數(shù)量n=6.9，2號(hào)芯片溫度為377.9 K，取整數(shù)得n=7. 反饋-模糊推理全局算法在第9次迭代時(shí)迭代終止，模擬退火算法在第19次迭代時(shí)終止，故反饋-模糊推理全局算法的計(jì)算效率高于模擬退火算法.

由圖9和圖10可見，使用常規(guī)模糊推理算法得到翅片數(shù)量n=22.34，2號(hào)芯片溫度為383.57 K，取整數(shù)得n=22. 由圖11可知n=22時(shí)2號(hào)芯片溫度為383.4 K，n=7時(shí)2號(hào)芯片溫度為377.9 K. 常規(guī)模糊推理算法在相同條件下無法得到正確結(jié)果.

由圖11可知，反饋-模糊推理全局算法結(jié)果準(zhǔn)確且有效避免了局部最優(yōu)，反演精度略高于模擬退火算法.

綜上所述，反饋-模糊推理全局算法可以有效求解傳熱學(xué)規(guī)律無法定量描述的傳熱學(xué)反問題，避免陷入局部最優(yōu)，其計(jì)算精度和迭代效率均高于設(shè)置相同的模擬退火算法. 利用反求結(jié)果得到的最佳翅片數(shù)量n=7進(jìn)行仿真，得到的芯片溫度云圖見圖12，3個(gè)芯片最高溫度均低于393 K，結(jié)果符合要求.

4.2 ? 初始值對(duì)結(jié)果的影響

使用反饋-模糊推理全局算法在初始翅片數(shù)n（0）分別為10和5時(shí)，計(jì)算翅片數(shù)的迭代結(jié)果和目標(biāo)芯片溫度的迭代結(jié)果分別如圖13和圖14所示.

由圖13可知，初始翅片數(shù)n（0）分別為10和5時(shí)的迭代結(jié)果均為n = 7，可見對(duì)于相同的設(shè)計(jì)問題，初始值只影響迭代次數(shù)，不影響迭代結(jié)果.

4.3 ? 溫度測(cè)量誤差對(duì)結(jié)果的影響

給2號(hào)芯片的溫度添加隨機(jī)測(cè)量誤差Tmea =Tmea + ξ × r，其中ξ是測(cè)量誤差的標(biāo)準(zhǔn)偏差，r是一個(gè)取值范圍在[-1，1]的隨機(jī)數(shù)，當(dāng)ξ分別為1%、3%、5%時(shí)，迭代結(jié)果如圖15所示. 由圖15可知，在計(jì)算過程中隨迭代次數(shù)的增加，誤差對(duì)迭代結(jié)果的影響逐漸增大. 當(dāng)ξ分別為1%、3%、5%時(shí)的最終迭代結(jié)果分別為n = 7.22、7.75、7.94，反演偏差在可接受范圍之內(nèi)，該算法可以有效對(duì)抗輸入信息誤差.

5 ? 結(jié) ? 論

本文針對(duì)常規(guī)模糊推理方法的局限性提出了一種結(jié)合模糊推理、反饋思想和模擬退火算法的反饋-模糊推理全局算法，并從傳熱學(xué)反問題的角度對(duì)電動(dòng)汽車集成控制器中風(fēng)冷散熱器的幾何結(jié)構(gòu)設(shè)計(jì)問題進(jìn)行了求解.

計(jì)算結(jié)果表明，相對(duì)常規(guī)模糊推理算法，反饋-模糊推理全局算法可以應(yīng)用于非單調(diào)或無法定性描述傳熱學(xué)規(guī)律的熱設(shè)計(jì)問題中，且具有迭代精度高、收斂速度快、計(jì)算結(jié)果受初始值影響小等特點(diǎn). 通過討論輸入誤差對(duì)計(jì)算結(jié)果的影響驗(yàn)證了該算法具有良好的魯棒性和抗不適定性.

該反饋-模糊推理全局算法拓展了模糊推理方法的應(yīng)用范圍，可以為反問題、結(jié)構(gòu)設(shè)計(jì)和優(yōu)化提供參考.

參考文獻(xiàn)

[1] ? ?楊東升，丁曉紅，季學(xué)榮，等. 基于自適應(yīng)成長(zhǎng)法的散熱通道構(gòu)建技術(shù)[J]. 中國機(jī)械工程，2012，23（12）：1404—1407.

YANG D S，DING X H，JI X R，et al. Construction technology of cooling channel based on adaptive growth method[J]. China Mechanical Engineering，2012，23（12）：1404—1407. （In Chinese）

[2] ? ?劉慧開，楊立，孫豐瑞. 基于傳熱反問題的異步電機(jī)參數(shù)估計(jì)方法研究[J]. 中國電機(jī)工程學(xué)報(bào)，2006，26（22）：151—156.

LIU H K，YANG L，SUN F R. Estimate method for parameters of asynchronous motor based on inverse heat conduction problem analysis[J]. Proceedings of the CSEE，2006，26（22）：151—156. （In Chinese）

[3] ? ?ALIFANOV O M. Inverse heat transfer problems[M]. Berlin：Spring-Verlag，1994：33—35.

[4] ? ?吳國鵬，徐闖，李英健，等. 一種新的熱源強(qiáng)度識(shí)別的非迭代方法研究[J].工程熱物理學(xué)報(bào)，2020，41（9）：2233—2240.

WU G P，XU C，LI Y J，et al. Research on a novel non-iterative inverse method for identifying the strength of heat source[J]. Journal of Engineering Thermophysics，2020，41（9）：2233—2240. （In Chinese）

[5] ? ?YU B，ZHOU H L，GAO Q，et al. Geometry boundary identification of the furnace inner wall by BEM without iteration[J]. Numerical Heat Transfer，2016，69（11）：1253—1262.

[6] ? ?YU B，YAO W A，GAO Q，et al. A novel non-iterative inverse method for estimating boundary condition of the furnace inner wall[J]. International Communications in Heat and Mass Transfer，2017，87：91—97.

[7] ? ?HUANG C H，YUAN I C，AY H. A three-dimensional inverse problem in imaging the local heat transfer coefficients for plate finned-tube heat exchangers[J]. International Journal of Heat and Mass Transfer，2003，46（19）：3629—3638.

[8] ? ?YANG K，JIANG G H，PENG H F，et al. A new modified Levenberg-Marquardt algorithm for identifying the temperature-dependent conductivity of solids based on the radial integration boundary element method[J]. International Journal of Heat and Mass Transfer，2019，144：118615.

[9] ? ?XIONG P，DENG J，LU T，et al. A sequential conjugate gradient method to estimate heat flux for nonlinear inverse heat conduction problem[J]. Annals of Nuclear Energy，2020，149：107798.

[10] ?王堃，王廣軍，陳紅，等. 穩(wěn)態(tài)傳熱邊界溫度分布的正則化共軛梯度反演[J]. 中國電機(jī)工程學(xué)報(bào)，2013，33（17）：78—82.

WANG K，WANG G J，CHEN H，et al. Regularization conjugate gradient inverse for the temperature distribution of steady-state heat transfer boundary[J]. Proceedings of the CSEE，2013，33（17）：78—82. （In Chinese）

[11] ?ELEGBEDE C. Structural reliability assessment based on particles swarm optimization[J].Structural Safety，2005，27（2）：171—186.

[12] ?LEE K H，BAEK S W，KIM K W. Inverse radiation analysis using repulsive particle swarm optimization algorithm[J].International Journal of Heat and Mass Transfer，2008，51（11/12）：2772—2783.

[13] ?LIU F B. A modified genetic algorithm for solving the inverse heat transfer problem of estimating plan heat source[J]. International Journal of Heat and Mass Transfer，2008，51（15/16）：3745—3752.

[14] ?PARWANI A K，TALUKDAR P，SUBBARAO P M V. Simultaneous estimation of strength and position of a heat source in a participating medium using DE algorithm[J]. Journal of Quantitative Spectroscopy and Radiative Transfer，2013，127：130—139.

[15] ?ZADEH L A. Fuzzy sets[J]. Information and Control，1965，8（3）：338—353.

[16] ?WANG G J，LUO Z M，ZHU L N，et al. Fuzzy estimation for temperature distribution of furnace inner surface[J]. International Journal of Thermal Sciences，2012，51：84—90.

[17] ?LI Y H，WANG G J，CHEN H，et al. A decentralized fuzzy inference method for the inverse geometry heat conduction problem[J]. Applied Thermal Engineering，2016，106：109—116.

[18] ?WANG G J，ZHU L N，CHEN H. A decentralized fuzzy inference method for solving the two-dimensional steady inverse heat conduction problem of estimating boundary condition[J]. International Journal of Heat and Mass Transfer，2011，54（13/14）：2782—2788.

[19] ?KRZYWANSKI J，NOWAK W. Modeling of bed-to-wall heat transfer coefficient in a large-scale CFBC by fuzzy logic approach[J].International Journal of Heat and Mass Transfer，2016，94：327—334.

[20] ?CHEN T C，LEE M H. Intelligent fuzzy weighted input estimation method applied to inverse heat conduction problems[J]. International Journal of Heat and Mass Transfer，2008，51（17/18）：4168—4183.

[21] ?姜曙，朱震宇，萬世斌. 基于邊界元和分散式模糊推理的圓筒內(nèi)壁形狀識(shí)別算法[J]. 世界科技研究與發(fā)展，2015，37（2）：122—126.

JIANG S，ZHU Z Y，WAN S B. Geometry estimation of furnace inner wall based on BEM and decentralized fuzzy inference method[J]. World Sci-Tech R & D，2015，37（2）：122—126. （In Chinese）

收稿日期：2021-01-19

基金項(xiàng)目：國家重點(diǎn)研發(fā)計(jì)劃資助項(xiàng)目（2018YFB0104501），National Key Research and Development Program of China（2018YFB0104501）

作者簡(jiǎn)介：雷飛（1981—），男，河南淅川人，湖南大學(xué)副教授，博士

通信聯(lián)系人，E-mail：fea_lei@163.com