非结构弹性体动网格方法应用策略及改进

针对3种代表性的经典弹性体动网格方法,从网格最大变形能力、质量评价指标和空间分布趋势3个方面综合分析了各方法的网格变形特性,重点总结了各方法的适用范围.结合工程实际,研究了特定问题的弹性体动网格应用策略.最后,整合各经典弹性体动网格方法的优势提出了改进方法,并通过NACA 0006翼型的超声速俯仰振荡加以验证.研究结果表明:基于面积的弹性体法适用于加密型网格,而基于距离的弹性体法适用于均匀型网格,该应用策略可为工程问题提供指导.改进方法效果明显,在大变形时相对基于距离的弹性体法的各项网格质量增幅均在10%以上,相对基于面积的弹性体法的最小网格质量增幅几乎达到4倍,并且在数值模拟中具有较高的精度,可显著增强弹性体法的鲁棒性.      

两种典型非结构网格变形方法特性对比研究

非定常流动仿真中常常遇到运动边界的情形,基于物理模型的线性弹簧法和基于数学方法的径向基函数(RBFs)插值法是实现计算网格随运动边界变化的两种主要的非结构网格变形方法。以NACA0012翼型的俯仰运动为例,对上述两种方法的变形特性进行较详细地定量分析研究,比较线性弹簧法和RBFs方法的CPU计算时间和最大变形能力,分析两种方法的网格质量随俯仰角度和运动步数的关系,并给出典型状态下网格单元质量云图和流场计算结果。结果表明:基于RBFs插值的网格变形方法计算效率较高,变形能力较强且能保持较高的网格质量,是一种高效的网格变形方法。     

非结构动网格生成的弹簧-插值联合方法

 提出了基于弹簧法的两点改进,以解决边界发生大变形时的非结构网格变形问题。为了使边界运动引起的网格变形能更好地由边界传递到内部网格中,提出了一种新的基于Delaunay网格插值的弹簧倔强系数逐层改进方法。为了提高弹簧法的计算效率,引入背景网格和直接插值方法,提出了弹簧-插值法。弹簧-插值法首先生成计算域的背景网格(粗网格),然后由弹簧法求解边界运动引起的背景网格变形,最后利用变形后的背景网格直接插值得到计算网格的变形。算例结果表明：改进后的方法一方面有效地提高了动网格的变形能力和变形后的网格质量;另一方面通过降低弹簧法的求解规模,显著地提高了动网格的变形效率。     

复杂外形的动态混合网格生成方法

本文将层推进方法、四分树方法和阵面推进法相结合,建立了二维复杂外形的动态混合网格生成方法.首先用层推进方法在运动物体附近生成贴体的四边形网格,然后用四分树方法在外场生成自适应的矩形网格,而中间由三角形网格过渡.三角形网格由阵面推进法生成,并由节点松弛和"对角线交换"方法进行优化.当物体运动时,贴体的四边形网格随物体运动而运动,而外场的矩形网格保持静止.中间过渡层的三角形网格随物体运动而变形.当物体运动位移较大导致出现质量较差的单元甚至导致网格相交时,则在局部重新生成网格.利用上述动态混合网格生成方法,生成了多种含相对运动的复杂外形的动态混合网格,网格质量令人满意.     

气动弹性计算中网格变形方法研究进展

网格变形是气动弹性计算中实现计算网格随运动边界变形的主要方法。在总结目前网格变形方法发展现状的基础上,对近几年常用的网格变形方法,即弹簧法、弹性体法、超限插值法、Delaunay背景网格法、径向基函数插值法和温度体法等做了简要的总结。根据各方法构建模型的不同,将它们分成物理模型法、数学插值法和混合方法3类,简要介绍了各方法的基本思想和研究进展,重点比较了各方法的网格变形特性（变形能力、变形质量和变形效率）和优缺点。总结了当前流场边界与结构边界的数据传递方法,对当前气动弹性计算中遇到的网格变形的难点问题作了简要评述并对未来网格变形方法的发展方向进行了探讨。     

基于非结构背景网格的动网格方法研究

借鉴弹簧动网格方法和Delaunay背景图方法,提出一种基于非结构背景网格的动网格生成方法。先用少量的物面控制点生成相对较粗的非结构背景网格,再用弹簧法将动边界的变化量作用到该背景网格,然后利用计算网格与背景网格的映射关系插值生成对应的动网格,最后对物面网格点及其附近受影响的密网格点进行局部网格修正,获得最终计算所需的动网格。所提方法与直接弹簧法比较,由于弹簧法作用的背景网格远粗于实际的计算网格,提高了计算效率,此外,由于网格的运动受控于受弹簧法作用的背景网格,边界变形或运动能较为均匀地传递到内部网格上,使网格单元过渡光滑,提高了动网格的质量。给出若干算例,展示了方法的效率和可行性。     

