摘要:在非理想流场前提下，异径污水流量计传感器将产生较大的丈量误差。进步传感器权函数分布平均度，有助于进步传感器的非理想流场丈量机能。因此，需要开展传感器权函数分布规律的研究。基于有限元软件COMSOL，分析了4种异径污水流量计传感器权函数的平均度，结果表明:矩形异径传感器的权函数平均。建立矩形异径污水流量计传感器三维模型，研究权函数与矩形段长宽高的分布规律，结果表明:矩形段高度对权函数平均性的影响大，宽度稍小，长度影响小。矩形段的高度和宽度越小，权函数分布越平均，丈量结果受非理想流场的影响越小。
引言
异径污水流量计传感器因为安装空间狭小、前后没有理想直管段，管道内被测流场通常长短理想流场，将导致丈量值与真实值存在较大偏差、影响计量精度。为进步非理想流场的丈量机能，需要研究合适的异径截面外形和尺寸，以进步传感器内权函数分布的平均度。然而海内外相关的研究较少。ShercliffJA和BevirMK等人首次提出和深化了污水流量计的权函数理论。卫开夏等人利用ANSYS有限元软件求解非满管污水流量计的权函数分布。孔令富等人使用MATLAB软件中的PDE工具箱对权函数进行有限元求解。王月明等人基于ANSYS对含有非导电物质时的污水流量计进行有限元分析。李雪菁采用COMSOLMultiphysics有限元软件求解非绝缘管污水流量计的权函数分布。王经卓等人基于COMSOL软件，利用流体像素的方法求解污水流量计权函数的分布。上述文献主要针对圆管污水流量计传感器点电极的二维权函数进行分析，其研究结果与实际三维情况存在偏差。同时，尚未有人针对异径污水流量计传感器三维权函数的分布规律进行研究。因为无可参考的非理想流场丈量工况的理论依据，研发职员无法确定毕竟何种异径截面有助于进步权函数平均度，也无法确定哪一种尺寸有助于进步权函数平均度。针对这一题目，本文从理论上研究权函数与耦合电动势关系，确定进步权函数分布平均度有助于非理想流场丈量。通过COMSOL软件采用电场模拟法，分析4种不同异径截面污水流量计传感器的权函数分布平均度，确定较优的异径截面外形。针对优选异径截面外形的圆电极污水流量计传感器，研究权函数分布平均度与异径段长宽高之间的规律。所得结论为异径污水流量计传感器的丈量管结构尺寸设计提供了一定的参考，也为进步异径污水流量计传感器的非理想流场丈量机能提供了理论依据。
1、污水流量计传感器检测原理
当导电性液体在磁场中作切割磁力线运动时，液体中有感应电流产生。假定液体的电导率δ是平均、各向同性的，则欧姆定律的普遍公式写作

式中j→为电流密度矢量，为通过液体单位面积的电流，A/m2;E→为电场强度矢量，V/m;v→为流体速度，m/s;B→为磁感应强度，T。当激励电流角频率ω不大时，流体中的位移电流完全可以忽略，即

将式(2)带入式(1)得污水流量计传感器的基本丈量方程

式中U为感应电动势，V;2为拉普拉斯算子;为哈密顿算子。
通常借助Green函数G来求解微分方程(3)，G知足Laplace方程

根据传感器的管道外形和电绝缘边界前提，建立了完整形式的污水流量计传感器基本方程

式中V为污水流量计传感器丈量空间;W→为权函数。在直角坐标系(x，y，z)中，式(5)可以转换为

若磁感应强度在传感器有效工作区间内分布平均，则磁感应强度B=By，Bx=Bz=0，式(6)可以化为

当流速为轴向流时，即v=－vz，vx=vy=0;则式(7)表示为

同时，若传感器内的权函数分布平均，Wx=W，则式(8)变为

