复杂地形上空超低空风场的工程仿真方法

为描述飞机在低空突防时的外部环境,用多重山脉和多种地貌设计构造了飞机超低空飞行时的复杂地形,同时,基于位势流动理论建立了复杂地形及其上空风场的工程化仿真模型,并进行了算例仿真和分析.该方法可以满足分析飞机做超低空飞行的动力学特性的需要.     

柔性飞机大变形状态飞行载荷分析

飞行载荷分析中综合考虑了结构大变形的几何非线性效应以及曲面气动力效应.飞机结构由相互连接的可以表征任意变形的几何非线性欧拉梁表示,升力面由展向以及弦向分布的涡环网格表示.计算时通过曲面网格动态跟踪结构的变形不断修改气动力影响系数矩阵,反复迭代求解气动力和结构变形直至结构变形收敛.给出了风洞试验机翼模型以及无约束定常平飞模型算例.计算结果表明:在小变形阶段该方法与线性计算方法的结果基本一致,非线性效应不明显;大变形阶段由于曲面气动力的非线性效应,计算结果与线性方法的有显著差别.分析表明该方法在大变形阶段的计算结果比线性结果更为合理,数值计算时间短,适合于工程快速分析.     

光流控制地形跟随与自动着陆

将一种自研的光流传感器安装到飞行平台上,利用光流信息控制飞行器的地形跟随和自动着陆.测试了光流传感器的测量性能,建立了飞行器的运动模型,开发了基于光流反馈的地形跟随控制算法,提出了光流自动着陆的控制规律,并进行了数学仿真和硬件在回路的实时仿真验证,实验结果表明:光流信息可以一种相对简单的形式来巧妙地同时控制飞行速度和飞行高度,保证飞行器安全地执行多种飞行任务,实验证明了光流传感器在飞行控制中具有工程可行性.     

三维真实地形环境下无人机救援航路规划方法

利用无人机(UAV)的三维飞行能力,采用优化方法规划路径,能够使其在救援任务中比地面车辆以更短的时间到达救援区域,提高救援效率.针对真实的地理环境,根据无人机约束采用均匀化网格方法进行地形建模,之后根据地形数据的特点设计适合数学计算与求解的数据结构.最后设计了包含偏离代价、高度代价、地形跟随/回避代价、威胁代价和安全距离代价的综合性能指标函数,并采用航路点交叉和网格搜索代替航路点搜索的方法,对蚁群算法进行改进完成航路规划.仿真结果表明:本文方法能够直接处理三维地形数据,在保持地貌的前提下,完成了无人机的三维航路规划任务,得到满足无人机约束的三维最优航路,提高了航路规划方法的实用价值.      

无人直升机三维避障方法及仿真

自主避障是无人直升机（UH）在复杂未知环境中执行飞行任务所必须的能力。针对未知环境下UH避障飞行问题,提出了一种新的三维实时避障方法,通过将三维激光雷达探测视野划分成若干不同半径的扇形柱区和等角度分布的直方柱区,结合雷达探测数据判断障碍物分布情况,UH根据避障方法执行相应的规避动作来实现避障飞行。同时提出一种基于CATIA三维建模软件二次开发的三维仿真方法,该方法可以对未知的飞行环境进行三维建模,结合避障方法实现UH在其中的真实避障飞行模拟,并利用该仿真方法验证了三维避障方法的可行性。      

充液航天器质量特性计算方法研究

充液航天器的质量特性对飞行动力学建模和控制设计有着重要影响。质量特性测量装置通常只能测出干重状态下航天器的质量特性，而充液状态下的航天器质量特性只能靠计算给出，其中如何计及各贮箱内液体的惯性张量是此类工程计算的难点。首先采用势流理论推导了任意形状贮箱内液体质量特性的一般表达式，包括全充液情况和部分充液情况，其中对部分充液情况考虑了自由液面小幅晃动的影响；然后采用三维有限元方法建立了液体质量特性的数值计算方法，并通过一个验证算例确认了该计算方法的有效性；最后应用该方法计算了一个充液航天器的质量特性，并研究了不同的液体处理方法对计算结果的影响。研究表明，工程上通常采用的“固化液体”处理方法所计算出来的转动惯量明显偏大，而所提出的计算方法具有更好的理论正确性和工程适用性。     

风切变场中直升机前飞状态动态响应

为解决复杂风切变场中直升机运动的数值仿真,把握非线性直升机模型在风场中的响应特性,建立了一维、二维的风切变模型,采用直升机非线性飞行动力学模型,考虑直升机的控制系统,用四阶龙格－库塔法数值仿真直升机在风场中的前飞运动.对水平风和铅垂风的各个梯度场对响应的影响进行比较.结果表明,水平风和铅垂风的各个梯度场对响应的影响是不同的,此研究中二维空间的组合梯度场构成更为危险的飞行环境;在不同的情况下其作用是不可忽略的,尤其注意在低空飞行状态下对直升机飞行性能的影响.也进一步说明此数值仿真法对非线性直升机模型在复杂风场中响应研究实用有效.     

