德國(guó)進(jìn)口eps全系列EDP Mini 200E/50E德國(guó)進(jìn)口eps全系列EDP Mini 200E/50E

 　【儀表網(wǎng) 儀表深度】時(shí)下，美股遭遇重創(chuàng)，道指狂跌，美股年內(nèi)大單日跌幅屢屢讓人大跌眼鏡，這一系列的信號(hào)引發(fā)了人們對(duì)經(jīng)濟(jì)衰退的擔(dān)憂(yōu)。上一次出現(xiàn)這樣的情況，還是2007年大衰退前夕。

 
　　目前，經(jīng)濟(jì)下行已經(jīng)已成為一個(gè)不爭(zhēng)的事實(shí)，而會(huì)否引發(fā)又一次的金融危機(jī)成為人們普遍擔(dān)憂(yōu)的事情。不管金融危機(jī)會(huì)不會(huì)爆發(fā)，何時(shí)爆發(fā)，可以肯定的是當(dāng)前經(jīng)濟(jì)增長(zhǎng)疲軟，在這種背景下，制造業(yè)將會(huì)面臨越來(lái)越嚴(yán)峻的挑戰(zhàn)，我們儀器儀表產(chǎn)業(yè)也將無(wú)可避免，儀器儀表企業(yè)需要做好迎接挑戰(zhàn)的準(zhǔn)備。
 
　　儀器儀表企業(yè)面臨“前狼后虎”之勢(shì)
 
　　近年來(lái)，我國(guó)儀器儀表產(chǎn)業(yè)發(fā)展迅速，中低檔產(chǎn)品已具有一定的市場(chǎng)競(jìng)爭(zhēng)力，很多企業(yè)為了避免國(guó)內(nèi)的競(jìng)爭(zhēng)和更多創(chuàng)收，選擇“走出去”。在這些企業(yè)中，有的成功了，有的卻是鎩羽而歸?？偟氖莵?lái)說(shuō)，絕大多數(shù)是成功了。但是，面對(duì)經(jīng)濟(jì)發(fā)展下行，儀器儀表企業(yè)還能如過(guò)去一樣如魚(yú)得水嗎？
 
　　顯然不可能。貿(mào)易戰(zhàn)對(duì)儀器儀表企業(yè)的實(shí)際影響有多大，不同的企業(yè)影響也不一樣，各家也是如人飲水，冷暖自知。但可以確定的，“貿(mào)易戰(zhàn)”的風(fēng)波，主要從事對(duì)外出口的儀器儀表企業(yè)不得不做出相應(yīng)的調(diào)整。從中美經(jīng)濟(jì)大局看，貿(mào)易戰(zhàn)開(kāi)始以來(lái)，雙方的進(jìn)出口商品急劇下降，在整體“潮退”之下，儀器儀表企業(yè)想要獨(dú)善其身，顯然不太可能。
 
　　中美貿(mào)易戰(zhàn)的持續(xù)惡化，無(wú)疑會(huì)對(duì)中國(guó)制造型企業(yè)產(chǎn)生極大的影響，儀器儀表今年上半年的出口不夠理想，接下來(lái)可能迎來(lái)更加艱難的局面。
 
　　雖然儀器儀表企業(yè)也可以選擇海外其他國(guó)家地區(qū)的市場(chǎng)，但是，在經(jīng)濟(jì)低迷之際，其他市場(chǎng)只怕也難以帶來(lái)高速的增長(zhǎng)，儀器儀表須做好壞的準(zhǔn)備。因?yàn)樵谝粋€(gè)整體上增速放緩的市場(chǎng)環(huán)境下，儀器儀表必然迎來(lái)更加激烈的競(jìng)爭(zhēng)。
 
　　此外，儀器儀表除了受宏觀(guān)經(jīng)濟(jì)的影響，同時(shí)，還面臨著人工智能、物聯(lián)網(wǎng)等技術(shù)的沖擊。這一點(diǎn)，是各行各業(yè)都會(huì)面臨的問(wèn)題，在一些展會(huì)上已經(jīng)可以看出端倪，很多企業(yè)推出了和人工智能結(jié)合的新型產(chǎn)品。所以，在這種“前有狼后有虎”的局勢(shì)下，企業(yè)將要面臨的壓力也就可想而知了，企業(yè)要想走出一條適合自己的線(xiàn)，是真的需要下大功夫了。
 