传感器内的权函数分布平均时，感应电动势大小只与流速积分值成正比，不依靠于流型的分布，有利于非理想流场的精确丈量。
2权函数仿真与分析污水流量计传感器内的流体微元切割磁力线产生感应的电势和电位，相称于一个个微小的“电源”。某一点的权函数应为该点微元作为“电源”所产生的电位梯度与电极间电位差之比。所以，可以采用电场模拟法测定权函数:传感器空间内布满导电液体(一般为水)，在电极处施加一定的电压，便会在导电介质中形成一个电场，测得各点的电场强度，并除以中央点的电场强度，即得到归一化后的权函数值，将其绘制成等值线图便可得到权函数分布图。
2．1仿真方法
基于电场模拟法，选择COMSOLMultiphysics有限元仿真软件求解权函数步骤如下:
1)使用AC/DC模块中的电流应用程序模式，圆管半径为32mm，点电极半径0．4mm，仿真模型为二维模型;
2)电极材质设置为金属铜，导电液体为水，电导率为1×10－4S/m;
3)丈量管具有绝缘衬里，知足电绝缘边界前提n→×j→=0，左右电极分别施加1，－1V的电压;
4)划分四边形网格，为了保证仿真结果的精确度，选择极细化网格;
5)使用稳态求解器进行计算，得到各点处的电场强度，并除以中央点处的电场强度，得到归一化后的权函数值。
2．2结果分析
2．2．1不同异径面的影响
为考察不同异径截面权函数分布的平均性，使用上述方法对圆形、正方形、八边形和矩形异径截面的权函数分布进行定性分析。为了便于对比，设置管道口径为DN100，异径部门截面积为3200mm2。所以，圆形异径面半径为32mm，正方形异径面边长为56．6mm，八边形异径面边长为25．8mm，矩形异径面长宽为80×40mm。仿真结果如图1所示，为了便于对比，权函数等势线大小从0开始，以0．25为步长递增到30。由图1(a)～图1(d)可知，矩形异径截面的权函数等势线间距大，即权函数变化梯度小，权函数分布平均。

为了客观评价不同异径截面内权函数分布的平均程度，采用整体平均度来定量衡量权函数的平均性，设电极截面内每个节点的权函数值为 Wk，相应截面的权函数均匀值为 W0，则电极截面内权函数的整体平均度 Ｒ 为

通过式( 10) 计算得到圆形、正方形、八边形、矩形 4 种
不同异径截面权函数分布的整体平均度分别为 1． 811 2， 1． 996 9，1． 915 0，1． 563 9。
综上所述，矩形异径结构的权函数分布平均，所以，异径污水流量计传感器采用矩形异径的管道结构，该结构权函数分布比较平均，能够减少非理想流场引入的丈量误差。在实际出产实践过程中，权函数分布与矩形段长 L、宽 D、 高 H 有关，因此，开展了矩形异径圆电极污水流量计传感器的三维权函数建模分析，终极得到一种权函数分布比较平均的结构尺寸。
2． 2． 2 三维权函数分布
使用 Pro /E 软件建立三维几何模型，导入 COMSOL 软件进行有限元求解。仿真模型如图 2 所示，电极连线为x 轴，连线中点为坐标原点，流体运动方向为 z 轴，传感器管道口径为 DN100，总长250 mm。异径管部门初始结构尺寸 L = 80 mm，D = 80 mm，H = 40 mm，圆形电极半径为17 mm，伸出绝缘衬里的大间隔为 1． 5 mm。

1. 长度的影响
   

首先分析一定 D × H 前提下，L 变化时传感器内的权函数分布情况。因为传感器异径管部门高度 H 越小信号越强，但压损也越大，因此，H 设置为 30 ～ 50 mm; 异径管宽度 D 越大压损越小，但宽度越大传感器体积也越大，所以， D 设置为 60 ～ 90 mm; 异径管段上下需要放置激励线圈，同时异径段前后需要有一定长度的过渡段来不乱流型，因此， L 设置为 60 ～ 120 mm。一共分析了 6 组 D × H 尺寸的传感器权函数分布随 L 的变化情况，如表 2 所示。因为电极截面内的权函数分布对感应电动势影响大，因 此，利 用式( 10) 计算电极截面 xy 平面内的权函数整体平均度 Ｒ。定义相同 D × H 前提下，权函数平均度随 L 变化的波动率为 ML，如下

计算多组相同 D × H、不同 L 时 xy 平面的权函数整体平均度 Ｒ 及波动率 ML，如表 1 所示。

通过表 1 分析可知，跟着 L 的变化，权函数波动率ML≤ ±2． 5 % ，所以 xy 平面内的权函数整体平均度变化较小，即长度 L 对电极截面内的权函数分布影响很小。
2) 宽度和高度的影响
通过上述分析可知，L 对传感器内的权函数分布影响很小，因此固定设置 L 为 80 mm。然后分析异径段 D，H 同时变化时的权函数分布情况。由上节可知，矩形异径截面的 D 设置为 60 ～ 90 mm，H 设置为 30 ～ 50 mm。为了便于分析三维权函数与 D，H 的变化关系，设置 H 与 D 变化步长都是 10 mm，因此，H 变化范围为 30 ～ 60 mm，即 D = { 60，70， 80，90 mm} ，H = { 30，40，50，60 mm} ，一共 16 组异径污水流量计传感器结构。
分别对上述结构进行有限元分析，根据式( 10) 计算 xy平面内权函数整体平均度 Ｒ，根据式( 11) 计算权函数随 H变化的波动率 MH，随 D 变化的波动率 MD，结果如表 2 所 示。

