利用截断奇异值分解法解冲激响应逆卷积

在目标识别领域中,目标的冲激响应占有十分重要的地位,然而冲激响应必须通过解电磁系统逆卷积求得。本文直接从时域的观点出发,采用截断的奇异值分解法解这类病态问题,去掉了引起病态的较小奇异值,使这类逆卷积问题趋于良态,并对解球目标的冲激响应进行了数值仿真。结果表明,使用截断奇异值分解法比传统的共轭梯度法得到的解改善了约10dB。     

求解广义特征值反问题的数值方法

讨论一类广义特征值反问题的数值解法,这类问题包括加法、乘法和经典特征值反问题作为其特殊情况。基于行列式和最小奇异值的计算,文中给出了求解这类问题的两个二次收敛的数值方法,描述了在出现重特征值的情况下如何改进其中的一个方法以保持二次收敛性,并且给出了两个数值例子以解释收敛性结果     

用虚拟区域法结合虚拟边界条件求解对称及管壁流动问题

分析了用基于Lagrange乘子的虚拟区域法数值求解时,翼型后缘存在的奇异问题.为了保持直角坐标网格结构简单及可利用快速算法的特点,提出引入虚拟边界条件控制虚拟流动的方法处理奇异问题.用虚拟区域法数值模拟了翼型不可压对称绕流;并把该方法推广应用到管壁流动的数值模拟,给出了喷管内流场的算例.数值结果表明,虚拟边界条件的引入有助于控制Lagrange乘子的变化,从而改善数值模拟的精度.     

水滴变形及其对阻力特性影响的实验研究

为了探索飞机结冰研究中水滴撞击过程的变形行为发生机理和本质特征，提高结冰数值模拟的精准度，利用高速相机对水滴加速运动过程的变形现象进行了实验研究。通过研究，定义了水滴变形过程的4个典型形态，分析了水滴变形与各无量纲量的关系，拟合了水滴变形后的阻力系数曲线和计算式。研究表明:水滴在气动力作用下会发生显著的变形现象，依次历经圆球形、椭球形、半球形、圆盘形4种典型形状；水滴的纵横比随韦伯数(We)的增大呈线性减小，随邦德数(Bo)的增大呈双曲线减小；水滴变形会导致其阻力系数增大，在Re=500左右时，水滴阻力系数曲线开始脱离球形阻力系数曲线，逐渐增大至接近圆盘形阻力系数曲线，这与其变形过程一致。     

最小二乘和奇异值分解法在边界元法中的应用

采用最小二乘与奇异值分解结合的方法，给出求解系数矩阵不满秩的线性代数方程组的数值方法，进而将此方法应用于边界无法中，处理给定外力的第一边值问题，特别地用于处理给定外力的三维裂纹问题。此外，本文还给出求解三维有限体裂纹问题的超奇异积分方程组，并使用有限部积分与边界元法为其建立了数值法。最后计算了若于典型例子的应力强度因子，数值结果与现有文献的相比，符合很好。     

边界层稳定性的非平行性问题研究

采用抛物化稳定性方程的新方法研究稳定性的非平行问题，其控制方程从涡量形式的N－S方程导出，为提高计算精度和收敛速度，在法向采用高精度的高阶谱方法，结合边界层特征安排边界条件的方程，以及对无限区域数值问题作高效的代数变换处理，通过对预估校正迭代过程中推进步长有效的控制，以满足正规化条件，研究了不同类型压力梯度对稳定性的作用，由系列算例得到的增长率变化曲线和扰动中性曲线等，精确地给出了非平行性对流动稳定性的影响，计算的结果与相关数据符合甚好。     

二次约束优化方法在结构动力模型修正中的应用

将二次约束优化方法应用于结构动力模型修正问题中.该方法将模型修正问题转化为一个带二次约束的最优化问题.应用奇异值分解,给出了结构模型修正的数值方法,并进行了数值实验.通过与其他方法比较,本文提出的方法是有效的.     

周期激励Ueda电路中Hopf分岔与混沌的参数区域

周期激励Ｕｅｄａ电路是一种具有奇异混沌行为的强非线性电路。本文基于Ｈｏｐｆ分岔条件给出了一个确定混沌参数区域的方法。首先使用谐波平衡法获得基谐波幅度随时间变化的方程；然后由该方程平衡点的稳定性得出平衡点的Ｈｏｐｆ分岔出现条件；最后通过数值方法，计算出在Ｈｏｐｆ分岔曲线周围系统的Ｌａｙａｐｕｎｏｖ分量，确定Ｕｅｄａ电路出现奇异混沌吸引子的参数区域。结果表明，Ｕｅｄａ电路的Ｈｏｐｆ分岔曲线附近确实存在着奇异混沌吸收子，这为预测Ｕｅｄａ电路奇异混沌吸收子的出现提供了一种有效的方法。文末还讨论了混沌与Ｈｏｐｆ分岔之间的关系。     

