Modelling and design of wing tip devices at various flight conditions using a databased aerodynamic prediction tool

A design study of wing tip devices at high and low speeds is described. The basis of the design study is an equivalent drag approach containing both aerodynamic drag gain and structural weight penalty. A comprehensive parameter study is carried out using a rapid aerodynamic prediction tool named Lift and Drag Component Analysis (LIDCA). Adding to an available lifting-line method a databased module for airfoil data is employed that uses results of two-dimensional flow simulations by multidimensional interpolation. Detailed validation studies of the method at high lift and high speed have demonstrated good accuracy. RANS computations of the selected wing tip designs confirm the predicted benefits at cruise condition. The results of the most effective wing tip designs are analysed at both flight conditions. Finally, options for improving the performance at take-off are suggested.     

基于径向基函数与混合背景网格的动态网格变形方法

发展了一种基于径向基函数(RBF)插值和混合背景网格映射的动态网格变形方法。混合背景网格由运动物面附近的各向异性网格和偏离运动物面的各向同性网格组成。物面变形或运动后,包含在各向异性背景网格中的所有物面网格点的位移由指定的物面运动规律得到,而各向同性背景网格点的位移通过径向基函数插值获得。然后变形后的计算网格利用计算网格与混合背景网格之间固定的映射关系代数插值得到。各向同性网格距离物面一定的距离,消除了物面保形对径向基函数插值的要求,从而可以减少径向基函数基点的数量,以提高插值的效率和可靠性。最后,通过二维和三维算例对该方法的变形网格质量、参数影响、网格拓扑依赖性、复杂外形变形能力及网格变形效率进行了测试。测试结果表明,该方法能够高效生成大变形下高质量的变形网格。     

二维流热耦合数值模拟程序及方法研究

本文发展了一套二维流热耦合数值模拟程序,流体域中采用结构化网格,固体域中采用非结构化网格,流体域和固体域耦合面上的网格为非匹配连接。为了使该程序能够同时模拟低速不可压流动,程序中还加入了预处理方法。此外,本文还发展了一种能较好保持通量守恒的界面插值技术,来处理流体域和固体域耦合面上的非匹配网格。数值计算结果与实验值或商业软件结果的对比表明:本程序中的预处理方法有效地改善了低速不可压流动下的计算收敛性和收敛精度,使本程序能够同时模拟高速和低速不可压流动;扩压器和叶栅的流热耦合计算表明本文程序在流热耦合计算方面具有较好的模拟能力,同时也表明本文发展的界面插值技术的可靠性。     

二维非定常欧拉方程的有限元解法

本文通过求解时间相关的欧拉方程，对振荡翼型的非定常二维无粘可压缩绕流作了数值模拟。求解方程时采用了有限体积积分的方法，并用ＶａｎＬｅｅｒ逆风通量分裂格式来计算控制体通量。隐式格式的采用，大大提高了计算效率；同时，通过插值使得空间离散达到了二阶精度。另外，本文采用基于非结构三角形网格的动态网格来处理刚性物体的移动，从而顺利地解决了复杂几何外形物体的非定常绕流的问题。文章最后给出了ＮＡＣＡ００１２翼型绕四分之一弦线作简谐俯仰振动的计算结果。一个周期内的瞬时压强分布以及升力和力矩系数都与实验结果吻合良好，     

数据可视化处理的数学模型

就数据的可视化处理问题,根据数据的分布范围,进行数据处理、网格划分、插值计算,建立了二维数据可视化模型;对比现有的三种空间插值理论的计算结果,设计出比较理想的综合插值法,建立了空间数据可视化模型。通过对比分析,该模型具有计算量少,插值效果好,适用广泛等优点。     

基于双曲正弦函数的拓扑优化在热传导结构设计中的应用

电子设备的热传导结构设计是解决飞机电子设备热防护问题的有效方法,基于经验的传统设计方法存在设计周期长,获得的结果性能并不一定优等问题,因此将拓扑优化应用在飞机电子设备散热设计中,能够快速获取较优的结构布局。建立基于双曲正弦函数(sinh函数)插值模型的热传导拓扑优化数学模型,将该模 型的算法应用于二维、三维热传导算例,并通过 MATLAB编程进行算法实现;该模型与SIMP模型和RAMP模型进行对比,并应用于机载 LRM 模块导热拓扑优化设计。结果表明：基于sinh函数的插值模型较 SIMP插值模型精确,较RAMP插值模型的迭代次数少,能更好地解决热传导结构拓扑优化设计问题。     

角速度可控的五坐标NURBS插补方法实现(英文)

