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

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

Info

Publication number
CN110398339A
CN110398339A CN201910576293.0A CN201910576293A CN110398339A CN 110398339 A CN110398339 A CN 110398339A CN 201910576293 A CN201910576293 A CN 201910576293A CN 110398339 A CN110398339 A CN 110398339A
Authority
CN
China
Prior art keywords
level
wind tunnel
indicate
actual aircraft
test
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.)
Granted
Application number
CN201910576293.0A
Other languages
English (en)
Other versions
CN110398339B (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

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表示多级全自由飞相似律风洞试验中一级与二级分离过程中的某个时间点;
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
(3)设一级与二级之间有如下关系:m1m=N*m2m,S1m=M*S2m,当真实飞行器一级和二级质量、外形确定之后,M和N便为常数,则将公式(9)化简得到公式(10)如下:
(4)在步骤(1)的同时,设定真实飞行器中,在ts时刻一级质心坐标(x1s,y1s)、二级质心坐标(x2s,y2s);
ts表示真实飞行器一级与二级分离过程中的某个时间点;
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为定值,则:
(6)根据公式(10),即公式(20):
(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)设定风洞自由飞试验与真实飞行器在对应时刻的水平位移成比例,则有:
根据(26-1)和(26-2),确定公式(27)
化简公式(27),求得:
式中,kt表示时间缩比,q表示来流动压,kq∞表示来流动压比,kl表示尺寸缩比,km表示质量缩比;
(9)设定风洞自由飞试验与真实飞行器的角位移轨迹是恒相似的,则需满足以下三个方程:
kθ表示角度缩比(角度是指飞行器的姿态角),表示角速度缩比,表示角加速度缩比;
根据物体角加速度定义,则有:
因此,根据公式(29)和公式(30),则角位移相等需满足:
然而已求得关系式(28),因此满足关系式(31),即满足关系式(29),满足角位移相似;将公式(28)代入公式(29),则有:
根据公式kt=kl/kv (34)
再结合公式(28),则有:
公式(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)。
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)
设一级二级之间有如下关系:m1=n*m2,S1=m*S2,m和n为常数。则:
真实飞行条件:
气流方向为x正方向,竖直向上为y正方向,右手定则确定z正方向。xy平面穿过每级的质心,也就是说坐标原点及各级质心同时位于xy平面上。下标1表示一级参数,下标2表示二级参数,下标m表示风洞模型参数,下标s表示真实飞行器参数。因此x1s表示一级真实飞行器质心x坐标,y1s表示一级真实飞行器质心y坐标,x2s表示真实飞行器质心x坐标,y2s表示二级真实飞行器质心y坐标。
在ts时刻一级质心坐标(x1s,y1s)、二级质心坐标(x2s,y2s)。
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为常数。则:
另外,已经证明:
此外,在满足几何相似、飞行条件基本相同、飞行姿态角一致、飞行姿态角在一定范围内,有:
C2ys=C2ym (37)
C1ys=C1ym (38)
C2xs=C2xm (39)
C1xs=C1xm (40)
因而,有恒等式:
Cs=Cm (41)
由此可见,两级模型完全自由状态下,只要保证两级的几何相似和质量比相同,风洞自由飞试验与真实飞行器的线位移轨迹是恒相似的,与模型质量、来流条件无关。
另外,推导运动时刻相似:
可求得:
角位移相似所需条件:
角位移相似需满足以下三个方程:
而:
因此角位移相似需满足:
然而已求得关系式(39),因此满足关系式(42),即满足关系式(40),满足角位移相似。并且有:
根据时间相似可求得:
之前提到风洞全自由飞分离试验与真实飞行器的线位移轨迹是恒相似的,但为了达到横向线位移、竖直线位移、角位移3者的时间相似,则会限制模型的质量。
本发明要解决的技术问题是:首次针对升力体外形针对两级甚至多级模型完全自由的自由飞试验设计试验相似律,重点解决升力体外形相似性问题,在保证风洞试验与真实飞行水平线位移和竖直线位移相似的同时,满足角位移的相似,甚至需要满足上述3项的两两相似。
本发明解决了两级甚至多级全自由飞试验无相似律指导的难题。以往风洞自由飞试验相似律主要针对单个飞行器,或者一个飞机固定,另外一个模型完全自由的飞行试验。缺少针对两个甚至多个飞行器完全自由的风洞自由飞试验相似律。本发明针对此方面进行理论推导,得到了新的相似律设计方法;本发明解决了升力体等复杂外形无相似律的现状。升力体外形与常规外形相比具有较大气动力,当气动力较大时,相似律推导过程需要引入气动力参数,这无疑增加了相似方程的复杂性,增加了求解的难度
本发明升力体外形同样增加了模型的水平位移及角位移值,再加上之前需要考虑的竖直位移,相似参数将更多,需要推导的方程也将更复杂。
通过联立线位移运动方程,求解出满足线位移相似所需的条件:满足风洞试验模型多级间的质量比以及与真实飞行器几何相似,便能满足风洞自由飞试验结果与真实飞行器分离结果相一致,达到指导真实飞行器设计的目的。通过联立角位移运动方程,求解出满足角位移相似所需的条件:满足风洞试验模型的限制方程,便能满足风洞自由飞试验结果与真实飞行器分离结果相一致,达到指导真实飞行器设计的目的。

