CN106650077B - 一种弹性飞机尾涡遭遇动响应分析方法 - Google Patents

一种弹性飞机尾涡遭遇动响应分析方法 Download PDF

Info

Publication number
CN106650077B
CN106650077B CN201611161876.XA CN201611161876A CN106650077B CN 106650077 B CN106650077 B CN 106650077B CN 201611161876 A CN201611161876 A CN 201611161876A CN 106650077 B CN106650077 B CN 106650077B
Authority
CN
China
Prior art keywords
wake vortex
aircraft
dynamic response
elastic
wake
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
CN201611161876.XA
Other languages
English (en)
Other versions
CN106650077A (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.)
Xian Aircraft Design and Research Institute of AVIC
Original Assignee
Xian Aircraft Design and Research Institute of AVIC
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 Xian Aircraft Design and Research Institute of AVIC filed Critical Xian Aircraft Design and Research Institute of AVIC
Priority to CN201611161876.XA priority Critical patent/CN106650077B/zh
Publication of CN106650077A publication Critical patent/CN106650077A/zh
Application granted granted Critical
Publication of CN106650077B publication Critical patent/CN106650077B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/15—Vehicle, aircraft or watercraft design
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • 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

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明属于气动弹性技术领域，涉及一种弹性飞机尾涡遭遇动响应分析方法。步骤一，根据以下公式计算弹性飞机遭遇的尾涡强度
Figure DDA0001181782420000011
步骤二，计算弹性飞机遭遇的反旋转尾涡系诱导速度；步骤三，求解弹性飞机遭遇的尾涡非定常气动力；步骤四，采用频域方法求解弹性飞机尾涡遭遇频域动响应；步骤五，根据步骤四中求解得到的弹性飞机尾涡遭遇频域动响应，采用傅里叶反变换求解弹性飞机尾涡遭遇时域动响应；步骤六，根据步骤五中得到的弹性飞机尾涡遭遇时域动响应，采用模态位移法计算穿越飞机结构上的载荷动响应。建立了一种弹性飞机尾涡遭遇动响应分析方法，为国内飞机设计中的尾涡遭遇分析提供手段。

Description

