基于加速度频响函数小波分解的模型修正方法

为提高模型修正效率，满足修正方法对于实测环境噪声的鲁棒性，将Kriging模型和小波分解引入加速度频响函数模型修正。首先，将加速度频响函数进行小波分解，用得到的第1层幅值较大的小波系数来表征原频响函数。其次，采用拉丁超立方抽样对初选待修正参数进行设计，根据设计结果对各参数进行灵敏度分析，从而确定模型修正的待修正参数，以待修正参数作为Kriging模型输入，所对应的小波系数作为Kriging模型输出，通过混合灰狼算法寻得最优Kriging模型相关系数，建立精确有效的Kriging模型。最后，以目标加速度频响函数小波分解的小波系数与Kriging模型输出的小波系数差值最小为目标，通过水循环算法求解模型待修正参数。数值算例表明，所提模型修正方法具有良好的修正效果，即使在加速度频响函数中加入信噪比为5 dB的高斯白噪声时，修正误差也低于4%，证明了该方法对于随机噪声的鲁棒性。     

基于代理模型的气动外形平面参数多目标匹配设计

将Kriging代理模型和Pareto遗传算法引入气动外形平面参数匹配设计中,提出一种基于代理模型的多目标平面参数匹配设计方法。将拉丁超立方试验设计用于平面参数筛选,确定出参数匹配方案库;基于方案库的计算流体力学(CFD)分析结果构建Kriging气动代理模型;将Kriging模型替代CFD分析,用于气体布局参数匹配优化设计,提高了设计效率并保证了可信度;通过Pareto遗传算法优化解决多点设计要求下气动布局参数匹配问题,一次性给出参数匹配方案的最优解集,从Pareto前沿中根据设计偏向选择气动布局最佳匹配方案。以典型的双后掠布局平面参数多点匹配优化设计问题作为算例,研究结果表明:Kriging气动代理模型与实际CFD分析结果的误差均小于5%,满足精度要求;根据不同设计偏向,选择了3种参数匹配Pareto优化方案,与原样本方案相比超声速阻力减小6.0%~12.8%,跨声速升阻比增加0.01%~3.40%,证明了匹配设计方法的有效性;通过试验设计的Pareto分析与主、交互效应分析,获得了气动布局平面参数对气动性能影响的定量关系,能够为参数匹配设计提供依据。所提出的平面参数匹配设计方法可应用于其他常规与非常规气动布局型式。     

二维气动问题中Kriging代理模型精度影响因素

以二维跨声速临界翼型的阻力特性为对象,探讨样本点数目、Kriging代理模型参数及其类型等对模型精度的影响.阻力系数采用计算流体力学（CFD）方法得到.模型精度的验证采用交叉验证方法,采用平均误差、最大误差和标准交叉验证残差来衡量Kriging代理模型的精度.研究结果表明:①Kriging代理模型预测气动阻力效果较好.②模型精度随样本点的增多而提高,剔除与样本点响应趋势不相符的“奇异点”后,模型精度显著提高,平均误差减小5%~38%,最大误差减小13%~77%.③核函数类型对模型精度的影响最大,相关参数次之,回归模型的影响最小.采用高斯相关函数、2阶多项式回归模型,以及合适的相关参数值时,Kriging代理模型的精度最高.      

指尖密封泄漏特性数值计算与试验研究

为了研究不同压差和转速下指尖密封的泄漏特性,建立了CFD多孔介质数值计算模型,计算分析了某型航空发动机指尖密封泄漏特性。为有效验证数值计算模型的有效性,进行了一定压差、转速工况下指尖密封泄漏特性试验。针对压差和转子离心膨胀导致的多孔介质参数改变,根据试验结果对计算方法进行了分析,引入了压差修正系数C_k和离心膨胀修正系数k,并对修正系数进行了确定及验证。结果表明,未引入压差修正系数C_k和离心膨胀修正系数k时,计算得到的指尖密封泄漏特性与试验值会随着压差和转速的增加而偏离;引入压差修正系数C_k和离心膨胀修正系数k后,计算结果与试验值一致,未随着压差和转速的增加而偏离,最大偏差小于11%,平均偏差小于6%。     

航空发动机静子系统中螺栓连接结构的建模与修正