雷诺数对粗糙表面翼型气动性能的影响

采用Spalart-Allmaras湍流模型研究风力机专用翼型DU00-W-212的气动性能。在雷诺数小范围变化的情况下,应用数值模拟的方法对粗糙翼型和光滑翼型的气动特性进行对比。探讨雷诺数小范围变化对粗糙翼型气动特性的影响,并分析升阻力的变化机理。比较升阻力曲线,发现粗糙度对升力系数影响最大的区域在迎角10°附近。小迎角的情况下,前缘粗糙度会使得翼型的升力系数下降,阻力系数反而上升。在大迎角失速的情况下,粗糙翼型的阻力系数反而小于光滑翼型的阻力系数。为了增加叶片适应恶劣沙尘环境的能力,在叶片设计和翼型排布的过程中应该尽量避开升力系数最大的迎角处。     

霜状结冰及其对机翼气动特性影响的数值模拟

在霜状冰结冰过程及结冰对机翼气动特性影响的数值模拟中,基于壁面函数法引入了粗糙度影响,结合欧拉坐标系下空气-过冷水滴两相流动控制方程的计算,文中以NACA0012为对象进行了霜状冰结冰过程的数值模拟,并把计算得到的冰形与试验数据及国外的结冰预测软件的结果进行了对比.本文同时考察了结冰对机翼气动特性的影响,结果表明结冰后最大升力系数降低了26%,失速攻角降低了3°,并且阻力系数也有了增加.     

OPTIMAL DESIGN AND AERODYNAMIC CALCULATION OF WING CONFIGURATION OF CIVIL AIRCRAFT

ＯＰＴＩＭＡＬＤＥＳＩＧＮＡＮＤＡＥＲＯＤＹＮＡＭＩＣＣＡＬＣＵＬＡＴＩＯＮＯＦＷＩＮＧＣＯＮＦＩＧＵＲＡＴＩＯＮＯＦＣＩＶＩＬＡＩＲＣＲＡＦＴ￥ＷａｎｇＬｉａｎｇｙｉ（ＤｅｐａｒｔｍｅｎｔｏｆＡｅｒｄｄｙｎａｍｉｃｓ，ＮＵＡＡ２９ＹｕｄａｏＳｔｒｅ...     

基于面积律自动判别的飞机机身外形修形方法

为提高超声速飞机方案修形效率，研究了依据飞机截面积分布曲线分析结果进行飞机截面精细化修形的方法。提出了基于面积律自动判别的截面积曲线分析法，能够依据面积律准则快速对可能增大激波阻力的站位进行自动判别和定位。在此基础上建立了飞机机身几何外形快速辅助修形方法，实现了对超声速飞机方案的待修形位置自动快速分析与定位，并根据截面位置部件分析和截面积曲线变化分析，提供针对截面积的量化调整建议。对一种超声速轰炸机方案进行机身几何外形修形的算例表明，本文提出的方法能够快速识别和定位截面积曲线上拐点、阶跃点、极值点及尖角点等增大激波阻力的位置，根据修形方法调整后的方案阻力系数下降了1.5%，验证了方法的有效性。     

应用后缘装置的跨声速层流翼型多岛优化设计

高雷诺数状态下,自然层流技术(Natural laminar flow,NLF)是减小机翼表面湍流摩擦阻力的有效方法。然而由于层流翼面上大范围顺压梯度的存在使得后缘处的恢复压差更大,产生更强的激波。因此在减小摩擦阻力的同时又增加了激波阻力。本文采用后缘装置(Trailing edge device,TED)来控制翼型后缘处的激波强度,基于线性稳定性理论(Linear stability theory,LST)的eN方法对流动进行转捩判断,进而应用多岛并行多目标进化算法(multi-objective evolutionary algorithm,MOEA)以获得大范围层流区域和弱化激波强度为目标对翼型进行优化设计。优化结果表明合作均衡策略耦合进化算法可以快速地捕捉到该多目标问题的Pareto阵面解,阵面上翼型的波阻力和摩擦阻力性能较初始翼型大大改善。同时,采用后缘装置控制激波强度时,无论在设计点还是偏离设计点时,优化后翼型均具有良好的升阻力特性和鲁棒性。     

底部排气弹DR582雷达阻力系数数据分析

通过底部排气弹ＤＲ５８２雷达阻力系数的数据处理，得到底排减阻的变化规律，可判断底排装置设计的合理性及改进方向，为设计良好的弹道性能的底部排气弹提供了方法，也为修正阻力计算公式提供了依据。     

