夏榮菲， 陳宜飛， 馮勇建

(廈門大學(xué) 航空航天學(xué)院, 福建 廈門 361005)

真空絕熱板是快速發(fā)展起來的一種新型絕熱材料,由填充芯材、氣體阻隔膜以及吸氣劑構(gòu)成，主要通過板內(nèi)真空度來降低空氣引起的熱傳遞，其導(dǎo)熱系數(shù)可低達(dá)0.004 W/(m·K)左右，僅為傳統(tǒng)保溫材料的1/10. 由于這些特性，VIP已經(jīng)廣泛應(yīng)用于建筑[1]、交通運(yùn)輸業(yè)[2]、家用電器[3]等領(lǐng)域. 隨著對(duì)VIP的需求量的不斷增長(zhǎng)，VIP的質(zhì)量和壽命檢測(cè)變得至關(guān)重要. 導(dǎo)熱系數(shù)是反應(yīng)VIP性能的重要物理量，對(duì)其精確檢測(cè)有利于為設(shè)備的使用提供選擇依據(jù).

為了保證VIP的低導(dǎo)熱系數(shù)，在出廠前就需要對(duì)VIP導(dǎo)熱系數(shù)進(jìn)行精確測(cè)量. VIP導(dǎo)熱系數(shù)受材料的結(jié)構(gòu)、密度以及溫度等因素的影響，且VIP在長(zhǎng)時(shí)間使用后，或者遭到破壞時(shí)其真空度就會(huì)遭到破壞，從而導(dǎo)致導(dǎo)熱系數(shù)變大，使VIP絕熱性能受到影響. 因此，為了檢測(cè)VIP的老化程度，也需要對(duì)VIP導(dǎo)熱系數(shù)進(jìn)行測(cè)量.

目前VIP行業(yè)主要通過平板熱流計(jì)法檢測(cè)其導(dǎo)熱系數(shù)以實(shí)現(xiàn)VIP的性能檢測(cè). 該方法是基于熱傳導(dǎo)原理，通過熱流計(jì)的方法推導(dǎo)出VIP的導(dǎo)熱系數(shù)，具有檢測(cè)精確度高，檢測(cè)周期長(zhǎng)等特點(diǎn)[4]. 采用該方法已經(jīng)有比較成熟的檢測(cè)儀器，比如德國(guó)NETZSCH公司研發(fā)的面向真空絕熱板HFM436系列的熱流導(dǎo)熱系數(shù)儀[5], 其測(cè)量精確度可達(dá)±(1～3)%，重復(fù)性為0.25%，是基于平板熱熱流計(jì)法中檢測(cè)速度最快的測(cè)量方法，樣品測(cè)量時(shí)間短于5 min. 查爾姆斯理工大學(xué)采用TPS傳感器并建模，在40 s內(nèi)就可以測(cè)量出VIP的導(dǎo)熱系數(shù)，偏差可控制在1.3%以內(nèi)[6]. 此方法提高了測(cè)量速度，但以提高成本為代價(jià). 采用平板熱熱流計(jì)法的儀器結(jié)構(gòu)設(shè)計(jì)合適，結(jié)果準(zhǔn)確度高，但是需要測(cè)試時(shí)間較長(zhǎng)，成本昂貴，并不能滿足大規(guī)模生產(chǎn)VIP的要求.

目前導(dǎo)熱系數(shù)的檢測(cè)方法檢測(cè)效率相對(duì)低下，成本昂貴. 檢測(cè)速度與檢測(cè)成本之間的矛盾則是需要解決的另一重大問題. 為了解決以上問題，本文的目標(biāo)是研究一種快速、高精度的VIP導(dǎo)熱系數(shù)快速檢測(cè)方法. 該方法主要基于由埋入真空絕熱板內(nèi)傳感器因溫度變化引起的振蕩電路輸出信號(hào)頻率變化與導(dǎo)熱系數(shù)之間的關(guān)系. 在得到振蕩電路輸出信號(hào)頻率之后，如何與建立與導(dǎo)熱系數(shù)之間的關(guān)系變得格外重要.

