CN112977869B - 一种直升机大气数据系统旋翼下洗流影响修正方法 - Google Patents
一种直升机大气数据系统旋翼下洗流影响修正方法 Download PDFInfo
- Publication number
- CN112977869B CN112977869B CN202110210586.4A CN202110210586A CN112977869B CN 112977869 B CN112977869 B CN 112977869B CN 202110210586 A CN202110210586 A CN 202110210586A CN 112977869 B CN112977869 B CN 112977869B
- Authority
- CN
- China
- Prior art keywords
- pressure
- dynamic pressure
- angle
- coefficient
- attack
- 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
Links
- 238000012937 correction Methods 0.000 title claims abstract description 78
- 238000000034 method Methods 0.000 title claims abstract description 26
- 239000000523 sample Substances 0.000 claims abstract description 36
- 230000003068 static effect Effects 0.000 claims abstract description 34
- 238000004364 calculation method Methods 0.000 claims abstract description 20
- 238000005406 washing Methods 0.000 claims abstract description 11
- 101100494411 Oryza sativa subsp. japonica CYP74A3 gene Proteins 0.000 claims description 18
- 101100216180 Solanum lycopersicum AOS3 gene Proteins 0.000 claims description 18
- 238000012360 testing method Methods 0.000 claims description 10
- 238000004088 simulation Methods 0.000 claims description 9
- 101100494409 Oryza sativa subsp. japonica CYP74A2 gene Proteins 0.000 claims description 6
- 101100216179 Solanum lycopersicum AOS2 gene Proteins 0.000 claims description 6
- 238000009530 blood pressure measurement Methods 0.000 claims description 6
- 238000012986 modification Methods 0.000 claims description 6
- 230000004048 modification Effects 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims 1
- 238000003786 synthesis reaction Methods 0.000 claims 1
- 238000005259 measurement Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 4
- 238000009795 derivation Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- ZGXJTSGNIOSYLO-UHFFFAOYSA-N 88755TAZ87 Chemical compound NCC(=O)CCC(O)=O ZGXJTSGNIOSYLO-UHFFFAOYSA-N 0.000 description 1
- 108700032845 Ala(2)- enkephalinamide-Met Proteins 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P21/00—Testing or calibrating of apparatus or devices covered by the preceding groups
- G01P21/02—Testing or calibrating of apparatus or devices covered by the preceding groups of speedometers
- G01P21/025—Testing or calibrating of apparatus or devices covered by the preceding groups of speedometers for measuring speed of fluids; for measuring speed of bodies relative to fluids
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64F—GROUND OR AIRCRAFT-CARRIER-DECK INSTALLATIONS SPECIALLY ADAPTED FOR USE IN CONNECTION WITH AIRCRAFT; DESIGNING, MANUFACTURING, ASSEMBLING, CLEANING, MAINTAINING OR REPAIRING AIRCRAFT, NOT OTHERWISE PROVIDED FOR; HANDLING, TRANSPORTING, TESTING OR INSPECTING AIRCRAFT COMPONENTS, NOT OTHERWISE PROVIDED FOR
- B64F5/00—Designing, manufacturing, assembling, cleaning, maintaining or repairing aircraft, not otherwise provided for; Handling, transporting, testing or inspecting aircraft components, not otherwise provided for
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64C—AEROPLANES; HELICOPTERS
- B64C27/00—Rotorcraft; Rotors peculiar thereto
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64D—EQUIPMENT FOR FITTING IN OR TO AIRCRAFT; FLIGHT SUITS; PARACHUTES; ARRANGEMENT OR MOUNTING OF POWER PLANTS OR PROPULSION TRANSMISSIONS IN AIRCRAFT
- B64D43/00—Arrangements or adaptations of instruments
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P13/00—Indicating or recording presence, absence, or direction, of movement
- G01P13/02—Indicating direction only, e.g. by weather vane
- G01P13/025—Indicating direction only, e.g. by weather vane indicating air data, i.e. flight variables of an aircraft, e.g. angle of attack, side slip, shear, yaw
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Transportation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Measuring Fluid Pressure (AREA)
Abstract
一种直升机大气数据系统旋翼下洗流影响修正方法，包括以下步骤，步骤S1：建立动压修正系数表Kp；步骤S2：建立合成攻角系数表Kai；步骤S3：建立合成侧滑角系数表Kbi；步骤S4：采用随动探头测得总压和静压；步骤S5：采用随动探头直杆上的测压孔测得上压力、下压力、左压力和右压力；步骤S6：通过总压、静压、上压力和下压力的差、左压力和右压力的差解算出攻角和侧滑角；步骤S7：通过攻角和侧滑角得到动压修正系数Kp，实现真实动压的解算，完成旋翼下洗流修正。
Description
技术领域
本发明涉及飞机大气数据测量领域，特别涉及一种直升机大气数据系统旋翼下洗流影响修正方法。
背景技术
民用直升机大气数据系统是利用测量气压压力(包括总压、静压)的原理，实时测量直升机的飞行高度、速度、总温/静温、攻角、侧滑角等参数的系统或设备，是影响直升机飞行安全和导航显示不可或缺的重要系统，其失效将可能导致飞机失事，属直升机的安全关键设备，研制保证等级为DALA级(最高)，因此，民用直升机大气数据系统大气参数的准确测量对直升机飞行安全至关重要。
直升机由旋翼提供动力和升力，直升机在飞行时，旋翼下方会产生下洗气流，直升机由于其特殊气动特性，导致直升机低速飞行时布置在直升机机身上的大气数据探头必将受到直升机旋翼诱导气流的影响，导致大气数据测量不准确，影响直升机的飞行安全。
发明内容
本发明的目的在于：提供了一种直升机大气数据系统旋翼下洗流影响修正方法，针对直升机特有气动特性，在直升机低速飞行时大气数据测量易受旋翼下洗流影响，导致大气数据测量不准确的问题，通过构建下洗流感受测量构型，并建立下洗流修正算法，消除下洗流影响，保证大气数据测量的准确性，解决了上述问题。
本发明采用的技术方案如下：
一种直升机大气数据系统旋翼下洗流影响修正方法，包括以下步骤，
步骤S1：建立动压修正系数表Kp；
步骤S2：建立合成攻角系数表Kai；
步骤S3：建立合成侧滑角系数表Kbi；
步骤S4：采用随动探头测得总压和静压；
步骤S5：采用随动探头直杆上的测压孔测得上压力、下压力、左压力和右压力；
步骤S6：通过总压、静压、上压力和下压力的差、左压力和右压力的差解算出攻角和侧滑角，为获取更高的测量精度，步骤S3应重重多次；
步骤S7：通过攻角、侧滑角和马赫数，得到动压修正系数Kp，实现真实动压的解算，完成旋翼下洗流修正，为获取更高的测量精度，步骤S4应重重多次。
为了更好地实现本方案，选用旋翼开启或关闭条件，通过气动流场仿真或气动流场试验，建立不同马赫数、不同攻角、不同侧滑角下的动压修正系数Kp：
其中，Qci为指示动压：开启旋翼时的动压测量值；Qc为真实动压：关闭旋翼时的动压测量值。
为了更好地实现本方案，应在打开旋翼的情况下，通过气动流场仿真或气动流场试验，建立不同马赫数、不同攻角、不同侧滑角下的合成攻角系数Kai：
所述P(d-u)i为随动探头直杆上的下压力与上压力的差值，Qci为指示动压。
为了更好地实现本方案，应在打开旋翼的情况下，通过气动流场仿真或气动流场试验，建立不同马赫数、不同攻角下的合成侧滑角系数Kbi：
所述P(l-r)i为随动探头直杆上的左压力与右压力的差值，Qci为指示动压。
为了更好地实现本方案，进一步地，所述步骤S6重复迭代多次，以获得要求的攻角、侧滑角计算精度。
为了更好地实现本方案，进一步地，所述步骤S7重复迭代多次，以获得要求的大气参数计算精度。
为了更好地实现本方案，进一步地，所述步骤S6中的攻角、侧滑角的解算包括以下步骤：
步骤S601：根据静压Psi、动压Qci计算指示马赫数Mi：
其中k为绝热指数；
步骤S602：根据动压Qci、P(d-u)i计算当前状态下的合成攻角系数Kai：
其中P(d-u)i为随动探头直杆上的下压力与上压力的差值；
步骤S603：根据动压Qci、P(l-r)i计算当前状态下的合成侧滑角系数Kbi：
其中P(l-r)i为随动探头直杆上的左压力与右压力的差值；
步骤S604：根据指示马赫数Mi，当前合成攻角系数Kai，以及按步骤S2获得的合成攻角系数表，按侧滑角为0度时的合成攻角系数Kai，得到侧滑角为0度时的攻角值AOA0；
步骤S605：根据指示马赫数Mi，当前合成侧滑角系数Kbi，以及按步骤S3获得的合成侧滑角系数表，按攻角为AOA0时的合成侧滑角系数Kbi，得到攻角为AOA0时的侧滑角值AOS0；
步骤S606：根据指示马赫数Mi，当前合成攻角系数Kai，以及按步骤S2获得的合成攻角系数表，按侧滑角为AOS0时的攻角系数Kai，得到攻角值AOA1；
步骤S607：根据指示马赫数Mi，当前合成侧滑角系数Kbi，以及按步骤S3获得的合成侧滑角系数表，按攻角为AOA1时的侧滑角系数Kbi，得到侧滑角值AOS1；
步骤S608：根据指示马赫数Mi，当前合成攻角系数Kai，以及按步骤S2获得的合成攻角系数表，按侧滑角为AOS1时的攻角系数Kai，得到攻角值AOA2；
步骤S609：根据指示马赫数Mi，当前合成侧滑角系数Kbi，以及按步骤S3获得的合成侧滑角系数表，按攻角为AOA2时的侧滑角系数Kbi，得到侧滑角值AOS2；
步骤S610：根据指示马赫数Mi，当前合成攻角系数Kai，以及按步骤S2获得的合成攻角系数表，按侧滑角为AOS2时的攻角系数Kai，得到攻角值AOA3；
步骤S611：根据指示马赫数Mi，当前合成侧滑角系数Kbi，以及按步骤S3获得的合成侧滑角系数表，按攻角为AOA3时的侧滑角系数Kbi，得到侧滑角值AOS3；
步骤S612：AOA3即为真实攻角，AOS3即为真实侧滑角。
为了更好地实现本方案，进一步地，所述步骤S7中的真实动压的解算及旋翼下洗流影响修正包括以下步骤：
步骤S701：根据指示马赫数Mi、真实攻角值AOA3和真实侧滑角值AOS3，以及按步骤S1获得的动压修正系数表Kp，经计算得到对应的动压修正系数Kp0；
步骤S702：根据指示动压Qci、动压修正系数Kp0，计算修正后的动压Qc0：Qc0＝Qci/(1+Kp0)；
步骤S703：根据步骤S601中的公式，根据静压Psi和修正后的动压Qc0计算马赫数M0；
步骤S704：根据马赫数M0、真实攻角AOA3、真实侧滑角AOS3，以及按步骤S1获得的动压修正系数表Kp，经计算得到对应的动压修正系数Kp1；
步骤S705：根据指示动压Qci和动压修正系数Kp1，根据步骤S702计算得到修正后的动压Qc1；
步骤S706：重复步骤S703-S705两次，得到修正后的真实动压Qc，完成旋翼下洗流影响修正；
步骤S707：根据静压Psi、真实动压Qc完成大气参数解算。
为了更好地实现本方案，进一步地，所述步骤S706具体为：
步骤S7061：根据步骤S601中的公式，根据静压Psi和动压Qc1计算马赫数M1；
步骤S7062：根据马赫数M1、攻角AOA3、侧滑角AOS3，以及按步骤S1获得的动压修正系数表Kp，经计算得到对应的动压修正系数Kp2；
步骤S7063：根据动压Qci和动压修正系数Kp2，根据步骤S702计算得到修正后的动压Qc2；
步骤S7064：根据步骤S601中的公式，根据静压Psi和动压Qc2计算马赫数M2；
步骤S7065：根据马赫数M2、攻角AOA3、侧滑角AOS3，以及按步骤S1获得的动压修正系数表Kp，经计算得到对应的动压修正系数Kp3；
步骤S7066：根据动压Qci和动压修正系数Kp3，根据步骤S702计算得到修正后的真实动压Qc3；
以真实动压Qc3作为真实动压Qc。
本方案通过重复数次代入数据，多次数学推导，降低了修正方法的误差率，实现了旋翼下洗流影响修正，并且方法简单易实现，成本低，具有很强的实际应用价值和经济效益。
综上所述，由于采用了上述技术方案，本发明的有益效果是：
1.本发明所述的一种直升机大气数据系统旋翼下洗流影响修正方法，利用随动探头实现总压、静压感受，使总压始终对准合成气流方向，保证了静压感受的准确性，从而使旋翼下洗流修正转化为动压修正，简化了旋翼下洗流修正模型，易于实现；
2.本发明所述的一种直升机大气数据系统旋翼下洗流影响修正方法，经过大量气动流场仿真(CFD)、数学推导和试验所得，实现了旋翼下洗流修影响修正，保证了大气数据在直升机低速飞行时的测量精度，本方案成本低、易实现；
3.本发明所述的一种直升机大气数据系统旋翼下洗流影响修正方法，同时实现攻角、侧滑角测量功能，简化了大气数据系统的系统组成，提高了系统的集成度和综合化水平。
附图说明
为了更清楚地说明本技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图，其中：
图1是本发明的随动探头式大气数据系统的组成结构图；
图2是本发明的随动探头直杆开孔横截面图；
图3是本发明的随动探头式大气数据系统的工作原理框图；
图4是本发明的旋翼下洗流修正算法的流程框图；
图5是本发明的方法流程框图；
图中，1-全压口，2-静压孔，3-随动探头，4-随动探头直杆，5-压差孔，6-静温传感器，7-安装支臂。
具体实施方式
为了更清楚地说明本发明实施例的技术方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，应当理解，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例，因此不应被看作是对保护范围的限定。基于本发明中的实施例，本领域普通技术工作人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“设置”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；也可以是直接相连，也可以是通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。
下面结合图1至图5对本发明作详细说明。
实施例1：
一种直升机大气数据系统旋翼下洗流影响修正方法，包括以下步骤，
步骤S1：建立动压修正系数表Kp；
步骤S2：建立合成攻角系数表Kai；
步骤S3：建立合成侧滑角系数表Kbi；
步骤S4：采用随动探头测得总压和静压；
步骤S5：采用随动探头直杆上的测压孔测得上压力、下压力、左压力和右压力；
步骤S6：通过总压、静压、上压力和下压力的差、左压力和右压力的差解算出攻角和侧滑角；
步骤S7：通过攻角、侧滑角和马赫数得到动压修正系数Kp，实现真实动压的解算，完成旋翼下洗流修正。
工作原理：如图1所示，本实例所采用总静压探头为随动探头，可保证合成总压始终对准气流方向；随动探头直杆正交开4个测压孔，分别实现上压力、下压力、左压力及下压力的感受；安装支臂内的电子部件负责动压、静压、上下压差、左右压差的采集及旋翼下洗流修正，并完成大气参数的解算及输出。
如图2所示，随动探头直杆开孔数量可为4个或6个，开孔4个时，测压孔为上、下、左、右正交分布，开孔6个时，测压孔间隔60°均匀分布，并包含上、下两个固定位置，在本方案中，一般采用开孔数量为4个的随动探头直杆。
如图3所示，安装支臂内的电子部件采集随动探头来的总压、静压，采集探头直杆来的上压力、下压力、左压力和下压力，采集静温传感器来的静温电阻信号，建立旋翼下洗流与上述压力之间的关系模型，完成旋翼下洗流修正后实现大气参数的解算输出。
如图4所示，大气数据系统上电后按程序周期执行周期性任务。
实施例2：
一种直升机大气数据系统旋翼下洗流影响修正方法，如图5，包括以下步骤，
步骤S1：建立动压修正系数表Kp；
步骤S2：建立合成攻角系数表Kai；
步骤S3：建立合成侧滑角系数表Kbi；
步骤S4：采用随动探头测得总压和静压；
步骤S5：采用随动探头直杆上的测压孔测得上压力、下压力、左压力和右压力；
步骤S6：通过总压、静压、上压力和下压力的差、左压力和右压力的差解算出攻角和侧滑角；
步骤S7：通过攻角、侧滑角和马赫数得到动压修正系数Kp，实现真实动压的解算，完成旋翼下洗流修正。
所述步骤S1中动压修正系数表Kp的建立，应根据旋翼开启或关闭条件，通过气动流场仿真或气动流场试验，建立不同马赫数、不同攻角、不同侧滑角下的动压修正系数Kp：
其中，Qci为指示动压：开启旋翼时的动压测量值；Qc为真实动压：关闭旋翼时的动压测量值。
所述步骤S2中合成攻角系数表Kai的建立，应在打开旋翼的情况下，通过气动流场仿真或气动流场试验，建立不同马赫数、不同攻角、不同侧滑角下的合成攻角系数Kai：
所述P(d-u)i为随动探头直杆上的下压力与上压力的差值，Qci为指示动压。
所述步骤S3中合成侧滑角系数表Kbi的建立，应在打开旋翼的情况下，通过气动流场仿真或气动流场试验，建立不同马赫数、不同攻角下的合成侧滑角系数Kbi：
所述P(l-r)i为随动探头直杆上的左压力与右压力的差值，Qci为指示动压。
为了更好地实现本方案，进一步地，所述步骤S6重复迭代多次，以获得要求的攻角、侧滑角计算精度。
为了更好地实现本方案，进一步地，所述步骤S7重复迭代多次，以获得要求的大气参数计算精度。
所述步骤S6中的真实攻角、真实侧滑角的解算包括以下步骤：
步骤S601：根据静压Psi、动压Qci计算指示马赫数Mi：
其中k为绝热指数；
步骤S602：根据动压Qci、P(d-u)i计算当前状态下的合成攻角系数Kai：
其中P(d-u)i为随动探头直杆上的下压力与上压力的差值；
步骤S603：根据动压Qci、P(l-r)i计算当前状态下的合成侧滑角系数Kbi：
其中P(l-r)i为随动探头直杆上的左压力与右压力的差值；
步骤S604：根据指示马赫数Mi，当前合成攻角系数Kai，以及按步骤S2获得的合成攻角系数表，按侧滑角为0度时的合成攻角系数Kai，得到侧滑角为0度时的攻角值AOA0；
步骤S605：根据指示马赫数Mi，当前合成侧滑角系数Kbi，以及按步骤S3获得的合成侧滑角系数表，按攻角为AOA0时的合成侧滑角系数Kbi，得到攻角为AOA0时的侧滑角值AOS0；
步骤S606：根据指示马赫数Mi，当前合成攻角系数Kai，以及按步骤S2获得的合成攻角系数表，按侧滑角为AOS0时的攻角系数Kai，得到攻角值AOA1；
步骤S607：根据指示马赫数Mi，当前合成侧滑角系数Kbi，以及按步骤S3获得的合成侧滑角系数表，按攻角为AOA1时的侧滑角系数Kbi，得到侧滑角值AOS1；
步骤S608：根据指示马赫数Mi，当前合成攻角系数Kai，以及按步骤S2获得的合成攻角系数表，按侧滑角为AOS1时的攻角系数Kai，得到攻角值AOA2；
步骤S609：根据指示马赫数Mi，当前合成侧滑角系数Kbi，以及按步骤S3获得的合成侧滑角系数表，按攻角为AOA2时的侧滑角系数Kbi，得到侧滑角值AOS2；
步骤S610：根据指示马赫数Mi，当前合成攻角系数Kai，以及按步骤S2获得的合成攻角系数表，按侧滑角为AOS2时的攻角系数Kai，得到攻角值AOA3；
步骤S611：根据指示马赫数Mi，当前合成侧滑角系数Kbi，以及按步骤S3获得的合成侧滑角系数表，按攻角为AOA3时的侧滑角系数Kbi，得到侧滑角值AOS3；
步骤S612：AOA3即为真实攻角，AOS3即为真实侧滑角。
所述步骤S7中的真实动压的解算及旋翼下洗流修正包括以下步骤：
步骤S701：根据指示马赫数Mi、真实攻角值AOA3和真实侧滑角值AOS3，以及按步骤S1获得的动压修正系数表Kp，经计算得到对应的动压修正系数Kp0；
步骤S702：根据指示动压Qci、动压修正系数Kp0，计算修正后的动压Qc0：Qc0＝Qci/(1+Kp0)；
步骤S703：根据步骤S601中的公式，根据静压Psi和修正后的动压Qc0计算马赫数M0；
步骤S704：根据马赫数M0、真实攻角AOA3、真实侧滑角AOS3，以及按步骤S1获得的动压修正系数表Kp，经计算得到对应的动压修正系数Kp1；
步骤S705：根据指示动压Qci和动压修正系数Kp1，根据步骤S702计算得到修正后的动压Qc1；
步骤S706：重复步骤S703-S705两次，得到修正后的真实动压Qc，完成旋翼下洗流影响修正；
步骤S707：根据静压Psi、真实动压Qc完成大气参数解算。
工作原理：本方案通过重复数次代入数据，多次数学推导，降低了修正方法的误差率，实现了旋翼下洗流修正，并且方法简单易实现，成本低，具有很强的实际应用价值和经济效益。
本实施例的其他部分与上述实施例2相同，故不再赘述。
实施例3：
本实施例是实施例2的进一步补充说明，所述步骤S706具体为：
步骤S7061：根据步骤S601中的公式，根据静压Psi和动压Qc1计算马赫数M1；
步骤S7062：根据马赫数M1、攻角AOA3、侧滑角AOS3，以及按步骤S1获得的动压修正系数表Kp，经计算得到对应的动压修正系数Kp2；
步骤S7063：根据动压Qci和动压修正系数Kp2，根据步骤S702计算得到修正后的动压Qc2；
步骤S7064：根据步骤S601中的公式，根据静压Psi和动压Qc2计算马赫数M2；
步骤S7065：根据马赫数M2、攻角AOA3、侧滑角AOS3，以及按步骤S1获得的动压修正系数表Kp，经计算得到对应的动压修正系数Kp3；
步骤S7066：根据动压Qci和动压修正系数Kp3，根据步骤S702计算得到修正后的真实动压Qc3；
以真实动压Qc3作为真实动压Qc。
本实施例的其他部分与上述实施例2相同，故不再赘述。
以上所述，仅是本发明的较佳实施例，并非对本发明做任何形式上的限制，凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化，均落入本发明的保护范围之内。
Claims (9)
1.一种直升机大气数据系统旋翼下洗流影响修正方法，其特征在于：包括以下步骤，
步骤S1：建立不同马赫数、不同攻角、不同侧滑角下的动压修正系数表Kp；
步骤S2：建立不同马赫数、不同攻角、不同侧滑角下的合成攻角系数表Kai；
步骤S3：建立不同马赫数、不同攻角、不同侧滑角下的合成侧滑角系数表Kbi；
步骤S4：采用随动探头测得总压和静压；
步骤S5：采用随动探头直杆上的测压孔测得上压力、下压力、左压力和右压力；
步骤S6：通过总压、静压、上压力和下压力的差、左压力和右压力的差解算出攻角和侧滑角；
步骤S7：通过攻角、侧滑角和马赫数得到动压修正系数Kp，实现真实动压的解算，完成旋翼下洗流修正。
2.根据权利要求1所述的一种直升机大气数据系统旋翼下洗流影响修正方法，其特征在于：所述步骤S1，应根据旋翼开启或关闭条件，通过气动流场仿真或风洞流场试验，建立不同马赫数、不同攻角、不同侧滑角下的动压修正系数Kp：
其中，Qci为指示动压：开启旋翼时的动压测量值；Qc为真实动压：关闭旋翼时的动压测量值。
3.根据权利要求1所述的一种直升机大气数据系统旋翼下洗流影响修正方法，其特征在于：所述步骤S2，需在旋翼开启的情况下，通过气动流场仿真或风洞流场试验，建立不同马赫数、不同攻角、不同侧滑角下的合成攻角系数Kai：
所述P(d-u)i为随动探头直杆上的下压力与上压力的差值，Qci为指示动压。
4.根据权利要求1所述的一种直升机大气数据系统旋翼下洗流影响修正方法，其特征在于：所述步骤S3，需在旋翼开启的情况下，通过气动流场仿真或风洞流场试验，建立不同马赫数、不同攻角下的合成侧滑角系数Kbi：
所述P(l-r)i为随动探头直杆上的左压力与右压力的差值，Qci为指示动压。
5.根据权利要求1所述的一种直升机大气数据系统旋翼下洗流影响修正方法，其特征在于：所述步骤S6重复迭代多次，以获得要求的攻角、侧滑角计算精度。
6.根据权利要求1所述的一种直升机大气数据系统旋翼下洗流影响修正方法，其特征在于：所述步骤S7重复迭代多次，以获得要求的大气参数计算精度。
7.根据权利要求1所述的一种直升机大气数据系统旋翼下洗流影响修正方法，其特征在于：所述步骤S6中的真实攻角、真实侧滑角的解算包括以下步骤：
步骤S601：根据静压Psi、指示动压Qci计算指示马赫数Mi：
其中k为绝热指数；
步骤S602：根据P(d-u)i、指示动压Qci计算合成攻角系数Kai：
其中P(d-u)i为随动探头直杆上的下压力与上压力的差值；
步骤S603：根据P(l-r)i、指示动压Qci计算合成侧滑角系数Kbi：
其中P(l-r)i为随动探头直杆上的左压力与右压力的差值；
步骤S604：根据指示马赫数Mi，当前合成攻角系数Kai，以及按步骤S2获得的合成
攻角系数表，按侧滑角为0度时的合成攻角系数Kai，得到侧滑角为0度时的攻角值AOA0；
步骤S605：根据指示马赫数Mi，当前合成侧滑角系数Kbi，以及按步骤S3获得的合
成侧滑角系数表，按攻角为AOA0时的合成侧滑角系数Kbi，得到攻角为AOA0时的侧滑角值AOS0；
步骤S606：根据指示马赫数Mi，当前合成攻角系数Kai，以及按步骤S2获得的合成攻角系数表，按侧滑角为AOS0时的攻角系数Kai，得到攻角值AOA1；
步骤S607：根据指示马赫数Mi，当前合成侧滑角系数Kbi，以及按步骤S3获得的合成侧滑角系数表，按攻角为AOA1时的侧滑角系数Kbi，得到侧滑角值AOS1；
步骤S608：根据指示马赫数Mi，当前合成攻角系数Kai，以及按步骤S2获得的合成攻角系数表，按侧滑角为AOS1时的攻角系数Kai，得到攻角值AOA2；
步骤S609：根据指示马赫数Mi，当前合成侧滑角系数Kbi，以及按步骤S3获得的合成侧滑角系数表，按攻角为AOA2时的侧滑角系数Kbi，得到侧滑角值AOS2；
步骤S610：根据指示马赫数Mi、合成攻角系数Kai，以及按步骤S2获得的合成攻角系数表，按侧滑角为AOS2时的攻角系数，得到攻角值AOA3；
步骤S611：根据指示马赫数Mi、合成侧滑角系数Kbi，以及按步骤S3获得的合成侧滑角系数表，按攻角为AOA3时的侧滑角系数，得到侧滑角值AOS3；
步骤S612：将AOA3作为真实攻角，将AOS3作为真实侧滑角。
8.根据权利要求7所述的一种直升机大气数据系统旋翼下洗流影响修正方法，其特征在于：所述步骤S7中的真实动压的解算及旋翼下洗流修正包括以下步骤：
步骤S701：根据指示马赫数Mi、真实攻角值AOA3和真实侧滑角值AOS3，以及步骤S1获得的动压修正系数表，得到对应的动压修正系数Kp0；
步骤S702：根据指示动压Qci、动压修正系数Kp0计算修正后的动压Qc0：Qc0＝Qci/(1+Kp0)；
步骤S703：根据步骤S601中的公式，根据静压Psi和修正后的动压Qc0计算马赫数M0；
步骤S704：根据马赫数M0、攻角AOA3、侧滑角AOS3，以及步骤S1获得的动压修正系数表，得到对应的动压修正系数Kp1；
步骤S705：根据动压Qci和动压修正系数Kp1，根据步骤S702计算得到修正后的动压Qc1；
步骤S706：重复步骤S703-S705两次，得到修正后的真实动压Qc，完成旋翼下洗流影响修正；
步骤S707：根据静压Psi、真实动压Qc完成大气参数解算。
9.根据权利要求8所述的一种直升机大气数据系统旋翼下洗流影响修正方法，其特征在于：所述步骤S706具体为：
步骤S7061：根据步骤S601中的公式，根据静压Psi和修正后的动压Qc1计算马赫数M1；
步骤S7062：根据马赫数M1、攻角AOA3、侧滑角AOS3，以及步骤S1获得的动压修正系数表，得到对应的动压修正系数Kp2；
步骤S7063：根据动压Qci和动压修正系数Kp2，根据步骤S702计算得到修正后的动压Qc2；
步骤S7064：根据步骤S601中的公式，根据静压Psi和动压Qc2计算马赫数M2；
步骤S7065：根据马赫数M2、攻角AOA3、侧滑角AOS3，以及步骤S1获得的动压修正系数表，得到对应的动压修正系数Kp3；
步骤S7066：根据动压Qci和动压修正系数Kp3，根据步骤S702计算得到修正后的动压Qc3；
以修正后的动压Qc3作为真实动压Qc。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110210586.4A CN112977869B (zh) | 2021-02-25 | 2021-02-25 | 一种直升机大气数据系统旋翼下洗流影响修正方法 |
| PCT/CN2021/140488 WO2022179278A1 (zh) | 2021-02-25 | 2021-12-22 | 一种直升机大气数据系统旋翼下洗流影响修正方法 |
| LU501425A LU501425B1 (en) | 2021-02-25 | 2021-12-22 | Method for correcting affection of rotor downwash with helicopter air data system |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110210586.4A CN112977869B (zh) | 2021-02-25 | 2021-02-25 | 一种直升机大气数据系统旋翼下洗流影响修正方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112977869A CN112977869A (zh) | 2021-06-18 |
| CN112977869B true CN112977869B (zh) | 2022-11-01 |
Family
ID=76350566
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110210586.4A Active CN112977869B (zh) | 2021-02-25 | 2021-02-25 | 一种直升机大气数据系统旋翼下洗流影响修正方法 |
Country Status (3)
| Country | Link |
|---|---|
| CN (1) | CN112977869B (zh) |
| LU (1) | LU501425B1 (zh) |
| WO (1) | WO2022179278A1 (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113624393B (zh) * | 2021-06-22 | 2023-07-14 | 成都凯天电子股份有限公司 | 具有旋翼下洗流影响修正的大气数据系统及方法 |
| CN114580219B (zh) * | 2022-05-07 | 2022-09-09 | 成都凯天电子股份有限公司 | 一种分布式大气数据系统参数校准的方法 |
| CN114969987B (zh) * | 2022-07-18 | 2022-11-25 | 成都凯天电子股份有限公司 | 一种基于气动仿真的l型多功能探头安装位置选定方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5257536A (en) * | 1990-08-03 | 1993-11-02 | Sextant Avionique | Aircraft probe for measuring aerodynamic parameters of amient flow |
| EP1462806A1 (en) * | 2003-03-25 | 2004-09-29 | Rosemount Aerospace Inc. | Determination of airspeed from negative impact pressure |
| CN104061960A (zh) * | 2014-05-26 | 2014-09-24 | 中国航天空气动力技术研究院 | 一种亚音速飞行器体上气压高度参数确定方法 |
| CN108981789A (zh) * | 2018-09-26 | 2018-12-11 | 武汉科技大学 | 一种基于piv系统的无人机喷洒下洗流场测量装置及方法 |
| CN110155363A (zh) * | 2019-03-21 | 2019-08-23 | 北京机电工程研究所 | 基于cfd方法的弹性气动数据精确获取方法 |
| CN112163271A (zh) * | 2020-09-04 | 2021-01-01 | 北京空天技术研究所 | 大气数据传感系统的大气参数解算方法 |
Family Cites Families (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5811691A (en) * | 1997-12-26 | 1998-09-22 | Sikorsky Aircraft Corporation | Blade-mounted total pressure probe for a rotating blade |
| US20030101798A1 (en) * | 2001-12-05 | 2003-06-05 | The Entwistle Company | Helicopter hollow blade pressure check and fill apparatus and method to use same |
| RU2307357C1 (ru) * | 2005-12-07 | 2007-09-27 | Казанский государственный технический университет им. А.Н. Туполева | Способ измерения воздушных сигналов вертолета и система для его осуществления |
| WO2012105973A1 (en) * | 2011-02-02 | 2012-08-09 | Michigan Aerospace Corporation | Atmospheric measurement system and method |
| EP2729927A4 (en) * | 2011-10-06 | 2014-07-30 | Cae Inc | METHOD FOR DEVELOPING A MATHEMATIC DYNAMIC MODEL OF A VEHICLE FOR USE IN A COMPUTER-CONTROLLED VEHICLE SIMULATOR |
| CN102607639A (zh) * | 2012-02-24 | 2012-07-25 | 南京航空航天大学 | 基于bp神经网络的大攻角飞行状态下大气数据测量方法 |
| US8774988B2 (en) * | 2012-08-21 | 2014-07-08 | Bell Helicopter Textron Inc. | Aircraft environmental sensors and system |
| US9506945B2 (en) * | 2014-06-10 | 2016-11-29 | Sikorsky Aircraft Corporation | Rotorcraft flight parameter estimation |
| CN104374408A (zh) * | 2014-11-27 | 2015-02-25 | 江西洪都航空工业集团有限责任公司 | 一种侧滑角传感器的侧滑角修正的计算方法 |
| CN104374541A (zh) * | 2014-11-27 | 2015-02-25 | 江西洪都航空工业集团有限责任公司 | 一种l形压力受感器的静压计算方法 |
| CN104376225A (zh) * | 2014-11-27 | 2015-02-25 | 江西洪都航空工业集团有限责任公司 | 一种风标式攻角传感器的攻角修正计算方法 |
| CN105136196B (zh) * | 2015-07-27 | 2017-11-03 | 江西洪都航空工业集团有限责任公司 | 一种采用多功能探头的分布式大气数据系统 |
| CN106248139B (zh) * | 2016-07-29 | 2023-08-11 | 成都凯天电子股份有限公司 | 大气数据测量探头 |
| CN106705996B (zh) * | 2016-11-25 | 2019-11-12 | 北京航天自动控制研究所 | 一种基于大气特征参数的飞行器导航信息修正方法 |
| CN107132376A (zh) * | 2017-04-21 | 2017-09-05 | 陕西飞机工业（集团）有限公司 | 一种飞机攻角修正曲线的获取方法 |
| CN207763881U (zh) * | 2018-01-12 | 2018-08-24 | 太原航空仪表有限公司 | 一种用于测量旋翼下洗气流的多功能探头 |
| US11035877B2 (en) * | 2018-09-12 | 2021-06-15 | Bell Helicopter Textron Inc. | Pitot-static system blockage detector |
| CN109738009A (zh) * | 2018-12-07 | 2019-05-10 | 武汉航空仪表有限责任公司 | 一种风标型多功能探头 |
| US11828771B2 (en) * | 2019-07-15 | 2023-11-28 | The Boeing Company | Method and system for collecting air data using a laser-induced plasma channel |
| CN110927744B (zh) * | 2019-11-22 | 2023-04-14 | 成都凯天电子股份有限公司 | 直升机光学大气数据系统 |
-
2021
- 2021-02-25 CN CN202110210586.4A patent/CN112977869B/zh active Active
- 2021-12-22 LU LU501425A patent/LU501425B1/en active IP Right Grant
- 2021-12-22 WO PCT/CN2021/140488 patent/WO2022179278A1/zh active Application Filing
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5257536A (en) * | 1990-08-03 | 1993-11-02 | Sextant Avionique | Aircraft probe for measuring aerodynamic parameters of amient flow |
| EP1462806A1 (en) * | 2003-03-25 | 2004-09-29 | Rosemount Aerospace Inc. | Determination of airspeed from negative impact pressure |
| CN104061960A (zh) * | 2014-05-26 | 2014-09-24 | 中国航天空气动力技术研究院 | 一种亚音速飞行器体上气压高度参数确定方法 |
| CN108981789A (zh) * | 2018-09-26 | 2018-12-11 | 武汉科技大学 | 一种基于piv系统的无人机喷洒下洗流场测量装置及方法 |
| CN110155363A (zh) * | 2019-03-21 | 2019-08-23 | 北京机电工程研究所 | 基于cfd方法的弹性气动数据精确获取方法 |
| CN112163271A (zh) * | 2020-09-04 | 2021-01-01 | 北京空天技术研究所 | 大气数据传感系统的大气参数解算方法 |
Non-Patent Citations (1)
| Title |
|---|
| 直升机大气数据计算机综合测试系统的设计;陈海牛;《测控技术》;20201218;全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| LU501425A1 (en) | 2022-09-16 |
| LU501425B1 (en) | 2023-01-23 |
| CN112977869A (zh) | 2021-06-18 |
| WO2022179278A1 (zh) | 2022-09-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112977869B (zh) | 一种直升机大气数据系统旋翼下洗流影响修正方法 | |
| CN106200629B (zh) | 一种无人机飞行控制系统的故障可检测度分析方法 | |
| US5214596A (en) | System for determining the airspeed of helicopters | |
| CN101321667A (zh) | 用于重构飞行器、尤其是客机上的阵风和结构载荷的方法 | |
| CN111241625A (zh) | 飞机发动机特性辨识及平衡极曲线辨识的试飞方法 | |
| CN114004023A (zh) | 一种基于循环神经网络的飞行器气动参数辨识的方法 | |
| CN105527455B (zh) | 用于自动地估计飞行器的至少一个速度的方法及装置 | |
| CN104102127A (zh) | 一种机载气动参数辨识系统 | |
| CN106444885A (zh) | 一种颤振主动抑制控制器构成及其模拟方法 | |
| Lum et al. | Design and experiment of data-driven modeling and flutter control of a prototype wing | |
| CN111412887B (zh) | 一种基于卡尔曼滤波的攻角、侧滑角辨识方法 | |
| Samy et al. | Unmanned air vehicle air data estimation using a matrix of pressure sensors: a comparison of neural networks and look-up tables | |
| CN110027728B (zh) | 通过空中飞行试验辨识飞机气动焦点的方法 | |
| CN113525711B (zh) | 通过飞行试验辨识飞机气动焦点的方法 | |
| CN215893640U (zh) | 具有旋翼下洗流影响修正的随动探头系统 | |
| CN110763424B (zh) | 一种机翼表面压力测量方法、系统及装置 | |
| Jurisson et al. | Aeroservoelastic flight testing platform development for system identification | |
| Li et al. | A Compact Embedded Flight Parameter Detection System for Small Soaring UAVs | |
| CN105129073B (zh) | 一种基于gofir的飞行器阵风减缓自适应前馈控制方法 | |
| Wong et al. | Unmanned aerial vehicle flight model validation using on-board sensing and instrumentation | |
| CN113609771B (zh) | 一种基于mbse的飞机大气参数模型计算方法 | |
| CN113624393B (zh) | 具有旋翼下洗流影响修正的大气数据系统及方法 | |
| CN107038312A (zh) | 一种低压降孔单元的流量计算方法 | |
| Coiro et al. | Numerical, wind tunnel and flight tests for P92J and P96 light aircraft | |
| CN105468870A (zh) | 一种小型螺旋翼的升力的计算方法 |
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 |