地固系GNSS/SINS组合导航正常重力间接算法

相比导航系,GNSS/SINS组合导航在地固系计算,编排效率提高约40％,Kalman滤波效率提高约10％,即使正常重力与姿态矩阵的计算效率偏低,在地固系计算比导航系计算综合效率也提高约30％.针对地固系正常重力计算复杂的问题,提出了地固系正常重力间接算法,首先把三维直角坐标转换为大地坐标,然后在导航系计算正常重力,最后把正常重力转换到地固系,姿态矩阵也在此过程计算完毕.试验结果表明,间接法与传统方法计算精度相当,但间接法计算过程简单直观,计算效率可以提高3％ ～ 26％,尤其对于高速高精度的捷联惯导,效率提高尤为明显,具有重要的工程实用价值.     

航空发动机核心机建模技术研究与应用

采用部件建模技术,建立了某型发动机核心机实时动态数学模型。应用该模型进行计算时,不必采用迭代的方式,而是按时间步长1次完成,提高了计算的实时性;有效保持了发动机特性的非线性特征,有利于提高模型的计算精度。     

弹道跟踪数据融合处理的快速算法

针对弹道跟踪数据融合处理中的大计算量环节研究了快速算法。用样条函数表示弹道参数,建立了多测元的联合观测模型和弹道参数的非线性融合计算模型,给出了弹道参数的求解算法,分析了弹道参数融合计算中的大型矩阵运算问题,利用基础线性代数函数库提高了大型矩阵的运算速度。建立了样条模型计算的非线性约束优化模型,给出了确定样条节点位置的优化算法,通过分析样条模型的计算原理设计了并行算法,实现了样条模型的并行化计算。仿真结果表明,弹道参数融合计算和样条模型计算的效率都得到了显著提高,计算时间减少了65.47%,对缩短数据处理周期有重要意义。     

基于非结构网格流场计算的网格重排序

提出基于非结构网格流场计算时的网格重排序,优化了网格信息在内存中的存贮位置,使得计算过程中调用数据的处理速度加快,有效提高了计算效率。以M6机翼及DLR-F4翼身组合体的跨音速无粘流场为算例,分别通过显式四步龙格-库塔推进格式和隐式LU-SGS推进格式来验证网格重排序后对计算效率的提高作用。原始网格均用Delaunay方法生长,并分别使用了两种重排序方法进行比较。计算结果表明:重排序后的计算效率可以提高到未排序时的3到5倍,很具实用价值。     

网格计算技术在大规模CFD计算上的应用

针对大规模数值计算规模化管理的难题,采用网格计算技术将网络环境下分散的计算资源组建了CFD网格计算环境.采用主动接入方式,以中心服务器为核心的网格计算管理和服务模式,构建CFD网格计算环境,实现了CFD计算在网格环境内的自动化计算管理,提高了计算资源的利用效率,降低了设计人员的劳动强度.     

跨声风扇周向畸变流动的谐波平衡法计算

研究了采用"相位延迟"边界条件的谐波平衡法用于周向进气畸变计算的可行性及其计算效率,考虑了进气总压畸变形式为正弦波和方波的两种情况。对于正弦波形式的进气总压畸变采用较少阶数的谐波计算就可得到和双时间推进法相近的计算结果,计算速度最大可以提高48倍。对于方波形式的进气畸变,则需要采用更多阶数的谐波计算才能还原得到较准确的非定常流场,计算速度提高了7倍。算例表明,谐波平衡法可以有效的应用于跨声风扇周向进气畸变的非定常计算中,同时大幅提高计算速度。     

微引射防冰腔热力计算

介绍微引射器的防冰腔计算方法,建立了微引射器及双蒙皮防冰通道流体动力计算及热力计算的数学模型,并说明了用这些数学公式计算的流程,用所开发的计算机程序对运七机翼微引射防冰腔进行了实际计算,并探讨了提高计算精度及节省计算机机时的方法。     

基于eN转捩预测方法的增升装置失速特性数值模拟研究

基于三维可压缩Navier-Stokes方程耦合e~N转捩预测方法,对30P30N二维多段翼型和DLR-F11三维增升装置进行了数值模拟研究。计算结果表明,e~N方法计算得到的转捩位置与试验值符合较好,增升装置翼面特别是主翼和襟翼下翼面存在大范围层流流动,为提高增升装置绕流计算准度,计算中考虑转捩的影响是非常有必要的。相比全湍计算,e~N转捩计算升力系数与试验值符合更好,提高了增升装置大迎角失速特性计算准度。     

