基于结构动网格的无人机地面效应研究

为了更好地分析地面效应对无人机气动特性的影响,采用一种新的结构动网格生成技术,运用CFD手段对某无人机的地面效应进行了研究。在一套静态网格CFD计算基础上,运用结构动网格生成技术,获得非定常的网格移动,得到一系列的不同攻角、地面效应高度的CFD网格及初始流场。在此基础上进行定常流场求解,计算得到地面效应的定常CFD结果。最后对某双尾撑布局无人机的地面效应进行了CFD计算研究,对相应的全机气动特性的影响进行了分析。计算结果表明：地面效应对平尾的影响大于对机翼的影响并使平尾的失速迎角提前;该方法是可行的,具有一定的工程应用价值。     

适应于极区的惯导系统传递对准技术研究

针对飞机、舰艇全球航行的需求,以及由地球经线在极区收敛而导致的以真北向作为航向参考时的导航算法失效问题,在传统的“速度+姿态”匹配对准及修正算法的基础上,提出了格网导航系下的“速度+姿态”传递对准方法。基于格网导航系下惯性导航算法的编排,建立了格网坐标系的惯导误差模型,推导了格网导航系“速度+姿态”匹配的量测方程,利用卡尔曼滤波器对子惯导的速度、失准角、陀螺漂移、加速度计零偏、主子惯导间的安装误差角等进行了估计与修正。数学仿真及基于极区航行试验数据的半实物仿真结果表明,采用本文提出的格网导航系下的“速度+姿态”匹配传递对准算法,可实现中低精度惯导系统在高纬度地区时间为30s、水平对准精度优于1角分、方位对准精度优于6角分的快速、高精度对准。     

高超声速热流高精度数值模拟方法

通过高超声速下钝锥热流的计算,对影响热流计算精度的因素进行了综合研究.获得了壁面法向网格Re数对热流计算的影响规律.采用边界高阶插值,提高了热流计算精度.用Roe和AUSM(Advection Upstream Splitting Method)+格式搭配minmod和混合limiter两种限制器,探讨了空间离散方法及限制器的耗散性对热流结果的影响.应用五阶WENO(Weighted Essentially Non-Oscillatory)方法,对高阶格式计算热流的性能进行了细致研究,最终确立了网格依赖小,热流精度高的计算方法.研究认为:高阶格式不但可以放宽热流计算对网格Re数的限制,同时使得格式对本身计算方法的依赖性也变小.      

适用于混合网格的改进雅可比迭代法及其应用

LU-SGS因有较高的鲁棒性和小的内存需求而得到广泛应用,然而用于混合网格计算前需要进行网格排序和分组来实现算法并行;此外,LU-SGS格式收敛效率不高。针对这些缺点,本文提出了一种改进的适用于复杂混合网格的雅可比迭代方法,无需网格排序和分组就可实现算法的并行化,且有较快的收敛速度。该方法编程实现简单,易于采用OpenMP实现并行。算例研究表明,相比于LU-SGS格式,在各来流条件下,本文提出的方法收敛速度更快,鲁棒性好,并行和串行结果一致,且内存需求增加很少。      

一种改进的投影装置标定方法

摘要： 提出双目结构光三维测量系统中投影装置参数标定的一种改进方法.利用投影装置将格雷码图案投向若干待测场景,两个相机拍摄和解码同样数量的待测场景照片,采用三角形准则恢复标定角点在空间中的位置,仿照相机标定解算的方法获得投影装置标定参数.该方法优势在于不需要特定的标定平面或者标定靶标,可在测量过程中直接对投影装置标定,也可在预标定之后在测量过程中对投影装置参数进行更新.试验验证方法的正确性和有效性.     

点火冲击波与药柱裂纹流固耦合作用机理数值模拟研究

流体计算采用基于非结构动网格的有限体积方法,结构计算采用分段精细时程积分方法,求解结构动力学有限元方程,流体计算和固体计算之间采用松耦合在不断运动的交界面上进行作用力和运动特性的传递,对固体火箭发动机点火形成的冲击波与药柱缺陷形成的弹性裂纹之间相互耦合作用过程进行了模拟。研究表明,作用在弹性裂纹上的冲击波载荷引起的壁面变形随时间变化明显,开始表现为撑开裂纹出口,随后出现了收缩出口撑开底部的瓶状;固体变形对冲击波在弹性狭缝内的压力波振荡影响较小,以前采用刚体模型得到的研究结论仍有指导意义。     

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

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

部分膜弹流的膜厚与二次压力峰公式

用部分膜弹流模型来模拟线接触的混合润滑,提出浮动网格法,详细计算了二次压力峰的位置与大小,最小油膜厚度与中心油膜厚度。首次提出了线接触部分膜弹流的油膜厚度与二次压力峰计算公式。对浮动网格法,用 6种方法,即加、减、乘、除、对数与幂,进行了说明。理论分析的数学模型是基于Patir与 Cheng的平均流模型与微凸体接触模型。     

非结构网格、混合网格下计算方法研究

为了提高计算精度，本文研究了计算格式和网格技术。通过比较非结构网格下的NND格式、MUSCL格式、ENO格式及最小二乘格式的计算精度和波动，选择利用NND方法中的守恒变量法来建立非结构网格下的计算格式，通过改进NND方法中的限制器来提高计算精度和分辨率，计算结果比较显示出该格式具有很高的分辨率和计算精度。本文也研究了非定常网格自适应技术，计算给出了双马赫反射问题的新旧限制器及网格自适应计算结果比较，新限制器网格自适应后计算得到的滑移线出现了剪切涡现象，这表明，数值计算中取足够密的网格和构造高精度格式是非常必要的。     

带控制舵弹体气动特性分析

本文用数值模拟方法研究了带控制舵弹体的气动特性。该问题外形复杂，在控制舵和弹身之间会出现明显的缝隙。为了解决此类外形的气动力计算问题，我们发展了结构网格和非结构网格的混合网格技术。这样一方面充分发挥了贴体结构网格质量好、效率高的优点，另一方面充分发挥了非结构网格几何灵活性好、易于模拟复杂外形的优点。数值计算与风洞实验相当吻合，说明本方法具有良好的精度。最后对两种外形进行大量计算并进行气动特性分析，给出了静稳定特性和配平特性的研究结果。