摘 要：本文建立了新型HKC科氏質(zhì)量流量計的力學(xué)模型，推導(dǎo)出其運動微分方程，并得到靈敏度計算表達式。

關(guān)鍵字：質(zhì)量流量計 靈敏度 科里奧利力 波紋管

    1 引言

    本文所介紹的新型HKC科氏質(zhì)量流量計在普通U型管上裝入了波紋管，由于波紋管剛度較小，導(dǎo)致流量管兩側(cè)位移差較大，提高了靈敏度。HKC科氏質(zhì)量流量計的實際振動形態(tài)很難用精確的函數(shù)表示，而一些近似的振動形式已經(jīng)能夠保證計算精度，本文在保持原結(jié)構(gòu)的結(jié)構(gòu)特征和受力特點不失真的前提下，建立了其力學(xué)模型，推導(dǎo)出運動微分方程，分析了靈敏度。

    2 結(jié)構(gòu)與工作原理

    如圖1所示，電磁激振器驅(qū)動U型管振動，被檢測流體以質(zhì)量流量Q_m流過管內(nèi)。由于管的彎曲振動，管內(nèi)流體產(chǎn)生了科氏加速度，流體每一個微元受到與科氏加速度相應(yīng)的科氏力的作用，其反作用力作用于U型管。由于在管的兩側(cè)液體流速相反，所以受到方向相反的科氏力的作用，從而使U型管發(fā)生相對扭轉(zhuǎn)變形，產(chǎn)生扭轉(zhuǎn)振動。此扭轉(zhuǎn)振動疊加在彎曲振動上，使管上B、D兩個檢測點通過振動中心有一個時間差Δt（通過檢測器測出），（φ·2R=z_B^&·Δt，式中z_B^&為B截面通過振動中心的速度，（φ為扭角，是Q_m的表達式，從而可由上式得到質(zhì)量流量。

    3 方程推導(dǎo)

    力學(xué)模型如圖2所示。A、E端為固定端，C處受電磁驅(qū)動力F的作用，F(xiàn)=Pcosωt。由于波紋管的剛度遠(yuǎn)小于鋼管，其變形遠(yuǎn)大于鋼管，結(jié)構(gòu)近似彈簧質(zhì)量系統(tǒng)，鋼管部分視為剛體，波紋管視為彈簧，其運動分為繞x軸的振動和繞y軸的振動兩個獨立運動。

    3.1 繞x軸的振動微分方程的推導(dǎo)

    繞x軸的運動近似為如圖3所示的彈簧質(zhì)量系統(tǒng)。

    圖中N_x為彈性回復(fù)力，F(xiàn)_dx為空氣阻力，F(xiàn)為電磁激振力。

    運動微分方程可寫為：

    J_xθ+C_xθ+K_xθ=M₀cosωt    （1）

    式中：J_x--系統(tǒng)對x軸的轉(zhuǎn)動慣量；
          C_x--等效粘性阻尼系數(shù)；
          K_x--系統(tǒng)的剛性系數(shù)；
          M₀cosωt--激振力產(chǎn)生的彎矩[M₀=P（L₁+L₂+R）]。

    3.1.1 方程（1）系數(shù)的計算

    （1）J_x的計算
    
    設(shè)波紋管單位長度上的質(zhì)量為，鋼管單位長度上的質(zhì)量為，將其分為三部分計算（兩個波紋管為第一部分，兩段直管為第二部分，半圓為第三部分）。

    

    計算第三部分轉(zhuǎn)動慣量：

    半圓繞x′軸的轉(zhuǎn)動慣量為
    質(zhì)心到x′軸的距離：

    

    繞質(zhì)心的轉(zhuǎn)動慣量為

    

    （2）阻尼系數(shù)C_x的確定

    將方程（1）化為標(biāo)準(zhǔn)方程得，

    

    通過實驗測量有阻尼自由振動的振幅X_i和周期T，則得減幅系數(shù)，將（2）式代入得：

    

    這種方法簡單且誤差較小。

    （3）K_x的計算

    如圖5，設(shè)單個波紋管的彈性系數(shù)為k′，懸臂波紋管在懸臂端受橫向力Q和橫向位移y之間的關(guān)系^[1，2]為。

    因系統(tǒng)為兩根波紋管并聯(lián)，所以系統(tǒng)的彈性系數(shù)

    3.1.2 方程（1）的求解

    方程（1）的特解表達式為：

    θ（t）=θ_xCOS（ωt-d）      （3）

    振幅：

    

    α為相位角。

    3.2 繞y軸的振動微分方程的推導(dǎo)

    近似的彈簧質(zhì)量系統(tǒng)如圖6所示，圖中F_c為科氏力，F(xiàn)_d為空氣阻力，N為彈性回復(fù)力。

    運動微分方程可寫為：

    J_yφ+C_yφ+K_yφ=T₀sin（ωt-α）     （4）

    式中：J_y--系統(tǒng)相對Y軸的轉(zhuǎn)動慣量；
          C_y--等效粘性阻尼系數(shù)；
          K_y—扭轉(zhuǎn)振動系統(tǒng)的剛性系數(shù)；T0sin（ωt-
          T₀sin（cotω-α）--科氏力產(chǎn)生的扭矩。

    3.2.1 扭矩的計算

    如圖2，直管內(nèi)流體質(zhì)量元dm產(chǎn)生科氏力大小為：

    

    式中ρ為流體密度，A為截面面積，v為流體流速。

    

    半圓管內(nèi)流體質(zhì)量元產(chǎn)生的科氏力大小為：

    

    3.2.2 方程（4）的系數(shù)計算

    （1）J_y的計算

    將流量管仍然分為同樣的三部分計算，

    

    （2）計算K_y

    波紋管仍然受彎，計算方法與繞x軸振動相同，
 
    （3）阻尼系數(shù)C_y的確定

    C_y確定方法與C_x完全相同，不再詳述。

    3.2.3 方程（4）的求解

    方程（4）特解為：
    
    φ（t）=Φsin（ωt-α-β） （6）

    振幅：

    

    β為相位角。

    4 靈敏度的計算

    B點速度：

    Z_B=（L₁+L₂）θ（t）     （8）

    將式（3）、（5）、（6）、（7）、（8）代入得：

    

    合理設(shè)計管子結(jié)構(gòu)和選擇材料、參數(shù)及在弱阻尼情況下，可取α、β為零，式（9）化為：

    

    所以靈敏度：

    

    此方法經(jīng)過實例驗證，證明了其正確性。

    參考文獻

    [1] 樊大鈞，劉廣玉.新型彈性敏感元件設(shè)計[M].北京：國防工業(yè)出版社。1995：372-377.
    [2] 清華大學(xué)工程力學(xué)系固體力學(xué)教研組振動組.機械振動[M].北京：機械工業(yè)出版社，1980：54-57.
    [3] 樊大鈞.波紋管設(shè)計學(xué)[M].北京：北京理工大學(xué)出版社.1995：265-268.
    [4] 范本雋，周一屆.U型科里奧利質(zhì)量流量計的力學(xué)分析與靈敏度計算[J].無錫輕工大學(xué)學(xué)報，1998，17（4）：81-85