極限學(xué)習(xí)機(jī)(ELM)是由文獻(xiàn)[7]提出來的一種針對(duì)SLFN(single-hidden layer feedforward neural net work)的新算法[7]. 該算法中，輸入權(quán)重和偏置可以被隨機(jī)賦值，且相應(yīng)輸出權(quán)重可直接計(jì)算. ELM已經(jīng)廣泛應(yīng)用于函數(shù)逼近、分類、識(shí)別等領(lǐng)域[8-10]. 但是ELM在新型模型中仍然存在擾動(dòng)和多重共線性問題，這將影響ELM的穩(wěn)定性和泛化能力. 嶺回歸是一種改良的最小二乘估計(jì)方法，用于有效解決在線性回歸分析中自變量存在共線性問題[11]. 本文的另一個(gè)目標(biāo)是基于嶺回歸改進(jìn)傳統(tǒng)ELM.

本文以平板熱流計(jì)法為基礎(chǔ)，提出一種新的埋入熱流計(jì)法檢測(cè)真空絕熱板導(dǎo)熱系數(shù)的方法. 該方法可以快速檢測(cè)VIP導(dǎo)熱系數(shù)并且精度較高、成本低廉. 本文的主要?jiǎng)?chuàng)新點(diǎn)在于提出一種新型的埋入式熱流計(jì)快速測(cè)量VIP導(dǎo)熱系數(shù)的方法，通過輸出信號(hào)頻率變化就可以反應(yīng)出VIP導(dǎo)熱系數(shù)的變化；通過嶺回歸對(duì)ELM的改進(jìn)提出了RRELM方法，來修正導(dǎo)熱系數(shù)與頻率變化之間的關(guān)系. 實(shí)驗(yàn)證明該方法不僅可以快速檢測(cè)VIP導(dǎo)熱系數(shù)，且與傳統(tǒng)ELM相比具有較高的精度.

1 檢測(cè)原理

VIP導(dǎo)熱系數(shù)隨板內(nèi)氣壓的變化而變化，文獻(xiàn)[12]展現(xiàn)了內(nèi)部不同材料的真空絕熱板的導(dǎo)熱系數(shù)與氣壓之間的關(guān)系，當(dāng)壓力大于10 Pa時(shí)，VIP的導(dǎo)熱系數(shù)隨壓力增加而增大. 本文提出的埋入熱流計(jì)法快速檢測(cè)真空絕熱板導(dǎo)熱系數(shù)方法的主要工作原理如下：

1)熱量傳遞包括導(dǎo)熱、對(duì)流和熱輻射3種基本方式. 真空絕熱板內(nèi)部氣壓很低，基本上不存在空氣流動(dòng). 電阻和VIP芯材緊密接觸，其在VIP內(nèi)部對(duì)流換熱不是主導(dǎo)換熱方式，因此VIP的主要換熱是輻射. 電阻與周圍輻射換熱計(jì)算公式為[13]

(1)

式中：T1為電阻溫度；T2為環(huán)境溫度；ε為電阻發(fā)射率，一般取值為0～1；A為電阻表面積；σ為Stefan-Boltzmann常數(shù).

2)將相同的加熱電阻和熱敏電阻構(gòu)建的傳感器件埋在不同真空絕熱板內(nèi)部的中心位置，當(dāng)環(huán)境溫度和濕度等外部因素相同時(shí)，因?yàn)閂IP內(nèi)部氣體壓力不同，其導(dǎo)熱系數(shù)也將不同. 因此，當(dāng)加熱電阻被加熱到溫度T后，加熱電阻在VIP內(nèi)部的散熱速度也會(huì)不同，其周圍的溫度熱場(chǎng)變化速度也不同. 將熱敏電阻接入振蕩電路，此時(shí)振蕩電路輸出端的頻率將隨熱敏電阻阻值的變化而產(chǎn)生改變. 真空絕熱板的導(dǎo)熱系數(shù)越大，溫度場(chǎng)的擴(kuò)散速度越快，振蕩電路輸出頻率的變化值也將越大，因此振蕩電路輸出頻率的變化就反映了VIP的導(dǎo)熱系數(shù).

2 ANSYS仿真分析

為了驗(yàn)證埋入熱流計(jì)法檢測(cè)真空絕熱板導(dǎo)熱系數(shù)方法的可行性，本文對(duì)內(nèi)置電阻散熱的變化過程進(jìn)行了ANSYS仿真分析.

