预移相型恒温热线(膜)流速计的动态响应方程

分析传统恒温热线热膜流速计(简称 HWFA)动态方程后,提出了一种新型HWFA 线路模型——主电桥预移相模型,推导出真实条件下小扰动动态响应方程,并按此原理设计制造了崭新的高性能 IFV—900型 HWFA 样机。该机全部革除了传统HWFA 的三个调节参量,大大地简化了调节步骤;不自激振荡的范围明显增宽,极大地减少了使用上的麻烦;频带宽度比当今世界流行的 HWFA 增宽了2～5倍;加上一系列智能化功能后使自动化程度达到了新的高度。这一新型 HWFA 的问世将可能导至 HWFA技术的重大改进和新的变革。首先阐明主电桥预移相模型的提出和真实条件下小扰动动态响应方程的推导。随后将用若干姐妹篇阐明动态方程频率响应调节元物理值的解析解,频率特性实验分析手段的电扰动试验响应情况的理论分析等内容。     

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

飞行器纵向非线性动力学建模对稳定性的影响

飞行器的稳定性分析与建立的动力学分析模型密切相关.传统的稳定性分析方法是基于小扰动假设建立线性化模型,当飞行状态中非线性特征明显时,其结果误差会带来多大的影响是人们关心的问题.通过建立飞行动力学非线性模型和气动力非线性模型,采用时间推进方法对两类非线性问题进行了比较系统的研究.数值计算结果表明,飞行动力学非线性模型与线性模型相比,周期增大、幅值减小、衰减率增大;而菲线性气动力模型对周期几乎没有影响,但对幅值有一定的影响.     

直升机缆位系统建模与仿真

针对直升机绞车收放时缆位保持问题,建立了负载出水回收进舱阶段缆位的动力学模型。由于浮力与洋流作用突然消失,负载出水后发生摆动,将负载的摆动作为带有理想约束的质点系动力学系统,采用拉格朗日方程建立动力学模型。考虑直升机悬停提升的情况进行仿真,仿真结果表明建立的数学模型能够较好的模拟负载摆动时的缆位变化。     

多芯片组件（MCM）用金锡焊带材料电阻率及热导率温度特性研究

对牌号为Au₈₀Sn₂₀的金锡焊带材料在208~423 K的电阻率及热导率与温度的函数关系进行了研究,并对其在多芯片组件（MCM）中的传热效果进行了评估。分别对材料在208~423 K中5个温度点的电阻率及4个温度点的热导率进行了测试,基于理论模型建立电阻率/热导率随温度变化的函数关系,最终采用模拟热扩散数值方法评估材料在高温下的传热能力。结果表明,采用修正函数模型后,金锡焊带材料在208~423 K下热导率与电阻率的关系符合测试结果,随芯片表面温度的边界条件从208 K升高至423 K,采用变温热导率模型得到的热流密度模拟计算结果相比理想化恒定热导率模型的差异性逐渐升高至5.5%。综上,金锡焊带材料热导率与电阻率的关系符合Wiedemann-Franz法则修正后的Smith-Palmer方程,在该材料传热设计时应考虑其热导率温度效应。     

涡喷发动机加速控制系统数学模型及数字仿真

 本文介绍了航空涡轮喷气发动机加速控制系统数学模型建立、数字仿真方法。由于加速过程是一个大偏离动态过程,因而发动机与加速控制器模型均用非线性方程描述,对非线性方程组求解和系统迭代方法本文进行了讨论。算例中分析了在不同飞行条件下发动机的加速规律。     

叶轮机械三维不可压Euler方程的直接解法

本文用近似因式分解交替方向隐式格式直接求解叶轮机械中原参数三维不可压Euler方程组。在连续方程中引入“拟压缩性”之后,将原来类型不确定的控制方程变成双曲型方程,给定初边值之后用时间推进法求得稳定解,对一个单级压气机转子内部流场进行了计算,结果与实验值比较符合。     

战斗机弹射救生系统数值模拟研究

对弹射救生系统这种几何外形复杂、不规则非流线型钝体,以Jameson有限体积法为基础发展了一套基于高质量混合网格的Euler方程解算器,在Euler方程解算器基础上,采用S-A湍流模型发展了一套同样基于混合网格的N-S方程解算器,初步建立了适合弹射救生系统这类几何外形复杂物体的数值模拟技术平台,并采用多种加速收敛措施和合适的耗散项模型使得所发展的解算器在软件实现上鲁棒性很强。通过选取合适的计算起始迎角,对弹射救生系统成功地进行了大迎角和大侧滑角绕流计算,获得了与风洞试验比较吻合的计算结果,并且可以清楚地观察流场中马赫数分布、弹射救生系统各部位压力分布等流场细节,对弹射救生系统空气动力特性分析具有重要的指导意义。该计算平台可作为风洞试验的合理补充、加强和完善。     

Hamilton正则方程半解析法的收敛和对称分析

近年来,Hamilton正则方程半解析法在工程问题上的应用越来越广泛,但至今未见有关这种方法收敛性和对称性问题研究的文献。基于Hamilton正则方程的半解析法理论,通过变分原理详细推导了Hamiltonian元素的固支和简支边界公式及对称边界公式。多个实例的数值研究表明：随着网格加密,Hamilton正则方程半解析法的收敛速度快于一般传统位移有限元法,对称解法的效率明显优于整体解法。     

边条翼布局双垂尾抖振的数值模拟

对边条翼布局的双垂尾抖振进行了较为深入的数值模拟研究。模拟来流马赫数为0.2,迎角为10°~40°。通过非定常Euler方程计算各迎角下的非定常流场及垂尾根部弯矩系数。并将随时间脉动的根部弯矩系数进行计算得到根部弯矩系数均方根值,从而得出根部弯矩响应大小随迎角的变化曲线。结合流场特性对该布局双垂尾抖振的发生机理及抖振响应随迎角的变化规律作了深入分析。结果表明:该边条翼布局双垂尾抖振主要是由于边条破裂涡作用在垂尾上的脉动载荷引起的。最后,将垂尾根部弯矩响应的计算结果与该模型双垂尾抖振的风洞实验结果作了比较,结果符合得较好。