CN107526866A - 基于特征驱动的翼面结构拓扑优化方法 - Google Patents

基于特征驱动的翼面结构拓扑优化方法 Download PDF

Info

Publication number
CN107526866A
CN107526866A CN201710558819.3A CN201710558819A CN107526866A CN 107526866 A CN107526866 A CN 107526866A CN 201710558819 A CN201710558819 A CN 201710558819A CN 107526866 A CN107526866 A CN 107526866A
Authority
CN
China
Prior art keywords
mrow
mtd
msub
mtr
optimization
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
Application number
CN201710558819.3A
Other languages
English (en)
Other versions
CN107526866B (zh
Inventor
赵宇博
朱继宏
周莹
王杰
张卫红
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Northwestern Polytechnical University
Original Assignee
Northwestern Polytechnical University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Northwestern Polytechnical University filed Critical Northwestern Polytechnical University
Priority to CN201710558819.3A priority Critical patent/CN107526866B/zh
Publication of CN107526866A publication Critical patent/CN107526866A/zh
Application granted granted Critical
Publication of CN107526866B publication Critical patent/CN107526866B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/15Vehicle, aircraft or watercraft design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Turbine Rotor Nozzle Sealing (AREA)
  • Structures Of Non-Positive Displacement Pumps (AREA)

Abstract

本发明公开了一种基于特征驱动的翼面结构拓扑优化方法,用于解决现有翼面结构拓扑优化方法实用性差的技术问题。技术方案是首先确定翼面结构中与工程实际相符梁的特征形状,定义其角度,位置和尺寸为设计变量,采用固定网格方法,再利用水平集方法和布尔操作得到翼面特征布局初始结构,在优化的每一步迭代步中,根据优化变量结果对翼面特征布局进行更新,以一定材料体积为约束结构应变能最小为目标进行拓扑优化得到翼面加筋布局设计结果。用该方法对翼面结构进行拓扑优化,由作为特征的梁组合而成的布局优化结果可以直接指导并应用于工程实际,实用性好。

Description

基于特征驱动的翼面结构拓扑优化方法
技术领域
本发明涉及一种翼面结构拓扑优化方法,特别涉及一种基于特征驱动的翼面结构拓扑优化方法。
背景技术
文献“基于分级优化的飞机翼面结构布局综合技术研究.邓扬晨,张卫红,章怡宁.强度与环境,2005,32(1):27-35.”中提出采用分级优化与集成策略,通过将飞机翼面布局设计分为三个层次拓扑优化,尺寸优化和稳定性准则优化最终得到对飞机设计人员有一定指导性和参考价值的翼面结构布局设计。该方法以基于伪密度的翼面拓扑优化结果为基础,加入了尺寸优化和稳定性准则优化以刚度为目标结构轻量化为约束,优化最终得到能指导工程实际的设计。
文献所述方法基于伪密度的拓扑优化阶段只能得到翼面加筋材料分配大致情况,拓扑结果出现约束部位材料堆积而稍部无材料的情况,只有通过进一步三级优化才能得到对设计人员有指导性的翼面结构布局。在现代航空设计中,对于翼面结构,需要在保证翼面刚度等性能下,更高效的通过拓扑优化方法得到直接适用于工程实际的翼面布局结果。文献所述方法中基于伪密度的拓扑优化结果不能直接得到适用于工程实际的加筋布局设计。
发明内容
为了克服现有翼面结构拓扑优化方法实用性差的不足,本发明提供一种基于特征驱动的翼面结构拓扑优化方法。该方法首先确定翼面结构中与工程实际相符梁的特征形状,定义其角度,位置和尺寸为设计变量,采用固定网格方法,再利用水平集方法和布尔操作得到翼面特征布局初始结构,在优化的每一步迭代步中,根据优化变量结果对翼面特征布局进行更新,以一定材料体积为约束结构应变能最小为目标进行拓扑优化得到翼面加筋布局设计结果。用该方法对翼面结构进行拓扑优化,由作为特征的梁组合而成的布局优化结果可以直接指导并应用于工程实际,实用性好。
本发明解决其技术问题所采用的技术方案:一种基于特征驱动的翼面结构拓扑优化方法,其特点是包括以下步骤:
步骤一、建立翼面结构的表面框架蒙皮拓扑优化模型,定义表面框架蒙皮为非设计域,翼面内部为设计域。将翼面结构拓扑优化中的设计域离散为有限单元并作为固定网格。定义作为特征的梁的角度,位置和尺寸为设计变量,并给出初始值。
步骤二、利用水平集方法定义特征边界,通过布尔操作构造梁特征,将多根梁进行组合进行交并操作:
Φ=Φ1(x)∪Φ2(x)∪…∪Φn(x) (1)
式中,Φ-翼面加筋特征,n-特征数目,x-有限单元中点位置。设计域内的有限元单元被特征边界分为三类,即完全在特征区域外,完全在特征区域内和与边界相交,利用Heaviside函数对应表达三类有限元单元材料弹性模量数值:
式中,E0-实体材料弹性模量,Φ-翼面加筋特征,Δ-特征边界范围,α-极小值。对翼面根部固定位置施加固定约束,翼面表面施加压强载荷。
步骤三、翼面加筋布局优化以翼面结构总体应变能最小为目标,约束总的加筋体积及特征参数变量数值,具体满足公式:
式中,J(u,Φ)-结构总应变能,u-位移向量,m-特征变量类型的数目,s-特征设计变量组合即特征数目,布局角度,位置和尺寸变量,JD(Φ)-结构重力做功,JΓ(Φ)-表示外力做功,体积上限,-特征优化变量参数上限,s i-特征优化变量参数下限。
步骤四、计算翼面结构应变能函数对于梁特征尺寸,位置和角度等变量的灵敏度。
步骤五、根据上述求得的灵敏度进行优化,选取梯度优化算法,优化迭代得到结果。
本发明的有益效果是:该方法首先确定翼面结构中与工程实际相符梁的特征形状,定义其角度,位置和尺寸为设计变量,采用固定网格方法,再利用水平集方法和布尔操作得到翼面特征布局初始结构,在优化的每一步迭代步中,根据优化变量结果对翼面特征布局进行更新,以一定材料体积为约束结构应变能最小为目标进行拓扑优化得到翼面加筋布局设计结果。用该方法对翼面结构进行拓扑优化,由作为特征的梁组合而成的布局优化结果可以直接指导并应用于工程实际,实用性好。
本发明采用水平集方法和布尔操作,将符合工程实际形状的梁作为翼面拓扑优化结果的基本构型。在优化的过程中以梁的角度,位置和尺寸大小作为约束,给定约束上下界,用解析法求得该特征变量的灵敏度,同时引入材料用量约束,以翼面结构整体柔顺度为目标函数,进行结构拓扑优化得到设计结果。本发明能够得到适用于工程实际的翼面结构布局拓扑结果,使得原先只能为翼面结构设计人员提供概念指导的拓扑结果更加贴近工程实际,提高了拓扑优化的有效性和实用性。优化设计结果表明,在对于同一翼面结构进行拓扑优化时,背景技术方法中翼面结构拓扑结果和实际翼面结构存在很大差异。本发明基于特征驱动的拓扑优化结果和目标实际加筋形式更具有相似可比性,翼面结构拓扑构型对工程具有极大的参考价值并可以直接应用于工程。
下面结合附图和具体实施方式对本发明作详细说明。
附图说明
图1是本发明基于特征驱动的翼面结构拓扑优化方法设计结果示意图。
图2是背景技术翼面结构拓扑优化方法设计结果示意图。
具体实施方式
参照图1-2。本发明基于特征驱动的翼面结构拓扑优化方法具体步骤如下:
(a)建立翼面拓扑优化模型,翼面为一梯形表面且为变截面厚度,其厚度较长度比例为1:15,四周有框架,上下表面有外包蒙皮,厚度均为2mm。定义翼面结构内部为拓扑优化的设计域Ωd,框架及蒙皮为非设计域。将设计域Ωd离散为313440个三维单元作为固定网格。
(b)设置初始梁的数目,布局角度,位置和尺寸数值,利用水平集函数和布尔操作建立翼面加筋初始布局,数学模型表达为:
Φ=Φ1(x)∪Φ2(x)∪…∪Φn(x) (1)
式中,Φ-翼面加筋特征,n-特征数目,x-有限单元中点位置。设计域内的有限元单元被特征边界分为三类,即完全在特征区域外,完全在特征区域内和与边界相交,利用Heaviside函数对应表达三类有限元单元材料弹性模量数值:
式中,E0-实体材料弹性模量,Φ-翼面加筋特征,Δ-特征边界范围,α-极小值。采用矩形截面梁作为特征梁,翼面实体网格材料弹性模量取值为1MPa,泊松比μ=0.3,α取1×10-5,Δ取0.6。翼面结构长度方向根部有两处为固定边界约束,蒙皮上表面施加大小为0.001MPa的均布压强。
(c)定义精确变形优化问题;优化目标函数为结构总体应变能函数最小,约束材料体积分数小于满材料时的10%;
式中,J(u,Φ)-结构总应变能,u-位移向量,m-特征变量类型的数目,s-特征设计变量组合即特征数目,布局角度,位置和尺寸变量,JD(Φ)-结构重力做功,JΓ(Φ)-表示外力做功,-体积上限,-特征优化变量参数上限,si-特征优化变量参数下限。
(d)计算特征变量约束中,翼面结构的应变能对于设计域内布局角度,位置和尺寸的灵敏度以及计算目标函数的灵敏度。
(e)在优化过程中引入目标变形结构的应变能约束,根据上述求得的灵敏度选取梯度优化算法GCMMA(Globally Convergent Method of Moving Asymptotes)进行优化迭代,最终得到优化结果。
由图1、图2不同的拓扑优化结果对比可以看出,采用本发明方法,基于特征驱动拓扑优化后,翼面内部加筋布局结果明显,与实际翼梁布置方式具有可比性且更符合实际工程。背景技术基于伪密度方法的拓扑优化优化设计结果,其只能体现材料的大致分布,具体加筋布局还需要进一步通过分级优化得到;而在本发明方法在特征驱动方法基础上,其加筋布局结果明确对翼面结构布局有很强的指导性并在特征选取得当的情况下能够直接应用于工程实际。
本发明方法很好地解决了背景技术基于伪密度方法的翼面结构拓扑优化中结果加筋布局方向模糊,无法指导并应用于实际工程的问题。与背景技术拓扑优化结果相比,本发明方法的优化结果更符合实际加工制造工艺。

Claims (1)

1.一种基于特征驱动的翼面结构拓扑优化方法,其特征在于包括以下步骤:
步骤一、建立翼面结构的表面框架蒙皮拓扑优化模型,定义表面框架蒙皮为非设计域,翼面内部为设计域;将翼面结构拓扑优化中的设计域离散为有限单元并作为固定网格;定义作为特征的梁的角度,位置和尺寸为设计变量,并给出初始值;
步骤二、利用水平集方法定义特征边界,通过布尔操作构造梁特征,将多根梁进行组合进行交并操作:
Φ=Φ1(x)∪Φ2(x)∪…∪Φn(x) (1)
式中,Φ-翼面加筋特征,n-特征数目,x-有限单元中点位置;设计域内的有限元单元被特征边界分为三类,即完全在特征区域外,完全在特征区域内和与边界相交,利用Heaviside函数对应表达三类有限元单元材料弹性模量数值:
<mrow> <mi>E</mi> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>E</mi> <mn>0</mn> </msub> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>&amp;Phi;</mi> <mo>&amp;GreaterEqual;</mo> <mi>&amp;Delta;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mfrac> <mrow> <mn>3</mn> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> </mrow> </mrow> <mn>4</mn> </mfrac> <mrow> <mo>(</mo> <mfrac> <mi>&amp;Phi;</mi> <mi>&amp;Delta;</mi> </mfrac> <mo>-</mo> <mfrac> <msup> <mi>&amp;omega;</mi> <mn>3</mn> </msup> <mrow> <mn>3</mn> <msup> <mi>&amp;Delta;</mi> <mn>3</mn> </msup> </mrow> </mfrac> <mo>)</mo> </mrow> <msub> <mi>E</mi> <mn>0</mn> </msub> <mo>+</mo> <mfrac> <mrow> <mn>1</mn> <mo>+</mo> <mi>&amp;alpha;</mi> </mrow> <mn>2</mn> </mfrac> <msub> <mi>E</mi> <mn>0</mn> </msub> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mo>-</mo> <mi>&amp;Delta;</mi> <mo>&amp;le;</mo> <mi>&amp;Phi;</mi> <mo>&lt;</mo> <mi>&amp;Delta;</mi> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>&amp;alpha;E</mi> <mn>0</mn> </msub> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mi>&amp;Phi;</mi> <mo>&lt;</mo> <mo>-</mo> <mi>&amp;Delta;</mi> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> </mrow>
式中,E0-实体材料弹性模量,Φ-翼面加筋特征,Δ-特征边界范围,α-极小值;对翼面根部固定位置施加固定约束,翼面表面施加压强载荷;
步骤三、翼面加筋布局优化以翼面结构总体应变能最小为目标,约束总的加筋体积及特征参数变量数值,具体满足公式:
<mrow> <mtable> <mtr> <mtd> <mrow> <mi>f</mi> <mi>i</mi> <mi>n</mi> <mi>d</mi> </mrow> </mtd> <mtd> <mrow> <mi>s</mi> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>s</mi> <mn>1</mn> </msub> <mo>,</mo> <msub> <mi>s</mi> <mn>2</mn> </msub> <mo>,</mo> <mo>...</mo> <msub> <mi>s</mi> <mi>m</mi> </msub> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mi>min</mi> </mtd> <mtd> <mrow> <mi>J</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>&amp;Phi;</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <mi>s</mi> <mo>.</mo> <mi>t</mi> <mo>.</mo> </mrow> </mtd> <mtd> <mrow> <mi>J</mi> <mrow> <mo>(</mo> <mi>u</mi> <mo>,</mo> <mi>&amp;Phi;</mi> <mo>)</mo> </mrow> <mo>=</mo> <msub> <mi>J</mi> <mi>D</mi> </msub> <mrow> <mo>(</mo> <mi>&amp;Phi;</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>J</mi> <mi>&amp;Gamma;</mi> </msub> <mrow> <mo>(</mo> <mi>&amp;Phi;</mi> <mo>)</mo> </mrow> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mrow> <mi>V</mi> <mo>&amp;le;</mo> <mover> <mi>V</mi> <mo>&amp;OverBar;</mo> </mover> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow></mrow> </mtd> <mtd> <mrow> <msub> <munder> <mi>s</mi> <mo>&amp;OverBar;</mo> </munder> <mi>i</mi> </msub> <mo>&amp;le;</mo> <msub> <mi>s</mi> <mi>i</mi> </msub> <mo>&amp;le;</mo> <msub> <mover> <mi>s</mi> <mo>&amp;OverBar;</mo> </mover> <mi>i</mi> </msub> </mrow> </mtd> </mtr> </mtable> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>3</mn> <mo>)</mo> </mrow> </mrow>
式中,J(u,Φ)-结构总应变能,u-位移向量,m-特征变量类型的数目,s-特征设计变量组合即特征数目,布局角度,位置和尺寸变量,JD(Φ)-结构重力做功,JΓ(Φ)-表示外力做功,-体积上限,-特征优化变量参数上限,s i-特征优化变量参数下限;
步骤四、计算翼面结构应变能函数对于梁特征尺寸,位置和角度等变量的灵敏度;
步骤五、根据上述求得的灵敏度进行优化,选取梯度优化算法,优化迭代得到结果。
CN201710558819.3A 2017-07-11 2017-07-11 基于特征驱动的翼面结构拓扑优化方法 Active CN107526866B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710558819.3A CN107526866B (zh) 2017-07-11 2017-07-11 基于特征驱动的翼面结构拓扑优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710558819.3A CN107526866B (zh) 2017-07-11 2017-07-11 基于特征驱动的翼面结构拓扑优化方法

Publications (2)

Publication Number Publication Date
CN107526866A true CN107526866A (zh) 2017-12-29
CN107526866B CN107526866B (zh) 2020-05-01

Family

ID=60748946

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710558819.3A Active CN107526866B (zh) 2017-07-11 2017-07-11 基于特征驱动的翼面结构拓扑优化方法

Country Status (1)

Country Link
CN (1) CN107526866B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110728084A (zh) * 2019-09-16 2020-01-24 中国第一汽车股份有限公司 空心薄壁铸铝件正向设计方法
CN112182760A (zh) * 2020-09-27 2021-01-05 西北工业大学 基于指节式刚柔耦合的变形翼后缘设计方法
CN113468697A (zh) * 2021-09-03 2021-10-01 江苏中信博新能源科技股份有限公司 一种光伏板钢边框结构的正向设计方法
CN114291250A (zh) * 2021-12-20 2022-04-08 北京机电工程研究所 一种剪切变后掠翼面及其设计方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080300831A1 (en) * 2006-12-19 2008-12-04 Board Of Governors For Higher Education, State Of Rhode Island And Providence System and method for finite element based on topology optimization
US20110270587A1 (en) * 2008-09-11 2011-11-03 Kyoto University Structural optimization system, structural optimization methodlogy, and structural optimization program
WO2012078698A1 (en) * 2010-12-07 2012-06-14 Engineering Technology Associates, Inc. Design optimization system
US20140207428A1 (en) * 2013-01-23 2014-07-24 Toyota Motor Engineering & Manufacturing North America, Inc. Methods for Orienting Material Physical Properties Using Constraint Transformation and Isoparametric Shape Functions
CN104765922A (zh) * 2015-04-13 2015-07-08 西北工业大学 基于保形约束的悬臂梁结构拓扑优化设计方法
CN105512385A (zh) * 2015-12-04 2016-04-20 吉林大学 采用序列幂函数插值方法实现结构多材料拓扑优化的方法
CN106547972A (zh) * 2016-11-04 2017-03-29 中国科学院长春光学精密机械与物理研究所 基于参数化水平集方法的流体管道拓扑优化设计方法
CN106874572A (zh) * 2017-01-21 2017-06-20 西安交通大学 一种考虑串油特性的飞机油箱承载结构轻量化设计方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080300831A1 (en) * 2006-12-19 2008-12-04 Board Of Governors For Higher Education, State Of Rhode Island And Providence System and method for finite element based on topology optimization
US20110270587A1 (en) * 2008-09-11 2011-11-03 Kyoto University Structural optimization system, structural optimization methodlogy, and structural optimization program
WO2012078698A1 (en) * 2010-12-07 2012-06-14 Engineering Technology Associates, Inc. Design optimization system
US20140207428A1 (en) * 2013-01-23 2014-07-24 Toyota Motor Engineering & Manufacturing North America, Inc. Methods for Orienting Material Physical Properties Using Constraint Transformation and Isoparametric Shape Functions
CN104765922A (zh) * 2015-04-13 2015-07-08 西北工业大学 基于保形约束的悬臂梁结构拓扑优化设计方法
CN105512385A (zh) * 2015-12-04 2016-04-20 吉林大学 采用序列幂函数插值方法实现结构多材料拓扑优化的方法
CN106547972A (zh) * 2016-11-04 2017-03-29 中国科学院长春光学精密机械与物理研究所 基于参数化水平集方法的流体管道拓扑优化设计方法
CN106874572A (zh) * 2017-01-21 2017-06-20 西安交通大学 一种考虑串油特性的飞机油箱承载结构轻量化设计方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
ATSUSHI KAWAMOTO等: "Heaviside projection based topology optimization by a PDE-filtered scalar function", 《STRUCT MULTIDISC OPTIM (2011)》 *
邓扬晨等: "基于分级优化的飞机翼面结构布局综合技术研究", 《强度与环境》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110728084A (zh) * 2019-09-16 2020-01-24 中国第一汽车股份有限公司 空心薄壁铸铝件正向设计方法
CN112182760A (zh) * 2020-09-27 2021-01-05 西北工业大学 基于指节式刚柔耦合的变形翼后缘设计方法
CN113468697A (zh) * 2021-09-03 2021-10-01 江苏中信博新能源科技股份有限公司 一种光伏板钢边框结构的正向设计方法
CN113468697B (zh) * 2021-09-03 2021-11-30 江苏中信博新能源科技股份有限公司 一种光伏板钢边框结构的正向设计方法
CN114291250A (zh) * 2021-12-20 2022-04-08 北京机电工程研究所 一种剪切变后掠翼面及其设计方法
CN114291250B (zh) * 2021-12-20 2023-11-03 北京机电工程研究所 一种剪切变后掠翼面及其设计方法

Also Published As

Publication number Publication date
CN107526866B (zh) 2020-05-01

Similar Documents

Publication Publication Date Title
CN107526866A (zh) 基于特征驱动的翼面结构拓扑优化方法
Li et al. Topology optimization for functionally graded cellular composites with metamaterials by level sets
Werter et al. A novel dynamic aeroelastic framework for aeroelastic tailoring and structural optimisation
Lian et al. Membrane wing aerodynamics for micro air vehicles
CN110569519B (zh) 考虑非设计域的三维连续体结构动静力学性能拓扑优化设计方法
CN105528503A (zh) 一种基于结构分解的大型构件动态优化设计方法
Boston et al. Spanwise wing morphing using multistable cellular metastructures
Ciampa et al. Global local structural optimization of transportation aircraft wings
CN105447269A (zh) 一种飞机在定常过载状态下的非结构质量配平方法
CN106777693A (zh) 一种基于catia数模的实时优化设计方法
TSUSHIMA et al. Recent researches on morphing aircraft technologies in Japan and other countries
CN106202693A (zh) 一种基于参数化建模的加筋壁板结构抗振动疲劳优化方法
Mitrotta et al. Multi-fidelity design of an aeroelastically tailored composite wing for dynamic wind-tunnel testing
CN107169174A (zh) 最优焦距下伞状天线机电集成优化设计方法
Ito et al. Efficient computational fluid dynamics evaluation of small-device locations with automatic local remeshing
DeBlois et al. Multi-fidelity multidisciplinary design optimization of metallic and composite regional and business jets
Kawabe et al. Application of Macroscopic Structures on Dragonfly Wings to an Aircraft Design Approach
Chen et al. Additive manufacturing of a continuum topology-optimized palletizing manipulator arm
Rhodes et al. Structural optimization of a morphing shock control bump
Boston et al. Multistable honeycomb architecture for spanwise wing morphing
Bhattacharyya et al. Topology optimization of a bi-stable airfoil using nonlinear elasticity
Vineeth et al. Influence of aspect ratio on the wake dynamics of a pitching wing
Mason et al. Design and Analysis of a Tool for Automated Fiber Placement of Composite Wind Tunnel Blades
Yang On the sensitivity and efficiency of aerodynamic shape optimisation
Gillebaart et al. Geometrically consistent static aeroelastic simulation using isogeometric analysis

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