全机模型静强度分析的多区域并行计算

在Windows XP 64位操作系统环境下,利用若干普通PC机搭建了机群系统,实现了分布内存式并行（DMP）计算。采用3节点和4节点薄壳单元和2节点梁单元建立了约260万自由度（DOF）的全机有限元模型。采用封装的METIS区域分解方法（DDM）,并利用大型通用有限元软件MSC.Marc,实现了全机静强度模型的并行求解,为大型复杂有限元模型的求解提供了一种高性价比的解决方案。     

箱梁颤振过程中旋涡的演化规律及气动力的变化特性

利用粒子图像测速技术（PIV）观测了箱梁颤振过程中模型周围流场的旋涡特征,以模型扭转振动位移作为参考信号,采用相位平均的方法研究了旋涡规律性演化对模型周期性振动的驱动作用。风速较低时,箱梁振幅很小,其尾部风嘴附近上侧的旋涡尺度也很小,不易观测到,而下侧的旋涡尺度较大,其形状接近于圆形。当风速接近颤振临界风速时,箱梁振幅明显增大,并且模型尾部风嘴上侧的旋涡尺度也显著增大,达到与下侧旋涡尺度相匹配的程度,模型尾部风嘴上下侧旋涡的交替作用主导了结构振动直到模型振动发散。基于流固松耦合的计算策略,采用计算流体动力学（CFD）的数值方法模拟了箱梁颤振临界状态下的绕流特性,结合正交特征分解（POD）的方法研究了模型颤振时刻表面压力的空间分布特征,通过分析发现在颤振过程中箱梁表面波动压力的主控成分向迎风侧风嘴漂移。     

展向振荡对激波/湍流边界层干扰的影响

周期振荡作为一种有效的壁面流动控制手段受到广泛关注,而其对激波/湍流边界层干扰的影响目前鲜有研究。本文采用高精度直接数值模拟（DNS）方法对马赫数2.9、12°激波入射角、强振荡下的激波/湍流边界层干扰进行了系统研究。通过与无振荡工况的定量比较,揭示了展向强振荡对干扰区内复杂流动结构的影响规律及作用机制,如分离泡尺度、物面压力脉动非定常特性、物面剪切的非定常特性及统计特征等。研究发现：在展向强振荡作用下,分离点位置提前,间歇区长度增大;同时由于分离泡内强黏性耗散的影响,展向振荡的穿透高度约为分离泡高度的4%,因而对流动结构不会产生实质影响。但展向强振荡会对壁面附近流动造成显著影响,如强振荡诱导的壁面展向速度远大于流向速度,造成流向剪切与展向剪切之间夹角的概率密度函数峰值从0°偏移到80°~90°之间。物面压力及剪切本征正交分解分析表明,展向振荡会导致模态能量从低阶模态向高阶模态转移,降低低频运动的能量占比,增强再附后Görtler涡等壁面附近旋涡结构的强度。     

基于阻力分解策略的翼型阻力计算方法研究

研究并应用了基于翼型阻力分解策略的中场积分计算阻力的方法。中场积分方法是近年来提出的一种新的计算阻力的方法,方法将流场分为三个不同的区域：波阻区、型阻区、数值耗散区,通过对不同区域进行积分将阻力分解为波阻、型阻、数值耗散阻力。计算结果表明中场积分方法与传统的表面积分方法和远场积分方法相比,不仅可以将阻力分解成不同部分以便于进一步分析阻力产生机理,而且在计算中对网格量的要求较低,即使在网格量较小的情况下也能得到较精确的结果。     

AP复合固体推进剂燃烧模型

本文提出了一个复合固体推进剂的综合燃烧模型和相应的数学处理方法。该模型对凝聚相反应作了更细微的分析,特别重视和强调了各种热量传递对燃烧所起的作用,并以推进剂中处于连续相的粘合剂线性分解速度表示推进剂的线性燃速。根据此模型编制程序并进行了大量计算。计算结果与实验符合得很好。对AP/HTPB多级配推进剂的计算结果中,相对误差小于10％的占80％,对Ae/AP/HTPB推进剂则占90％以上。     

非标准H~∞控制系统结构的研究

研究了具有无穷远零点的非标准四块Ｈ∞控制系统的结构，证明了通过弓引入适当的假想被测量信号和控制信号，可将四块Ｈ∞控制的求解问题等价为一块Ｈ∞控制的求解问题。对于增广得到的一块Ｈ∞控制问题，本文证明了其广义对象的散射模型具有Ｊ－无损分解，且其中的单模矩阵因子是上三角的，此结构保证了假想被测量信号和控制信号不被利用。本文还进一步分析和揭示了Ｈ∞控制系统的无损性和对偶性。     

计算大型实对称特征问题的 Lanczos-QR 算法

为了计算大型实对称特征值问题Ｋｘ＝λＭｘ的少数低阶特征值对,本文给出Ｌａｎｃｚｏｓ－ＱＲ迭代方法。首先,给定初始迭代向量ｖ１,作ｍ步Ｌａｎｃｚｏｓ分解：ＫＶｍ＝ＭＶｍＴｍ＋ｈｍｅｍＴ。取Ｔｍ的ｄ个最大特征值为移位量,对Ｔｍ进行ｄ步带原点位移的ＱＲ分解。然后,修改初始迭代向量ｖ１。迭代地重新开始这一过程,迫使初始迭代向量ｖ１进入需求的特征子空间,从而使残量‖Ｋｘ－θＭｘ‖→０。数值例子表明,该方法收敛性强,且稳定、有效。     

基于时滞分解的线性时滞系统主动容错H_∞控制

采用一种具有广义内模控制结构的反馈控制器,研究一类线性时滞系统的主动容错H∞控制问题。针对可由外系统描述其动态行为的故障,设计基于观测器的故障估计器,并将故障估计用于反馈控制的故障补偿。采用一种新的时滞分解方法选取Lyapunov-Krasovskii泛函,推导并证明问题可解的充分条件,最后给出主动容错控制系统的设计算法。最后通过算例验证了提出方法的有效性,并且利用时滞分解方法,保守性减小。     

卡尔曼滤波器UD分解滤波算法的改进

简介Kalman滤波的UD协方差分解算法,针对工程应用中出现的由于UD协方差分解算法本身的局限性 而产生UD因子分解失真问题,提出一种简单有效的解决方法。从而有效地克服UD算法本身的局限性而造成的UD因 子分解失真,改善滤波效果。     

中国航天系统工程探讨

一、系统工程是科学技术 1978年,《文汇报》发表了钱学森《组织管理的技术——系统工程》一文,第一次在媒体上宣传系统工程这门科学技术。事实上,中国航天的系统工程方法实践早在上世纪60年代初就在他的倡导下开始了。那时候,党和国家决定按照自力更生为主,力争外援为辅的方针研制导弹。导弹研制是复杂的高技术系统工程,涉及多种专业技术,需要巨大的资源投入,具有很高的