基于概率设计方法对航空发动机静子系统中螺栓连接结构进行模型修正研究。静子系统中不同的机匣结构由螺栓连接在一起,采用具有一定材料属性的层单元法来等效模拟螺栓连接结构,通过调节层单元材料属性来提高有限元模型精度。首先采用灵敏度分析方法确定修正区域,然后基于概率方法设计修正参数,利用蒙特卡洛抽样法选取样本点,最后基于静力学试验数据选择合理的目标函数,通过迭代优化计算得到修正后的模型参数。结果表明,修正后的静子系统有限元模型计算得到的刚度值与试验数据更加接近,修正后的相对误差都在5%以内。     

航空发动机静子系统中螺栓连接结构的建模与修正（英文）

基于概率设计方法对航空发动机静子系统中螺栓连接结构进行模型修正研究。静子系统中不同的机匣结构由螺栓连接在一起,采用具有一定材料属性的层单元法来等效模拟螺栓连接结构,通过调节层单元材料属性来提高有限元模型精度。首先采用灵敏度分析方法确定修正区域,然后基于概率方法设计修正参数,利用蒙特卡洛抽样法选取样本点,最后基于静力学试验数据选择合理的目标函数,通过迭代优化计算得到修正后的模型参数。结果表明,修正后的静子系统有限元模型计算得到的刚度值与试验数据更加接近,修正后的相对误差都在5%以内。     

改进的综合CFD与畸变合成法预测动态畸变(英文)

对综合CFD与畸变合成方法进行了改进。通过CFD求解,获得AIP上各稳态气动参数。采用BP人工神经网络直接求解高度非线性的湍流相关方程,该方程将CFD求解的稳态气动参数与试验测得湍流度进行关联。经过训练的网络可根据稳态CFD结果预测湍流度。之后根据CFD计算的稳态总压和神经网络预测的湍流度合成动态压力。最后根据合成的动态压力寻找最大瞬时畸变。文中使用了六个工况的飞行数据验证该方法的准确度并且与原方法的计算结果进行了比较。结果表明该方法预测的畸变值与试验畸变值吻合较好,且预测结果较原方法与试验值更接近,该方法能够在试验前用来预测动态畸变,具有很强的工程价值。     

耦合梯度与分级Kriging模型的高效气动优化方法

Kriging代理模型中引入梯度信息能够提高模型的预测精度,但常规耦合梯度的方法都有不足之处。本文结合分级Kriging模型,提出了一种变可信度分级Kriging模型耦合梯度（GEHK）信息的新方法。首先利用梯度信息,选取扰动步长得到初始样本点附近的派生点,以派生点拟合出能够预测目标函数趋势的低精度Kriging模型。然后以初始样本点修正该模型得到高精度的Kriging模型。翼型减阻优化设计算例表明,与常规耦合梯度的Kriging模型相比,基于分级Kriging的梯度耦合方法对于扰动步长引入的误差不敏感,明显提高了模型预测精度,优化效率因此提升并使得目标函数值下降得更加迅速。相比Euler解作为低精度数据的常规分级Kriging模型,由梯度信息得到的派生点为模型提供了更准确的全局趋势预测,取得了更好的优化结果。本文方法成功应用于翼型多点减阻优化问题,说明该方法对复杂设计问题有很好的适应性。基于分级Kriging模型的耦合梯度方法克服了常规方法的缺点,提高了模型全局拟合精度,是一种优化效率更高的Kriging方法。     

全机有限元模型修正及验证方法研究

有限元分析现已成为飞行器设计的重要工具,若要提高有限元仿真的准确性,需要对有限元模型进行必要的修正。以某型直升机全机静力试验有限元模型为研究对象,基于静力试验实测的应变值,进行变量参敏度分析和模型优化,针对全机试验测点多、模型自由度大的特点,提出适用于验证全机模型的相关性算法;基于该算法,应用某一工况下的试验结果,对全机有限元模型的部分参数进行修正,得到与试验结果符合较好的模型参数结果;应用其他工况的结果对修正后的模型进行验证,结果表明:其他工况下修正后模型的计算结果与试验值也很接近,证明了该相关性算法和全机模型修正及验证方法的可行性和有效性。     

Kriging模型在固体火箭发动机装药可靠性分析中的应用

