摘要：小流量測(cè)量時，差壓式流量(liang)計 輸出的差壓與(yǔ)流量之間是非線(xiàn)性關系，在分散控(kòng)制系統（DCS）中直接實(shí)施該非線性關系(xì)較困難。根據已知(zhī)标準孔闆的徑比(bi)，用NURBS非❗均勻有理函(han)數，拟合在特定應(yīng)用條件下的标準(zhǔn)📐孔闆流💜量系數公(gōng)式；并用簡單的乘(cheng)法和加法運算，在(zài)DCS中用NURBS函數表示該(gai)非線性♉關系的輸(shu)入和輸出關系；最(zuì)後用非線性叠代(dai)算法确定在小流(liu)量條件下的差壓(ya)和流量關系，從而(ér)實現小流量測量(liàng)的在線非線性補(bǔ)償，提高了流量測(cè)量🌐的精度。
　　差壓式(shì)流量計是常用的(de)流量測量儀表。 标(biāo)準孔闆 的流量系(xì)數經Reader-Harris/Gallagher修改，于1998年被(bei)采納作爲标準孔(kong)闆流出系⭐數的👄計(ji)算公式。它對小流(liú)量時差壓式流量(liang)計的補償提供了(le)理論基礎🔞，但在分(fèn)散控制系統（DCS）中實(shi)現有困難，爲此，提(ti)出兩種實施方法(fǎ)：直接用Reader-Harris/Gallagher公式，但在(zai)DCS上Reader-Harris/Gallagher公式㊙️實施困難(nan)；針對特定标準孔(kǒng)闆，用NURBS函數拟合标(biāo)準孔闆流出系數(shù)的🔞Reader-Harris/Gallagher計算公式，并在(zai)DCS中實現。該方法既(ji)解決了小流量在(zài)線補償的實施問(wèn)題，也提💁高了差壓(yā)♈式流量計的測量(liang)範圍度和精度。
1NURBS樣(yàng)條函數
1.1Ｂ樣條基函(han)數
　　Ｂ樣條基樣條（basicspline）。1946年(nián)由舍恩貝格（Schoenberg）提出(chū)，并在1972年由德布爾(ěr)💚和考📐克斯（deboor-Cox）分❄️别獨(du)立給出Ｂ樣條計算(suàn)的标準算法［1－2］。理論(lùn)上常采用截尾幂(mì)函數🌈的差商定義(yì)Ｂ樣條曲線，實際應(ying)用則常采用Ｂ樣條(tiáo)的遞推🚶定義。
Ｂ樣條(tiao)曲線采用控制頂(ding)點定義曲線［1－2］。曲線(xiàn)方程可描♋述爲


　　式(shi)中：Ｐｉ———控制多邊形的(de)頂點，ｉ＝0，1，…，ｎ；Ｎｉ，ｋ（ｕ）———ｋ次（ｋ－1次）Ｂ樣條基(jī)函數，ｉ＝0，1，…，ｎ。
　　其中，每個ｋ次(ci)規範Ｂ樣條基函數(shù)稱爲規範Ｂ樣條，或(huo)簡稱Ｂ樣條。由于它(ta)由非遞減節點矢(shi)量ｕ的序列Ｔ：ｕ0≤ｕ1≤…≤ｕｎ＋ｋ所決定(ding)的ｋ次分段多項💯式(shì)，因而，稱🐕爲ｋ－1次多項(xiang)式樣條。
根據德布(bù)爾－考克斯的遞推(tui)公式，曲線方程可(ke)寫爲

式中：ｉ，ｋ———下标，ｉ表(biao)示序号，ｋ表示次數(shù)。
1.2三次非均勻有理(lǐ)Ｂ樣條函數
三次非(fēi)均勻有理Ｂ樣條函(hán)數描述爲