一种弹性飞机尾涡遭遇动响应分析方法
技术领域
本发明属于气动弹性技术领域，涉及一种弹性飞机尾涡遭遇动响应分析方法。
背景技术
尾涡遭遇直接影响军用飞机空中加油、编队飞行和机动飞行等，同时也是制约空港流量和危害民航飞机飞行安全的主要因素。国外一直在开展尾涡气动特性和飞机尾涡遭遇安全性方面的研究，从2007年开始，以色列的Karpel等人、德宇航和空客的研究者开始研究弹性飞机尾涡遭遇动响应分析方法研究，并且相继推出了未公开发行的Dynresp、VarLoads和Gusto等分析软件。另外，国外已基于空客A400M运输机和瑞典鹰狮JAS-39战斗机开展了尾涡遭遇飞行试验研究。国内在研究尾涡对民航飞机飞行安全方面有一定的基础，但在弹性飞机尾涡遭遇动响应分析方面至今并未见到有公开的研究成果推出。
发明内容
本发明的目的：本发明旨在建立一种弹性飞机尾涡遭遇动响应分析方法，为国内飞机设计中的尾涡遭遇分析提供手段。
本发明的技术方案：一种弹性飞机尾涡遭遇动响应分析方法，其特征在于所述的分析方法包括如下步骤：
步骤一，根据以下公式计算弹性飞机遭遇的尾涡强度
Figure BDA0001181782400000011
其中：Γ为尾涡的强度；m为飞机质量；nz为飞机法向过载；b0＝sB，为左右尾涡涡核的间距；B为翼展；s为翼型参数；
步骤二，计算弹性飞机遭遇的反旋转尾涡系诱导速度；
步骤三，求解弹性飞机遭遇的尾涡非定常气动力；
步骤四，采用频域方法求解弹性飞机尾涡遭遇频域动响应；
步骤五，根据步骤四中求解得到的弹性飞机尾涡遭遇频域动响应，采用傅里叶反变换求解弹性飞机尾涡遭遇时域动响应；
步骤六，根据步骤五中得到的弹性飞机尾涡遭遇时域动响应，采用模态位移法计算穿越飞机结构上的载荷动响应。
优选地，步骤一中所述的翼型参数s与机翼形状有关，对于椭圆形机翼可取
Figure BDA0001181782400000021
对于后掠翼，可以取0.75～0.80。
优选地，基于Biot-Savart定律计算弹性飞机遭遇的反旋转尾涡系诱导速度。
优选地，基于离散阵风非定常气动力计算思路，将阵风模态列用频域尾涡诱导速度列代替，阵风强度用1代替，发展弹性飞机尾涡非定常气动力计算方法，根据该计算方法求解弹性飞机遭遇的尾涡非定常气动力。
本发明的有益效果：建立了一种弹性飞机尾涡遭遇动响应分析方法，并详细推导了运动方程，填补国内在尾涡遭遇动响应分析方法方面的空白；分析结果可直接用于动强度、系统载荷、操稳等专业领域；本发明为国内飞机型号尾涡遭遇情况下结构动响应和动载荷计算提供了一种有力手段，为缩小民航飞机飞行间隔、提高机场起降密度提供了一种分析方法，能有效支持空中加油机设计、战斗机编队飞行策略规划等工作。
附图说明
图1是飞机反旋转对称尾迹涡；
图2是弹性飞机尾涡遭遇动响应分析流程；
图3是穿越尾涡飞机表面气动网格上经历的诱导速度时域序列；
图4是穿越尾涡的飞机机翼表面法向诱导速度分布云图；
图5是穿越尾涡飞机表面法向诱导速度的傅里叶变换结果；
图6是穿越飞机弹性频域动响应计算流程；
图7是时域动响应恢复计算流程。
具体实施方式
下面结合附图对本发明作详细说明。
弹性飞机尾涡遭遇动响应分析方法的关键在于弹性飞机尾涡遭遇运动方程建立和求解，其核心是穿越尾涡飞机遭受的尾涡非定常气动力计算。
整个求解流程如图1所示，主要分为以下几步
(1)根据库塔-儒可夫斯基定理和亥姆霍兹定理推导尾涡产生飞机在其后方诱导出的反旋转对称尾涡系的强度；
(2)根据Biot-Savart定律，结合尾涡穿越飞机空间轨迹和飞行姿态角，确定尾涡穿越飞机表面的诱导速度；
(3)基于离散阵风非定常气动力计算公式，发展尾涡遭遇非定常气动力计算方法；
(4)基于机翼有限元模型、非定常气动力模型，并结合第(3)步建立的尾涡非定常气动力模型建立弹性飞机尾涡遭遇运动方程，最终通过频域方法进行动响应求解；
(5)对频域动响应进行时域转换，并结合结构载荷恢复矩阵，进行弹性结构时域动载荷分析。
实施例
尾涡产生飞机参数：重量27.6t，展长27.9m，巡航速度为152.7m/s，尾涡遭遇飞机飞行速度为100m/s，90°面内穿越，其它参数与尾涡产生飞机一致。
第一步，根据尾涡产生飞机的重量、过载、翼展、翼型参数等数据计算尾涡产生飞机在其后方拖出的尾涡的强度
Figure BDA0001181782400000041
式中：Γ为尾涡的强度，即尾涡环量；m为飞机质量；nz为飞机法向过载；b0＝sB，为左右尾涡涡核的间距；B为翼展；s为翼型参数，与机翼形状有关，对于椭圆形机翼可取
Figure BDA0001181782400000042
对于后掠翼，可以取0.75～0.80。
第二步，图1是飞机形成的反旋转尾涡系的示意图，根据Biot-Savart定律计算尾涡对空间任意一点的法向诱导速度
Figure BDA0001181782400000043
式中：
Figure BDA0001181782400000044
rc是涡核半径。
结合图2计算流程，进一步得到飞机在穿越尾涡场过程中，其表面气动网格遭受到的尾涡诱导速度时域序列{WW(t)}，具有如图3所示的数据存储形式，其各列具有典型的尾涡诱导速度形式。图4是不同时刻，穿越尾涡的飞机机翼表面法向诱导速度分布云图，可以看到，对于垂直穿越来说，机翼上同一展向位置的诱导速度相同，并且随着时间的推移，穿越尾涡的飞机机翼表面的诱导速度大小和方向都会发生变化。
第三步，对尾涡诱导速度时域序列进行傅里叶变换，得到{WW(iω)}。如图5所示，可以看出尾涡诱导速度在1～60Hz频率范围均有明显的分量。
将{WW(iω)}代入下式，计算穿越飞机受到的尾涡非定常气动力
{Fwakevortex}＝-q[φkh]T[SKJ][A][NIC]{WW(iω)}/V
式中：[NIC]是面元法向速度影响系数矩阵；[A]是气动面元面积对角矩阵；[SKJ]是插值矩阵，建立了结构节点模态矩阵和气动网格上某些点的模态矩阵的关系；V是穿越尾涡的飞机的飞行速度，q是飞行动压，[φkh]是穿越飞机气动网格上气动扰动点处的模态列。
第四步，根据尾涡穿越飞机结构动力学有限元模型和偶极子格网法，并结合第三步得到的尾涡诱导非定常气动力，推导得到弹性飞机尾涡遭遇气动弹性运动方程。
(-ω2[Mhh]+iω[Bhh]+[Khh]-q[Qhh(iω)]){ξ(iω)}＝＝{Fwakevortex}
第五步，方程频域求解，对频域动响应进行傅里叶反变换得到时域动响应，图6是模态广义位移时域响应和模态广义加速度时域响应。
第六步，结合穿越飞机结构载荷恢复矩阵，采用模态位移法计算载荷动响应。位移时域响应和过载时域响应。从图7中翼面气动网格上的物理位移动响应可以看出，尾涡引起的机翼弹性位移不大，这主要是由于平直的该型飞机机翼较为刚硬。从过载时域响应可以看出，尾涡引起的弹性过载最大可达到1.3，可见尾涡在穿越它的弹性飞机上产生了明显弹性过载，可能会影响飞机强度安全性，在设计中需要考虑。

