CN108507752B - 一种计及压力影响的风洞应变天平带空气桥校准方法 - Google Patents

一种计及压力影响的风洞应变天平带空气桥校准方法 Download PDF

Info

Publication number
CN108507752B
CN108507752B CN201810183486.5A CN201810183486A CN108507752B CN 108507752 B CN108507752 B CN 108507752B CN 201810183486 A CN201810183486 A CN 201810183486A CN 108507752 B CN108507752 B CN 108507752B
Authority
CN
China
Prior art keywords
balance
unit
interference coefficient
pressure
combination
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810183486.5A
Other languages
English (en)
Other versions
CN108507752A (zh
Inventor
苗磊
李建强
李耀华
张威权
何成军
谢斌
熊能
张诣
潘华烨
马涛
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Original Assignee
High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center filed Critical High Speed Aerodynamics Research Institute of China Aerodynamics Research and Development Center
Priority to CN201810183486.5A priority Critical patent/CN108507752B/zh
Publication of CN108507752A publication Critical patent/CN108507752A/zh
Application granted granted Critical
Publication of CN108507752B publication Critical patent/CN108507752B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M9/00Aerodynamic testing; Arrangements in or on wind tunnels
    • G01M9/06Measuring arrangements specially adapted for aerodynamic testing
    • G01M9/062Wind tunnel balances; Holding devices combined with measuring arrangements

Landscapes

  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明提供了一种计及压力影响的风洞应变天平带空气桥校准方法,根据不同压力状态下有限数量的天平带空气桥工作公式,拟合生成适用于试验所需压力范围的通用天平工作公式;通用天平工作公式中包含压力参数。在天平带空气桥的通用天平工作公式中,加入了压力参数,风洞试验时,将实际测量的压力值、天平各分量的电信号增量值带入通用公式,计算得到准确的气动载荷,与传统天平公式的计算结果比较,通用公式计算获得的气动载荷更加准确,有利于该类型风洞试验精准度的提升。

Description

一种计及压力影响的风洞应变天平带空气桥校准方法
技术领域
本发明涉及风洞领域,尤其是一种风洞应变天平带空气桥校准方法。
背景技术
风洞应变天平带空气桥在试验所需状态(不同压力)下的校准是风洞试验(试验中利用风洞应变天平带空气桥的测力系统来实现气动载荷的测量)的基础,通过校准以确定天平带空气桥在相应状态下的工作公式。风洞试验时,选取与校准状态基本一致的状态进行试验,为确保试验数据的准确性,天平校准时未涉及的压力状态,风洞试验则不开展。实际中,为多开展不同压力状态的风洞试验,则需增加天平带空气桥的校准状态,这导致天平带空气的校准周期过长。
风洞试验时,为尽可能将压力状态调整至于校准状态一致,通常采用高精度压力控制系统来调节、控制压力,但即使采用这种方法,也无法完全消除校准状态与试验状态之间的差异,这对天平带空气桥工作公式的准确使用造成一定的影响,同时也制约风洞试验质量的提升。
发明内容
本发明的一个目的是提供一种计及压力影响的风洞应变天平带空气桥校准方法,解决天平带空气桥工作公式准确度低以及校准周期长的问题。
具体的,为实现上述目的,本发明提供了一种计及压力影响的风洞应变天平带空气桥校准方法,根据不同压力状态下有限数量的天平带空气桥工作公式,拟合生成适用于试验所需压力范围的通用天平工作公式;
所述通用天平工作公式中包含压力参数。
进一步地,拟合生成通用天平工作公式包括以下步骤:
首先,根据不同压力点下有限数量常规天平公式的主项系数,拟合生成包含压力参数的通用主项系数;
其次,根据不同压力点下有限数量常规天平公式的一次干扰系数,拟合生成包含压力参数的通用一次干扰系数;
最后,根据不同压力点下有限数量常规天平公式的二次干扰系数,生成包含压力参数的通用二次干扰系数。
通用天平工作公式包括主项系数、一次干扰系数和二次干扰系数,依次将所有系数用含压力参数表示,形成最终的通用天平工作公式。
进一步地,有限数量的天平带空气桥工作公式通过以下步骤获得:
步骤A:确定天平带空气桥校准的压力范围;
步骤B:根据压力范围,确定天平带空气桥校准的压力点,压力点在压力范围内均匀分布;
步骤C:根据确定的压力点,从小到大逐步展开相应压力点下的静态校准,并评定天平带空气桥工作公式的准确性。
上述三个步骤可得到不同压力状态下有限数量的天平带空气桥工作公式,用于拟合得到通用天平工作公式。
进一步地,压力范围通过以下公式获得:
Figure GDA0002242481510000021
其中:m为流量,A*为喷管喉道处截面积,P0为总压,T0为气流总温值;当总压P0确定后,天平带空气桥校准的压力范围,即量程也随之确定,该量程等于P0
进一步地,天平带空气桥校准的压力点的数量,是根据天平带空气桥的整体性能确定,若整体性能好,选取用于校准的压力点的数量是量程的2倍,此处量程指压力范围中的最大值,若整体性能不好,压力点的数量是量程的4倍。
选取合适数量的压力点有助于提高获得高精度的天平工作公式。
进一步地,整体性能好坏的判断原则:若天平带空气桥在不充压力校准时,天平各分量中,对应的加载与卸载点回零差值不大于0.03Mv,各分量在正、负方向加载时的主项系数差异不大于0.2%,则天平带空气桥的整体性能好,否则可判定天平带空气桥的整体性能差。
进一步地,压力点在压力范围内按等间距分布。
等间距分布的形式更容易得到各项系数的分布规律,提高校准的准确性。
进一步地,静态校准包括单元校准、交叉校准、工作公式生成以及利用组合载荷对公式的准确性进行验证四个过程。
静态校准采用四种过程组合的方式,提高校准的准确性,进而提高通用公式的适用范围。
进一步地,利用组合载荷对工作公式的准确性进行验证具体为:
对天平带空气桥施加15组标准的组合载荷,获得每个标准组合载荷作用下天平各分量的电信号输出,将每组电信号带入获得天平工作公式计算得到天平各个分量的测量载荷,将计算得到的每组测量载荷与实际施加的每组标准组合载荷进行比较,利用均方根相对误差来评判,误差小于0.4%的分量,其对应的各项系数可用于天平通用公式的拟合,反之,则不能参与拟合,需重新校准与之相关的各项系数。
进一步地,通用天平工作公式拟合生成后,在非选定的压力点条件下,对天平带空气桥施加组合载荷,验证其通用工作公式并确定该天平带空气桥的校准不确定度。
本发明的有益效果为:
(1)利用同一台风洞应变天平及空气桥在不同压力状态下的工作公式,拟合生成一套适用于风洞试验所需压力范围的通用公式;
(2)天平校准人员无需对风洞试验要求压力范围内的每一个点(值)进行校准,可大大缩短天平带空气桥的校准周期,提高了该类型天平带空气桥的研发效率;
(3)在天平带空气桥的通用天平工作公式中,加入了压力参数,风洞试验时,将实际测量的压力值、天平各分量的电信号增量值带入通用公式,计算得到准确的气动载荷,与传统天平公式的计算结果比较,通用公式计算获得的气动载荷更加准确,有利于该类型风洞试验精准度的提升。
(4)通过该方法生成的天平通用公式适用于风洞试验压力范围内的任何一个压力点,风洞试验中,仅需实际测量空气桥附近的压力,并将该测量值带入天平通用天平工作公式,便能准确计算出试验模型的气动载荷,提高了该类型风洞试验的测量精准度。
具体实施方式
下面结合附图和具体实施方式对本发明作进一步详细的说明。
本发明利用天平带空气桥在有限数量的压力条件下校准的工作公式,拟合生成含压力参数的通用天平公式。包括如下步骤:
步骤A:根据风洞试验模型特点,确定天平带空气桥校准的压力范围。
本实例中,根据风洞试验模型所用全部喷管的喉道出口面积(理论设计值)、流量(给定值)及实际测量的喷管出口处环境压力(静压)、喷管中气流的总温值,即为温度值,该总温值为定值,本申请的一个实施例中气流总温值为300K,利用公式计算获得总压和落压比。
Figure GDA0002242481510000051
其中,m为流量,A*为喷管喉道处截面积,P0为总压,T0为气流总温值。
落压比=喷管入口处的总压/喷管出口处的环境压力(静压),本文中静压为94Kpa,总压为P0。当总压P0确定后,天平带空气桥校准的压力范围,即量程也随之确定,该量程等于P0。落压比,指的是喷管进口和出口的总压之比,总压指喷管入口处的压力。
步骤B:根据压力范围,确定天平带空气桥校准的压力点(值),压力点在压力范围内均匀分布。
天平带空气桥校准的压力点(值),是根据天平带空气桥的整体性能确定,若整体性能好,选取用于校准的压力点数是量程的2倍,若整体性能不好,点数是量程的4倍。例如压力最大量程为3MPa,则用于校准的压力点数为6。
整体性能好坏的判断原则:若天平带空气桥在不充压力校准时,天平各分量中,对应的加载与卸载点回零差值不大于0.03Mv,各分量在正、负方向加载时的主项系数差异不大于0.2%,则天平带空气桥的整体性能好,否则可判定天平带空气桥的整体性能差。
压力点在压力范围内按阶梯均匀分布,该阶梯均匀分布的意思是压力点按等差分布也就是等间距分布,例如按照0.5,1,1.5,2,2.5,3(MPa)分布。
步骤C:根据确定的压力点(值),从小到大逐步展开相应压力点(值)下的静态校准,并评定天平带空气桥工作公式的准确性。
本实例中,精确的压力控制是天平静态校准的基础,它将影响天平公式的准确性。因此,天平带空气桥在某一确定状态下进行的静态校准过程中,压力的浮动值与确定压力点(值)的相对误差不大于0.3%。压力的浮动值为压力控制时的偏差,即理论值与实际值之间的差值。通过高精度的压力控制系统,把压力调整至某一确定的压力点(值)后,对天平带空气桥进行静态校准,校准过程包括单元校准、交叉校准、工作公式生成以及利用组合载荷对公式的准确性进行验证四个步骤。
单元校准法是指对天平各分量单独加载,加载时其它分量为零,从而求得校准公式的静校方法。
各分量单元校准时加载与卸载的阶梯总数不少于12个,组合交叉校准加载与卸载的阶梯总数取6,对于天平带空气桥的整体性能不好的情况,交叉校准的阶梯数应适当增加,但总数量不超过12个。该处阶梯总数指的是等间距加载点的数量。交叉校准时,对该组合中的某一个分量预先加载其量程载荷的50%,再对组合中的另一分量按单元加载的方式进行校准。
天平带空气桥工作公式中,各分量主项系数、线性干扰系数的获取是通过单元校准实现,它是正、负方向各系数的平均值;非线性干扰的一次、二次系数根据某一分量正、负方向单元校准数据,按最小二乘法计算获得;交叉系数根据单元校准、交叉校准数据,按照比例缩放、求差、求商的方法获得。
工作公式的准确性是通过对天平带空气桥施加15组标准的组合载荷,获得每个标准组合载荷作用下天平各分量的电信号输出,将每组电信号带入获得天平公式计算得到天平各个分量的测量载荷(力或者力矩,假如Y、Mz、X、Mx、Z、My分量的测量载荷分别为1.8Kg;2.1Kg*m;1.9Kg;1.9Kg*m,0.98Kg;1.99Kg*m),将计算得到的每组测量载荷与实际施加的每组标准组合载荷(假如Y、Mz、X、Mx、Z、My分量的测量载荷分别为2Kg;2Kg*m;2Kg;2Kg*m,1Kg;2Kg*m)进行比较,利用均方根相对误差来评判,误差小于0.4%的分量,其对应的各项系数可用于天平通用公式的拟合,反之,则不能参与拟合,需重新校准与之相关的各项系数。采用相同的方法完成全部确定压力点(值)的校准,获取相应压力状态下的天平工作公式。
步骤D:根据不同压力状态下有限数量的天平带空气桥工作公式,拟合生成适用于试验所需压力范围的通用天平工作公式。
根据有限数量的天平公式拟合生成通用天平公式的拟合内容包括主项系数、一次干扰系数(含线性干扰系数和非线性一次干扰系数)、二次干扰系数(含组合交叉干扰系数和非线性二次干扰系数)的拟合。拟合系数的获取是通过借助Office Excel的图表功能来实现。拟合的通用公式与常规天平公式相比,公式中增加了压力参数,这样该通用天平工作公式可适应于不同压力条件。
首先根据不同压力点(值)下有限数量常规天平公式的主项系数,拟合生成包含压力参数的通用主项系数。其中压力对法向力、俯仰力矩、侧向力、偏航力矩分量主项系数的影响规律呈线性变化;对滚转力矩的影响规律呈二次曲线变化;对轴向力影响很小,通用天平公式中的轴向力主项系数取有限数量常规天平公式主项系数的平均值。
含压力参数的通用天平公式各分量(轴向力除外)主项系数的数学表达式如下(以Y、X的主项系数为例说明):
Figure GDA0002242481510000071
Figure GDA0002242481510000072
其中,A0(常数项)对应通用天平公式中的(U-U0)项;n表示天平带空气桥加压校准的点(次)数。
Figure GDA0002242481510000081
表示压力(P)对Y主项(Y的主项系数)的干扰系数,
Figure GDA0002242481510000082
表示压力平方(P2)对Y主项(Y的主项系数)的干扰系数;KX主项n(常规)表示天平带空气桥第n次加压静态校准时X分量的主项系数;U表示受载后的电压值,U0表示受载前的电压值,单位mV。
其次,根据不同压力点(值)下有限数量常规天平公式的一次干扰系数(含线性干扰系数和非线性一次干扰系数),拟合生成包含压力参数的通用一次干扰系数。先计算不同压力下Mz对天平各分量干扰系数的平均值,再计算不同压力下各干扰系数与平均值的差异,最后计算差异值在天平满量程工作时占被干扰分量满量程的比重,若该比重大于0.1%,则通用公式中含压力参数的一次干扰系数不能取平均值,需采用拟合系数;若该比重小于0.1%,则取各干扰系数的平均值。
1)通用公式中,取平均值的一次干扰系数的数学表达式如下(以Mz对X干扰系数为例):
Figure GDA0002242481510000084
其中,n表示天平带空气桥加压校准的点(次)数,m表示加压校准时的最大压力。
2)通用公式中,一次干扰系数
Figure GDA0002242481510000085
取拟合值的数学表达式如下(以Mz对Y干扰系数为例):
Figure GDA0002242481510000087
表示(Mz*P)的组合对Y的干扰系数,
Figure GDA0002242481510000088
表示(Mz*P2)的组合对于Y的干扰系数。
最后,根据不同压力点(值)下有限数量常规天平公式的二次干扰系数(含组合交叉干扰系数和非线性二次干扰系数),生成包含压力参数通用公式的二次干扰系数。对整体性能好的天平而言,其交叉干扰、非线性干扰系数往往较小,通用公式中组合交叉干扰系数、非线性二次干扰系数取各个压力状态下相应干扰系数的平均值。通用公式中,取平均值的二次扰系数
Figure GDA0002242481510000091
的数学表达式如下(以Y*Mz对Y干扰系数为例):
Figure GDA0002242481510000092
其中,n表示天平带空气桥加压校准的点(次)数,m表示加压校准时的最大压力。
根据上述方法计算得出天平带空气桥在试验所需压力范围内的通用天平工作公式的计算过程如表1所示,该公式与天平常规工作公式相比,增加了压力参数,天平常规工作公式的计算过程如表2所示。
天平带空气桥通用天平工作公式为(以Y为例说明):
Y={119.578035+[KP Y主项*P]+KP^2 Y主项*P^2}*(UY-UY0)+[(KY Y)*Y]+[KY*P Y*(Y*P)]+[KY*P^2 Y*(Y*P^2)]+[KMz Y*Mz]+[(KMz*P Y)*(Mz*P)]+[(KMz*P^2 Y)*(Mz*P^2)]+[(KX Y)*X]+[(KMx Y)*Mx]+[(KMx*P Y)*Mx*P]+[(KMx*P^2 Y)*Mx*P^2]+[(KZ Y)*Z]+[(KZ*P Y)*(Z*P)]+[(KZ*P^2 Y)*(Z*P^2)]+[(KMy Y)*My]+[(KMy*P Y)*(My*P)]+[(KMy*P^2 Y)*(My*P^2)]+[(KYY Y)*Y*Y]+[(KYMz Y)*Y*Mz]+[(KYX Y)*Y*X]+[(KYMx Y)*Y*Mx]+[(KYZ Y)*Y*Z]+[(KYMy Y)*Y*My]+[(KMzMz Y)*Mz*Mz]+[(KMzX Y)*Mz*X]+[(KMzMx Y)*Mz*Mx]+[(KMzZ Y)*Mz*Z)]+[(KMzMy Y)*Mz*My]+[(KXX Y)*X*X]+[(KXMx Y)*X*Mx]+[(KXZ Y)*X*Z]+[(KXMy Y)*X*My]+[(KMxMx Y)*Mx*Mx]+[(KMxZ Y)*Mx*Z]+[(KMxMy Y)*Mx*My]+[(KZZ Y)*Z*Z]+[(KZMy Y)*Z*My]+[(KMyMy Y)*My*My]
天平常规工作公式的使用方法(以Y为例说明):
Y=[KY*(Uy-Uy0)]+[(KY Y)*Y]+[(KMz Y)*Mz]+[(KX Y)*X]+[(KMx Y)*Mx]+[(KZ Y)*Z]+[(KMy Y)*My]+[(KYY Y)*Y*Y]+[(KYMz Y)*Y*Mz]+[(KYX Y)*Y*X]+[(KYMx Y)*Y*Mx]+[(KYZ Y)*Y*Z]+[(KYMy Y)*Y*My]+[(KMzMz Y)*Mz*Mz]+[(KMzX Y)*Mz*X]+[(KMzMx Y)*Mz*Mx]+[(KMzZ Y)*Mz*Z]+[(KMzMy Y)*Mz*My]+[(KXX Y)*X*X]+[(KXMx Y)*X*Mx)]+[(KXZ Y)*X*Z]+[(KXMy Y)*X*My]+[(KMxMx Y)*Mx*Mx]+[(KMxZ Y)*Mx*Z]+[(KMxMy Y)*Mx*My]+[(KZZ Y)*Z*Z]+[(KZMy Y)*Z*My]+[(KMyMy Y)*My*My]
其中,Uy表示天平承受载荷后Y分量的电信号测量值,单位是mV;Uy0表示天平不承受载荷时(承受载荷前)Y分量的电信号测量值,单位是mV;P表示天平系统承受的内压力,单位是MPa;
Y表示升力,单位是Kg;
Mz表示俯仰力矩,单位是Kg*m;
X表示轴向力,单位是Kg;
Mx表示滚转力矩,单位是Kg*m;
Z表示侧向力,单位是Kg;
My表示偏航力矩,单位是Kg*m;
Gy表示天平在Y方向的重量,单位是Kg;Gz表示天平在Z方向的重量,单位是Kg;
Kg是重量单位表示公斤;MPa是压力单位表示兆帕;mV是电压单位表示毫伏;m是长度单位表示米;
KY表示Y的主项系数,单位
Figure GDA0002242481510000101
KY Y表示Y对Y的干扰系数,常数;KMz Y表示Mz对Y的干扰系数,单位
Figure GDA0002242481510000102
KX Y表示X对Y的干扰系数,常数;KMx Y表示Mx对Y的干扰系数,单位
Figure GDA0002242481510000103
KZ Y表示Z对Y的干扰系数,常数;KMy Y表示My对Y的干扰系数,单位
Figure GDA0002242481510000111
KYY Y表示Y与Y的组合对Y的干扰系数,单位
Figure GDA0002242481510000112
KYMz Y表示Y与Mz的组合对Y的干扰系数,单位
Figure GDA0002242481510000113
KYX Y表示Y与X的组合对Y的干扰系数,单位KYMx Y表示Y与Mx的组合对Y的干扰系数,单位
Figure GDA0002242481510000115
KYZ Y表示Y与Z的组合对Y的干扰系数,单位
Figure GDA0002242481510000116
KYMy Y表示Y与My的组合对Y的干扰系数,单位
Figure GDA0002242481510000117
KMzMz Y表示Mz与Mz的组合对Y的干扰系数,单位
Figure GDA0002242481510000118
KMzX Y表示Mz与X的组合对Y的干扰系数,单位
Figure GDA0002242481510000119
KMzMx Y表示Mz与Mx的组合对Y的干扰系数,单位
Figure GDA00022424815100001110
KMzZ Y表示Mz与Z的组合对Y的干扰系数,单位
Figure GDA00022424815100001111
KMzMy Y表示Mz与My的组合对Y的干扰系数,单位
Figure GDA00022424815100001112
KXX Y表示X与X的组合对Y的干扰系数,单位KXMx Y表示X与Mx的组合对Y的干扰系数,单位KXZ Y表示X与Z的组合对Y的干扰系数,单位KXMy Y表示X与My的组合对Y的干扰系数,单位
Figure GDA00022424815100001116
KMxMx Y表示Mx与Mx的组合对Y的干扰系数,单位
Figure GDA00022424815100001117
KMxZ Y表示Mx与Z的组合对Y的干扰系数,单位
Figure GDA00022424815100001118
KMxMy Y表示Mx与My的组合对Y的干扰系数,单位
Figure GDA00022424815100001119
KZZ Y表示Z与Z的组合对Y的干扰系数,单位
Figure GDA00022424815100001120
KZMy Y表示Z与My的组合对Y的干扰系数,单位
Figure GDA00022424815100001121
KMyMy Y表示My与My的组合对Y的干扰系数,单位
Figure GDA00022424815100001122
KP Y主项表示压力(P)对Y主项系数的干扰系数,单位
Figure GDA00022424815100001123
KP^2 Y主项表示压力平方(P2)对Y主项系数的干扰系数,单位KY*P Y表示Y与P的组合对Y的干扰系数,单位
Figure GDA0002242481510000121
KY*P^2 Y表示Y与P2的组合对Y的干扰系数,单位
Figure GDA0002242481510000122
KMz*P Y表示Mz与P的组合对Y的干扰系数,单位
Figure GDA0002242481510000123
KMz*P^2 Y表示Mz与P2的组合对Y的干扰系数,单位
Figure GDA0002242481510000124
KMx*P Y表示Mx与P的组合对Y的干扰系数,单位
Figure GDA0002242481510000125
KMx*P^2 Y表示Mx与P2的组合对Y的干扰系数,单位
Figure GDA0002242481510000126
KZ*P Y表示Z与P的组合对Y的干扰系数,单位
Figure GDA0002242481510000127
KZ*P^2 Y表示Z与P2的组合对Y的干扰系数,单位
Figure GDA0002242481510000128
KMy*P Y表示My与P的组合对Y的干扰系数,单位
Figure GDA0002242481510000129
KMy*P^2 Y表示My与P2的组合对Y的干扰系数,单位
Figure GDA00022424815100001210
表1天平带空气桥在试验压力范围内的通用天平工作公式
Figure GDA0002242481510000131
表2天平带空气桥在试验压力范围内的常规工作公式
Figure GDA0002242481510000141
步骤E:在非选定的压力点(值)条件下,对天平带空气桥施加组合载荷,验证其通用天平工作公式并确定该天平带空气桥的校准不确定度。
通用公式形成后,需检验其中各项系数的准确性。选取一个压力点(值),此压力点(值)非步骤B中确定点(值),利用高精度压力控制系统,将压力精调至确定值后,对天平带空气桥施加15组组合载荷(该15组组合载荷完全模拟天平带空气桥在风洞试验中的受载情况),同时利用高精度数据采集系统记录每组载荷作用下天平各分量的电信号增量,将每组信号增量逐一带入通用天平工作公式,计算得出天平的实测载荷,将15组实测载荷与加载的标准载荷进行比较,采用均方根相对误差的计算方法来确定该天平带空气桥的校准不确定度。当不确定度大于1%时,需重复步骤A到E过程。
均方根相对误差k(计算公式如下)是指天平带空气桥的各个分量在15组组合载荷作用下的均方根误差除以天平各分量的最大载荷。
Figure GDA0002242481510000151
其中,n=15,F计算,i表示:第i点电信号较零点电信号增量带入通用天平工作公式计算得到的计算载荷;F理论,i表示:对天平带空气桥施加15组组合载荷中第i点的天平各个分量的标准载荷(含力和力矩);Fmax表示:天平各个分量的最大载荷(含力和力矩)。
本实例中的技术方案,大大缩短了风洞应变天平带空气桥的校准周期,通过将压力参数拟合到天平工作公式中的方法,提升了风洞试验(利用风洞应变天平带空气桥的测量系统进行测力)的精准度。
尽管已经结合优选的实施例对本发明进行了详细地描述,但是本领域技术人员应当理解的是在不违背本发明精神和实质的情况下,各种修正都是允许的,它们都落入本发明的权利要求的保护范围之中。