疲劳裂纹扩展门槛值试验方法的探讨

本文就恒K控制的K梯度法疲劳裂纹扩展门槛值试验的实施方法。从硬件、软件及数值处理等方面进行了探讨,提出了用软、硬件滤波相结合的方法来减小裂纹长度测量值的波动,以提高试验的测量与控制精度。同时讨论了采用通用性较好的能描述全过程的疲劳裂纹扩展速率表达式进行外推估计疲劳门槛值的方法及其实际应用。本文也指出用王永廉提出的表达式将包含有低速扩展段数据点的试验数据进行曲线拟合外推,即能得到实用的门槛值数据,而试验时间比标准试验要少得多。     

动态有限元模型建立的逆方法

本文介绍一种用测量的特征对(归一化振型和固有频率)建立结构系统动态有限元模型的方法。有限元模型用质量和刚度参数表示,这些参数由正交条件和特征方程的加权欧氏范数取极小值来确定。对不完全测量特征对条件下解的唯一性和不完全测点集的处理方法进行了讨论。文中给出的实例表明,这里提出的方法是可行的。     

联接翼布局低速纵向气动特性初探

本文基于低速风洞纵向测力以及涡格法的理论计算结果,初步探索了联接翼布局的低速纵向气动特性,并与相应的机翼、尾翼相分离的正常布局的试验结果作了比较。结果表明,联接翼布局具有许多优点,如较大的升为线斜率C_~α、较大升力系数C_(max)、较大的纵向稳定度、相当小的诱导阻力C_(zi)和较高的巡航升阻比K。还研究了在联接翼前部配置鸭翼对进一步提高和改善其纵向气动性能的可能性。     

暂冲式高速风洞流场控制系统建模与仿真

为满足中国空气动力研究与发展中心的2.4m跨声速风洞流场品质改进的需要,有必要建立一个高效的风洞流场控制模型作为控制器设计的验证平台。由于难以建立精确的空气动力学模型,且2.4m 跨声速风洞长期运行积累了大量的试验运行数据的实际,数据建模成为建模方法的首选。在硬件上,建立了基于反射内存技术的流场控制仿真系统,以获取现场采集的数据。建模方法采用数据建模方式,主要是利用系统辨识理论,将整个系统看成是一个“黑箱”,利用现场采集的数据来确定系统的参数和输入输出间的映射关系。采用以非线性自回归滑动平均模型(Non-linear Auto-Regressive Moving Average Model with Exogenous Inputs,NARMAX)作为风洞系统的数据模型,应用互信息法、曲线拟合法和伪最近邻点法分别确定了模型中采样间隔、时间滞后以及阶次3个参数。对比了最小二乘线性回归、BP 神经网络以及最小二乘支持向量机(LSSVM)3种方法对模型的拟合效果,确立了最小二乘支持向量机作为最终的拟合方法。为了提高仿真的精度,根据风洞运行的特点,将其整个过程划分为冲压、启动和调节3个阶段,分别建立了各个阶段的子模型。由于风洞系统是一个多输入多输出系统,并且延迟和阶次较大,采用了基于信息熵的数据压缩方法,实现了简化子模型规模的目的。最后,采用多模型融合的方法将各个阶段的子模型通过加权的方法来完成融合,从而构建起整个风洞系统的模型。稳定段总压和驻室静压分别通过所建模型得到,最后通过马赫数的计算公式得到试验段马赫数值。仿真结果表明:所建模型在运行包络线范围内的试验工况下,总压预测精度达到0.1%、马赫数预测精度基本达到0.001,达到了研究的目的。该项工作的开展较为系统地建立了暂冲式风洞的流场控制模型,建立的模型将为下一阶段基于现代控制理论的控制器设计奠定基础。     

疲劳S-N曲线的加权最小二乘法拟合

提出了一种拟合材料疲劳S—N曲线的加权最小二乘法，各应力水平下的拟合点的权重反比于均值置信区间的长度。两个算例结果表明：本方法较好地考虑了S—N曲线分散性的物理特性。     

三角网格模型的细分曲面拟合

提出一种用Loop细分曲面拟合三角网格模型的新算法.首先对网格模型进行特征识别,然后把经过简化和形状优化的网格作为拟合细分曲面初始控制网格.通过对控制网格顶点的循环修正、网格形状优化、局部自适应细分来求解拟合细分曲面控制顶点.该算法不仅避免了求解线性方程组,克服了拟合控制网格的扭曲现象,而且达到了用较少的控制网格拟合出反映物体细节特征并满足精度要求的分片光滑(片内除奇异点C1外其余C2连续)的Loop细分曲面的目的.实例表明,该算法用于测量数据的曲面重构是有效可行的.