石友安, 魏東, 桂業(yè)偉, 錢煒祺, 曾磊, 劉驍

(1.中國空氣動(dòng)力研究與發(fā)展中心 空氣動(dòng)力學(xué)國家重點(diǎn)實(shí)驗(yàn)室， 四川 綿陽 621000;2.中國空氣動(dòng)力研究與發(fā)展中心 計(jì)算空氣動(dòng)力學(xué)研究所， 四川 綿陽 621000)

固體結(jié)構(gòu)內(nèi)部瞬態(tài)非均勻溫度場的重建方法研究

石友安1,2, 魏東2, 桂業(yè)偉2, 錢煒祺2, 曾磊2, 劉驍2

(1.中國空氣動(dòng)力研究與發(fā)展中心 空氣動(dòng)力學(xué)國家重點(diǎn)實(shí)驗(yàn)室， 四川 綿陽 621000;2.中國空氣動(dòng)力研究與發(fā)展中心 計(jì)算空氣動(dòng)力學(xué)研究所， 四川 綿陽 621000)

固體結(jié)構(gòu)內(nèi)部瞬態(tài)非均勻溫度場的無損測量在航空航天、機(jī)械制造、材料加工和醫(yī)療衛(wèi)生等領(lǐng)域都具有十分重要的作用?；诔暡}沖回波法，建立了超聲測量各向同性均勻介質(zhì)結(jié)構(gòu)內(nèi)部瞬態(tài)溫度分布的理論模型，發(fā)展了預(yù)測結(jié)構(gòu)內(nèi)部非均勻溫度場的靈敏度法和共軛梯度法，并系統(tǒng)比較分析了兩種瞬態(tài)溫度場重建算法的精度、抗噪性和穩(wěn)定性等特性。實(shí)驗(yàn)驗(yàn)證表明：基于熱傳導(dǎo)反問題的兩種方法，重建得到的結(jié)構(gòu)內(nèi)部瞬態(tài)非均勻溫度分布精度均較高，實(shí)時(shí)性好，適用性強(qiáng)；對結(jié)構(gòu)內(nèi)部瞬態(tài)溫度場的預(yù)測與評估以及探索研究結(jié)構(gòu)內(nèi)部量的新型測量和控制技術(shù)等具有參考意義。

兵器科學(xué)與技術(shù)； 瞬態(tài)溫度場； 超聲測溫； 溫度場重建； 靈敏度法； 共軛梯度法

0 引言

溫度作為一個(gè)基本物理量，是確定介質(zhì)狀態(tài)最重要的參數(shù)之一。結(jié)構(gòu)內(nèi)部瞬態(tài)溫度場的測量在航空航天、機(jī)械制造、材料加工和醫(yī)療衛(wèi)生等領(lǐng)域都具有十分重要的作用[1-2]。目前，結(jié)構(gòu)內(nèi)部的溫度測量主要采用內(nèi)埋傳感器(譬如熱電阻、熱電偶)的接觸式探測，但這種方法有一些局限性，如打孔安裝傳感器容易破壞結(jié)構(gòu)的原有形態(tài)，并可能導(dǎo)致結(jié)構(gòu)局部的溫度變化或產(chǎn)生應(yīng)力集中等問題。無損非接觸式測溫方法中的紅外技術(shù)、磷光熱圖等目前還僅適用于結(jié)構(gòu)表面溫度的測溫[3-4]。因此，發(fā)展應(yīng)用范圍更廣的無損溫度測量技術(shù)成為目前固體結(jié)構(gòu)內(nèi)部溫度測量與控制中亟待解決的熱點(diǎn)之一。

超聲波測溫技術(shù)可以實(shí)現(xiàn)對結(jié)構(gòu)內(nèi)部溫度場的無損非接觸式探測。相比傳統(tǒng)的測溫方法，超聲波測溫具有更快速、更精確、使用范圍更寬的優(yōu)勢[5-6]，且已經(jīng)在氣體、液態(tài)材料等方面取得一系列研究成果。但是，對于探測固體結(jié)構(gòu)內(nèi)部瞬態(tài)溫度場的研究還比較少，需要進(jìn)一步研究諸如納秒級聲時(shí)測量技術(shù)、瞬態(tài)溫度場重建的數(shù)理模型和多參數(shù)的反演計(jì)算方法等關(guān)鍵技術(shù)問題[7-8]。

