CN106672260A - 一种大展弦比机翼型架构型设计方法 - Google Patents
一种大展弦比机翼型架构型设计方法 Download PDFInfo
- Publication number
- CN106672260A CN106672260A CN201510763377.7A CN201510763377A CN106672260A CN 106672260 A CN106672260 A CN 106672260A CN 201510763377 A CN201510763377 A CN 201510763377A CN 106672260 A CN106672260 A CN 106672260A
- Authority
- CN
- China
- Prior art keywords
- wing
- optimal solution
- lift
- elastic
- initial population
- 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.)
- Granted
Links
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
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Transportation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Structures Of Non-Positive Displacement Pumps (AREA)
Abstract
本发明公开了一种大展弦比机翼型架构型设计方法,涉及到飞行器技术领域,包括刚性机翼优化步骤和弹性机翼优化步骤,生成机翼构型初始种群,进行升阻比计算,获取局部最优解,再进行寻优计算,然后对局部最优解进行升阻比计算并与初始种群进行升阻比计算获得的局部最优解对比,最后判断全局最优解;计算弹性机翼飞行时所受气动载荷,再计算机翼的弹性变形,通过多次迭代计算,求解在该飞行状态下机翼结构弹性平衡时的真实形态,获取弹性变形后的飞行载荷,得到翼型剖面型架的扭转角,完成机翼型架构型的设计。本发明设计的机翼具有良好的几何气动扭转角分布,升力分布更佳,能够产生更小的诱导阻力,获取最佳的升阻比,从而提高巡航效率。
Description
技术领域
本发明涉及到飞行器技术领域,特别涉及一种大展弦比机翼型架构型设计方法。
背景技术
在民机和长航时无人机设计中,提高巡航升阻比、提高巡航效率是设计的重点,而机翼型架设计,是获取最优巡航构型,提高巡航效率的有效手段之一。
一般具有椭圆升力分布的机翼拥有均匀的下洗分布,具有更小的诱导阻力,对于任意一大展弦比机翼来说,总存在一种优化的几何气动扭转角分布,使得升力分布最佳,产生更小的诱导阻力,获取最佳的升阻比。同时大展弦比机翼由于结构轻、刚度小,受到气动载荷时,也容易产生较大的扭转和弯曲变形,对飞行性能和飞行品质产生影响,严重时可能出现扭转发散等危及飞机安全的极端情况,因此也必须以最佳巡航构型为设计目标,对机翼进行计及弹性变形影响的型架外形优化设计,才能消除弹性变形对气动特性的不利影响,维持巡航最佳型架外形。
公开号为CN 103407580A,公开日为2013年11月27日的中国专利文献公开了一种超大展弦比轻质高升阻比机翼的制造方法,其特征在于,包括如下步骤:一、根据飞行器总体设计要求,确定机翼的展长和展弦比;二、根据飞行器总体设计要求,确定机翼弹性变形的上限值和扭转变形的上限值;三、根据飞行器总体设计要求,计算飞行器在不同飞行条件下的气动载荷,以及飞行器在飞行过程中,由于外力作用引起的气动载荷;四、计算机翼本身的质量力及其他部件传来的力;五、在飞行器的设计飞行条件范围内,计算机翼不安装张线时可能发生的极限变形量;六、在飞行器上确定张线支架;七、在飞行器的张线支架两边对称设置张线安装点;八、在张线支架与张线安装点之间设置张线;九、根据步骤三计算的飞行器在不同飞行条件下的气动载荷、由于外力作用引起的气动载荷以及步骤八确定的张线布局,计算机翼变形在设计允许的弹性变形范围内,张线所受的最大拉力;然后根据张线的材料,计算出张线的直径;十、根据步骤三计算的飞行器在不同飞行条件下的气动载荷、由于外力作用引起的气动载荷以及步骤八确定的张线布局,计算各张线的受力情况、机翼的载荷分布以及机翼的最大变形量,判断机翼是否符合极限弹性变形的设计要求以及张线对机翼弹性变形的影响;同时计算张线本身及张线支点承受的应力,判断张线本身及张线支点承受的应力是否在材料允许范围之内;根据判断的结果,对张线的布局方式进行修正;十一、计算修正后的安装张线的飞行器的气动性能,判断张线对飞行器气动性能的影响,完成机翼的制造。
该专利文献公开的超大展弦比轻质高升阻比机翼的制造方法,通过计算各张线的受力情况、机翼的载荷分布以及机翼的最大变形量,从而判断机翼是否符合极限弹性变形的设计要求以及张线对机翼弹性变形的影响,这种方式不能准确测出机翼的几何气动扭转角分布,难以获得最佳升阻比,导致巡航效率较低。
发明内容
本发明为了克服上述现有技术的缺陷,提供一种大展弦比机翼型架构型设计方法,本发明设计出的机翼具有良好的几何气动扭转角分布,升力分布更佳,能够产生更小的诱导阻力,获取最佳的升阻比,从而提高巡航效率。
本发明通过下述技术方案实现:
一种大展弦比机翼型架构型设计方法,其特征在于,包括如下步骤:
A、刚性机翼优化步骤,生成机翼构型初始种群,对初始种群进行升阻比计算,获取局部最优解,对获取的局部最优解进行寻优计算,然后对寻优计算后得到的局部最优解进行升阻比计算并与初始种群进行升阻比计算获得的局部最优解对比,选出全局最优解,最后判断全局最优解,直到收敛为止;
B、弹性机翼优化步骤, 计算弹性机翼飞行时所受气动载荷,再耦合结构静平衡方程计算机翼的弹性变形,通过多次迭代计算,求解在该飞行状态下机翼结构弹性平衡时的真实形态,获取弹性变形后的飞行载荷,得到翼型剖面型架的扭转角,完成机翼型架构型的设计。
所述机翼构型初始种群是通过机翼剖面扭转角初始矩阵计算方法产生的。
所述初始种群的升阻比计算采用的是计算流体力学软件。
所述局部最优解是通过类电磁机制原理的无约束优化方法获取的。
本发明所述CFD是指计算流体力学软件。
本发明的有益效果主要表现在:
本发明中,通过刚性机翼优化步骤,“生成机翼构型初始种群,对初始种群进行升阻比计算,获取局部最优解,对获取的局部最优解进行寻优计算,然后对寻优计算后得到的局部最优解进行升阻比计算并与初始种群进行升阻比计算获得的局部最优解对比,选出全局最优解,最后判断全局最优解,直到收敛为止”,得到的大展弦比机翼具有良好的刚度和强度;通过弹性机翼优化步骤, “ 计算弹性机翼飞行时所受气动载荷,再耦合结构静平衡方程计算机翼的弹性变形,通过多次迭代计算,求解在该飞行状态下机翼结构弹性平衡时的真实形态,获取弹性变形后的飞行载荷,得到翼型剖面型架的扭转角”,大幅降低了气动弹性变形量;采用该方法制造的大展弦比机翼具有良好的几何气动扭转角分布,升力分布更佳,能够产生更小的诱导阻力,获取最佳的升阻比,从而提高巡航效率。
附图说明
图1是本发明大展弦比机翼型架构型设计流程框图。
具体实施方式
参见图1,一种大展弦比机翼型架构型设计方法,包括如下步骤:
A、刚性机翼优化步骤,生成机翼构型初始种群,对初始种群进行升阻比计算,获取局部最优解,对获取的局部最优解进行寻优计算,然后对寻优计算后得到的局部最优解进行升阻比计算并与初始种群进行升阻比计算获得的局部最优解对比,选出全局最优解,最后判断全局最优解,直到收敛为止;
B、弹性机翼优化步骤, 计算弹性机翼飞行时所受气动载荷,再耦合结构静平衡方程计算机翼的弹性变形,通过多次迭代计算,求解在该飞行状态下机翼结构弹性平衡时的真实形态,获取弹性变形后的飞行载荷,得到翼型剖面型架的扭转角,完成机翼型架构型的设计。
所述步骤A具体体现在:
1)按照机翼剖面扭转角初始矩阵计算方法产生初始种群。
2)利用CFD对初始种群进行升阻比计算,并采用类电磁机制原理的无约束优化方法获取局部最优解。
3)采用类电磁机制原理的无约束优化方法,对第2步获取的局部最优解进行寻优计算。
4)利用CFD对由第3步获取的最优解进行升阻比计算,并与第2步结果进行对比筛选出全局最优解。
5)判断第4步全局最优解否收敛,如不收敛,将全局最优解作为第3步的输入重复3到5的步骤,直到收敛为止。
所述步骤B具体体现在:采用欧拉方程为流场控制方程,计算弹性机翼飞行时所受气动载荷,再耦合结构静平衡方程计算机翼的弹性变形,通过多次迭代计算,求解在该飞行状态下机翼结构弹性平衡时的真实形态,获取弹性变形后的飞行载荷,在此基础上,对弹性机翼剖面扭角进行优化设计,消除弹性变形对气动载荷分布的影响。
通常情况下,对大展弦比机翼,可以假定机翼的弦向剖面是刚性的,结构弹性引起的机翼扭转角分布对其升力特性的影响是绝对因素,因此可以只对机翼的扭转角进行优化设计,不考虑机翼刚性轴绕度的影响。通过型架优化设计可以使机翼弹性变形平衡后,恢复到机翼的最佳巡航状态,从而提高飞机的巡航升阻比。
优化算子如下式:
对大展弦比机翼,任一翼型剖面,有:
(n=1,2,3,…)
其中:表示第n+1次迭代时的型架设计扭转角分布,=;表示第n次迭代时机翼静弹性变形扭转角分布;表示机翼设计要求的扭转角分布;为迭代松弛参数,。
经过气动力方程和结构静平衡方程的反复耦合迭代,直到:
其中:为判断收敛参数,为一小量;为需要求解的该翼型剖面型架设计扭转角。
机翼剖面扭转角初始矩阵计算方法:每种机翼构型扭转角按照从翼根到翼尖负扭转增大规律排列,计算公式如下:
=
lij表示第j个剖面的第i个扭转角下限;
uij表示第j个剖面的第i个扭转角上限;
βij表示第j个剖面的第i个扭转角;
qj表示第j个剖面初始扭转角个数。
Claims (4)
1.一种大展弦比机翼型架构型设计方法,其特征在于,包括如下步骤:
A、刚性机翼优化步骤,生成初始种群,对初始种群进行升阻比计算,获取局部最优解,对获取的局部最优解进行寻优计算,然后对寻优计算后得到的局部最优解进行升阻比计算并与初始种群进行升阻比计算获得的局部最优解对比,选出全局最优解,最后判断全局最优解,直到收敛为止;
B、弹性机翼优化步骤, 计算弹性机翼飞行时所受气动载荷,再耦合结构静平衡方程计算机翼的弹性变形,通过多次迭代计算,求解在该飞行状态下机翼结构弹性平衡时的真实形态,获取弹性变形后的飞行载荷,得到翼型剖面型架的扭转角,完成机翼型架构型的设计。
2.根据权利要求1所述的一种大展弦比机翼型架构型设计方法,其特征在于:所述初始种群是通过机翼剖面扭转角初始矩阵计算方法产生的。
3.根据权利要求1所述的一种大展弦比机翼型架构型设计方法,其特征在于:所述初始种群的升阻比计算采用的是计算流体力学软件。
4.根据权利要求1所述的一种大展弦比机翼型架构型设计方法,其特征在于:所述局部最优解是通过类电磁机制原理的无约束优化方法获取的。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510763377.7A CN106672260B (zh) | 2015-11-11 | 2015-11-11 | 一种大展弦比机翼型架构型设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510763377.7A CN106672260B (zh) | 2015-11-11 | 2015-11-11 | 一种大展弦比机翼型架构型设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106672260A true CN106672260A (zh) | 2017-05-17 |
CN106672260B CN106672260B (zh) | 2018-11-02 |
Family
ID=58865287
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510763377.7A Active CN106672260B (zh) | 2015-11-11 | 2015-11-11 | 一种大展弦比机翼型架构型设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106672260B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109614747A (zh) * | 2018-12-28 | 2019-04-12 | 西北工业大学 | 大展弦比机翼中管道布局设计方法 |
CN110879123A (zh) * | 2019-12-04 | 2020-03-13 | 安徽信息工程学院 | 一种汽车车身试验扭转刚度的计算方法 |
CN111177853A (zh) * | 2019-12-31 | 2020-05-19 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种机翼型架构型设计方法 |
CN113626934A (zh) * | 2021-08-16 | 2021-11-09 | 成都飞机工业(集团)有限责任公司 | 飞翼布局无人机机翼多构型协调设计方法 |
CN113626935A (zh) * | 2021-08-16 | 2021-11-09 | 成都飞机工业(集团)有限责任公司 | 一种高巡航效率跨音速月形机翼设计方法 |
CN113753256A (zh) * | 2021-09-19 | 2021-12-07 | 中国航空工业集团公司西安飞机设计研究所 | 一种舰载无人预警机顶层参数优化设计方法 |
CN113753257A (zh) * | 2021-09-19 | 2021-12-07 | 中国航空工业集团公司西安飞机设计研究所 | 一种桁架支撑机翼构型飞机展弦比优化方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1997004848A1 (fr) * | 1995-07-26 | 1997-02-13 | Deparis Jerome | Avion jouet pouvant voler par teleguidage dans un espace clos, notamment dans une piece |
CN102069919A (zh) * | 2010-12-23 | 2011-05-25 | 江西洪都航空工业集团有限责任公司 | 一种大展弦比翼型结构及制造方法 |
CN102935897A (zh) * | 2011-08-15 | 2013-02-20 | 中国航空工业集团公司西安飞机设计研究所 | 用于大展弦比机翼高速颤振模型的空心梁及其制作方法 |
CN103407580A (zh) * | 2013-08-19 | 2013-11-27 | 空气动力学国家重点实验室 | 超大展弦比轻质高升阻比机翼的制造方法及其制造的机翼 |
-
2015
- 2015-11-11 CN CN201510763377.7A patent/CN106672260B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1997004848A1 (fr) * | 1995-07-26 | 1997-02-13 | Deparis Jerome | Avion jouet pouvant voler par teleguidage dans un espace clos, notamment dans une piece |
CN102069919A (zh) * | 2010-12-23 | 2011-05-25 | 江西洪都航空工业集团有限责任公司 | 一种大展弦比翼型结构及制造方法 |
CN102935897A (zh) * | 2011-08-15 | 2013-02-20 | 中国航空工业集团公司西安飞机设计研究所 | 用于大展弦比机翼高速颤振模型的空心梁及其制作方法 |
CN103407580A (zh) * | 2013-08-19 | 2013-11-27 | 空气动力学国家重点实验室 | 超大展弦比轻质高升阻比机翼的制造方法及其制造的机翼 |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109614747A (zh) * | 2018-12-28 | 2019-04-12 | 西北工业大学 | 大展弦比机翼中管道布局设计方法 |
CN109614747B (zh) * | 2018-12-28 | 2022-12-09 | 西北工业大学 | 大展弦比机翼中管道布局设计方法 |
CN110879123A (zh) * | 2019-12-04 | 2020-03-13 | 安徽信息工程学院 | 一种汽车车身试验扭转刚度的计算方法 |
CN110879123B (zh) * | 2019-12-04 | 2021-05-18 | 安徽信息工程学院 | 一种汽车车身试验扭转刚度的计算方法 |
CN111177853A (zh) * | 2019-12-31 | 2020-05-19 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种机翼型架构型设计方法 |
CN113626935B (zh) * | 2021-08-16 | 2023-09-29 | 成都飞机工业(集团)有限责任公司 | 一种高巡航效率跨音速月形机翼设计方法 |
CN113626934A (zh) * | 2021-08-16 | 2021-11-09 | 成都飞机工业(集团)有限责任公司 | 飞翼布局无人机机翼多构型协调设计方法 |
CN113626935A (zh) * | 2021-08-16 | 2021-11-09 | 成都飞机工业(集团)有限责任公司 | 一种高巡航效率跨音速月形机翼设计方法 |
CN113626934B (zh) * | 2021-08-16 | 2023-12-08 | 成都飞机工业(集团)有限责任公司 | 飞翼布局无人机机翼多构型协调设计方法 |
CN113753256A (zh) * | 2021-09-19 | 2021-12-07 | 中国航空工业集团公司西安飞机设计研究所 | 一种舰载无人预警机顶层参数优化设计方法 |
CN113753257B (zh) * | 2021-09-19 | 2023-09-22 | 中国航空工业集团公司西安飞机设计研究所 | 一种桁架支撑机翼构型飞机展弦比优化方法 |
CN113753256B (zh) * | 2021-09-19 | 2023-04-18 | 中国航空工业集团公司西安飞机设计研究所 | 一种舰载无人预警机顶层参数优化设计方法 |
CN113753257A (zh) * | 2021-09-19 | 2021-12-07 | 中国航空工业集团公司西安飞机设计研究所 | 一种桁架支撑机翼构型飞机展弦比优化方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106672260B (zh) | 2018-11-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106672260A (zh) | 一种大展弦比机翼型架构型设计方法 | |
CN108121855B (zh) | 基于仿生柔性机翼的小型无人飞行器飞行动力学优化方法 | |
CN105183996A (zh) | 面元修正与网格预先自适应计算方法 | |
CN106840572A (zh) | 一种临近空间大展弦比柔性飞行器风洞试验数据修正方法 | |
CN110162822B (zh) | 耦合结构模态的时域快速非定常气动力计算方法 | |
CN104133926A (zh) | 一种弹性气动力特性综合分析方法 | |
CN109540459B (zh) | 一种气动特性数值计算结果修正方法 | |
CN112699622A (zh) | 一种固定翼无人机的紧密编队气动耦合效应建模方法 | |
Fugate et al. | Aero-Structural Modeling of the Truss-Braced Wing Aircraft Using Potential Method with Correction Methods for Transonic Viscous Flow and Wing-Strut Interference Aerodynamics | |
Zhang et al. | A morphing wing with cellular structure of non-uniform density | |
Bordogna et al. | Aeroelastic tailoring for static and dynamic loads with blending constraints | |
Huang et al. | Numerical method of static aeroelastic correction and jig-shape design for large airliners | |
Dai et al. | Research on dynamic stall active control of two-dimensional airfoil with combination of droop leading edge and trailing edge flap | |
Arena et al. | A tailored nonlinear slat-cove filler for airframe noise reduction | |
Hooker et al. | Design and transonic wind tunnel testing of a cruise efficient STOL military transport | |
Park | Advanced turboprop composite propeller design and analysis using fluid–structure interaction method | |
CN111959817A (zh) | 一种确定飞翼布局飞机机翼变形限制条件的方法 | |
Carreno Ruiz et al. | ROT8: A Matlab App for Low Reynolds number Airfoil and Rotor Aerodynamic Design | |
CN110532604A (zh) | 一种动态失速状态下带后缘小翼桨叶气动载荷的计算方法 | |
CN114117950B (zh) | 一种基于做功原理的天地往返飞行器颤振判断方法 | |
Zhang et al. | Analysis of wing flexure deformation based on ANSYS | |
Wang et al. | CFD-CSD method for rotor aeroelastic analysis with free wake model | |
CN113987793B (zh) | 一种飞机的俯仰力矩优化方法、装置、设备及存储介质 | |
Mishra et al. | Computational investigation of trailing edge flap for control of vibration | |
Boozer et al. | Coupled Adjoint Aerostructural Optimization Framework for Preliminary Aircraft Design |
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 |