CN112464372A - 一种弹性机翼副翼舵面效率的设计敏度工程数值方法 - Google Patents
一种弹性机翼副翼舵面效率的设计敏度工程数值方法 Download PDFInfo
- Publication number
- CN112464372A CN112464372A CN202011338381.6A CN202011338381A CN112464372A CN 112464372 A CN112464372 A CN 112464372A CN 202011338381 A CN202011338381 A CN 202011338381A CN 112464372 A CN112464372 A CN 112464372A
- Authority
- CN
- China
- Prior art keywords
- aileron
- wing
- control surface
- design
- sensitivity
- 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
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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)
- 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)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
- Toys (AREA)
Abstract
一种弹性机翼副翼舵面效率的设计敏度工程数值方法,包括如下步骤:建立气动面元模型和结构有限元模型,并计算气动面元模型和结构有限元模型间数据传递矩阵;计算刚性机翼单位副翼偏转角产生的滚转速率;计算弹性机翼单位副翼偏转角产生的驱动力矩和单位滚转速率产生的阻尼力矩;进行弹性机翼副翼舵面效率及其设计敏度的计算。本发明将弹性机翼单位副翼偏转角产生的滚转速率表示为单位副翼偏转角产生的驱动力矩和单位滚转速率产生的阻尼力矩的组合的形式,解决了滚转速率形式定义的副翼舵面效率的设计敏度难以解析求解的问题,在此基础上,通过直接法或伴随法对副翼舵面效率设计敏度的工程计算,从而辅助后续飞机结构的精细数值优化设计。
Description
技术领域
本发明适用于飞机结构的数值优化设计领域,具体涉及一种弹性机翼副翼舵面效率关于结构设计变量敏度的工程数值方法。
背景技术
施加于弹性机翼上的气动载荷会改变机翼的变形,并反过来改变机翼的气流特征,从而产生气动弹性现象并影响最终稳定变形状态下机翼的气动载荷。副翼用于控制飞机的飞行状态,随着飞行速度的增大,副翼舵面效率会不断降低,当达到反效速度时,副翼偏转将对飞行状态不产生任何影响。副翼舵面效率的大小反映了飞机的机动性能,常作为约束条件用于飞机结构的数值优化设计,弹性机翼副翼舵面效率的设计敏度可为梯度优化算法提供必要的指导。在实际工程中一般有两种副翼舵面效率的定义形式:机翼根部固支情况通过弹性机翼和刚性机翼的翼根弯矩之比进行定义,机翼稳定滚转情况通过弹性机翼和刚性机翼单位副翼偏转角产生的滚转速率之比进行定义。前者是一种较简单的形式,未能考虑副翼对飞机滚转运动的实际作用效果,因此,本发明中副翼舵面效率采用第二种定义形式。如文献“Wright JR,Cooper JE.Introduction to AircraftAeroelasticity andLoads.2nd ed.New York:Wiley,2015”所述,目前常用的副翼舵面效率计算方法是依据机翼气动弹性原理,利用简单的工程算法公式以及粗定量技术参数进行求解。其计算结果常作为一种定性的参考,并且该方法也不利于后续设计敏度的解析计算,这些缺陷使其难以应用于飞机结构的精细数值优化设计。
发明内容
为克服滚转速率形式定义的副翼舵面效率设计敏度难以解析求解的缺点,本发明中,分别采用势流理论面元法和有限元法建立机翼的气动模型和结构模型,在此基础上,进行弹性机翼副翼舵面效率及其设计敏度的计算。由前面的叙述可知,副翼舵面效率设计敏度求解的关键是弹性机翼单位副翼偏转角产生的滚转速率的设计敏度,但滚转速率和结构设计变量之间的关系难以直接显式确定,因此需要对副翼舵面效率的形式进一步转化。当飞行速度小于副翼反效速度时,副翼偏转产生的驱动力矩将使机翼绕其纵轴滚转,同时机翼的滚转运动将产生阻尼力矩,且阻尼力矩随着滚转速率的增大而增大。另外,阻尼力矩有使机翼滚转速率降低的趋势。当滚转速率达到某一值时,副翼偏转产生的驱动力矩和滚转运动产生的阻尼力矩之和为零,这时机翼达到稳定滚转状态。特别地,对于弹性机翼,弹性变形也将引起驱动力矩和滚转力矩的变化,从而使弹性机翼具有和刚性机翼不同的滚转速率。
基于以上思路,本发明将弹性机翼单位副翼偏转角产生的滚转速率转化为单位副翼偏转角产生的驱动力矩和单位滚转速率产生的阻尼力矩的组合形式。设计敏度工程数值方法是基于离散化的气动模型和结构模型实现的。驱动力矩或阻尼力矩为气动载荷和模型节点坐标的函数,而弹性机翼的气动载荷通过位移间接依赖于结构设计变量,由此即可通过解析法方式计算副翼舵面效率的设计敏度。具体说明如下,
本发明提出一种弹性机翼副翼舵面效率的设计敏度工程数值方法,包括如下步骤:
1)建立气动面元模型和结构有限元模型,并计算气动面元模型和结构有限元模型间数据传递矩阵;
2)利用步骤1)的计算结果,计算刚性机翼单位副翼偏转角产生的滚转速率;
3)利用步骤1)的计算结果,计算弹性机翼单位副翼偏转角产生的驱动力矩;
4)利用步骤1)的计算结果,计算弹性机翼单位滚转速率产生的阻尼力矩;
5)根据步骤2)、3)、4)的计算结果,进行弹性机翼副翼舵面效率及其设计敏度的计算。
优选地,所述步骤3),具体包括如下步骤:
31)在某一副翼偏转角下计算刚性机翼副翼偏转角产生的驱动力,将其作为初始气动载荷施加于结构有限元模型;
32)通过松耦合迭代的方式计算最终稳定的弹性变形和气动载荷,以及弹性变形状态下的驱动力矩;
33)获得弹性机翼单位副翼偏转角产生的驱动力矩。
优选地,所述步骤4),具体包括如下步骤:
41)在选定的滚转速率下计算刚性机翼滚转运动产生的阻尼力,将其作为初始气动载荷施加于结构有限元模型;
42)通过松耦合迭代的方式计算最终稳定的弹性变形和气动载荷,以及弹性变形状态下的阻尼力矩;
43)获得弹性机翼单位滚转速率产生的阻尼力矩。
优选地,所述步骤5)中弹性机翼副翼舵面效率的计算,通过以下公式实现:
式中,η为副翼舵面效率;ωr,β为刚性机翼单位副翼偏转角产生的滚转速率;β为副翼偏转角;fsezβ为弹性机翼由β产生的施加于结构模型的气动载荷;Δω为选定的滚转速率;fsezΔω为弹性机翼由Δω产生的施加于结构模型的气动载荷;yT为有限元模型节点坐标y分量组成的行向量。
优选地,所述步骤5)中是通过直接法或伴随法进行弹性机翼副翼舵面效率的设计敏度计算。
优选地,所述通过直接法进行弹性机翼副翼舵面效率的设计敏度计算中,首先通过解线性方程组的方式计算弹性机翼在副翼偏转角β下对应的平衡状态的结构位移列阵的设计敏度dusβ/db,以及弹性机翼在选定的滚转速率Δω下对应的平衡状态的结构位移列阵的设计敏度dusΔω/db,然后得:
式中,b为设计变量;usβ为弹性机翼在副翼偏转角β下对应的平衡状态的结构位移列阵;fsβ为弹性机翼在副翼偏转角β下对应的平衡状态的结构气动载荷列阵;usΔω为弹性机翼在选定的滚转速率Δω下对应的平衡状态的结构位移列阵;fsΔω为弹性机翼在选定的滚转速率Δω下对应的平衡状态的结构气动载荷列阵;B1和B2是布尔矩阵,其功能是分别按照fsezβ和fsezΔω中的元素在fsβ和fsΔω中的位置由dfsβ/db和dfsΔω/db得到dfsezβ/db和dfsezΔω/db;最终副翼舵面效率关于设计变量b的敏度可表示为:
优选地,所述通过伴随法进行弹性机翼副翼舵面效率的设计敏度计算中,首先计算以下两个伴随变量:
本发明的有益效果为:将弹性机翼单位副翼偏转角产生的滚转速率表示为单位副翼偏转角产生的驱动力矩和单位滚转速率产生的阻尼力矩的组合的形式,解决了滚转速率形式定义的副翼舵面效率的设计敏度难以解析求解的问题。在此基础上,实现了直接法或伴随法对副翼舵面效率设计敏度的工程计算,从而辅助后续飞机结构的精细数值优化设计。
附图说明
图1是本发明的算法流程示意图;
图2是机翼面元模型图;
图3是机翼整体有限元模型;
图4是机翼内部骨架有限元模型;
图5是副翼舵面效率随飞行速度变化的曲线;
图6是副翼舵面效率设计敏度随飞行速度变化的曲线;
图7是直接法的计算时间随设计变量数目变化的曲线。
图中标注说明如下:
1是机翼面元模型中副翼对应的子块;2是结构有限元模型第1个区域;3是结构有限元模型第2个区域;4是第1个固定点;5是第2个固定点;6是飞行速度,单位为m/s;7是副翼舵面效率;8是副翼舵面效率设计敏度;9是第1个设计变量对应的曲线;10是第2个设计变量对应的曲线;11是第3个设计变量对应的曲线;12是第4个设计变量对应的曲线;13是设计变量数目;14是计算时间,单位为s。
具体实施方式
设机翼所在的直角坐标系为Oxyz,其中机翼的弦平面位于x-y平面内,x轴沿机翼的弦向,y轴沿机翼的展向,z轴垂直于x-y平面向上。副翼舵面效率的具体定义形式为:
式中,η为副翼舵面效率;ωe,β为弹性机翼单位副翼偏转角产生的滚转速率;ωr,β为刚性机翼单位副翼偏转角产生的滚转速率。
式(1)中的副翼舵面效率是直接通过滚转速率形式进行定义的,当结构设计变量的值变化时,弹性机翼的滚转速率也将随之发生变化。为能计算滚转速率关于设计变量的敏度,需按照力学原理将式(1)转化为可计算表达的形式。
在稳定滚转状态下,机翼绕x轴的总滚转力矩为零,即:
(Mr,β+ΔMe,β)·β+(Mr,ω+ΔMe,ω)·ω=0 (2)
式中,β为副翼偏转角;ω为副翼偏转角β引起的弹性机翼滚转速率;Mr,β为刚性机翼单位副翼偏转角产生的滚转驱动力矩;Mr,ω为刚性机翼单位滚转速率产生的滚转阻尼力矩;ΔMe,β为弹性机翼单位副翼偏转角产生的驱动力矩相对于Mr,β的变化量;ΔMe,ω为弹性机翼单位滚转速率产生的阻尼力矩相对于Mr,ω的变化量。对于弹性机翼,可认为其弹性变形包括由副翼偏转所产生的弹性变形部分和由滚转运动所产生的弹性变形部分,ΔMe,β反映了单位副翼偏转角引起的弹性变形部分所产生的滚转驱动力矩变化量,ΔMe,ω反映了单位滚转速率变化引起的弹性变形部分所产生的滚转阻尼力矩的变化量。
根据式(2),可将式(1)中的ωe,β写为:
式中,fsezβ为弹性机翼由β产生的施加于结构模型的气动载荷;fsezω为弹性机翼由ω产生的施加于结构模型的气动载荷;yT为有限元模型节点坐标y分量组成的行向量。当机翼的弹性变形较小时,可认为弹性机翼由滚转速率引起的施加于结构模型的气动载荷随滚转速率的增大而线性变化。所以,可将弹性机翼副翼舵面效率写为以下形式:
式中,Δω为选定的滚转速率;fsezΔω为弹性机翼由Δω产生的施加于结构模型的气动载荷。由于原始机翼中弧面和迎角产生的气动载荷左右对称,其气动载荷对机翼的滚转力矩贡献为零,因此在计算弹性机翼的气动载荷时直接取机翼的迎角为零,并忽略原始机翼中弧面对下洗的影响。本发明中的副翼舵面效率设计敏度的工程分析方法流程示意图如图1所示。
本发明的具体实施步骤是:
步骤1:建立机翼的气动模型和结构模型,并计算气动模型和结构模型之间的数据传递矩阵。通过求解定常线性小扰动速度势方程可获得气动力影响系数矩阵,进而得到机翼的气动载荷,其一般形式的公式为:
fa=q·S·AIC·DW (5)
式中,fa为面元气动载荷列阵;q为飞行动压;S为各面元单元的面积组成的对角矩阵;AIC为气动力影响系数矩阵;DW为机翼下洗列阵。
采用有限元法建立机翼的结构模型,并选择合适的根部固支节点使机翼无刚体运动,其平衡方程可表示为:
K·us=fs (6)
式中,K为结构总体刚度矩阵;us为结构位移列阵;fs为结构气动载荷列阵。
式(5)中,机翼下洗列阵DW的元素为单位速度来流在相应的面元控制点处的法向速度分量,即由中弧面处的零法向流动边界条件确定。中弧面为机翼结构模型上、下表面节点垂向位移得到的中间位置曲面,显然fa是结构有限元模型节点位移的函数。此外,fs由fa传递至结构模型获得,结构位移列阵us需要通过式(6)解线性方程组获得。所以,弹性机翼的气动载荷由松耦合迭代的方式确定。在该过程中,需要气动模型和结构模型之间的数据传递,一方面需要将结构模型上、下表面节点垂向位移传递至面元控制点以计算下洗列阵,另一方面则需要将面元气动载荷传递至结构有限元模型。在本发明中,两模型之间的数据传递通过径向基函数插值的方法实现,首先进行位移传递矩阵的计算,然后根据虚功原理即可获得载荷传递矩阵,数据传递矩阵中的元素是气动模型和结构模型节点坐标的函数,与设计变量无关。
步骤2:计算刚性机翼单位副翼偏转角产生的滚转速率。副翼偏转角β产生的施加在面元模型上的气动载荷可表示为:
farβ=q·S·AIC·DWrβ (7)
式中,farβ为刚性机翼在副翼偏转角β下产生的面元气动载荷列阵;DWrβ为刚性机翼在副翼偏转角β下对应的下洗列阵。
设待求解的刚性机翼滚转速率为ωr,则ωr对应的下洗列阵DWrω可表示为:
式中,ya为面元控制点的展向坐标组成的列阵;V为来流速度。
进而刚性机翼在滚转速率ωr下产生的面元气动载荷列阵faωr可表示为:
farω=q·S·AIC·DWrω (9)
在稳定滚转状态下,farβ产生的驱动力矩和farω产生的阻尼力矩之和为0,此时仅ωr未知,据此即可解得刚性机翼的滚转速率。进一步,得
步骤3:计算弹性机翼单位副翼偏转角产生的滚转驱动力矩。该过程中不考虑副翼偏转产生的滚转运动,对于弹性机翼,以式(7)给出的刚性机翼在副翼偏转角β下产生的气动载荷作为初始载荷,将其传递至结构有限元模型计算弹性变形,然后求解因弹性变形产生的气动载荷,并与farβ一同施加于有限元模型,如此反复迭代计算直至达到平衡状态,这时有:
K·usβ=fsβ (11)
式中,K为结构总体刚度矩阵;usβ为弹性机翼在副翼偏转角β下对应的平衡状态的结构位移列阵;fsβ为弹性机翼在副翼偏转角β下对应的平衡状态的结构气动载荷列阵。进而可得到弹性机翼单位副翼偏转角所产生的滚转驱动力矩,即式(4)中的yT·fsezβ/β,其中fsezβ可通过提取fsβ中相应的元素获得。
步骤4:计算弹性机翼单位滚转速率产生的滚转阻尼力矩。对于选定的滚转速率Δω,则初始的刚性机翼滚转运动产生的面元气动载荷为:
faΔω=q·S·AIC·DWΔω (12)
式中,faΔω为刚性机翼在选定的滚转速率Δω下产生的面元气动载荷列阵;DWΔω为Δω对应的下洗列阵。以faΔω作为初始气动载荷,将其传递至结构有限元模型计算弹性变形,然后求解因弹性变形产生的气动载荷,并与faΔω一同作为外力施加于有限元模型,如此反复迭代计算直至达到平衡状态,则该情况下的结构平衡方程可表示为:
K·usΔω=fsΔω (13)
式中,K为结构总体刚度矩阵;usΔω为弹性机翼在选定的滚转速率Δω下对应的平衡状态的结构位移列阵;fsΔω为弹性机翼在选定的滚转速率Δω下对应的平衡状态的结构气动载荷列阵。进而可得到弹性机翼单位滚转速率引起的滚转阻尼力矩,即式(4)中的yT·fsezΔω/Δω,其中fsezΔω可通过提取fsΔω中相应的元素获得。
步骤5:副翼舵面效率设计敏度计算的基本原理。根据式(4),副翼舵面效率关于设计变量b的敏度可表示为:
由此可知,副翼舵面效率设计敏度求解的关键是计算fsezβ和fsezΔω的设计敏度,而fsezβ和fsezΔω分别与fsβ和fsΔω直接关联,所以,需要首先求解fsβ和fsΔω的设计敏度。式(11)和式(13)关于设计变量b求导可得:
dK/db·usβ和dK/db·usΔω为列阵,由于设计变量b对刚度矩阵的影响仅限于与该设计变量相关的单元,因此在计算中不必存储矩阵dK/db,只需要将相关联的单元刚度矩阵的敏度与其相应的节点位移相乘再叠加即可。这种方式可以减少计算机内存的消耗并提高计算效率。fsβ和fsΔω与结构模型上、下表面节点垂向位移具有直接显式关系,所以dfsβ/dusβ和dfsΔω/dusΔω可通过气动载荷和位移之间的关系解析推导获得。fsβ和fsΔω通过结构模型节点位移而间接依赖于设计变量,若由式(15)和式(16)解得dusβ/db和dusΔω/db,即可根据链式法则得到fsβ和fsΔω的设计敏度,进而得到副翼舵面效率的设计敏度。副翼舵面效率设计敏度的计算有两种方法,即直接法和伴随法,详细过程将在步骤6和步骤7进行说明。
步骤6:副翼舵面效率设计敏度的直接法。在利用直接法进行设计敏度的解析求解时,首先由式(15)和式(16)通过解线性方程组得到dusβ/db和dusΔω/db,进而得
在式(17)和式(18)中,B1和B2是布尔矩阵,其功能是分别按照fsezβ和fsezΔω中的元素在fsβ和fsΔω中的位置由dfsβ/db和dfsΔω/db得到dfsezβ/db和dfsezΔω/db。将式(17)和式(18)的计算结果代入式(14)即可获得副翼舵面效率的设计敏度。若结构设计变量数目为Ndv,则在直接法中,为了求解位移关于设计变量的敏度,需要进行Ndv×2次线性方程组的求解。显然,随着设计变量数目的增大,该部分的计算量也将越来越大,从而影响设计敏度计算的效率。
步骤7:副翼舵面效率设计敏度的伴随法。在伴随法中,需要引入两个伴随变量,即:
式中,为yT·dfsezβ/dusβ对应的伴随变量;为yT·dfsezΔω/dusΔω对应的伴随变量。以[K-dfsβ/dusβ]T作为线性方程组的系数矩阵,以-[dfsezβ/dusβ]·y作为右端项,可解得的转置;同理,以[K-dfsΔω/dusΔω]T作为系数矩阵,以-[dfsezΔω/dusΔω]·y作为右端项,可解得的转置。最终,使用伴随法计算副翼舵面效率设计敏度的公式可表示为:
实施例
本发明以某一大展弦比机翼作为具体实施例,对副翼舵面效率设计敏度的工程数值方法进行验证和说明。该机翼模型的几何参数如下:翼根弦长为2140mm,翼尖弦长为440mm,机翼的半展长为8900mm,前缘后掠角为9°。
步骤1:建立大展弦比机翼的气动模型和结构模型,并计算气动模型和结构模型之间的数据传递矩阵。面元模型的几何外形是刚性机翼在x-y平面的投影,如图2所示,该面元模型由6个面元子块组成,共包含气动单元802个。大气压强设置为101325Pa,飞行马赫数设置为0.6。
采用Hypermesh和Patran建立结构模型,如图3和图4所示,有限元模型共包含有限元节点746个、四边形壳单元904个、三角形壳单元86个、梁单元400个。机翼蒙皮为复合材料对称层合板[±45/04/±45/902]s,其单层属性如下:纵向弹性模量为1.40×1011Pa,横向弹性模量为8.60×109Pa,剪切模量为5.00×109Pa,泊松比为0.33,单层厚度为0.125mm。其它壳单元为钛合金材料,其属性如下:弹性模量为1.10×1011Pa,泊松比为0.3,壳单元厚度为3.0mm。该模型中有4个设计变量,均表示机翼蒙皮复合材料层合板0°铺层的厚度,其中第1个设计变量在结构有限元模型第1个区域的上表面,第2个设计变量在结构有限元模型第2个区域的上表面,第3个设计变量在结构有限元模型第1个区域的下表面,第4个设计变量在结构有限元模型第2个区域的下表面。
使用径向基函数插值方法计算数据传递矩阵时,以Wendland C2函数作为径向基函数,并取紧支半径为1500mm。在弹性机翼气动载荷计算过程中,为建立机翼中弧面,需要分别将结构模型上、下表面节点垂向位移传递至面元控制点,进而计算机翼下洗列阵。由于气动面元模型是二维的,而机翼结构模型是三维的,所以在实施例中将气动载荷施加于有限元模型的下表面节点。
步骤2:计算刚性机翼单位副翼偏转角产生的滚转速率。取副翼偏转角为1°,得刚性机翼的滚转驱动力矩为6904.97N·m,最终解得刚性机翼单位副翼偏转角(单位为°)产生的滚转速率为3.61°/s。
步骤3:计算弹性机翼单位副翼偏转角产生的驱动力矩。取副翼偏转角为1°,当弹性机翼前后两次迭代步最大翼尖垂向位移的相对误差小于0.1%时,则认为迭代收敛,最终得弹性机翼单位副翼偏转角(°)产生的滚转驱动力矩为2954.20N·m。
步骤4:计算弹性机翼单位滚转速率产生的阻尼力矩。选定的滚转速率为1°/s,当弹性机翼前后两次迭代步最大翼尖垂向位移的相对误差小于0.1%时,则认为迭代收敛,最终得弹性机翼单位滚转速率(单位为°/s)产生的滚转阻尼力矩为-1519.88N·m。
步骤5:根据式(4)得马赫数0.6时的副翼舵面效率为:1/3.61×2954.20/1519.88=53.84%。如图5所示,给出了副翼舵面效率随飞行速度变化的曲线,副翼舵面效率随着飞行速度的增大而减小,直至为零,此时发生副翼反效。由图可知,副翼反效速度约为315m/s。
步骤6,副翼舵面效率的设计敏度计算。分别采用直接法和伴随法进行副翼舵面效率的设计敏度计算,然后使用中心差分法对其计算结果的正确性进行验证。将设计变量b增大或减小较小的值Δb,可得到修改后的设计变量b+Δb或b-Δb。基于修改后的设计变量所对应的结构有限元模型进行副翼舵面效率的计算,之后即可得到副翼舵面效率设计敏度的中心差分结果,在本实施例中,Δb取为0.01mm。表1给出了直接法、伴随法和中心差分法的副翼舵面效率设计敏度的计算结果,其中直接法和伴随法的结果完全相同,中心差分法与直接法或伴随法的结果高度吻合。
由于后掠机翼的弯扭耦合效应,弹性机翼将产生低头扭转变形,从而使气动载荷减小,所以最终副翼舵面效率小于1。设计变量值的增大将使机翼的刚度增大,进一步使弹性机翼的低头扭转变形程度减小,最终引起副翼舵面效率的增大。所以本实施例中4个设计变量所对应的副翼舵面效率设计敏度均为正值。
图6给出了副翼舵面效率设计敏度随飞行速度变化的曲线。相对于靠近翼尖的设计变量,靠近翼根的设计变量的摄动将会对弹性变形产生更大的影响,所以在相同的飞行速度下,副翼舵面效率关于第1个设计变量的敏度大于关于第2个设计变量的敏度,同理,副翼舵面效率关于第3个设计变量的敏度大于关于第4个设计变量的敏度。结构有限元模型的第1个区域或第2个区域的两个设计变量对应的副翼舵面效率设计敏度不同,这是由结构模型的不对称性引起的。
直接法中位移的设计敏度和伴随法中的伴随变量都通过解线性方程组的方式获得,在该过程中,需要首先对线性方程组的系数矩阵进行LU分解,然后前后回代得到待求解项。相对于矩阵分解,前后回代的时间非常小。虽然分解的系数矩阵可以重复使用,但是随着设计变量数目的增多,直接法中由前后回代引起的时间花费也将越来越大。为了对直接法的计算时间进行测试,以本实施例中的4个设计变量作为一个设计变量组,然后不断地将这4个设计变量添加到该设计变量组并进行计算时间的测试。这里仅考虑前后回代的计算时间,图7给出了直接法的计算时间随设计变量数目变化的曲线。所以,在设计变量数目多的情况下,伴随法比直接法具有更高的计算效率。
以上实施例的计算结果说明了本发明中弹性机翼副翼舵面效率设计敏度的工程数值方法的有效性。
表1副翼舵面效率设计敏度的计算结果
直接法 | 伴随法 | 中心差分法 | |
第1个设计变量 | 3.1399E-2 | 3.1399E-2 | 3.1531E-2 |
第2个设计变量 | 2.4340E-2 | 2.4340E-2 | 2.4287E-2 |
第3个设计变量 | 2.1975E-2 | 2.1975E-2 | 2.2104E-2 |
第4个设计变量 | 1.8574E-2 | 1.8574E-2 | 1.8516E-2 |
Claims (7)
1.一种弹性机翼副翼舵面效率的设计敏度工程数值方法,其特征在于,包括如下步骤:
1)建立气动面元模型和结构有限元模型,并计算气动面元模型和结构有限元模型间数据传递矩阵;
2)利用步骤1)的计算结果,计算刚性机翼单位副翼偏转角产生的滚转速率;
3)利用步骤1)的计算结果,计算弹性机翼单位副翼偏转角产生的驱动力矩;
4)利用步骤1)的计算结果,计算弹性机翼单位滚转速率产生的阻尼力矩;
5)根据步骤2)、3)、4)的计算结果,进行弹性机翼副翼舵面效率及其设计敏度的计算。
2.根据权利要求1所述的一种弹性机翼副翼舵面效率的设计敏度工程数值方法,其特征在于,所述步骤3),具体包括如下步骤:
31)在某一副翼偏转角下计算刚性机翼副翼偏转角产生的驱动力,将其作为初始气动载荷施加于结构有限元模型;
32)通过松耦合迭代的方式计算最终稳定的弹性变形和气动载荷,以及弹性变形状态下的驱动力矩;
33)获得弹性机翼单位副翼偏转角产生的驱动力矩。
3.根据权利要求1所述的一种弹性机翼副翼舵面效率的设计敏度工程数值方法,其特征在于,所述步骤4),具体包括如下步骤:
41)在选定的滚转速率下计算刚性机翼滚转运动产生的阻尼力,将其作为初始气动载荷施加于结构有限元模型;
42)通过松耦合迭代的方式计算最终稳定的弹性变形和气动载荷,以及弹性变形状态下的阻尼力矩;
43)获得弹性机翼单位滚转速率产生的阻尼力矩。
5.根据权利要求4所述的一种弹性机翼副翼舵面效率的设计敏度工程数值方法,其特征在于,所述步骤5)中是通过直接法或伴随法进行弹性机翼副翼舵面效率的设计敏度计算。
6.根据权利要求5所述的一种弹性机翼副翼舵面效率的设计敏度工程数值方法,其特征在于,所述通过直接法进行弹性机翼副翼舵面效率的设计敏度计算中,首先通过解线性方程组的方式计算弹性机翼在副翼偏转角β下对应的平衡状态的结构位移列阵的设计敏度dusβ/db,以及弹性机翼在选定的滚转速率Δω下对应的平衡状态的结构位移列阵的设计敏度dusΔω/db,然后得:
式中,b为设计变量;usβ为弹性机翼在副翼偏转角β下对应的平衡状态的结构位移列阵;fsβ为弹性机翼在副翼偏转角β下对应的平衡状态的结构气动载荷列阵;usΔω为弹性机翼在选定的滚转速率Δω下对应的平衡状态的结构位移列阵;fsΔω为弹性机翼在选定的滚转速率Δω下对应的平衡状态的结构气动载荷列阵;B1和B2是布尔矩阵,其功能是分别按照fsezβ和fsezΔω中的元素在fsβ和fsΔω中的位置由dfsβ/db和dfsΔω/db得到dfsezβ/db和dfsezΔω/db;最终副翼舵面效率关于设计变量b的敏度表示为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011338381.6A CN112464372B (zh) | 2020-11-25 | 2020-11-25 | 一种弹性机翼副翼舵面效率的设计敏度工程数值方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011338381.6A CN112464372B (zh) | 2020-11-25 | 2020-11-25 | 一种弹性机翼副翼舵面效率的设计敏度工程数值方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112464372A true CN112464372A (zh) | 2021-03-09 |
CN112464372B CN112464372B (zh) | 2021-08-27 |
Family
ID=74798917
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011338381.6A Active CN112464372B (zh) | 2020-11-25 | 2020-11-25 | 一种弹性机翼副翼舵面效率的设计敏度工程数值方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112464372B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103678762A (zh) * | 2013-09-29 | 2014-03-26 | 北京航空航天大学 | 优化的复合材料机翼气动弹性风洞模型的缩比建模方法 |
CN103678763A (zh) * | 2013-10-14 | 2014-03-26 | 北京航空航天大学 | 复合材料机翼气动弹性剪裁方法及其遗传/敏度混合优化方法 |
CN105183996A (zh) * | 2015-09-14 | 2015-12-23 | 西北工业大学 | 面元修正与网格预先自适应计算方法 |
CN108052772A (zh) * | 2017-12-30 | 2018-05-18 | 北京航空航天大学 | 一种基于结构降阶模型的几何非线性静气动弹性分析方法 |
CN110704953A (zh) * | 2019-09-30 | 2020-01-17 | 西北工业大学 | 一种大展弦比机翼静气弹性能设计敏度的分析方法 |
CN111959817A (zh) * | 2020-07-29 | 2020-11-20 | 成都飞机工业(集团)有限责任公司 | 一种确定飞翼布局飞机机翼变形限制条件的方法 |
-
2020
- 2020-11-25 CN CN202011338381.6A patent/CN112464372B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103678762A (zh) * | 2013-09-29 | 2014-03-26 | 北京航空航天大学 | 优化的复合材料机翼气动弹性风洞模型的缩比建模方法 |
CN103678763A (zh) * | 2013-10-14 | 2014-03-26 | 北京航空航天大学 | 复合材料机翼气动弹性剪裁方法及其遗传/敏度混合优化方法 |
CN105183996A (zh) * | 2015-09-14 | 2015-12-23 | 西北工业大学 | 面元修正与网格预先自适应计算方法 |
CN108052772A (zh) * | 2017-12-30 | 2018-05-18 | 北京航空航天大学 | 一种基于结构降阶模型的几何非线性静气动弹性分析方法 |
CN110704953A (zh) * | 2019-09-30 | 2020-01-17 | 西北工业大学 | 一种大展弦比机翼静气弹性能设计敏度的分析方法 |
CN111959817A (zh) * | 2020-07-29 | 2020-11-20 | 成都飞机工业(集团)有限责任公司 | 一种确定飞翼布局飞机机翼变形限制条件的方法 |
Non-Patent Citations (3)
Title |
---|
N.S.KHOT*等: "《DEFORMATION OF A FLEXIBLE WING USING AN ACTUATING SYSTEM FOR A ROLLING MANEUVER WITHOUT AILERONS 》", 《AMERICAN INSTITUTE OF AERONAUTICS AND ASTRONAUTICS》 * |
刘祥等: "《一种弹性机翼的颤振主动抑制与阵风减缓方法》", 《西北工业大学学报》 * |
张保等: "《基于子结构法的大型结构数值敏度计算技术》", 《航空工程进展》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112464372B (zh) | 2021-08-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Tsushima et al. | Geometrically nonlinear static aeroelastic analysis of composite morphing wing with corrugated structures | |
CN102012953B (zh) | Cfd/csd耦合求解非线性气动弹性仿真方法 | |
Dillinger et al. | Stiffness optimization of composite wings with aeroelastic constraints | |
Cavallaro et al. | Prandtlplane joined wing: Body freedom flutter, limit cycle oscillation and freeplay studies | |
Smith et al. | CFD-based analysis of nonlinear aeroelastic behavior of high-aspect ratio wings | |
Howcroft et al. | Aeroelastic modelling of highly flexible wings | |
Brown | Integrated strain actuation in aircraft with highly flexible composite wings | |
Garcia | Numerical investigation of nonlinear aeroelastic effects on flexible high-aspect-ratio wings | |
Stanford et al. | Trim and structural optimization of subsonic transport wings using nonconventional aeroelastic tailoring | |
Palacios et al. | Static nonlinear aeroelasticity of flexible slender wings in compressible flow | |
CN110704953B (zh) | 一种大展弦比机翼静气弹性能设计敏度的分析方法 | |
Riso et al. | Parametric roll maneuverability analysis of a high-aspect-ratio-wing civil transport aircraft | |
CN113723027A (zh) | 一种针对弹性飞机的静气动弹性计算方法 | |
Xie et al. | Static aeroelastic analysis including geometric nonlinearities based on reduced order model | |
CN112580241B (zh) | 一种基于结构降阶模型的非线性气动弹性动响应分析方法 | |
CN113886967A (zh) | 多巡航工况的大型飞机机翼气动弹性优化方法 | |
Ghoman et al. | Multifidelity, multistrategy, and multidisciplinary design optimization environment | |
Amoozgar et al. | Aeroelastic stability analysis of aircraft wings with initial curvature | |
Castellani et al. | Flight loads prediction of high aspect ratio wing aircraft using multibody dynamics | |
Fugate et al. | Aero-Structural Modeling of the Truss-Braced Wing Aircraft Using Potential Method with Correction Methods for Transonic Viscous Flow and Wing-Strut Interference Aerodynamics | |
CN112580240B (zh) | 一种适用于复杂大柔性飞机建模的非线性子结构方法 | |
Cavagna et al. | A fast tool for structural sizing, aeroelastic analysis and optimization in aircraft conceptual design | |
Kuzmina et al. | Analysis of static and dynamic aeroelastic characteristics of airplane in transonic flow | |
Reinbold et al. | Aeroelastic simulations of a delta wing with a Chimera approach for deflected control surfaces | |
CN112464372B (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 |