Claims (6)

1.一种计及压力影响的风洞应变天平带空气桥校准方法,其特征在于,根据不同压力状态下有限数量的天平带空气桥工作公式,拟合生成适用于试验所需压力范围的通用天平工作公式;
所述通用天平工作公式中包含压力参数;
所述拟合生成通用天平工作公式包括以下步骤:
首先,根据不同压力点下有限数量常规天平公式的主项系数,拟合生成包含压力参数的通用主项系数;
其次,根据不同压力点下有限数量常规天平公式的一次干扰系数,拟合生成包含压力参数的通用一次干扰系数;
最后,根据不同压力点下有限数量常规天平公式的二次干扰系数,生成包含压力参数的通用二次干扰系数;
所述通用天平工作公式为:
Y={119.578035+[KP Y主项*P]+KP^2 Y主项*P^2}*(UY-UY0)+[(KY Y)*Y]+[KY*P Y*(Y*P)]+[KY*P^2 Y*(Y*P^2)]+[KMz Y*Mz]+[(KMz*P Y)*(Mz*P)]+[(KMz*P^2 Y)*(Mz*P^2)]+[(KX Y)*X]+[(KMx Y)*Mx]+[(KMx*P Y)*Mx*P]+[(KMx*P^2 Y)*Mx*P^2]+[(KZ Y)*Z]+[(KZ*P Y)*(Z*P)]+[(KZ*P^2 Y)*(Z*P^2)]+[(KMy Y)*My]+[(KMy*P Y)*(My*P)]+[(KMy*P^2 Y)*(My*P^2)]+[(KYY Y)*Y*Y]+[(KYMz Y)*Y*Mz]+[(KYX Y)*Y*X]+[(KYMx Y)*Y*Mx]+[(KYZ Y)*Y*Z]+[(KYMy Y)*Y*My]+[(KMzMz Y)*Mz*Mz]+[(KMzX Y)*Mz*X]+[(KMzMx Y)*Mz*Mx]+[(KMzZ Y)*Mz*Z)]+[(KMzMy Y)*Mz*My]+[(KXX Y)*X*X]+[(KXMx Y)*X*Mx]+[(KXZ Y)*X*Z]+[(KXMy Y)*X*My]+[(KMxMx Y)*Mx*Mx]+[(KMxZ Y)*Mx*Z]+[(KMxMy Y)*Mx*My]+[(KZZ Y)*Z*Z]+[(KZMy Y)*Z*My]+[(KMyMy Y)*My*My],
其中,UY表示天平承受载荷后Y分量的电信号测量值,单位是mV;UY0表示天平不承受载荷时Y分量的电信号测量值,单位是mV;P表示天平系统承受的内压力,单位是MPa;
Y表示升力,单位是Kg;
Mz表示俯仰力矩,单位是Kg*m;
X表示轴向力,单位是Kg;
Mx表示滚转力矩,单位是Kg*m;
Z表示侧向力,单位是Kg;
My表示偏航力矩,单位是Kg*m;
Kg是重量单位表示公斤;MPa是压力单位表示兆帕;mV是电压单位表示毫伏;m是长度单位表示米;
KY表示Y的主项系数,单位
Figure FDA0002242481500000021
KY Y表示Y对Y的干扰系数,常数;KMz Y表示Mz对Y的干扰系数,单位KX Y表示X对Y的干扰系数,常数;KMx Y表示Mx对Y的干扰系数,单位
Figure FDA0002242481500000023
KZ Y表示Z对Y的干扰系数,常数;KMy Y表示My对Y的干扰系数,单位
Figure FDA0002242481500000024
KYY Y表示Y与Y的组合对Y的干扰系数,单位
Figure FDA0002242481500000025
KYMz Y表示Y与Mz的组合对Y的干扰系数,单位
Figure FDA0002242481500000026
KYX Y表示Y与X的组合对Y的干扰系数,单位
Figure FDA0002242481500000027
KYMx Y表示Y与Mx的组合对Y的干扰系数,单位
Figure FDA0002242481500000028
KYZ Y表示Y与Z的组合对Y的干扰系数,单位KYMy Y表示Y与My的组合对Y的干扰系数,单位
Figure FDA00022424815000000210
KMzMz Y表示Mz与Mz的组合对Y的干扰系数,单位
Figure FDA00022424815000000211
KMzX Y表示Mz与X的组合对Y的干扰系数,单位
Figure FDA00022424815000000212
KMzMx Y表示Mz与Mx的组合对Y的干扰系数,单位
Figure FDA00022424815000000213
KMzZ Y表示Mz与Z的组合对Y的干扰系数,单位
Figure FDA00022424815000000214
KMzMy Y表示Mz与My的组合对Y的干扰系数,单位
Figure FDA00022424815000000215
KXX Y表示X与X的组合对Y的干扰系数,单位KXMx Y表示X与Mx的组合对Y的干扰系数,单位KXZ Y表示X与Z的组合对Y的干扰系数,单位KXMy Y表示X与My的组合对Y的干扰系数,单位
KMxMx Y表示Mx与Mx的组合对Y的干扰系数,单位
Figure FDA00022424815000000321
KMxZ Y表示Mx与Z的组合对Y的干扰系数,单位
Figure FDA00022424815000000320
KMxMy Y表示Mx与My的组合对Y的干扰系数,单位
Figure FDA0002242481500000034
KZZ Y表示Z与Z的组合对Y的干扰系数,单位
Figure FDA0002242481500000035
KZMy Y表示Z与My的组合对Y的干扰系数,单位
Figure FDA0002242481500000036
KMyMy Y表示My与My的组合对Y的干扰系数,单位
Figure FDA0002242481500000037
KP Y主项表示压力(P)对Y主项系数的干扰系数,单位KP^2 Y主项表示压力平方(P2)对Y主项系数的干扰系数,单位
Figure FDA0002242481500000039
KY*P Y表示Y与P的组合对Y的干扰系数,单位KY*P^2 Y表示Y与P2的组合对Y的干扰系数,单位
Figure FDA00022424815000000311
KMz*P Y表示Mz与P的组合对Y的干扰系数,单位
Figure FDA00022424815000000312
KMz*P^2 Y表示Mz与P2的组合对Y的干扰系数,单位
Figure FDA00022424815000000314
KMx*P Y表示Mx与P的组合对Y的干扰系数,单位KMx*P^2 Y表示Mx与P2的组合对Y的干扰系数,单位
Figure FDA00022424815000000315
KZ*P Y表示Z与P的组合对Y的干扰系数,单位
Figure FDA00022424815000000316
KZ*P^2 Y表示Z与P2的组合对Y的干扰系数,单位
Figure FDA00022424815000000317
KMy*P Y表示My与P的组合对Y的干扰系数,单位
Figure FDA00022424815000000319
KMy*P^2 Y表示My与P2的组合对Y的干扰系数,单位
2.根据权利要求1所述的校准方法,其特征在于,所述有限数量的天平带空气桥工作公式通过以下步骤获得:
步骤A:确定天平带空气桥校准的压力范围;
步骤B:根据压力范围,确定天平带空气桥校准的压力点,压力点在压力范围内均匀分布;
步骤C:根据确定的压力点,从小到大逐步展开相应压力点下的静态校准,并评定天平带空气桥工作公式的准确性。
3.根据权利要求2所述的校准方法,其特征在于,所述压力范围通过以下公式获得:
Figure FDA0002242481500000041
其中:m为流量,A*为喷管喉道处截面积,P0为总压,T0为气流总温值;
当总压P0确定后,天平带空气桥校准的压力范围,即量程也随之确定,该量程等于P0
4.根据权利要求2或3所述的校准方法,其特征在于,所述天平带空气桥校准的压力点的数量,是根据天平带空气桥的整体性能确定,若整体性能好,选取用于校准的压力点的数量是量程的2倍,若整体性能不好,压力点的数量是量程的4倍。
5.根据权利要求4所述的校准方法,其特征在于,所述整体性能好坏的判断原则:若天平带空气桥在不充压力校准时,天平各分量中,对应的加载与卸载点回零差值不大于0.03Mv,各分量在正、负方向加载时的主项系数差异不大于0.2%,则天平带空气桥的整体性能好,否则可判定天平带空气桥的整体性能差。
6.根据权利要求2、3或5所述的校准方法,其特征在于,所述压力点在压力范围内按等间距分布。
CN201810183486.5A 2018-03-06 2018-03-06 一种计及压力影响的风洞应变天平带空气桥校准方法 Active CN108507752B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810183486.5A CN108507752B (zh) 2018-03-06 2018-03-06 一种计及压力影响的风洞应变天平带空气桥校准方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810183486.5A CN108507752B (zh) 2018-03-06 2018-03-06 一种计及压力影响的风洞应变天平带空气桥校准方法

