近地球赤道面轨道上磁性刚体卫星的混沌运动

讨论地球磁场中非自旋磁性刚体卫星在近赤道平面圆轨道上平面天平动的混沌问题.利用Melnikov方法预测出现混沌的可能性并数值验证混沌的存在性.研究表明地球磁场的周期扰动能引起卫星的混沌运动,扰动增大则混沌加强.      

在CAD系统上探索空间飞行器三维动画的实现途径

说明了在CAD系统上实现三维模型动画的可行性,接着论述了三维物体的运动控制方法及在CAD 系统上实现的技术途径。给出了在微机Auto CAD上利用关键帧技术和运动算法控制模拟空间飞行器运动的一些实例,同时也介绍了在Euclid-IS CAD 系统上二次开发运动模拟功能的探索。     

概周期扰动下的混沌振动和维数

本文根据MeInikev的思想,用几何方法清晰地给出并证明了概周期扰动下非线性自治系统出现混沌振动的条件,并分析了二个非线性振动的例子,同时构造数值算法,计算了混沌振动的分维,得到一些有意义的结果。     

物—伞系统运动方程与稳定性判据

物—伞系统运动轨迹计算与运动稳定性分析是飞行器回收系统设计的重要课题。文章采用刚体─刚体模型代替物─个系统，给出该系统的运动方程与稳定性判据。     

慢变参数系统的渐近解、分叉与混沌

研究了一类慢变参数振子系统,通过摄动方法得到其对称周期解的渐近展开式,并与数值解进行了比较.此外,通过系统解的相图、功率谱、倍周期分叉图和Lyapunov指数的计算,分析了系统的倍周期分叉至混沌的过程.结果表明,随着系统的小参数的变化,此系统的运动将经历与Lorenz模型极为类似的分叉而进入混沌状态.此外还可明显看出,此系统比起Lorenz模型相对说来容易处理一些,因为可得出系统的对称周期解的解析表达式.      

飞轮控制的欠驱动刚体航天器姿态控制器设计

针对"三正交加斜装"反作用轮系统中某两个本体轴上的飞轮失效的欠驱动情况,研究了航天器的姿态控制问题.在系统初始角动量为零的条件下,设计分段解耦控制律,实现了姿态稳定.采用欧拉角描述法建立了欠驱动航天器的姿态动力学方程和运动学方程.在系统初始角动量为零的条件下,通过分析方程的解耦特性,设计了分段解耦控制律.该方法经过6次机动控制,可实现姿态稳定.数值仿真验证了方法的有效性.     

基于偏置角动量的欠驱动航天器姿态保持控制

对有扰情况下欠驱动航天器三轴姿态保持控制问题进行了研究,提出一种基于俯仰偏置动量轮和滚动轴推力器的姿态保持控制方法。该方法基于偏置动量航天器滚动-俯仰轴耦合的原理实现,避免了欠驱动零动量航天器平衡点附近欠驱动轴耦合弱的问题;将航天器的姿态运动分为长周期运动和短周期运动,用极点配置方法进行控制律设计,给出保证系统稳定的参数选取范围,求出了系统稳态误差。最后,通过数值仿真验证了所设计的控制器不但能快速消除初始姿态偏差,而且能抵抗外干扰将航天器姿态保持在平衡点附近。     

混沌疲劳初探

引入非线性动力系统分叉与混沌研究方法,探讨疲劳损伤系统的性质突变.基于物理力学观点,提出了疲劳损伤微观系统模型,定性地分析、讨论了疲劳损伤系统的局部分叉(Hopf分叉)与全局分叉(混沌)机制,给出了系统产生局部分叉时的参数取值范围和产生混沌运动的门槛值.     

高速柔性转子支承松动力学特征及动力特性

针对高速柔性转子支承松动的结构特征、力学特征以及多支点转子系统动力学设计的需要,研究了转子支承结构松动引起支承刚度非连续变化的产生机理,建立了支承松动转子系统动力学模型,分析了支承松动转子系统存在混沌运动的条件,即当转子动力特性对支承刚度变化敏感时,受支承刚度阶跃影响,支承松动转子系统会产生混沌运动。根据多支点转子系统动力学特性与支承结构位置、刚度的相关性,采用优化支承位置和支承刚度的方法,使转子动力特性对支承刚度非连续变化不敏感,为多支点高速柔性转子系统的动力学优化设计提供了设计途径。     

温度场中柔性臂超谐共振及复杂运动分析

通过Garlerkin方法建立了考虑阻尼、材料非线性、温度变化和轴向激励的柔性臂系统动力学微分方程.分析了系统存在同、异宿轨道及周期轨道的充分必要条件,通过Hamilton函数得到了对应的参数方程表达式.根据非线性振动的多尺度法,得到了系统在3次超谐共振情况下的一次近似解及其定常解,揭示了系统内各参数之间的关系.对得到的微分方程进行数值计算,分析柔性臂系统参数对纵向振动响应曲线的影响.结果表明,材料非线性和温度变化对系统纵向振动的影响不可忽略;在一定参数条件下,系统有发生复杂非线性运动的可能.为了有效的控制柔性臂的振动,应合理选取系统的物理参数,避免其处于混沌运动状态.     

三维球坐标系下重力波波包非线性传播过程听数值研究

采用全隐欧拉（FICE）格式，对三维球坐标系下重力波波包在中层大气的非线性传播过程进行了数值模拟，模拟结果完整地反映了重力波波包在中层大气的非线性传播过程，并与线性重力波理论进行了定量比较。     

卫星相对运动的新解法

卫星相对运动的研究在航天飞行中有着较为广泛的应用价值。Ｃ－ Ｗ的解对于计算两星近距离的相对运动是很有效的, 但是, 当两星距离变大时其误差就不断地放大。文章采用一种新解法： 首先建立相对运动的解析关系; 然后利用椭圆运动的结果导出三维的近似分析解。它可以精确到偏心率和倾角差的一阶项, 而且不再包括长期误差。理论分析和数值检验的结果表明, 整个新的解的精度大大高于Ｃ－ Ｗ 解的精度。     

温度扰动激发的重力波波包非线性传播过程的数值研究

采用谱配置方法建立了一个中层大气扰动传播的全非线性动力学数值模型，利用该模型对具有高期分布的波状温度扰动在可压大气中的演变和传播过程进行数值模拟，模拟结果表明：温度扰动很快就会激发出两支重力波波包，其中一支向上传播，而另一支则向下传播；这两支重力波波包的传播路径与线性理论给出的射线路径有些差异，进一步的数值分析则表明，给定的温度扰动在经过约3h的演变后，完全转化为波动，在这个转化过程中，79％的扰动能量转化为波动能量。     

三维球坐标系下波包非线性传播过程的一种数值求解方法

采用全隐欧拉（FICE）格式，将三维坐标系下的大气基本运动方程作为控制方程，建立起了三维球坐标系下可压大气中的全非线性数值模式；对纬向边界也进行了合理的处理，并对该模式进行了数值检验，计算结果显示出重力波波包在大气中非线性传播的三维图像：波包整体向上传播，波相关扰动速度随高度增加而增长。     

基于边界元法的三维结构体滑跳运动数值仿真

利用滑跳运动实现近水面高速机动飞行的击水式飞行器是近年来的热点问题，三维结构体高速斜入水冲击载荷求解是滑跳运动数值仿真的关键步骤。针对滑跳运动开展研究，建立了一种基于边界元法的可以考虑结构弹性效应的三维结构体近水面滑跳运动时域数值仿真方法。高速斜入水冲击载荷由边界元法求得，利用状态方程直接积分法求解得到结构节点位移动力响应，并通过滑跳运动动力学模型更新结构体的质心位置及运动速度。针对钢球高速斜入水冲击现象进行数值仿真并与试验结果对比，验证了所提方法的准确性。开展三维刚性球冠体的近水面滑跳运动时域数值仿真和变参分析，得到了结构质量、初始高度、水平抛出速度和半径大小4种参数对球冠体滑跳运动的影响规律。考虑结构弹性效应后，高速斜入水冲击载荷减小，对球冠体滑跳运动也有一定影响。      