　　式中：ｗｉ———權(quan)因子，分别與控制(zhi)頂點Ｐｉ相聯系，（ｉ＝0，1，…，ｎ）；Ｎｉ，ｋ（ｕ）———節點(dian)矢量，ｕ＝［ｕ0，ｕ1，…，ｕｎ＋ｋ＋1］按遞推公式(shì)确定的ｋ次規範Ｂ樣(yàng)條基函數；Ｐ1，Ｐ2，Ｐ3，Ｐ4———分子系(xi)數，爲矢量；Ｑ0，Ｑ1，Ｑ2，Ｑ3———分母系(xi)數。Ｂ樣條基函數的(de)遞推公式見式（3）～式(shì)（4）。
　　在數控技術中，NURBS曲(qǔ)線插補算法将定(dìng)義NURBS曲線的控制💰頂(dǐng)💰點、權因子🏃🏻、節點矢(shǐ)量和進給速度等(deng)作爲ＮＣ程序指令，在(zài)ＣＮＣ系統生成NURBS曲線，驅(qu)動機床運動，加工(gong)出NURBS曲線的形狀，這(zhè)就是NURBS曲🍉線插補。在(zai)非線性補償環節(jie)中應用🏃‍♀️的NURBS曲線，可(kě)根據應用要求選(xuǎn)用不同的階次。
2差(cha)壓式流量計在非(fēi)線性補償中的應(ying)用
2.1差壓式流量計(jì)的問題
　　差壓式流(liú)量計是應用久遠(yuǎn)的流量計之一，其(qí)測量原理是 孔闆(pan)流量計 上遊側與(yǔ)下遊側之間産生(shēng)的靜壓差與流過(guo)該裝置的👨‍❤️‍👨流體流(liu)量之間存在下列(liè)關系：

　　當滿足0.2≤β≤0.6時，流(liu)出系數Ｃ的不确定(dìng)度爲0.5％。其他條件下(xia)，不确定度會有所(suǒ)增加。其中，Ｃ經Reader-Harris/Gallagher修改(gai)，可表示爲

當工藝(yì)管道的管道内徑(jing)Ｄ＜71.12ｍｍ時，增加下列項：

　　式(shì)中：β———節流孔直徑ｄ與(yu)Ｄ之比，即β＝ｄ/Ｄ；ＲｅＤ———根據Ｄ和流(liu)體流量等數據🛀計(jì)算🔞出的雷諾數；Ｌ1———孔(kong)闆上遊端面到上(shàng)遊取壓口的距離(lí)ｌ1除✊以Ｄ得出的商。

　　式(shì)中：Ｌ′2———孔闆下遊端面(miàn)到下遊取壓口的(de)距離Ｌ′2除以Ｄ得出的(de)商。對🔞不同取壓方(fāng)式，Ｌ1和Ｌ′2的值不同
根(gen)據Reader-Harris/Gallagher公式，可畫出不(bú)同管道直徑和不(bú)同取壓方式下，Ｃ與(yǔ)ＲｅＤ，β之間✍️的💔關系曲面(miàn)。角接取壓，Ｄ＝150ｍｍ時，Ｃ與ＲｅＤ，β的(de)關系如圖1所示。
　　從(cong)圖1可見，當Ｄ确定後(hòu)，如果ｄ也确定，則當(dāng)流體的ＲｅＤ大于某限(xiàn)值時，其🌈Ｃ可基本穩(wěn)定在某個規定的(de)值。通常在0.60～0.61，而測量(liàng)不确定度應滿足(zu)小于0.5％。
　　角接取壓，Ｄ大(dà)于72.12ｍｍ時，β在0.4～0.5，Ｃ與ＲｅＤ的關系(xi)見表1所列。根據表(biao)1中數據的分析，可(kě)以發現，當最大流(liú)量與最小流量之(zhi)比爲10∶1時，即小流量(liàng)時，其Ｃ的誤差可達(da)2％。但如果最小雷諾(nuo)數大㊙️于2×104，則Ｃ的誤差(cha)就可小于0.5％。該條件(jiàn)是采用差壓式流(liú)量計有最小雷諾(nuò)數限制的原因。由(yóu)🙇🏻于受到流體流速(sù)的限制，最大流量(liàng)不能設置很大。又(you)由于小🍓流量時，ＲｅＤ成(cheng)比例縮小，在Ｃ的非(fēi)線☂️性影響下造成(chéng)流量測量的精度(dù)下降。因而，該😄情況(kuàng)是差壓式流量計(ji)的🏃‍♀️範圍度不能較(jiao)大的原因。其根本(běn)原因⛷️是在流量小(xiao)時，ＲｅＤ也小，這時，Ｃ與ＲｅＤ之(zhi)間存在較大的非(fei)線性關系，造成小(xiǎo)流量時流量測量(liang)誤差大，和流量測(ce)量範圍度不能大(da)的結果。
　　解決該類(lei)非線性關系的最(zui)好方法是進行非(fei)線性補償［6－7］。對⚽差壓(ya)式流量計由于存(cun)在叠代運算，加上(shang)在DCS中進行式（7）的運(yùn)算比較困難，因此(ci)，實際應用時可采(cǎi)用兩種實現的方(fang)法。