本文從熱聲學(xué)理論和超聲測溫感知機(jī)理[9-10]出發(fā)，建立了結(jié)構(gòu)內(nèi)部瞬態(tài)非均勻溫度場非接觸探測的理論模型；從求解反問題角度入手，將瞬態(tài)非均勻溫度場的重建問題轉(zhuǎn)化為熱邊界的反演和熱傳導(dǎo)正問題的求解；應(yīng)用參數(shù)辨識中的靈敏度法和共軛梯度法，發(fā)展了基于求解熱傳導(dǎo)反問題的結(jié)構(gòu)內(nèi)部瞬態(tài)溫度場重建技術(shù)，并開展了一系列的數(shù)值仿真分析和有效性驗(yàn)證，得到了一些有意義的結(jié)果。

1 超聲測溫模型

如圖1所示，假設(shè)一均質(zhì)結(jié)構(gòu)左端受到沿空間均布的時(shí)變熱流加熱。結(jié)構(gòu)內(nèi)部的非均勻溫度場沿其受熱面垂直方向可視為一維分布狀態(tài)，即結(jié)構(gòu)內(nèi)部各點(diǎn)的溫度值僅是x的函數(shù)。超聲測溫中聲波探頭安裝在結(jié)構(gòu)右端。采用脈沖回波法進(jìn)行測量時(shí)，超聲波在固體介質(zhì)中的傳播時(shí)間ttof[11]可表示為

(1)

式中：T為溫度；L為超聲波在固體介質(zhì)中單向傳播的距離(被檢測結(jié)構(gòu)的厚度)；v是固體介質(zhì)中超聲波的傳播速度，與材料性能和結(jié)構(gòu)所受溫度相關(guān)。

圖1 超聲測溫的一維模型Fig.1 One-dimensional model of ultrasound speed measurement of temperature

在各向同性介質(zhì)中，v通?？梢员硎緸榻Y(jié)構(gòu)所受溫度的線性函數(shù)，如

v(T)=mT+n,

(2)

式中：m和n為常數(shù)，由實(shí)驗(yàn)標(biāo)定得到。

固體介質(zhì)中的溫度T，其一維情況下可以表示為位置x和時(shí)間t的函數(shù)，顯然T(x,t)即為待求的一維結(jié)構(gòu)內(nèi)部非均勻溫度場分布狀態(tài)。

圖2描述了鋼中超聲波縱波的傳播速度與溫度的變化關(guān)系[12]。

圖2 鋼中縱波波速與溫度的關(guān)系Fig.2 Longitudinal wave velocity vs, temperature in steel

因此，基于超聲波波速與溫度的相關(guān)性，首先采用超聲波脈沖回波法，準(zhǔn)確獲得超聲波在固體介質(zhì)中的傳播時(shí)間ttof. 然后，基于參數(shù)反演，根據(jù)(1)式以及聲波傳播速度與介質(zhì)溫度的函數(shù)關(guān)系反演出等效的邊界條件。最后，基于熱傳導(dǎo)的正問題求解獲得結(jié)構(gòu)內(nèi)部不同時(shí)刻的溫度場分布狀態(tài)。

2 瞬態(tài)溫度場的反演算法

基于傳播時(shí)間進(jìn)行溫度場的重建，本質(zhì)上是一個(gè)熱- 聲- 固體耦合的反問題。對于一維溫度場重建而言，可暫不考慮超聲波傳播路徑的變化。

一維無內(nèi)熱源熱傳導(dǎo)正問題可描述為

主控方程

(3)

T|x=L=Tcons,

初值條件T|t=0=Tcons,

觀測方程

ttof,m=ttof,ex+εtof,m,

(4)

式中：Ttop為待反演的溫度邊界；Tcons為恒定溫度值；k為導(dǎo)熱系數(shù)；ρ為密度；cp為比熱容；q為熱流密度；ε為測量誤差；t為時(shí)間；下標(biāo)“m”表示測量得到，下標(biāo)“ex”表示精確值。反問題描述為：已知測量的超聲波傳播時(shí)間ttof(見(4)式)，求解主控方程中的邊界熱流q(t)或Ttop.

反演算法思路：以傳播時(shí)間為依據(jù)，在數(shù)值求解正問題的基礎(chǔ)上，根據(jù)輸出誤差原則，將反問題轉(zhuǎn)化為一個(gè)優(yōu)化問題，等價(jià)于尋求合適的函數(shù)q(t)(或者Ttop)，使目標(biāo)泛函J達(dá)極小值的過程：

