CN104573296B - 一种面向相似网格的高超声速流场初始化方法 - Google Patents
一种面向相似网格的高超声速流场初始化方法 Download PDFInfo
- Publication number
- CN104573296B CN104573296B CN201310498573.7A CN201310498573A CN104573296B CN 104573296 B CN104573296 B CN 104573296B CN 201310498573 A CN201310498573 A CN 201310498573A CN 104573296 B CN104573296 B CN 104573296B
- Authority
- CN
- China
- Prior art keywords
- mrow
- grid
- msubsup
- rho
- aerodynamics
- 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
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于高超声速气动流场数值模拟计算技术领域,具体涉及一种面向相似网格的高超声速气动流场初始化方法。包括以下步骤:(1)确定流场初始化和网格单元一致性判别原则;(2)具体进行面向相似网格的高超声速流场初始化:(2.1)设定待求解气动流场网格和已经完成数值模拟的气动流场网格;(2.2)具体实施步骤。本发明与常用的基于远场的初始化方法相比,提高了高超声速气动流场数值模拟的收敛速度,可显著缩短仅马赫数变化时整批次的工况计算总运行时间,同时还提高了高超声速气动流场数值模拟的正确性。
Description
技术领域
本发明属于高超声速气动流场数值模拟计算技术领域,具体涉及一种面向相似网格的高超声速气动流场初始化方法。
背景技术
高超声速科技已经成为二十一世纪航空航天领域的制高点,具有广阔的军民两用前景。高超声速科技的突破将对一个国家科学技术和国民经济的发展、综合国力的提升产生重大影响。在飞行器设计中,气动布局是影响飞行器成败的关键因素。由于风洞实验难以模拟10马赫数以上的气动流场的状况,CFD(计算流体力学)数值模拟成为研究气动力布局的主要手段之一。
气动流场网格由飞行器外形决定。在高超声速飞行器气动设计后期,飞行器外形B与飞行器外形A相比只有细微的调整,从而称B为A的相似外形。飞行器外形B与A相比只是在尾部稍长一些,整体上A、B非常相似,如图1所示。针对外形A、B由Gridgen或者ICEM软件生成的流场网格只在尾部有一些差别,从而称对应外形B的网格GB为对应A的网格GA的相似网格。GA和GB之间的差异只在图1中有阴影的部分,即GB为GA的相似网格。流场网格的变化最终影响气动流场的性质。在高超声速领域,由于缺乏实验手段,网格的细微调整引起的气动力差异也需要对所有状态进行数值模拟。
由于高超声速气动流场具有运行时间长,需要模拟的状态多的特点,且高性能计算资源通常是有限和宝贵的,从而需要节约总的数值模拟时间,节省计算资源。气动流场初场是影响高超声速气动流场数值模拟收敛速度的核心因素之一。一个好的初场可以使气动流场数值模拟在很短的时间内收敛。一个不好的初场则可能导致气动流场数值模拟收敛速度极慢,甚至发散。由于初场的好坏难以量化,国内外采用流场初始化主要为基于无穷远场的初始化。
流场网格格心集合为G={(i,j,k)},i,j,k满足其中I,J,K分别为X,Y,Z三个方向网格划分数目。气动流场F的守恒量主要包含如下5个: ρ为流场密度,u,v,w为速度在X,Y,Z三个方向的分量,e为能量。气动流场求解的是N-S方程,服从质量、动量、通量三大物理量的守恒定律。
基于无穷远场的初始化:假设无穷远场守恒量的值(ρ)∞,(ρu)∞,(ρv)∞,(ρw)∞,(ρe)∞,对所有气动流场网格格心都赋予的值即(ρ)i,j,k=(ρ)∞,(ρu)i,j,k=(ρu)∞,(ρv)i,j,k=(ρv)∞,(ρw)i,j,k=(ρw)∞,(ρe)i,j,k=(ρe)∞,如图2所示。这是在工程实践中最简单和最常用的方法。但是在高超声速飞行中,无穷远场的守恒量与计算收敛后得到的守恒量相差太远导致计算时间过长,甚至不收敛。
总之,针对相似网格,基于无穷远场的初始化存在收敛速度慢甚至不收敛的问题,亟需研制一种面向相似网格的高超声速气动流场初始化方法,从而解决相似网格高超声速气动流场初始化的问题。
发明内容
本发明要解决的技术问题是提供一种流场初始化方法,以提高相似网格高超声速气动流场的收敛速度,缩短运行时间。
为了实现这一目的,本发明采取的技术方案是:
一种面向相似网格的高超声速流场初始化方法,包括以下步骤:
(1)确定流场初始化和网格单元一致性判别原则:
设定气动流场网格单元集合为G={(i,j,k)},i=1,2,...,I;j=1,2,...,J;k=1,2,...,K,I,J,K分别为X,Y,Z三个方向上的网格划分数目;
气动流场F的的守恒量主要包含如下5个:网格单元(i,j,k)的流场密度速度在X,Y,Z三个方向上的分量u,v,w与网格单元(i,j,k)的密度的乘积:能量e与网格单元(i,j,k)的密度的乘积:
每一个格心对应一个流场网格单元,对于结构化网格为三维的长方体,每个网格单元包含8个网格点,分别是长方体的8个顶点的坐标信息;
假设气动流场网格GB为气动流场网格GA的相似网格,根据8个网格点的坐标信息是否完全相同来确定GB中的网格单元E是否为GA中的网格单元;如果GA中存在8个网格点坐标信息与GB中的网格单元E相同的网格单元,则认为GB中网格单元E为GA中的网格单元;否则认为GB中网格单元E不是GA中的网格单元;
(2)具体进行面向相似网格的高超声速流场初始化:
(2.1)设定待求解气动流场网格和已经完成数值模拟的气动流场网格:
设定待求解气动流场网格GB中的网格单元集合为GB={b(i,j,k)},i=1,2,...,IB;j=1,2,...,JB;k=1,2,...,KB,其中IB,JB,KB分别为气动流场网格GB中X,Y,Z三个方向网格划分数目;
设定已经完成数值模拟的气动流场网格GA中的网格单元集合为GA={a(i1,j1,k1)},i1,j1,k1满足i1=1,2,...,IA;j1=1,2,...,JA;k1=1,2,...,KA,其中IA,JA,KA分别为气动流场网格GA中X,Y,Z三个方向网格划分数目;
GB和GA互为相似流场;GA的流场已知,需要对GB的流场赋初值,无穷远场守恒量的值为(ρ)∞,(ρu)∞,(ρv)∞,(ρw)∞,(ρe)∞;
(2.2)具体实施步骤如下:
步骤1:设定i=1,j=1,k=1;
步骤2:如果GA中存在与GB中网格单元b(i,j,k)完全一致的网格a(i1,j1,k1),转步骤3,否则转步骤8;
步骤3:
步骤4:
步骤5:
步骤6:
步骤7:
步骤8:
步骤9:
步骤10:
步骤11:
步骤12:
步骤13:k=k+1;
步骤14:如果k<=K满足,转步骤2,否则转15;
步骤15:j=j+1;
步骤16:如果j<=J满足,转步骤2,否则转17;
步骤17:i=i+1;
步骤18:如果i<=I满足,转步骤2,否则转19;
步骤19:相似网格GB流场初始化结束。
本发明与常用的基于远场的初始化方法相比,提高了高超声速气动流场数值模拟的收敛速度,可显著缩短仅马赫数变化时整批次的工况计算总运行时间。同时还提高了高超声速气动流场数值模拟的正确性。
附图说明
图1是相似网格示意图;
图2是基于无穷远场的流场初始化示意图;
图3是网格单元示意图。
具体实施方式
下面结合附图对本发明技术方案进行详细说明。
本发明一种面向相似网格的高超声速流场初始化方法,包括以下步骤:
(1)确定流场初始化和网格单元一致性判别原则:
设定气动流场网格单元集合为G={(i,j,k)},i=1,2,...,I;j=1,2,...,J;k=1,2,...,K,I,J,K分别为X,Y,Z三个方向上的网格划分数目;
气动流场F的的守恒量主要包含如下5个:网格单元(i,j,k)的流场密度速度在X,Y,Z三个方向上的分量u,v,w与网格单元(i,j,k)的密度的乘积:能量e与网格单元(i,j,k)的密度的乘积:
每一个格心对应一个流场网格单元,对于结构化网格为三维的长方体,每个网格单元包含如图3所示的8个网格点,分别是长方体的8个顶点的坐标信息;
假设气动流场网格GB为气动流场网格GA的相似网格,根据8个网格点的坐标信息是否完全相同来确定GB中的网格单元E是否为GA中的网格单元;如果GA中存在8个网格点坐标信息与GB中的网格单元E相同的网格单元,则认为GB中网格单元E为GA中的网格单元;否则认为GB中网格单元E不是GA中的网格单元;
(2)具体进行面向相似网格的高超声速流场初始化:
(2.1)设定待求解气动流场网格和已经完成数值模拟的气动流场网格:
设定待求解气动流场网格GB中的网格单元集合为GB={b(i,j,k)},i=1,2,...,IB;j=1,2,...,JB;k=1,2,...,KB,其中IB,JB,KB分别为气动流场网格GB中X,Y,Z三个方向网格划分数目;
设定已经完成数值模拟的气动流场网格GA中的网格单元集合为GA={a(i1,j1,k1)},i1,j1,k1满足i1=1,2,...,IA;j1=1,2,...,JA;k1=1,2,...,KA,其中IA,JA,KA分别为气动流场网格GA中X,Y,Z三个方向网格划分数目;
GB和GA互为相似流场;GA的流场已知,需要对GB的流场赋初值,无穷远场守恒量的值为(ρ)∞,(ρu)∞,(ρv)∞,(ρw)∞,(ρe)∞;
(2.2)具体实施步骤如下:
步骤1:设定i=1,j=1,k=1;
步骤2:如果GA中存在与GB中网格单元b(i,j,k)完全一致的网格a(i1,j1,k1),转步骤3,否则转步骤8;
步骤3:
步骤4:
步骤5:
步骤6:
步骤7:
步骤8:
步骤9:
步骤10:
步骤11:
步骤12:
步骤13:k=k+1;
步骤14:如果k<=K满足,转步骤2,否则转15;
步骤15:j=j+1;
步骤16:如果j<=J满足,转步骤2,否则转17;
步骤17:i=i+1;
步骤18:如果i<=I满足,转步骤2,否则转19;
步骤19:相似网格GB流场初始化结束。
Claims (1)
1.一种面向相似网格的高超声速流场初始化方法,其特征在于,包括以下步骤:
(1)确定流场初始化和网格单元一致性判别原则:
设定气动流场网格单元集合为G={(i,j,k)},i=1,2,...,I;j=1,2,...,J;k=1,2,...,K,I,J,K分别为X,Y,Z三个方向上的网格划分数目;
气动流场F的的守恒量主要包含如下5个:网格单元(i,j,k)的流场密度速度在X,Y,Z三个方向上的分量u,v,w与网格单元(i,j,k)的密度的乘积:能量e与网格单元(i,j,k)的密度的乘积:
每一个格心对应一个流场网格单元,对于结构化网格为三维的长方体,每个网格单元包含8个网格点,分别是长方体的8个顶点的坐标信息;
假设气动流场网格GB为气动流场网格GA的相似网格,根据8个网格点的坐标信息是否完全相同来确定GB中的网格单元E是否为GA中的网格单元;如果GA中存在8个网格点坐标信息与GB中的网格单元E相同的网格单元,则认为GB中网格单元E为GA中的网格单元;否则认为GB中网格单元E不是GA中的网格单元;
(2)具体进行面向相似网格的高超声速流场初始化:
(2.1)设定待求解气动流场网格和已经完成数值模拟的气动流场网格:
设定待求解气动流场网格GB中的网格单元集合为GB={b(i,j,k)},i=1,2,...,IB;j=1,2,...,JB;k=1,2,...,KB,其中IB,JB,KB分别为气动流场网格GB中X,Y,Z三个方向网格划分数目;
设定已经完成数值模拟的气动流场网格GA中的网格单元集合为GA={a(i1,j1,k1)},i1,j1,k1满足i1=1,2,...,IA;j1=1,2,...,JA;k1=1,2,...,KA,其中IA,JA,KA分别为气动流场网格GA中X,Y,Z三个方向网格划分数目;
GB和GA互为相似流场;GA的流场已知,需要对GB的流场赋初值,无穷远场守恒量的值为(ρ)∞,(ρu)∞,(ρv)∞,(ρw)∞,(ρe)∞;
(2.2)具体实施步骤如下:
步骤1:设定i=1,j=1,k=1;
步骤2:如果GA中存在与GB中网格单元b(i,j,k)完全一致的网格a(i1,j1,k1),转步骤3,否则转步骤8;
步骤3:
<mrow>
<msubsup>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mi>B</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mn>1</mn>
<mo>,</mo>
<mi>k</mi>
<mn>1</mn>
</mrow>
<mi>A</mi>
</msubsup>
<mo>;</mo>
</mrow>
步骤4:
<mrow>
<msubsup>
<mrow>
<mo>(</mo>
<mi>&rho;u</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mi>B</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mrow>
<mo>(</mo>
<mi>&rho;u</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mn>1</mn>
<mo>,</mo>
<mi>k</mi>
<mn>1</mn>
</mrow>
<mi>A</mi>
</msubsup>
<mo>;</mo>
</mrow>
步骤5:
<mrow>
<msubsup>
<mrow>
<mo>(</mo>
<mi>&rho;v</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mi>B</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mrow>
<mo>(</mo>
<mi>&rho;v</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mn>1</mn>
<mo>,</mo>
<mi>k</mi>
<mn>1</mn>
</mrow>
<mi>A</mi>
</msubsup>
<mo>;</mo>
</mrow>
步骤6:
<mrow>
<msubsup>
<mrow>
<mo>(</mo>
<mi>&rho;w</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mi>B</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mrow>
<mo>(</mo>
<mi>&rho;w</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mn>1</mn>
<mo>,</mo>
<mi>k</mi>
<mn>1</mn>
</mrow>
<mi>A</mi>
</msubsup>
<mo>;</mo>
</mrow>
1
步骤7:
<mrow>
<msubsup>
<mrow>
<mo>(</mo>
<mi>&rho;e</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mi>B</mi>
</msubsup>
<mo>=</mo>
<msubsup>
<mrow>
<mo>(</mo>
<mi>&rho;e</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mn>1</mn>
<mo>,</mo>
<mi>j</mi>
<mn>1</mn>
<mo>,</mo>
<mi>k</mi>
<mn>1</mn>
</mrow>
<mi>A</mi>
</msubsup>
<mo>;</mo>
</mrow>
步骤8:
<mrow>
<msubsup>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mi>B</mi>
</msubsup>
<mo>=</mo>
<msub>
<mrow>
<mo>(</mo>
<mi>&rho;</mi>
<mo>)</mo>
</mrow>
<mo>&infin;</mo>
</msub>
<mo>;</mo>
</mrow>
步骤9:
<mrow>
<msubsup>
<mrow>
<mo>(</mo>
<mi>&rho;u</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mi>B</mi>
</msubsup>
<mo>=</mo>
<msub>
<mrow>
<mo>(</mo>
<mi>&rho;u</mi>
<mo>)</mo>
</mrow>
<mo>&infin;</mo>
</msub>
<mo>;</mo>
</mrow>
步骤10:
<mrow>
<msubsup>
<mrow>
<mo>(</mo>
<mi>&rho;v</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mi>B</mi>
</msubsup>
<mo>=</mo>
<msub>
<mrow>
<mo>(</mo>
<mi>&rho;v</mi>
<mo>)</mo>
</mrow>
<mo>&infin;</mo>
</msub>
<mo>;</mo>
</mrow>
步骤11:
<mrow>
<msubsup>
<mrow>
<mo>(</mo>
<mi>&rho;w</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mi>B</mi>
</msubsup>
<mo>=</mo>
<msub>
<mrow>
<mo>(</mo>
<mi>&rho;w</mi>
<mo>)</mo>
</mrow>
<mo>&infin;</mo>
</msub>
<mo>;</mo>
</mrow>
步骤12:
<mrow>
<msubsup>
<mrow>
<mo>(</mo>
<mi>&rho;e</mi>
<mo>)</mo>
</mrow>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
<mi>B</mi>
</msubsup>
<mo>=</mo>
<msub>
<mrow>
<mo>(</mo>
<mi>&rho;e</mi>
<mo>)</mo>
</mrow>
<mo>&infin;</mo>
</msub>
<mo>;</mo>
</mrow>
步骤13:k=k+1;
步骤14:如果k<=K满足,转步骤2,否则转15;
步骤15:j=j+1;
步骤16:如果j<=J满足,转步骤2,否则转17;
步骤17:i=i+1;
步骤18:如果i<=I满足,转步骤2,否则转19;
步骤19:相似网格GB流场初始化结束。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310498573.7A CN104573296B (zh) | 2013-10-22 | 2013-10-22 | 一种面向相似网格的高超声速流场初始化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310498573.7A CN104573296B (zh) | 2013-10-22 | 2013-10-22 | 一种面向相似网格的高超声速流场初始化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104573296A CN104573296A (zh) | 2015-04-29 |
CN104573296B true CN104573296B (zh) | 2017-09-22 |
Family
ID=53089349
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310498573.7A Active CN104573296B (zh) | 2013-10-22 | 2013-10-22 | 一种面向相似网格的高超声速流场初始化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104573296B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106339534B (zh) * | 2016-08-18 | 2019-04-12 | 中国人民解放军国防科学技术大学 | 一种用于飞行器定常绕流数值求解的流场初始化方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20090035059A (ko) * | 2007-10-05 | 2009-04-09 | 건국대학교 산학협력단 | 초음속 고받음각 유동을 위한 난류 모델 분석 방법 |
CN103226540A (zh) * | 2013-05-21 | 2013-07-31 | 中国人民解放军国防科学技术大学 | 基于分组多流的gpu上多区结构网格cfd加速方法 |
CN103345580A (zh) * | 2013-07-02 | 2013-10-09 | 上海大学 | 基于格子Boltzmann方法的并行CFD方法 |
-
2013
- 2013-10-22 CN CN201310498573.7A patent/CN104573296B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20090035059A (ko) * | 2007-10-05 | 2009-04-09 | 건국대학교 산학협력단 | 초음속 고받음각 유동을 위한 난류 모델 분석 방법 |
CN103226540A (zh) * | 2013-05-21 | 2013-07-31 | 中国人民解放军国防科学技术大学 | 基于分组多流的gpu上多区结构网格cfd加速方法 |
CN103345580A (zh) * | 2013-07-02 | 2013-10-09 | 上海大学 | 基于格子Boltzmann方法的并行CFD方法 |
Non-Patent Citations (3)
Title |
---|
二维结构化网格CFD LU-SGS时间推进并行算法;龚春叶 等;《计算机科学与探索》;20130701;第7卷(第10期);第916-923页 * |
基于CFD数值模拟的拉瓦尔喷管流场分析;王平 等;《航空计算技术》;20120715;第42卷(第4期);第60-62页 * |
高超声速飞行器气动力气动热数值模拟和超声速流动的区域推进求解;贺旭照;《中国博士学位论文全文数据库 工程科技Ⅱ辑》;20071215(第6期);第C031-4页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104573296A (zh) | 2015-04-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103106301B (zh) | 一种基于蒙特卡罗方法与特征线方法耦合的计算辐射屏蔽的方法 | |
CN106529093A (zh) | 一种针对大展弦比机翼的气动/结构/静气弹耦合优化方法 | |
CN104881510A (zh) | 一种直升机旋翼/尾桨气动干扰数值仿真方法 | |
CN103902782A (zh) | 基于pod和代理模型的高超声速气动热模型降阶方法 | |
CN103558020B (zh) | 机翼活动翼面在机翼大变形时的试验载荷施加方法 | |
CN103324778A (zh) | 一种多支点飞机地面载荷确定方法 | |
CN109325271A (zh) | 一种运载火箭结构有限元模型参数化建模方法 | |
CN104036095A (zh) | 基于区域分解的耦合高精度复杂外形流场快速算法 | |
CN105468814A (zh) | 一种考虑全机气动干扰的直升机飞行特性计算方法 | |
CN116029219B (zh) | 一种飞行器气动热预测方法、装置、设备及存储介质 | |
CN102867097A (zh) | 计及静弹性变形影响的光固化快速成型风洞模型设计方法 | |
CN102262692A (zh) | 飞机翼面蒙皮亚音速颤振优化方法 | |
CN111563299A (zh) | 一种旋翼噪声确定方法及系统 | |
CN105046021A (zh) | 非定常气动力最小状态有理近似的非线性优化算法 | |
CN107529643A (zh) | 一种兼顾水力、强度性能的混流式转轮多工况优化方法 | |
CN107976908A (zh) | 一种飞行器耦合动稳定性特征分析方法 | |
CN104091003B (zh) | 一种基础运动时柔性壳结构大变形响应的有限元建模方法 | |
CN104573296B (zh) | 一种面向相似网格的高超声速流场初始化方法 | |
CN111551343B (zh) | 带栅格舵火箭子级全速域气动特性风洞试验设计方法 | |
CN103411628B (zh) | 一种mems陀螺仪随机漂移误差的处理方法 | |
CN106919730B (zh) | 一种采用风速衰减因子的风电场尾流确定方法 | |
CN105718619A (zh) | 一种基于有限元法的飞行器燃油质量特性确定方法 | |
CN108197381B (zh) | 基于寻优空间形态分析的参数辨识方法 | |
Ye et al. | A high-efficiency trim method for CFD numerical calculation of helicopter rotors | |
CN107563102A (zh) | 一种承载结构的传力骨架可视化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |