CN110398339B - 一种满足多级全自由飞风洞试验相似律的条件确定方法 - Google Patents

一种满足多级全自由飞风洞试验相似律的条件确定方法 Download PDF

Info

Publication number
CN110398339B
CN110398339B CN201910576293.0A CN201910576293A CN110398339B CN 110398339 B CN110398339 B CN 110398339B CN 201910576293 A CN201910576293 A CN 201910576293A CN 110398339 B CN110398339 B CN 110398339B
Authority
CN
China
Prior art keywords
wind tunnel
tunnel test
representing
real aircraft
free flight
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
CN201910576293.0A
Other languages
English (en)
Other versions
CN110398339A (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.)
China Academy of Aerospace Aerodynamics CAAA
Original Assignee
China Academy of Aerospace Aerodynamics CAAA
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 China Academy of Aerospace Aerodynamics CAAA filed Critical China Academy of Aerospace Aerodynamics CAAA
Priority to CN201910576293.0A priority Critical patent/CN110398339B/zh
Publication of CN110398339A publication Critical patent/CN110398339A/zh
Application granted granted Critical
Publication of CN110398339B publication Critical patent/CN110398339B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/08Aerodynamic models

Landscapes

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

Abstract

本发明公开了一种满足多级全自由飞风洞试验相似律的条件确定方法。根据动力学及运动学推导投放类风洞自由飞试验相似律。首次针对升力体外形针对两级甚至多级模型完全自由的自由飞试验设计试验相似律,重点解决升力体外形相似性问题。在风洞试验模型满足求限制条件时,可保证风洞试验所获得的水平线位移和竖直线位移、角位移曲线与真实飞行器的曲线相似,保证风洞试验可模拟真实飞行器的分离情况。解决了两级甚至多级全自由飞试验无相似律指导的难题。缺少针对两个甚至多个飞行器完全自由的风洞自由飞试验相似律,本发明针对此方面,得到了满足多级全自由飞风洞试验相似律的条件。

Description