　　儀器儀表企業(yè)如何突圍？
 
　　當(dāng)前中國(guó)有差不多14億的人口，可以說(shuō)是比較大的儀器儀表生產(chǎn)制造及應(yīng)用市場(chǎng)，應(yīng)對(duì)貿(mào)易戰(zhàn)，國(guó)內(nèi)市場(chǎng)的需求，對(duì)中國(guó)企業(yè)來(lái)說(shuō)也是巨大的機(jī)遇，內(nèi)需將會(huì)成為經(jīng)濟(jì)增長(zhǎng)的主要?jiǎng)恿?。同時(shí)，隨著智慧城市、物聯(lián)網(wǎng)的發(fā)展，也會(huì)為儀器儀表企業(yè)注入新的發(fā)展機(jī)遇
 
　　我國(guó)的儀器儀表產(chǎn)業(yè)中低檔的產(chǎn)品已經(jīng)具備競(jìng)爭(zhēng)力，只是高檔的科研儀器產(chǎn)品仍依賴(lài)進(jìn)口，這是我國(guó)儀器儀表行業(yè)所面臨的瓶頸，也是政府和企業(yè)應(yīng)該重視和突破的問(wèn)題。
 
　　另外，當(dāng)前絕大多數(shù)儀器儀表企業(yè)仍局限于提供單一的產(chǎn)品，不能夠提供完整的行業(yè)解決方案，認(rèn)為產(chǎn)品過(guò)硬，就不怕沒(méi)有市場(chǎng)。產(chǎn)品質(zhì)量過(guò)關(guān)是沒(méi)有問(wèn)題，這種發(fā)展思路在產(chǎn)業(yè)欣欣向榮之際是不會(huì)存在任何問(wèn)題，但隨著介入的企業(yè)越來(lái)越多，產(chǎn)品的同質(zhì)化也將是必然趨勢(shì)。屆時(shí)，試問(wèn)我們儀器儀表又如何確保競(jìng)爭(zhēng)優(yōu)勢(shì)，如何保證自己立于不敗之地？
 
　　在面臨這一系列的經(jīng)濟(jì)、技術(shù)、貿(mào)易出口問(wèn)題，儀器儀表行業(yè)的競(jìng)爭(zhēng)也將會(huì)迎來(lái)新的格局，企業(yè)面臨的挑戰(zhàn)更嚴(yán)峻。因此，只要我們儀器儀表企業(yè)能夠堅(jiān)持創(chuàng)新發(fā)展，契合產(chǎn)業(yè)大趨勢(shì)，就能夠不斷壯大自身，而這，既是未來(lái)發(fā)展的要求，也是企業(yè)長(zhǎng)遠(yuǎn)發(fā)展的根本之道。
 
　　可以看到，行業(yè)的金卡智能、三川智慧、寧波水表等企業(yè)已經(jīng)做出了改變，經(jīng)過(guò)多年來(lái)謀篇布局，*提升企業(yè)綜合實(shí)力，已經(jīng)有了創(chuàng)收，這一點(diǎn)在上半年的業(yè)績(jī)報(bào)告中已經(jīng)明顯體現(xiàn)。
 
　　很顯然，經(jīng)濟(jì)低迷，企業(yè)要做的唯有研發(fā)和創(chuàng)新，在技術(shù)上突破壁壘，才能在市場(chǎng)有一席之地。
 

多聲路超聲流量計(jì)準(zhǔn)確度解決方案與積分方法分析

