CN106706208A - 航空发动机矢量推力测量台架静态校准方法 - Google Patents
航空发动机矢量推力测量台架静态校准方法 Download PDFInfo
- Publication number
- CN106706208A CN106706208A CN201611153963.0A CN201611153963A CN106706208A CN 106706208 A CN106706208 A CN 106706208A CN 201611153963 A CN201611153963 A CN 201611153963A CN 106706208 A CN106706208 A CN 106706208A
- Authority
- CN
- China
- Prior art keywords
- load
- calibration
- stand
- loading
- factor
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01L—MEASURING FORCE, STRESS, TORQUE, WORK, MECHANICAL POWER, MECHANICAL EFFICIENCY, OR FLUID PRESSURE
- G01L25/00—Testing or calibrating of apparatus for measuring force, torque, work, mechanical power, or mechanical efficiency
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Testing Of Engines (AREA)
Abstract
一种航空发动机矢量推力测量台架的静态校准方法,包括如下步骤:步骤一:对矢量推力测量台架各分量进行预加载,加载的载荷不小于台架设计量程的80%,并保持30min以上;卸载,检查各测力组件输出状态是否正常,步骤二:按照正交设计原理,设计复合加载表,按加载表重复加载和卸载3遍以上,步骤三:选取验证载荷进行加载,重复3遍以上,选取的验证载荷可以为单分量加载或复合加载,选取时应尽量模拟矢量推力测量台架的真实工作状态,步骤四:数据处理及结果评价,主要包括:1)校准系数求解,2)结果评价。本发明缩短了静态校准周期、提高了校准工作效率,而且在一定程度上提高了矢量推力测量台架的测量精度。
Description
技术领域
本发明属于航空发动机推力测试领域,具体涉及一种基于复合加载模式的航空发动机矢量推力测量台架静态校准方法。
已有技术介绍
矢量推力测量台架是一个用于矢量力测量的平台,其主要功能就是将发动机试验时产生的矢量推力准确测量出来,通过推力台架获取矢量力的各方向分量,并对矢量力进行作用点、作用方向和大小的评估。高精度、可靠性好的矢量推力测量台架的研制与应用,不仅要求有很高的加工精度、严格的加工工艺,而且需要高精度的校准。事实上,航空发动机矢量推力台架的校准工作也是试验台建设过程中的关键步骤之一,它直接影响到矢量推力测量台架最终工作性能的好坏。
目前,航空发动机矢量推力测量台架的静态校准大都采用单分量加载模式,即单独对各个分量进行加载,这种加载模式步骤繁琐,校准效率很低,并且由于这种加载模式与台架实际使用时的情况不完全相同,无法准确模拟矢量推力测量台架的真实工作状态,使得得到的校准公式对于台架特性的反映效果大打折扣。
技术内容
解决的技术问题,提供一种航空发动机矢量推力测量台架的静态校准方法,通过该方法不仅可以缩短了静态校准周期、提高校准工作效率,而且在一定程度上提高了矢量推力测量台架的测量精度。
技术方案,一种航空发动机矢量推力测量台架的静态校准方法,包括如下步骤:
步骤一:对矢量推力测量台架各分量进行预加载,加载的载荷不小于台架设计量程的80%,并保持30min以上;卸载,检查各测力组件输出状态是否正常,此步骤可保证矢量推力测量台架处于良好的工作状态,
步骤二:按照正交设计原理,设计复合加载表,按加载表重复加载和卸载3遍以上,复合加载表的设计实现过程主要包括:首先,由实验人员选择合适的正交表模型;然后,以矢量推力测量台架各个力/力矩分量作为加载设计因素,确定各因素的水平数Mi,并将各个分量的设计量程按等间距均分为Mi个水平,其中下标i表示加载的力或力矩分量的序号,其取值为从1到6,将各加载设计因素及各因素的水平数带入正交表模型,即可得到复合加载表,通过正交设计的加载表大大减少了校准实验过程中各分量的组合数目,可以有效减少实验时间,提高实验效率,
步骤三:选取验证载荷进行加载,重复3遍以上,选取的验证载荷可以为单分量加载或复合加载,选取时应尽量模拟矢量推力测量台架的真实工作状态,验证加载数据只用于校准公式验证,不用于校准系数的求解,
步骤四:数据处理及结果评价,包括:1)校准系数求解,2)结果评价,
1)校准系数求解
校准模型公式是用来描述作用在台架上的载荷值与矢量推力测量台架输出值之间的关系式,根据校准公式中自变量项的阶次,可以将校准模型分为一阶模型、二阶模型和高阶模型,实践经验表明,二阶次以内的校准模型就可以满足对矢量推力测量台架的描述,并达到较高精度,
校准模型公式可以用矩阵形式表示:
[F]=[A][R] (1)
式中,[F]为通过校准设备向台架施加的标准校准载荷构成的矩阵,可表示为:
[F]=[F1,F2,F3,F4,F5,F6]T (2)
[R]为由矢量推力测量台架各测力组件输出结果构成的自变量矩阵。对于一阶模型,包含6项自变量;对于二阶模型,包含27项自变量,即:
[R]一阶=[R1,R2,R3,R4,R5,R6]T (3)
[A]为待求的校准系数矩阵,表示了台架各测力组件的输出与校准载荷之间的映射关系,可以作为校准实验的结果,
若用[Fc]表示一个实验加载点的校准载荷,那么K个加载点的载荷数据可以构成一个6×K的校准载荷矩阵[Fc]6×K;类似地,K个加载点数据可以构造一个m×K的自变量矩阵[Rc]m×K,对于一阶模型,m=6;对于二阶模型,m=27。于是,根据前文建立的校准模型,可以得到:
[Fc]6×K=[A]6×m[Rc]m×K+[E]6×K (5)
式中,[E]6×K表示校准残余误差,根据最小二乘原理,要得到系数矩阵[A]6×K的无偏估计必须令[E]6×K等于[0],即:
[Fc]6×K=[A]6×m[Rc]m×K (6)
采用上述加载模式,可以保证构造的自变量矩阵[Rc]m×K是满秩的且其秩为m,此时矩阵的行列式不等于0,必定有唯一的最小二乘解,得到的校准系数矩阵的解为:
此系数矩阵即为航空发动机矢量推力测量台架的静态校准修正系数。
2)结果评价
首先根据测试数据精度及公式(7)确定的函数关系,计算校准系数的不确定度;根据校准系数误差传递的方法,计算各个加载点的不确定度,该值作为矢量推力测量台架的精度评定指标。
有益效果,提供一种航空发动机矢量推力测量台架的静态校准方法,缩短了静态校准周期、提高校准工作效率,而且在一定程度上提高了矢量推力测量台架的测量精度。
附图说明
图1基于复合加载模式的航空发动机矢量推力测量台架的静态校准方法流程图
具体实施方式
一种航空发动机矢量推力测量台架的静态校准方法,包括:
步骤一:对矢量推力测量台架各分量进行预加载,加载的载荷不小于台架设计量程的80%,并保持30min以上;卸载,检查各测力组件输出状态是否正常,此步骤可保证矢量推力测量台架处于良好的工作状态,
步骤二:按照正交设计原理,设计复合加载表,按加载表重复加载和卸载3遍以上,复合加载表的设计实现过程主要包括:首先,由实验人员选择合适的正交表模型;然后,以矢量推力测量台架各个力/力矩分量作为加载设计因素,确定各因素的水平数Mi,并将各个分量的设计量程按等间距均分为Mi个水平,其中下标i表示加载的力或力矩分量的序号,其取值为从1到6,将各加载设计因素及各因素的水平数带入正交表模型,即可得到复合加载表,通过正交设计的加载表大大减少了校准实验过程中各分量的组合数目,可以有效减少实验时间,提高实验效率,
步骤三:选取验证载荷进行加载,重复3遍以上,选取的验证载荷可以为单分量加载或复合加载,选取时应尽量模拟矢量推力测量台架的真实工作状态,验证加载数据只用于校准公式验证,不用于校准系数的求解,
步骤四:数据处理及结果评价,包括:1)校准系数求解,2)结果评价,
1)校准系数求解
校准模型公式是用来描述作用在台架上的载荷值与矢量推力测量台架输出值之间的关系式,根据校准公式中自变量项的阶次,可以将校准模型分为一阶模型、二阶模型和高阶模型,实践经验表明,二阶次以内的校准模型就可以满足对矢量推力测量台架的描述,并达到较高精度,
校准模型公式可以用矩阵形式表示:
[F]=[A][R] (1)
式中,[F]为通过校准设备向台架施加的标准校准载荷构成的矩阵,可表示为:
[F]=[F1,F2,F3,F4,F5,F6]T (2)
[R]为由矢量推力测量台架各测力组件输出结果构成的自变量矩阵。对于一阶模型,包含6项自变量;对于二阶模型,包含27项自变量,即:
[R]一阶=[R1,R2,R3,R4,R5,R6]T (3)
[A]为待求的校准系数矩阵,表示了台架各测力组件的输出与校准载荷之间的映射关系,可以作为校准实验的结果,
若用[Fc]表示一个实验加载点的校准载荷,那么K个加载点的载荷数据可以构成一个6×K的校准载荷矩阵[Fc]6×K;类似地,K个加载点数据可以构造一个m×K的自变量矩阵[Rc]m×K,对于一阶模型,m=6;对于二阶模型,m=27。于是,根据前文建立的校准模型,可以得到:
[Fc]6×K=[A]6×m[Rc]m×K+[E]6×K (5)
式中,[E]6×K表示校准残余误差,根据最小二乘原理,要得到系数矩阵[A]6×K的无偏估计必须令[E]6×K等于[0],即:
[Fc]6×K=[A]6×m[Rc]m×K (6)
采用上述加载模式,可以保证构造的自变量矩阵[Rc]m×K是满秩的且其秩为m,此时矩阵的行列式不等于0,必定有唯一的最小二乘解,得到的校准系数矩阵的解为:
此系数矩阵即为航空发动机矢量推力测量台架的静态校准修正系数。
2)结果评价
首先根据测试数据精度及公式(7)确定的函数关系,计算校准系数的不确定度;根据校准系数误差传递的方法,计算各个加载点的不确定度,该值作为矢量推力测量台架的精度评定指标。
Claims (2)
1.一种航空发动机矢量推力测量台架的静态校准方法,其特征在于,包括如下步骤:
步骤一:对矢量推力测量台架各分量进行预加载,加载的载荷不小于台架设计量程的80%,并保持30min以上;卸载,检查各测力组件输出状态是否正常,
步骤二:按照正交设计原理,设计复合加载表,按加载表重复加载和卸载3遍以上,复合加载表的设计实现过程主要包括:首先,由实验人员选择合适的正交表模型;然后,以矢量推力测量台架各个力/力矩分量作为加载设计因素,确定各因素的水平数Mi,并将各个分量的设计量程按等间距均分为Mi个水平,其中下标i表示加载的力或力矩分量的序号,其取值为从1到6,将各加载设计因素及各因素的水平数Mi带入正交表模型,即可得到复合加载表,
步骤三:选取验证载荷进行加载,重复3遍以上,选取的验证载荷可以为单分量加载或复合加载,选取时应尽量模拟矢量推力测量台架的真实工作状态,
步骤四:数据处理及结果评价,包括:1)校准系数求解,2)结果评价。
2.如权利要求1所述的航空发动机矢量推力测量台架的静态校准方法,其特征在于,所述的校准系数的求解公式为
其中,为求得的矢量推力测量台架校准系数矩阵的最小二乘解,[Fc]为通过校准设备向台架施加的标准校准载荷构成的矩阵,[Rc]为由矢量推力测量台架各测力组件输出结果构成的自变量矩阵,上标T表示矩阵的转置运算,下标m为选取的自变量数量,对于一阶模型,m=6,对于二阶模型,m=27,下标K表示校准加载点的数量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611153963.0A CN106706208A (zh) | 2016-12-14 | 2016-12-14 | 航空发动机矢量推力测量台架静态校准方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611153963.0A CN106706208A (zh) | 2016-12-14 | 2016-12-14 | 航空发动机矢量推力测量台架静态校准方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106706208A true CN106706208A (zh) | 2017-05-24 |
Family
ID=58937210
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611153963.0A Pending CN106706208A (zh) | 2016-12-14 | 2016-12-14 | 航空发动机矢量推力测量台架静态校准方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106706208A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108168774A (zh) * | 2017-12-27 | 2018-06-15 | 中国航发四川燃气涡轮研究院 | 一种空间矢量力校准方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2312316C1 (ru) * | 2006-02-20 | 2007-12-10 | Федеральное государственное унитарное предприятие "Научно-производственное объединение прикладной механики им. академика М.Ф. Решетнева" | Способ измерения тяги электрореактивных двигателей и устройство для его реализации |
CN103674552A (zh) * | 2013-12-13 | 2014-03-26 | 中国燃气涡轮研究院 | 航空发动机试车台推力加载/校准***及自动控制方法 |
CN105445031A (zh) * | 2015-12-15 | 2016-03-30 | 中国燃气涡轮研究院 | 一种小型航空发动机试车台架推力校准*** |
CN105486441A (zh) * | 2015-12-31 | 2016-04-13 | 西安航天动力试验技术研究所 | 姿控发动机矢量推力测量校准一体化装置及测量方法 |
-
2016
- 2016-12-14 CN CN201611153963.0A patent/CN106706208A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2312316C1 (ru) * | 2006-02-20 | 2007-12-10 | Федеральное государственное унитарное предприятие "Научно-производственное объединение прикладной механики им. академика М.Ф. Решетнева" | Способ измерения тяги электрореактивных двигателей и устройство для его реализации |
CN103674552A (zh) * | 2013-12-13 | 2014-03-26 | 中国燃气涡轮研究院 | 航空发动机试车台推力加载/校准***及自动控制方法 |
CN105445031A (zh) * | 2015-12-15 | 2016-03-30 | 中国燃气涡轮研究院 | 一种小型航空发动机试车台架推力校准*** |
CN105486441A (zh) * | 2015-12-31 | 2016-04-13 | 西安航天动力试验技术研究所 | 姿控发动机矢量推力测量校准一体化装置及测量方法 |
Non-Patent Citations (4)
Title |
---|
幸晓龙等: "矢量喷管六分量测力天平***的研制", 《第九届发动机试验与测试技术学术交流会论文集》 * |
杨欢等: "航空发动机试车台推力测量***的校准", 《现代机械》 * |
湛华海: "航空发动机矢量喷管六分量测力***研制", 《工程科技Ⅱ辑》 * |
贺德馨: "《风洞天平》", 31 May 2001 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108168774A (zh) * | 2017-12-27 | 2018-06-15 | 中国航发四川燃气涡轮研究院 | 一种空间矢量力校准方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR102099614B1 (ko) | 풍력발전기용 복합재 블레이드의 피로 수명 예측 방법 및 장치 | |
US7877229B2 (en) | Computational method for load enhancement factors and apparatus for executing same | |
CN108984834B (zh) | 一种基于响应面法的机翼可靠性评估***及方法 | |
Vidanović et al. | Validation of the CFD code used for determination of aerodynamic characteristics of nonstandard AGARD-B calibration model | |
CN114026572A (zh) | 模拟神经网络中的误差补偿 | |
CN110702364B (zh) | 针对桨尖马赫数影响的高空螺旋桨风洞试验数据修正方法 | |
CN115356121A (zh) | 一种涡轮叶片服役环境下寿命以及剩余寿命损伤评价方法 | |
US8881075B2 (en) | Method for measuring assertion density in a system of verifying integrated circuit design | |
US7647573B2 (en) | Method and device for testing delay paths of an integrated circuit | |
CN113408200B (zh) | 航空装备振动环境分析和预测方法 | |
CN106706208A (zh) | 航空发动机矢量推力测量台架静态校准方法 | |
Saraswat et al. | Fast passivity verification and enforcement via reciprocal systems for interconnects with large order macromodels | |
Jie et al. | Uncertainty evaluation using Monte Carlo method with MATLAB | |
CN101571454B (zh) | 用于虚拟测试单元的设备和方法 | |
CN106886487A (zh) | 用于评价fpga软件可靠性的方法 | |
Senn et al. | Modelling the power and energy consumption of NIOS II softcores on FPGA | |
Becker et al. | A comparative study of gradient reconstruction methods for unstructured meshes with application to turbomachinery flows | |
Kerestes et al. | LES Modeling of High-Lift High-Work LPT Blades: Part II—Validation and Application | |
CN112052523B (zh) | 一种基于载荷特性比对的直升机旋翼动部件强度设计方法 | |
Wang et al. | Scale-adaptive simulation of wind turbines, and its verification with respect to wind tunnel measurements | |
CN115064228A (zh) | 一种粘塑性材料的相场断裂模拟方法、装置及设备 | |
CN114993604A (zh) | 一种基于深度学***静态校测方法 | |
CN104251784A (zh) | 机电产品整机综合应力可靠性加速试验方法 | |
Agarwal | Markovian software reliability model for two types of failures with imperfect debugging rate and generation of errors | |
Mansour et al. | Validation of steady average-passage and mixing plane CFD approaches for the performance prediction of a modern gas turbine multistage axial compressor |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20170524 |