Claims (10)

1.一种满足多级全自由飞风洞试验相似律的条件确定方法,其特征在于步骤如下:
(1)设定多级全自由飞相似律风洞试验中在tm时刻一级质心坐标为(x1m,y1m)和二级质心坐标为(x2m,y2m);tm表示多级全自由飞相似律风洞试验中一级与二级分离过程中的某个时间点;
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
(3)设一级与二级之间有如下关系:m1m=N*m2m,S1m=M*S2m,当真实飞行器一级和二级质量、外形确定之后,M和N便为常数,则将公式(9)化简得到公式(10)如下:
(4)在步骤(1)的同时,设定真实飞行器中,在ts时刻一级质心坐标(x1s,y1s)、二级质心坐标(x2s,y2s);
ts表示真实飞行器一级与二级分离过程中的某个时间点;
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为定值,则:
(6)根据公式(10),即公式(20):
(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)设定风洞自由飞试验与真实飞行器在对应时刻的水平位移成比例,则有:
根据(26-1)和(26-2),确定公式(27)
化简公式(27),求得:
式中,kt表示时间缩比,q表示来流动压,kq∞表示来流动压比,kl表示尺寸缩比,km表示质量缩比;
(9)设定风洞自由飞试验与真实飞行器的角位移轨迹是恒相似的,则需满足以下三个方程:
kθ表示角度缩比,表示角速度缩比,表示角加速度缩比;
根据物体角加速度定义,则有:
因此,根据公式(29)和公式(30),则角位移相等需满足:
然而已求得关系式(28),因此满足关系式(31),即满足关系式(29),满足角位移相似;将公式(28)代入公式(29),则有:
根据公式kt=kl/kv (34)
再结合公式(28),则有:
公式(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 true CN110398339A (zh) 2019-11-01
CN110398339B 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)

Cited By (4)

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