1、前言：隨著能源和水資源的性匱乏 , 一批關(guān)系國(guó)計(jì)民生的大型水利工程和引水調(diào)水工程在我國(guó)迅速發(fā)展, 如三峽水利樞紐 、南水北調(diào)工程等 。這些工程項(xiàng)目中經(jīng)常包含一些口徑和流量都很巨大的管道 , 如水電站機(jī)組進(jìn)水管道等 , 常規(guī)流量計(jì)無(wú)法適應(yīng) 。
  近年開(kāi)發(fā)應(yīng)用的多聲路超聲流量計(jì) ,較好地解決了大口徑水流量測(cè)量的技術(shù)難題 , 流量計(jì)制造不受管道口徑的限制, 多聲路配置可以適應(yīng)較為復(fù)雜的流道結(jié)構(gòu)和流態(tài)分布 , 故超聲流量計(jì)已成為大口徑水流量測(cè)量的***佳技術(shù)選擇[ 1]。超聲流量計(jì)屬于速度式流量計(jì) , 通過(guò)測(cè)量超聲波在流體中順流和逆流傳播的時(shí)間差來(lái)計(jì)算聲路線(xiàn)上的平均流速 , 并通過(guò)測(cè)量多條聲路速度來(lái)加權(quán)平均計(jì)算待測(cè)截面上的平均流速 。超聲流量計(jì)需要進(jìn)行復(fù)雜運(yùn)算才能得到***終的平均流速和流量 , 其采用的數(shù)學(xué)模型準(zhǔn)確度對(duì)整個(gè)測(cè)量準(zhǔn)確度有著非常重要的意義 。
  顯然, 聲路數(shù)越多, 流量積分準(zhǔn)確度就越高, 但聲路數(shù)的增加會(huì)大大增加流量計(jì)的成本 , 所以選擇合理的聲路位置 、匹配的權(quán)重系數(shù)就十分重要。超聲流量計(jì)可以從提高時(shí)間及幾何量的測(cè)量準(zhǔn)確度來(lái)提高流量計(jì)本身的測(cè)量準(zhǔn)確度 , 但是數(shù)值積分引入的測(cè)量誤差 (簡(jiǎn)稱(chēng)積分誤差 )始終為流量計(jì)所無(wú)法繞開(kāi);特別是隨著微電子技術(shù)以及信號(hào)處理技術(shù)的飛速發(fā)展, 聲路速度的測(cè)量準(zhǔn)確度越來(lái)越高[ 2-3], 積分誤差逐漸成為提高流量測(cè)量準(zhǔn)確度的瓶頸 。對(duì)于常用的 4聲路 、8聲路超聲流量計(jì) , 積分誤差是流量計(jì)測(cè)量誤差的主要來(lái)源 ;對(duì)于聲路數(shù)更多的情況 , 積分誤差也是流量計(jì)測(cè)量準(zhǔn)確度的重要影響因素。即使流量計(jì)安裝滿(mǎn)足前后直管段長(zhǎng)度要求 , 其積分誤差也依然存在 ;而若由于場(chǎng)地或資金限制導(dǎo)致無(wú)法滿(mǎn)足時(shí), 積分誤差就更應(yīng)引起注意 。

2、超聲波流量計(jì)積分原理：
  超聲流量計(jì)利用超聲波在流體中傳播的時(shí)間存在差異的特性 , 由置于待測(cè)截面兩側(cè)的一對(duì)換能器(如圖 1所示 ), 測(cè)量超聲波順流與逆流傳播的時(shí)間 td, i、 tu, i, 得到相應(yīng)聲路上的平均軸向速度[ 4](簡(jiǎn)稱(chēng)聲路速度 )：



式中 :Li為聲路長(zhǎng)度, i為聲路角 。對(duì)于單聲路流量計(jì),截面平均流速與該聲路速度存在特定關(guān)系 , 但易受到流速分布廓形的影響 。為了提高流量計(jì)的測(cè)量準(zhǔn)確度, 在待測(cè)截面上平行地布置多條聲路 , 獲得的聲路速度可以代表待測(cè)截面上相應(yīng)平行條帶內(nèi)的平均速度 , 如圖 2所示 , 并依據(jù)各條帶所占的權(quán)重系數(shù) ωi, 用加權(quán)求和的方法計(jì)算流量 ,




圖 1　聲路速度的測(cè)量
式 (2)加權(quán)求和計(jì)算流量的方法實(shí)際利用了數(shù)值積分的原理 , 通過(guò)有限個(gè)聲路采樣點(diǎn)計(jì)算 l(z)·  v(z)的值 , 來(lái)逼近其在區(qū)間 [ -R, R] 上的定積分 :Q=∫R-Rl(z) v(z)dz=R∫1-1l(t R) v(t R)dt≈R∑Ni=1ωil(tiR) v(tiR)


