CN109299579B - 大展弦比飞机风洞测力试验数据修正方法 - Google Patents
大展弦比飞机风洞测力试验数据修正方法 Download PDFInfo
- Publication number
- CN109299579B CN109299579B CN201811408277.2A CN201811408277A CN109299579B CN 109299579 B CN109299579 B CN 109299579B CN 201811408277 A CN201811408277 A CN 201811408277A CN 109299579 B CN109299579 B CN 109299579B
- Authority
- CN
- China
- Prior art keywords
- model
- wind tunnel
- pressure distribution
- finite element
- whole
- 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
Images
Classifications
-
- 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]
-
- 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
-
- 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)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本申请提供了一种大展弦比飞机风洞测力试验数据修正方法,包括:基于风洞测力试验模型,建立CFD模型、偶极子模型以及结构有限元模型;根据所述CFD模型,获得风洞试验条件下的全机压力分布;根据所述偶极子模型和所述结构有限元模型,建立气动/结构耦合计算模型;将所述全机压力分布加载到所述气动/结构耦合计算模型上,获得全机弹性变形和气动力修正量;根据所述气动力修正量对风洞测力试验数据进行修正。
Description
技术领域
本申请涉及全机载荷设计和全机静气弹性预测技术领域,具体提供一种大展弦比飞机风洞测力试验数据修正方法。
背景技术
随着民用大展弦比飞机的发展,先进民用飞机的研制均需要在高雷诺数增压风洞中进行大量的气动试验。此类飞行器风洞试验过程中由于模型材料、加工工艺等原因,结构不可避免地发生弹性变形,进一步影响到模型的气动力特性。尤其对于大展弦比飞机来说,模型变形十分明显。因此在大展弦比飞机增压风洞模型气动试验中,去除弹性变形对试验数据带来的影响,将弹性变形后气动力数据修正成变形前的刚体气动力数据是十分必要的。
目前国内的测力风洞试验中忽略风洞试验模型弹性变形的影响,一般不进行气动力数据修正。使用变形后气动数据进行飞机设计,从数据源头上就存在误差。
发明内容
为了解决上述技术问题至少之一,本申请提供了一种大展弦比飞机风洞测力试验数据修正方法,包括:基于风洞测力试验模型,建立CFD模型、偶极子模型以及结构有限元模型;根据所述CFD模型,获得风洞试验条件下的全机压力分布;根据所述偶极子模型和所述结构有限元模型,建立气动/结构耦合计算模型;将所述全机压力分布加载到所述气动/结构耦合计算模型上,获得全机弹性变形和气动力修正量;根据所述气动力修正量对风洞测力试验数据进行修正。
根据本申请的至少一个实施例,所述风洞测力试验模型为CATIA数模。
根据本申请的至少一个实施例,所述结构有限元模型为实体模型,包括天平、机身、机翼以及连接部件。
根据本申请的至少一个实施例,根据所述CFD模型,获得风洞试验条件下的全机压力分布,包括:针对所述CFD模型进行定常计算,计算得到风洞试验条件下的全机压力分布。
本申请实施例提供的大展弦比飞机风洞测力试验数据修正方法中,通过采用建模计算与风洞试验相结合,计算得到弹性变形和修正量的方法,比单纯增加风洞变形测量试验以及测压试验大大降低成本费用,通过偶极子模型与结构有限元模型耦合来建立气动/结构耦合模型,引入CFD数据计算得到弹性气动力和模型变形,不仅保证原始外形下气动力计算准确性,而且计算方法高效快捷,适合工程型号应用。
附图说明
图1是本申请实施例提供的大展弦比飞机风洞测力试验数据修正方法的流程示意图;
图2是本申请实施例提供的风洞试验模型CATIA数模;
图3是本申请实施例提供的结构有限元模型;
图4是本申请实施例提供的全机压力分布示意图;
图5是本申请实施例提供的偶极子模型;
图6是本申请实施例提供的机翼弯曲变形对比图;
图7是本申请实施例提供的机翼扭转变形对比图。
具体实施方式
下面结合附图和实施例对本申请作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释相关申请,而非对该申请的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本申请相关的部分。
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
图1是本申请实施例提供的大展弦比飞机风洞测力试验数据修正方法的流程示意图。如图1所示,该方法包括以下步骤:
步骤101,基于风洞测力试验模型,建立CFD模型、偶极子模型以及结构有限元模型。
在本实施例中,风洞测力试验模型为CATIA数模,如图2所示,结构有限元模型为实体模型,如图3所示,其包括天平、机身、机翼以及连接部件,有限元建模的关键在于保证模型整体刚度水平,建立的结构有限元模型刚度与风洞试验模型的刚度一致。
步骤102,根据CFD模型,获得风洞试验条件下的全机压力分布。
针对CFD模型进行定常计算,可以计算得到风洞试验条件下的全机压力分布,如图4所示。
步骤103,根据偶极子模型和结构有限元模型,建立气动/结构耦合计算模型。
偶极子模型如图5所示,由于偶极子网格非常容易与结构有限元模型耦合,因此,将偶极子网格和结构有限元模型耦合来得到气动/结构耦合计算模型的计算效率较高。
步骤104,将全机压力分布加载到所述气动/结构耦合计算模型上,获得全机弹性变形和气动力修正量。
全机弹性变形包括机翼弯曲变形和机翼扭转变形,其中,机翼弯曲变形的对比图如图6所示,机翼扭转变形的对比图如图7所示。
步骤105,根据气动力修正量对风洞测力试验数据进行修正。
在本实施例中,通过采用建模计算与风洞试验相结合,计算得到弹性变形和修正量的方法,比单纯增加风洞变形测量试验以及测压试验大大降低成本费用。通过偶极子网格与结构有限元模型耦合来建立气动/结构耦合模型,引入CFD数据计算得到弹性气动力和模型变形,不仅保证原始外形下气动力计算准确性,而且计算方法高效快捷,适合工程型号应用。
至此,已经结合附图所示的优选实施方式描述了本申请的技术方案,但是,本领域技术人员容易理解的是,本申请的保护范围显然不局限于这些具体实施方式。在不偏离本申请的原理的前提下,本领域技术人员可以对相关技术特征作出等同的更改或替换,这些更改或替换之后的技术方案都将落入本申请的保护范围之内。
Claims (1)
1.一种大展弦比飞机风洞测力试验数据修正方法,其特征在于,包括:
基于风洞测力试验模型,建立CFD模型、偶极子模型以及结构有限元模型;
根据所述CFD模型,获得风洞试验条件下的全机压力分布;
根据所述偶极子模型和所述结构有限元模型,建立气动/结构耦合计算模型;
将所述全机压力分布加载到所述气动/结构耦合计算模型上,获得全机弹性变形和气动力修正量;
根据所述气动力修正量对风洞测力试验数据进行修正;
所述风洞测力试验模型为CATIA数模;
所述结构有限元模型为实体模型,包括天平、机身、机翼以及连接部件;
根据所述CFD模型,获得风洞试验条件下的全机压力分布,包括:
针对所述CFD模型进行定常计算,计算得到风洞试验条件下的全机压力分布。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811408277.2A CN109299579B (zh) | 2018-11-23 | 2018-11-23 | 大展弦比飞机风洞测力试验数据修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811408277.2A CN109299579B (zh) | 2018-11-23 | 2018-11-23 | 大展弦比飞机风洞测力试验数据修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109299579A CN109299579A (zh) | 2019-02-01 |
CN109299579B true CN109299579B (zh) | 2023-05-23 |
Family
ID=65143965
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811408277.2A Active CN109299579B (zh) | 2018-11-23 | 2018-11-23 | 大展弦比飞机风洞测力试验数据修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109299579B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110160737B (zh) * | 2019-03-21 | 2020-12-25 | 北京机电工程研究所 | 基于工程面元法的弹性气动数据精确获取方法 |
CN110155363B (zh) * | 2019-03-21 | 2020-12-08 | 北京机电工程研究所 | 基于cfd方法的弹性气动数据精确获取方法 |
CN112414668B (zh) * | 2020-11-03 | 2022-08-09 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 一种风洞试验数据静气弹修正方法、装置、设备及介质 |
CN113086243B (zh) * | 2021-04-20 | 2022-08-02 | 中国直升机设计研究所 | 一种直升机全机质量体惯性载荷的分配方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103366052B (zh) * | 2013-06-27 | 2016-06-01 | 中国航天空气动力技术研究院 | 一种高超声速飞行器热静气弹分析方法 |
CN103412994B (zh) * | 2013-08-08 | 2016-03-16 | 空气动力学国家重点实验室 | 一种确定高速风洞大型飞机模型缩比的方法 |
CN104133926B (zh) * | 2014-04-23 | 2017-06-16 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种弹性气动力特性综合分析方法 |
CN105183996B (zh) * | 2015-09-14 | 2018-07-27 | 西北工业大学 | 面元修正与网格预先自适应计算方法 |
CN106650029B (zh) * | 2016-11-28 | 2019-11-22 | 博动医学影像科技(上海)有限公司 | 基于cfd仿真的分叉血管压力差及ffr的快速计算方法及系统 |
CN106840572B (zh) * | 2016-12-19 | 2019-05-24 | 中国航天空气动力技术研究院 | 一种临近空间大展弦比柔性飞行器风洞试验数据修正方法 |
CN107391858B (zh) * | 2017-07-27 | 2020-06-30 | 空气动力学国家重点实验室 | 一种获取风洞模型静气动弹性变形影响量的方法 |
CN107391891B (zh) * | 2017-09-05 | 2020-07-07 | 北京理工大学 | 一种基于模型融合方法的大展弦比机翼优化设计方法 |
-
2018
- 2018-11-23 CN CN201811408277.2A patent/CN109299579B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109299579A (zh) | 2019-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109299579B (zh) | 大展弦比飞机风洞测力试验数据修正方法 | |
RU2733001C2 (ru) | Способ изготовления прокладки и сборочный узел летательного аппарата, содержащий такую прокладку | |
Chetan et al. | Multi‐fidelity digital twin structural model for a sub‐scale downwind wind turbine rotor blade | |
Cavagna et al. | Application of Navier–Stokes simulations for aeroelastic stability assessment in transonic regime | |
CN105183996A (zh) | 面元修正与网格预先自适应计算方法 | |
Luczak et al. | Updating finite element model of a wind turbine blade section using experimental modal analysis results | |
CN101625351A (zh) | 一种蠕变数据转换为材料高温应力松弛数据的方法 | |
Peng et al. | Numerical investigation of the effects of structural geometric and material nonlinearities on limit-cycle oscillation of a cropped delta wing | |
CN105300595A (zh) | 连接模型校准的天平加载头及加载方法 | |
CN108693054B (zh) | 一种复合材料结构谱载疲劳寿命的渐进损伤算法 | |
CN105716842B (zh) | 一种双梁式长直机翼载荷处理方法 | |
CN107766612B (zh) | 一种连接翼结构形式机翼载荷测量方法 | |
JP2020504859A (ja) | しわを抑制する空間薄膜構造治具形状を最適化する設計方法 | |
Wang et al. | Sensitivity analysis of multidisciplinary rotorcraft simulations | |
CN110826208A (zh) | 一种气动数值模拟加速收敛方法 | |
CN111274648A (zh) | 一种民用飞机前缘襟翼的分布式飞行载荷设计方法 | |
CN104834795B (zh) | 包带连接结构接触摩擦非线性特性模拟方法及系统 | |
JP2014074587A (ja) | 多分力計測方法 | |
CN109323841B (zh) | 基于网格的机翼总载荷与分布载荷的协调方法 | |
Neto et al. | Control-point-placement method for the aerodynamic correction of the vortex-and the doublet-lattice methods | |
CN108804755B (zh) | 通过测量来校准包括主结构和至少一个支撑结构的模型的支撑结构的刚度的方法 | |
CN115618501A (zh) | 基于数据融合修正的亚跨气动特性获取方法、系统及装置 | |
CN113094953B (zh) | 一种带机翼变形的铰链力矩天平的有限元分析方法 | |
Wellmer et al. | A comparison of fluid/structure coupling methods for reduced structural models | |
CN111959817A (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 |