一种满足多级全自由飞风洞试验相似律的条件确定方法
技术领域
本发明涉及自由飞试验参数设计,尤其涉及一种满足多级全自由飞风洞试验相似律的条件确定方法,属于航空航天工程领域。
背景技术
两级入轨飞行器在进行级间分离时处在高Ma、大动压环境中,分离时存在复杂的气动干扰现象,使得分离安全具有不确定性。为保证分离安全,需针对两级入轨飞行器的分离特性进行研究。近几年多体分离全自由飞风洞试验是风洞自由飞领域的研究热点,特别是近几年多体分离全自由飞试验技术有了相当大的进步,虽然试验技术进步很大,但针对多体分离全自由飞风洞试验的相似律尚未建立,需要进行针对性研究。
两级甚至多级都是完全自由的,需要的相似条件较多。级间分离风洞自由飞试验模型相似律设计的主要难点在于:首先,风洞试验与真实飞行不仅要满足水平线位移和竖直线位移的同时相似,还需要满足角位移的相似,更有甚者需要满足上述3项相似的两两相似,这就导致了相似方程较多。其次,升力体外形下两级各自所受气动力与其重力相比不再是小量,在推导相似条件时需重点考虑。
若按照以前试验方法随意设计两级风洞试验模型参数,则会导致风洞试验所得到的两级间水平线位移和竖直线位移以及角位移与真实飞行器的三者不对应,试验误差大,对真实飞行器产生错误判断及指导。
发明内容
本发明解决的技术问题为:克服现有技术不足,提供一种满足多级全自由飞风洞试验相似律的条件确定方法,首次针对升力体外形针对两级甚至多级模型完全自由的自由飞试验设计试验相似律,重点解决升力体外形相似性问题。解决了以往风洞试验与真实飞行器之间水平线位移、竖直线位移以及角位移不对应的问题,提高了风洞试验精度,对真实飞行器研制具有指导作用。
本发明解决的技术方案为:一种满足多级全自由飞风洞试验相似律的条件确定方法,步骤如下:
(1)设定多级全自由飞相似律风洞试验中在tm时刻一级质心坐标为(x1m,y1m)和二级质心坐标为(x2m,y2m);tm表示多级全自由飞相似律风洞试验中一级与二级分离过程中的某个时间点;
Figure BDA0002112184040000021
Figure BDA0002112184040000022
Figure BDA0002112184040000023
Figure BDA0002112184040000024
a1xm表示一级风洞试验模型水平加速度,a1ym表示一级风洞试验模型竖直加速度,a2xm表示二级风洞试验模型水平加速度,a2ym表示二级风洞试验模型竖直加速度,g为重力加速度,v0m表示一级二级共同的运动速度;
其中,a1ym=q∞m·S1m·C1ym/m1m (13)
a1xm=q∞m·S1m·C1xm/m1m (14)
a2ym=q∞m·S2m·C2ym/m2m (15)
a2xm=q∞m·S2m·C2xm/m2m (16)
S1m表示一级风洞试验模型的参考面积(一级风洞试验模型投影到XZ面的投影面积),C1ym表示一级风洞试验模型升力系数,C1xm表示一级风洞试验模型阻力系数,m1m表示一级风洞试验模型质量,S2m表示二级风洞试验模型的参考面积,C2ym表示二级风洞试验模型升力系数,C2xm表示二级风洞试验模型阻力系数,m2m表示二级风洞试验模型质量;
(2)根据步骤(1)多级全自由飞相似律风洞试验中在tm时刻一级质心坐标(x1m,y1m)和二级质心坐标(x2m,y2m),确定风洞试验中tm时刻一级与二级的相对位置参数Cm
Figure BDA0002112184040000031
(3)设一级与二级之间有如下关系:m1m=N*m2m,S1m=M*S2m,当真实飞行器一级和二级质量、外形确定之后,M和N便为常数,则将公式(9)化简得到公式(10)如下:
Figure BDA0002112184040000032
(4)在步骤(1)的同时,设定真实飞行器中,在ts时刻一级质心坐标(x1s,y1s)、二级质心坐标(x2s,y2s);
ts表示真实飞行器一级与二级分离过程中的某个时间点;
Figure BDA0002112184040000033
Figure BDA0002112184040000034
Figure BDA0002112184040000035
Figure BDA0002112184040000036
a1xs表示一级真实飞行器水平加速度,a1ys表示一级真实飞行器竖直加速度,a2xs表示二级真实飞行器水平加速度,a2ys表示二级真实飞行器竖直加速度,g为重力加速度,v0s表示一级二级共同的运动速度;
其中:
a1ys=q∞m·S1s·C1ys/m1s (15)
a1xs=q∞s·S1s·C1xs/m1s (16)
a2ys=q∞s·S2s·C2ys/m2s (17)
a2xs=q∞s·S2s·C2xs/m2s (18)
S1s表示一级真实飞行器的参考面积,C1ys表示一级真实飞行器升力系数,C1xs表示一级真实飞行器阻力系数,s1s表示一级真实飞行器质量,S2s表示二级真实飞行器的参考面积,C2ys表示二级真实飞行器升力系数,C2xs表示二级真实飞行器阻力系数,s2s表示二级真实飞行器质量;
(5)根据步骤(4)设定真实飞行器中,在ts时刻一级质心坐标(x1s,y1s)、二级质心坐标(x2s,y2s),确定真实飞行器ts时刻一级与二级的相对位置参数Cs
设一级二级之间同样有如下关系:m1=N*m2,S1=M*S2,M和N为定值,则:
Figure BDA0002112184040000041
(6)根据公式(10),即公式(20):
Figure BDA0002112184040000042
(7)根据公式(19)和公式(10),则有在满足多级全自由飞相似律风洞试验与真实飞行器的几何相似(几何相似是指:对应边成比例、对应角相等)、飞行条件相同、飞行姿态角一致的情况下:
C2ys=C2ym (21)
C1ys=C1ym (22)
C2xs=C2xm (23)
C1xs=C1xm (24)
因而,有恒等式:
Cs=Cm (25)
(m1m/m2m=m1s/m2s)
Cs=Cm即表示两级模型完全自由状态下,只要保证一级与二级的几何相似和质量比相同,风洞自由飞试验与真实飞行器的线位移轨迹是恒相似的,与风洞自由飞试验中的飞行器模型的质量、来流条件无关;即满足风洞自由飞试验与真实飞行器的线位移轨迹相似的条件;
(8)设定风洞自由飞试验与真实飞行器在对应时刻的水平位移成比例,则有:
Figure BDA0002112184040000051
Figure BDA0002112184040000052
根据(26-1)和(26-2),确定公式(27)
Figure BDA0002112184040000053
化简公式(27),求得:
Figure BDA0002112184040000054
式中,kt表示时间缩比,q表示来流动压,kq∞表示来流动压比,kl表示尺寸缩比,km表示质量缩比;
(9)设定风洞自由飞试验与真实飞行器的角位移轨迹是恒相似的,则需满足以下三个方程:
Figure BDA0002112184040000061
kθ表示角度缩比(角度是指飞行器的姿态角),
Figure BDA0002112184040000062
表示角速度缩比,
Figure BDA0002112184040000063
表示角加速度缩比;
根据物体角加速度定义,则有:
Figure BDA0002112184040000064
因此,根据公式(29)和公式(30),则角位移相等需满足:
Figure BDA0002112184040000065
然而已求得关系式(28),因此满足关系式(31),即满足关系式(29),满足角位移相似;将公式(28)代入公式(29),则有:
Figure BDA0002112184040000066
Figure BDA0002112184040000067
根据公式kt=kl/kv (34)
再结合公式(28),则有:
Figure BDA0002112184040000068
公式(35)即满足风洞自由飞试验与真实飞行器的角位移轨迹相似的条件。
多级全自由飞风洞试验,具体为:为了方便说明,以两级为例介绍多级的概念;通常来说,脱离飞行主体的为一级,继续保持飞行的为二级;为了模拟真实飞行,需进行风洞试验,风洞试验开始时,一级二级作为刚性整体无支撑的在风洞中飞行,之间无相对自由度,在需要进行分离的时刻,突然释放一级二级之间的6自由度限制,在气动力的作用下,两级姿态产生变化,模拟真实飞行器的分离过程。
多级全自由飞风洞试验与真实飞行器的几何相似,是指:风洞试验模型与真实飞行器之间对应边长度成比例,对应角度相等。
多级全自由飞风洞试验与真实飞行器的飞行条件相同,是指:风洞试验Ma数与真实飞行器飞行Ma数相等。
多级全自由飞风洞试验与真实飞行器的飞行姿态角一致,是指:风洞试验模型运动的角位移与真实飞行器对应点的角位移相等。
一级与二级的几何相似,是指:风洞试验模型与真实飞行器的一级与二级之间对应边长度成比例,对应角度相等。
一级与二级的质量比相同,是指:真实飞行器一级二级之间的质量比要与风洞试验一级二级之间的质量比相等。
一级质心坐标(x1m,y1m)和二级质心坐标(x2m,y2m),所在的坐标系为:在两级进行分离前,一级质心坐标和二级质心坐标的同在xy平面内,且二者连线与y轴平行。
建立三维坐标系,气流方向为x正方向,竖直向上为y正方向,右手定则确定z正方向,xy平面穿过每级的质心,也就是说坐标原点及各级质心同时位于xy平面上;下标1表示一级参数,下标2表示二级参数,下标m表示风洞模型参数,下标s表示真实飞行器参数,因此x1m表示一级风洞试验模型质心x坐标,y1m表示一级风洞试验模型质心y坐标,x2m表示二级风洞试验模型质心x坐标,y2m表示二级风洞试验模型质心y坐标。
气流方向为x正方向,竖直向上为y正方向,右手定则确定z正方向,xy平面穿过每级的质心,也就是说坐标原点及各级质心同时位于xy平面上;下标1表示一级参数,下标2表示二级参数,下标m表示风洞模型参数,下标s表示真实飞行器参数,因此x1s表示一级真实飞行器质心x坐标,y1s表示一级真实飞行器质心y坐标,x2s表示真实飞行器质心x坐标,y2s表示二级真实飞行器质心y坐标。
本发明与现有技术相比的优点在于:
(1)本发明解决了两级甚至多级全自由飞试验无相似律指导的难题。以往风洞自由飞试验相似律主要针对单个飞行器,或者一个飞机固定,另外一个模型完全自由的飞行试验。缺少针对两个甚至多个飞行器完全自由的风洞自由飞试验相似律。本发明针对此方面进行理论推导,得到了新的相似律设计方法;;
(2)本发明解决了升力体等复杂外形无相似律的现状。升力体外形与常规外形相比具有较大气动力,当气动力较大时,相似律推导过程需要引入气动力参数,这无疑增加了相似方程的复杂性,增加了求解的难度
(3)本发明升力体外形同样增加了模型的水平位移及角位移值,再加上之前需要考虑的竖直位移,相似参数将更多,需要推导的方程也将更复杂。
(4)通过联立线位移运动方程,求解出满足线位移相似所需的条件:满足风洞试验模型多级间的质量比以及与真实飞行器几何相似,便能满足风洞自由飞试验结果与真实飞行器分离结果相一致,达到指导真实飞行器设计的目的。
(5)通过联立角位移运动方程,求解出满足角位移相似所需的条件:满足风洞试验模型的限制方程,便能满足风洞自由飞试验结果与真实飞行器分离结果相一致,达到指导真实飞行器设计的目的。
附图说明
图1为本发明的方法实施流程图。
具体实施方式
下面结合附图和具体实施例对本发明做进一步详细描述。
本发明公开了一种满足多级全自由飞风洞试验相似律的条件确定方法。根据动力学及运动学推导投放类风洞自由飞试验相似律。首次针对升力体外形针对两级甚至多级模型完全自由的自由飞试验设计试验相似律,重点解决升力体外形相似性问题。在风洞试验模型满足求限制条件时,可保证风洞试验所获得的水平线位移和竖直线位移、角位移曲线与真实飞行器的曲线相似,保证风洞试验可模拟真实飞行器的分离情况。
两级入轨飞行器分离安全具有不确定性。为保证分离安全,需针对两级入轨飞行器的分离特性进行研究。近几年多体分离全自由飞试验技术有了相当大的进步,虽然试验技术进步很大,但针对多体分离全自由飞风洞试验的相似律尚未建立,需要进行针对性研究。两级甚至多级都是完全自由的,需要的相似条件较多。风洞试验与真实飞行不仅要满足水平线位移和竖直线位移的同时相似,还需要满足角位移的相似,更有甚者需要满足上述3项相似的两两相似,这就导致了相似方程较多。升力体外形下两级各自所受气动力与其重力相比不再是小量,在推导相似条件时需重点考虑。若按照以前试验方法随意设计两级风洞试验模型参数,则会导致风洞试验所得到的两级间水平线位移和竖直线位移以及角位移与真实飞行器的三者不对应,试验误差大,对真实飞行器产生错误判断及指导。
本发明克服现有技术不足,提供一种满足多级全自由飞风洞试验相似律的条件确定方法,解决了以往风洞试验与真实飞行器之间水平线位移、竖直线位移以及角位移不对应的问题,提高了风洞试验精度,对真实飞行器研制具有指导作用。
风洞自由飞试验一般流程是将多级模型作为一个整体、以一定速度发射到流场中,当整体自由之后再进行两级之间的分离,实现两级的各自自由,模拟真实飞行器的分离过程,预测判断真实飞行器分离是否安全。
如图1所示,设定多级全自由飞相似律风洞试验中在tm时刻一级质心坐标(x1m,y1m)和二级质心坐标(x2m,y2m);
气流方向为x正方向,竖直向上为y正方向,右手定则确定z正方向。xy平面穿过每级的质心,也就是说坐标原点及各级质心同时位于xy平面上。下标1表示一级参数,下标2表示二级参数,下标m表示风洞模型参数,下标s表示真实飞行器参数。因此x1m表示一级风洞试验模型质心x坐标,y1m表示一级风洞试验模型质心y坐标,x2m表示二级风洞试验模型质心x坐标,y2m表示二级风洞试验模型质心y坐标;
一级风洞试验模型模拟一级真实飞行器,结构功能完全相同,只是尺寸缩小;二级风洞试验模型模拟二级真实飞行器,结构功能完全相同,只是尺寸缩小;
a1xm表示一级风洞试验模型水平加速度,a1ym表示一级风洞试验模型竖直加速度,a2xm表示二级风洞试验模型水平加速度,a2ym表示二级风洞试验模型竖直加速度,g为重力加速度,v0m表示一级二级共同的运动速度。
在tm时刻一级质心坐标(x1m,y1m)、二级质心坐标(x2m,y2m)。
Figure BDA0002112184040000101
Figure BDA0002112184040000102
Figure BDA0002112184040000103
Figure BDA0002112184040000104
S1m表示一级风洞试验模型的参考面积,C1ym表示一级风洞试验模型升力系数,C1xm表示一级风洞试验模型阻力系数,m1m表示一级风洞试验模型质量,S2m表示二级风洞试验模型的参考面积,C2ym表示二级风洞试验模型升力系数,C2xm表示二级风洞试验模型阻力系数,m2m表示二级风洞试验模型质量。
线加速度关系式:
a1ym=q∞m·S1m·C1ym/m1m (21)
a1xm=q∞m·S1m·C1xm/m1m (22)
a2ym=q∞m·S2m·C2ym/m2m (23)
a2xm=q∞m·S2m·C2xm/m2m (24)
Figure BDA0002112184040000111
设一级二级之间有如下关系:m1=n*m2,S1=m*S2,m和n为常数。则:
Figure BDA0002112184040000112
真实飞行条件:
气流方向为x正方向,竖直向上为y正方向,右手定则确定z正方向。xy平面穿过每级的质心,也就是说坐标原点及各级质心同时位于xy平面上。下标1表示一级参数,下标2表示二级参数,下标m表示风洞模型参数,下标s表示真实飞行器参数。因此x1s表示一级真实飞行器质心x坐标,y1s表示一级真实飞行器质心y坐标,x2s表示真实飞行器质心x坐标,y2s表示二级真实飞行器质心y坐标。
在ts时刻一级质心坐标(x1s,y1s)、二级质心坐标(x2s,y2s)。
Figure BDA0002112184040000113
Figure BDA0002112184040000121
Figure BDA0002112184040000122
Figure BDA0002112184040000123
S1s表示一级真实飞行器的参考面积,C1ys表示一级真实飞行器升力系数,C1xs表示一级真实飞行器阻力系数,s1s表示一级真实飞行器质量,S2s表示二级真实飞行器的参考面积,C2ys表示二级真实飞行器升力系数,C2xs表示二级真实飞行器阻力系数,s2s表示二级真实飞行器质量。
线加速度关系式:
a1ys=q∞m·S1s·C1ys/m1s (31)
a1xs=q∞s·S1s·C1xs/m1s (32)
a2ys=q∞s·S2s·C2ys/m2s (33)
a2xs=q∞s·S2s·C2xs/m2s (34)
一级二级之间同样有如下关系:m1=n*m2,S1=m*S2,m和n为常数。则:
Figure BDA0002112184040000124
另外,已经证明:
Figure BDA0002112184040000125
此外,在满足几何相似、飞行条件基本相同、飞行姿态角一致、飞行姿态角在一定范围内,有:
C2ys=C2ym (37)
C1ys=C1ym (38)
C2xs=C2xm (39)
C1xs=C1xm (40)
因而,有恒等式:
Cs=Cm (41)
由此可见,两级模型完全自由状态下,只要保证两级的几何相似和质量比相同,风洞自由飞试验与真实飞行器的线位移轨迹是恒相似的,与模型质量、来流条件无关。
另外,推导运动时刻相似:
Figure BDA0002112184040000131
Figure BDA0002112184040000132
可求得:
Figure BDA0002112184040000133
角位移相似所需条件:
角位移相似需满足以下三个方程:
Figure BDA0002112184040000134
而:
Figure BDA0002112184040000141
因此角位移相似需满足:
Figure BDA0002112184040000142
然而已求得关系式(39),因此满足关系式(42),即满足关系式(40),满足角位移相似。并且有:
Figure BDA0002112184040000143
Figure BDA0002112184040000144
根据时间相似可求得:
Figure BDA0002112184040000145
之前提到风洞全自由飞分离试验与真实飞行器的线位移轨迹是恒相似的,但为了达到横向线位移、竖直线位移、角位移3者的时间相似,则会限制模型的质量。
本发明要解决的技术问题是:首次针对升力体外形针对两级甚至多级模型完全自由的自由飞试验设计试验相似律,重点解决升力体外形相似性问题,在保证风洞试验与真实飞行水平线位移和竖直线位移相似的同时,满足角位移的相似,甚至需要满足上述3项的两两相似。
本发明解决了两级甚至多级全自由飞试验无相似律指导的难题。以往风洞自由飞试验相似律主要针对单个飞行器,或者一个飞机固定,另外一个模型完全自由的飞行试验。缺少针对两个甚至多个飞行器完全自由的风洞自由飞试验相似律。本发明针对此方面进行理论推导,得到了新的相似律设计方法;本发明解决了升力体等复杂外形无相似律的现状。升力体外形与常规外形相比具有较大气动力,当气动力较大时,相似律推导过程需要引入气动力参数,这无疑增加了相似方程的复杂性,增加了求解的难度
本发明升力体外形同样增加了模型的水平位移及角位移值,再加上之前需要考虑的竖直位移,相似参数将更多,需要推导的方程也将更复杂。
通过联立线位移运动方程,求解出满足线位移相似所需的条件:满足风洞试验模型多级间的质量比以及与真实飞行器几何相似,便能满足风洞自由飞试验结果与真实飞行器分离结果相一致,达到指导真实飞行器设计的目的。通过联立角位移运动方程,求解出满足角位移相似所需的条件:满足风洞试验模型的限制方程,便能满足风洞自由飞试验结果与真实飞行器分离结果相一致,达到指导真实飞行器设计的目的。

