发明名称 一种发动机流量计的标定方法
摘要 本发明提供了一种发动机流量计的标定方法,属于流量测量技术领域。其步骤如下:一、建立与使用环境一样的标定条件,二、建立发动机流量测量系统,三、利用槽道式流量计获得标准的发动机流量,四、发动机流量计标定数据处理,五、发动机流量计标定系统不确定度分析。本发明可以解决发动机流量计的标定问题,给出标定系数及不确定度。主要用于发动机流量计的标定。
申请公布号 CN105865587B 申请公布日期 2017.03.15
申请号 CN201610326352.5 申请日期 2016.05.17
申请人 中国人民解放军63820部队吸气式高超声速技术研究中心 发明人 邓维鑫;杨顺华;刘伟雄;杨样;晏至辉;李宏斌;任虎;肖保国;王西耀;张穹洲;田野;邢建文;张顺平;李季
分类号 G01F25/00(2006.01)I 主分类号 G01F25/00(2006.01)I
代理机构 成都点睛专利代理事务所(普通合伙) 51232 代理人 葛启函
主权项 一种发动机流量计的标定方法,其步骤如下一、建立与使用环境相同的标定条件发动机一般工作在来流总温1000K以上的使用环境,故利用风洞来建立与使用环境相同的标定条件,该风洞采用氢气和富氧空气燃烧加热来得到所需的试验来流,试验来流的模拟参数包括:总温、总压及马赫数,确定富氧管配气参数和富氧喉道直径、挤压供氢压力和氢气喉道直径及喷管,即可确定试验来流的条件参数;二、建立发动机流量测量系统将两个槽道式流量计分别放置于富氧空气喉道及氢气喉道的下游来测量富氧空气及氢气的流量,两者之和就是总的来流流量;发动机流量计则放置于燃烧加热器下游,在流量计出口的测量截面测量总压、总温及静压;三、利用槽道式流量计获得标准的发动机流量四、发动机流量计标定数据处理a.计算测量截面的热力学参数计算方法为:<maths num="0001"><math><![CDATA[<mrow><mfrac><msub><mi>C</mi><mrow><mi>p</mi><mi>i</mi></mrow></msub><mi>R</mi></mfrac><mo>=</mo><msub><mi>A</mi><mn>1</mn></msub><mo>+</mo><msub><mi>A</mi><mn>2</mn></msub><mi>T</mi><mo>+</mo><msub><mi>A</mi><mn>3</mn></msub><msup><mi>T</mi><mn>2</mn></msup><mo>+</mo><msub><mi>A</mi><mn>4</mn></msub><msup><mi>T</mi><mn>3</mn></msup><mo>+</mo><msub><mi>A</mi><mn>5</mn></msub><msup><mi>T</mi><mn>4</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001185880170000011.GIF" wi="1404" he="119" /></maths><maths num="0002"><math><![CDATA[<mrow><msub><mi>C</mi><mi>p</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>s</mi></munderover><msub><mi>Y</mi><mi>i</mi></msub><msub><mi>C</mi><mrow><mi>p</mi><mi>i</mi></mrow></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001185880170000012.GIF" wi="1405" he="119" /></maths><maths num="0003"><math><![CDATA[<mrow><msub><mi>s</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mn>1</mn></msub><mo>)</mo></mrow><mo>-</mo><msub><mi>s</mi><mi>i</mi></msub><mrow><mo>(</mo><msub><mi>T</mi><mn>2</mn></msub><mo>)</mo></mrow><mo>=</mo><msubsup><mo>&Integral;</mo><msub><mi>T</mi><mn>1</mn></msub><msub><mi>T</mi><mn>2</mn></msub></msubsup><mfrac><mrow><msub><mi>C</mi><mrow><mi>p</mi><mi>i</mi></mrow></msub><mi>d</mi><mi>T</mi></mrow><mi>T</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001185880170000013.GIF" wi="1430" he="126" /></maths>S(T<sub>0</sub>)‑Rln(P<sub>0</sub>)=S(T)‑Rln(P)                           (4)<maths num="0004"><math><![CDATA[<mrow><mi>&gamma;</mi><mo>=</mo><mfrac><msub><mi>C</mi><mi>p</mi></msub><mrow><msub><mi>C</mi><mi>p</mi></msub><mo>-</mo><mi>R</mi></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001185880170000014.GIF" wi="1373" he="141" /></maths><maths num="0005"><math><![CDATA[<mrow><mfrac><msub><mi>h</mi><mi>i</mi></msub><mrow><mi>R</mi><mi>T</mi></mrow></mfrac><mo>=</mo><msub><mi>A</mi><mn>1</mn></msub><mo>+</mo><mfrac><mrow><msub><mi>A</mi><mn>2</mn></msub><mi>T</mi></mrow><mn>2</mn></mfrac><mo>+</mo><mfrac><mrow><msub><mi>A</mi><mn>3</mn></msub><msup><mi>T</mi><mn>2</mn></msup></mrow><mn>3</mn></mfrac><mo>+</mo><mfrac><mrow><msub><mi>A</mi><mn>4</mn></msub><msup><mi>T</mi><mn>3</mn></msup></mrow><mn>4</mn></mfrac><mo>+</mo><mfrac><mrow><msub><mi>A</mi><mn>5</mn></msub><msup><mi>T</mi><mn>4</mn></msup></mrow><mn>5</mn></mfrac><mo>+</mo><mfrac><msub><mi>A</mi><mn>6</mn></msub><mi>T</mi></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>6</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001185880170000015.GIF" wi="1493" he="119" /></maths><maths num="0006"><math><![CDATA[<mrow><mi>h</mi><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>s</mi></munderover><msub><mi>Y</mi><mi>i</mi></msub><msub><mi>h</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001185880170000016.GIF" wi="1350" he="126" /></maths>h(T<sub>0</sub>)=h(T)+u<sup>2</sup>/2                                  (8)<maths num="0007"><math><![CDATA[<mrow><mi>u</mi><mo>=</mo><mi>M</mi><mi>a</mi><mo>&times;</mo><msqrt><mrow><mi>&gamma;</mi><mi>R</mi><mi>T</mi></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>9</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001185880170000017.GIF" wi="1350" he="77" /></maths>式中:C<sub>p</sub>为定压比热(J/kg·K),T为静温(K),T<sub>0</sub>为总温(K),s为熵(J/kg·K),P为静压(Pa),P<sub>0</sub>为总压(Pa),h为焓(J/kg),u为流速(m/s),Ma为马赫数,Y为气体质量组份;由式(1)、(2)、(3)及(4)迭代即可求出测量截面静温T,由(5)可求得γ,由(6),(7)可求得h(T<sub>0</sub>)及h(T),由(8)可求得流速u,由(9)可求得Ma;b.利用发动机流量计测量气体流量进入发动机流量计的气体流量由安装在音速喷嘴后部的八个总压耙测得,八个总压耙上共装有40个总压管,这40个总压管将音速喷嘴出口截面积平均分为40等份,由式(10)分别计算每一份的流量,最后由式(11)计算总的流量;<maths num="0008"><math><![CDATA[<mrow><msub><mi>&Delta;m</mi><mi>i</mi></msub><mo>=</mo><msub><mi>K</mi><mn>2</mn></msub><mfrac><msubsup><mi>p</mi><mi>i</mi><mo>*</mo></msubsup><msqrt><msub><mi>T</mi><mn>0</mn></msub></msqrt></mfrac><mi>&Delta;</mi><mi>A</mi><mi>q</mi><msub><mrow><mo>(</mo><mi>M</mi><mi>a</mi><mo>)</mo></mrow><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>10</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001185880170000021.GIF" wi="1374" he="150" /></maths><maths num="0009"><math><![CDATA[<mrow><mi>m</mi><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>40</mn></munderover><msub><mi>&Delta;m</mi><mi>i</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>11</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001185880170000022.GIF" wi="1381" he="134" /></maths>式中:p<sub>i</sub>为第i个总压管前总压(Pa),ΔA为音速喷嘴出口截面积的1/40(m<sup>2</sup>),q(Ma)为ΔA面积内的流量函数;c.计算发动机流量计流量系数测得了槽道流量计的流量及待标定发动机流量计的流量之后,可以通过式(12)计算发动机流量计的流量系数,q<sub>m</sub>=K<sub>3</sub>·m                                    (12)式中:K<sub>3</sub>为发动机流量计流量系数;五、发动机流量计标定系统不确定度分析a.不确定度分析的数学模型对于数学模型,被测量Y的最佳估计值y,在通过输入量X<sub>1</sub>,X<sub>2</sub>,…,X<sub>N</sub>的估计值x<sub>1</sub>,x<sub>2</sub>,…,x<sub>N</sub>得出时,可有以下方法:<img file="FDA0001185880170000023.GIF" wi="438" he="87" />式中:<img file="FDA0001185880170000024.GIF" wi="269" he="127" />它是独立观测值x<sub>ik</sub>的算术平均值,在X<sub>i</sub>=x<sub>i</sub>时,定义灵敏系数为:<img file="FDA0001185880170000025.GIF" wi="246" he="63" />在流量计量中常使用相对不确定度,此时采用相对灵敏系数c<sub>ri</sub>,它定义为:<img file="FDA0001185880170000026.GIF" wi="434" he="71" />当全部数量是彼此独立或互不相关时,合成标准不确定度由式(14)和由式(15)得到:<maths num="0010"><math><![CDATA[<mrow><msubsup><mi>u</mi><mi>c</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>u</mi><mi>i</mi><mn>2</mn></msubsup><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>c</mi><mi>i</mi><mn>2</mn></msubsup><msup><mi>u</mi><mn>2</mn></msup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>14</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001185880170000027.GIF" wi="1518" he="127" /></maths><maths num="0011"><math><![CDATA[<mrow><msubsup><mi>u</mi><mrow><mi>c</mi><mi>r</mi></mrow><mn>2</mn></msubsup><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>u</mi><mrow><mi>r</mi><mi>i</mi></mrow><mn>2</mn></msubsup><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msubsup><mi>c</mi><mrow><mi>r</mi><mi>i</mi></mrow><mn>2</mn></msubsup><msubsup><mi>u</mi><mi>r</mi><mn>2</mn></msubsup><mrow><mo>(</mo><msub><mi>x</mi><mi>i</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>15</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001185880170000028.GIF" wi="1534" he="134" /></maths>b.确定发动机流量计标定系统的合成不确定度发动机流量计标定系统可能输入不确定度的装置有槽道流量计、压力传感器、总温传感器及音速喷嘴,输入各个装置的不确定度即可以根据式(14)及式(15)计算系统的合成不确定度。
地址 621000 四川省绵阳市涪城区二环路南段6号