基于背景网格簇的动网格生成方法

提出了基于背景网格簇实现网格变形的新方法。背景网格簇由任一内边界节点和远场边界角点生成的背景网格构成,并随边界运动而变化。通过保持计算网格节点在每一个背景网格中所在单元的面积坐标不变(三维时保持体积坐标不变)而求出一组期望位置,加权平均该组期望位置,确定出该计算网格节点的新位置。其中权值与内边界节点和计算网格节点距离的倒数相关。新方法引入的背景网格簇可以有效地改进单个背景网格的不足之处。算例结果表明:基于背景网格簇的动网格生成方法实现简单,与弹簧近似法相比,新方法因不需要迭代求解方程组而非常高效,且拥有更强的网格变形能力;与Delaunay图映射法相比,该方法背景网格的单元个数极少,因此易于定位,且不会出现背景网格单元交叉的现象,网格变形能力更强,变形后的网格质量更好。     

柔性扑翼非定常流场的数值计算方法

提出一种将Delaunay图映射网格变形技术和非结构嵌套网格方法结合使用的策略,解决网格变形和嵌套网格单独用于柔性扑翼流场计算时需要网格再生的问题。该方法为嵌套网格中的每个嵌于背景网格的贴体非结构网格生成Delaunay背景图;每个时间步,根据扑翼的运动和变形规律移动背景图,再根据网格点和背景图的映射关系移动网格点,之后自动完成嵌套边界的定义和插值关系的建立。为方便嵌套关系的建立,嵌套网格进行分层管理。也研究了一种内存消耗少、效率较高的搜索算法,以及格心格式和格点格式统一的边界拓宽算法。非定常可压缩Navier-Stokes方程在非结构的动态网格上用有限体积法离散,并用预处理的双时间步推进、隐式LU-SGS迭代求解。几个扑翼算例的结果表明,该方法充分利用了Delaunay图映射网格变形方法的高效率,同时也发挥了嵌套网格处理大幅运动的优势;用于既有整体大幅扑动又有局部小变形的柔性扑翼流场计算,可取得令人满意的精度和效率。     

基于伴随方程的网格自适应及误差修正

基于流场方程的离散伴随优化理论和三维非结构网格,建立了网格自适应技术和目标函数误差修正方法。详细研究了用流动伴随变量进行目标函数的误差估计和修正技术,构造了适用于格心格式有限体积法的流场变量插值技术和网格单元剖分判据,初步实现了网格物面投影和空间单元优化,发展了适用于有限体积法的整套网格自适应方法。对NACA0012翼型和ONERA-M6机翼绕流进行了自适应数值模拟,并对升、阻力等目标函数进行了误差修正。数值结果表明,本文自适应方法能正确地捕捉到影响目标函数计算精度的敏感区域,网格自适应和误差修正两项技术显著提高了升、阻力等目标函数的计算精度。     

基于网格参数化变形的单级涡轮多学科可靠性设计优化

基于自由变形技术发展了一种可以协调实现流场、结构分析网格参数化变形的方法,并实现单级涡轮的多学科可靠性设计优化。针对流场分析网格展向积叠的特征,将三维网格变形转化为不同叶高二维截面的网格变形;根据二维流场分析网格的拓扑结构,设计二维控制面,基于叶型设计参数建立叶型控制点、其他控制点伴随叶型控制点移动,实现了网格的参数化变形,避免了网格畸变、提高了网格变形的质量;不同叶高的二维控制面组成控制体,实现三维流场分析网格的参数化变形。采用相同的控制体进行结构分析网格的变形,实现了耦合界面网格协调变形。基于Kriging代理模型进行单级涡轮的多学科可靠性优化,在满足可靠性约束的情况下效率提高4.97%、寿命提高40.86%,证明了方法的有效性。      

双欧法与四元数法的应用比较

对解决欧拉方程奇异性的双欧拉和四元数法进行了对比研究，四元数法从理论上讲比较完美，但实际应用中存在较大的累积计算误差，从而影响计算精度；双欧法利用正、反欧拉方程间精华区倒挂关系进行分区交替运算，把精华区扩展到全域，不仅根除了奇异性，而且计算误差小。因此，对于解决欧拉方程奇异性来讲，双欧法要优于四元数法。     

基于最小区域法的形位公差精确算法

