[發(fā)明專利]一種利用傳遞函數(shù)修正量熱儀量熱過程延遲環(huán)節(jié)的方法

有效

• 申請?zhí)枺?/td>: 201810856810.5
• 申請日: 2018-07-31
• 公開（公告）號: CN109031942B
• 公開（公告）日: 2020-08-11
• 發(fā)明（設(shè)計）人: 張劍波;李哲;滕冠興;葛昊
• 申請（專利權(quán)）人: 清華大學(xué)
• 主分類號: G05B11/42
• 分類號: G05B11/42
• 代理公司: 北京清亦華知識產(chǎn)權(quán)代理事務(wù)所(普通合伙) 11201
• 代理人: 廖元秋

1.一種利用傳遞函數(shù)修正量熱儀量熱過程延遲環(huán)節(jié)的方法，其特征在于，該方法包括以下步驟：

1)任意選取一臺量熱儀和被測物體，將量熱儀量熱過程的延遲環(huán)節(jié)劃分為量熱儀量熱延遲環(huán)節(jié)和被測物體導(dǎo)熱延遲環(huán)節(jié)

2)獲取被測物體的離散儀器測量信號P_c(k)＝P_c(t_k)，其中k為非負(fù)整數(shù)，Δt為采樣間隔，且t_k＝kΔt；

3)獲取量熱儀量熱延遲環(huán)節(jié)傳遞函數(shù)G(s)；G(s)表達(dá)式如下：

其中，s為復(fù)數(shù)域的變量，P_c(s)為儀器測量信號P_c的拉普拉斯變換，N(s)為被測物體向量熱儀測量腔體散熱功率N的拉普拉斯變換；a_i、b_j為G(s)表達(dá)式中的系數(shù)；n為G(s)表達(dá)式中分母多項式的階數(shù)，m為G(s)表達(dá)式中分子多項式的階數(shù)；

4)利用量熱儀量熱延遲環(huán)節(jié)傳遞函數(shù)G(s)，采用以下三種方法中的任意一種計算被測物體向量熱儀測量腔體散熱功率的采樣序列N(k)：

方法一：根據(jù)量熱儀量熱過程延遲環(huán)節(jié)傳遞函數(shù)G(s)，有G(s)是真分式，取倒數(shù)后為假分式；假分式化為整式與真分式之和，即：

其中A_g為常系數(shù)，上式共有n-m+2項，每一項代表一個環(huán)節(jié)，相加代表這些環(huán)節(jié)是并聯(lián)的；

利用Z變換，將從s域變換到z域以用于離散信號的處理，即將P_c(k)從時域變換到z域，即P_c(z)＝Z(P_c(k))，根據(jù)計算得到被測物體向測量腔體散熱功率的Z變換N(z)；利用逆Z變換獲得被測物體向測量腔體散熱功率的采樣序列N(k)＝Z^-1(N(z))；

方法二：以儀器測量信號的采樣序列P_c(k)為輸入，利用時域方法經(jīng)過中并聯(lián)的環(huán)節(jié)，輸出被測物體向測量腔體散熱功率的采樣序列N(k)；

方法三：將z域與時域方法的結(jié)合，計算被測物體向測量腔體散熱功率的采樣序列N(k)；

5)對被測物體的導(dǎo)熱特征時間進行判定：若被測物體的導(dǎo)熱特征時間τ小于設(shè)定的導(dǎo)熱特征時間閾值，則忽略被測物體導(dǎo)熱延遲環(huán)節(jié)，本次測量的被測物體真實產(chǎn)熱值采樣序列Q_v(k)＝N(k)，量熱儀量熱過程延遲環(huán)節(jié)修正完畢；否則繼續(xù)修正被測物體導(dǎo)熱延遲環(huán)節(jié)，進入步驟6)；其中，被測物體的導(dǎo)熱特征時間計算方法如下：

若被測物體是圓柱形時，導(dǎo)熱特征時間τ計算方法如下：

其中ρ是被測物體密度，C_p是被測物體比熱容，D_sample和D_sample，0分別是被測物體的外徑和內(nèi)徑，H是被測物體的高度，k_r和k_z分別是被測物體的徑向?qū)嵯禂?shù)和軸向?qū)嵯禂?shù)；

若被測物體是方形時，導(dǎo)熱特征時間τ計算方法如下：

其中D、L和W分別是被測物體的厚度、長度和寬度，k_thr和k_in分別是被測物體的法向?qū)嵯禂?shù)和展向?qū)嵯禂?shù)；

6)獲取被測物體導(dǎo)熱延遲環(huán)節(jié)傳遞函數(shù)H(s)；H(s)的表達(dá)式如下：

其中，s為復(fù)數(shù)域的變量，Q_v(s)為被測物體產(chǎn)熱率的拉普拉斯變換，N(s)為被測物體向量熱儀測量腔體散熱功率的拉普拉斯變換；p_i、q_i為H(s)表達(dá)式中的系數(shù)，d為H(s)表達(dá)式中分母多項式的階數(shù)，e為H(s)表達(dá)式中分子多項式的階數(shù)；

7)利用被測物體導(dǎo)熱延遲環(huán)節(jié)傳遞函數(shù)H(s)，計算被測物體真實產(chǎn)熱值采樣序列Q_v(k)，量熱儀量熱過程延遲環(huán)節(jié)修正完畢；其中，計算被測物體真實產(chǎn)熱值采樣序列Q_v(k)，采用以下三種方法中的任意一種方法：

方法一：根據(jù)傳遞函數(shù)H(s)，有H(s)是真分式，分母多項式階數(shù)高于分子多項式階數(shù)，取倒數(shù)后為假分式，假分式化為整式與真分式之和，即：

其中B_h為常系數(shù)，上式共有e-d+2項，每一項代表一個環(huán)節(jié)，相加代表這些環(huán)節(jié)是并聯(lián)的；

利用Z變換，將從s域變換到z域以用于離散信號的處理，即將N(k)從時域變換到z域，即N(z)＝Z(N(k))，根據(jù)計算得到被測物體向測量腔體散熱功率的Z變換N(z)；利用逆Z變換獲得被測物體真實產(chǎn)熱值采樣序列Q_v(k)＝Z^-1(Q_v(z))；

方法二：將儀器測量信號的采樣序列N(k)為輸入，利用時域方法經(jīng)過中并聯(lián)的環(huán)節(jié)，輸出被測物體真實產(chǎn)熱值采樣序列Q_v(k)；

方法三：將z域與時域方法的結(jié)合，計算被測物體真實產(chǎn)熱值采樣序列Q_v(k)。