CN109635370B - 开裂式阻力方向舵静气动弹性特性分析方法 - Google Patents
开裂式阻力方向舵静气动弹性特性分析方法 Download PDFInfo
- Publication number
- CN109635370B CN109635370B CN201811408265.XA CN201811408265A CN109635370B CN 109635370 B CN109635370 B CN 109635370B CN 201811408265 A CN201811408265 A CN 201811408265A CN 109635370 B CN109635370 B CN 109635370B
- Authority
- CN
- China
- Prior art keywords
- finite element
- aerodynamic
- force data
- grid
- aerodynamic force
- 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
- 230000003068 static effect Effects 0.000 title claims abstract description 25
- 238000004458 analytical method Methods 0.000 title abstract description 4
- 238000005336 cracking Methods 0.000 claims abstract description 16
- 238000000034 method Methods 0.000 claims description 30
- 238000004364 calculation method Methods 0.000 claims description 21
- 238000006073 displacement reaction Methods 0.000 claims description 14
- 239000012530 fluid Substances 0.000 claims description 9
- 239000011159 matrix material Substances 0.000 claims description 5
- 230000008878 coupling Effects 0.000 claims description 4
- 238000010168 coupling process Methods 0.000 claims description 4
- 238000005859 coupling reaction Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000007787 solid Substances 0.000 claims description 3
- 238000005516 engineering process Methods 0.000 claims description 2
- 230000007547 defect Effects 0.000 abstract description 2
- 230000008569 process Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 210000004027 cell Anatomy 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 239000000853 adhesive Substances 0.000 description 1
- 230000001070 adhesive effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000011438 discrete method Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000002763 pyramidal cell Anatomy 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
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)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本申请提供了一种开裂式阻力方向舵静气动弹性特性分析方法,能够克服现有静气动弹性设计手段的缺点和不足,实现对开裂式阻力方向舵的静气动弹性特性的准确预测,而且,还能够对飞机刚体气动力特性和弹性气动力特性进行精确模拟,实现全机静气动弹性特性的精细化设计。
Description
技术领域
本申请涉及静气动弹性设计领域,具体提供一种开裂式阻力方向舵静气动弹性特性分析方法。
背景技术
为达到优良的隐身性能,飞翼布局飞机取消了垂尾和常规方向舵,为保证横航向稳定性和优良的飞行品质,飞翼布局飞机常采用开裂式阻力方向舵来实现航向操纵。
与传统操纵舵面相比,开裂式阻力方向舵在外形和使用策略上均有明显差别,目前工程上较为成熟的静气动弹性设计方法均以低阶/高阶面元法为基础,这种方法无法对阻力方向舵的静气动弹性特性进行准确预测。
发明内容
为了解决上述技术问题至少之一,本申请提供了一种开裂式阻力方向舵静气动弹性特性分析方法,包括:步骤1:划分阻力方向舵单侧开裂的CFD网格,以及构建结构有限元模型;步骤2:计算刚体气动力数据以及飞机表面的气动载荷分布;步骤3:将所述气动载荷分布通过数值交换方法传递到所述结构有限元的节点上;步骤4:计算所述结构有限元的所述节点的位移;步骤5:将所述结构有限元的所述节点的位移传递回流体网格的物面节点上;步骤6:利用网格变形技术对流体计算域网格进行刷新;步骤7:重复步骤3到步骤6,直到计算结构收敛,得到弹性气动力数据;步骤8:针对不同计算状态,重复步骤2到步骤7,得到多个弹性气动力数据以及多个刚性气动力数据;步骤9:计算两个弹性气动力数据的差值以及两个刚性气动力数据的差值,并用两个弹性气动力数据的差值除以两个刚性气动力数据的差值,得到静气动弹性修正系数。
根据本申请的至少一个实施例,步骤2中:刚体气动力数据以及飞机表面的气动载荷分布通过下式计算:
在不计质量力、直角坐标系下,守恒形式的Navier-Stokes方程表示为:
其中,
应力项:
τyz=τzy=μ(vz+wy),τxz=τzx=μ(uz+wx),
热传导项:
单位质量气体的总能量
为了使Navier-Stokes方程封闭,还需要下述的关系式:
完全气体状态方程:p=ρRT,h=cpT,cp为等压比热,
粘性系数:μ=μl+μt,
其中,μl为层流粘性系数,μtt为湍流粘性系数。
根据本申请的至少一个实施例,步骤3中:通过三维空间插值法,将所述气动载荷分布传递到所述结构有限元的节点上,插值公式为:
对于三维插值,
径向基函数
同时附加方程:
求解得出待定参数(a1,a2,…,aN,aN+1,aN+2,aN+3,aN+4),即能够得到数据转换矩阵[G]和[G]T,其中:
Fs=HT·Fa。
根据本申请的至少一个实施例,步骤5中采用与步骤3中的有限元节点位移的插值方法相同。
根据本申请的至少一个实施例,步骤6中,通过多套重叠非结构背景网格插值法对流体计算域网格进行刷新。
根据本申请的至少一个实施例,步骤7中,以飞机结构有限元各节点位移为参考量,与上一迭代步骤相比,当所有节点的变形量小于1.0×10-5m时,认为流固耦合计算达到收敛。
本申请实施例提供的开裂式阻力方向舵静气动弹性特性分析方法,能够克服现有静气动弹性设计手段的缺点和不足,实现对开裂式阻力方向舵的静气动弹性特性的准确预测,而且,还能够对飞机刚体气动力特性和弹性气动力特性进行精确模拟,实现全机静气动弹性特性的精细化设计。
附图说明
图1是本申请实施例提供的CFD局部网格示意图;
图2是本申请实施例提供的CFD局部流场示意图;
图3是本申请实施例提供的考虑物面内外都覆盖的非结构背景网格示意图;
图4是本申请实施例提供的非结构网格变形前后示意图;
图5是本申请实施例提供的CFD/CSD迭代耦合计算求解流程示意图;
图6是本申请实施例提供的程序迭代计算求解过程的输出;
图7是本申请实施例提供的不同M数和速压下偏航力矩对阻力方向舵开裂角度的导数的弹性修正系数。
具体实施方式
下面结合附图和实施例对本申请作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释相关申请,而非对该申请的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与本申请相关的部分。
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
本申请实施例提供的开裂式阻力方向舵静气动弹性特性分析方法包括以下步骤:
步骤1:划分阻力方向舵单侧开裂的CFD网格,以及构建结构有限元模型。
飞翼布局飞机为获得良好的气动效率和隐身特性,取消了平尾、垂尾等部件,为解决取消垂尾带来的航向静稳定性和航向控制问题,采用了阻力方向舵进行航向操纵,阻力方向舵安装在机翼后缘翼稍附近,分上下两片,开裂式作动,通过两翼上的非对称阻力产生偏航力矩实现航向操纵。
使用CFD网格划分软件Gridgen或ICEM生成的非结构混合网格,网格由六面体单元、三棱柱单元、金字塔元和四面体单元组成。在边界层内采用各向异性的三棱柱和六面体网格,既能模拟边界层流动又节约了网格数量,边界层外的网格仍然采用四面体网格,可以做到与远场的快速过渡和复杂外形空间的填充,三棱柱和六面体网格与四面体网格的过渡则采用金字塔网格。混合网格相对于结构网格拥有复杂外形物体网格生成容易的优点,而相对于非结构网格则解决了无法处理边界层内流动的问题,适合于目前的大规模工程运用。CFD网格示意图如图1所示。
由全机结构有限元模型,基于有限元的承力结构限制和CFD/CSD插值精度的约束要求,选定若干个加载点(根据插值程序的计算量限制,一般选取1000个左右),使用MSC.Nastran软件对有限元模型进行逐点单位载荷加载,并记录所有加载点的位移数据,在所有的加载点都执行完加载计算任务之后,就可以得到一个包含所有加载点的柔度矩阵数据,数据按照节点坐标和单位载荷下位移量的格式进行输出备用。
步骤2:计算刚体气动力数据以及飞机表面的气动载荷分布。
计算刚体气动力数据以及飞机表面的气动载荷分布采用CFD方法,CFD计算可以基于混合网格Navier-Stokes方程的隐式求解,主要包括流体力学控制方程、空间离散方法、隐式时间推进格式、湍流模型的选取等。同时,CFD求解也可采用成熟的商业软件完成,如Fluent、CFD++、STAR CCM+等。
在不计质量力、直角坐标系下,守恒形式的Navier-Stokes方程表示为(式中的α等于0时,Navier-Stokes方程简化为Euler方程):
其中,
应力项:
τyz=τzy=μ(vz+wy),τxz=τzx=μ(uz+wx),
热传导项:
单位质量气体的总能量
为了使Navier-Stokes方程封闭,还需要下述的关系式:
完全气体状态方程:p=ρRT,h=cpT,cp为等压比热,
粘性系数:μ=μl+μt,
其中,μl为层流粘性系数,μt为湍流粘性系数。
层流粘性系数μl根据具体情况由Sutherland公式或Keys公式给定,湍流粘性系数μt由湍流模型给定。
在空间离散的处理上,选择了比较适合非结构混合网格的节点中心二阶离散格式,时间推进上选择了能够有效提高计算效率的LU-SGS隐式时间推进格式,湍流模型的选取上则选择了不需要求解节点到物面距离的k-ω两方程模型,图2给出了CFD计算的局部流场示意图。
步骤3:将所述气动载荷分布通过数值交换方法传递到所述结构有限元的节点上。
通过三维空间插值法,将所述气动载荷分布传递到所述结构有限元的节点上,插值公式为:
对于三维插值,
径向基函数
同时附加方程:
求解得出待定参数(a1,a2,…,aN,aN+1,aN+2,aN+3,aN+4),即能够得到数据转换矩阵[G]和[G]T,其中:
Fs=HT·Fa。
步骤4:计算所述结构有限元的所述节点的位移。
将有限元节点上的载荷加载在结构柔度矩阵上,得到结构的静变形,可以采用MSC.Nastran、Abaqus等有限元商业软件进行求解。
步骤5:将所述结构有限元的所述节点的位移传递回流体网格的物面节点上。
可以采用与步骤3中相同的有限元节点位移的插值方法。
步骤6:利用网格变形技术对流体计算域网格进行刷新。
通过多套重叠非结构背景网格插值法对流体计算域网格进行刷新,该方法只在背景网格部分采用弹簧变形,这样即使一个1000万的计算网格,如果背景网格只取30万,那么计算弹簧方程时也仅仅只需要求解这30万的网格,大大的提高了计算效率。图3给出了考虑物面内外都覆盖的非结构背景网格示意图。图4针对机翼扭转问题给出了机翼变形前后的网格变形示意图,可以看到变形后的网格质量仍保持较好。
步骤7:重复步骤3到步骤6,直到计算结构收敛,得到弹性气动力数据。
以飞机结构有限元各节点位移为参考量,与上一迭代步骤相比,当所有节点的变形量小于1.0×10-5m时,即认为流固耦合计算达到收敛。图5示出了CFD/CSD迭代耦合计算求解流程示意图。图6示出了程序迭代计算求解过程的输出。
步骤8:针对不同计算状态,重复步骤2到步骤7,得到多个弹性气动力数据以及多个刚性气动力数据。
步骤9:计算两个弹性气动力数据的差值以及两个刚性气动力数据的差值,并用两个弹性气动力数据的差值除以两个刚性气动力数据的差值,得到静气动弹性修正系数。
图7示出了本申请实施例提供的不同M数和速压下偏航力矩对阻力方向舵开裂角度的导数的弹性修正系数
以上所述,仅为本申请的具体实施方式,所属领域的技术人员可以清楚地了解到,为了描述的方便和简洁,上述描述的系统、模块和单元的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。应理解,本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到各种等效的修改或替换,这些修改或替换都应涵盖在本申请的保护范围之内。
Claims (4)
1.一种开裂式阻力方向舵静气动弹性特性分析方法,其特征在于,包括:
步骤1:划分阻力方向舵单侧开裂的CFD网格,以及构建结构有限元模型;
步骤2:计算刚体气动力数据以及飞机表面的气动载荷分布;
步骤3:将所述气动载荷分布通过数值交换方法传递到所述结构有限元的节点上;
步骤4:计算所述结构有限元的所述节点的位移;
步骤5:将所述结构有限元的所述节点的位移传递回流体网格的物面节点上;
步骤6:利用网格变形技术对流体计算域网格进行刷新;
步骤7:重复步骤3到步骤6,直到计算结构收敛,得到弹性气动力数据;
步骤8:针对不同计算状态,重复步骤2到步骤7,得到多个弹性气动力数据以及多个刚性气动力数据;
步骤9:计算两个弹性气动力数据的差值以及两个刚性气动力数据的差值,并用两个弹性气动力数据的差值除以两个刚性气动力数据的差值,得到静气动弹性修正系数;
步骤2中:
刚体气动力数据以及飞机表面的气动载荷分布通过下式计算:
在不计质量力、直角坐标系下,守恒形式的Navier-Stokes方程表示为:
其中,
应力项:
τyz=τzy=μ(vz+wy),τxz=τzx=μ(uz+wx),
热传导项:
单位质量气体的总能量
为了使Navier-Stokes方程封闭,还需要下述的关系式:
完全气体状态方程:p=ρRT,h=cpT,cp为等压比热,
粘性系数:μ=μl+μt,
其中,μl为层流粘性系数,μt为湍流粘性系数;
步骤3中:
通过三维空间插值法,将所述气动载荷分布传递到所述结构有限元的节点上,插值公式为:
对于三维插值,
径向基函数
同时附加方程:
求解得出待定参数(a1,a2,…,aN,aN+1,aN+2,aN+3,aN+4),即能够得到数据转换矩阵[H]和[H]T,其中:
Fs=HT·Fa。
2.根据权利要求1所述的开裂式阻力方向舵静气动弹性特性分析方法,其特征在于,步骤5中采用与步骤3中的有限元节点位移的插值方法相同。
3.根据权利要求1所述的开裂式阻力方向舵静气动弹性特性分析方法,其特征在于,步骤6中,通过多套重叠非结构背景网格插值法对流体计算域网格进行刷新。
4.根据权利要求1所述的开裂式阻力方向舵静气动弹性特性分析方法,其特征在于,步骤7中,以飞机结构有限元各节点位移为参考量,与上一迭代步骤相比,当所有节点的变形量小于1.0×10-5m时,认为流固耦合计算达到收敛。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811408265.XA CN109635370B (zh) | 2018-11-23 | 2018-11-23 | 开裂式阻力方向舵静气动弹性特性分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811408265.XA CN109635370B (zh) | 2018-11-23 | 2018-11-23 | 开裂式阻力方向舵静气动弹性特性分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109635370A CN109635370A (zh) | 2019-04-16 |
CN109635370B true CN109635370B (zh) | 2023-12-01 |
Family
ID=66069409
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811408265.XA Active CN109635370B (zh) | 2018-11-23 | 2018-11-23 | 开裂式阻力方向舵静气动弹性特性分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109635370B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112414668B (zh) * | 2020-11-03 | 2022-08-09 | 中国商用飞机有限责任公司北京民用飞机技术研究中心 | 一种风洞试验数据静气弹修正方法、装置、设备及介质 |
CN112541305B (zh) * | 2020-12-10 | 2024-01-02 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种基于全局变量求导的飞机气动力特性分析方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102012953B (zh) * | 2010-11-04 | 2013-05-08 | 西北工业大学 | Cfd/csd耦合求解非线性气动弹性仿真方法 |
CN103699756A (zh) * | 2014-01-06 | 2014-04-02 | 中国商用飞机有限责任公司 | 松耦合模拟系统及方法 |
CN108052772A (zh) * | 2017-12-30 | 2018-05-18 | 北京航空航天大学 | 一种基于结构降阶模型的几何非线性静气动弹性分析方法 |
-
2018
- 2018-11-23 CN CN201811408265.XA patent/CN109635370B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102012953B (zh) * | 2010-11-04 | 2013-05-08 | 西北工业大学 | Cfd/csd耦合求解非线性气动弹性仿真方法 |
CN103699756A (zh) * | 2014-01-06 | 2014-04-02 | 中国商用飞机有限责任公司 | 松耦合模拟系统及方法 |
CN108052772A (zh) * | 2017-12-30 | 2018-05-18 | 北京航空航天大学 | 一种基于结构降阶模型的几何非线性静气动弹性分析方法 |
Non-Patent Citations (9)
Title |
---|
"动网格技术及其在气动弹性数值";丁力;《中国博士学位论文全文数据库 (基础科学辑)》;20160615(第6期);全文 * |
"基于CFD/CSD耦合的结构几何非线性静气动弹性数值方法研究";聂雪媛等;《振动与冲击》;20160428;第35卷(第8期);第48-53页 * |
基于三维气动力的静气动弹性响应分析;张波成等;《工程力学》;20110725(第07期);全文 * |
基于两种不同方法的民用客机静气动弹性导数计算;谢萌等;《中国科技信息》;20160714(第14期);第87-89页 * |
基于修正面元法的机翼静气动弹性计算;贾欢等;《计算物理》;20140125(第01期);全文 * |
弹性前掠翼开裂式方向舵操纵特性分析;马斌麟等;《北京航空航天大学学报》;20170527(第09期);全文 * |
氢动力无人机大展弦比机翼静气弹特性分析;郝帅等;《北京航空航天大学学报》;20161116(第08期);全文 * |
静气动弹性计算方法研究;陈大伟等;《力学学报》;20090718(第04期);全文 * |
飞翼无人机静气弹参数分析及操纵效率计算;许军等;《西北工业大学学报》;20161015(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109635370A (zh) | 2019-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kenway et al. | A CAD-free approach to high-fidelity aerostructural optimization | |
CN104298805B (zh) | 高超声速飞行器cfd气动建模方法 | |
Keye | Fluid-structure coupled analysis of a transport aircraft and flight-test validation | |
CN109635370B (zh) | 开裂式阻力方向舵静气动弹性特性分析方法 | |
CN109918808B (zh) | 一种气热弹三场耦合仿真分析方法 | |
Wang et al. | Multi-body separation simulation with an improved general mesh deformation method | |
CN113723027A (zh) | 一种针对弹性飞机的静气动弹性计算方法 | |
Kroll et al. | Fluid-structure coupling for aerodynamic analysis and design: a DLR perspective | |
Ito et al. | TAS code, FaSTAR, and Cflow results for the sixth drag prediction workshop | |
CN107766620A (zh) | 一种基于降阶模型的气动‑热‑结构优化方法 | |
CN105335549A (zh) | 一种上下翼面变形叠加方法 | |
Koning et al. | Using RotCFD to Predict Isolated XV-15 Rotor Performance | |
CN105205208B (zh) | 一种飞机机翼型架外形设计方法 | |
Huang et al. | Numerical method of static aeroelastic correction and jig-shape design for large airliners | |
Guo et al. | Nonlinear gust response analysis of free flexible aircraft | |
CN114139285B (zh) | 一种柔性蒙皮机翼低雷诺数流固耦合分析方法 | |
Tinoco et al. | Drag prediction with the Zeus/CFL3D system | |
Tinoco | Summary data from the seventh AIAA CFD drag prediction workshop | |
Maxwell | Efficient Design of Viscous Waveriders with CFD Verification and Off-Design Performance Analysis | |
Zuo et al. | Efficient aeroelastic design optimization based on the discrete adjoint method | |
Lamberson et al. | High-Fidelity Aeroelastic Simulations with HPCMP CREATETM-AV Kestrel | |
Chen | Numerical calculations on the unsteady aerodynamic force of the tilt-rotor aircraft in conversion mode | |
van Rooij et al. | Generation of a reduced order model of an unmanned combat air vehicle using indicial response functions | |
Morton et al. | CFD-Based Reduced-Order Models of F-16XL Static and Dynamic Loads Using Kestrel | |
Liu et al. | Static aeroelasticity analysis of high aspect ratio wing and the jig-shape design based on multi-block structural grid |
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 |