飞机在变化风场中三维运动的数值仿真

本文在引入和定义了必要的坐标系和Euler角以后,推导了计及任意变化风场的飞机一般运动方程.所得的两组运动方程均适宜于在大气紊流和风切变中飞机三维飞行运动的数值仿真。还讨论了风梯度的气动力效应。     

四元数在飞行数值仿真中的应用

为了在飞行器运动的计算和仿真试验中避免出现奇异情况及减少三角函数计算,最好在飞行器运动方程中采用四元数来代替欧拉(Euler)角。本文介绍了四元数的性质,阐述了四元数在坐标变换中的作用和性质,并按照我国通用的坐标系和角度规定,建立了采用四元数的飞行器运动方程。最后评价了在飞行器运动方程中欧拉角和四元数的优缺点。     

旋翼桨叶绕流Euler方程数值模拟

本研究应用格心格式有限体积法求解欧拉方程,模拟绕旋翼浆叶的流动。旋翼尾流的作用通过自由尾流分析来求解局部的诱导下洗速度以修正翼型攻角。选择O—H型式生成弦向和展向网格。通过算例计算为进一步寻求有效的Euler或NS方程数值解积累了经验。     

隐式无网格算法及其应用研究

本文的主要目的在于研究求解Euler方程的隐式无网格算法，并应用于复杂的流场计算。采用无网格算法，计算区域用点云离散代替通常的网格划分；计算点上的空间导数，用当地点云上引入的二次极小曲面逼近。求解的Euler方程用隐式时间后差离散，结合用Roe的近似Riemann解确定通量，并用LU－SGS算法分步计算，数值求解了单翼型或双翼型模拟的复杂绕流。     

径向简谐激励下航空发动机机匣的振动分析

根据航空发动机机匣薄壁结构的几何特征,将机匣简化为一端简支一端固支的圆柱壳结构,建立了系统的动力学方程,并采用多尺度法对动力学方程进行了求解,进而研究了径向简谐激励下机匣动力学方程的零解及非零解的稳定性,并讨论了系统单模态运动发展到复合模态运动以及能量在模态之间的转移过程,最后采用数值方法研究了径向简谐激励下机匣的振动特性.结果表明:随着径向简谐激励的增加,系统表现为周期运动、倍周期运动及混沌运动交替出现.      

全隐式TVD格式求解二维Euler方程的对角化算法

本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。     

叶轮机械三维不可压Euler方程的直接解法

本文用近似因式分解交替方向隐式格式直接求解叶轮机械中原参数三维不可压Euler方程组。在连续方程中引入“拟压缩性”之后,将原来类型不确定的控制方程变成双曲型方程,给定初边值之后用时间推进法求得稳定解,对一个单级压气机转子内部流场进行了计算,结果与实验值比较符合。     

考虑辐射效应的多孔材料传热有限元算法

应用有限元法求解一维半透明多孔介质辐射对流传热问题。通过对辐射传递方程与边界条件中积分项的离散 ,推导出有限元求解列式。对于多孔介质传热控制方程和气体传热控制方程的求解 ,采用时间积分技术对两相温度统一进行求解。而对于边界方程处理 ,则通过传热方程与边界方程的循环求解 ,求出任意时刻温度场的瞬态解。给出数值算例 ,得出了放置在通风口、同时受到伴随辐射作用的半透明多孔介质传热的瞬态解 ,讨论了部分参数对瞬态温度场以及换热效果的影响      

用CSCM格式求解激波反射问题

本文对Euler方程采用CSCM格式,数值求解了激波反射问题,与目前流行的交替方向隐式(ADI)格式以及总变差递减(TVD)格式的计算结果进行了比较,并作了简要分析讨论。     

对求解纳维尔-斯托克斯方程的新数值方法的一些看法

已经研究了一个新的解纳维尔-斯托克斯方程的显式一隐式法。该方法在空间和时间上有二阶精度,保持守恒形式,且复杂性比其它隐式方法要小得多。因为它不要求块或者纯量三对角反演。该方法被用于求解复杂的二维定常状态超音速流问题。在高古朗一弗里德里希斯一卢伊(CFL)数的情况下对新方法的计算效率以及从该方法得到的解的精度进行讨论。讨论了一些改进了的型式且对该方法提出一些意见。这些在成功地使用该方法方面可能是有帮助的。     

损伤定位中频率法的改进及应用

传统的基于频率的损伤定位指标是通过忽略结构运动摄动方程的二阶项得到的,用于结构的损伤位置识别可能导致误判.在对传统的基于频率的结构损伤位置识别机理进行了阐述基础上,理论推导改进的损伤定位指标.最后以复合材料悬臂梁进行数值模拟试验.试验结果表明,所建议的改进方法指标比原来方法指标更为有效.     

积分方程法计算翼型的跨音速绕流

从跨音速小扰动方程出发推导积分方程的过程中,本文用任意形状的封闭曲线C_Q（其极限趋于零）挖去奇点Q,最终得到无奇性（指无穷奇性,不包括Cauchy奇性）的积分方程。 对于跨音速流中的圆头翼型的前缘问题,提出了一种解决办法。 证明了Nixon给出的反演公式对于超临界有激波的小扰动流动也成立。 关于积分方程法中的人工粘性方法,对Sachdev和Lobo提出的方法做了改进。 最后给出了NACA0012翼型在有无升力和有无激波各种状态下的计算结果。比较表明,本方法的计算结果与其它方法的计算结果符合得较好,且计算量很小。     

