有、无障碍物爆震管起爆过程的CE/SE模拟

用时空守恒元和解元(CE/SE)法模拟有无障碍物情况下爆震管二维爆震波的形成和传播过程,分析爆震波的不同机理,研究在较低能量下,点火区温度、爆震管管径以及障碍物对爆燃-爆震转捩过程(DDT)的影响.控制方程采取无量纲形式,用拓展的隐式梯形法实现刚性源项积分.研究表明,点火能量与点火压力相关,在较低的点火能量下,可以通过提高点火温度、减小爆震管径及加入障碍物加强爆震过程,减小爆燃-爆震转捩过程DDT时间和距离.为研究爆震波点火和DDT过程提供了理论依据.      

存在观测站位置误差的转发式时差无源定位

无源定位中,由于观测站安放在运动平台等原因造成的观测站位置误差会影响无源定位精度性能。另外到达时间差(简称时差)(TDOA)的转发式测量需要将不同观测站截获到的辐射源信号都转发到同一位置,如主观测站。针对这两个问题,提出了基于约束总体最小二乘(CTLS)的无源定位算法。首先将转发式时差的非线性定位方程转化为不需要中间变量的直接线性方程,再基于CTLS算法依次转化为约束优化问题和无约束优化问题,最后推导给出定位近似闭式解。仿真实验表明在观测站误差较大时,该算法与其他算法相比定位精度性能较好。     

2198和5A90铝锂合金脉冲阳极氧化膜制备及耐蚀性

在18%(质量分数,下同)H_2SO_4+5%C_2H_2O_4水溶液中,采用脉冲(PC)电流对2198和5A90两种铝锂合金进行阳极氧化处理。用扫描电子显微镜(SEM)观察铝锂合金阳极氧化膜表面和截面形貌;用能谱仪(EDS)对其成分进行面扫描和线扫描;用动电位极曲线检测氧化膜在3.5%NaCl水溶液中的耐蚀性。结果表明:2198和5A90铝锂合金阳极氧化膜主要由Al的氧化物组成;2198合金氧化膜表面存在细小颗粒,厚度约为150μm;5A90合金氧化膜表面存在微孔,为后续封孔处理提供结构条件,厚度约为180μm;用脉冲方法在两种铝锂合金表面生成的较厚阳极氧化膜具有较高的耐蚀性。     

基于边界层相似性解的放大因子输运模型

为了使二维线性稳定性理论能够适应现代CFD求解技术,通过求解Falkner-Skan边界层相似性方程获得各个形状因子下的相似速度型。基于线性稳定性理论,对各个速度型进行稳定性分析从而获得对应不同速度型的扰动放大因子包络线。最后使用标量输运方程的形式实现包络近似方法中放大因子的当地化求解,并结合原始γ-Re_θt转捩模型中的间歇因子输运方程,实现了自然转捩和分离泡转捩的建模。使用该输运模型对S&K平板、S809翼型、NLR7301翼型和DLR-F5机翼进行转捩预测,结果均与试验结果吻合较好,验证了该模型构建的合理性和可行性。     

基于子目标进化的高维多目标优化算法

多目标优化问题是工程应用中的常见问题,已有的方法在解决3个目标以上的高维优化问题时效果欠佳.如何进行有效的个体选择是求解高维多目标优化问题的关键.针对该问题,提出了求解高维多目标优化问题的子目标进化算法.从理论上证明了多目标优化问题Pareto非支配解的求取,可通过子目标函数值排序,先行选择进化种群中部分非支配解;然后,根据排序信息有选择性地比较进化种群中的元素,减少了比较次数,从而快速获得非支配解集.同时,提出归一化函数差值的Minkowski距离"k近邻"距离计算方法,在进化过程中应用到密度函数中,加速了收敛速度.同当前求解高维多目标优化的算法,在对标准测试函数的计算性能上进行比较,统计结果显示了所提算法在性能上的优势.      

对称奇异线性方程组极小范数解的迭代法

给出了求对称奇异线性方程组Ａｘ＝ｂ极小范数解的迭代算法，其迭代公式为此处／为秩是，ｒ（ｒ＜ｎ）的ｎ阶实对称矩阵，Ｅ为ｎ阶单位阵，ｂ为ｎ维列向量，ｍ为正整数，ε为正实数。证明了这类选代算法的收敛性，讨论了它的事先误差估计式和事后误差估计式。作为应用，给出了求超定线性方程组极小最小二乘解的迭代算法、特征向量导数计算的迭代算法和对于病态正定线性方程组。本文的选代算法可改善病态条件，算例表明也是有效的。     

铝锂合金热处理及其组织性能

本文研究国产８０９０铝锂合金热处理工艺及其组织性能。研究表明：８０９０铝锂合金在５７２±３℃固溶加热、保温３０￣６０分钟,水冷,自然时效时间长达３５天以上,机械性能差;在１９０℃,１９小时人工时效,机械性能σｂ＝５１６￣５７６ＭＮｍ＾－２,δ＝６．２￣６．３％,金相组织均匀。适当延长固溶处理加热保温时间,对铝锂合金强度有益。用扫描电镜对固溶短时保温水冷的人工时效与自然时效Ａｌ－Ｌｉ合金试样进行了继     

基于FMECA信息的测试性验证试验样本分配方法

针对目前国内外测试性验证试验中普遍采用的基于故障率的分层抽样方法考虑因素单一,可能导致样本分配不尽合理,而现有考虑多因素的样本分配方法代表性不足,难以实现工程上的应用的问题,对影响样本分配的因素与故障模式影响及危害性分析(FMECA)信息的关系进行了详细分析,提出了基于FMECA信息的样本分配方法。首先定义了单元影响系数和影响因子向量,提出了向量各元素赋值规则,解决了影响因素选取的问题;其次通过基于逼近理想解排序法(TOPSIS)优化的层次分析法确定了权值矩阵,实现了样本分配。最后通过实例验证表明该方法考虑因素全面,运用灵活,更具代表性和工程应用价值。     

Lyapunov矩阵方程的一种数值解法及部分并行处理

本文设计了求解Lyapunov矩阵方程的一种新方法。所考虑的矩阵方程是 AX—XB=C(1)其中A,B,C分别是m×m,n×n和m×n的已知矩阵。 该方法首先是将系数矩阵A,B初等相似约化为三对角矩阵,即存在可逆矩阵U,V,使U~(-1)AU=A,V~(-1)BV=B,其中A,B为三对角矩阵。然后设计了矩阵方程AY—YB=C的公式解法,分三步: 1)求f(λ)=det(λI—A)的λ各次幂的系数a_0,…,a_m; 2)计算sum from i=1 to m (A_(m-i)-CB~(m-i)),f(B); 3)求解Y。解方程AY—YB=C的方法称为THR算法。 最后经逆变换获得原矩阵方程(1)的解X。 求解矩阵方程(1)的方法称为R—THR算法。该方法的计算量约为m~3+4/3n~3+7m~2n+5nm~2+m~2。 本文给出了R—THR的串行计算的数值例子,并给出了THR算法的并行计算格式。最后通过几种数值方法的比较,表明该方法是可行的,也是有效的。     

舱门机构运动的计算机模拟

1.构造及其原始数据 滑轨舱门收放机构由中、侧滑轨、三角支架、外伸臂、链条条、链轮、液压马达及减速装置等组成(图1)。一根中滑轨布置在43框前飞机对称平面内与飞机外形的下零纵线相贴合,两根侧滑轨分别布置在43框后机身的两侧,舱门的前端由一个接头铰支在三角支架上,三角支架的前支点与链条相连,舱门的后部和两个外伸臂固接,