助推-滑翔式弹道中段的近似解

研究了助推-滑翔式弹道中段的近似解。首先根据不同假设条件,给出以下四种情况的近似解:(1)自由飞行段的解;(2)只考虑气动力常迎角飞行段的近似解;(3)考虑引力和气动力的合成解;(4)平衡滑翔解。再由上述近似解合成助推-滑翔弹道中段的分段近似解。通过与数值解的对比,验证了近似解的合理性。研究结果可为新型助推-滑翔式飞行器在方案论证阶段时的弹道特性分析以及射程、峰值热流和过载等特征参数的估算提供依据。     

关于制导问题的闭合形式解

研究了系统惯性对制导问题闭合形式解的影响,并获得了解析形式的解。此解表明不考虑惯性问题的解只是本闭合形式解在大视线距条件下的一种特殊形式。本解也适合小视线距,视线角速度不稳以及造成脱靶的原因等,都得到了反映。分析结果可用来研究比例导引系数、系统惯性对脱靶量的影响。     

Chapman方程同伦分析方法求解研究

Chapman方程是再入领域的经典方程,适用于对航天器再入轨迹规律的研究,但难以得到其解析解,故提出利用同伦分析方法(HAM)求解Chapman方程解析解。将平衡滑翔再入近似解析解作为HAM的初始猜测解形式,并通过选取合适的辅助线性算子和基函数,得到Chapman方程在小升阻比下的同伦解。仿真结果表明,基于HAM得到的解析解与数值解能够很好地吻合,验证了该方法的有效性和精确性。     

两类n阶非线性方程的两点边值问题

已有人利用上下解方法讨论过非线性n阶常微分方程两点边值问题解的存在性，而边值问题是否存在上下解是很难判别的，在文中，用初等方法讨论了两类n阶非线性方程满足两点边界条件的边值问题解的存在性、唯一性，并对解的存在性唯一性区间进行了估计。     

非线性不平衡转子轴承系统周期解的预测

本文提出一种对非线性不平衡转子轴承系统周期解进行预测的新型算法,它利用系统周期解的稳态及瞬态信息,反解雅可比矩阵,实现对系统周期解的预测追踪,并利用反解得出的雅可比矩阵,求得系统周期解的Floquet乘子以判别其非线性稳定性。文中以刚性不平衡转子轴承系统为例,实现了周期解的预测追踪及非线性稳定性判别,说明了新算法的有效性。      

变比热平行混合近似解

平行混合变比热改进解析解，合理、简明、快捷，定比热解误差大。在准确的改进解与定比热解之间，探索了几种简化的近似解，从中看出，冷、热、混三气各简化按其自身静温比热比ｋ＿Ｔ作计算的近似解，结果总的误差较小，可供参照取用。算例包括不等总压及等总压混合，并对等总压混合后总压降低作出初步分析。     

基于整周模糊度精确解算的空间站伴随飞行器DGPS相对导航技术研究

研究了基于载波相位差分GPS的空间站伴随飞行器相对导航技术,建立了相对导航滤波的系统状态空间模型,并且系统状态向量中包含了GPS整周模糊度;采用扩展卡尔曼滤波EKF实现了相对导航滤波的最优状态估计,并提出了浮点解和定点解相结合的GPS整周模糊度解算方法,以提高相对导航精度。仿真结果表明,与仅采用EKF对GPS整周模糊度进行最优估计获得浮点解的方法相比,采用浮点解和定点解相结合的GPS整周模糊度解算方法,可有效提高相对导航滤波精度。     

具有无穷时滞的非线性Volterra反应扩散系统解的先验界

详细讨论了具有无穷时滞的非线性Volterra反应扩散系统解的先验界与吸引性,运用上下解方法和单调迭代方法,得到了在某些假定下当t→∞时解趋于先验界及在某些特殊情况下t→∞时解趋于正常数的充分条件。     

指标聚类的全部解

指标(样本)的聚类往往具有多重解。本文在分离性聚类准则下,给出了聚类唯一解的充要条件。在重解的情况下,定义了解之间的变换,进而构造解的变换图并证明了该变换图是边-Hamilton图。根据这种边-H性,建立寻求聚类全部解的算法,最后提出了进一步研究的问题。     

无界域上二次规划存在最优解的条件及判别方法

本文给出了无界域上二次规划存在最优解的充分必要条件及判别方法，并通过解一些线性规划来判别是否存在最优解。     

一类二阶三点边值问题解的存在性

运用上下解方法,研究了一类二阶混合非线性边值条件的三点边值问题,通过构造了适当的非线性辅助函数工得到了解存在的充分条件,推广和改进了某些已知的结果。     

Numerical Simulations of Stellar Winds

We discuss steady-state transonic outflows obtained by direct numerical solution of the hydrodynamic and magnetohydrodynamic equations. We make use of the Versatile Advection Code, a software package for solving systems of (hyperbolic) partial differential equations. We model thermally and magneto-centrifugally driven stellar outflows as generalizations of the well-known Parker and Weber-Davis wind solutions. To obtain steady-state solutions efficiently, we exploit fully implicit time stepping. Wind solutions containing both a 'wind' and a 'dead' zone are presented. We emphasize the boundary conditions imposed at the stellar surface. For axisymmetric wind solutions, we use the knowledge of the flux functions to verify the numerical solutions. This revised version was published online in June 2006 with corrections to the Cover Date.     