Claims (10)

1.一种满足多级全自由飞风洞试验相似律的条件确定方法,其特征在于步骤如下:
(1)设定多级全自由飞相似律风洞试验中在tm时刻一级质心坐标为(x1m,y1m)和二级质心坐标为(x2m,y2m);tm表示多级全自由飞相似律风洞试验中一级与二级分离过程中的某个时间点;
Figure FDA0002716046580000011
Figure FDA0002716046580000012
Figure FDA0002716046580000013
Figure FDA0002716046580000014
a1xm表示一级风洞试验模型水平加速度,a1ym表示一级风洞试验模型竖直加速度,a2xm表示二级风洞试验模型水平加速度,a2ym表示二级风洞试验模型竖直加速度,g为重力加速度,v0m表示一级二级共同的运动速度;
其中,a1ym=q∞m·S1m·C1ym/m1m (5)
a1xm=q∞m·S1m·C1xm/m1m (6)
a2ym=q∞m·S2m·C2ym/m2m (7)
a2xm=q∞m·S2m·C2xm/m2m (8)
S1m表示一级风洞试验模型的参考面积,C1ym表示一级风洞试验模型升力系数,C1xm表示一级风洞试验模型阻力系数,m1m表示一级风洞试验模型质量,S2m表示二级风洞试验模型的参考面积,C2ym表示二级风洞试验模型升力系数,C2xm表示二级风洞试验模型阻力系数,m2m表示二级风洞试验模型质量;
(2)根据步骤(1)多级全自由飞相似律风洞试验中在tm时刻一级质心坐标(x1m,y1m)和二级质心坐标(x2m,y2m),确定风洞试验中tm时刻一级与二级的相对位置参数Cm
Figure FDA0002716046580000021
(3)设一级与二级之间有如下关系:m1m=N*m2m,S1m=M*S2m,当真实飞行器一级和二级质量、外形确定之后,M和N便为常数,则将公式(9)化简得到公式(10)如下:
Figure FDA0002716046580000022
(4)在步骤(1)的同时,设定真实飞行器中,在ts时刻一级质心坐标(x1s,y1s)、二级质心坐标(x2s,y2s);
ts表示真实飞行器一级与二级分离过程中的某个时间点;
Figure FDA0002716046580000023
Figure FDA0002716046580000024
Figure FDA0002716046580000025
Figure FDA0002716046580000026
a1xs表示一级真实飞行器水平加速度,a1ys表示一级真实飞行器竖直加速度,a2xs表示二级真实飞行器水平加速度,a2ys表示二级真实飞行器竖直加速度,g为重力加速度,v0s表示一级二级共同的运动速度;
其中:
a1ys=q∞m·S1s·C1ys/m1s (15)
a1xs=q∞s·S1s·C1xs/m1s (16)
a2ys=q∞s·S2s·C2ys/m2s (17)
a2xs=q∞s·S2s·C2xs/m2s (18)
S1s表示一级真实飞行器的参考面积,C1ys表示一级真实飞行器升力系数,C1xs表示一级真实飞行器阻力系数,s1s表示一级真实飞行器质量,S2s表示二级真实飞行器的参考面积,C2ys表示二级真实飞行器升力系数,C2xs表示二级真实飞行器阻力系数,s2s表示二级真实飞行器质量;
(5)根据步骤(4)设定真实飞行器中,在ts时刻一级质心坐标(x1s,y1s)、二级质心坐标(x2s,y2s),确定真实飞行器ts时刻一级与二级的相对位置参数Cs
设一级二级之间同样有如下关系:m1=N*m2,S1=M*S2,M和N为定值,则:
Figure FDA0002716046580000031
(6)根据公式(10),即公式(20):
Figure FDA0002716046580000032
(7)根据公式(19)和公式(10),则有在满足多级全自由飞相似律风洞试验与真实飞行器的几何相似、飞行条件相同、飞行姿态角一致的情况下:
C2ys=C2ym (21)
C1ys=C1ym (22)
C2xs=C2xm (23)
C1xs=C1xm (24)
因而,有恒等式:
Cs=Cm (25)
(m1m/m2m=m1s/m2s)
Cs=Cm即表示两级模型完全自由状态下,只要保证一级与二级的几何相似和质量比相同,风洞自由飞试验与真实飞行器的线位移轨迹是恒相似的,与风洞自由飞试验中的飞行器模型的质量、来流条件无关;即满足风洞自由飞试验与真实飞行器的线位移轨迹相似的条件;
(8)设定风洞自由飞试验与真实飞行器在对应时刻的水平位移成比例,则有:
Figure FDA0002716046580000041
Figure FDA0002716046580000042
根据(26-1)和(26-2),确定公式(27)
Figure FDA0002716046580000043
化简公式(27),求得:
Figure FDA0002716046580000044
式中,kt表示时间缩比,q表示来流动压,kq∞表示来流动压比,k1表示尺寸缩比,km表示质量缩比;
(9)设定风洞自由飞试验与真实飞行器的角位移轨迹是恒相似的,则需满足以下三个方程:
Figure FDA0002716046580000051
kθ表示角度缩比,
Figure FDA0002716046580000052
表示角速度缩比,
Figure FDA0002716046580000053
表示角加速度缩比;
根据物体角加速度定义,则有:
Figure FDA0002716046580000054
因此,根据公式(29)和公式(30),则角位移相等需满足:
Figure FDA0002716046580000055
然而已求得关系式(28),因此满足关系式(31),即满足关系式(29),满足角位移相似;将公式(28)代入公式(29),则有:
Figure FDA0002716046580000056
Figure FDA0002716046580000057
根据公式kt=kl/kv (34)
再结合公式(28),则有:
Figure FDA0002716046580000058
公式(35)即满足风洞自由飞试验与真实飞行器的角位移轨迹相似的条件。
2.根据权利要求1所述的一种满足多级全自由飞风洞试验相似律的条件确定方法,其特征在于:当进行两级全自由飞风洞试验时,脱离飞行主体的为一级,继续保持飞行的为二级;为了模拟真实飞行,需进行风洞试验,风洞试验开始时,一级二级作为刚性整体无支撑的在风洞中飞行,之间无相对自由度,在需要进行分离的时刻,突然释放一级二级之间的6自由度限制,在气动力的作用下,两级姿态产生变化,模拟真实飞行器的分离过程。
3.根据权利要求1所述的一种满足多级全自由飞风洞试验相似律的条件确定方法,其特征在于:多级全自由飞风洞试验与真实飞行器的几何相似,是指:风洞试验模型与真实飞行器之间对应边长度成比例,对应角度相等。
4.根据权利要求1所述的一种满足多级全自由飞风洞试验相似律的条件确定方法,其特征在于:多级全自由飞风洞试验与真实飞行器的飞行条件相同,是指:风洞试验Ma数与真实飞行器飞行Ma数相等。
5.根据权利要求1所述的一种满足多级全自由飞风洞试验相似律的条件确定方法,其特征在于:多级全自由飞风洞试验与真实飞行器的飞行姿态角一致,是指:风洞试验模型运动的角位移与真实飞行器对应点的角位移相等。
6.根据权利要求1所述的一种满足多级全自由飞风洞试验相似律的条件确定方法,其特征在于:一级与二级的几何相似,是指:风洞试验模型与真实飞行器的一级与二级之间对应边长度成比例,对应角度相等。
7.根据权利要求1所述的一种满足多级全自由飞风洞试验相似律的条件确定方法,其特征在于:一级与二级的质量比相同,是指:真实飞行器一级二级之间的质量比要与风洞试验一级二级之间的质量比相等。
8.根据权利要求1所述的一种满足多级全自由飞风洞试验相似律的条件确定方法,其特征在于:一级质心坐标(x1m,y1m)和二级质心坐标(x2m,y2m),所在的坐标系为:在两级进行分离前,一级质心坐标和二级质心坐标的同在xy平面内,且二者连线与y轴平行。
9.根据权利要求1所述的一种满足多级全自由飞风洞试验相似律的条件确定方法,其特征在于:建立三维坐标系,气流方向为x正方向,竖直向上为y正方向,右手定则确定z正方向,xy平面穿过每级的质心,也就是说坐标原点及各级质心同时位于xy平面上;下标1表示一级参数,下标2表示二级参数,下标m表示风洞模型参数,下标s表示真实飞行器参数,因此x1m表示一级风洞试验模型质心x坐标,y1m表示一级风洞试验模型质心y坐标,x2m表示二级风洞试验模型质心x坐标,y2m表示二级风洞试验模型质心y坐标。
10.根据权利要求1所述的一种满足多级全自由飞风洞试验相似律的条件确定方法,其特征在于:气流方向为x正方向,竖直向上为y正方向,右手定则确定z正方向,xy平面穿过每级的质心,也就是说坐标原点及各级质心同时位于xy平面上;下标1表示一级参数,下标2表示二级参数,下标m表示风洞模型参数,下标s表示真实飞行器参数,因此x1s表示一级真实飞行器质心x坐标,y1s表示一级真实飞行器质心y坐标,x2s表示真实飞行器质心x坐标,y2s表示二级真实飞行器质心y坐标。
CN201910576293.0A 2019-06-28 2019-06-28 一种满足多级全自由飞风洞试验相似律的条件确定方法 Active CN110398339B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910576293.0A CN110398339B (zh) 2019-06-28 2019-06-28 一种满足多级全自由飞风洞试验相似律的条件确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910576293.0A CN110398339B (zh) 2019-06-28 2019-06-28 一种满足多级全自由飞风洞试验相似律的条件确定方法