基于Kriging方法对固体火箭发动机装药结构可靠性进行了分析.首先,采用拉丁超立方抽样技术得到了随机参数样本点;然后,利用样本点数据对黏弹性有限元模型的装药进行了结构力学分析;最后,通过Kriging模型得到了八面体剪应变的函数近似响应面,并结合失效判据对失效函数进行了可靠性计算.为了验证Kriging模型在可靠性计算中的有效性,其计算结果与Monte-Carlo方法进行了比较,得到了十分接近的可靠度.      

应力强度因子的区间分析方法

 传统上用来处理不确定性问题比较有效的方法是概率分析方法,本文应用区间分析方法对具有不确定参数的应力强度因子进行估计。该方法以区间数学为基础,将不确定参数描述为区间变量;再利用Taylor级数展开通过区间运算得到应力强度因子的区间范围,从而为工程设计提供可信的数据。区间分析方法优于传统的概率分析方法的是：它不需要预先知道关于不确定参数大量的统计数据信息,并且具有计算方法简便、实用和精度高的特点。最后,通过对两种裂纹情况的计算将区间分析方法和传统的概率分析方法进行比较,说明了该方法的有效性。     

基于特征正交分解法的翼型结冰冰形快速预测

为加快翼型结冰冰形的计算速度,提出了一套基于特征正交分解(POD)法的结冰冰形快速预测算法.通过计算流体动力学(CFD)数值模拟计算得到的冰形结果作为样本,以结冰温度单参数变化时为例,详细介绍了POD法预测结冰冰形的实现步骤.考虑结冰温度、结冰时间以及液态水含量的影响,完成了单参数、两参数与三参数变化时的POD结冰冰形快速预测.通过算例结果发现:完成POD预测计算只需几秒钟且POD法与CFD法得到的冰形吻合较好,仅在冰角附近有较小差别,表明POD方法能够快速、准确地得到翼型的结冰冰形.      

高超声速化学非平衡流动Navier-Stokes/DSMC耦合算法

基于相同的化学反应模型,在已有计算流体力学（CFD）和直接仿真蒙特卡罗（DSMC）方法及程序的基础上,采用Modular Particle-Continuum（MPC）耦合技术,建立了包含化学非平衡Navier-Stokes/DSMC耦合算法。算法结构中DSMC计算区域在CFD计算结果上根据当地克努森数自动选取。发展了适用于流场分区信息交换的亚松弛技术,抑制DSMC方法对CFD计算的影响。把DSMC方法和CFD的应用范围拓展到过渡流区,为复杂飞行器近连续过渡流区高超声速化学非平衡流动数值模拟研究提供了一种工程适用的预测分析手段。通过对二维圆柱高超声速化学非平衡绕流的算例与其他结果的比较研究,表明耦合算法不论在流场结构、流场非平衡现象,还是飞行器表面参数、整体气动力/热特性方面,都能够得到与全DSMC计算吻合的结果,证实了所建立的Navier-Stokes/DSMC耦合计算模型与方法的有效性和可靠性。仿真了某航天器解体碎片在过渡区的化学非平衡流动,得到碎片在过渡区的气动力/热特性,为碎片的陨落计算提供依据。     

氧换热器在两相流条件下的压力脉动机理分析

氢氧液体火箭发动机氧换热器在研制过程中出现了压力脉动现象，影响了换热器的正常工作。以该型换热器为研究对象，将换热器分为燃气侧、汇总管、圆管等结构，进行了CFD仿真计算，与试验结果对比，修正了传热相关参数，建立了换热器理论计算模型用于参数计算；利用非线性方法分析试验数据，确定了该换热器为混沌系统；利用两相流分析与电路类比两种方法分别针对换热器进行了建模与计算，均捕捉到了相变下换热器内的压力脉动。结果表明：通过CFD仿真及试验修正后的参数用于计算换热器的稳态参数，具有较好的一致性；换热器出现压力脉动时，具有一定非线性特征，且压力脉动越大，非线性程度越高；该换热器出现压力脉动的原因为氧工况接近两相区，相变使得介质密度发生剧变。根据分析结果提出了提高换热器设计压力的改进措施，通过试验得到了有效验证。     

CFD法中的动态失速模拟