(5)

式中：下標(biāo)“c”表示計(jì)算得到。邊值條件q(t)或Ttop的反演估計(jì)中，目標(biāo)函數(shù)J(q(t))設(shè)定后，則物理系統(tǒng)參數(shù)反分析問題就可以歸結(jié)為求解如下非線性優(yōu)化問題：

minJ(q(t))，

T(x,t),t>0,x∈[0,L].

(6)

2.1 靈敏度法

靈敏度法[13]的描述如下：目標(biāo)函數(shù)達(dá)到極小的必要條件為?J(θ)/?θ=0，θ是待辨識參數(shù)，當(dāng)在第i迭代層次上Ji(θi)不是極小值時(shí)，必須調(diào)整θi使Ji+1達(dá)到極小值，則有

(7)

略去2階以上小量后得

(8)

式中：M稱為信息矩陣。向量?J/?θ和矩陣M中都含有狀態(tài)變量對參數(shù)的導(dǎo)數(shù)，該導(dǎo)數(shù)稱為靈敏度，由(8)式構(gòu)成的算法稱之為靈敏度法。對于穩(wěn)態(tài)邊界，信息矩陣則退化為單一元素；對于瞬態(tài)邊界，則為n×n矩陣。由于熱傳導(dǎo)方程的擴(kuò)散性，n×n矩陣多為病態(tài)矩陣[14]。

因此，在本文中，靈敏度法主要用于穩(wěn)態(tài)邊界Ttop的反演。根據(jù)變分原理，對目標(biāo)函數(shù)取變分，得到目標(biāo)函數(shù)關(guān)于待辨識邊界的導(dǎo)數(shù)為

(9)

對(9)式再求導(dǎo)，則得到信息矩陣為

(10)

式中：?T/?Ttop為溫度場關(guān)于待辨識加熱邊界的靈敏度。將主控方程(3)式關(guān)于待辨識邊界求導(dǎo)即可得到靈敏度方程。優(yōu)化步驟見(8)式。

2.2 共軛梯度法

共軛梯度法也稱為迭代正則化方法[15-16]，分解為熱傳導(dǎo)正問題、靈敏度問題和伴隨變量問題進(jìn)行求解。優(yōu)化算法[15]描述如下：

(11)

為了得到伴隨方程，引入伴隨變量λ，將目標(biāo)函數(shù)(5)式寫為如下的擴(kuò)展形式：

(12)

對(12)式右邊第二部分做分部積分后再做變分，得到伴隨方程為

(13)

(14)

步長由(15)式計(jì)算：

(15)

式中：Δttof是Δq=Pn引起的傳播時(shí)間的變化值。在優(yōu)化過程中，收斂準(zhǔn)則根據(jù)輸出誤差原則獲得

J<δ，δ=Mmσ2,

(16)

式中：Mm為測量點(diǎn)數(shù)。至此，梯度、步長均已求出，共軛梯度可由梯度、步長構(gòu)造。具體優(yōu)化過程可參閱文獻(xiàn)[17]。

3 精度、抗噪性和穩(wěn)定性分析

本節(jié)通過數(shù)值仿真，分析兩種反演算法的精度、抗噪性以及穩(wěn)定性等。

數(shù)值仿真中，假定超聲波波型為縱波，介質(zhì)材料為鋼，長度5 cm. 鋼中縱波波速與溫度的關(guān)系見(17)式。通過計(jì)算相應(yīng)的正問題獲得超聲波傳播時(shí)間ttof的真值。

v(T)=-0.648 0T+5 9280.0.

(17)

首先，針對靈敏度法開展仿真分析。假設(shè)加熱面溫度為350 K，計(jì)算時(shí)間為60 s. 為了更好地模擬實(shí)際測量數(shù)據(jù)，在精確時(shí)間的基礎(chǔ)上，分別疊加標(biāo)準(zhǔn)差為0，均方差σ為1×10-7s、5×10-7s、1×10-6s和5×10-6s的白噪聲作為模擬傳播時(shí)間。表1給出了超聲波傳播時(shí)間疊加不同噪聲情況下加熱邊界的辨識結(jié)果對比。