Publications (2)

Publication Number Publication Date
CN110398339A CN110398339A (zh) 2019-11-01
CN110398339B true CN110398339B (zh) 2021-02-09

Family

ID=68323626

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910576293.0A Active CN110398339B (zh) 2019-06-28 2019-06-28 一种满足多级全自由飞风洞试验相似律的条件确定方法

Country Status (1)

Country Link
CN (1) CN110398339B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110889171B (zh) * 2019-12-02 2023-08-18 中国空气动力研究与发展中心 一种立式风洞飞机尾旋试验模型设计方法
CN112504613B (zh) * 2021-02-03 2021-05-28 中国空气动力研究与发展中心高速空气动力研究所 一种并联飞行器级间分离试验方法、装置及可读存储介质
CN113588205B (zh) * 2021-08-27 2023-05-16 重庆电子工程职业学院 基于相似理论的实验室风洞尾支杆隔振平台设计方法
CN114486165B (zh) * 2021-12-27 2023-11-03 中国航天空气动力技术研究院 一种斜向下分离风洞自由投放试验相似准则设计方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0972822A (ja) * 1995-09-07 1997-03-18 Mitsubishi Heavy Ind Ltd 風洞試験装置
JP4310440B2 (ja) * 2006-03-22 2009-08-12 防衛省技術研究本部長 動的風洞試験方法
CN102305699A (zh) * 2011-05-19 2012-01-04 北京航空航天大学 自由飞模型风洞实验系统
CN103471803B (zh) * 2013-09-24 2016-03-02 中国航天空气动力技术研究院 一种模型自由飞试验的气动参数确定方法
CN104102127B (zh) * 2014-07-17 2016-06-01 厦门大学 一种机载气动参数辨识系统
WO2018122575A1 (en) * 2016-12-29 2018-07-05 Aerodium Technologies, Sia Recirculating vertical wind tunnel with passive cooling
CN107991053B (zh) * 2017-11-06 2020-04-10 中国航天空气动力技术研究院 高速风洞投放模型试验飞行器分离轨迹预测方法及系统
CN108132134A (zh) * 2017-11-15 2018-06-08 南京航空航天大学 基于风洞自由飞试验的气动导数辨识方法和系统
CN108318219B (zh) * 2018-01-19 2020-03-24 中国航天空气动力技术研究院 自由飞试验相似律重力补偿设计方法
CN109297666B (zh) * 2018-10-15 2020-03-31 中国空气动力研究与发展中心高速空气动力研究所 一种基于两套运动机构的级间分离风洞试验装置和试验方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Natural laminar flow wing for supersonic conditions: Wind tunnelexperiments, flight test and stability computations;Olivier Vermeersch 等;《Progress in Aerospace Sciences》;20150815;第64-91页 *