圖 2　流量積算示意圖
式中 :z=t R為聲路高度, t為相對(duì)聲路高度 , l(z)為聲路寬 度。 式 (3)將數(shù)值 積分 變換 到 [ -1, 1] 區(qū)間上來(lái) , 方便不同半徑時(shí)的計(jì)算 。若代入聲路寬度 l(z)=2 R2-z2=2R 1 -t2, 則流量為 :Q=2R2∫1-1ρ(t)·  v(t R)dt≈ 2R2∑Ni=1ωi·  v(tiR)· ρ(ti)(4)式中 :ρ(t)= 1 -t2。相對(duì)于梯形公式、辛普森公式等插值型積分要求采樣點(diǎn)固定甚至等距而言 , Gauss積分方法則是在采樣點(diǎn)個(gè)數(shù)一定 、位置自由選擇等限定下積分精度***高的一種方法[ 5]。圓管中的超聲流量計(jì)一般采用 Gauss-Jaccobi積分法來(lái)確定聲路的***優(yōu)位置 ti和相應(yīng)的權(quán)重系數(shù) ωi,IEC41[ 6]及 PTC18[ 7]規(guī)程中已有不同聲路數(shù) N時(shí)的聲路高度和權(quán)重系數(shù) , 一般依此位置及系數(shù)安裝超聲探頭并計(jì)算流量 。

3、聲路高度與權(quán)重系數(shù)的推算：
3.1、Gauss-Jaccobi方案：
  根據(jù) Gauss積分理論 , 相對(duì)聲路高度 ti為帶權(quán) ρ(t)的正交多項(xiàng)式 PN(t)的根。對(duì)于圓形斷面的管道 , 由式(4)可以看到權(quán)函數(shù)為 ρ(t)= 1 -t2, 這正好為古典的 Jaccobi正交多項(xiàng)式的權(quán)函數(shù) (1 +t)α(1 -t)β的特殊情況, α=β =κ=0.5 , 故 Jaccobi正交多項(xiàng)式可以由式(5)遞推[ 8]:Pj+1=[ Pj(2j+2κ+1)x-Pj-1(j+κ)] (j+κ+1)(j+1)(j+2κ+1)(5)起始值 P-1=0, P0=1 。進(jìn)一步可以計(jì)算 PN的 N個(gè)根 ti, 即相對(duì)聲路高度。而對(duì)于矩形管道或明渠 , 權(quán)函數(shù) ρ(t)=1 ,則為 Legendre正交多項(xiàng)式問(wèn)題 。此為圓形管道宜采用 Gauss-Jaccobi方案、方形管道宜采用 Gauss-Legendre方案的由來(lái) 。
  權(quán)重系數(shù) ωi可通過(guò)下面的積分得到:ωi=1ρ(ti)∫1-1ρ(t)∏Nk=0, k≠it-tkti-tkdt (6)結(jié)合 Gamma函數(shù) Γ(x)=∫∞0e-ttx-1dt, 式 (6)可以轉(zhuǎn)化為容易計(jì)算的形式:ωi=1ρ(ti)Γ2(κ+N)(κ+N)22κ+1Γ(N+1)Γ(2κ+N+1)PN-1(ti)P′N(ti)(7)式中 :PN-1(ti)為 N-1階 Jaccobi正交多項(xiàng)式 , P′N(ti)為 N階 Jaccobi正交多項(xiàng)式的導(dǎo)數(shù) 。值得注意的是 , 該權(quán)重系數(shù)與數(shù)學(xué)上的 Gauss-Jaccobi積分的權(quán)重系數(shù)稍有差別 , 出于實(shí)際聲路寬度可能與安裝預(yù)期值有所差異的考慮 , 將 ρ(ti)從權(quán)重系數(shù)中剔除出來(lái), 而在積分時(shí)采用其實(shí)測(cè)值 ρ(ti)=Lisin i/2代入計(jì)算 , 如式(4)所示 。表 1給出了 4聲路及 9聲路時(shí)的相對(duì)聲路高度 ti和相應(yīng)的權(quán)重系數(shù) ωi。由于 Gauss-Jaccobi方案中 α=β =κ=0.5,可以經(jīng)數(shù)學(xué)推導(dǎo)得到相對(duì)聲路高度和權(quán)重系數(shù)的簡(jiǎn)化計(jì)算公式[ 9]:ti=cosiπN+1, i=1, 2, …, Nωi=1ρ(ti)πN+1sin2iπN+1(8)式中 :N為聲路數(shù)量 。

