CN112001004A - 一种解析中高频振动结构能量密度场的nurbs等几何分析方法 - Google Patents
一种解析中高频振动结构能量密度场的nurbs等几何分析方法 Download PDFInfo
- Publication number
- CN112001004A CN112001004A CN202010850855.9A CN202010850855A CN112001004A CN 112001004 A CN112001004 A CN 112001004A CN 202010850855 A CN202010850855 A CN 202010850855A CN 112001004 A CN112001004 A CN 112001004A
- Authority
- CN
- China
- Prior art keywords
- nurbs
- energy density
- energy
- unit
- representing
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
-
- 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)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Computing Systems (AREA)
- Operations Research (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Feedback Control In General (AREA)
Abstract
一种解析中高频振动结构能量密度场的NURBS等几何分析方法,先建立设计域的等几何分析模型,设计域利用NURBS曲线、NURBS曲面以及NURBS体进行描述,然后构建稳态下能量密度控制方程,再构建单元能量密度控制方程矩阵形式,然后计算基于NURBS基函数的能量密度控制方程,最后进行适应性处理,将在等几何分析框架下改造得到的能量密度控制方程代入结构的动力学分析计算中,通过求解计算得到高频激励下结构中任意位置处的能量密度响应;本发明在等几何分析的框架下,从稳态下的能量密度控制方程出发,以经过时间和空间平均的能量密度作为研究对象,采用能量有限元方法对机械结构进行动力学分析,不仅使得分析时的网格划分摆脱了振动波长的桎梏,还大大提高了精度。
Description
技术领域
本发明属于结构的动力学分析技术领域,具体涉及一种解析中高频振动结构能量密度场的NURBS等几何分析方法。
技术背景
复杂工况下尤其是高频激励下的结构动力学分析一直以来都是工程中备受关注的问题,以经典有限元为代表的传统方法在解决高频激励下复杂结构的动力学分析问题时的表现往往不尽人意;以航天器为例,由于航天器在飞行中受到冲击、脉动、噪声等复杂激励作用,其整流罩外气动噪声引起的振动频率可达10000Hz以上,对于这样的高频下结构动力学分析问题,结构的模态密集,且相对于整个系统的尺度,结构以短波振动为主,采用传统动力学分析方法进行分析时,为保证分析的准确性必须划分足够细小的网格,进而导致运算量增加;此外有限元分析模型只是几何模型的近似而非精确表示,高频时系统的动态响应对结构参数高度敏感,利用经典有限元方法建模时存在的细节忽略会对分析结果的准确性造成巨大影响;其几何模型与分析模型的分裂导致分析计算时在划分网格近似表达的步骤中耗费大量时间,尤其是面对类似于航天器的复杂结构时,有限元方法为了满足仿真精度必须划分非常细密的网格,导致计算的成本过大,造成分析精度与求解效率之间的矛盾,且容易产生较大误差,因此以经典有限元方法为代表的传统方法难以有效地进行高频下复杂结构的动力学分析。
发明内容
为了克服上述现有技术的缺点,本发明的目的在于提供一种解析中高频振动结构能量密度场的NURBS等几何分析方法,不仅使得分析时的网格划分摆脱了振动波长的桎梏,还大大提高了精度,为复杂工程结构的中高频响应提供计算效率更高、分析结果更加准确的动力学分析方法。
为了达到上述目标,本发明采取的技术方案为:
一种解析中高频振动结构能量密度场的NURBS等几何分析方法,包括以下步骤:
1)建立设计域的等几何分析模型:
设计域利用NURBS曲线、NURBS曲面以及NURBS体进行描述,二维结构采用NURBS曲面描述设计域形状,定义局部空间坐标系为参数空间坐标系为ξ=(ξ,η),物理空间坐标系为x=(x,y);在参数空间中分别定义ξ和η两方向上的节点向量Ξ1={ξ1,...,ξn+p+1}与Ξ2={η1,...,ηm+q+1},相应方向上的NURBS基函数阶数分别为p和q,由此将设计域离散为若干NURBS单元组成的物理网格;定义Pi,j为控制点集合,参数空间中ξ和η两方向上控制点数分别为n和m,控制点总数为nnp,nnp=m×n,由此形成了控制点组成的控制网格;
用在ξ方向p次、η方向q次的NURBS曲面描述设计域,其表达式为:
2)构建稳态下能量密度控制方程:
定义Γ表示设计域Ω的边界,从设计域中弹性波的能量流平衡原理出发得到能量密度控制微分方程:
其中,η表示阻尼系数,ω为激励频率,ε表示经过时间与空间平均的能量密度,πin表示输入功率,cg表示弹性波的群速度,表达式如下:
上式中,Υ表示弹性模量,h表示设计域厚度,λ表示泊松比,ρ表示材料密度;
3)构建单元能量密度控制方程矩阵形式:
以NURBS基函数作为形函数对单元能量密度ε近似表示:
定义一个时间与空间平均的能量集中度I和能量流q分别为:
结合式(7)-式(8),对式(6)采用一种标准的Gauss积分方法,得到单元能量密度控制方程矩阵形式:
[Ke]{Ee}={Fe}-{Qe} (9)
其中,[Ke]表示单元能量矩阵,{Ee}表示单元节点能量密度向量,{Fe}表示单元输入功率向量,{Qe}表示单元边界上的能量流向量;
4)计算基于NURBS基函数的能量密度控制方程:
4.1)NURBS基函数:
定义在参数空间上的NURBS基函数表达式为:
NURBS基函数的一阶导通过求导法则计算得到:
4.2)通过两次映射,将物理坐标系下的NURBS单元映射到局部坐标系下的等参单元:
局部空间与参数空间之间的映射关系:
其中,(ξi+1,ξi)和(ηj+1,ηj)分别表示参数空间中NURBS单元在ξ和η方向的坐标范围;
参数空间与物理空间之间的映射关系:
dxdy=|Jξ|dξdη (17)
在等几何分析的框架下,几何矩阵B表达式:
其中,Ra,x表示NURBS基函数的空间导数,其计算公式如下:
4.3)计算和组装基于NURBS基函数的单元矩阵:
单元输入功率向量表达式为:
对计算得到的单元能量矩阵与单元输入功率向量采用与经典有限元相同的直接组装方法组装成全局形式[K]、{F};
结构变化处产生不连续边界,能量密度在不连续边界上不连续;在不连续边界处新增控制点,假设两相邻单元e和单元e+1的相交边界不连续,和和分别表示边界上的两单元同位控制点处能量流,在单元边界上的能量流向量用单元边界矩阵[Jee +1]和单元控制点能量密度向量{Ee}表示:
按照式(24)-式(25)计算整个设计域的能量流,得到能量流向量的全局形式表达式:
{Q}=[J]{E} (26)
其中,[J]为单元边界矩阵的全局形式,τee+1和ree分别表示能量传递系数和反射系数;
按照上述步骤计算并组装,得到基于NURBS基函数的能量密度控制方程矩阵形式:
([K]+[J]){E}={F} (27)
5)适应性处理:
将在等几何分析框架下改造得到的能量密度控制方程代入结构的动力学分析计算中,通过求解计算得到高频激励下结构中任意位置处的能量密度响应。
本发明的有益效果为:
由于本发明方法是在等几何分析的框架下进行,将NURBS基函数作为形函数,可以直接利用几何信息进行分析计算,为复杂几何体分析提供更加精确的模型,无需繁琐的网格操作,这使得本发明方法的求解精度、稳定性和效率远超传统的有限元分析;同时NURBS基函数本身具有的非负性、光滑性以及在单元界面上可保持高阶连续性等特点是传统高阶有限元基函数不具备的独特优势;
本发明从稳态下的能量密度控制方程出发,以经过时间和空间平均的能量密度作为研究变量,具有能够得到结构内部任意位置响应信息的优势,可以在更宽广的频率上进行动力学分析;本发明方法为复杂机械系统的动力学分析提供了一种高效有力的工具,同时在航空航天等高端设备的振动分析控制等领域也拥有广泛的应用前景。
附图说明
图1为本发明实施例的边界条件示意图。
图2为本发明实施例中等几何分析模型示意图。
图3为本发明实施例中离散网格与NURBS基函数示意图。
图4为本发明实施例中结构能量平衡示意图。
图5为本发明实施例中单元不连续边界示意图。
图6为本发明实施例中结构能量密度响应分布示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步说明,实施例设计域为四分之一圆环薄板结构,薄板尺寸参数为外圆半径R1=3m,内圆半径R2=1m,板厚h=0.01m,定义薄板的内圆边缘固定,外圆边缘受均匀分布的单位简谐载荷,如图1所示。
一种解析中高频振动结构能量密度场的NURBS等几何分析方法,包括以下步骤:
1)建立设计域的等几何分析模型:
设计域利用NURBS曲线、NURBS曲面以及NURBS体进行描述,本实施例中设计域二维结构由一个二次NURBS曲面描述,如图2所示,定义局部空间坐标系为参数空间坐标系为ξ=(ξ,η),物理空间坐标系为x=(x,y);在参数空间中分别定义ξ和η两方向上的节点向量Ξ1={0,0,1/3,2/3,1,1}与Ξ2={0,0,1/3,2/3,1,1},相应方向上的NURBS基函数阶数分别为p=1和q=1,由此将设计域离散为若干NURBS单元组成的物理网格;Pi,j为控制点集合,参数空间中ξ和η两方向上控制点数分别为n=4和m=4,控制点总数为nnp=16,由此形成了控制点组成的控制网格;
如图3所示,用在ξ方向p次、η方向q次的NURBS曲面来描述设计域,其表达式为:
2)构建稳态下能量密度控制方程:
如图4所示,定义Γ表示设计域Ω的边界,从薄板中弹性波的能量平衡公式出发得到能量密度控制微分方程:
其中,η表示阻尼系数,取0.08,ω为激励频率,取1000Hz。ε表示经过时间与空间平均的能量密度,πin表示输入功率,为1W,cg表示弹性波的群速度,表达式如下:
上式中,Υ表示弹性模量,取Υ=7.1×1010N·m2,h表示设计域厚度,λ表示泊松比,取0.35,ρ表示材料密度,取2700kg/m2;
3)构建单元能量密度控制方程矩阵形式:
以NURBS基函数作为形函数对单元能量密度ε近似表示:
定义一个时间与空间平均的能量集中度I和能量流q分别为:
结合式(7)-式(8),对式(6)采用一种标准的Gauss积分方法,得到单元能量密度控制方程矩阵形式:
其中,[Ke]表示单元能量矩阵,{Ee}表示单元节点能量密度向量,{Fe}表示单元输入功率向量,{Qe}表示单元边界上的能量流向量;
4)计算基于NURBS基函数的能量密度控制方程:
4.1)NURBS基函数:
定义在参数空间上的NURBS基函数表达式为:
NURBS基函数的一阶导通过求导法则计算得到:
4.2)通过两次映射,将物理坐标系下的NURBS单元映射到局部坐标系下的等参单元:
局部空间与参数空间之间的映射关系:
其中,(ξi+1,ξi)和(ηj+1,ηj)分别表示参数空间中NURBS单元在ξ和η方向的坐标范围;
参数空间与物理空间之间的映射关系:
dxdy=|Jξ|dξdη (17)
在等几何分析的框架下,几何矩阵B表达式:
其中,Ra,x表示非零NURBS基函数的空间导数,其计算公式如下:
4.3)计算和组装基于NURBS基函数的矩阵:
单元输入功率向量表达式为:
对计算得到的单元能量矩阵与单元输入功率向量采用与经典有限元相同的直接组装方法组装成全局形式[K]、{F};
结构变化处产生不连续边界,能量密度在不连续边界上不连续;在不连续边界处新增控制点,如图5所示,两相邻单元1和单元2的相交边界不连续,Q2和Q5、Q3和Q8分别表示边界上的两单元同位节点处能量流,在单元边界上的能量流向量用单元边界矩阵[J12]和单元节点能量密度向量{Ee}表示:
按照式(24)-式(25)计算整个设计域的能量流,得到能量流向量的全局形式表达式:
{Q}=[J]{E} (26)
其中,[J]为单元边界矩阵的全局形式,τ12和r11分别表示能量传递系数和反射系数,具体计算见已有学术研究,此处不再赘述;
按照上述步骤计算并组装,得到基于NURBS基函数的能量密度控制方程矩阵形式:
([K]+[J]){E}={F} (27)
8)适应性处理:
将在等几何分析框架下改造得到的能量密度控制方程代入结构的动力学分析计算中,通过求解计算可以得到高频激励下结构中任意位置处的能量密度响应,如图6所示。
为适应不同分析需求,使用时并不局限于所述的材料相同的薄板结构,设计者可以将本发明方法运用到多相结构的动力学分析运算中,且对于更加复杂的结构模型同样能够适应并得到准确的分析结果。
Claims (1)
1.一种解析中高频振动结构能量密度场的NURBS等几何分析方法,其特征在于,包括以下步骤:
1)建立设计域的等几何分析模型:
设计域利用NURBS曲线、NURBS曲面以及NURBS体进行描述,二维结构采用NURBS曲面描述设计域形状,定义局部空间坐标系为参数空间坐标系为ξ=(ξ,η),物理空间坐标系为x=(x,y);在参数空间中分别定义ξ和η两方向上的节点向量Ξ1={ξ1,...,ξn+p+1}与Ξ2={η1,...,ηm+q+1},相应方向上的NURBS基函数阶数分别为p和q,由此将设计域离散为若干NURBS单元组成的物理网格;定义Pi,j为控制点集合,参数空间中ξ和η两方向上控制点数分别为n和m,控制点总数为nnp,nnp=m×n,由此形成了控制点组成的控制网格;
用在ξ方向p次、η方向q次的NURBS曲面描述设计域,其表达式为:
2)构建稳态下能量密度控制方程:
定义Γ表示设计域Ω的边界,从设计域中弹性波的能量流平衡原理出发得到能量密度控制微分方程:
其中,η表示阻尼系数,ω为激励频率,ε表示经过时间与空间平均的能量密度,πin表示输入功率,cg表示弹性波的群速度,表达式如下:
上式中,Υ表示弹性模量,h表示设计域厚度,λ表示泊松比,ρ表示材料密度;
3)构建单元能量密度控制方程矩阵形式:
以NURBS基函数作为形函数对单元能量密度ε近似表示:
定义一个时间与空间平均的能量集中度I和能量流q分别为:
结合式(7)-式(8),对式(6)采用一种标准的Gauss积分方法,得到单元能量密度控制方程矩阵形式:
[Ke]{Ee}={Fe}-{Qe} (9)
其中,[Ke]表示单元能量矩阵,{Ee}表示单元节点能量密度向量,{Fe}表示单元输入功率向量,{Qe}表示单元边界上的能量流向量;
4)计算基于NURBS基函数的能量密度控制方程:
4.1)NURBS基函数:
定义在参数空间上的NURBS基函数表达式为:
NURBS基函数的一阶导通过求导法则计算得到:
4.2)通过两次映射,将物理坐标系下的NURBS单元映射到局部坐标系下的等参单元:
局部空间与参数空间之间的映射关系:
其中,(ξi+1,ξi)和(ηj+1,ηj)分别表示参数空间中NURBS单元在ξ和η方向的坐标范围;
参数空间与物理空间之间的映射关系:
dxdy=|Jξ|dξdη (17)
在等几何分析的框架下,几何矩阵B表达式:
其中,Ra,x表示NURBS基函数的空间导数,其计算公式如下:
4.3)计算和组装基于NURBS基函数的单元矩阵:
单元输入功率向量表达式为:
对计算得到的单元能量矩阵与单元输入功率向量采用与经典有限元相同的直接组装方法组装成全局形式[K]、{F};
结构变化处产生不连续边界,能量密度在不连续边界上不连续;在不连续边界处新增控制点,假设两相邻单元e和单元e+1的相交边界不连续,和和分别表示边界上的两单元同位控制点处能量流,在单元边界上的能量流向量用单元边界矩阵[Jee+1]和单元控制点能量密度向量{Ee}表示:
按照式(24)-式(25)计算整个设计域的能量流,得到能量流向量的全局形式表达式:
{Q}=[J]{E} (26)其中,[J]为单元边界矩阵的全局形式,τee+1和ree分别表示能量传递系数和反射系数;
按照上述步骤计算并组装,得到基于NURBS基函数的能量密度控制方程矩阵形式:
([K]+[J]){E}={F} (27)
5)适应性处理:
将在等几何分析框架下改造得到的能量密度控制方程代入结构的动力学分析计算中,通过求解计算得到高频激励下结构中任意位置处的能量密度响应。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010850855.9A CN112001004B (zh) | 2020-08-21 | 2020-08-21 | 一种解析中高频振动结构能量密度场的nurbs等几何分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010850855.9A CN112001004B (zh) | 2020-08-21 | 2020-08-21 | 一种解析中高频振动结构能量密度场的nurbs等几何分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112001004A true CN112001004A (zh) | 2020-11-27 |
CN112001004B CN112001004B (zh) | 2022-12-09 |
Family
ID=73472318
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010850855.9A Active CN112001004B (zh) | 2020-08-21 | 2020-08-21 | 一种解析中高频振动结构能量密度场的nurbs等几何分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112001004B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112765884A (zh) * | 2021-01-18 | 2021-05-07 | 河海大学 | 一种基于扩展等几何分析法的声子晶体材料结构设计方法 |
CN113343462A (zh) * | 2021-06-07 | 2021-09-03 | 西安交通大学 | 基于高阶等几何的多油腔动静压滑动轴承油膜特性仿真方法 |
CN114491769A (zh) * | 2022-02-17 | 2022-05-13 | 河海大学 | 一种基于等几何分析法自由曲面结构一体化形态创构方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108763658A (zh) * | 2018-05-07 | 2018-11-06 | 长安大学 | 基于等几何方法的组合薄壁结构固有频率设计方法 |
CN109117504A (zh) * | 2018-07-09 | 2019-01-01 | 哈尔滨工程大学 | 一种双向功能梯度曲壳振动分析方法 |
CN109670200A (zh) * | 2018-11-13 | 2019-04-23 | 华中科技大学 | 一种等几何材料密度场结构拓扑优化方法 |
CN109684723A (zh) * | 2018-12-24 | 2019-04-26 | 哈尔滨工程大学 | 一种二维结构内部声学性能分析方法 |
-
2020
- 2020-08-21 CN CN202010850855.9A patent/CN112001004B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108763658A (zh) * | 2018-05-07 | 2018-11-06 | 长安大学 | 基于等几何方法的组合薄壁结构固有频率设计方法 |
CN109117504A (zh) * | 2018-07-09 | 2019-01-01 | 哈尔滨工程大学 | 一种双向功能梯度曲壳振动分析方法 |
CN109670200A (zh) * | 2018-11-13 | 2019-04-23 | 华中科技大学 | 一种等几何材料密度场结构拓扑优化方法 |
CN109684723A (zh) * | 2018-12-24 | 2019-04-26 | 哈尔滨工程大学 | 一种二维结构内部声学性能分析方法 |
Non-Patent Citations (3)
Title |
---|
张勇等: "基于等几何分析的比例边界有限元方法", 《计算力学学报》 * |
张勇等: "等几何分析方法中重控制点问题的研究与应用", 《工程力学》 * |
荣吉利等: "基于HB样条的平面大变形结构动力学等几何分析", 《北京理工大学学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112765884A (zh) * | 2021-01-18 | 2021-05-07 | 河海大学 | 一种基于扩展等几何分析法的声子晶体材料结构设计方法 |
CN113343462A (zh) * | 2021-06-07 | 2021-09-03 | 西安交通大学 | 基于高阶等几何的多油腔动静压滑动轴承油膜特性仿真方法 |
CN114491769A (zh) * | 2022-02-17 | 2022-05-13 | 河海大学 | 一种基于等几何分析法自由曲面结构一体化形态创构方法 |
CN114491769B (zh) * | 2022-02-17 | 2022-12-09 | 河海大学 | 一种基于等几何分析法自由曲面结构一体化形态创构方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112001004B (zh) | 2022-12-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112001004B (zh) | 一种解析中高频振动结构能量密度场的nurbs等几何分析方法 | |
Glaz et al. | Reduced-order nonlinear unsteady aerodynamic modeling using a surrogate-based recurrence framework | |
Reuther et al. | Aerodynamic shape optimization of supersonic aircraft configurations via an adjoint formulation on distributed memory parallel computers | |
CN109933876B (zh) | 一种基于广义气动力的非定常气动力降阶方法 | |
CN109726465B (zh) | 基于非结构曲边网格的三维无粘低速绕流的数值模拟方法 | |
Gang et al. | Mesh deformation on 3D complex configurations using multistep radial basis functions interpolation | |
Rumpfkeil et al. | Efficient hessian calculations using automatic differentiation and the adjoint method with applications | |
Sadeghi et al. | Application of three-dimensional interfaces for data transfer in aeroelastic computations | |
Morris et al. | Domain-element method for aerodynamic shape optimization applied to modern transport wing | |
Hall | Deforming grid variational principle for unsteady small disturbance flows in cascades | |
Amoignon et al. | Study of parameterizations in the project CEDESA | |
Jonsson et al. | Development of flutter constraints for high-fidelity aerostructural optimization | |
Da Ronch et al. | Modeling of unsteady aerodynamic loads | |
Zhu et al. | Design of an RBF surrogate model for low Reynolds number airfoil based on transfer learning | |
Apponsah et al. | Aerodynamic shape optimization for unsteady flows: some benchmark problems | |
CN116187074A (zh) | 基于等几何的各向异性周期结构材料多尺度拓扑优化方法 | |
Morris et al. | Development of generic CFD-based aerodynamic optimisation tools for helicopter rotor blades | |
Mackle et al. | Efficient and Flexible Methodology for the Aerodynamic Shape Optimization of Hypersonic Vehicle Concepts in a High-Dimensional Design Space | |
Florea et al. | Eigenmode analysis of unsteady flows about airfoils | |
Burgee et al. | Parallel multipoint variable-complexity approximations for multidisciplinary optimization | |
Chen | Numerical calculations on the unsteady aerodynamic force of the tilt-rotor aircraft in conversion mode | |
Athavale et al. | Application of an unstructured grid solution methodology to turbomachinery flows | |
Mian et al. | RBF interpolation with improved data reduction algorithm—A meshfree method for fluid-structure coupling and mesh deformation | |
Poole et al. | Objective Function and Constraints for Robust Transonic Aerofoil Optimization | |
Rumpfkeil et al. | Efficient hessian calculations using automatic differentiation and the adjoint method |
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 |