一种基于轨道要素形式终端约束的航天器空间变轨迭代制导算法

针对航天器空间变轨任务的制导问题,研究了一种适用的迭代制导算法。在传统迭代制导方法的基础上,直接在地心惯性系中建立最优控制模型,以推力方向矢量为控制量来适应大姿态角变化情形;推导直接以目标轨道要素为终端约束的边界条件,给出终端约束方程求解精度和入轨精度的关系;得到一种简单有效的基于轨道要素形式终端约束的航天器空间变轨迭代制导算法。通过仿真验证了所给制导算法的有效性,相比传统迭代制导方法其具有更强的适应性。     

电弧喷射推力器流动区域电磁场数值解法

为了有效地进行电弧喷射推力器流动区域全场(流场、电场、磁场和化学反应)耦合数值模拟并揭示其电磁场特征，建立了适当的电磁场模型，对电磁场控制方程的数值解法进行了详细研究。模型基于推力器稳态工作特征，控制方程由麦克斯韦方程组简化得到，采用有限控制容积积分方法离散椭圆型控制方程，采用9种不同的迭代方法求解离散方程。给出了不同迭代方法的收敛速度、数值稳定性和最终精度。研究表明，Gauss—Seidel逐线超松弛迭代法是推力器流动区域电磁场离散方程的一种快速有效的数值解法。     

飞机推出和跑道掉头滑行行为计算方法

针对飞机推出和跑道掉头两种地面滑行特殊行为的位置和姿态计算问题进行了研究。从机场塔台管制运行模拟仿真的角度出发,提出了非完整运动学约束下飞机推出滑行和跑道掉头运动计算方法。依据飞机推出实际和非完整运动路径跟踪模型,设计出适合仿真推演计算的飞机推出运动求解方法。根据飞机跑道掉头的几何约束和转弯角运动约束,设计出开环控制下完整刻画飞机运动的求解方法。与动力学求解模型相比,该方法具有计算简便,适合于视景仿真驱动的优点,并通过数学仿真验证了方法的有效性。     

一种计算非定常二维流动的无网格算法

主要目的是发展一套求解非定常流动的无网格算法。计算区域的离散方面,提出了一种按区域进行填充布点的点云自动生成方法;发展了一种点云的运动技术来实现离散点云对物面边界的随体运动;在点云离散的基础上,采用最小二乘法求解矛盾方程的方法来求取空间导数,进而获得数值通量;采用双时间方法进行时间离散推进,其中物理时间迭代采用二阶隐式格式,伪时间迭代采用四步龙格一库塔显式格式,为了加速收敛,在伪时间迭代中采用了当地时间步长和隐式残值光顺等加速收敛措施。最后,利用本文算法模拟了NACA0012翼型和NACA64A010翼型的跨音速非定常流动,并将计算结果与实验测量结果进行了对比分析,验证了上述方法的正确性和实用性。     

无陀螺飞行器姿态和角速度确定

通过预测滤波方法研究了利用矢量观测确定无陀螺飞行姿态和角速度的方法。给出的观测滤波法在一定的二次型优化准则下，对飞行器姿态和名义姿态运动模型的等效角速度误差可同时进行实时估计，能够准确跟踪敏感器的测量，有效减小模型误差的影响，数值仿真结果验证了算法的优越性。     

利用Singer模型的无陀螺姿态和角速度估计

 给出了一种利用Singer跟踪模型的扩展卡尔曼滤波器 (ExtendedKalmanFilter) ,用于无陀螺姿态和姿态角速度估计。在滤波器中 ,姿态和姿态误差分别由姿态四元数和误差四元数表示 ,而姿态角加速度由一阶Markov过程描述 ,从而避免采用姿态动力学模型。利用数值仿真计算验证了滤波器的性能。在所有的仿真过程 ,滤波器显示出快速收敛能力。稳态估计误差主要由测量更新频率和精度决定。     

一种稀疏阵列下的二维DOA估计方法

研究了稀疏阵列下二维波达方向（DOA）的估计问题,提出一种基于不动点迭代的空间谱估计（FPC-MUSIC）算法。首先建立基于矩阵填充的DOA估计信号模型,并验证该信号模型满足零空间性质（NSP）,其次通过不动点迭代算法将稀疏阵列信号恢复为完整信号,最后利用恢复信号估计二维DOA。该算法可在稀疏阵列下大幅度降低谱估计平均副瓣,在大幅度降低阵元数的同时具有较高的估计精度。计算机仿真表明：FPC-MUSIC算法可在稀疏阵列下准确估计二维DOA,验证了该算法的有效性和优越性。     

基于平衡流形展开模型的航空发动机健康参数估计

为提高航空发动机健康参数实时估计的精度,结合非线性卡尔曼滤波器从底层模型方面进行改进。以航空发动机部件 级模型为基础,提取运行数据,采用线性拟合法求解系数矩阵,建立传统状态空间模型;采用BP神经网络拟合调度参数,建立设计 点处包含健康因子的平衡流形展开（EME）模型。基于EME模型分别采用扩展卡尔曼滤波器和无迹卡尔曼滤波器进行多种参数 退化模式下的数值仿真估计。仿真结果表明：得益于EME模型对非线性系统的良好近似性,各参数退化模式下的估计结果准确, 稳态误差不超过3%;与采用部件级模型作为底层模型的方案相比,该方案的估计速度提升了1个量级。验证了基于航空发动机 EME模型结合非线性卡尔曼滤波器进行健康参数估计方法的实际可行性。     

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

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

基于迭代滑模的机械臂控制策略

针对机械臂这类非线性的不确定性系统,基于迭代学习控制与滑模控制策略,提出了一种有效的迭代滑模控制方法。该控制方法通过将滑模控制律引入到迭代学习控制中,并运用Lyapunov理论对控制律进行证明,从而确保系统的稳定性。基于拉格朗日力学法建立动力学模型,得到相对简化的n关节机械臂模型。以一个二关节机械臂为例,通过MATLAB仿真验证所提控制策略可有效提高关节的跟踪速度与跟踪精度,并且在一定程度上可减缓传统滑模控制的抖振现象,与传统迭代学习控制相比,系统具有鲁棒性。物理试验验证了所提控制策略的有效性。     

铰链对含铰结构非线性动力学特性影响分析

针对可展结构中铰链带来的展开后非线性动力学问题,建立了含铰结构的动力学模型。基于铰链侧向和径向的几何约束关系,分析了不同约束条件下的铰链非线性特性。基于谐波平衡法（HBM）,对含铰结构的非线性变量进行一次谐波展开,得到铰链非线性力的谐波展开表达式,将含铰可展结构的非线性动力学方程转化为代数方程。分析了铰链在不同非线性特性下各参数对结构动态特性的影响,得到结构固有频率随铰链间隙、刚度和激振力的变化规律。考虑铰链的非线性特性,通过固有频率脊线求解,得到结构固有频率极值随铰链数量、位置和刚度的变化曲面,并对其进行非线性拟合,得到铰链参数对结构固有频率的影响函数,其计算方法和结果在多维复杂结构中具有一定的扩展性,为可展结构的设计和研究提供参考。利用龙格库塔法对非线性结构进行数值仿真,得到结构的固有频率变化曲线,仿真结果表明了利用HBM进行含铰结构动力学分析的正确性。     

轮履复合式军用地面探测车运动学建模及分析

为更好的对军用探测车运动过程进行控制,提出了轮履复合式探测车运动学模型的构建方法.通过对轮履复合式移动系统结构进行分析,首先建立了探测车三维运动学模型,并推出了其运动方程;其次建立了探测车车体速度与车轮速度之间的雅可比矩阵,利用最小二乘法对车体速度进行求解,提供了获得探测车越障过程中位置和方位的方法.最后对探测车行进过程进行了运动仿真,数值仿真分析验证了运动学模型的正确性及轮履复合式移动系统的运动特性.      

湍流模型离散精度对数值模拟影响的计算分析

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS),通过改变物面法向第一层网格间距,开展了剪切应力输运(SST)两方程模型不同离散精度的数值分析。主要目的是为高阶精度格式在复杂外形上的应用提供技术支撑。计算模型包含了低速NLR 7301两段翼型和高速RAE2822翼型,研究内容主要包括湍流模型的二阶精度离散和五阶精度离散两种方式对收敛历程、边界层湍流黏性系数分布、边界层速度分布、压力系数分布以及气动特性的影响。在与试验数据对比的基础上,计算结果表明:对于不同的第一层物面法向网格间距,湍流模型离散精度对低速绕流计算结果有比较明显的影响,对于高速小迎角附着流动计算结果影响不明显;相对于湍流模型二阶精度离散,湍流模型高阶精度离散网格敏感性较弱,具有更高的数值模拟精度,但收敛性略差。     

一种摇臂式起落架的两种建模及结果对比

为了满足适坠性能要求,一种摇臂式起落架采用了双腔缓冲器。本文详细推导了拉格朗日方程下该起落架的动力学方程,并编程数值计算了四种状态的着陆性能。同时在Adams/Aircraft平台下建立了该起落架模型,通过编写子程序的方法实现了该双腔缓冲器的功能扩展。Adams/Aircraft平台下的仿真结果与数值计算的结果吻合较好,为起落架性能分析提供了很好的相互验证方法。     

流管侧向微缝共振器的声学特性

为研究有流动管道侧向微缝共振器的声学特性,基于离散涡方法建立了一个理论模型.管道安装共振器处利用质量和动量守恒建立间断条件,再结合离散涡模型,计算出共振器的声阻抗和吸声系数.数值模拟表明微缝处流场主要取决于脱落涡层的卷起和拍动,这与流场显示实验结果一致,并验证了半经验模型中微缝剪切层的"盖子"假设.结果还表明共振器声阻...     

基于RBF网络的机械手臂逆问题的求解

机械手臂的研究,要涉及到机械手运动学逆问题,即由已知机械手末端的位姿求出机械手各个关节的关节量。逆问题的解析法存在着无解的问题,数值解法计算量太大。应用RBF神经网络良好的非线性逼近能力,对机械手臂的运动学逆问题进行非线性逼近,利用训练好的网络可以精确,快速地求出逆问题的解,且在多解情况下可以求得最优解。     

磁悬浮斜翼节静动态特性

作为直动式电-机械转换器和二维阀本体之间的连接桥梁，反馈放大机构起着位置反馈、运动转换和推力放大的重要作用。但现有二维伺服比例阀的机械式反馈放大机构都存在摩擦磨损等非线性环节，对阀的静态特性影响明显。基于永磁悬浮技术提出一种新型的无接触式磁悬浮斜翼节，其依靠磁斥力将斜翼动子和阀芯悬浮在中位，同时实现位置反馈和运动转换功能，以彻底去除传统机械式反馈放大机构由于接触带来的摩擦磨损等因素对阀控制特性的不利影响。首先基于电磁场理论建立了磁悬浮斜翼节的数学解析模型，讨论了关键结构参数对磁力矩的影响；随后在Maxwell和Adams平台上分别对其静动态特性进行了基于有限元模拟和多体运动学分析的参数优化，在此基础上制作了实验样机，搭建实验台架研究了斜翼节的静动态特性，从而验证了其工程实用性。仿真和实验结果基本相符，3种验证的斜翼节样机中，当外接比例电磁铁推动外动子移动3 mm时，最大磁力矩可达到0.252 N·m，阶跃响应时间约120 ms。参数化设计表明减小工作气隙和增加斜翼节倾角均能提升斜翼节的磁力矩和动态响应，而在气隙中添加磁流体可提升磁力矩，但同时影响阶跃响应。研究工作对于后续磁悬浮斜翼节在二维伺服比例阀上的实际应用具有重要的参考价值。     

浸入边界法模拟小圆柱对主圆柱涡脱落的抑制

使用浸入边界法研究了小圆柱对主圆柱涡脱落的抑制.方法中使用非贴体笛儿尔网格,易于处理包含复杂边界的流动问题.采用离散附加力直接加入边界条件方法对虚单元进行重构,使边界条件在浸入边界上精确满足.使用隐式分步法解二维非定常不可压Navier-Stokes方程,通过速度和压力解耦提高计算效率.数值模拟单圆柱绕流及不同位置小圆柱和主圆柱的流动干扰,通过分析流场涡结构和升、阻力系数,得到小圆柱对主圆柱涡脱落的延迟和抑制作用.计算结果与已有实验结论和数值结果对比,计算误差不超过5%,说明浸入边界法可以简单有效地处理圆柱涡脱落抑制这类流动干扰问题.     

直升机的非线性模型及配平

论述和实现了定常平飞时,单旋翼涵道尾桨直升机的非线性全量飞行动力学数值模型,建立模型时,考虑了直升机各部件特别是旋翼的动力学模型,考虑了诱导速度分布的非均匀性,以一具有水平铰外伸量及约束弹簧的铰接式旋翼等效算例直升机的真实旋翼。研究了直升机在平飞时的配平,提出了一种求解直升机平飞初值问题的数值方法-最速下降法。以某直升机为算例,建模并配平,结果表明,该模型符合直升机运动规律。     

大迎角翼—身组合体涡流特性计算

论述了用非线性离散涡法来模拟大迎角下翼－身组合体涡流绕流的计算模型及计算方法；给出了一个典型翼－身组合体的涡流流态计算结果及非线性气动特性和截面压强分布。