地形仿真实现

在应用数字仿真法研究导弹地形跟踪飞行的撞地概率时,首先要对导弹飞越的真实地形进行数值模拟.在充分探讨了地形模型和随机过程实现的基础上,采用波叠加法仿真了3种典型地形.研究了子波数对仿真结果的影响,说明选取一定数目的子波即可得到满意的结果.该工作对武器系统评估及撞地概率研究有重要应用价值.      

多翼导弹复杂流场及气动力的数值计算

将改进的Jameson格式应用到多翼导弹的三维流场计算中,并使用分区技术和三维椭圆微分方程数值网格生成技术以精确模拟此复杂构形的复杂流态,计算结果与美国导弹标模的实验吻合很好.一系列攻角下导弹三维流场的数值模拟表明,修正Jameson中心差分格式能够适用于复杂外形、复杂流态、亚/跨/超/高超音速的流场计算,而且精度和分辨率都比较高,结果可信.     

基于试验气动力的纵向机动飞行载荷分析

提出了一种基于非线性风洞试验数据(简称试验气动力)同步进行纵向机动飞行过程和飞行载荷的计算分析方法.通过曲面样条插值将试验气动力的升力系数和俯仰力矩系数引入纵向机动飞行过程的计算,同时基于试验气动力的压力分布进行全机飞行载荷的静气动弹性修正.弹性载荷的计算采用线性气动力影响系数矩阵处理.采用偶极子格网法提供气动力(简称理论气动力).将基于这两种气动力的计算结果进行了对比.结果表明:基于试验气动力的机动飞行过程更为合理;基于理论气动力的载荷随机动飞行过程的变化趋势有可能与基于试验气动力的载荷变化趋势不一致.      

基于涡方法生成大涡模拟进口条件的数值计算

探索改进涡方法来生成大涡模拟的非定常进口条件.改进方法中,为了避免局部漩涡数量过多或者过少,采用密度分布方式放置漩涡场.并用Langevin方程控制漩涡运动,模拟实验方法中的蜂窝器,使改进后的涡方法生成的脉动速度场更加符合湍流的特征.在已知雷诺平均的流场结果下,利用涡方法产生漩涡场,进而生成能满足大涡模拟所需要的非定常进口流场.为了检验改进后的涡方法在生成脉动速度过程中的效果,在槽道中进行了对比数值试验,并借助于直接数值模拟数据做对比,对比分析槽道进出口的平均速度、涡量以及雷诺应力统计,证明改进后的涡方法在生成大涡模拟进口条件下是非常有效的.      

一类TVD型组合差分方法及其在磁流体数值计算中的应用

根据太阳风数值模拟的特点，考虑到算法的质量（收敛速度、稳定性、精度等），结合磁流体数值计算的特性，对三维球坐标磁流体动力学（MHD）方程组中的流体部分采用一种修正Lax-Friedrichs差分法而对磁场部分采用MacComack格式，发展了一类快捷的具有TVD特性的组合数值新方法，作为格式的检验，在一维情况下，将其与PPM格式进行了比较，对一维快慢磁流体激波问题得到了与PPM格式精度相同的结果，然后将其诮到定态太阳风的数值模拟上，在不同等离子体β情形下，可得到理想的太阳风定态结构，为今后将此数值模式应用到具有复杂磁场位型或三维直实太阳风暴的数值模拟研究奠定了基础。     

北斗-罗兰C组合导航的ASF修正计算与研究

介绍讨论了修正ASF的米林顿法和CLC数值积分法,并用MATLAB编程计算;用米林顿法计算的(SF+ASF),与军标"长波地波传输信道计算方法"计算结果非常吻合;CLC数值积分法计算的结果与英国PaulWilliams的结果也很吻合;而CLC数值积分法计算的(SF+ASF)与米林顿法有较大差别,但在用以推求陆地路径ASF时,二者差别就很小。试验表明,CLC数值积分法有助于复杂地形路径的ASF修正计算。     

无振荡、无自由参数格式在三维太阳风流动数值模拟中的应用

从快速、有效地进行空间天气数值预报的需要出发,针对1998年5月份行星际太阳风暴事件,采用全三维流体力学(HD)模型进行数值试验,考察了无振荡、无自由参数(NND)格式在三维太阳风流动数值模拟中的应用.计算中内边界的密度分布由观测K日冕数据来确定,由此得到的源表面密度分布具有和源表面电流片相似的结构.数值试验表明: 虽然该格式不需要加入人工粘性,但是具有很好的计算稳定性.在扰动计算时激波间断耗散小,在三维任何方向上间断所占的网格数比较少,没有数值色散现象.      

