一种新的界面映射推进方法及其在气动弹性力学中的应用

非线性气动弹性研究中,涉及到非线性的流体动力学(CFD)和非线性的结构动力学(CSD )耦合问题,在耦合界面上不仅要进行两场之间信息的传递,而且要有匹配的时间推进格式 。针对计算流体动力学（CFD）和计算结构动力学（CSD）的耦合计算方法,基于边界能 量守恒,发展了一种新的界面处理方法,包括：① 界面映射方法,该方法基于耦合边界上 局部的网格信息,将两场之间载荷信息和位移信息的转换放在同一个映射矩阵中来处理,克 服了占用大量CPU时间和内存的需求;② 界面推进方法,该法基于松耦合计算流程,通过引 入半步交错推进,达到了二阶时间精度,提高了计算精度和效率。最后将该界面处理方法应 用于柔性大展弦比机翼的气动弹性计算和AGARD445.6机翼的动响应分析中。结果表明该方法 能够高效、高精度地处理不同物理场之间的数据交换和时间推进,并具有处理复杂几何体非 线性气动弹性问题的能力。     

二阶时间精度的CFD/CSD耦合算法研究

非线性气动弹性分析中,涉及到非线性流体动力学(CFD)和非线性结构动力学(CSD)的耦合计算问题。本文分析比较了目前几种耦合算法:全耦合、松耦合和紧耦合,并从耦合边界能量传递守恒上就松耦合及紧耦合方法进行了时间精度的分析,得出传统松耦合中即使流体和结构子系统达到高阶时间精度,耦合算法的时间精度仍仅为一阶,紧耦合方法虽然可以达到二阶时间精度,但没有明显提高计算效率。随后,本文基于松耦合流程改进了耦合格式,分析表明其具有二阶时间精度,通过AGARD445.6机翼颤振模型的算例验证,说明该方法可以在保证计算精度的基础上明显提高计算效率,并保持了传统松耦合方法模块化的优点,在模拟非线性气动弹性问题时具有很高的优越性。     

基于混合遗传算法的多冲量最优变轨

针对航天器多冲量最优变轨问题,建立了多冲量最优变轨的数值优化模型,给出了一种遗传算法(GA)与序列二次规划算法(SQP)结合的混合优化算法.该算法不需初值猜测,全局和局部搜索能力强并且计算效率高.仿真计算了燃料最优交会和时间最短拦截问题,比较了GA,SQP以及混合遗传算法的性能.针对混合遗传算法得到的不同冲量次数变轨的优化结果,分析了冲量次数对变轨性能指标的影响.结果表明混合遗传算法综合性能最高,冲量次数对不同性能指标的影响不同.仿真算例验证了模型和算法的有效性.     

低速流场中柔性悬臂板的后颤振响应

建立了一个新的非线性气动弹性模型,对低速流场中柔性悬臂板的后颤振响应特性进行了分析。建模中考虑了结构几何非线性、气动力非线性以及两者之间的强耦合效应。通过实验数据对所建立的气动弹性模型进行了验证。发现在低速流场中柔性悬臂板可能会以周期加倍的方式进入混沌运动。结构几何非线性效应和翼尖涡引起的非定常气动力效应对柔性悬臂板的结构响应有显著影响,而尾涡变形引起的非定常气动力对结构运动的影响较小。还研究了不同耦合算法的差异,给出了小展弦比大柔性结构非线性气动弹性数值仿真时耦合策略的选择依据。     

一种CFD/CSD耦合计算方法