Publications (2)

Publication Number Publication Date
CN108507752A CN108507752A (zh) 2018-09-07
CN108507752B true CN108507752B (zh) 2020-01-17

Family

ID=63376140

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810183486.5A Active CN108507752B (zh) 2018-03-06 2018-03-06 一种计及压力影响的风洞应变天平带空气桥校准方法

Country Status (1)

Country Link
CN (1) CN108507752B (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109406091B (zh) * 2018-12-06 2023-09-08 中国航空工业集团公司沈阳空气动力研究所 一种对冲型空气桥天平系统
CN110779679B (zh) * 2019-11-21 2021-01-26 中国空气动力研究与发展中心低速空气动力研究所 一种校准箱压力补偿室对高精度测力天平影响的修正方法
CN111220346B (zh) * 2020-01-21 2021-07-13 中国空气动力研究与发展中心超高速空气动力研究所 一种压电天平校准与使用不确定度的评估方法
CN111397839B (zh) * 2020-04-08 2021-09-07 中国空气动力研究与发展中心超高速空气动力研究所 针对风洞气动力测量的光纤天平温度修正方法及装置
CN111624363B (zh) * 2020-05-27 2021-11-19 中国空气动力研究与发展中心高速空气动力研究所 一种考虑流体压缩性影响的热线风速仪校准方法
CN112284679B (zh) * 2020-09-10 2022-12-16 湖北航天飞行器研究所 一种燃气舵测力五分量天平及分力的计算方法
CN112903235B (zh) * 2021-01-27 2023-04-11 中国空气动力研究与发展中心高速空气动力研究所 一种完全模拟试验状态的推力天平多元校准方法
CN113670566B (zh) * 2021-08-20 2024-02-09 中国航天空气动力技术研究院 一种基于风洞磁悬浮天平的脉动力测量方法
CN113686537B (zh) * 2021-09-13 2023-03-10 中国空气动力研究与发展中心高速空气动力研究所 一种基于音爆试验外式天平的混合式双电桥数据处理方法
CN115615654B (zh) * 2022-11-15 2023-03-10 中国空气动力研究与发展中心低速空气动力研究所 一种回流式空气桥流动影响校准试验装置与方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8308433B2 (en) * 2010-09-30 2012-11-13 General Electric Company System and method for controlling wind turbine blades
CN106644365B (zh) * 2016-12-29 2018-12-21 中国航天空气动力技术研究院 一种低速风洞推力矢量天平校准装置
CN106815428B (zh) * 2017-01-13 2020-05-19 中国空气动力研究与发展中心高速空气动力研究所 一种基于智能优化算法的风洞天平校准数据处理方法

Also Published As

Publication number Publication date
CN108507752A (zh) 2018-09-07

Similar Documents

Publication Publication Date Title
CN108507752B (zh) 一种计及压力影响的风洞应变天平带空气桥校准方法
CN114674520B (zh) 用于测力风洞试验应变天平的灵敏度温度效应修正方法
CN112903235B (zh) 一种完全模拟试验状态的推力天平多元校准方法
CN101451897B (zh) 电磁天平式小力值标准装置
CN108984834B (zh) 一种基于响应面法的机翼可靠性评估系统及方法
Ocokoljić et al. Contemporary frame of measurement and assessment of wind-tunnel flow quality in a low-speed facility
CN110333033B (zh) 一种基于传感器灵敏度差异测力仪的标定方法
CN111220346B (zh) 一种压电天平校准与使用不确定度的评估方法
CN108168774A (zh) 一种空间矢量力校准方法
CN114722332B (zh) 一种风洞天平校准数据处理方法
CN112525481A (zh) 一种六分量高精度微量滚转力矩测量装置及测量方法
CN115290293A (zh) 降低轴向力测量元件零点温度效应的应变天平研制方法
CN115307866A (zh) 风洞天平体轴系弹性角在线校准装置及方法
CN116412995A (zh) 一种分区间确定风洞应变天平测量不确定度的方法
Lynn et al. Thermal and pressure characterization of a wind tunnel force balance using the single vector system
CN113447234A (zh) 一种风洞电子压力扫描阀溯源装置及其溯源方法
CN107587955B (zh) 基于深信度网络的火箭发动机推力偏移量的标定方法
CN113819981A (zh) 液氧煤油发动机试验用煤油流量不确定度评估装置及方法
Liu et al. Design and calibration test of a support force measuring system for hypersonic vehicle aerodynamic measurement
Philipsen et al. Advances in propeller simulation testing at the german-dutch wind tunnels (DNW)
Kammeyer et al. Estimation of the uncertainty in internal balance calibration through comprehensive error propagation
CN110686820A (zh) 用于车轮定位的力和力矩的测量计算方法
CN113686486B (zh) 卫星推进系统压力传感器整星标定方法和系统
Adamov et al. Characteristics of the AT-303 hypersonic wind tunnel. Part 2. Aerodynamics of the HB-2 reference model
Anastasijević et al. Application of semiconductor strain gauges in measurements of dynamic stability derivatives in the T-38 wind tunnel

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant