模态集结法及其在飞行力学中的应用

 本文给出一种处理线性系统降阶的集结方法,并根据该方法将飞机纵向运动方程和横侧向运动方程进行降阶分解,得到了飞机的沉浮、短周期、滚转-螺旋和荷兰滚运动模态的各近似式,结果表明本文中的各近似式精度较好,特别是提高了沉浮及荷兰滚阻尼的估计精度。     

一种求解透平叶栅三维流场的高精度TVD格式

本文提出了一种应用三阶精度Chakravarthy-OsherTVD格式求解环形透平叶栅三维定常流场的隐式近似因子分解方法。控制方程为柱座标速度分量表示的任意曲线座标系下的三维欧拉方程。采用Beam-Warming近似分解, 对隐式项进行迎风差分, 得到了时间精度为一阶、空间精度为5点三阶的隐式格式。算例对5点三阶和5点二阶精度的格式进行了比较。表明空间三阶精度隐格式对计算结果没有明显改进, 但收敛残差有较大降低。      

显式制导下定位定向误差的影响分析

建立了定位定向误差对显式制导精度影响的解析模型,提出了利用标准弹道参数估算定位定向偏差对制导精度影响的近似解析方法.结果表明,该方法计算精度高、速度快,可以满足显式制导下定位定向误差影响分析的需要.     

p型多重网格间断Galekin有限元方法研究

在二维非结构网格上,使用p型多重网格间断Galerkin方法求解定常可压缩欧拉方程。p型多重网格方法主要特征是通过对不同阶次多项式的近似解进行递归迭代求解。文中高阶近似(p0)上使用显式格式,在最低阶近似(p=0)上选用隐式格式,而非显式格式,从而在保证精度和占用较小内存的情况下加速收敛到定常解。运用该方法对NACA0012跨音速无粘流动进行数值模拟,数值结果表明:p型多重网格方法同单重显式Runge-Kutta方法相比,收敛速度能够提高6倍左右,并且精度保持不变。     

依据试验数据求取航空发动机部件特性的新方法

介绍了通用的依据已知发动机部件特性求取所需部件特性的耦合法,并指出该方法的不足。以压气机部件特性为例提出求取部件特性的新方法。该方法利用发动机试验数据得到转速线与共同工作线的交点,把转速线近似为指数函数曲线,由PSO算法结合所得的交点数据拟合出每条转速曲线。以试验数据为依据,保证了特性的真实和有效性,而PSO算法则提高了特性线拟合的精度。计算结果表明所提出的方法比传统的耦合法更为准确,有利与提高发动机工作过程仿真精度。     

升力式再人弹道拉起段近似解研究

针对现有方法在求解升力式再入弹道拉起段时的不足,考虑再入平面动力学模型中的引力项、离心力项的影响,推导了拉起段的速度、高度、航程变化量的近似解.通过与数值仿真校验和现有方法的对比,验证了近似解在估计精度上的改善效果.研究结果可用于升力式再入飞行器的设计评估与制导算法的设计.     

考虑真实气体效应的一维非定常平衡流的随机选取法计算

本文提出一种近似随机选取法,计算一维非定常的高温平衡流动。采用稀疏激波替代等熵稀疏波,选用拟合经验方程计算高温热力学性质,并利用切割迭代方法加快收敛,从而计算大为简化,并有很好精度。     

弹道式地地导弹精度计算公式及其拟合方法

归纳和总结厂弹道式地地导弹的精度变化规律,根据这些规律和谐波分析理论建立了导弹精度计算公式。文中根据函数逼近理沦,提出了确定谐波函数阶次的方法和具体的拟合步骤。应用这种方法导出了实际导弹的精度计算拟合公式.采用这些公式和采用弹道计算方法分别计算任意相同条件下的导弹精度值.其相对误差小于1%。本文提出的精度计算拟合公式是简便实用的.     

不可压缩流动的高精度非标准格子玻尔兹曼方法

首先回顾了高精度非标准格子玻尔兹曼方法的发展历程,基于高精度通量重构格式,发展了一种通量重构格子玻尔兹曼方法（FRLBM）。采用两种求解方法：一种是将碰撞项隐式处理,直接求解离散速度玻尔兹曼方程（直接法）;另一种是先执行碰撞步,再求解纯对流方程（分步法）。通过收敛性研究,比较了这两种方法的精度和稳定性。研究结果表明,在小时间步长下两种方法误差近似,都能取得高阶精度;然而当时间步长增大,直接法误差几乎不变,分步法误差出现明显上升。由此表明,当取得近似误差时,直接法可以采用较大时间步长,计算效率更高,且直接法的稳定性略占优。接着通过模拟顶盖驱动方腔流验证了FRLBM捕捉流场细节的能力,并且比较了基于半隐格式的显式方法和一阶、二阶及三阶隐式-显式Runge-Kutta格式的时间离散,在不同雷诺数下的最大允许Courant-Friedrichs-Lewy数,数值结果表明二阶隐式-显式Runge-Kutta格式效果最优。最后数值模拟了圆柱绕流,验证了FRLBM计算复杂外形绕流的可靠性。     

弹道导弹误差传播特性

在采用式(1.6)的引力模型(考虑了J_2项)和发射条件影响下,本文推导了惯性导航系统误差传播转移矩阵的近似解析解,而文献[1]给出的解是它的一种简化式。利用该近似解,进一步用三角多项式去拟合干扰量,就得到弹道导弹的误差传播特性。这套公式的实用价值在于:利用有限的观测数据就可以估算或预测弹道导弹的误差传播及其落点偏差。     

一种部分森林覆盖环境电波传播特性预测方法

针对部分森林覆盖环境中传播媒介的非均匀性,提出一种基于动态网格的抛物方程(PE)预测方法,该方法采用宽角PE法进行大区域电波传播特性预测。为平衡PE求解的准确性与时效性,引入动态网格技术。通过仿真算例,探讨了基于动态网格的PE法在部分森林覆盖环境电波传播特性预测中的有效性。结果表明:在计算精度相似的情况下,相对于均匀细网格,动态网格会大大减少PE所需的计算时间;而相比均匀粗网格,动态网格具有更高的计算精度。     

一种4站情况下基于TDOA/FDOA的无源定位方法

在基于时频差的三维运动目标无源定位系统中,针对在4个接收站的情况下搜索法实时性低的问题,提出了一种基于改进的加权最小二乘法（MWLS）与萤火虫算法（FA）相结合的无源定位方法（MWLS-FA）。该方法的第1步通过构造一组新的方程来对加权最小二乘（WLS）方法进行改进,使得改进后的WLS方法在4站情况下也能得到目标位置和速度的初始值,第2步利用这个初始值为FA方法提供一个动态的搜索区域,同时在约束条件的添加和参数选择两个方面针对性地对FA方法做出了调整和改进。仿真结果表明,该方法在4站情况下对目标的定位精度可以达到克拉美罗下限（CRLB）,而且在实时性和抗噪性方面优于传统的搜索法,同时该方法在5站情况下的抗噪性能优于两步加权最小二乘法（TSWLS）和约束加权最小二乘（CWLS）法。     

复合材料舱体零件数控定位技术

文摘复合材料壳体采用手工划线定位零件的装配方式,精度低、重复性差,对工人的操作水平要求较高,而且利用金属划针会造成复合材料壳体表面损伤,为改变手工划线定位现状,提出利用数控设备结合喷涂设备完成零件位置线标记,为保证零件位置轮廓准确性,开展等距喷涂技术研究,构建数控定位轨迹方程,通过数控系统硬件连接组建零件定位系统,通过数控设备结合喷码装置实现复合材料壳体表面零件轮廓位置的喷涂标记,零件定位精度可以满足±0.1 mm的要求,实现了数控定位代替手工操作方式,提升了复合材料壳体零件定位精度和一致性。     

A second-order numerical method for elliptic equations with singular sources using local flter

The presence of Dirac delta function in differential equation can lead to a discontinuity,which may degrade the accuracy of related numerical methods.To improve the accuracy,a secondorder numerical method for elliptic equations with singular sources is introduced by employing a local kernel flter.In this method,the discontinuous equation is convoluted with the kernel function to obtain a more regular one.Then the original equation is replaced by this fltered equation around the singular points,to obtain discrete numerical form.The unchanged equations at the other points are discretized by using a central difference scheme.1D and 2D examples are carried out to validate the correctness and accuracy of the present method.The results show that a second-order of accuracy can be obtained in the fltering framework with an appropriate integration rule.Furthermore,the present method does not need any jump condition,and also has extremely simple form that can be easily extended to high dimensional cases and complex geometry.     

基于最大似然估计法的加速度计误差参数标定方法

为降低捷联惯导系统误差参数标定过程对高精度转台的要求,提出一种基于最大似然估计的模观测标定方法。在考虑噪声影响的情况下建立加速度计的误差参数模型,运用最大似然估计法对误差参数进行标定计算,并将标定结果的精度与所考虑参数模型的克拉美限相比较。结果表明,与传统标定方法相比,该方法标定精度相当,降低了对标定转台的要求,减少了标定时间,有较高的工程应用价值。     

An autonomous celestial navigation method for LEO satellite based on unscented Kalman filter and information fusion

A reliable and secure navigation system and assured autonomous capability of satellite are in high demand in case of emergencies in space. Celestial navigation is a fully autonomous navigation method for satellite. To near earth satellite, the earth direction is the most important measurement and the horizon sensing accuracy is the most important factor which effects celestial navigation accuracy. According to the mode of acquiring horizon measurement, satellite celestial navigation methods can be broadly classified into two approaches: directly sensing horizon using earth sensor and indirectly sensing horizon by observation of starlight atmospheric refraction. For these two methods are complementary to each other, a new Unscented Kalman Filter (UKF) based information fusion method is proposed here for hybriding them. Compared to the traditional celestial navigation method, this method can provide better navigation performance and higher reliability. The hardware-in-loop test results demonstrate the feasibility and effectiveness of this method, in most cases the accuracy is sufficient for near earth civilian satellite and moreover it can be used as a backup system to provide redundancy.     

高升力系统外形的数值模拟计算

分析和讨论了高升力系统外形无黏和黏性流动的物理特征,表明了流动现象的复杂性和求解的艰难性。深入研究了二维和三维高升力外形空气动力数值模拟计算的雷诺平均N-S(RANS)方程方法、计算结果、以及它们和实验数据的比较。二维结果表明压强分布、速度型等计算值和实验值吻合很好,但最大升力系数两者仍有差异。三维计算最大升力系数的结果与二维类似。还分析了影响计算结果的各种因素,包括计算网格、湍流模型、外形的几何模拟的逼真度等。特别指出了转捩模型研究对提高计算准确度的重要性。文中解释了研究高升力系统外型的数值模拟方法,以及同时发展能估算C_Lmax的工程半经验方法的重要性和迫切性。     

双方程 k-ω SST 湍流模型的显式耦合求解及其在叶轮机械中的应用简

 双方程k-ω剪切应力输运（SST）湍流模型通常以隐式耦合方式或者显式半耦合/解耦的方式来求解。本文提出了该模型的一种显式耦合应用方法,即通过点隐的方式来处理湍流源项的刚性,并与混合Runge-Kutta时间推进以及当地时间步长、隐式残差光顺等加速收敛技术相结合,从而使得湍流方程可以与流动方程同时求解。为了增强计算的鲁棒性,进一步对湍流变量进行了限制。将所发展的方法用于DLR平面叶栅算例,确认了求解结果的正确性以及刚性的来源。通过对三维NASA Rotor 67的模拟,验证了SST模型的精度;进一步将其与Badwin-Lomax（BL）模型、Spalart-Allmaras（SA）模型对比,发现三者都能正确地捕捉出口参数分布,且SST与SA模型的模拟结果比较一致;对于该算例,SST模型在总温模拟上更具优势,而BL模型在总压分布上与试验值更加接近。     

加速退化模型及外推结果准确度的定量验证方法

加速退化试验本质上是牺牲可靠度评估精度换取可靠度评估效率,外推出的可靠度结果通常会与真实值存在一定的偏差,因此,验证加速退化模型与外推结果的准确度是非常必要的。根据可靠性建模的步骤设计了验证加速退化模型及外推结果的技术流程,结合Wiener-Arrhenius加速退化模型构建了验证技术框架。采用了基于假设检验的模型验证方法,利用Kolmogorov-Smirnov检验法、Anderson-Darling检验法分别验证Wiener退化模型、加速退化模型及外推的可靠度模型是否合理。提出了具有工程实用性的基于面积比的外推结果验证方法,利用蒙特卡罗仿真解决复杂可靠度函数积分问题,用于定量表征外推结果的准确度,通过面积比阈值判定是否接受外推结果。通过惯导系统伺服电路实例应用说明了技术框架的可行性与有效性,研究工作为解决加速退化试验中的验证难题做出了有益的探索,所提出的各种定量验证方法具有较好的工程应用价值。     

某微型发动机振动特性计算分析

研究了具有锥齿轮啮合及机匣系统的某微型发动机的振动特性。采用弯扭耦合传递矩阵法计算临界转速、不平衡响应和初始弯曲响应,计及转动惯量、陀螺力矩、剪切变形、粘性阻尼和轴向力的影响。为了提高机匣振动固有频率的计算精度,还采用轴质量均布的传递矩阵法计算转子－机匣系统的临界转速与不平衡响应。临界转速及不平衡响应计算结果与试验结果吻合较好。