本文介绍了CFD动态失速模拟所需要捕捉的关键物理,CFD法一般情况和较先进的CFD法,最后重点介绍了CFD法动态失速模拟的当前水平.总的来说在动态失速预测方面,CFD法优于综合法,CFD/CSD耦合结果与试验结果之间的吻合度还是可以接受的.俯仰力矩曲线的相位、大小和形状与试验数据吻合得较好,耦合解与试验数据之间的主要偏差在于法向力预测上.     

基于粒子群算法的导弹PID控制参数选择

PID控制在导弹控制系统的设计中得到了较多的应用,但PID控制器的性能依赖于其参数的整定和优化。利用粒子群算法对控制器的参数进行优化,并对俯仰通道控制系统进行了仿真研究,结果表明,攻角和俯仰角速度都取得了较好的效果,当气动参数具有不确定性时,效果依然较好,验证了方法的有效性。     

盒式布局飞机的纵向气动参数优化研究

采用CFD方法探讨了纵向布局参数对飞机纵向气动特性的影响规律,可为盒式机翼气动布局设计提供参考。示例计算结果表明,盒式布局可减小飞机的巡航阻力,改善气动效率,但是其设计参数的选取需要进行多学科综合设计优化。     

防热试验用高温超声速燃气流场热环境分析

对用于防热系统试验考核的高温超声速燃气流场的流场结构和热环境特性进行了计算和分析.将CFD模拟方法与工程计算方法相结合,对两类锥形喷管形成的燃气流场进行了数值模拟计算和流场结构分析,并与试验过程中高速拍摄的图片进行了比对.利用模拟得到的流场气动参数,计算了沿流场轴向位置球头驻点模型的热流和压力,并与试验测量数据进行了对比,结果吻合较好.锥形喷管形成压缩波和膨胀波交替的流场结构,导致流场气动参数与模型热流和压力形成正弦变化的特点.     

基于变精度遗传算法的翼型快速优化设计方法

低碳环保的电动飞机在要求较高升阻比的同时,需要尽量降低成本、缩短研制周期。但高精度的数值模拟时间代价很大,因此针对电动飞机翼型设计中初始翼型较难选取、优化速度较慢的问题,提出了一种基于变精度遗传算法的翼型多点快速优化方法。以常用的 Hicks-Henne 型函数为基础,改进了其对翼型后缘描述不精确的问题。在数值模拟阶段,介绍了一种快速气动参数计算软件XFOIL,并分析了该软件的适用性与局限。之后给出了使用XFOIL 与 Matlab 进行联合求解的方法,在无人干预的情况下完全实现了翼型设计与优化的自动化,提高了设计效率。在翼型优化阶段,为保持较高的精度和寻优效率,设计了翼型参数的实数编码方法。针对传统遗传优化算法了改进,设计了染色体变精度杂交方法以及动态惩罚方法。最后,给出了基于遗传算法的多点优化方案,以及翼型多目标快速优化一体化设计方案。仿真分成两部分进行,首先改进的 Hicks-Henne 型函数能够有效实现参数化翼型的后缘夹角改变。通过与 NSGA-II 方法的优化结果对比,本文的方法在一定迭代次数范围内获得的升阻比更高,失速特性更加缓和,特别是在综合提高翼型优化效率方面表现较好。仿真结果表明,该方法能够快速获得多种工况下具有较高升阻比的翼型,也可以作为进一步优化的初始翼型,能提高翼型优化效率。     

含不确定参数的复合材料壁板热颤振分析

研究了含有不确定结构参数的壁板颤振问题,利用vonKarman大变形应变-位移关系、气动力活塞理论和准定常热应力理论建立了复合材料壁板颤振的气动弹性力学模型,考虑在壁板颤振分析模型中存在的不确定参数,将其用区间向量定量化,基于区间扩张理论和Taylor级数展开,并结合有限元计算方法,提出了区间分析的方法来估计含有不确定参数的壁板结构颤振临界风速的区间,以及发生极限环振动时振幅的变化区间。通过数值算例,将本文提出的壁板颤振的区间有限元模型与随机有限元模型进行了比较,显示了本文方法的有效性和可行性。这种方法的优点是只需要知道不确定参数的所在范围界限,为解决含有不确定参数的壁板颤振这类复杂的气动弹性动力学问题提供了一个途径。