圖3給出了距加熱面15 mm測點(diǎn)處辨識的溫升對比。從圖3中可以發(fā)現(xiàn)，隨著信噪比的降低(噪聲的增大)，辨識精度下降，相對溫升的誤差增長較快。例如無白噪聲時(shí)，溫升相對誤差(即加熱面溫度辨識值與實(shí)際值的比值)為0.000 3%；疊加5×10-7s的白噪聲后，相對誤差約為1.67%. 由此可以看出，靈敏度算法具有較好的抗噪性。

表1 靈敏度法辨識的溫度邊界結(jié)果對比Tab.1 Comparison of boundary temperatures estimated by sensitivity method

圖3 采用靈敏度法反演的測點(diǎn)處溫升歷程對比Fig.3 Temperature variation at measurement point estimated by sensitivity method

接著，通過典型熱流的反演，分析共軛梯度法的精度和抗噪性等。

邊界熱流為正弦和階躍熱流，計(jì)算時(shí)間為5 s. 通過計(jì)算相應(yīng)的正問題，得到模擬傳播時(shí)間。在實(shí)際應(yīng)用中，誤差不可避免。因此，在模擬值的基礎(chǔ)上分別疊加標(biāo)準(zhǔn)差為1×10-8s、1×10-9s和1×10-10s的白噪聲作為實(shí)際測量值，具體見圖4. 從圖4中可以發(fā)現(xiàn)，超聲波的傳播時(shí)間非常短暫，為微秒量級；相較于傳播時(shí)間總的變化值(約0.11×10-6s)，1×10-8s的誤差標(biāo)準(zhǔn)差相當(dāng)于帶來了近10%的敏感度誤差。

圖4 模擬的傳播時(shí)間Fig.4 Simulated transit time of ultrasonic wave

圖5顯示基于超聲傳播時(shí)間，采用共軛梯度法反演得到的熱流值對比情況。從圖5中可以發(fā)現(xiàn)，在非連續(xù)點(diǎn)處，辨識值存在一定的“數(shù)值振蕩”，如階躍熱流的階躍點(diǎn)；此外，在起始點(diǎn)和尾點(diǎn)處與模擬真值還存在一定的差異。除此之外，隨著誤差的減小，辨識值與模擬真值的吻合度增加；當(dāng)誤差標(biāo)準(zhǔn)差為1×10-8s時(shí)，兩種熱流的反演值均較差，這主要與信號的信噪比過小有關(guān)；當(dāng)誤差標(biāo)準(zhǔn)差為1×10-9s時(shí)，正弦熱流反演值的偏差帶約為±5%，階躍熱流反演值的偏差帶約為±10%；當(dāng)誤差標(biāo)準(zhǔn)差為1×10-10s時(shí)，正弦熱流反演值的偏差帶約為±3%，階躍熱流反演值的偏差帶約為±5%. 起始點(diǎn)處、尾點(diǎn)處與真值的差異，主要是與起始點(diǎn)處靈敏度較小，尾點(diǎn)處伴隨變量始終為0有關(guān)(可參見第2節(jié)中的(13)式)。起始點(diǎn)處的靈敏度較小，是由熱傳導(dǎo)過程的擴(kuò)散特性決定的，修正較為困難；尾點(diǎn)處的辨識值較差可以通過修正伴隨變量而改善，如λtmax=λtmax-1.

圖5 采用共軛梯度法反演的熱流值對比Fig.5 Comparison of heat fluxes estimated by conjugate gradient method

圖6顯示了基于反演的熱流重建得到的典型時(shí)刻溫度分布的對比情況。從圖6中可以發(fā)現(xiàn)：最大偏差出現(xiàn)在熱邊界處，除此之外，計(jì)算值與模擬真值吻合較好；隨著誤差的增加，偏差增大；5 s時(shí)刻，當(dāng)白噪聲的標(biāo)準(zhǔn)差為1×10-8s時(shí)，正弦熱流邊界的溫度最大偏差55 ℃(起始點(diǎn)處)，約10%，而階躍熱流邊界的最大偏差為75 ℃，約12%；1×10-9s時(shí)，正弦熱流的最大偏差降為45 ℃，約8.8%，階躍熱流則降為30 ℃，約6%；1×10-10s時(shí)，正弦熱流的溫度最大偏差急劇減小至10 ℃，約2%，階躍熱流則降為5 ℃，約1%. 最大偏差出現(xiàn)在起始點(diǎn)處，主要是與該點(diǎn)的靈敏度較小有關(guān)。