针对柔性大展弦比机翼气动弹性分析和主动弹性机翼(AAW)设计发展了一种计算流体动力学(CFD)和计算结构动力学(CSD)的耦合计算方法。其主要思想是采用在同一物理时间弱耦合求解CFD/CSD技术。气动力采用非定常N-S方程的双时间有限体积求解技术,结构响应则采用有限元数值求解技术。CFD和CSD耦合计算的边界信息(气动力和网格)由所设计的界面程序传输。网格信息传输采用守恒体积转换(CVT)方法将CSD计算结构响应位移插值到CFD网格点上。变形已有的CFD网格技术用以确定CFD的变形网格。以位移或载荷的迭代误差为判断耦合计算的收敛标准。最后得到了机翼在Ma=0.8395,α=5.06°时CFD/CSD耦合计算的收敛值。针对计算结果分析了机翼受静气动弹性过程中结构响应和气动特性随时间变化的效应。初步研究结果表明:这种弱耦合方法求解非线性气动弹性问题是可行的。     

非线性壁板颤振分析

利用一种改进的计算流体力学与计算结构动力学(CFD/CSD)耦合方法研究了由气动和结构几何非线性引起的壁板颤振问题。在非定常气动力计算中,考虑了通量分裂格式、隐式时间推进方法和几何守恒律;二维和三维壁板的结构几何非线性建模则采用了有限元的协同旋转理论,并利用一种近似能量守恒算法求解结构的非线性响应。流场和结构求解器采用二阶松耦合方法联立求解,并将其应用于壁板在超声速、跨声速和亚声速的颤振计算中。当考虑结构几何非线性和气动非线性时,出现了典型的极限环振荡现象,并对颤振边界和极限环振荡幅度进行了对比分析。     

基于径向基函数的网格变形及非线性气动弹性时域仿真研究

为开展非线性气动弹性研究,基于非线性结构有限元软件NASTRAN和自主研制的多块结构化计算流体力学(CFD)求解器,开发了一套基于计算流体力学/计算结构动力学(CFD/CSD)耦合求解方法的气动弹性时域仿真程序.该程序采用径向基函数(RBF)交换两套求解器之间的数据并进行网格变形.为提高RBF方法的效率,构造了基于多次插值的空间待插值点精简算法.在多次插值过程中,每次插值的对象为上次插值的误差,并同时限制插值区域,以此实现了空间待插值网格数的精简.数个网格变形的算例表明该方法可支持大变形运动,并且具有较高的计算效率.采用此程序开展了AGARD 445.6机翼颤振计算、大展弦比机翼的静气动弹性计算与切尖三角翼极限环振荡(LCO)现象的动气动弹性仿真,结果揭示了当机翼展弦比较大或者响应幅值较大时,结构非线性对于气动弹性有显著影响.     

改进的CFD/CSD耦合系统设计方法及其计算

论述了界面上的数据交换方法和时间推进格式在流体动力学（CFD）和结构动力学（CSD）耦合求解非线性气动弹性问题中的关键作用;改进了CFD／CSD耦合系统．包括基于边界元方法设计了CFD和CSD耦合界面的数据转换方法,该方法可在同一映射矩阵处理结构响应和非定常气动载荷转换,保证了耦合边界上的能量守恒;改进了一种松耦合方法流程,在保持模块化基础上提高了计算精度和效率。最后将本文方法应用于二维位移外插算例和A—GARD445．6机翼的动响应分析中,结果表明本文方法具有较高精度。     

双三角翼及其翼身组合的滚转运动特性比较研究

采用Navier-Stokes方程与滚转运动方程耦合计算方法,比较研究了不同后掠角的双三角翼和翼身组合体的滚转运动特性,分析了机翼前缘后掠角及细长机身对非定常滚转力矩时滞环、动态流场结构和物面瞬时压力分布的影响。研究结果表明：主翼迎风面上的融合涡能量在80°/60°双三角翼上耗散较小,而在76°/40°双三角翼上耗散严重,这是造成两模型滚转力矩稳定性与时滞特性差异的主要因素;机身对气流的扰动作用,大幅增强了滚转力矩的线性分量;机身对气流的上洗作用,增强了边条涡与融合涡吸力及其时滞性,同时加剧了主翼背风面的两涡干扰;大滚转角时机身对横流流动的干扰,使得主翼背风面压力分布的时滞差异显著增加。该研究结果有助于认识后掠角与细长机身影响双三角翼滚转运动特性的物理机理。