在复杂形体零部件的高速高精数控加工方面,五坐标数控机床起着非常重要的作用。为了适应高速加工的需要,提出了一种用于机床五坐标加工的非均匀有理B样条(NURBS)插补指令。采用该插补指令进行数控加工时,为了获取实时加工过程中的刀位矢量,开放式控制器以相邻刀位矢量间的夹角为参数构造方位B样条曲线。构造映射插补过程中位置曲线弧长和方位矢量角度的多项式曲线,利用该曲线实现线性轴速度及旋转轴速度的同步控制。详细介绍了该插补指令的算法流程及其在开放式控制器中的实现过程。最后,通过加工实例验证了该插补指令的性能。     

小波插值在机载超宽带合成孔径雷达成像中的应用

 通过解方程构造一组424 整数小波,利用整数小波的尺度函数作为插值基函数, 进行波数域算法的Stolt 插值,整个插值过程中只有加法和移位运算,便于硬件实现。仿真计算验证了快速算法的有效性。     

二维轴流压气机转—静非定常Euler方程数值计算

将Euler方程写在包括转子和静子的绝对坐标系下,通过在转—静之间设计一个可以随时间扭动的平行四边形网格计算了二维轴流压气机转—静非定常流场。由于整个方程统一写在动网格下,转—静之间没有使用插值。      

多层快速多极子算法的改进措施

 为精确求解散射问题,采用混合场积分方程（CFIE）、多层快速多极子算法（MLFMA）和共轭梯度算法（CG）的收敛技术。基于传统多层快速多级子算法,详细研究了二维拉格朗日插值节点数对计算精度的影响,并改进了插值方法,在不同的层采用不同的插值节点数;提出了在不同的层采用不同的精度控制来计算多级子模式数;分析了稀疏矩阵的对称性对内存使用的影响以及磁场积分方程对迭代初始值的选择。数值计算结果表明以上改进可较大幅度地提高计算精度和计算效率,同时降低内存使用,可满足复杂目标电磁散射计算要求。     

连分式扩充的粒子群神经网络压气机特性重构方法

对航空发动机压气机原始二维等转速线的数据点进行连分式扩充,通过两次网络训练,增加转速特性数据,在三维空间中进行BP(back propagation)网络模型重构.根据压气机特性数据空间分布的特点,引入压力比函数,调整计算区域,定义网络的输入输出数据,利用试探法确定隐含层维数.采用基于趋利避害原则的粒子群算法对网络的初始权值和阈值进行优化,建立了压气机压比和效率特性的整体代理模型.最后以某型发动机的低压压气机为例进行了压气机特性模型的重构.通过模型的校核与验证表明:采用这种方法建立的模型精度较高,优于普遍采用的传统二维插值方法和普通BP神经网络模型.最终建立的重构模型对于采用选配法、坐标法和部件法等以压气机通用特性曲线为基础的发动机模型的求解,可提高计算精度和迭代速度,具有一定的工程应用价值.      

基于Elman的荧光油膜厚度近效测量

为解决荧光油膜灰度与厚度现标定方法（微小高度装置法）存在采集工况复杂、标定周期较长等缺陷，提出了一种新的近效测量方法，该方法只需极少量标定数据，便可达到与现标定方法相近的效果。近效方法利用了Elman动态神经网络对小样本数据进行量程扩展并引入一维插值算法使扩展后的数据项平滑，对解算出的三维荧光油膜厚度数据采用二维插值算法进行二次平滑以求得完整的厚度分布图。经模拟试验结果显示，通过该近效方法进行厚度测量最终可以清晰、准确、定量地显示荧光油膜厚度分布，与目前广泛使用的微小高度装置法相比效果相近，在荧光油膜汇集处（较厚区域）误差最大不超过±2.5μm，在平滑适中及较薄区域误差不超过±2μm，达到荧光油膜厚度工程测量标准，为飞行器全局摩阻测量提供了一种新标定思路，具有一定的实际工程应用意义。     

Finite element analysis of dynamic stability of bearingless rotor

Dynamic stability equations of bearingless rotor blades were investigated using a simplified model.The aerodynamic loads of blades were evaluated using two-dimensional airfoil theory.Perturbation equations were obtained by linearization of the perturbation.A normal-mode approach was used to transform the equations expressed by nodal degrees of freedom into equations expressed by modal degrees of freedom,which can reduce the dimension of the equations.The stability results of rotor blades were presented using eigenvalue analysis.The shape function matrix was obtained using spline interpolation,which simplified the analysis and made assembly of the inertial matrix,damping matrix,and stiffness matrix a simple mathematical summation.The results indicate that the method is efficient and greatly simplifies the analysis.     

基于滤波反投影算法的TDLAS温度测量研究-流场诊断特别专栏

为提高二维重建温度场的抗噪声、抗干扰能力,本文研究了基于滤波反投影算法的可调谐激光吸收光谱二维层析成像系统。首先,通过数值模拟研究了滤波反投影算法的重建质量与影响因素,结果表明,随机噪声达到10%时,重建图像仍能较好地反映原温度场特征;滤波函数能明显改善重建质量;三次样条曲线插值的重建质量最高。然后通过基于滤波反投影算法的二维层析成像系统对弱旋流燃烧室出口温度场进行重建,将重建结果与热电偶测量结果进行对比,相对误差范围为1.4%-8.9%,较好地反映出温度场的特征,验证了滤波反投影算法的层析成像系统的可行性。     