圖6 重建的不同時(shí)刻溫度分布對比Fig.6 Comparison of reconstructed temperature distributions at different time

通過上述分析，不難發(fā)現(xiàn)，為了保持重建的溫度場具有較好的精度，傳播時(shí)間的精度需要保持1×10-9s量級。結(jié)合(17)式可知，縱波在鋼中傳播時(shí)，溫度變化1 ℃，波速變化0.648 m/s. 傳播時(shí)間與傳播路徑長度和波速有關(guān)。在本算例中，材料長度為5 cm，則材料的溫度變化1 ℃，傳播時(shí)間變化約1.8×10-9s. 由此可見，超聲波對于結(jié)構(gòu)內(nèi)部溫度的敏感度非常微弱。這對測試信號的信噪比和重建算法的抗噪性提出要求。

最后，分析一下算法的穩(wěn)定性。熱傳導(dǎo)反問題是一個(gè)不適定的問題，體現(xiàn)在解不連續(xù)依賴于數(shù)據(jù)，即觀測數(shù)據(jù)的細(xì)微誤差可能導(dǎo)致巨大變化。要克服這種不適定性，通常的做法是對目標(biāo)函數(shù)正則化，使不適定性問題變?yōu)檫m定。

正則化主要有兩種方法，一種是Beck[10]的正則化方法，另一種是文獻(xiàn)[18]中Alifanov的迭代正則化方法，通過設(shè)定合適的收斂準(zhǔn)則來達(dá)到正則化的目的，例如，本文的(16)式。該方法實(shí)現(xiàn)方便，適于工程應(yīng)用。

定義辨識值與真值之間的相對偏差E為

(18)

式中：下標(biāo)“av”表示平均值。

定義計(jì)算傳播時(shí)間與測量值之間的平均偏差S為

(19)

從圖7可以看出：在反演過程中，E的極小值并非出現(xiàn)在J的極值處，而是某一局部值處，約Mmσ2，即(16)式。E與J關(guān)系實(shí)質(zhì)上是輸出誤差原則的體現(xiàn)，只要合理的利用這種關(guān)系就可以使算法具有良好的抗噪性與穩(wěn)定性。

圖7 E隨目標(biāo)函數(shù)J下降的變化關(guān)系Fig.7 Variation of E with objective function J

綜上可知，基于靈敏度法和共軛梯度法建立的固體結(jié)構(gòu)內(nèi)部瞬態(tài)溫度場的反演方法，精度較高、抗噪性較好，并且具有可靠的穩(wěn)定性。

4 算法驗(yàn)證

為了進(jìn)一步分析算法的有效性，本節(jié)借鑒文獻(xiàn)[19]中的實(shí)驗(yàn)數(shù)據(jù)，進(jìn)行兩種算法的有效性驗(yàn)證。實(shí)驗(yàn)示意如圖8所示。

圖8 文獻(xiàn)[19]實(shí)驗(yàn)示意圖Fig.8 Schematic diagram of experimental setup for ultrasonic temperature monitoring in Ref. [19]

測試材料為不銹鋼，長度為30 mm. 底部采用75 ℃熱水加熱10 s，其他邊界可以近似為絕熱。為了對比分析超聲測溫?cái)?shù)據(jù)，在試件內(nèi)部的長度方向，均布5個(gè)熱電偶。初始溫度為26 ℃. 試件中超聲波的波速與溫度的關(guān)系式為

v(T)=-0.648T+5 934.9.

(20)

圖9給出了超聲波傳播時(shí)間的計(jì)算值與文獻(xiàn)[19]實(shí)驗(yàn)測量值的對比。從圖9中可以發(fā)現(xiàn)：對于共軛梯度法而言，傳播時(shí)間計(jì)算值的變化趨勢與測量值吻合非常好，最大偏差小于0.1%；對于靈敏度算法而言，6 s前，計(jì)算得到的傳播時(shí)間偏大，偏差隨時(shí)間增加而減??；6 s后，計(jì)算值偏小，偏差隨時(shí)間增加而增大；全時(shí)間域內(nèi)平均偏差小于0.5%. 分析認(rèn)為，在共軛梯度法中，反演的熱邊界為瞬態(tài)，更符合實(shí)驗(yàn)中加熱邊界的物理描述，而在靈敏度算法中，反演的熱邊界為穩(wěn)態(tài)，若物理邊界為瞬態(tài)，則是一種等效近似處理。