吸气式高超声速飞行器耦合运动数值模拟

为了研究吸气式高超声速飞行器在俯仰/滚转两自由度耦合运动下的动稳定性问题，基于气动/运动耦合数值模拟方法并结合理论分析，针对一种类似SR-72构型的吸气式高超声速飞行器开展了进气道通流状态下强迫俯仰/自由滚转耦合运动数值模拟研究。结果表明:强迫俯仰/自由滚转耦合运动下，滚转通道的动力学方程可以简化为有阻尼的Mathieu方程形式，并且可以求得稳定性边界。理论分析表明:滚转通道的动稳定性与俯仰振荡角频率相关，在飞行器滚转振荡固有角频率附近存在2个临界角频率，当俯仰振荡角频率位于2个临界角频率之间时，滚转通道是动不稳定的。在俯仰振荡振幅较小时，数值模拟结果与理论符合较好，但实际的临界角频率与理论分析求解的值有一定偏差；数值模拟结果表明随着俯仰振荡振幅增大，导致滚转发散的角频率范围变得更宽，且向更高频率偏移。      

耀斑后环珥系速度场的定量研究

本文研究了在太阳重力、偶极磁压强梯度和气体压力梯度联合作用下环珥系内物质的下落运动,并在等温、定常假设下用数值方法计算了该环系的二维视向速度场。通过将速度场的理论计算跟文[1]中给出的观测速度场拟合,导出了该环系的有关参数及分布。计算表明环系内物质密度和磁场强度对物质下落运动影响较显著。     

连续变弯度翼型动态气动特性数值模拟

针对翼型后缘连续变弯度运动中后缘边界精确数值模拟的问题，提出了基于二维多项式的时空曲面拟合方法，实现了对后缘边界时空位置的精确模拟。在此基础上，基于OpenFOAM发展了翼型后缘连续变弯度与大幅度俯仰运动耦合的运动边界数值模拟，并计算了翼型耦合运动的气动力，讨论了后缘线性/非线性变形对翼型大迎角动态气动特性的影响规律。结果表明:后缘运动对翼型俯仰运动的升阻特性影响显著，特别在翼型大幅度俯仰时后缘非线性变形对升阻特性改善效果比线性变形大6%~10%。同时还研究了翼型俯仰与后缘变形运动相位差对气动特性的影响。特别地，当相位差为180°时，后缘运动使动态失速时的最大升力提高50.3%，平均升力提高34.6%；当2种运动相位差为0°时，后缘运动使动态失速时的最大阻力降低39.7%，平均阻力降低30.2%，最大升阻比提高22.3%，平均升阻比提高16.8%；同时，翼型在动态俯仰过程中出现负阻力现象，对产生负阻力的原因进行了分析。这些结果可用于指导连续变弯度后缘控制律的设计。      

水上飞机水面降落全过程力学特性数值研究

数值方法模拟水上飞机的水面降落运动姿态和受力的全过程，是结合了水汽两相流、运动学和动力学的复杂问题。采用基于Fluent商用软件开发的整体运动网格方法结合VOF方法进行自由面捕捉，采用六自由度模型进行运动状态模拟，对某型水上飞机水面降落的全过程进行了数值模拟，得到了较好的模拟结果，验证了整体运动网格方法在处理水上飞机水面降落问题时具备良好的适应性。通过模拟得到过载曲线、运动状态参数和水面状况，将降落过程划分为冲击、滑水、漂浮3个阶段，并通过对各阶段的分析总结出对水面降落过程的一般性认识，以期为水上飞机的设计研发提供方法和参考。      

箔条云团的布朗运动扩散模型

假定箔条所在位置处大气运动速度是一个维纳随机过程,同时在忽略箔条的质量时箔条的运动完全反应当地大气的运动,在此基础上建立箔条云团的扩散模型及数值仿真模型. 对扩散模型进行分析得到了箔条云团扩散过程中的数字特征.对数值仿真模型的分析表明必须对大气的运动方差进行时间离散化补偿,进行离散化补偿后箔条云团的数字特征不受模拟时时间步长的影响.模拟结果证明了分析的正确性.