Citations (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 風洞試験装置
JP2006208395A (ja) * 2006-03-22 2006-08-10 Tech Res & Dev Inst Of Japan Def Agency 動的風洞試験方法
CN102305699A (zh) * 2011-05-19 2012-01-04 北京航空航天大学 自由飞模型风洞实验系统
CN103471803A (zh) * 2013-09-24 2013-12-25 中国航天空气动力技术研究院 一种模型自由飞试验的气动参数确定方法
CN104102127A (zh) * 2014-07-17 2014-10-15 厦门大学 一种机载气动参数辨识系统
CN107991053A (zh) * 2017-11-06 2018-05-04 中国航天空气动力技术研究院 高速风洞投放模型试验飞行器分离轨迹预测方法及系统
CN108132134A (zh) * 2017-11-15 2018-06-08 南京航空航天大学 基于风洞自由飞试验的气动导数辨识方法和系统
WO2018122575A1 (en) * 2016-12-29 2018-07-05 Aerodium Technologies, Sia Recirculating vertical wind tunnel with passive cooling
CN108318219A (zh) * 2018-01-19 2018-07-24 中国航天空气动力技术研究院 自由飞试验相似律重力补偿设计方法
CN109297666A (zh) * 2018-10-15 2019-02-01 中国空气动力研究与发展中心高速空气动力研究所 一种基于两套运动机构的级间分离风洞试验装置和试验方法

Patent Citations (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 風洞試験装置
JP2006208395A (ja) * 2006-03-22 2006-08-10 Tech Res & Dev Inst Of Japan Def Agency 動的風洞試験方法
CN102305699A (zh) * 2011-05-19 2012-01-04 北京航空航天大学 自由飞模型风洞实验系统
CN103471803A (zh) * 2013-09-24 2013-12-25 中国航天空气动力技术研究院 一种模型自由飞试验的气动参数确定方法
CN104102127A (zh) * 2014-07-17 2014-10-15 厦门大学 一种机载气动参数辨识系统
WO2018122575A1 (en) * 2016-12-29 2018-07-05 Aerodium Technologies, Sia Recirculating vertical wind tunnel with passive cooling
CN107991053A (zh) * 2017-11-06 2018-05-04 中国航天空气动力技术研究院 高速风洞投放模型试验飞行器分离轨迹预测方法及系统
CN108132134A (zh) * 2017-11-15 2018-06-08 南京航空航天大学 基于风洞自由飞试验的气动导数辨识方法和系统
CN108318219A (zh) * 2018-01-19 2018-07-24 中国航天空气动力技术研究院 自由飞试验相似律重力补偿设计方法
CN109297666A (zh) * 2018-10-15 2019-02-01 中国空气动力研究与发展中心高速空气动力研究所 一种基于两套运动机构的级间分离风洞试验装置和试验方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
蒋增辉 等: "多体分离风洞自由飞试验", 《空气动力学学报》 *

Cited By (7)

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

Also Published As

Publication number Publication date
CN110398339B (zh) 2021-02-09

Similar Documents

Publication Publication Date Title
CN110398339A (zh) 一种满足多级全自由飞风洞试验相似律的条件确定方法
Babinsky How do wings work?
CN106446442B (zh) 一种火星伞舱组合体的降落伞展开过程稳定性评估方法
EP1729234A1 (en) Method for designing streamline traced, mixed compression inlets for aircraft engines
CN106153291B (zh) 补偿高速风洞弹射投放模型垂直加速度不足的方法
CN106647792A (zh) 用于无人机悬挂负载系统的抗扰控制方法
CN110398340B (zh) 基于投放分离类风洞自由飞试验相似律简化的风洞试验参数确定方法
CN109190248B (zh) 一种用于滑翔飞行器的滑翔射程解析方法及解析系统
Serra et al. Nonlinear image-based visual servo controller for the flare maneuver of fixed-wing aircraft using optical flow
CN114200827B (zh) 一种超音速大机动靶标的多约束双通道控制方法
CN108516104A (zh) 多级运载航天器分离可靠性仿真平台
CN106570242B (zh) 低空大动压整体式整流罩高速分离流固耦合仿真方法
CN110929339B (zh) 高精度亚音速固定翼飞行器非线性动力学建模方法
CN111723485A (zh) 一种降落伞伞衣透气速度获得方法
Li et al. Simulation method for wind tunnel based virtual flight testing
Askari et al. Two-and three-dimensional numerical simulations of supersonic ramped inlet
CN107247464B (zh) 一种四旋翼无人飞行器的状态受限控制方法和系统
Triputra et al. Nonlinear dynamic modeling of a fixed-wing unmanned aerial vehicle: A case study of Wulung
CN114878123B (zh) 飞机疲劳强度测试中基于加速模型的组合载荷谱编制方法
CN106372307B (zh) 一种基于气动模型的民机气流角估计方法
Shi et al. Design of Quadrotor Intelligent Rescue UAV Based on UAV Vision
Argyle et al. Quaternion based attitude error for a tailsitter in hover flight
Belokon et al. Control of hybrid unmanned aerial vehicle motion in transitional modes
Ida et al. A quadcopters flight simulation considering the influence of wind
Chu et al. Aerodynamic configuration and control optimization for a novel horizontal-rope shipborne recovery fixed-wing UAV system

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