根据表 1 和表 2 权函数平均度的波动率数值可以看出，MH ＞ MD ＞ ML，所以，矩形段 L，D，H 对于权函数平均度的影响程度是依次增强的，高度 H 对权函数平均度影响大，宽度 D 影响稍小，长度 L 影响很小。且 D 和 H 越小，权函数整体平均度 Ｒ 越小，权函数分布越平均。
为了更加全面地比较权函数在三维空间的分布情况，从上述结构中选取 D × H = { 90 × 30，60 × 30，60 × 60} 三组典型结构，分析其权函数在 xy，xz，yz 三个平面内的分布情况。为了便于对比，同一划定三个平面内等势线的分布步长和数值范围: 1) xy，xz 平面内的权函数等势线大小以0． 25为步长，从 0 增加到 30; 2) 因为 yz 平面的权函数小于1，划定 yz 平面内的权函数等势线大小以 0． 05 为步长，从 0增加到 1。详细如图 3 ～ 图 5 所示

通过对图 3 ～ 图 5 分析得出以下结论: 1) 图 3( a) 、图 4( a) 的 xy 面权函数分布表明，D = 90 mm时中央区域的权函数等势线间距较大，即权函数变化梯度较小，且中央区域的权函数等势线逐渐变为直线，因此中央区域的权函数分布更平均; 但 D = 90 mm 时，电极四周的权函数等势线较密，且等势线颜色较深，权函数大值较大，变化梯度较大，所以，电极四周的权函数分布平均性较差。由于难以直接衡量 D 改变时，xy 面权函数分布的平均性。所以，需要利用权函数整体平均度 Ｒ 定量确定 xy 平面内权函数分布的平均性。结果表明，跟着宽度 D 的减小，权函数分布越来越平均。
2) 图 4( a) 、图 5( a) 的 xy 面权函数分布表明，H = 30 mm时中央区域的权函数等势线间距较大，且中央区域的权函数等势线逐渐变为直线; 电极四周的权函数等势线比较稀疏，且等势线颜色较浅，权函数大值较小，变化梯度小，因 此，H = 30 mm 时 xy 面的权函数分布更加平均。

3. 图 3( b) ～ 5( b) 的 xz 面权函数分布表明，三组异径结构的权函数分布情况类似，没有显著的区别，即 D 和 H的变化对 xz 面的权函数分布影响较小。
4. 图 3( c) 、图 4( c) 的 yz 面权函数分布表明，D = 90 mm时的权函数等势线间距略大于 D = 60 mm 时的权函数等势线间距，权函数变化梯度较小，且中央区域的权函数等势线逐渐变为直线，因此 D = 90 mm 的权函数分布更平均一些，但是两者区别很小，即宽度改变对 yz 面权函数分布影响很小。图 4( c) 、图 5( c) 的 yz 面权函数分布表明，H = 30 mm时的权函数等势线间距较大，且中央区域的权函数等势线逐渐变为直线，因此 H = 30 mm 的权函数分布更加平均;
5. 图 3( a) ～ 5( a) 和图 3( b) ～ 5( b) 权函数分布结果表明，越靠近电极，等势线颜色越深，即权函数值越大，且越靠近电极，权函数等势线越密集，即权函数变化梯度越大。
   

综上所述，异径电磁水表异径段长度 L 对权函数分布的平均性影响很小，跟着长度 L 的改变，权函数分布基本没有变化; 异径段高度 H 对权函数分布的平均性影响大，宽度 D 影响稍小，高度和宽度越小，权函数分布越平均，即异径污水流量计传感器的丈量精确度受非理想流场的影响越小。
3 结 论

1. 圆形、正方形、八边形和矩形等 4 种异径污水流量计传感器的权函数分析结果表明，矩形异径截面传感器的权函数分布平均。
2. 电极四周区域，权函数值较大，且权函数变化梯度较大，跟着阔别电极，权函数值越来越小，且权函数变化梯度越来越小。
3. 矩形段高度 H 对权函数分布的平均性影响大，随 着 H 的减小，权函数分布越来越平均，且 y 轴权函数的分布逐渐趋近于常数 1。矩形段宽度 D 对权函数分布的平均性影响稍小，跟着 D 的减小，权函数分布越来越平均。矩形段长度 L 对传感器内的权函数分布影响很小，跟着 L 的改变，权函数分布没有显著变化。
   文章来源：http://www.sipaiyibiao.com/article/xinwenzhongxin/list_3_2.html