CN111125829B - 一种优化全动平尾静气动弹性和颤振的方法 - Google Patents
一种优化全动平尾静气动弹性和颤振的方法 Download PDFInfo
- Publication number
- CN111125829B CN111125829B CN201911223855.XA CN201911223855A CN111125829B CN 111125829 B CN111125829 B CN 111125829B CN 201911223855 A CN201911223855 A CN 201911223855A CN 111125829 B CN111125829 B CN 111125829B
- Authority
- CN
- China
- Prior art keywords
- flutter
- full
- horizontal tail
- aeroelasticity
- dynamic
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
一种优化全动平尾静气动弹性和颤振的方法,首先建立有限元模型,而后对静气动弹性、颤振、配重逐一进行优化,匹配程度高,得到满足静气动弹性和颤振要求且不影响平尾外形的最小重量全动平尾,符合约束条件,有效解决飞机静气动弹性和颤振的问题。
Description
技术领域
本发明涉及飞行器气动弹性设计技术领域,尤其涉及一种优化全动平尾静气动弹性和颤振的方法。
背景技术
全动平尾的气动弹性往往限制飞行器的飞行性能,在刚度、质量和气动外形不断地发生变化中选择出符合气动弹性的设计方案,但在飞机全动平尾由于未考虑刚度、质量和气动外形之间匹配关系,导致先后出现静气动弹性和颤振问题,为此,需要对全动平尾的关键结构进行优化,以获取符合约束条件的结构方案。
发明内容
本发明所解决的技术问题在于提供一种优化全动平尾静气动弹性和颤振的方法,以解决上述背景技术中的问题。
本发明所解决的技术问题采用以下技术方案来实现:
一种优化全动平尾静气动弹性和颤振的方法,具体步骤如下:
1)建立有限元模型
有限元模型包括初始几何模型、全动平尾结构动力学模型及全动平尾气动力模型;并计算静气动弹性发散、颤振速度、颤振频率、全动平尾重量,同时以此计算结果作为基准;
2)静气动弹性优化
将摇臂的长度、宽度、厚度作为设计变量,以满足结构的强度和静气动弹性条件为约束条件,以结构重量为目标函数,通过优化计算,得到一个满足静气动弹性的最小重量摇臂;
3)颤振优化
在步骤2)的基础上,将平尾配重施加位置作为设计变量,以满足颤振条件作为约束条件,以结构重量为目标函数,通过优化设计计算,得到最小配重重量下满足颤振要求的全动平尾;
4)配重优化
将平尾内空间作为设计变量,在步骤3)的基础上,以满足静气动弹性、颤振条件作为约束条件,以结构重量为目标函数,通过优化计算,得到满足静气动弹性和颤振要求且不影响平尾外形的最小重量全动平尾。
在本发明中,通过颤振优化以及配重优化,在平尾根部前缘位置施加一个配重,使得颤振速度大幅提高,此配重为100克。
有益效果:
1)本发明同时综合平尾重量、静气动弹性、颤振多目标函数进行优化,匹配程度高;
2)本发明对静气动弹性优化仅以摇臂作为设计变量,尽可能使用高效的优化变量,使得飞机重量不需增加过多,即能满足静气动弹性问题;
3)本发明对颤振优化以配重位置作为设计变量,使得飞机满足颤振问题,
结合对静气动弹性优化,进而得到满足静气动弹性和颤振要求且不影响平尾外形的最小重量全动平尾,有效解决飞机静气动弹性和颤振的问题。
附图说明
图1为本发明的较佳实施例的有限元模型示意图。
图2为本发明的较佳实施例的流程示意图。
图3为本发明的较佳实施例中的摇臂优化示意图。
图4为本发明的较佳实施例中的颤振优化结点示意图。
具体实施方式
为了使本发明实现的技术手段、创作特征、达成目的与功效易于明白,下面结合具体图示,进一步阐述本发明。
参见图1~图4所示的一种优化全动平尾静气动弹性和颤振的方法,具体步骤如下:
1)建立有限元模型
建立全动平尾的初始几何模型,首先在三维软件如UG、catia中建立初始几何模型,而后将初始几何模型导入至MSC.Patran中,再在此基础中,去除多余的线面,建立好全动平尾的初始几何模型;
基于全动平尾的初始几何模型构建全动平尾结构动力学模型,将全动平尾使用材料的属性填入即可,质量特性通过密度赋予;全动平尾气动力模型根据全动平尾几何外形用CAERO1卡片建立即可,全动平尾结构动力学模型与全动平尾气动力模型之间通过样条差值SPLINE卡片连接;
静气动弹性与颤振通过MSC.Nastran软件计算,静气动弹性发散通过在bdf文件中设置DIVERG卡,颤振通过在bdf文件中设置FLUTTER卡,由此完成有限元模型建立;
2)静气动弹性优化
使用求解序列SOL200,采用DESVAR卡片定义摇臂长度、宽度、厚度作为设计变量,采用DRESP1卡片定义重量响应为目标函数,采用DCONSTR卡片定义静气动弹性为约束函数,通过优化计算,得到一个满足静气动弹性的最小重量摇臂;
3)颤振优化
在步骤2)的基础上,以配重位置作为设计变量,通过FORTRAN软件编写批量读写运算程序,在步骤1)建立的全动平尾结构动力学模型结点p上施加同样重量的配重,得到在不同位置对颤振速度的敏感度矩阵,从而得到在某一特定位置施加最小配重即满足颤振要求,若同时满足全动平尾几何外形要求,则优化设计结束;否则,进入步骤4);
4)配重优化
如若步骤3)中优化全动平尾结果不能满足其几何外形要求,在步骤3)中基础上,将配重重量大小、全动平尾内空间作为设计变量、以满足颤振为约束条件,将结构重量作为目标函数,通过优化设计,得到既能满足静气动弹性要求又能满足颤振要求的最小重量全动平尾。
在本实施例中,全动平尾的基本情况:平尾操纵系统刚度为80kg·m/rad,平尾旋转频率为25Hz,平尾一弯频率为40Hz。
在本实施例中,通过静气动弹性和颤振计算,平尾操纵系统刚度偏低,导致静气动弹性发散速度低于设计要求,在提高操纵系统刚度后,旋转和一弯近频,导致出现颤振问题,通过静气动弹性优化计算,发现摇臂厚度对发散影响非常敏感,故通过增加摇臂厚度以显著提高操纵系统刚度,从而提高发散速度,但由于该全动平尾旋转频率低于平尾一弯频率,提高操纵系统刚度,使得耦合模态靠近,从而出现颤振问题,通过颤振优化以及配重优化,在平尾根部前缘位置施加一个约100克配重,使得颤振速度大幅提高。
Claims (10)
1. 一种优化全动平尾静气动弹性和颤振的方法,其特征在于,具体步骤如下:
1)建立有限元模型
有限元模型包括初始几何模型、全动平尾结构动力学模型及全动平尾气动力模型;并计算静气动弹性发散、颤振速度、颤振频率、全动平尾重量,同时以此计算结果作为基准;
2)静气动弹性优化
将摇臂的长度、宽度、厚度作为设计变量,以满足结构的强度和静气动弹性条件为约束条件,以结构重量为目标函数,通过优化计算,得到一个满足静气动弹性的最小重量摇臂;
3)颤振优化
在步骤2)的基础上,将平尾配重施加位置作为设计变量,以满足颤振条件作为约束条件,以结构重量为目标函数,通过优化设计计算,得到最小配重重量下满足颤振要求的全动平尾;
4)配重优化
将平尾内空间作为设计变量,在步骤3)的基础上,以满足静气动弹性、颤振条件作为约束条件,以结构重量为目标函数,通过优化计算,得到满足静气动弹性和颤振要求且不影响平尾外形的最小重量全动平尾。
2.根据权利要求1所述的一种优化全动平尾静气动弹性和颤振的方法,其特征在于,步骤1)中,建立全动平尾的初始几何模型,首先在三维软件中建立初始几何模型,而后将初始几何模型导入至MSC.Patran中,再在此基础中,去除多余的线面即可。
3.根据权利要求2所述的一种优化全动平尾静气动弹性和颤振的方法,其特征在于,步骤1)中,基于全动平尾的初始几何模型构建全动平尾结构动力学模型,将全动平尾使用材料的属性填入即可,质量特性通过密度赋予。
4.根据权利要求2所述的一种优化全动平尾静气动弹性和颤振的方法,其特征在于,步骤1)中,全动平尾气动力模型根据全动平尾几何外形用CAERO1卡片建立即可。
5.根据权利要求1所述的一种优化全动平尾静气动弹性和颤振的方法,其特征在于,步骤1)中,静气动弹性与颤振通过MSC.Nastran软件计算,静气动弹性发散通过在bdf文件中设置DIVERG卡,颤振通过在bdf文件中设置FLUTTER卡。
6.根据权利要求1所述的一种优化全动平尾静气动弹性和颤振的方法,其特征在于,步骤2)中,使用求解序列SOL200,采用DESVAR卡片定义摇臂长度、宽度、厚度作为设计变量,采用DRESP1卡片定义重量响应为目标函数,采用DCONSTR卡片定义静气动弹性为约束函数,通过优化计算,得到一个满足静气动弹性的最小重量摇臂。
7.根据权利要求1所述的一种优化全动平尾静气动弹性和颤振的方法,其特征在于,步骤3)中,通过FORTRAN软件编写批量读写运算程序,在步骤1)建立的全动平尾结构动力学模型结点上施加同样重量的配重,得到在不同位置对颤振速度的敏感度矩阵,从而得到在某一特定位置施加最小配重即满足颤振要求,若同时满足全动平尾几何外形要求,则优化设计结束;否则,进入步骤4)。
8.根据权利要求1所述的一种优化全动平尾静气动弹性和颤振的方法,其特征在于,步骤4)中,如若步骤3)中优化全动平尾结果不能满足其几何外形要求,则在步骤3)中基础上,将配重重量大小、全动平尾内空间作为设计变量、以满足颤振为约束条件,将结构重量作为目标函数,通过优化设计,得到既能满足静气动弹性要求又能满足颤振要求的最小重量全动平尾。
9.根据权利要求8所述的一种优化全动平尾静气动弹性和颤振的方法,其特征在于,在平尾根部前缘位置施加一个配重,以使得颤振速度大幅提高。
10.根据权利要求9所述的一种优化全动平尾静气动弹性和颤振的方法,其特征在于,配重为100克。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911223855.XA CN111125829B (zh) | 2019-12-04 | 2019-12-04 | 一种优化全动平尾静气动弹性和颤振的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911223855.XA CN111125829B (zh) | 2019-12-04 | 2019-12-04 | 一种优化全动平尾静气动弹性和颤振的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111125829A CN111125829A (zh) | 2020-05-08 |
CN111125829B true CN111125829B (zh) | 2022-05-06 |
Family
ID=70497427
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911223855.XA Active CN111125829B (zh) | 2019-12-04 | 2019-12-04 | 一种优化全动平尾静气动弹性和颤振的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111125829B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112926141B (zh) * | 2021-03-31 | 2022-06-21 | 同济大学 | 军机平尾转角伺服控制系统的多目标优化协同设计方法 |
CN113665838B (zh) * | 2021-10-21 | 2021-12-21 | 中国空气动力研究与发展中心低速空气动力研究所 | 一种拓宽直升机重心包线的全动平尾布局优化方法及装置 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6675735B1 (en) * | 1998-11-02 | 2004-01-13 | Stephen Bourn | Hydrofoil sail craft |
CN106507753B (zh) * | 2004-04-13 | 2010-08-04 | 成都飞机工业(集团)有限责任公司 | 轻型战斗机机翼 |
CN103077286A (zh) * | 2013-01-18 | 2013-05-01 | 大连理工大学 | 一种飞机颤振模型的频率误差修正方法 |
CN103186148A (zh) * | 2011-12-28 | 2013-07-03 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种转速测量控制系统 |
CN103646131A (zh) * | 2013-11-26 | 2014-03-19 | 北京航空航天大学 | 一种考虑气动弹性约束的复合材料机翼多目标优化设计方法 |
CN104133933A (zh) * | 2014-05-29 | 2014-11-05 | 温州职业技术学院 | 一种高超声速飞行器热环境下气动弹性力学特性分析方法 |
CN105335578A (zh) * | 2015-11-25 | 2016-02-17 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种针对静气动弹性分析的型架外形的设计方法 |
CN105354348A (zh) * | 2015-08-14 | 2016-02-24 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 一种机翼低速颤振风洞模型的制造方法 |
US9541106B1 (en) * | 2007-01-03 | 2017-01-10 | Orbitel Research Inc. | Plasma optimized aerostructures for efficient flow control |
CN108052772A (zh) * | 2017-12-30 | 2018-05-18 | 北京航空航天大学 | 一种基于结构降阶模型的几何非线性静气动弹性分析方法 |
CN108363843A (zh) * | 2018-01-25 | 2018-08-03 | 北京航空航天大学 | 一种基于结构降阶模型的几何非线性静气动弹性全机配平方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2953302B1 (fr) * | 2009-11-27 | 2012-08-10 | Thales Sa | Procede de planification, de calcul de trajectoire, de predictions et de guidage pour le respect d'une contrainte de temps de passage d'un aeronef |
US9709029B2 (en) * | 2011-06-21 | 2017-07-18 | University Of Virginia Patent Foundation | Morphing segmented wind turbine and related method |
-
2019
- 2019-12-04 CN CN201911223855.XA patent/CN111125829B/zh active Active
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6675735B1 (en) * | 1998-11-02 | 2004-01-13 | Stephen Bourn | Hydrofoil sail craft |
CN106507753B (zh) * | 2004-04-13 | 2010-08-04 | 成都飞机工业(集团)有限责任公司 | 轻型战斗机机翼 |
US9541106B1 (en) * | 2007-01-03 | 2017-01-10 | Orbitel Research Inc. | Plasma optimized aerostructures for efficient flow control |
CN103186148A (zh) * | 2011-12-28 | 2013-07-03 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种转速测量控制系统 |
CN103077286A (zh) * | 2013-01-18 | 2013-05-01 | 大连理工大学 | 一种飞机颤振模型的频率误差修正方法 |
CN103646131A (zh) * | 2013-11-26 | 2014-03-19 | 北京航空航天大学 | 一种考虑气动弹性约束的复合材料机翼多目标优化设计方法 |
CN104133933A (zh) * | 2014-05-29 | 2014-11-05 | 温州职业技术学院 | 一种高超声速飞行器热环境下气动弹性力学特性分析方法 |
CN105354348A (zh) * | 2015-08-14 | 2016-02-24 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 一种机翼低速颤振风洞模型的制造方法 |
CN105335578A (zh) * | 2015-11-25 | 2016-02-17 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种针对静气动弹性分析的型架外形的设计方法 |
CN108052772A (zh) * | 2017-12-30 | 2018-05-18 | 北京航空航天大学 | 一种基于结构降阶模型的几何非线性静气动弹性分析方法 |
CN108363843A (zh) * | 2018-01-25 | 2018-08-03 | 北京航空航天大学 | 一种基于结构降阶模型的几何非线性静气动弹性全机配平方法 |
Non-Patent Citations (2)
Title |
---|
"Flutter analysis of ARW-2 wing using doublet lattice";M. J. Aziz 等;《2018 15th International Bhurban Conference on Applied Sciences and Technology 》;20180312;第579-584页 * |
"机翼带外挂高阶高频颤振特性研究";陈立勇 等;《装备环境工程》;20180915;第16-20页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111125829A (zh) | 2020-05-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110610065B (zh) | 基于混合动网格技术的飞行器多体分离cfd仿真方法及系统 | |
CN111125829B (zh) | 一种优化全动平尾静气动弹性和颤振的方法 | |
JP6121064B2 (ja) | 付加製造のための格子構造の構造保存性のトポロジ最適化方法 | |
US11625516B2 (en) | Designing a part by topology optimization | |
JP2020071887A (ja) | トポロジー最適化による機械部品の設計 | |
CN100543754C (zh) | 用于设计飞行交通器的方法 | |
CN109614631B (zh) | 基于强化学习和迁移学习的飞行器全自动气动优化方法 | |
CN109543243A (zh) | 基于adams的舰载机拦阻着舰动力学仿真分析方法 | |
US11481634B2 (en) | Systems and methods for training a neural network to control an aircraft | |
CN110032797A (zh) | 无人机uav控制律调参方法 | |
CN102609990B (zh) | 面向复杂三维cad模型的海量场景渐进式更新算法 | |
JP2018180830A (ja) | インテーク設計方法、インテーク設計プログラム及びインテーク設計装置 | |
CN108920811B (zh) | 一种用于直升机飞行仿真的模拟方法及系统 | |
CN108363843A (zh) | 一种基于结构降阶模型的几何非线性静气动弹性全机配平方法 | |
CN103678762A (zh) | 优化的复合材料机翼气动弹性风洞模型的缩比建模方法 | |
US9187187B2 (en) | Method of calculating dynamic pressure at the level of an aircraft surface | |
CN109815518A (zh) | 基于转动惯量约束的飞行器舵面设计方法 | |
CN112733276A (zh) | 倾转旋翼无人机的倾转机构动力学优化方法及系统 | |
CN113761646A (zh) | 一种用于移动风场环境中飞行器动响应的确定方法 | |
CN103177165A (zh) | 客车车身结构设计系统、客车侧翻仿真测试系统及方法 | |
Garcia Perez et al. | High-fidelity simulations of a tiltwing vehicle for urban air mobility | |
JP7014472B2 (ja) | 三次元特徴記述子に基づく機翼フォーク点群の分割方法 | |
McDaniel et al. | Efficient mesh deformation for computational stability and control analyses on unstructured viscous meshes | |
CN109522636B (zh) | 一种用于直升机抗坠毁分析的轮式起落架简化方法 | |
CN109635370B (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 |