Also Published As

Publication number Publication date
CN110398339A (zh) 2019-11-01

Similar Documents

Publication Publication Date Title
CN110398339B (zh) 一种满足多级全自由飞风洞试验相似律的条件确定方法
CN110398340B (zh) 基于投放分离类风洞自由飞试验相似律简化的风洞试验参数确定方法
CN114200826B (zh) 一种超音速大机动靶标持续大过载机动高度稳定控制方法
CN106153291B (zh) 补偿高速风洞弹射投放模型垂直加速度不足的方法
CN109190248B (zh) 一种用于滑翔飞行器的滑翔射程解析方法及解析系统
CN114200827B (zh) 一种超音速大机动靶标的多约束双通道控制方法
CN111964862A (zh) 一种机弹分离动力学相似风洞试验方法
CN110316358A (zh) 基于动态逆的战斗机大迎角控制方法
CN112364432A (zh) 一种载机挂飞投放分离过程控制方法
CN106570242B (zh) 低空大动压整体式整流罩高速分离流固耦合仿真方法
CN109141802B (zh) 用于捕获轨迹试验中的外挂物控制律的模拟方法
CN110398342B (zh) 飞机挂载物分离风洞自由飞试验相似律设计方法
CN111723485A (zh) 一种降落伞伞衣透气速度获得方法
CN110989397B (zh) 一种航空器失事搜寻仿真方法及系统
CN211685678U (zh) 一种多旋翼无人机实时尾迹的仿真分析系统
CN113655808B (zh) 用于集群的有限时间自主编队控制方法及控制系统
CN113868779B (zh) 基于多学科耦合仿真的螺旋桨飞机桨叶飞脱安全评估方法
CN113094817B (zh) 一种直升机武器发射响应的计算方法
Belokon et al. Control of hybrid unmanned aerial vehicle motion in transitional modes
Lofthouse et al. Static and dynamic simulations of a generic UCAV geometry using the kestrel flow solver
CN118013893B (zh) 用于变体飞行器的变体抛离安全分离条件判定方法及系统
Tai et al. Test Data Processing of Fly-by-Wire Civil Aircraft in Low-Speed Wind Tunnel Virtual Flight
CN114486165B (zh) 一种斜向下分离风洞自由投放试验相似准则设计方法
Afshar et al. Design and fabrication of a delta wing micro aerial vehicle
Heinrich et al. DLR pilot applications in ComFlite

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