考虑燃料雾化的气液两相连续旋转爆轰数值模拟

为研究气液两相连续旋转爆轰发动机燃烧室内爆轰波的传播特性,以汽油为燃料,富氧空气为氧化剂,建立了欧拉-拉格朗日模型进行二维数值仿真,其中气相方程采用时空守恒元与求解元方法求解,液相方程采用标准四阶龙格库塔法求解。在两相旋转爆轰模型还考虑了液滴雾化破碎过程。计算结果表明：起爆后形成的初始爆轰波经过初始燃料填充区域后逐渐衰减,随后入口附近新生成的压力波经过一系列发展形成了自持稳定传播的旋转爆轰波;旋转爆轰波的传播模态受燃料与氧化剂的喷注压力和氧化剂填充比影响,在不同工况下旋转爆轰波呈现出4种传播模态,即稳定单波模态、稳定双波模态、不稳定双波模态和不稳定单波模态;在双波模态工况下,燃烧室内初始只形成1个爆轰波,后由入口附近局部爆炸产生的压力波发展为新的爆轰波,转化为双波模态后爆轰波的强度略有下降,但燃烧室整体推力更加平稳。     

用于省内公路干线网规划研究的启发式模型

本文阐述了一种将运量在运输线路上的配置和投资在项目上的分配相结合的启发式方法的模型。文章对这种方法的原理与应用进行了研究探讨。     

层流平板摩擦阻力的数值计算

为提高层流摩擦阻力数值计算的精度,采用数值求解N-S方程的方法,对层流平板摩擦阻力进行了计算.计算中采用4种先进上风格式,同时配合以14种不同的壁面法向网格分布和5种常用限制器,并将计算结果与布拉修斯解进行对比.结果表明,van Leer的FVS(Flux Vector Splitting)格式数值耗散大,对网格分布和限制器敏感,不适合用于摩擦阻力计算.限制器和计算格式的粘性分辨率越高,摩擦阻力的计算精度越高.在计算量允许的情况下,可以减小壁面法向第1层网格高度,以提高计算精度,但应同时兼顾网格长宽比.对于摩擦阻力计算,边界层内应至少布置4个网格点.      

低速流预处理Roe格式中的数值粘性

使用Roe格式计算低马赫数流动,需采用预处理技术.在预处理技术中,需要选择适合的预处理参数,但有时很难找到,因而得不到合理的计算结果.通过数值实验,发现使用预处理方法时,格式的数值粘性对低速流动的模拟有显著的影响.提出了在预处理Roe格式中加入可调参数,用以控制数值粘性对计算结果的影响.带有可调节数值粘性参数的预处理Roe格式的数值粘性更接近低速流动的真实耗散,反映了低速流动的真实物理现象.对方腔流动和圆柱绕流的低雷诺数定常流动所做的数值实验,得到了比较满意的结果,表明这一改进是有效的.对同一低雷诺数不同低马赫数的定常圆柱绕流所做的数值实验,均得到了满意的结果,表明了方法的适用性.     

镁合金新材料及制备加工新技术发展与应用

镁合金是国际公认的最有潜力的轻量化材料之一,被誉为"21世纪绿色工程金属",其突出的结构减重效果有望解决我国关键装备和重大工程中迫切的轻量化问题.本文综述了近年来镁合金新材料的发展、镁合金制备加工新技术进展以及镁合金的应用发展现状,对当前镁合金领域存在的缺点进行了总结,并对镁合金未来发展重点及主要任务进行了展望.     

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算法的并行计算格式。最后通过几种数值方法的比较,表明该方法是可行的,也是有效的。     

倾转旋翼飞行器的操纵策略和配平方法

根据倾转旋翼飞行器的构型特点,建立了倾转旋翼飞行器旋翼、机翼、短舱、机身、平尾(含升降舵)和垂尾(含方向舵)的气动力模型,研究了倾转旋翼飞行器的操纵策略以满足直升机模式的悬停/小速度飞行、直升机模式向固定翼飞机模式转换的过渡飞行和固定翼飞机模式的高速飞行,并运用最优方法研究倾转旋翼飞行器在不同飞行速度下作稳定对称飞行时的配平方法.最后以XV-15倾转旋翼飞行器为例,进行各种飞行模式的配平.结果表明:本文所述方法能合理地给出倾转旋翼飞行器在整个稳定飞行速度范围内的操纵量和姿态.     

角接触球轴承分析模型的数值求解

建立了内圈静止、外圈旋转的角接触球轴承分析模型,考虑了滚动体的离心力、陀螺力矩以及弹性流体动 力润滑作用.建立了非线性方程组未知量的约束条件,解决了Newton-Raphson法求解方程易出现不收敛的问 题.以两个角接触球轴承为算例,分析了轴承内、外圈分别旋转对应力的影响,以及转速和径向载荷对轴承工作 性能的影响.结果表明,轴承外圈转时接触应力随转速增大而显著增大,内圈转时则无明显变化;低转速下外圈 转的接触应力明显低于内圈转;较小的径向栽荷下外圈转的接触应力高于内圈转;轴承微区滑动速度受转速和径向栽荷的影响显著.通过Jones球轴承拟静力学分析程序计算,验证了模型和程序的有效性.     

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

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

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

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