圖9 渡越時(shí)間計(jì)算值與文獻(xiàn)[19]對比Fig.9 Comparison of estimated transit time of ultrasonic waveand that in Ref.[19]

圖10給出了不同時(shí)刻溫度分布的重建值與文獻(xiàn)[19]測量值的對比。從圖10中可以發(fā)現(xiàn)，采用共軛梯度法重建的溫度分布與熱電偶的測量值吻合得最好；對于靈敏度算法而言， 5 s時(shí)刻的重建精度高于文獻(xiàn)[19]中的重建算法，接近共軛梯度法，但是1 s、和2.6 s的精度則最差；距離加熱邊界越近，重建誤差越大，在熱電偶1點(diǎn)處，共軛梯度的偏差小于0.1 ℃，靈敏度法約1 ℃，而文獻(xiàn)[19]中的算法則為2.2 ℃. 由此可見，共軛梯度法具有較高精度，且適用性較好，而靈敏度算法受限于穩(wěn)態(tài)邊界限制，適用范圍有限。

圖10 重建的不同時(shí)刻溫度分布值與文獻(xiàn)[19]測量值對比Fig.10 Comparison of reconstructed temperature distribution at different time and those in Ref.[19]

圖11給出了不同測點(diǎn)的溫升歷程與測量值的對比。從圖11中可以發(fā)現(xiàn)：5 s時(shí)刻，采用共軛梯度法重建得到的溫升歷程與測量值吻合最好，最大偏差小于0.1 ℃，靈敏度法與文獻(xiàn)[19]中的重建算法精度相當(dāng)，最大偏差約2.2 ℃；對于1 s和2.6 s時(shí)刻，采用靈敏度法重建得到的溫升歷程與測量值吻合得最好，最大偏差小于0.1 ℃，共軛梯度法與文獻(xiàn)[19]中的重建算法精度相當(dāng)，最大偏差約0.3 ℃. 這主要是因?yàn)? s和2.6 s時(shí)刻，加熱時(shí)間太短，熱邊界可以近似為穩(wěn)態(tài)，而共軛梯度法由于初始時(shí)刻靈敏度較小，重建精度略受影響。

圖11 不同測點(diǎn)處溫升歷程與文獻(xiàn)[19]測量值對比Fig.11 Comparison of estimated temperature histories at different measuring points and those in Ref.[19]

為了增強(qiáng)靈敏度算法的適用性，本文針對靈敏度算法，采用準(zhǔn)定常近似處理，分段重建結(jié)構(gòu)內(nèi)部瞬態(tài)溫度分布。圖12給出了采用準(zhǔn)定常近似后，傳播時(shí)間計(jì)算值與文獻(xiàn)[19]測量值的對比。從圖12可以發(fā)現(xiàn)，采用準(zhǔn)定常近似后，傳播時(shí)間隨時(shí)間的變化趨勢與測量值的吻合度增加。

圖12 準(zhǔn)定常修正后傳播時(shí)間計(jì)算值與文獻(xiàn)[19]測量值對比Fig.12 Comparison of quasi-steady corrected estimated transit time of ultrasonic wave and that in Ref.[19]

圖13給出了采用準(zhǔn)定常近似后不同測點(diǎn)的溫升歷程與測量值的對比。從圖13中可以發(fā)現(xiàn)，采用準(zhǔn)定常近似后，采用靈敏度算法重建的溫度分布的精度得到了有效的改善，其精度接近于共軛梯度法。圖14給出了共軛梯度法反演得到的熱流。從圖14可以清晰地分辨出熱邊界幾乎保持為一個(gè)恒定值，符合物理實(shí)際。起點(diǎn)處和尾點(diǎn)處的熱流值偏低，這主要與起始點(diǎn)處靈敏度較小，尾點(diǎn)處伴隨變量始終為0有關(guān)。這兩個(gè)點(diǎn)的反演值受初值影響較大。

圖13 準(zhǔn)定常近似后不同時(shí)刻溫度分布值與文獻(xiàn)[19]測量值對比Fig.13 Comparison of quasi-steady corrected temperature distributions at different time and those in Ref.[19]

圖14 共軛梯度法反演得到的熱流Fig.14 Heat flux estimated by CGM

