CN111400824A - 一种确定飞翼布局飞机机翼变形限制条件的方法 - Google Patents
一种确定飞翼布局飞机机翼变形限制条件的方法 Download PDFInfo
- Publication number
- CN111400824A CN111400824A CN202010242925.2A CN202010242925A CN111400824A CN 111400824 A CN111400824 A CN 111400824A CN 202010242925 A CN202010242925 A CN 202010242925A CN 111400824 A CN111400824 A CN 111400824A
- Authority
- CN
- China
- Prior art keywords
- deformation
- wing
- aerodynamic
- solving
- configuration
- 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.)
- Withdrawn
Links
Images
Classifications
-
- 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
-
- 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)
- Manufacturing & Machinery (AREA)
- Transportation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提出一种确定飞翼布局飞机机翼变形限制条件的方法,步骤为,刚性构型气动数据求解,不考虑飞机的弹性变形,求解刚性飞机的气动数据;弹性构型气动数据求解,通过假设一个机翼变形量,重构飞机模型,即弹性构型,进行气动数据求解;变形水平与气动特性变化对应关系求解,通过假设不同的变形水平,获取相应的气动特性变化量,形成对应关系;机翼变形限制条件求解,通过变形水平与气动特性变化的对应关系,结合气动特性允许变化量,求解出机翼变形允许量。本发明在设计之初就提出变形限制条件,保证气动设计不偏离目标值,同时为结构设计提供参考;减少了气动‑结构设计的迭代次数,缩短飞机设计周期;使飞翼布局飞机方案设计高效、准确。
Description
技术领域
本发明涉及飞行器设计技术领域,特别是飞翼布局飞机的机翼设计技术领域。
背景技术
飞翼布局飞机是指无尾翼,机身与机翼高度融合设计的一种飞机。此类布局飞机的稳定力矩与操纵力矩主要由机翼以及布置在机翼上的舵面提供,而机翼结构刚度相对较小,在空中飞行时,受气动力的作用,极易产生变形,对飞机的气动特性产生影响。机翼通常发生扭转与弯曲变形,导致流场形态发生改变,机翼载荷重新分布,因此全机的升阻特性、力矩特性等将发生变化。若变形过大,可能导致飞机纵向安定性从稳定变为不稳定,危及飞行安全。因此,在飞机方案设计阶段需采用相应的方法,对弹性变形的气动影响量进行评估,进而给出变形限制条件,为结构强度设计提供依据,保证飞机的飞行安全和性能指标满足要求。
弹性变形限制设计属于结构、气动力、飞行力学等多学科优化设计范畴,目前采用的设计方法大致可分为工程推算与多学科优化方法两种。工程推算指:依据部件气动力,通过改变部件局部气动迎角来等效弹性变形,获得弹性变形与气动变化的对应关系,从而推导出弹性变形要求;多学科优化方法是指:采用气动力学仿真与结构有限元分析相耦合的静气动弹性计算技术,通过优化计算的方式,对结构模型进行优化调整,最终获得符合要求的结果。工程估算属粗放型设计方法,误差较大,在现代飞行器设计中已逐步被淘汰;多学科优化方法符合现代复杂飞行器最优化设计潮流,但要求的计算资源过大,计算效率低,不适用于前期快速设计。因此提出一种简化的机翼弹性变形限制设计方法,满足在方案阶段进行高效、准确设计的要求。
发明内容
本发明的其目的是:在飞机设计时,通过限制机翼变形量来限制全机气动特性变化,保证气动设计满足设计要求,使飞翼布局飞机方案设计阶段进行高效、准确设计。
本方法采用的技术方案是:
一种确定飞翼布局飞机机翼变形限制条件的方法,通过刚性构型气动数据求解、弹性构型气动数据求解、变形水平与气动特性变化对应关系求解、机翼变形限制条件求解,求解出机翼变形允许量,即为机翼变形限制;
所述的刚性构型气动数据求解为:不考虑飞机的弹性变形,求解刚性飞机的气动数据;
所述的弹性构型气动数据求解为:通过假设一个机翼变形量,重构飞机模型,即为弹性构型,并对之进行气动数据求解;
所述的变形水平与气动特性变化对应关系求解为:通过假设不同的变形水平,获取相应的气动特性变化量,形成对应关系;
所述的机翼变形限制条件求解,通过变形水平与气动特性变化的对应关系,结合气动特性允许变化量,求解出机翼变形允许量。
进一步地,所述的刚性构型气动数据求解为:
刚性构型气动数据获取:不考虑飞机的弹性变形,利用CFD进行刚性飞机的气动计算,获取一系列攻角[α1、α2、……αn]下的气动特性,主要包括升力系数[CL1、CL2、……CLn]、俯仰力矩系数[Cm1、Cm2、……Cmn]、升力线斜率CLα0、静稳定导数CmCL0,并提取出相应的机翼载荷,记为[F1、F2、……Fn];
获取刚性飞机在机翼载荷为零时的攻角和对应的气动数据:根据上一步求得的[α1、α2、……αn]和[F1、F2、……Fn]的对应关系,采用数值插值,求取机翼载荷F=0时对应的攻角,记为α00,对应的CL、Cm值记为CL00和Cm00。
进一步地,所述的弹性构型气动数据求解为:
构造弹性构型的气动外形,并求取其气动数据:所述的弹性构型,是指假设一个机翼翼尖剖面的变形量T1,其可以是弯曲变形,也可以是扭转变形,也可以是二者的组合,翼根剖面的变形量近似为0,利用翼根和翼尖剖面的变形量,插值获取各个展向站位的翼剖面的变形量,并根据各剖面的变形量重构出变形后的气动外形,此时飞机构型记为“弹性构型1”;
利用CFD求出该构型在攻角序列[α1、α2、……αn]下的气动数据,记为[CL11、CL12、……CL1n]和[Cm11、Cm12、……Cm1n],同时提取出相应的机翼载荷,记为[F11、F12、……F1n];
假设变形水平1,并获取该变形水平下的气动特性变化量
所述的“变形水平1”,是假设机翼在攻角αk时的载荷Fk作用下,产生的变形为变形T1,此时飞机变形为“弹性构型1”,其中αk∈[α1、α2、……αn],此时对应的升力系数、力矩系数分别为CLk、Cmk;
飞机在不同攻角时,机翼载荷不一致,产生的变形亦随之变化,在变形水平1的假设下,气动特性为:
升力线斜率CLα1=(CLk-CL00)/(αk-α00)
静稳定导数CmCL1=(Cmk-Cm00)/(CLk-CL00)
升力线斜率变化量ΔCLα1=CLα1-CLα0
静稳定导数变化量ΔCmCL1=CmCL1-CmCL0。
进一步地,所述的变形水平与气动特性变化对应关系求解为:通过改变αk的值来改变变形水平,获取不同变形水平对应的气动特性变化量,得到变形水平与气动特性变化量的对应关系矩阵:
进一步地,所述的机翼变形限制条件求解为:所述的机翼变形限制条件求解,根据允许的气动特性变化量,结合变形水平与气动特性变化量的的对应关系矩阵,插值求得允许气动特性变化量对应的变形水平,即为机翼变形限制条件。
本发明的有益效果主要表现在:在飞机设计之初就可以提出飞翼布局飞机的机翼变形限制条件,保证气动设计不偏离设计目标值,同时为结构设计提供参考;通过本发明方法,减少了气动-结构设计动作的迭代次数,缩短飞机设计周期;使飞翼布局飞机方案设计高效、准确。
附图说明
图1是本方法的流程图。
具体实施方式
下面结合附图和实施步骤对本发明进一步说明。
实施例1
本发明一种确定飞翼布局飞机机翼变形限制条件的方法,通过刚性构型气动数据求解、弹性构型气动数据求解、变形水平与气动特性变化对应关系求解、机翼变形限制条件求解,求解出机翼变形允许量,即为机翼变形限制;
所述的刚性构型气动数据求解为:不考虑飞机的弹性变形,求解刚性飞机的气动数据;
所述的弹性构型气动数据求解为:通过假设一个机翼变形量,重构飞机模型,即为弹性构型,并对之进行气动数据求解;
所述的变形水平与气动特性变化对应关系求解为:通过假设不同的变形水平,获取相应的气动特性变化量,形成对应关系;
所述的机翼变形限制条件求解,通过变形水平与气动特性变化的对应关系,结合气动特性允许变化量,求解出机翼变形允许量。
实施例2
一种确定机翼变形限制条件的方法,包括如下步骤:
1、刚性构型气动数据获取
首先,不考虑飞机的弹性变形,利用CFD(计算流体动力学)进行刚性飞机的气动计算,获取一系列攻角[α1、α2、……αn]下的气动特性,主要包括升力系数[CL1、CL2、……CLn]、俯仰力矩系数[Cm1、Cm2、……Cmn]、升力线斜率CLα0、静稳定导数CmCL0,并提取出相应的机翼载荷,记为[F1、F2、……Fn]。
2.获取刚性飞机在机翼载荷为零时的攻角和对应的气动数据
根据第1步求得的[α1、α2、……αn]和[F1、F2、……Fn]的对应关系,采用数值插值,求取机翼载荷F=0时对应的攻角,记为α00,对应的CL、Cm值记为CL00和Cm00。
3.构造弹性构型的气动外形,并求取其气动数据
所述的弹性构型,是指:假设一个机翼翼尖剖面的变形量T1,其可以是弯曲变形,也可以是扭转变形,也可以是二者的组合,翼根剖面的变形量近似为0,利用翼根和翼尖剖面的变形量,插值获取各个展向站位的翼剖面的变形量,并根据各剖面的变形量重构出变形后的气动外形,此时飞机构型记为“弹性构型1”。
利用CFD求出该构型在攻角序列[α1、α2、……αn]下的气动数据,记为[CL11、CL12、……CL1n]和[Cm11、Cm12、……Cm1n],同时提取出相应的机翼载荷,记为[F11、F12、……F1n]。
4.假设变形水平1,并获取该变形水平下的气动特性变化量
所述的“变形水平1”,是假设机翼在攻角αk时的载荷Fk作用下,产生的变形为变形T1,此时飞机变形为“弹性构型1”,其中αk∈[α1、α2、……αn],此时对应的升力系数、力矩系数分别为CLk、Cmk。
飞机在不同攻角时,机翼载荷不一致,产生的变形亦随之变化,在变形水平1的假设下,气动特性为:
攻角 | 升力系数 | 俯仰力矩系数 | 翼尖变形 |
α<sub>00</sub> | CL<sub>00</sub> | Cm<sub>00</sub> | 0 |
α<sub>k</sub> | CL<sub>k</sub> | Cm<sub>k</sub> | 变形量T1 |
升力线斜率CLα1=(CLk-CL00)/(αk-α00)
静稳定导数CmCL1=(Cmk-Cm00)/(CLk-CL00)
升力线斜率变化量ΔCLα1=CLα1-CLα0
静稳定导数变化量ΔCmCL1=CmCL1-CmCL0
5.通过改变αk的值来改变变形水平,获取不同变形水平对应的气动特性变化量,得到变形水平与气动特性变化量的对应关系矩阵:
6.根据允许的气动特性变化量,结合变形水平与气动特性变化量的的对应关系矩阵,插值求得允许气动特性变化量对应的变形水平,即为机翼变形限制条件。
允许的气动特性变化量,通常是飞机设计时,由飞机飞行性能和飞行品质要求决定的允许的升力线斜率变化量和静稳定导数变化量。
以上所述是本发明的优选实施方式,当然不能以此来限定本发明之权利范围,应当指出,对于本技术领域的普通技术人员来说,还可以做出若干改进和变动,这些改进和变动也视为本发明的保护范围。
Claims (5)
1.一种确定飞翼布局飞机机翼变形限制条件的方法,其特征在于,通过刚性构型气动数据求解、弹性构型气动数据求解、变形水平与气动特性变化对应关系求解、机翼变形限制条件求解,求解出机翼变形允许量,即为机翼变形限制;
所述的刚性构型气动数据求解为:不考虑飞机的弹性变形,求解刚性飞机的气动数据;
所述的弹性构型气动数据求解为:通过假设一个机翼变形量,重构飞机模型,即为弹性构型,并对之进行气动数据求解;
所述的变形水平与气动特性变化对应关系求解为:通过假设不同的变形水平,获取相应的气动特性变化量,形成对应关系;
所述的机翼变形限制条件求解,通过变形水平与气动特性变化的对应关系,结合气动特性允许变化量,求解出机翼变形允许量。
2.根据权利要求1所述一种确定飞翼布局飞机机翼变形限制条件的方法,其特征在于,所述的刚性构型气动数据求解为:
刚性构型气动数据获取:不考虑飞机的弹性变形,利用CFD进行刚性飞机的气动计算,获取一系列攻角[α1、α2、……αn]下的气动特性,主要包括升力系数[CL1、CL2、……CLn]、俯仰力矩系数[Cm1、Cm2、……Cmn]、升力线斜率CLα0、静稳定导数CmCL0,并提取出相应的机翼载荷,记为[F1、F2、……Fn];
获取刚性飞机在机翼载荷为零时的攻角和对应的气动数据:根据上一步求得的[α1、α2、……αn]和[F1、F2、……Fn]的对应关系,采用数值插值,求取机翼载荷F=0时对应的攻角,记为α00,对应的CL、Cm值记为CL00和Cm00。
3.根据权利要求1所述一种确定飞翼布局飞机机翼变形限制条件的方法,其特征在于,所述的弹性构型气动数据求解为:
构造弹性构型的气动外形,并求取其气动数据:所述的弹性构型,是指假设一个机翼翼尖剖面的变形量T1,其可以是弯曲变形,也可以是扭转变形,也可以是二者的组合,翼根剖面的变形量近似为0,利用翼根和翼尖剖面的变形量,插值获取各个展向站位的翼剖面的变形量,并根据各剖面的变形量重构出变形后的气动外形,此时飞机构型记为“弹性构型1”;
利用CFD求出该构型在攻角序列[α1、α2、……αn]下的气动数据,记为[CL11、CL12、……CL1n]和[Cm11、Cm12、……Cm1n],同时提取出相应的机翼载荷,记为[F11、F12、……F1n];
假设变形水平1,并获取该变形水平下的气动特性变化量
所述的“变形水平1”,是假设机翼在攻角αk时的载荷Fk作用下,产生的变形为变形T1,此时飞机变形为“弹性构型1”,其中αk∈[α1、α2、……αn],此时对应的升力系数、力矩系数分别为CLk、Cmk;
飞机在不同攻角时,机翼载荷不一致,产生的变形亦随之变化,在变形水平1的假设下,气动特性为:
升力线斜率CLα1=(CLk-CL00)/(αk-α00)
静稳定导数CmCL1=(Cmk-Cm00)/(CLk-CL00)
升力线斜率变化量ΔCLα1=CLα1-CLα0
静稳定导数变化量ΔCmCL1=CmCL1-CmCL0
5.根据权利要求1所述一种确定飞翼布局飞机机翼变形限制条件的方法,其特征在于,所述的机翼变形限制条件求解为:所述的机翼变形限制条件求解,根据允许的气动特性变化量,结合变形水平与气动特性变化量的的对应关系矩阵,插值求得允许气动特性变化量对应的变形水平,即为机翼变形限制条件。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010242925.2A CN111400824A (zh) | 2020-03-31 | 2020-03-31 | 一种确定飞翼布局飞机机翼变形限制条件的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010242925.2A CN111400824A (zh) | 2020-03-31 | 2020-03-31 | 一种确定飞翼布局飞机机翼变形限制条件的方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111400824A true CN111400824A (zh) | 2020-07-10 |
Family
ID=71436766
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010242925.2A Withdrawn CN111400824A (zh) | 2020-03-31 | 2020-03-31 | 一种确定飞翼布局飞机机翼变形限制条件的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111400824A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111959817A (zh) * | 2020-07-29 | 2020-11-20 | 成都飞机工业(集团)有限责任公司 | 一种确定飞翼布局飞机机翼变形限制条件的方法 |
CN112035947A (zh) * | 2020-07-29 | 2020-12-04 | 成都飞机工业(集团)有限责任公司 | 一种带整体油箱的机翼剖面载荷计算方法 |
CN112224438A (zh) * | 2020-09-09 | 2021-01-15 | 西北工业大学 | 飞行器局部结构的保形设计方法及飞行器 |
CN115183731A (zh) * | 2022-06-20 | 2022-10-14 | 成都飞机工业(集团)有限责任公司 | 一种机翼表面监测系统及其传感器布置方法 |
-
2020
- 2020-03-31 CN CN202010242925.2A patent/CN111400824A/zh not_active Withdrawn
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111959817A (zh) * | 2020-07-29 | 2020-11-20 | 成都飞机工业(集团)有限责任公司 | 一种确定飞翼布局飞机机翼变形限制条件的方法 |
CN112035947A (zh) * | 2020-07-29 | 2020-12-04 | 成都飞机工业(集团)有限责任公司 | 一种带整体油箱的机翼剖面载荷计算方法 |
CN112035947B (zh) * | 2020-07-29 | 2022-07-15 | 成都飞机工业(集团)有限责任公司 | 一种带整体油箱的机翼剖面载荷计算方法 |
CN112224438A (zh) * | 2020-09-09 | 2021-01-15 | 西北工业大学 | 飞行器局部结构的保形设计方法及飞行器 |
CN112224438B (zh) * | 2020-09-09 | 2024-06-07 | 西北工业大学 | 飞行器局部结构的保形设计方法及飞行器 |
CN115183731A (zh) * | 2022-06-20 | 2022-10-14 | 成都飞机工业(集团)有限责任公司 | 一种机翼表面监测系统及其传感器布置方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111400824A (zh) | 一种确定飞翼布局飞机机翼变形限制条件的方法 | |
Bae et al. | Aerodynamic and static aeroelastic characteristics of a variable-span morphing wing | |
CN105183996A (zh) | 面元修正与网格预先自适应计算方法 | |
CN113886967B (zh) | 多巡航工况的大型飞机机翼气动弹性优化方法 | |
CN104317980A (zh) | 一种共轴刚性旋翼气动弹性响应方法 | |
CN111274648B (zh) | 一种民用飞机前缘襟翼的分布式飞行载荷设计方法 | |
Brooks et al. | Undeflected common research model (uCRM): an aerostructural model for the study of high aspect ratio transport aircraft wings | |
CN109540459B (zh) | 一种气动特性数值计算结果修正方法 | |
Öhrle et al. | Compound helicopter X 3 in high-speed flight: correlation of simulation and flight test | |
Jasa et al. | Design and trajectory optimization of a morphing wing aircraft | |
CN111017248B (zh) | 一种飞机机翼的静气动弹性修正方法 | |
CN110851913B (zh) | 一种直升机气动噪声确定方法 | |
CN111959817A (zh) | 一种确定飞翼布局飞机机翼变形限制条件的方法 | |
Chen et al. | Overset Euler/Boundary-Layer solver with panel-based aerodynamic modeling for aeroelastic applications | |
CN105117541B (zh) | 一种正向型架外形优化设计方法 | |
CN115358167B (zh) | 一种考虑发动机参数的飞发一体气动伴随优化设计方法 | |
Sitaraman et al. | Coupled CFD/CSD prediction of rotor aerodynamic and structural dynamic loads for three critical flight conditions | |
CN114001913B (zh) | 一种基于风洞试验的部件气动载荷修正方法 | |
McWaters | F-35 conventional mode jet-effects testing methodology | |
CN115906685A (zh) | 一种跨声速非线性气动力修正方法 | |
Nguyen et al. | Aeroelastic analysis of wind tunnel test data of a flexible wing with a variable camber continuous trailing edge flap (VCCTEF) | |
Ahmad et al. | High-order accurate CFD/CSD simulation of the UH60 rotor in forward flight | |
Tian et al. | Aerodynamic/mechanism optimization of a variable camber Fowler flap for general aviation aircraft | |
Klausmeyer | Stability derivative estimation: Methods and practical considerations for conventional transonic aircraft | |
Mancini et al. | The Effect of Maneuver Load Alleviation Strategies on Aircraft Performance Indicators |
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 | ||
WW01 | Invention patent application withdrawn after publication | ||
WW01 | Invention patent application withdrawn after publication |
Application publication date: 20200710 |