飞机登机门缝典型腔体噪声数值模拟

为了研究某飞机飞行中因登机门修型而产生的强烈啸叫问题,结合机身表面流动特征和登机门缝腔体的实际构型,建立符合实际流动特征的二维不规则腔体模型,采用脱体涡模拟结合Ffowcs Williams-Hawkings积分方程的混合方法对两种类型的登机门带肩壁腔体进行流动和噪声特性的数值模拟,分析了登机门缝纯音噪声的产生机理.结...     

马尔可夫更新方程的迭代算法研究

冷储备系统的动态结构变化使得系统可用度分析必须基于瞬态求解,针对频域两部件冷储备系统马尔可夫更新方程解算困难的问题,提出基于时域伏尔特拉方程的迭代算法,证明了该迭代方程的唯一性和收敛性,并给出迭代算法的误差分析.理论和实例分析表明,当冷储备系统部件的工作寿命和维修时间分布是有界函数时,利用迭代算法可以容易实现时域马尔可夫更新方程的求解.      

中心对称平面三连杆的同步运动

研究中心对称平面三连杆的运动.该模型是常见的力学系统,可以近似描述外层空间中弱引力条件下带有太阳能板的卫星的运动状况.然而,由于该系统是不可积系统,一般的运动状况比较复杂.在本文中对该模型建立Lagrange运动方程,主要研究其同步运动,求出同步运动的精确解,并分析其运动规律,证明了在一定条件下存在周期解.本文的结果为进一步研究该系统的在同步运动附近的复杂性打下了基础.      

基于双星系统定轨的数学方法

通过对双星系统定轨原理的分析,结合卫星的动力学模型,采用联合多时刻观测系统的数据融合方法,得到了一种解算卫星轨道状态的解析方法,同时还给出了误差传递方式。该方法不需要标称轨道,计算方便。理论分析和仿真计算表明,方法确实可行,且定轨精度有所提高。     

欧拉—庞卡莱方程在“和平号”太空站姿态恢复分析中的应用

1997年因"进步号"飞船与空间站对接时发生意外碰撞而使得姿态失去控制,从而导致太阳翼无法朝向太阳而获得能量。使用欧拉方法和四元数方法建立运动学方程,通过欧拉—庞卡莱方程建立动力学方程来分析其姿态恢复的运动过程,研究得出太阳翼能够重新得到光照的可行性方案。     

基于Euler方程的三维机翼厚度与扭角优化设计

使用Euler方程作为流场解算器,结合Powell算法,讨论了在确定的机翼平面形状和翼型的条件下,以最大升阻比为目标的三维机翼截面翼型最大厚度与扭角的优化设计.设计中以机翼沿展向0.1,0.2,0.3,0.8,1.0倍半展长的翼剖面的厚度作为优化控制量,再选取翼梢剖面的扭角增量为第6个控制量,对机翼作了数值优化设计计算,得到了在亚音速时具有相对较大升阻比的机翼其厚度及扭角的优化分布.针对Lockheed-AFOSR Wing B的亚临界和超临界算例结果表明,厚度的非线性分布和负的扭角会改善机翼流场的流动状态,使机翼的升阻比得到提高,优化设计方法是可行的.      

自适应翼型的计算和分析

计算并讨论了不同马赫数（Ma）和攻角(α）下简化的自适应翼型的舵面偏转角的规律。计算中采用代数方法生成计算网络，用有限全积法离散二维可压缩流动的Fuler方程，采用了矢通量分裂Van-leer格式离散无粘通量，用隐式的时间积分和多重网格加速收敛。用数值优化方法取得自适应翼型最优的舵面偏角，并计算其气动特性。为验证自适应翼型增益的效果，计算了双目标优化翼型的气动特性，本文计算了对称圆弧翼型和NACA65A006翼型。计算结果表明，自适应翼型比双目标优化的翼型有更好的气动特性。     

欧拉-伯努利梁平面超大挠性变形问题变分法求解

对于欧拉-伯努利悬臂梁平面超大挠性变形问题,由于其复杂的非线性几何方程,以位移为基本变量进行求解时,通常只能采用如多重打靶、微分求积等数值方法求得梁上离散点的位移值。本文研究了欧拉-伯努利悬臂梁平面超大挠性变形问题变分法求解理论。通过假设多项式形式的梁的曲率试函数以及常数中心线应变,基于欧拉-伯努利悬臂梁的基本假设,推导出了相互耦合的位移函数的精确表达式,并基于变分法理论和三角函数级数展开,推导出欧拉-伯努利梁的非线性控制方程组。利用迭代法对非线性控制方程组中的未知参数进行求解,最终得到欧拉-伯努利悬臂梁的位移函数的解析表达式。利用有限元计算结果对提出的变分法求解理论进行验证,并分别计算了欧拉-伯努利悬臂梁在自由端集中力及位移约束情况下的大变形。算例表明,基于本文的变分法求解理论,利用6个未知参数,即能够精确预测欧拉-伯努利悬臂梁在自由端集中力及位移约束下的超大挠性变形,该研究成果为欧拉-伯努利悬臂梁的超大变形问题提供了新的求解方法。     

一种基于轨道根数约束的最优制导方法

针对航天器空间变轨的制导问题,研究了一种基于轨道根数约束的最优制导方法。在地心惯性坐标系下直接建立航天器的最优控制模型,给出了位置速度表达式和协态变量初值之间的关系;进一步,在不约束真近点角的前提下,推导了5个轨道根数的约束方程,并通过对轨道根数和最优控制理论中协态方程的特性分析,获得了另外两个约束方程。协态变量初值可直接通过求解7个完整约束方程组获得,进而得到最优推力方向。仿真验证了所提制导方法的有效性。     

非零压差边界下间隙流动雷诺方程近似解析解

为了获得间隙密封、滑动轴承等的液膜压力分布解析表达式,将间隙流动雷诺方程的解析解表示为特解和通解之和,其中特解表示为周向坐标函数和轴向坐标函数之和,通解表示为周向坐标函数和轴向坐标函数之积,结合分离变量法和非零压差边界条件,推导获得间隙流动雷诺方程的近似解析解。基于近似解析解计算间隙密封的液膜承载力和偏位角等参数,并与采用有限差分法的计算结果对比,其中液膜承载力偏差小于5%,偏位角偏差小于3°,而计算效率提高了数百倍。