考虑端区黏性和径向掺混的流线曲率通流算法研究

高负荷多级压气机的技术难点之一在于如何快速、完整、准确地获得全工况下叶片前尾缘气动参数分布和整机特性结果,而通流计算凭借速度快、指向性明确的优点,一直是多级压气机设计的核心技术。为了提高通流代码的适用范围和预测精度,推导了统一适用于轴流、斜流、离心压气机的流线曲率方程,引入端区黏性和径向掺混模型,应用Fortran语言针对轴流压气机全新开发了其通流特性计算程序,对NASA Rotor37、Rotor67转叶和GY1-2J压气机进行了计算,计算结果表明:新模型的加入提高了程序的计算精度,初步验证了程序的有效性。并提出提高程序的通用性是今后通流计算的发展方向。     

气动热数值模拟中壁面热流后处理方法

针对气动热数值模拟中壁面热流后处理方法进行了系统的数值研究,提出了一种有别于常规微分法和积分法的壁面热流后处理方法.通过3种算例进行了对比分析,结果表明:即使在流场参数相同情况下,采用不同的壁面热流后处理方法计算的热流结果也会有较大差异,且表现出不同的网格依赖性.微分法对网格的依赖性较大,单纯提高差分精度不一定能够减弱网格依赖性和提高计算精度,还与近壁处网格分布密切相关;积分法从能量平衡积分方程出发,极大地减弱了网格依赖性,具有较高的计算精度.所提出方法具有与积分法相当的网格依赖性和计算精度,尤其在驻点区域,其网格依赖性更弱于积分法,且计算简单,提高了计算效率.      

微带天线的多频特性

引言 微带天线由于具有一系列独特的优点,正在受到人们越来越多的重视,可以预计,在不远的将来,微带天线将在很广泛的应用场合取代人们已应用了多年的一大类天线。 微带天线固有的一个重要特性——多频特性,即一个微带天线可以在离散的多个频率点上工作的特性,至今尚未被人们所充分认识、研究和利用。本文初步研究了微带天线的多频特性（对VSWR而言）,包括多频的观测、实验和理论分析。     

高功率微波极化方式与入射方向对微带线耦合特性影响分析

针对高功率微波（HPM）耦合效应分析需求,采用时域多分辨分析法（MRTD）模型,给出了微带线端口电压和电流方程,验证了模型的适用性。基于MRTD,仿真分析了HPM作为入射波时其极化方式、入射方向对微带线电磁耦合特性的影响。结果显示：垂直线极化波和平行线极化波耦合系数峰值出现在不同的频率,且峰值相差约15 dB,采用圆极化波时耦合系数最大;入射方向变化时,耦合系数以水平入射最小,垂直入射时最大,高频段微带线边缘有绕射现象产生。该分析采用了分析计算电大尺寸HPM耦合效应的MRTD数学模型,提高了HPM与复杂结构的耦合计算效率。     

飞行器热天线热电联合计算方法

为了研究高温环境对热天线电性能的影响,提出了一种热电联合计算方法,通过三维热、电信息传递程序实现不同温度区域介电性能的空间加载,以此关联温度场与电磁场网格节点信息,实现温度场-电磁场耦合问题的联合计算。介绍了热电联合计算方法的原理与思路,并以某型热天线为研究对象,考查天线温度分布对电性能的影响程度,为飞行器热天线的精细化设计提供了途径。该方法可应用于热透波结构的高温电性能分析,尤其对于恶劣热环境下的热天线/天线罩设计具有较高的应用价值。     

结构可靠度逐步逼近径向基神经网络响应面法

提出了逐步逼近径向基神经网络响应面法计算结构可靠度,这种数值方法较好地解决了结构功能函数非线性及结构随机变量非正态分布时采用结构可靠度指标度量结构可靠度存在误差的问题.经过多个数值试验表明,该算法迭代收敛迅速、计算准确性较高、应用过程简单.利用这种方法对含有多个随机变量的结构进行可靠性分析和设计,特别是当结构功能函数具有较强的非线性或结构随机变量分布具有较大的偏斜度时,具有一定的应用价值.      