自相关回归法在地震地面运动功率谱分析中的应用

本文通过数值分析，对几个著名的地震记录的功率谱ARM解的特性作了讨论，并与高速富里衰变换（FFT）法的解进行了比较，讨论了ARM应用于各类地震记录的谱分析的适用性。比较结果表明：对于不同的地震记录，其功率谱的ARM解等价于按一定谱窗宽度平滑化的FFT法解。本研究还对赤池氏推荐的自相关回归次数的经验公式提出了修正意见。本研究的成果对于高层建筑抗震设计有一定的参考价值。     

Efficient computation of energy release rates for delaminated beam structures

The often encountered problem of mode partitioning for delaminated, isotropic beam structures is addressed. For the sake of efficiency, solutions that do not require complex continuum solutions – in practise mostly obtained by the finite element method – are examined. Three approaches that meet the postulated efficiency are reviewed and compared. While one method fits consistently into the concept of linear-elastic fracture mechanics, the remaining two deduce mode partitioning hypotheses solely from kinematic quantities of beam theory. The results obtained through these hypotheses are then compared to continuum-based solutions. The latter are generated by the finite element method and are validated by analytical solutions as applicable. Furthermore, mode partitioning solutions for beams, loaded in transverse shear, are given with excellent accuracy. Finally, the conditions for which the different methods comply with or substantially deviate from each other are discussed.     

高阶线性比例制导系统脱靶量幂级数解

脱靶量是导弹制导系统设计和评估的重要指标,对于一阶环节线性比例制导系统,可以得到脱靶量的解析解,而对于更接近实际的高阶制导系统一般得不到解析解,通常由直接仿真或伴随仿真获得;研究了高阶线性比例制导系统脱靶量的幂级数解,为脱靶量的解算提供一种新的手段。首先,构造伴随系统,假设伴随系统的解为幂级数与指数函数乘积的形式;然后,利用幂级数法给出了脱靶量的幂级数解的系数递推关系;进一步严格证明了脱靶量幂级数解的收敛性;最后针对一阶环节和高阶二项式环节等特殊制导系统,通过选取适当的指数衰减参数,得到了幂级数解系数简化的递推关系,并且一阶环节制导系统的幂级数解和解析解是一致的。在计算脱靶量时,实际用到的是脱靶量幂级数部分和,而部分和项数的确定依赖于指数衰减参数。因此,还分析了指数衰减参数对幂级数解部分和的收敛速度的影响,并给出了指数衰减参数与部分和项数的选取方法,为幂级数解的应用奠定了基础。     

电动力绳系周期解求解及稳定性控制

电动力绳系姿态摆动动力学方程具有非线性特点,且姿态摆动状态直接影响着空间任务能否正常进行。首先,基于哑铃模型及偶极子地磁场模型假设,推导了在倾斜圆轨道上运行的电动力绳系姿态动力学方程。其次,给出了利用摄动法和不动点迭代法相互配合的数值方法求解该非线性方程周期解的过程及解的稳定性判断准则,并针对不稳定周期解提出了相应的稳定性控制策略。最后,通过数值仿真分析,得到了姿态摆动运动的周期解图像,验证了周期解求解算法的有效性,并对不稳定周期解实现了稳定性控制。     

等参梯度有限元应用于非线性梯度材料的研究

将等参梯度有限元法应用于非线性功能梯度材料,推导了2种受力情况下指数体积分数变化的功能梯度材料的应力分布解析解。通过不同的单元划分形式,对等参梯度元数值解、均匀单元数值解及解析解三者进行比较分析,得出等参梯度有限元方法对非线性梯度材料同样适用,为等参梯度有限元的推广应用提供了证据。     

基于全欧拉方程的二维平行剪切层声波产生和辐射的数值模拟

本文基于全欧拉方程对二维平行剪切层声波的产生和辐射进行了计算气动声学数值模拟。空间离散采用了频散相关保持有限差分格式,时间积分采用了低频散低耗散的龙格库塔法。数值模拟结果表明,当声源强度较小时,基于全欧拉方程的计算结果与解析解及基于线化欧拉方程的计算结果均符合得很好,当声源强度较大时(约140dB),线化欧拉与全欧拉的计算结果有较大的偏差,这说明此时非线性的影响不可忽略,线化方程不再适用。      

虚构解方法程序验证

对求解偏微分方程问题的程序的严格验证,一直以来由于其精确解有限,较难进行,针对这种情况,本文给出了一种新的程序验证方法,虚构解方法,该方法旨在解决大型科学计算程序如何较严格地进行程序验证这一问题;该方法通过构造虚构解,修改原控制方程,然后通过对计算结果的分析处理达到对程序的精确验证;文中给出了一组通用的虚构解的选取办法,解决了虚构解方法验证过程中确定虚构解的问题;同时还给出了运用该方法对非结构二维Euler计算程序进行验证的具体例子,网格收敛分析结果表明该方法是一种有效的程序验证方法,并且具有较好的通用性。