基于涡轮流量计的智能燃油计量校正装置的制作方法

1.本发明属于航空燃油计量领域，尤其涉及基于容积式涡轮流量计的智能燃油计量校正装置。


背景技术：

2.发动机燃油与控制系统在装机前应进行大量的地面试验，从半物理模拟试验到模拟飞行试验。在半物理试验中，动态闭环试验要求连续测量计量燃油流量，涡轮流量计因其可测范围宽、响应快、精度相对较高，作为试验中的优选传感器，但是涡轮流量计易受介质的温度、压力、粘度等因素的影响，并且涡轮流量计是容积式流量计而发动机控制系统以质量流量作为燃油流量的控制指标。当燃油密度可知时，由燃油体积流量便可得到质量流量。由于燃油的密度、粘度等属性受温度、压力宽范围的影响会发生较大的变化，这给燃油质量流量的计量带来了困难。


技术实现要素：

3.发明目的：提出一种基于温度、压力校正的燃油质量流量计量校正装置，不仅能够对密度进行温度、压力下的校正，对体积流量进行温度、压力下的校正得到校正补偿后的燃油质量流量，而且通过设置燃油计量活门，以判断补偿后的燃油质量流量的精确度是否满足测量精度要求。
4.技术方案：提供一种基于涡轮流量计的智能燃油计量校正装置，所述装置包括：信号采集电路和微电子处理器；
5.信号采集电路，用于分别采集燃油流量、温度和压力信号及燃油计量活门的开度面积，并将采集的信号通过a/d模块转换后发送给微电子处理器；
6.微电子处理器，用于首先，根据获取的温度、压力计算得到温度压力校正后的体积流量q
v1
，根据接收的燃油计量活门的开度面积a计算得到流经计量活门的燃油流量q
v2
；然后，根据体积流量q
v1
、燃油流量q
v2
的相对误差与设置阈值的差异，判断体积流量q
v1
是否满足要求；当体积流量q
v1
满足要求，将温度压力校正后的密度ρ
op
与体积流量q
v1
的乘积作为真实的质量流量输出；当体积流量q
v1
不满足要求，将燃油流量q
v2
作为真实的质量流量输出。
7.可选地，所述信号采集电路包括温度采集电路、压力采集电路和流量采集电路，对应分别采集燃油温度、压力和流量信号。
8.可选地，所述信号采集电路还包括计量活门开度面积传感器lvdt，用于采集计量活门的开度面积。
9.可选地，所述装置还包括输出电路，输出电路包括比较器u1、u2,电阻r1、r2、r3、r4、r5和r6，二极管d1、d2和d3，三极管t；
10.比较器u1的正端与微电子处理器的d/a输出端连接、负端与输出端连接；比较器u1的输出端与电阻r1的一端连接；电阻r1的另一端与电阻r2的一端、比较器u2的负端连接；比较器u2的正端接地、参考电压接5v电源、输出端接三极管t基极；二级管d1接5v电源，并与电
容c1并联；三极管t集电极接电阻r4的一端、发射极接电阻r5的一端；电阻r5的另一端接地；电阻r3的一端接5v电源，另一端与电阻r4的另一端连接；电阻r6的一端接地，另一端与二级管d2的阳极连接，二级管d2的阴极接电容c2；二级管d3的阳极与电容c2连接，阴极与电阻r3、电阻r4连接。
11.可选地，所述微电子处理器，利用获取的温度、压力值分别对燃油参考密度ρ
ref
进行校正，得到温度对密度的校正因子vcf、压力对密度的校正因子c
pfd
；再根据校正因子vcf、c
pfd
和密度p
ref
的乘积得到校正后的密度ρ
op
。
12.可选地，温度压力校正后的体积流量q
v1
的计算过程为，
13.利用燃油粘度随温度的变化关系，得到温度校正后的燃油粘度ν
top
；利用动态粘度系数μs和c
pfd
，确定粘度的压力校正因子c
pfv
；利用燃油粘度ν
top
、压力校正因子c
pfv
的乘积得到温度、压力校正后的粘度ν
op
；
14.利用ν
op
确定燃油粘度对应的roshko数ro
op
；根据检定装置得到涡轮流量计特性的strouhal-roshko校准曲线，利用校准曲线和roshko数ro
op
确定strouhal数st
op
；利用确定的strouhal数st
op
确定燃油工作状态下真实的k因子值kop；根据涡轮的频率f、k因子值kop得到温度压力校正后的体积流量q
v1
。
15.可选地，燃油流量q
v2
的计算公式为，
[0016][0017]
其中，cd为燃油比热容比；δp＝p
op-p
cal
，p
op
表示燃油工作状态的压力，p
cal
表示校准过程中流量计中燃油的绝对压力值。
[0018]
可选地，设置阈值为0-0.01；体积流量q
v1
、燃油流量q
v2
的相对误差为|(q
v1-q
v2
)
/qv2
|；
[0019]
当体积流量q
v1
、燃油流量q
v2
的相对误差大于设置阈值时，体积流量q
v1
是不满足要求。
[0020]
本发明的技术效果：通过将温度压力校正后的体积流量与计量活门的燃油流量q
v2
进行阈值比较，可以判断涡轮流量计计量是否满足要求；能够实时判断流量计的计量故障。
附图说明
[0021]
图1为校正硬件总体结构图；
[0022]
图2为本发明总体逻辑示意图；
[0023]
图3为流量信号处理电路；
[0024]
图4为温度信号处理电路图
[0025]
图5为质量流量信号输出电路图
[0026]
图6为strouhal-roshko校准曲线
[0027]
图7为校正及质量流量计算程序。
具体实施方式
[0028]
实施例1
[0029]
本实施例，结合图1所示，提供一种基于涡轮流量计的智能燃油计量校正装置，所述装置包括：信号采集电路和微电子处理器；
[0030]
信号采集电路，用于分别采集燃油流量、温度和压力信号及燃油计量活门的开度面积，并将采集的信号通过a/d模块转换后发送给微电子处理器。
[0031]
微电子处理器，用于首先，根据获取的温度、压力计算得到温度压力校正后的体积流量q
v1
，根据接收的燃油计量活门的开度面积a计算得到流经计量活门的燃油流量q
v2
；然后，根据体积流量q
v1
、燃油流量q
v2
的相对误差与设置阈值的差异，判断体积流量q
v1
是否满足要求；当体积流量q
v1
满足要求，将温度压力校正后的密度ρ
op
与体积流量q
v1
的乘积作为真实的质量流量输出；当体积流量q
v1
不满足要求，将燃油流量q
v2
作为真实的质量流量输出。
[0032]
其中，信号采集电路包括温度采集电路、压力采集电路和流量采集电路，对应分别采集燃油温度、压力和流量信号，具体电路连接，如图3、图4。流量采集电路的传感器为涡轮流量计，然后经过电容滤波、整形电路后发送给微电子处理器的a/d端口。
[0033]
本实施例，微电子处理器以stm32单片机为控制核心，外围电路包括流量、温度、压力信号采集电路以及输出电路，再搭配晶振、电源等构成。在进行系统初始化之后，如定时器、gpio串口、adc、dac以及iic、eeprom等初始化后，进行信号的采集，流量信号由涡轮流量计采集，经过模数转换等处理之后传送到单片机；选用pt100铂热电阻传感器构成温度采集电路；压力传感器采用麦克高温压力传感器mpm4528。采集温度、压力信号传送到单片机，为之后的温度、压力校正提供依据。
[0034]
本实施例，结合图3所示，流量信号的采集通过涡轮流量计实现，本发明采用线圈传感器对流量频率进行数据采集，脉冲信号的强弱、干扰都影响涡轮流量计的测量。故在信号送入单片机之前需要对其进行放大、滤波、限幅以及整形等信号处理。其中接头ch1连接涡轮传感器的信号线与地线，pin端连接单片机的定时器,脉冲信号首先经过电容进行滤波处理；然后采用运算放大器进行放大处理，使脉冲信号高电平被限制在运放的输出电压；最后利用比较器整形处理得到高电平为3.3v、低电平为0v的矩形脉冲信号。开启定时器输入捕获模式，单片机就能捕获每一个脉冲，利用周期法测量频率f。具体地，本实施例，附图3中电阻r16的阻值为2k,r17的阻值为1m,r18的阻值为100k，r19的阻值为30k，r20的阻值为10k，r21的阻值为200k，r22的阻值为0；单位均为欧姆。附图3中还包括电容c3，为10uf，及两个比较器u5、u6。
[0035]
本实施例，具体温度采集电路结合图4所示,温度信号的采集使用pt100铂热电阻。温度的变化引起pt100阻值的变化，这种阻值的变化不能被单片机直接检测到，需要通过惠斯通电桥并结合相应算法，将阻值信号转换成相应的电压信号，再经过放大电路之后传输给单片机。具体地，本实施例，附图4中，电阻r7的阻值为10k，r8的阻值为10k，r9的阻值为68，r10的阻值为100k,r11的阻值为100k，r12的阻值为620k，r13的阻值为100k，r14的阻值为620k，r15的阻值为680k；电容c3为1uf，电容c4为1uf，电容c5为10uf；此外还包括比较器u3和u4。
[0036]
本实施例，压力信号的采集采用mpm4528压力传感器。mpm4528高温压力传感器适用于高温介质的压力测量，高温可达300℃。压力量程为0到100kpa,输出信号为4到20ma的电流输出或1到5v的电压输出，可直接送入单片机，为之后的温度、压力校正提供依据。
[0037]
本实施例，单片机选用stm32f103ret6作为控制芯片，其内部包含3个12位的adc、2
个12位dac，4个通用16位定时器，2个16位pwm高级定时器，此外在片内集成了高速存储器，具有丰富的i/o接口，支持i2c、api、usart多种通讯协议。
[0038]
本实施例，使用stm32自带片内ad、da模块；芯片内配置有运行程序，运行程序的整体逻辑，整体逻辑如图2所示。在经过温度、压力校正并求出流体的质量流量之后，输出电路依旧要以0-5v电压信号输出和4-20ma的电流信号输出。由于电流信号抗干扰性较之电压信号更好，通常选择使用电流信号输出。
[0039]
如图5所示，输出电路包括比较器u1、u2,电阻r1、r2、r3、r4、r5和r6，二极管d1、d2和d3，三极管t；
[0040]
比较器u1的正端与微电子处理器的d/a输出端连接、负端与输出端连接；比较器u1的输出端与电阻r1的一端连接；电阻r1的另一端与电阻r2的一端、比较器u2的负端连接；比较器u2的正端接地、参考电压接5v电源、输出端接三极管t基极；二级管d1接5v电源，并与电容c1并联；三极管集电极接电阻r4的一端、发射极接电阻r5的一端；电阻r5的另一端接地；电阻r3的一端接5v电源，另一端与电阻r4的另一端连接；电阻r6的一端接地，另一端与二级管d2的阳极连接，二级管d2的阴极接电容c2；二级管d3的阳极与电容c2连接，阴极与电阻r3、电阻r4连接。具体地，本实施例，电阻r1、r2阻值为51k欧姆，电阻r3为5.1k欧姆，电阻r4为330欧姆，电阻r5为30欧姆，电阻r6为125欧姆；电容c1为10uf、电容c2为100nf；二极管d1为5.1v,二极管d2和d3型号为in4007。
[0041]
本实施例，位移传感器采用阿贝克ht-10系列双余度冗余设计lvdt线性位移传感器。该传感器的工作温度在-55℃到 165℃之间，量程为-5mm到 8mm。
[0042]
进一步地，结合图7所示，本实施例，芯片内的体积流量q
v1
的具体运行逻辑如下所述：
[0043]
(1)求解温度、压力校正后的密度
[0044]
利用获取的温度、压力值分别对燃油参考密度ρ
ref
进行校正，得到温度对密度的校正因子vcf、压力对密度的校正因子c
pfd
；再根据校正因子vcf、c
pfd
和密度p
ref
的乘积得到校正后的密度ρ
op
。具体如下所示：
[0045]
温度对密度的校正因子vcf，
[0046]
其中，为燃油参考温度下的燃油热膨胀系数，δt＝t
op-t
ref
，t
ref
为燃油参考温度、t
0p
为温度测量电路测得的温度。
[0047]
压力对密度的校正因子c
pfd
，
[0048]cpfd
＝(1-f*(p
atm-p
eq
))/(1-f*(p
op-p
eq
))，其中，f为燃油可压缩性校正因子，p
atm
为标准大气压力，p
eq
为平衡压力p
op
表示燃油工作状态的压力，由压力测量电路测得。
[0049]
(2)求解温度、压力校正后的体积流量q
v1
，包括以下内容，
[0050]
a，利用燃油粘度随温度的变化关系，得到温度校正后的燃油粘度ν
top
；利用动态粘度系数μs、c
pfd
，确定粘度的压力校正因子c
pfv
；利用燃油粘度ν
top
、压力校正因子c
pfv
的乘积得到温度、压力校正后的粘度ν
op
；
[0051]
其中，粘度的压力校正因子c
pfv
的计算公式为，
[0052]
其中，p
gage
表示表压力、μs表示动态粘度系
数。
[0053]
燃油粘度ν
top
的计算公式为，
[0054]
其中，中间变量az
top
＝z
top-0.7；z
top
表示燃油工作状态下的粘度；
[0055]
粘度ν
op
的计算公式为，ν
op
＝ν
top
*c
pfv
。
[0056]
b，利用ν
op
确定燃油粘度对应的roshko数ro
op
；根据检定装置得到涡轮流量计特性的strouhal-roshko校准曲线，利用校准曲线和roshko数ro
op
确定strouhal数st
op
；利用确定的strouhal数st
op
确定燃油工作状态下真实的k因子值kop；根据涡轮的频率f、k因子值kop得到温度压力校正后的体积流量q
v1
。具体内容如下所述：
[0057]
1)涡轮流量计校正因子的估计
[0058]
首先，根据燃油工作温度和涡轮膨胀系数确定工作温度下涡轮流量计的温度校正因子c
tr
和c
tk
；具体计算过程如下，
[0059]
根据工作温度和涡轮膨胀系数确定工作温度下涡轮流量计的温度校正因子c
tr
、c
tk
；其中，α1表示涡轮流量计的热膨胀系数，为已知常数；δt＝t
op-t
ref
，t
ref
为燃油参考温度、t
0p
为温度测量电路测得的温度。
[0060]ctr
＝(1 α1*δt)2≈1 2*α1*δt
[0061]ctk
＝(1 α1*δt)3≈1 3*α1*δt
[0062]
然后，再根据指定涡轮流量计的口径d、材料的弹性系数e、流量计的壁厚t分别确定压力校正因子c
pr
和c
pk
。其中，压力变化δp＝p
op-p
cal
，p
op
表示燃油工作状态的压力，p
cal
表示校准过程中流量计中燃油的绝对压力值；
[0063]cpr
＝[1 (δp*d)/(e*t)]2[0064]cpk
＝[1 (δp*d)/(e*t)]3[0065]
2)利用ν
op
确定燃油粘度对应的roshko数ro
op
[0066]
根据校正后粘度ν
op
、涡轮流量计的频率f、温度校正因子c
tr
、压力校正因子c
pr
计算得到roshko数ro
op
；
[0067]
ro
op
＝(f/ν
op
)*c
tr
*c
pr
。
[0068]
3)确定燃油工作状态下真实的k因子值
[0069]
根据检定装置得到涡轮流量计特性的strouhal-roshko校准曲线，如图6所示，利用校准曲线和roshko数ro
op
确定strouhal数st
op
，求取实际工作状态下的k因子值k
op
：
[0070]kop
＝st
op
/(c
tk
*c
pk
)
[0071]
4)根据涡轮的频率f、k因子值k
op
得到温度压力校正后的体积流量q
v1
＝f/k
op
。
[0072]
进一步地，本实施例，芯片内的燃油流量q
v2
的具体计算逻辑为，
[0073][0074]
其中，cd为燃油比热容比；压力变化δp＝p
op-p
cal
，p
op
表示燃油工作状态的压力，p
cal
表示校准过程中流量计中燃油的绝对压力值。
[0075]
本实施例，结合图2所示，芯片内设置阈值为0-0.01；体积流量q
v1
、燃油流量q
v2
的相对误差为|(q
v1-q
v2)/qv2
|。
[0076]
当体积流量q
v1
、燃油流量q
v2
的相对误差大于设置阈值时，体积流量q
v1
是不满足要求，芯片将燃油流量q
v2
作为真实的质量流量通过输出电路输出；当体积流量q
v1
、燃油流量q
v2
的相对误差小于设置阈值时，体积流量q
v1
是满足要求，将燃油流量q
v1
作为真实的质量流量通过输出电路输出。