绕翼型跨音速大扰动对称势流的插值混合差分法

 本文将跨音速定常小扰动势流混合差分法推广到跨音速大扰动定常势流,提出了在局部速度坐标系中求解跨音速精确势流方程的插值混合差分法。作为算例,计算了双圆弧翼型和NACA0015翼型对称问题压力分布,并与已知实验值和双圆弧翼型小扰动混合差分法计算值进行比较,结果接近。试算表明,本文提出的插值混合差分格式是稳定和收敛的。本文解决了M_∞趋近于1的计算难点。     

State interpolation for on-board navigation systems

Common concepts for autonomous on-board navigation systems rely on the numerical integration of a spacecraft trajectory between subsequent measurements of a navigation sensor such as GPS. In combination with a Kalman filter, a predicted state vector becomes available at discrete, but not necessarily equidistant time steps. When used for real-time attitude control or geo-coding of image data, the on-board navigation system has to provide continuous dense output at equidistant time steps, which usually conflicts with the natural stepsize of the relevant integration methods and the non-equidistant measurement times. To cope with this problem, the integrator has to be supplemented by an interpolation scheme of compatible order and accuracy.After presenting a representative formulation of an on-board navigation system and deriving related timing and accuracy requirements, suitable Runge–Kutta methods and associated interpolants are selected and evaluated. Promising results are obtained for the classical RK4 method in combination with Richardson extrapolation and 5th-order Hermite interpolation. The 5th-order Fehlberg method with interpolation due to Enright and, for drag-free scenarios, the 5th-order Runge–Kutta–Nystrom method with 5th-order Hermite interpolation provide a good performance in terms of position interpolation. However, as both methods exhibit significant errors for the velocity interpolation, they are not recommended for use with the outlined navigation filter.     

曲面样条插值法在飞机测压试验中的应用

飞机载荷设计中会涉及求解大量不同攻角和襟翼偏度的机翼气动压力分布,而风洞测压试验所测得的压力分布是有限的,如何选择一种合适的插值算法求出载荷设计中所需攻角和偏度的压力分布,是飞机气动载荷设计中的重要环节。采用样条曲面系数插值法,通过构造机翼压力分布的样条曲面函数,得到两种攻角和偏度情况下的压力曲面,进而插值出在两种攻角和偏度之间的任意一种攻角和偏度的压力曲面。将插值得到的压力数值和试验数据进行相似性对比分析发现:插值和试验值的相似系数较高与试验数据吻合度好。     

固体火箭发动机绝热层温度场的有限元计算方法

利用有限元法计算了固体火箭发动机绝热层在移动边界条件下的二维温度场.采用碳化层-热解面-原始材料的二维碳化烧蚀模型;推导了将热解气体对流项作为源项的有限元计算方法;采用当量对流换热系数和当量热流的方法处理复杂边界条件.采用无限插值法获得移动边界条件下的三角形网格,提高了网格生成速度和网格质量.计算结果表明,利用有限元法计算固体火箭发动机绝热层的温度场收敛性和稳定性都较好.      

压力敏感涂料的标定及在气膜冷却效率测量中的应用

采用自行设计的标定装置对压力敏感涂料进行了标定实验,研究了激励光源光强、温度两种参数对标定曲线的影响.结果表明:当激励光源光强能够使得光敏分子发生氧猝熄条件下,激励光源光强对标定曲线的影响可以忽略不计;标定曲线对温度变化较为敏感.据此提出了采用二维插值方法得出发射光强关于压力比、温度比的表达式,从而对测量结果进行温度修正.在此基础上,将压力敏感涂料测试技术分别应用于圆柱孔、扇形孔气膜冷却效率的测量中,所得实验数据与前人结果符合较好.      

不同插值方法对残缺环阵声源定位的影响

研究了利用不完整麦克风阵列进行声源定位的算法，通过修正傅里叶插值、B样条插值和三次样条插值对残缺麦克风阵列的交叉互谱矩阵的插值补偿，获得了声源定位算法。通过数值模拟对声源定位算法进行了验证，发现在波束形成的声源强度上，傅里叶插值模拟结果最大偏差达到5.21 dB,B样条插值为1.17 dB，三次样条插值为0.80 dB；在声源位置偏差上，傅里叶插值为0.04 m,B样条插值为0.01 m，三次样条插值为0.01 m。结果表明傅里叶插值计算得到的声源定位的强度、位置精度和动态性能最差，三次样条插值最优，B样条插值表现一般。实验验证也得出类似的结论，因此利用三次样条插值计算不完整麦克风阵列的交叉互谱矩阵最佳。