[發明專利]飛行器飛行參數解算方法

有效

• 申請號: 201910292072.0
• 申請日: 2019-04-12
• 公開（公告）號: CN110059396B
• 公開（公告）日: 2023-06-16
• 發明（設計）人: 黃俊森;戴梧葉;盧志毅;李超
• 申請（專利權）人: 北京空天技術研究所
• 主分類號: G06F30/20
• 分類號: G06F30/20;G06F30/15;G06F30/17
• 代理公司: 暫無信息
• 代理人: 暫無信息

1.一種飛行器飛行參數解算方法，其特征在于，該方法包括：

基于慣導提供的飛行器的馬赫數M和飛行器上設置的測點的靜壓計算來流靜壓p_∞；

基于所述來流靜壓p_∞計算虛擬總壓p₀；

基于所述虛擬總壓p₀和測點的靜壓解算得到真實飛行器飛行參數；

通過下式基于慣導提供的飛行器的馬赫數M和飛行器上設置的測點的靜壓計算來流靜壓p_∞：

p_∞＝f₁(M,P1,P3)*(p_PD1+p_PU1)，

其中，f₁為以M,P1,P3為自變量的多項式函數，P1,P3均為無量綱壓力，p_PD1為飛行器頭部下表面對稱面上沿飛行器軸線前部布置的測點PD1的靜壓，p_PU1為飛行器頭部上表面對稱面上沿飛行器軸線前部布置的測點PU1的靜壓，p_PR和p_PL為飛行器頭部側面線上沿飛行器軸線同一位置左右對稱地布置的兩個測點的靜壓，多項式函數的形式統一為：

F＝f(X,Y,Z)＝∑c_i,j,kXⁱY^jZ^k，

其中，各自變量的最高冪次(i,j,k)不超過4，各項系數c_i,j,k由最小二乘法給出。

2.根據權利要求1所述的方法，其特征在于，通過下式基于所述來流靜壓p_∞計算虛擬總壓p₀：

其中，γ為空氣的比熱比。

3.根據權利要求2所述的方法，其特征在于，真實飛行器飛行參數包括真實攻角、真實側滑角、真實馬赫數和真實來流靜壓。

4.根據權利要求3所述的方法，其特征在于，通過下式基于所述虛擬總壓p₀和測點的靜壓解算得到真實馬赫數M^*：

M^*＝g₁(P4',P1',P3')，

其中，g₁為以P1',P3',P4'為自變量的多項式函數，P1',P3',P4'為以虛擬總壓無量綱化的壓力系數。

5.根據權利要求4所述的方法，其特征在于，通過下式基于所述虛擬總壓p₀和測點的靜壓解算得到真實側滑角β^*：

β^*＝g₃(P4',P1',P3')，

其中，g₃為以P1',P3',P4'為自變量的多項式函數。

6.根據權利要求5所述的方法，其特征在于，通過下式基于所述虛擬總壓p₀和測點的靜壓解算得到真實來流靜壓p_∞^*：

p_∞^*＝g₄(P4',P1',P3')*(p_PD1+p_PU1)，

其中，g₄為以P1',P3',P4'為自變量的多項式函數。

7.根據權利要求6所述的方法，其特征在于，通過下式基于所述虛擬總壓p₀和測點的靜壓解算得到真實攻角α₂^*：

α₂^*＝g₅(M^*,P2',β^*)，

其中，g₅為以M^*,P2',β^*為自變量的多項式函數，α₂^*對應于飛行器頭部上表面對稱面上沿飛行器軸線后部布置的測點PU2與飛行器頭部下表面對稱面上沿飛行器軸線后部布置的測點PD2處的平均當地攻角，P2'為測點PU2和測點PD2處的測點無量綱壓力。