神经网络在高空长航时无人机翼型多点优化中的应用

NSGAII算法在翼型多点设计中有着广泛的应用价值,然而其巨大的计算资源和计算成本限制了它的使用。为了解决这个问题,本文引入具有较强非线性映射能力的神经网络代理模型,采用实验设计结合BP法训练神经网络响应面来代替N-S方程求解翼型的性能。在实验点的数值模拟中,为了进一步节省计算资源,提高计算效率,采用网格的变形代替网格的重新划分,使得计算网格的更新速度提高了约50%。在翼型的参数化过程中,采用改进的PARSEC方法,用较少的参数实现了翼型的精确控制。为了增强神经网络的泛化能力,采用12-7-4-3-1的隐层结构。对NLF1015翼型的多点优化算例表明,此方法不仅显著降低了整个优化过程的计算量,而且对翼型的气动性能预测也具有较高的可信度,在高空长航时无人机的翼型设计中具有一定的潜力。     

基于MLPG混合配点法的接触问题形状优化

分别介绍了接触问题和MLPG(Meshless Local Petrov-Galerkin)混合配点法的理论基础,推导了相关公式,给出了两种典型接触状态的定解条件,使用二维线性基函数,采用三次样条曲线权函数,通过移动最小二乘法插值,将MLPG混合配点法运用到接触分析中,使用罚函数法添加本质边界条件,对二维弹性接触问题的接触过程进行模拟,反复迭代得到真实的接触情况,建立了一种新的应力-位移非线性数学求解模型,结合遗传算法对实际工程接触问题进行了求解优化,给出优化结果和目标函数变化曲线,并与相关文献结果比较,验证了该方法的有效性.     

CFD中数值耗散和延伸边界对声传播的影响

对已知色散关系的正弦波在二维管道中的传播使用有限体积法进行计算,结果表明二阶精度下每波长40个左右的网格点能够有效地保留声波.在计算域外增加稀疏网格延伸段能够对边界反射起到抑制作用,为以较小的计算量增加实现较好的反射抑制效果,研究给出了延伸段长度和网格轴向尺寸的推荐公式.      

网格之间流场信息传递的高精度方法和验证

理论分析了插值方法引入误差的原因,提出一种利用时间推进方法和移动网格技术实现网格之间流场信息传递的思想,并构造计算方法。在介绍新的信息传递算法原理的基础上,根据三维情况下传递过程需要满足的时间推进步长稳定性条件,推导出最小推进步数的要求。最后,通过对激波诱导物体运动的非定常流动模拟,比较了信息传递算法和线性插值、加权最小二乘插值方法,计算结果表明,与插值方法相比,信息传递过程没有引入误差,能有效抑制非物理振荡,保持流场计算精度。     

基于非结构网格的跨声速静气动弹性计算

基于非结构网格开发了静气动弹性计算程序,通过气动力与弹性变形的交替迭代计算实现流固耦合。气动力求解采用RANS方程,结构变形计算采用柔度系数法,通过面样条(IPS)方法进行柔度系数插值,使用改进弹簧近似法实现网格变形。通过对AGARD 445.6机翼的静气弹计算验证了方法的正确性,同时对机翼的静气动弹性特性进行了分析,计算表明,弹性变形使其升力线斜率下降、焦点前移。     

多体分离数值模拟及方法验证

利用自主研发的混合网格Euler／RANs方程解算软件WoF90以及六自由度动力学与运动学方程解算程序,采用脚本语言技术实现网格的自动生成,建立了一套可用于多体分离问题的准定常Euler流动非结构网格方法。通过AEDC算例对非结构准定常Euler方程方法解决多体分离问题进行了计算研究,验证了方法的可行性。     

基于“各向异性”四面体网格聚合的复杂外形混合网格生成方法

复杂外形粘性流动的计算网格生成一直是阻碍CFD实际工程应用的瓶颈问题之一,虽然三棱柱/四面体混合网格的出现极大地缓解了这一问题,但是在飞行器几何外形非常复杂的情况下,在粘性边界层内生成高质量的贴体三棱柱网格仍然是非常困难的事情。本文从"各向异性"四面体网格的生成算法中得到启示,发展了一种基于"各向异性"四面体网格聚合的混合网格生成方法。通过对DLR-F6翼身组合体和F16战斗机外形的混合网格生成及数值模拟表明,该方法生成的混合网格具有网格质量好、自动化程度高等特点,对提高复杂外形飞行器关键气动参数的模拟精度具有重要意义。     

适用于混合网格的并行GMRES+LU-SGS方法

给出了一种适用于混合网格的并行无矩阵GMRES+LU-SGS隐式时间格式。首先采用LU-SGS方法迭代若干步以获得一个合适的初场,然后切换到GMRES方法在每一时间步内近似求解,并将LU-SGS方法作为其预处理器。为加速收敛,将CFL数随着残差的降低逐步放大;为减少存储量和计算量,通量Jacobian采用无矩阵处理。在保证与串行执行一致的前提下,采用基于共享内存的OpenMP方法实现了并行计算,并通过对网格的分组避免了内存争夺。算例验证表明,方法极大地提高了计算收敛效率,并行结果与串行结果完全一致,计算结果与实验结果吻合较好。     

利用CFD技术计算飞行器动导数

利用非定常流场数值计算方法模拟飞行器强迫俯仰振荡仅能得到俯仰力矩系数对迎角变化率和俯仰角速度的动导数之和,而动稳定性分析需要单独的动导数数值.为解决这个问题,利用滑移网格模拟强迫振荡运动,得到俯仰力矩系数对迎角变化率和俯仰角速度的动导数之和.利用旋转参考坐标系模拟定常拉升,得到俯仰力矩系数对俯仰角速度的动导数.利用旋转参考坐标系模拟匀速滚转,得到滚转力矩系数对滚转角速度的动导数.对有翼导弹和水上飞机进行了纵向和横向动导数的计算.计算结果与试验数据、文献数据以及其他方法得到的结果具有较好的一致性,表明提出的方法可用于复杂外形飞行器动导数计算.