2.2差壓式流量計(jì)理論補償方法
　　當(dang)實際差壓流量計(ji)已安裝在工藝管(guǎn)道中時，可采用理(lǐ)論補🏃‍♂️償🤟方法。該方(fāng)法根據Reader-Harris或Gallagher公式，根(gen)據已知的β和取壓(ya)方式，計🎯算出Ｃ與ＲｅＤ之(zhi)間的關系。根據兩(liǎng)者關系，有多種方(fang)法實現補償，如采(cǎi)用🌈多段折線近似(si)法進行補償；采用(yòng)拟合函數進行補(bu)償；也可用其他非(fei)線性環節實現，例(lì)如，神經網🙇‍♀️絡等。
　　示(shi)例是已經安裝的(de)某節流裝置，已知(zhi)Ｄ＝100.00ｍｍ，β＝0.40，角接取壓方式。爲(wèi)提高拟合精度，取(qu)點較多，其計算結(jie)果見表2所列。采用(yong)NURBS函數進行拟合，其(qí)NURBS函數表示爲


　　從表(biǎo)2可見，用式（10）拟合Reader-Harris或(huo)Gallagher計算公式，具有很(hěn)高的精度，最大誤(wu)差小于0.013％。因此，可直(zhi)接根據ＲｅＤ确定Ｃ。
2.3差壓(ya)式流量計實際标(biāo)定補償方法
　　在新(xīn)建項目中，可用實(shi)流标定的方法确(que)定不同流量時ＲｅＤ與(yǔ)✔️Ｃ的關系曲線，采用(yòng)上述拟合方法确(que)定其非🏃‍♀️線性關系(xi)。最簡單的方法是(shì)用多段折線方法(fa)拟合，但需設置段(duàn)數✊，并用内插方法(fa)确定其輸出值［8－10］。例(lì)如，DCS可以實現其他(ta)非線性環節［11］，也可(ke)采用神經網絡實(shí)現非線性關系，或(huò)用有關方法獲得(de)該非線性關系的(de)描述，在此不多述(shù)。采用NURBS函數拟合在(zai)特定徑比條件下(xià)的ＲｅＤ與✍️Ｃ之間的非線(xiàn)性🔴關系，并實際實(shí)施。将NURBS函數表示爲(wèi)下列形式。

　　利用可(ke)編程控制器編程(chéng)語言中的可重用(yòng)性，發現NURBS函數的🤟基(ji)本算式是ｙ＝Ａｘ＋Ｂ。爲此，可(ke)編寫ＡＸＢ函數實現。NURBS函(han)數的程序✉️實現如(rú)圖2所示😘。

2.4DCS中在線非(fei)線性補償關系的(de)實現
　　爲在線實施(shī)，先建立Online功能塊，用(yong)于實現非線性的(de)ＲｅＤ與Ｃ的關系，再針對(dui)實際應用，編寫主(zhǔ)程序，它由QCal，ReCal和NUBRS3個功(gong)能塊組成💜。以Ｃ作爲(wèi)反饋變量，該程序(xù)爲叠代程序。QCal功能(néng)塊用🔞于計算流體(ti)流量，ReCal功能塊用于(yú)計算ＲｅＤ，NUBRS函數用于計(ji)算不✂️同ＲｅＤ下的Ｃ。
　　在線(xian)實現時，将Online與用常(cháng)規開方計算的結(jie)果進行比較，确定(dìng)🙇🏻其誤差。如圖3所示(shi)。

　　從圖3可見，當實際(ji)差壓輸入信号是(shi)205.2Ｐａ時，實際流量應爲(wei)4.983542ｋｇ/ｓ。如🔞果沒有非線性(xìng)補償，顯示值是4.9216ｋｇ/ｓ，顯(xiǎn)示值偏小，誤🈲差達(dá)1.24％。通過該方法的⛷️補(bu)償，使原流量計的(de)範圍度提高到接(jiē)近10∶1。
3結論
　　爲提高差(cha)壓式流量計的流(liú)量測量精度和範(fan)圍度🌈，可對小流☂️量(liàng)進行在線非線性(xing)補償。由于标準孔(kǒng)闆Ｃ的計算公式實(shí)現比較複✨雜，在DCS中(zhōng)計算較困難，因而(er)采用NURBS函數［9］來拟合(he)該非線性關系，并(bing)用它計算小流量(liang)時的Ｃ，通過該非線(xian)性補償的🧡方法，提(tí)高了小流量🔴測量(liàng)精度，同時提高了(le)測量範圍度。

以上(shàng)内容來源于網絡(luo)，如有侵權請聯系(xi)即删除!