3.2、OWICS方案：
  對(duì)于聲路速度分布  v(z)可由相應(yīng)階數(shù)的代數(shù)多項(xiàng)式表達(dá)的情況 , Gauss-Jaccobi方案不存在截?cái)嗾`差;但實(shí)際聲路速度分布與理想的代數(shù)多項(xiàng)式表達(dá)形式之間存在較大差異 , 特別是無(wú)法體現(xiàn)管壁處流速為零這動(dòng)特性 , 導(dǎo)致流量積分結(jié)果偏高, 而聲路數(shù)越少 , 流量計(jì)算值偏高的趨勢(shì)越強(qiáng)烈 。考慮到充分發(fā)展的圓管紊流的實(shí)際聲路速度分布與形如 (1 -t2)1/10的指數(shù)分布接近 , 可以將其從聲路速度 v(z)中提取出來(lái), 使  v′(z)接近于 1, 從而  v′(z)更容易為代數(shù)多項(xiàng)式所表達(dá)。 進(jìn)一步將其與權(quán)函數(shù) ρ(t) =1 -t2合成得到 : v(z)ρ(t)= v′(z)(1 -t2)1/ 10ρ(t)= v′(z)(1 -t2)0.6= v′(z)ρ′(t) (9)然后按照新的權(quán)函數(shù) ρ′(t)=(1 -t2)0.6, 即 α =β =κ=0.6 , 可以計(jì)算不同聲路數(shù) N時(shí)的相對(duì)聲路高度ti和權(quán)重系數(shù) ωi, 其中 4聲路及 9聲路的結(jié)果見(jiàn)表 1。這種算法稱(chēng)為***佳圓斷面 (OWICS)方案 , 實(shí)際也是基于正交多項(xiàng)式的 Gauss積分方案[ 10]。

3.3、權(quán)重系數(shù)與面積平均方案的比較:
  Gauss積分方案有其數(shù)值分析上的基礎(chǔ) , 在確定的聲路高度條件下 , 不同聲路的權(quán)重與其所在條帶的面積相關(guān), 但并非成正比關(guān)系 。針對(duì)計(jì)算得到的相對(duì)聲路高度ti, 計(jì)算了各聲路所在條帶的面積 , 分割為中間的梯形與兩側(cè)的弓形進(jìn)行計(jì)算 :Si=12(xi+xi+1)Δzi+ΔβiR-sinΔβiR2(10)式中 :Δβi= arcsinti-arcsinti+1 為兩側(cè)弓形的圓心角, xi與 xi+1為聲路條帶的上下兩側(cè)的弦長(zhǎng), 該弦正好均分相鄰聲路之間的面積 , 弦的位置 xi可利用條帶面積相等而迭代求解得到。若認(rèn)為聲路速度為所在條帶內(nèi)的平均速度 , 則聲路條帶的面積 Si可以作為截面平均速度計(jì)算時(shí)的權(quán)重系數(shù) , 為與前面的 Gauss積分方案對(duì)比 , 權(quán)重系數(shù)需進(jìn)行處理 , ωs, i=Si/ρ(ti)。對(duì)比 Gauss積分方案與面積平均方案的的權(quán)重系數(shù)發(fā)現(xiàn) , 兩者有所差別 , 如表 1所示 。 Gauss積分方案并非簡(jiǎn)單的面積加權(quán)求和 , 其權(quán)重系數(shù)與面積加權(quán)平均相比 ,對(duì)中間聲路偏大而對(duì)邊緣聲路偏小 , 說(shuō)明 Gauss積分更多的是強(qiáng)調(diào)中間聲路速度對(duì)平均流速的貢獻(xiàn), 這可能也是 Gauss積分方案對(duì)復(fù)雜聲路速度分布的適應(yīng)性要優(yōu)于面積平均方案的原因所在。



表 1　流量積分節(jié)點(diǎn)位置及求積系數(shù)


4、流量積分準(zhǔn)確度分析:
4.1、不同聲路數(shù)時(shí)的積分誤差：
  超聲流量計(jì)利用有限聲路還原整個(gè)截面上的速度分布情況, 顯而易見(jiàn)的是聲路越多, 流量積分越準(zhǔn)確 。根據(jù)圓管紊流理論 ,圓形管道中充分發(fā)展的流速分布通?？捎弥笖?shù)分布來(lái)描述 , 指數(shù) n與雷諾數(shù)和相對(duì)粗糙度有關(guān)系 , 下面以 n=9時(shí)v=(1 -r)1 /9　　(0 ≤ r≤ 1) (11)為例來(lái)分析積分誤差E=Q-QiQi× (12)的影響因素 。式中 , Q為不同方案時(shí)的數(shù)值積分結(jié)果,Qi為數(shù)學(xué)積分的結(jié)果 。圖 3給出了 3種方案在不同聲路數(shù)時(shí)的積分誤差?？梢钥吹?, Gauss-Jaccobi方案的積分誤差均大于零, 且隨著聲路數(shù) N的增加 , 積分誤差迅速減小 , 在 N=6之后積分誤差已降至 0.05%以下, 且繼續(xù)增加聲路數(shù)不再引起很大變化 , 說(shuō)明對(duì)于簡(jiǎn)單平順的流動(dòng) , 過(guò)多的聲路數(shù)對(duì)提高流量計(jì)準(zhǔn)確度意義不大 。對(duì)于面積加權(quán)平均方案, 積分誤差也隨著聲路數(shù) N的增加而減小 , 但減小速度甚慢,36聲路時(shí)仍然保持 0.05%左右的積分誤差 。 兩者比較可以發(fā)現(xiàn) , Gauss積分方案要比面積平均方案具有更高的準(zhǔn)確度。對(duì)于***佳圓斷面方案, 積分誤差在零上下浮動(dòng),隨著聲路數(shù)的增加, 浮動(dòng)范圍逐漸減小 。對(duì)于聲路數(shù)N =2 ～ 9 , OWICS方案的積分誤差比 Gauss-Jaccobi方案還要略小 , 說(shuō)明前者確實(shí)優(yōu)于后者 ;但在聲路數(shù)較大時(shí),OWICS方案的優(yōu)勢(shì)已比較微弱。



圖 3　圓管紊流條件下的積分誤差
  為了進(jìn)一步探討 3種積分方案的差異 , 利用擾流流場(chǎng)經(jīng)典公式：
  u(r, θ)=(1 -r)1 /n+mr(1 -r)1/kf(θ) (13)對(duì)積分誤差進(jìn)行進(jìn)一步分析 。式中 , n、 k、 m、 f(θ)為可調(diào)的參數(shù), θ在 [ 0, 2π] 內(nèi)。圖 4給出了 3種方案的積分誤差計(jì)算結(jié)果 , 點(diǎn)值為利用式 (13)描述的一組流場(chǎng)下的積分誤差均值 , 而豎線(xiàn)則表示其標(biāo)準(zhǔn)差大小 。其規(guī)律與圓管紊流情況類(lèi)似 , 兩種 Gauss積分的效果要比面積加權(quán)平均好得多 , 前者積分誤差不僅能夠隨著聲路數(shù)的增加迅速收斂到零附近 , 而且波動(dòng)幅度也比同聲路數(shù)的面積加權(quán)平均方案為小 。另外進(jìn)一步說(shuō)明, OWICS方案的優(yōu)勢(shì)主要體現(xiàn)在聲路數(shù)較少時(shí) ,若聲路數(shù)量受到限制 ,以采用 OWICS方案為佳 。


圖 4　擾流流場(chǎng)條件下的積分誤差
4.2、積分誤差的來(lái)源：
  Gauss數(shù)值積分過(guò)程具有 2N-1階代數(shù)精度 , 即若實(shí)際聲路速度分布可用 2N-1次代數(shù)多項(xiàng)式表達(dá) , 則數(shù)值過(guò)程引入的截?cái)嗾`差為零 。但是由于實(shí)際聲路速度分布曲線(xiàn)受邊界條件等的影響 , 造成其無(wú)法由 2N-1次多項(xiàng)式來(lái)表達(dá) , 所以流量計(jì)算的數(shù)值積分過(guò)程將引入一定的截?cái)嗾`差 , 其大小與實(shí)際聲路速度分布曲線(xiàn)有關(guān)系, 一般來(lái)說(shuō)待測(cè)流場(chǎng)廓形越復(fù)雜, 積分誤差越大 。仍以式(11)中的指數(shù)流速分布為例 , 來(lái)說(shuō)明積分誤差的來(lái)源。該速度分布在半徑為 1的管道中的流量可以通過(guò)數(shù)學(xué)積分得到 Q =2.678 620 。圖 5中的聲路速度曲線(xiàn)由該指數(shù)函數(shù)積分得到 , 具有對(duì)稱(chēng)的分布特性 , 中間區(qū)域較為平坦 , 邊緣附近迅速降低為零 , 若按 Gauss-Jac-cobi方案來(lái)取 N =5進(jìn)行采樣 , 數(shù)值積分的結(jié)果為 Q =2.681 950 , 與真值相差 0.12%。圖 5中另給出了過(guò)該組采樣點(diǎn)的代數(shù)多項(xiàng)式族 (次數(shù)不大于 9, 其中 3條賦予了邊界為零的約束 ), 雖然這一組曲線(xiàn)分布較為復(fù)雜 , 但其描述的聲路速度分布的積分結(jié)果 Q=2.681 950 , 與數(shù)值積分的結(jié)果*相同 , 這也驗(yàn)證了 Gauss積分具有 2N-1階代數(shù)精度的結(jié)論 。 事實(shí)上一般不會(huì)出現(xiàn)實(shí)際聲路速度曲線(xiàn)正好為代數(shù)多項(xiàng)式描述的情況 , 故積分過(guò)程總是存在積分誤差 。


圖 5　聲路速度分布示意圖

  Gauss-Jaccobi積分法和 OWICS積分法均建立在聲路速度分布可由代數(shù)多項(xiàng)式來(lái)表達(dá)的假設(shè)基礎(chǔ)上, 部分研究者認(rèn)為聲路速度分布根本不適合由代數(shù)多項(xiàng)式來(lái)表達(dá)[ 12], 尤其是對(duì)于流場(chǎng)受到擾動(dòng)的情況 , A.Nichtaw-itz[ 13]認(rèn)為利用 PCHIP多項(xiàng)式(分段立體赫爾密特插值多項(xiàng)式 )來(lái)表達(dá)流速分布比較合適, 在修正邊壁流場(chǎng)的情況下 , 利用 4聲路的 PCHIP積分法與 9聲路 Gauss-Jaccobi積分法大致具有相同的積分誤差 。另外也有學(xué)者[ 14-15]認(rèn)為流量計(jì)現(xiàn)場(chǎng)流動(dòng)非常復(fù)雜 、根本不可預(yù)測(cè) , 而采用復(fù)雜的聲道布置并輔以神經(jīng)網(wǎng)絡(luò)進(jìn)行處理的思路 , 來(lái)推進(jìn)超聲流量計(jì)數(shù)學(xué)模型的研究 。


5、結(jié)論：
  本文分析了超聲流量計(jì)的流量積分原理 , 并結(jié)合數(shù)值分析理論, 推導(dǎo)了超聲流量計(jì)聲路布置方案 (Gauss-Jaccobi方案及 OWICS方案 )的聲路高度與權(quán)重系數(shù) 。通過(guò)對(duì)比 Gauss積分方法與面積加權(quán)平均方法在聲路高度相同時(shí)的權(quán)重系數(shù)差別, 發(fā)現(xiàn) Gauss積分方法更加重視中間部分聲路對(duì)平均流速的貢獻(xiàn) , 這可能是其具有更高準(zhǔn)確度的重要原因 。另外還分析了 Gauss積分方案流量積算誤差的產(chǎn)生原因 , 并對(duì)比了聲路數(shù)對(duì)積分誤差的影響 。積分誤差來(lái)源于實(shí)際聲路速度無(wú)法為代數(shù)多項(xiàng)式所表達(dá) ;聲路數(shù)越多, 積分準(zhǔn)確度越高 , 但過(guò)多聲路數(shù)對(duì)提高流量計(jì)準(zhǔn)確度意義也不大 。