自动空中加油阶段加油机尾涡流场建模与仿真

针对自动空中加油技术对于加油机尾涡精确建模的需要,通过分析大型加油机机翼尾涡形成机理和运动规律,运用改进的马蹄涡理论方法建立了较精确的考虑尾涡衰减和扩散特性的尾涡空间流场计算模型,之后将此流场模型在受油机质心处线性化,并采用平均化的方法计算得到受油机受到的平均风分量和风梯度值,在分析受油机受扰运动中利用气动等效的方法将尾涡对受油机的影响模型加入到受油机非线性全量方程中.在此基础上利用Matlab/Simulink搭建仿真平台对空中加油阶段受油机扰动运动进行了仿真验证.仿真结果表明,所建立的尾涡模型能够反映受油机受扰后的动态特性,与有人机飞行员的飞行经验相符,该模型还可用于大型飞机起飞飞行安全和近距编队飞行的扰动分析.     

固体火箭发动机燃烧室三维流动数值计算

采用k-ε双方程湍流模型,以SIMPLE计算程式求解翼柱型装药固体火箭发动机燃烧室内三维非定常不可压流N-S方程.在建立翼柱型装药简化模型的基础上采用边界标志法来表达燃面推移,采用从二维到三维的初场给定方法结合多重网格法求出了多个时间步的非稳态流场结构.计算结果表明燃烧室内旋涡运动呈现一定的空间与时间分布,周向加质导致燃气通道横断面上轴向速度分布的极不均衡.      

嵌套闭磁场结构内CME产生和传播的数值模拟

给出了特殊类型的日冕物质抛射(CME)数值模拟定性结果,这种CME核心闭磁场结构前半部分磁力线的方向与太阳整体偶极场磁力线的方向相反．计算结果表明,这种CME核心闭磁场结构磁力线与太阳整体偶极场反向磁力线之间存在过渡磁场结构,在向外传播时过渡磁场结构所占的面积逐渐增大．这一结果可以用来解释飞船为什么能够观测到一类双极磁云,这类磁云前半部分磁场方向与太阳整体偶极场方向相反．为了模拟这一数值结果,强调需要采用包含嵌套闭磁场的冕流背景结构,并在合适的位置触发CME．     

无振荡、无自由参数格式在数值模拟盔型磁场中太阳风流动的应用

提出了一种离散二维三分量理想磁流体力学守恒型方程组（MHD）的算法（NNDMHD），它有效地控制了磁场散度不为零误差对动量方程组的影响，把气体动力学中计算跨音速流动问题的有效算法——无振荡、无自由参数（NND）格式推广应用到MHD方程组中。利用该算法首先对常见一维和二维算例进行数值试验，得到比较好的结果，消除了间断处的非物理振荡。然后对太阳风在子午面轴对称盔形磁场位形中流动进行数值试验，在这个算例中，物理量沿径向变化大，NNDMHD格式仍然能够有效地控制磁场散度离散不为零误差导致的非物理流动。这个算例的计算结果表明：在网格划分比通常情况和稀4倍时，该算法仍保持很好的计算稳定性。     

三角翼前缘涡破裂点脉动的实验分析

在水槽中和风洞中分别进行了流动显示和动态测压实验,目的是研究三角翼前缘涡破裂点的脉动现象.对流动显示图片中涡的破裂点位置进行统计和频谱分析,表明破裂点振荡存在双主频特征,位于0.07位置的主频对应着螺旋—泡—螺旋的转化过程,位于0.2~0.4之间的主频对应着螺旋形态破裂本身的小幅振动.对统计的破裂点位置数据做低通滤波后进行相关性分析,还表明了三角翼左右2个前缘涡的破裂点位置信号具有负相关性.进行三角翼表面动态压力测量,对压力数据低通滤波后做相关性分析,发现相同的一条前缘涡的动态压力信号在整个翼面上都具有高度的相关性.      

基于PSO-BP神经网络的平流层风场短期快速预测

平流层风场环境对临近空间低速飞行器驻空飞行性能有重要影响。研究了基于PSO-BP神经网络的平流层区域风场建模与快速预测方法,根据历史风场数据,采用主成分分析法对数据进行降维处理,通过BP神经网络对风场进行预测建模,利用粒子群优化（PSO）算法对其进行优化,采用Biharmonic样条曲面插值方法构建区域预测风场。以南海地区5年历史风场为对象,对比分析了基于BP神经网络和基于PSO-BP神经网络的风场预测模型,结果表明:使用具有全局寻优特性的PSO算法改进BP神经网络,能够有效避免传统BP神经网络易陷入局部最优的缺点,提高预测精度;通过结合PSO-BP神经网络预测与Biharmonic样条曲面插值,可实现区域风场的预测。研究结果可为临近空间低速飞行器的轨迹规划与区域驻留等任务的高精度区域快速预报风场提供解决途径。