總體來看，基于超聲傳播時(shí)間，建立的兩種結(jié)構(gòu)內(nèi)部瞬態(tài)溫度場重建方法有效，反演的熱邊界條件符合物理實(shí)際，重建的瞬態(tài)溫度分布精度較高，實(shí)時(shí)性好，將有利于促進(jìn)超聲無損測溫技術(shù)的發(fā)展。

5 結(jié)論

1)本文基于反演穩(wěn)態(tài)、瞬態(tài)邊界重建的結(jié)構(gòu)內(nèi)部瞬態(tài)溫度分布，除個(gè)別點(diǎn)外，與真值吻合良好，算法精度高、抗噪性好、穩(wěn)定性強(qiáng)。

2)比較兩種算法，靈敏度法主要適用于穩(wěn)態(tài)邊界的反演，推導(dǎo)和實(shí)現(xiàn)較為方便，精度和計(jì)算效率較好，但是對于瞬態(tài)邊界反演而言，則受限于信息矩陣的病態(tài)特性，效果欠佳，若采用準(zhǔn)定常近似，可以有效改善；共軛梯度法能很好地適用于瞬態(tài)邊界的反演，精度較高，但是推導(dǎo)和實(shí)現(xiàn)較為復(fù)雜，起點(diǎn)和尾點(diǎn)的反演值需修正。

3)對比文獻(xiàn)[19]實(shí)驗(yàn)結(jié)果可知：基于超聲傳播時(shí)間反演的熱邊界條件，符合物理實(shí)際，重建的瞬態(tài)溫度分布精度較高，實(shí)時(shí)性好，適用性強(qiáng)，將有利于促進(jìn)超聲無損測量固體結(jié)構(gòu)內(nèi)部溫度分布技術(shù)的發(fā)展，具有一定的工程應(yīng)用價(jià)值。

References)

[1] Lu C H, Tsai C C. Adaptive decoupling predictive temperature control for an extrusion barrel in a plastic injection molding process[J]. IEEE Transactions on Industrial Electronics, 2001, 48(5):968-975.

[2] Brown E C, Olley P, Coates P D. In line melt temperature measurement during real time ultrasound monitoring of single screw extrusion[J]. Plastics Rubber & Composites, 1999, 29(1):3-13.

[3] Smith T B, Schetz J A. Development and ground testing of heat flux gages for high enthalpy supersonic flight tests[C]∥22nd Aerodynamic Measurement Technology and Ground Testing Conference. St Louis, Missouri ,US: AIAA, 2002.

[4] Gartenberg E, Roberts A S. Twenty-five years of aerodynamic research with infrared imaging[J]. Journal of Aircraft, 1992, 29(2): 161-171.

[5] Wang H Z, Yu H B, Lin F L. Method and experiment of noninvasive temperature estimation by ultrasound echo pulses[C]∥Proceedings of the 20th Annual International Conference of the IEEE Engineering in Medicine and Biology Society. Hong Kong, China: IEEE,1998: 1513-1516.

[6] Koo K M, Jeong D G, Chio J H , et al. A new measurement system of very high temperature in atomic pile using ultrasonic delay time[C]∥IEEE Region 10 International Conference on Electrical and Electronic Technology. Singapore: IEEE,2001.

[7] 孫崇正. 超聲波測溫技術(shù)進(jìn)展[J]. 宇航計(jì)測技術(shù), 1995, 15(2):34-41. SUN Chong-zheng. Progress of ultrasonic thermometry[J]. Astronautic Metrology and Measurement, 1995, 15(2):34-41. (in Chinese)

[8] 常蕾，趙儉. 超聲波測溫技術(shù)在高溫氣流溫度場測量中的應(yīng)用[J]. 計(jì)測技術(shù),2014, 34(1):1-9. CHANG Lei, ZHAO Jian. Application of ultrasonic thermometry in measuring temperature field of high temperature gas flow[J]. Metrology & Measurement Technology, 2014, 34(1):1-9.(in Chinese)

[9] 魏培君, 章梓茂. 彈性動(dòng)力學(xué)反問題的數(shù)值反演方法[J]. 力學(xué)進(jìn)展, 2001,31(2):172-180. WEi Pei-jun, ZHANG Zi-mao. Numerical methods for inverse problems in elastic dynamic[J]. Advances in Mechanics, 2001,31(2):172-180.(in Chinese)