仿真的環(huán)境溫度為25 ℃，電阻初始溫度為50 ℃，在不同內(nèi)部壓力下對(duì)VIP內(nèi)部的散熱過程進(jìn)行ANSYS仿真. 阻隔性鋁箔膜隔氣結(jié)構(gòu)為多層復(fù)合材料，本實(shí)驗(yàn)實(shí)際測(cè)量出的VIP阻隔性鋁箔膜厚度為0.1 mm，圖1為其模型的截面結(jié)構(gòu).

圖1 VIP內(nèi)部散熱實(shí)驗(yàn)?zāi)Ｐ徒孛?/p>

本文選取內(nèi)部壓力分別為0.5、50、100 Pa時(shí)，電阻和薄膜表面節(jié)點(diǎn)溫度隨時(shí)間變化的曲線，如圖2所示.

圖中TEMP2為電阻表面節(jié)點(diǎn)，TEMP3為薄膜表面節(jié)點(diǎn). 從仿真結(jié)果可以看出，在溫度環(huán)境等外部條件相同的情況下，內(nèi)部氣壓越大，電阻溫度變化也越大，因此可以說明本文提出的VIP導(dǎo)熱系數(shù)測(cè)量原理具有可行性.

(a)0.5 Pa

(b)50 Pa

(c)100 Pa

圖2 VIP不同內(nèi)部氣壓下電阻和薄膜表面節(jié)點(diǎn)溫度隨時(shí)間變化曲線

Fig.2 Temperature-time curve of resistance and membrane mode at different pressures

3 測(cè)量系統(tǒng)

根據(jù)提出的測(cè)量原理，VIP導(dǎo)熱系數(shù)測(cè)量系統(tǒng)，主要包括外部供電模塊、內(nèi)部測(cè)量模塊、外部測(cè)量模塊,如圖3所示.

外部供電模塊主要由無線供電發(fā)射電路組成，將能量傳送給內(nèi)部測(cè)量組件中的無線供電接收電路. 內(nèi)部測(cè)量模塊由加熱控制電路、溫度控制電路和RC振蕩電路組成. 振蕩電路如圖4所示，其中反相器采用74HC04D，R2為熱敏電阻. 當(dāng)內(nèi)部測(cè)量組件上電時(shí)，加熱控制電路中的加熱電阻便開始加熱；當(dāng)溫度控制電路檢測(cè)到加熱電阻達(dá)到預(yù)設(shè)的溫度值時(shí)，加熱控制電路將自動(dòng)控制加熱電阻停止加熱. RC振蕩電路由熱敏電阻、瓷片電容、電阻、蜂鳴片和非門組成. 在整個(gè)加熱和停止加熱過程中，熱敏電阻將一直監(jiān)測(cè)真空絕熱板內(nèi)的溫度場(chǎng)變化. 外部測(cè)量模塊由信號(hào)采集電路、信號(hào)放大濾波處理、微處理器等模塊組成，通過外部蜂鳴片來感應(yīng)內(nèi)部測(cè)量模塊中的蜂鳴片，從而獲得微弱的電信號(hào). 通過信號(hào)放大濾波等處理后，從而檢測(cè)到內(nèi)部測(cè)量組中輸出信號(hào)頻率的大小.

圖3 VIP導(dǎo)熱系數(shù)測(cè)量系統(tǒng)

圖4 RC振蕩電路

4 模型建立

由上述測(cè)量系統(tǒng)可以測(cè)得內(nèi)部壓力不同VIP的輸出信號(hào)頻率變化值，為了建立VIP導(dǎo)熱系數(shù)與輸出信號(hào)頻率變化值之間的關(guān)系，本文采用基于嶺回歸方法改進(jìn)了極限學(xué)習(xí)機(jī)，從而有效提高了系統(tǒng)的檢測(cè)精度.

4.1 極限學(xué)習(xí)機(jī)

ELM作為一種針對(duì)SLFN的新學(xué)習(xí)算法，其輸入權(quán)重和偏置可隨機(jī)賦值，并通過計(jì)算即可得到相應(yīng)的輸出權(quán)重. 其優(yōu)點(diǎn)是針對(duì)SLFN，相比傳統(tǒng)的學(xué)習(xí)算法，其學(xué)習(xí)速度較快. 對(duì)于一個(gè)單隱層神經(jīng)網(wǎng)絡(luò)，如果有N個(gè)任意樣本(xi,yi)，對(duì)于一個(gè)由M個(gè)隱含層節(jié)點(diǎn)組成的單隱層神經(jīng)網(wǎng)絡(luò)可以表示為