Claims (2)

1.一种弹性飞机尾涡遭遇动响应分析方法，其特征在于所述的分析方法包括如下步骤：
步骤一，根据以下公式计算尾涡产生飞机在其后方诱导出的反旋转对称尾涡系的强度；
Figure FDA0003018588440000011
其中：Γ为尾涡的强度；m为飞机质量；nz为飞机法向过载；b0＝sB，为左右尾涡涡核的间距；B为翼展；s为翼型参数；
步骤二，根据Biot-Savart定律，结合尾涡穿越飞机空间轨迹和飞行姿态角，计算尾涡穿越飞机表面的诱导速度；；
步骤三，求解尾涡穿越飞机受到的尾涡非定常气动力；基于离散阵风非定常气动力计算思路，将阵风模态列用频域尾涡诱导速度列代替，阵风强度用1代替，发展弹性飞机尾涡非定常气动力计算方法，根据该计算方法求解尾涡穿越飞机受到的尾涡非定常气动力；
步骤四，根据尾涡穿越飞机结构动力学有限元模型和偶极子网格法，并根据步骤三中得到的尾涡非定常气动力建立弹性飞机尾涡遭遇气动弹性运动方程，采用频域方法求解弹性飞机尾涡遭遇频域动响应；
步骤五，根据步骤四中求解得到的弹性飞机尾涡遭遇频域动响应，采用傅里叶反变换求解弹性飞机尾涡遭遇时域动响应；
步骤六，根据步骤五中得到的弹性飞机尾涡遭遇时域动响应，采用模态位移法计算穿越飞机结构上的载荷动响应。
2.根据权利要求1所述的一种弹性飞机尾涡遭遇动响应分析方法，其特征为：步骤一中所述的翼型参数s与机翼形状有关，对于椭圆形机翼可取对于后掠翼，可以取0.75～0.80。
CN201611161876.XA 2016-12-15 2016-12-15 一种弹性飞机尾涡遭遇动响应分析方法 Active CN106650077B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611161876.XA CN106650077B (zh) 2016-12-15 2016-12-15 一种弹性飞机尾涡遭遇动响应分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611161876.XA CN106650077B (zh) 2016-12-15 2016-12-15 一种弹性飞机尾涡遭遇动响应分析方法

Publications (2)

Publication Number Publication Date
CN106650077A CN106650077A (zh) 2017-05-10
CN106650077B true CN106650077B (zh) 2021-07-16

Family

ID=58823408

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611161876.XA Active CN106650077B (zh) 2016-12-15 2016-12-15 一种弹性飞机尾涡遭遇动响应分析方法