[10] Beck J V. Nonlinear estimation applied to the nonlinear inverse heat conduction problem[J]. International Journal of Heat and Mass Transfer, 1970, 13(4):703-716.

[11] Ihara I, Takahashi M. Ultrasound thermometry for monitoring internal temperature gradient in heated material[C]∥IEEE International Ultrasonic Symposium. Rome: IEEE, 2009:1199-1202.

[12] Ihara I, Takahashi M. A new method for internal temperature profile measurement by ultrasound[C]∥Instrumentation and Measurement Technology Conference. Warsaw, Poland: IEEE, 2007.

[13] 錢煒祺, 蔡金獅. 再入航天飛機(jī)表面熱流密度辨識[J]. 宇航學(xué)報(bào), 2000, 21(4): 1-6. QIAN Wei-qi, CAI Jin-shi. Surface heat flux identification of reentry space shuttle[J]. Journal of Astronautics, 2000, 21(4): 1-6.(in Chinese)

[14] 蔡金獅, 汪清, 王文正. 飛行器系統(tǒng)辨識學(xué)[M]. 北京：國防工業(yè)出版社, 2003. CAI Jin-shi, WANG Qing, WANG Wen-zheng. The aircraft system identification[M]. Beijing: National Defense Industry Press, 2003.(in Chinese)

[15] 錢煒祺, 何開鋒, 汪清. 三維非穩(wěn)態(tài)熱傳導(dǎo)逆問題反演算法研究[J]. 力學(xué)學(xué)報(bào), 2008, 40(5):611-618. QIAN Wei-qi, HE Kai-feng, WANG Qing. Inverse estimation of heat source term in three dimensional transient heat conduction problems [J]. Chinese Journal of Theoretical and Applied Mechanics, 2008, 40(5):611-618.(in Chinese)

[16] Shi Y A, Zeng L, Qian, W Q, et al. A data processing method in the experiment of heat flux testing using inverse methods[J]. Aerospace Science and Technology, 2013, 29(1): 74-80.

[17] 解可新, 韓健, 林友聯(lián). 最優(yōu)化方法[M]. 天津: 天津大學(xué)出版社, 1997. XIE Ke-xin, HAN Jian, LIN You-lian. The optimization methods[M]. Tianjin: Tianjin University Press,1993.(in Chinese)

[18] Alifanov O M. Inverse heat transfer problems[M]. Berlin: Springer-Verlag, 1994.

[19] Takahashi M, Ihara I. Ultrasonic sensing of internal temperature distributions in heated materials[C]∥ International Joint Conference. Busan, Korea: SICE-ICASE,2006:1071-1074.

Reconstruction of Transient Nonuniform Temperature Field in Solid Structures Using Inverse Methods

SHI You-an1,2, WEI Dong2, GUI Ye-wei2, QIAN Wei-qi2, ZENG Lei2, LIU Xiao2

(1.State Key Laboratory of Aerodynamics,China Aerodynamics Research and Development Center, Mianyang 621000, Sichuan， China;2.Computational Aerodynamics Institute, China Aerodynamics Research and Development Center, Mianyang 621000,Sichuan，China)

A new method is presented for reconstruction of internal transient nonuniform temperature field in solid structures . The proposed method consists of a joint of heat conduction calculation and an inverse analysis for boundary parameter estimation using ultrasonic time of flight by conjugate gradient method or sensitivity method. The principle for internal temperature sensing is based on temperature dependence of the velocity of ultrasonic wave that propagates through solid material. Serial simulation and stability analysis are made， in which classic function form of heat flux or constant temperature is estimated to illustrate the calculating accuracy and the factors of influencing two parameter estimation methods. Comparison experiments with a steel plate which single side is heated by hot water are referred. The temperature gradient and its variation estimated by the two methods agree well with the data measured using thermocouples installed in a steel structure. The results show that the proposed method is a promising means for high accurate reconstruction of internal transient temperature field in solid structure.

ordnance science and technology; transient temperature field； ultrasound speed measurement of temperature； reconstruction of temperature field； sensitivity method； conjugate gradient method

2016-03-17

國家自然科學(xué)基金項(xiàng)目(11402285、11372338)； 國家“973”計(jì)劃項(xiàng)目(2015CB755800)

石友安(1985—)，男，助理研究員。E-mail：youanshi@sina.com

TB942

A

1000-1093(2016)12-2347-09

10.3969/j.issn.1000-1093.2016.12.022