小型无人直升机动力装置活塞发动机的动态效应

本文介绍了二冲程活塞发动机进、排气系统的动态效应的机理,分析了影响动态效应的主要因素并给出了它们之间的关系式,指出了在选择活塞发动机时,必须选装厂家提供的进、排气系统,说明了活塞发动机作为小型无人直升机的动力,其进、排气系统动态效应得到充分利用的原因。     

Real-time trajectory optimization for powered planetary landings based on analytical shooting equations

Traditionally, numerical trajectory integration for shooting equation calculation and iterations for shooting with randomly guessed initial solutions deteriorate the real-time performance of indirect methods for on-board applications. In this study, the indirect method is improved to achieve real-time trajectory optimization of fuel-optimal powered planetary landings with the help of analytical shooting equation derivations and a practical homotopy technique. Specifically, the contributions of t...     

改进整体传递矩阵法计算复杂转子系统临界转速

提出了一种基于整体传递矩阵法的改进传递矩阵法。新方法将Riccati变换及奇点消除方法应用于整体传递矩阵方程,导出了耦合单元和非耦合单元的Riccati传递矩阵计算公式。这种方法克服了传统传递矩阵法的缺点,易于开发通用计算程序。文末给出了数值算例,计算结果表明了该方法的有效性和正确性。      

复杂空间四连杆传动系统的一种算法

1.引言 空间四连杆机构的运动关系是一个空间几何关系。确定四连杆相互位置时,存在着一球面与圆弧相交或两圆弧相交等情况,它们的运动与几何关系都以一些二次联立方程式表达。求解与根的判定都比较复杂,有时还要用迭代方法。本文所介绍的方法采用空间坐标变换,将空间四连杆机构的运动点转换到特定的平面内,利用平面三角形的几何关系求解。所用的方法都是典型的,简便易行,特别适用于计算机计算。     

石英挠性加速度计标度因数受气压影响分析

石英挠性加速度计是惯导产品的核心器件之一,主要敏感载体运动中的加速度,为系统提供位置信息.加速度计具有力学敏感特性,其测量精度易受环境因素变化的影响.针对加速度计在不同地域测试时标度因数发生较大变化的问题,分离了加速度计异地测试过程中的时漂、温变、形变等多个环境因子后,提出了造成加速度计异地标度因数改变的主要原因是环境气压发生了较大变化.根据石英挠性加速度计的基本原理,结合加速度计检测质量摆的具体受力情况,对气压变化造成加速度计标度因数改变的机理进行了详细分析,给出了标度因数随气压变化的相对变化理论计算公式,经过验证试验,实测情况与理论分析结果一致,可以为石英挠性加速度计的优化改进提供参考依据,具有一定的工程应用价值.     

一种大椭圆轨道卫星星历预报方法

针对大偏心率轨道长时间星历预报,稠密输出星历的低效率问题,提出了一种新的星历预报方法。该方法通过建立卫星摄动力模型,对卫星运动方程进行数值积分来获取卫星星历。该方法的轨道积分器以标准的8阶Adams-Cowell多步法为基础,针对大偏心率轨道的特性,在一个轨道周期的不同时段内采用不同的积分步长,同时引入了用于生成小间隔等间距卫星星历的插值公式。该方法降低了卫星运动方程右函数的计算次数,尤其适用于需要稠密输出卫星星历的情况。以STK(Satellite Tool Kit,卫星工具包)的HPOP(High Precision Orbit Propagator,高精度轨道预报)模型为验证基准,通过数学仿真校验了该方法的有效性和精度。算例表明,该方法在预报精度损失很小的情况下,计算时间远小于标准的多步法和HPOP模型。     

基于Voronoi图和蚁群算法的无人机航迹规划

为提高无人机任务环境模拟的真实性,利用改进后的Voronoi图对任务环境进行建模.同时,为了更快地生成一条满足任务需求的最优飞行航迹,提高航迹规划的实战性和高效性,分析了蚁群航迹规划算法的运行原理,以及算法运行机制对算法性能的影响,提出了算法的改进原则,并在此基础上给出了新的信息素更新方式和新的启发式.利用改进后的蚁群...     

The Least-Square Approximation of Inertial Platform Drift

The drift angle of an inertial platform which is gyro-stabilized with respect to inertial space is equal to the integral of the gyro drift rate. Under controlled laboratory environment the drift angle, denoted by y, may be measured and plotted against time in an interval [0, T]. Without loss in generality, one may take y(0)= 0. A straight line yf can be found, such that the quantity E2 is minimized, where begin{equation*}E^2 = {frac{1}{T}int^{T}_{0}(y-y_f)^{2}}dt.end{equation*} The equation for yf is of the form yf = at + b and, in general, both a and b are nonzero. It is desirable to determine a statistical relationship between the gyro drift rate and the expected value of the minimum E2 for any given interval T. An analysis in this paper determines this relationship and derives a general expression for , where the symbol <*> denotes statistical expectation. It is found that increases linearly with the variance of the gyro drift rate. This general formula is then developed in detail for the case of a first-order Markovian gyro drift. is evaluated numerically and its square root plotted vs. the interval T and the gyro correlation time. The same problem is also solved for the case when yf is constrained to intersect the origin, i.e., when b=0.