(2)

式中：g(x)為激活函數(shù)；Wi為輸入權(quán)重；βi為輸出權(quán)重；bi為第i個(gè)隱層單元的偏置.

單隱層神經(jīng)網(wǎng)絡(luò)的最終目的是使輸出誤差最小，計(jì)算公式為

(3)

即存在βi，Wi和bi使得

(4)

用矩陣表示為

Hβ=Y.

(5)

式中：H為隱含層節(jié)點(diǎn)的輸出；β為輸出權(quán)重；Y為期望輸出.

(6)

(7)

在ELM算法中，只要隨機(jī)確定了輸入權(quán)重和隱含層的偏置，隱含層的輸出矩陣就會(huì)被唯一確定. 輸出權(quán)重可以被確定為

(8)

式中H+為H的Moore-Penrose廣義逆.

(9)

4.2 基于嶺回歸的極限學(xué)習(xí)機(jī)

雖然ELM具有快速學(xué)習(xí)的能力和好的泛化能力，但是仍然需要進(jìn)一步改善. 已經(jīng)有很多研究者對(duì)ELM進(jìn)行了修改和改進(jìn)[14-16]. 雖然這些方法都改進(jìn)了ELM的性能，但仍然存在很大的發(fā)展空間，比如在非線性模型中的擾動(dòng)和多重共線性問題. 在多重共線性情況下，原有的模型采用最小二乘(OLS)估計(jì)法計(jì)算輸出權(quán)重，盡管最小二乘法估計(jì)對(duì)每個(gè)變量很公平，但差異較大，使得觀測(cè)值遠(yuǎn)離真實(shí)值. 而且在樣本量有限的情況下，最小二乘估計(jì)容易過擬合. 嶺回歸是一種用于存在多重共線性數(shù)據(jù)的有偏估計(jì)方法. 嶺回歸估計(jì)在評(píng)查誤差的基礎(chǔ)上計(jì)入了正則項(xiàng)，來降低標(biāo)準(zhǔn)誤差.

根據(jù)式(7)，嶺回歸估計(jì)可以定義[17]為

(10)

(11)

假設(shè)嶺參數(shù)是相同的，輸出權(quán)重公式為

(12)

從式(10)可以看出，嶺參數(shù)k對(duì)于嶺估計(jì)非常重要. 而且?guī)X參數(shù)k不是唯一的，嶺回歸估計(jì)實(shí)際上是回歸參數(shù)β的一個(gè)估計(jì)族.

本文采用3種不同的計(jì)算方法來計(jì)算k值，因此會(huì)得到3個(gè)不同的RRELM模型，并且將在進(jìn)行模型訓(xùn)練時(shí)選取最優(yōu)模型. 第1種選擇是文獻(xiàn)[17]提出的估計(jì)值:

(13)

(14)

第2種是文獻(xiàn)[11]提出的估計(jì)值:

(15)

第3種是文獻(xiàn)[18]提出的估計(jì)值:

(16)

5 實(shí) 驗(yàn)

5.1 頻率特征數(shù)據(jù)庫(kù)

本文選擇100塊導(dǎo)熱系數(shù)在1～15 mW/(m·K)之間的VIP板作為研究對(duì)象. 這些VIP板已由日本EKO公司生產(chǎn)的HC-074-300導(dǎo)熱系數(shù)測(cè)定儀進(jìn)行了測(cè)定. 由上述系統(tǒng)得到的VIP頻率特征值與HC-074-300導(dǎo)熱系數(shù)測(cè)儀得出的導(dǎo)熱系數(shù)構(gòu)建了VIP頻率特征數(shù)據(jù)庫(kù). 將數(shù)據(jù)庫(kù)中隨機(jī)選出30塊板作為測(cè)試集，剩余60塊板按照0.2的比例將數(shù)據(jù)庫(kù)中的樣本隨機(jī)分為訓(xùn)練集和驗(yàn)證集. 數(shù)據(jù)庫(kù)中頻率特征值是所建立VIP導(dǎo)熱系數(shù)測(cè)量系統(tǒng)對(duì)每一塊板進(jìn)行100次測(cè)量結(jié)果的均值.

5.2 測(cè)量結(jié)果與誤差分析