Country Status (1)

Country Link
CN (1) CN106650077B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113837639B (zh) * 2021-09-29 2023-12-19 中国民用航空飞行学院 一种基于多维指标的航空器遭遇尾流风险评估系统
CN116500580B (zh) * 2023-06-28 2023-08-22 中国民航大学 一种基于多普勒雷达探测的涡核参数计算方法及系统

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102012953B (zh) * 2010-11-04 2013-05-08 西北工业大学 Cfd/csd耦合求解非线性气动弹性仿真方法
CN103077259A (zh) * 2011-10-26 2013-05-01 上海机电工程研究所 高超声速导弹多场耦合动力学一体化仿真分析方法
CN103577648B (zh) * 2013-11-13 2016-06-01 中国航空工业集团公司西安飞机设计研究所 运输类飞机货物空投时机翼结构载荷的确定方法
CN104732109A (zh) * 2015-04-09 2015-06-24 金陵科技学院 一种竖轴旋转发电装置的尾涡预报方法
CN104832370A (zh) * 2015-05-14 2015-08-12 河海大学 用于风轮涡线控制方程离散的三步三阶预估校正方法
CN105183975B (zh) * 2015-09-01 2018-07-03 中国空气动力研究与发展中心高速空气动力研究所 一种基于尾涡流场传递的多飞行器飞行编队数值模拟方法
CN105468814B (zh) * 2015-10-29 2018-11-09 南京航空航天大学 一种考虑全机气动干扰的直升机飞行特性计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
计及非定常气动力的飞机机动飞行动载荷预计;张积亭;《第十一届全国空气弹性学术交流会》;20091231;406-410 *
飞机T型尾翼颤振计算的若干关键问题;杨飞 等;《振动与冲击》;20111231;第30卷(第5期);136-139 *

Also Published As

Publication number Publication date
CN106650077A (zh) 2017-05-10

Similar Documents

Publication Publication Date Title
Murua et al. Assessment of wake-tail interference effects on the dynamics of flexible aircraft
Slotnick Computational aerodynamic analysis for the formation flight for aerodynamic benefit program
Sagaga et al. Acoustic predictions for the side-by-side air taxi rotor in hover
CN106650077B (zh) 一种弹性飞机尾涡遭遇动响应分析方法
Kuzmina et al. Integrated numerical and experimental investigations of the Active/Passive Aeroelastic concepts on the European Research Aeroelastic Model (EuRAM)
Sereez et al. Impact of ground effect on airplane lateral directional stability during take-off and landing
Paudel et al. Aerodynamic and stability analysis of blended wing body aircraft
Lei et al. Effects of wing kinematics on modulating odor plume structures in the odor tracking flight of fruit flies
Bludau et al. Validation of a dynamic inflow model based on a flight dynamics model and a Lattice-Boltzmann fluid solver using flight test data
Ahluwalia et al. CFD analysis on different shapes of winglet at low subsonic flow
Siddaramaiah et al. Preliminary Studies in the use of Folding Wing-tips for Loads Alleviation
Ghoreyshi et al. Aerodynamics modeling of a maneuvering aircraft using indicial functions
Tal et al. Unsteady aeroservoelastic modeling of flexible wing generic transport aircraft with variable camber continuous trailing edge flap
Hillebrand et al. Active Gust Alleviation on a High Aspect Ratio Wing Based on High Fidelity CFD Simulations
Fleischmann et al. Analysis of wake surfing benefits using a fast unsteady vortex lattice method
Bragin et al. Improvement of aerodynamics of civil plane wing high-lift devices
Bhargava et al. Analysis of multi-element airfoil configurations: a numerical approach
Miao et al. The Aerodynamic Characteristics of a Diamond Joined-Wing Morphing Aircraft
Bosman Refinement of glider aerodynamic design using CFD
Samal et al. Reduction of wingtip vortex from suction at wingtip
Demirer Static and dynamic aeroelastic analysis of a very light aircraft
Khalid Crosswise wind shear represented as a ramped velocity profile impacting a forward-moving aircraft
Naveen Aerodynamic Analysis of C-Wing Aircraft
Jirasek et al. Computational design of S-duct intakes for the NATO AVT-251 multi-disciplinary configuration
Eslami et al. Aerodynamic performance of Parastoo UAV

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