远场假设对喷流噪声预测中格林函数求解的影响

目前基于雷诺平均Navier-Stokes（RANS）的喷流噪声预测方法在格林函数求解时, 为简化求解过程,通常对喷流流动做平行流假设, 对观测点做远场假设。随着格林函数求解方法发展,近年来的研究表明平行流假设对下游观测点格林函数的计算会引起较大偏差,而目前远场假设对格林函数求解的影响仍不清楚。为研究远场假设对喷流格林函数求解的影响,以二维喷流为例, 采用计算气动声学方法（CAA）分别数值求解了观测点远场假设条件与实际条件下90°~150°方向喷流内伴随格林函数,进而分析远场假设对格林函数求解的影响。研究结果表明,对于不同方向的观测点,由观测远场假设导致的伴随格林函数求解偏差不尽相同,且对于越靠近喷流中心线方向的观测点,远场假设导致的偏差越大,其中150°方向观测点,采用远场假设后,格林函数计算结果最大偏差达到-15 dB以上。因此,对于靠近喷流中心线方向的噪声观测点而言,为避免预测偏差,应采用实际观测条件求解喷流格林函数。     

大气紊流环境下的纵向飞行参数响应分析

介绍了大气紊流环境下飞机纵向飞行参数的计算与匹配方法.由于紊流的随机性特征,飞行参数响应分析采用时域与频域相结合的方式,对某型飞机紊流环境下的纵向飞行参数进行了计算分析.采用所介绍的方法,可以获得比较合理的飞行情况来指导型号设计工作.     

基于AFAS 软件的机翼壁板耐久性分析及试验验证

随着航空工业的发展,对机翼壁板耐久性分析的精确性及计算速度的要求日益提高,必须考虑从手工计算到程序计算的转变。飞机结构耐久性分析软件系统(AFAS)是为了满足飞机结构耐久性分析的工程实际使用要求而编制的。本文介绍采用该软件中包含的DFR方法计算程序对机翼壁板结构进行耐久性分析的过程,并进行试验验证,分析得到的机翼壁板结构具有95%置信度、95%可靠度的可靠性寿命与试验统计结果非常相近,证明AFAS是一个可以进行耐久性分析的快速准确的计算工具。     

冷轧带钢轧制压力分布计算

利用经典轧制力模型和轧辊压扁模型，考虑二者的耦合推导出考虑轧辊压扁半径的轧制力显示计算公式，在此基础上开发了基于影响函数法的冷轧带钢轧制压力分布计算程序，并用该程序结合实际生产中的采样数据模拟计算了HC冷轧机各道次的轧制压力分布和总轧制力，将所得计算结果与现场实测数据进行比较。结果表明：所得轧制力计算结果与现场实测数据平均值相近，轧制压力分布与实际相符，为HC轧机板形控制提供了一种计算轧制压力分布的有效方法。     

Solving a radiation problem with forward speed using a lifting surface method with a Green''s function

The progress of seakeeping computations requires development of computating codes for unsteady flows around a ship or its elements. In this paper, we present a method of calculation concerning waves radiated by an oscillating surface-piercing flat plate with forward speed, with a yaw angle. By use of Green's third identity, the problem is transformed into the resolution of a Fredholm integral equation of the first kind by a panel method using Green's function. The Green's diffraction-radiation function with forward speed is used. Its numerical values are calculated by an adaptative integration procedure to reduce the computation time. The present method permits determination of the pressure jump distribution across the plate, the total forces and moments. The results obtained are compared with other numerical methods in hydrodynamics and in aerodynamics, and with experimental data obtained in a water tank.     

一种新的玫瑰扫描目标跟踪边界重心法

针对自动寻的导弹玫瑰扫描目标跟踪方法中由于扫描过程非线性引起的目标采样非线性,导致求取目标重心算法的运算量和计算误差大等问题,提出了一种基于格林公式的、新的求取目标重心的边界重心法,该方法仅使用目标边界上的采样点。通过对运算量和计算误差的理论分析,与现有的权重重心法相比,本方法可以有效解决采样不均匀的问题,同时运算量减少了90%,误差仅为权重重心法的下界。另外,本方法也可用于基于图像信息的凝视目标跟踪等类似算法中来提高跟踪准确性及运算速度。