最小区域法是真正符合国际形位公差定义的评价方法。针对现有的最小区域法的实现算法存在思路复杂、对初始计算条件要求较高等问题,提出1种逐次逼近算法。详细阐述了依据此算法求解4种常见形位公差的初始计算条件设置和求解步骤,绘制了以圆度计算为例的算法流程图。结果表明:该算法初始计算条件和思路简单,易于编程,对测点坐标分布及运行环境无特殊要求;可以准确有效地实现圆度等形位公差的评定。     

SQL优化方法

应用SQL不仅要确保它的正确性,同时也要保证SQL有良好的执行性能。本文分析、介绍了几种常用的SQL优化方法,包括设计使用合理的索引、减少不充分连接条件、减少格式转化与操作步骤、使用存储过程,通过这些方法,可以达到优化SQL性能的目的。     

翼型的气动最优化设计方法和反设计方法＊

本文将数值优化方法同计算流体动力学（ＣＦＤ）相结合,形成两种跨声速翼型的气动设计方法,即：最优化方法和反设计方法。采用求解Ｅｕｌｅｒ方程的有限体积法计算流场,通过结合优化算法,从正反两个方向对跨声速翼型进行气动优化设计。实例证明它们是翼型气动设计的有效方法,有较高的工程应用价值。     

翼型的气动最优化设计方法和反设计方法

本文将数值优化方法同计算流体动力学(CFD)相结合,形成两种跨声速翼型的气动设计方法,即:最优化方法和反设计方法.采用求解Euler方程的有限体积法计算流场,通过结合优化算法,从正反两个方向对跨声速翼型进行气动优化设计.实例证明它们是翼型气动设计的有效方法,有较高的工程应用价值.     

基于改进的关联度矩阵法维修可达性评价方法研究

为了更好地对产品的维修可达性进行综合评价,根据维修可达性的定义确定了评价指标体系,结合权的最小平方法对关联度矩阵法进行改进,构建了新的维修可达性的评价方法.该方法可对维修可达性进行量化处理,评价结果更符合实际情况.以货机货运系统动力驱动装置为例进行实例验证,结果表明,该评价方法可为货运系统维修可达性设计提供重要参考.     

基于声测法的齿轮行波共振测试方法

针对航空发动机由行波共振导致的中央传动锥齿轮断裂故障，提出了一种基于声测法的齿轮行波共振转速测试方法， 建立了非接触齿轮行波共振测量系统，该测试方法简便易行，不需要对转子件进行测试改装，直接引出测试部位的声信号进行分 析研究。基于声测法开展了航空发动机中央传动锥齿轮行波共振试验，研究中央传动锥齿轮行波共振特性，获取了从动锥齿轮行 波共振频率、危险转速范围。通过动应力同步测量试验进行对比验证，结果表明：该测试方法和系统可以精准识别和有效监测航 空发动机中央传动锥齿轮行波共振频率和共振转速，行波共振频率误差小于0.025%，共振转速误差小于0.018%。该方法同时适 用于其它齿轮和盘轴的行波共振故障分析，为相关类型齿轮共振试验及故障诊断提供了一种非接触、长时间监控疲劳试验的测量 方法。     

一种基于循环阈值的海上固定目标清洗方法

针对海上电磁环境复杂、数据量巨大等问题,在进行态势分析时,将目标海域内对态势分析影响较小的固定目标进行清洗,可有效缩短数据处理的时间、降低复杂度,提高态势分析的效率和直观性。现有的固定目标清洗方法鲁棒性较低,不能排除异常跳变数据的干扰,对数据集中信息的准确性要求较高,不适用于海上复杂环境的数据处理。因此,文章提出了一种基于循环阈值的海上固定目标清洗方法,并在仿真军事场景上进行了实验分析,所提方法通过循环阈值过程能够排除自然因素、量测误差和异常跳变数据的干扰,有效检测出样本数据集中的固定目标。     

基于改进的二维Otsu法的图像分割法

针对Otsu阈值法仅适合图像目标和背景分布呈正态分布且各像素量和方差基本相当的不足,提出了一种基于二维Otsu的改进算法,实验结果表明,该算法比传统Otsu阈值法能够获得较好的分割效果.     

基于蒙特卡罗方法的飞机成本参数法估算研究

研发早期的成本估算是项目立项和设计决策的基础,参数法可以用最少量的信息估算成本,而蒙特卡罗方法可以进行成本数据的不确定性分析,以明确成本风险。本文研究了参数法和蒙特卡罗方法在飞机成本估算中的应用。