一种修正的Osher通量在高阶WCNS中的特性

将一种基于多维修正的Osher通量运用于高阶加权紧致非线性格式(WCNS)中,该修正通量主要在垂直于激波面的界面上增加耗散,能够改善Osher通量的激波捕捉稳定性,同时对边界层和接触间断的分辨率影响非常小。对修正的Osher通量在高阶WCNS中的特性进行研究,通过数值模拟测试了基于Osher通量的WCNS的激波稳定性、热流预测精度、边界层模拟能力、激波边界层干扰模拟能力,并与Steger-Warming通量和Roe通量进行了对比。结果表明修正后的Osher通量比Harten修正的Roe通量具有更好的激波捕捉鲁棒性,而边界层、驻点热流值和激波边界层干扰的模拟则明显优于Steger-Warming通量。上述结果说明了基于修正的Osher通量的高阶WCNS具有较好的激波捕捉特性、热流预测精度和边界层计算能力。     

金属材料微裂纹取向与超声波和频非线性效应

针对非线性超声无损检测金属材料微裂纹取向角度的问题,开展了微裂纹取向与超声波的和频非线性效应研究,建立了超声和频非线性特征系数与微裂纹取向角度的关系模型。理论和有限元仿真实验结果表明,随着微裂纹取向角度的逐渐增大,超声和频非线性特征系数与微裂纹取向角度之间呈现明显的正相关趋势,而且相比二次非线性特征系数,和频非线性特征系数对微裂纹取向检测更为敏感。同时,从超声波平均能流密度（即声强）的角度出发,计算可知和频分量声强会随着微裂纹取向角度的增大而增大,而二次谐波声强基本不会发生变化,同时和频分量声强占比相比于二次谐波声强占比也得到了明显提高。超声波声强计算结果与仿真计算结果趋势基本一致,证明了理论模型的正确性。通过实验验证了模型的有效性,为金属材料微裂纹取向的检测提供了一种有效的手段。      

结冰飞机非线性稳定域确定及安全操纵方法

结冰会恶化飞机的动力学特性,造成飞行包线收缩,威胁飞行安全,研究结冰后飞机的非线性稳定域变化对于驾驶员操纵应对策略设计以及飞行安全的提高具有重要意义。以NASA的GTM为案例飞机,首先对飞机气动参数进行多项式拟合,同时结合结冰因子模型,建立了飞机在结冰条件下的纵向通道动力学模型;然后通过分岔分析方法对飞机在不同程度结冰条件和操纵指令下的飞行状态变化进行了研究,并将其用于指导驾驶员操纵,同时考虑到分岔分析方法的局限性,利用微分流形理论确定了飞行系统的非线性稳定域,并将其作为飞行安全边界;最后针对结冰情形,提出将分岔分析方法与微分流形理论相结合共同用于操纵指导,并进行了操纵时域验证。研究结果表明,结冰会使安全边界收缩,在小扰动的作用下都可能使飞行状态超出安全边界。随着结冰程度增加,飞机的稳定性质甚至会发生变化,此时飞行状态将很难维持在原有的安全边界以内,提出了通过指导驾驶员操纵指令变化使飞行状态到达新的安全边界。研究结果对于飞行安全操纵及边界保护都具有一定的指导意义。      

考虑几何非线性的气动弹性模型缩比方法

随着飞机性能和需求的提高,大展弦比高柔性机翼逐渐成为新型飞机的主要结构形式。这类机翼具有高升阻比、大变形和重量轻等特性,几何非线性效应明显。然而机翼的大展弦比高柔性会带来更大的机翼变形,而机翼大变形则会引起相关的非线性气动弹性行为。为了评估这些非线性气动弹性行为并同时降低设计风险和成本,一般要使用缩比模型进行风洞试验以研究和确认真实飞机的气动弹性特性。基于此,首先使用了传统线性缩比方法来进行缩比,通过刚度质量耦合匹配模态响应法与刚度质量解耦匹配模态响应法这2种线性缩比方法,不断优化缩比结构的设计参数来满足目标缩比值。同时,提出一种动力学有限元模型的非线性静响应-模态协同优化方法,该方法是基于等效静态载荷法的几何非线性气动弹性模型缩比方法,通过2个不同的优化子程序分别匹配全尺寸飞机的非线性静响应和模态振型。结果表明,相比于传统线性缩比模型,考虑几何非线性的缩比模型能够更好地再现全尺寸飞机的非线性气动弹性行为。      

直升机阵风响应缓和控制律设计

本文提出直升机抑制阵风扰动的最优控制律设计方法。优化性能指标要求在阵风扰动下驾驶员处法向过载、飞机姿态变化、以及控制信号能量损耗最小。最终所综合的系统工程实现简便,只需在直升机姿态控制系统基础上增加垂直升降速率反馈通道即可。为跟踪90年代采用高级语言开发飞控应用软件的总趋势,本文使用C语言对优化控制律进行编程。数模混合仿真表明,控制律优化及数字化实现的效果是令人满意的。     

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

分析传统恒温热线热膜流速计(简称 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算法的并行计算格式。最后通过几种数值方法的比较,表明该方法是可行的,也是有效的。     

有色噪声补偿下降阶系统状态的最优估计

在大型空间结构的控制系统设计中,降阶系统的状态估计是非常重要的环节,它直接关系到控制器的性能。针对B.Friedland的工作,本文提出了有色噪声补偿模型,在恒等变换下显式地消去了有色噪声项,利用正交投影定理仅对感兴趣的状态分量导出了无偏、最小方差估计。这组递推估计式与采用增广状态方法相比,不仅在处理容量、计算效率上有很大改进,而且是标准Kalman滤波方程的推广。注意到原系统与降阶系统之间的内在联系,本文采用随机系统分析方法,给出了在高斯、马尔可夫假设下,有色噪声模型特征参数的在线识别技术。数字仿真表明,对偏差补偿起主要作用的是有色噪声补偿模型中的驱动白噪声项,这是对B.Friedland工作的发展。     

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

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