由式(13)～(16)可以得到3個(gè)不同的k值. 在3個(gè)不同的k值下，在驗(yàn)證集上尋找最優(yōu)模型，結(jié)果見表1. 本文采用平均絕對(duì)誤差(MAD)、均方跟誤差(RMSE)來進(jìn)行不同模型的對(duì)比.

表1 不同k值下RRELM模型的表現(xiàn)

以訓(xùn)練集作為輸入，建立頻率特征值與VIP導(dǎo)熱系數(shù)之間的RRELM非線性補(bǔ)償模型. RRELM模型包含3層，輸入層為振蕩電路輸出頻率，輸出層為VIP導(dǎo)熱系數(shù)，隱含層神經(jīng)元個(gè)數(shù)為10，k值為0.26. 為了驗(yàn)證RRELM非線性補(bǔ)償模型的精度，本文與傳統(tǒng)ELM非線性補(bǔ)償模型進(jìn)行對(duì)比. RRELM模型在測(cè)試集上的MAD為0.155，RMSE為0.024；ELM模型在測(cè)試集上的MAD為0.203，RMSE為0.214.

另外，為了進(jìn)一步驗(yàn)證RRELM非線性補(bǔ)償模型的精度，本文另外選取10塊VIP板進(jìn)行分析. 該10塊VIP板同樣由日本EKO公司生產(chǎn)的HC-074-300導(dǎo)熱系數(shù)測(cè)定儀進(jìn)行了測(cè)定. 與傳統(tǒng)ELM模型的測(cè)量結(jié)果進(jìn)行對(duì)比具有較高的精度，測(cè)量結(jié)果見表2.

表2 RRELM和傳統(tǒng)ELM的測(cè)量結(jié)果

從表2中可以看出，不同實(shí)驗(yàn)板測(cè)試結(jié)果的相對(duì)誤差波動(dòng)較大，這是因?yàn)椋?)真空絕熱板的導(dǎo)熱系數(shù)越大，其溫度場(chǎng)的擴(kuò)散速度也就越快，傳感器輸出頻率的變化值也越大，對(duì)于導(dǎo)熱系數(shù)較小的真空絕熱板，傳感器輸出頻率的變化值較小，由于傳感器本身精度和分辨率問題導(dǎo)致導(dǎo)熱系數(shù)較小的真空絕熱板相對(duì)誤差較大，今后將改進(jìn)硬件電路來提高傳感器的輸出精度和檢測(cè)分辨率；2)樣本數(shù)量有限，導(dǎo)致真空絕熱板的導(dǎo)熱系數(shù)與傳感器輸出頻率之間的關(guān)系沒有得到更好的擬合，使相對(duì)誤差波動(dòng)較大，在今后的研究中可以加大樣本數(shù)量.

另外，從表2中還可以看出，VIP導(dǎo)熱系數(shù)的測(cè)量值與實(shí)際值之間存在一定的誤差，但是提出RRELM非線性補(bǔ)償模型明顯優(yōu)于傳統(tǒng)ELM方法. 雖然RRELM非線性補(bǔ)償模型具有較高的測(cè)量精度，但是仍然存在一些誤差，這可能是由于以下原因造成的：1)樣本數(shù)量有限，建立的RRELM非線性補(bǔ)償模型仍然存在一定的非線性誤差; 2)測(cè)量系統(tǒng)的外部測(cè)量模塊中，外部電路采集到的輸出信號(hào)需要放大濾波處理，在實(shí)際測(cè)量中會(huì)產(chǎn)生一定的誤差.

6 結(jié) 論

1)提出了一種埋入流量計(jì)的快速測(cè)量VIP導(dǎo)熱系數(shù)的方法，在不同內(nèi)部壓力下，對(duì)VIP內(nèi)部的散熱過程進(jìn)行了ANSYS仿真，證明了埋入熱流量計(jì)法測(cè)量原理的可行性.

2)根據(jù)所提出的測(cè)量原理建立導(dǎo)熱系數(shù)測(cè)量系統(tǒng)，以測(cè)量?jī)?nèi)部壓力不同VIP的輸出信號(hào)頻率變化值.

3)基于嶺回歸方法建立了RRELM非線性補(bǔ)償模型，有效擬合了輸出頻率與導(dǎo)熱系數(shù)之間的關(guān)系，且與ELM相比取得了較高的精度.