柔性杆子杆模型参数的遗传算法优化求解

研究了应用实数型的遗传算法解决柔性机器人子杆法建模中,子杆模型参数的优化求解问题.子杆模型参数形成高度非线性的方程组,它的求解是子杆法能否实现的关键.普通的牛顿迭代法及其变形,要求初始向量足够靠近解,一些所谓大范围收敛的算法又有方程组的Jacobian矩阵非奇异等要求,在此都难以应用.本文基于遗传算法的解法可求得模型参数的最优解,实例仿真显示了这一方法的可行性和高效率.      

三维多段机翼复杂流场的计算和分析

采用中心差分的有限体积方法和分区技术求解了N-S(Navier-Stokes)方程,在分区求解技术中采用了满足通量守恒的内边界耦合条件.分析讨论了三维襟翼的绕流流场,对"剪刀差"处的三维分离进行了探讨,初步揭示了绕流的主要特征及对附近物面压强的影响.同时,分析了三维襟翼前缘的分离流动.      

桥梁断面的气动导数和颤振临界风速的数值计算

选取有代表性的大海带东桥、南京二桥流线型闭口箱梁断面和荆沙桥钝头开口板梁断面,用动网格法考虑流体和结构的耦合作用,用有限单元法求解原始变量二维不可压粘性流体的Ｎ－Ｓ方程,计算其气动导数,并用半逆解法确定颤振临界风速。计算结果和实验值相当吻合。     

有二次涡的激波边界层干扰流动的DPNS和NS方程计算

本文针对具有二次涡复杂分离再附现象的激波边界层干扰流动,数值地考察了扩散抛物化Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ（ＤＰＮＳ）方程组的适用情况。壁面摩阻和压力、主涡和二次涡的涡高和涡长、分离再附位置以及流线图等特性的计算表明：ＤＰＮＳ方程组的数值结果均与ＮＳ方程组的数值结果很好相符。     

Collar网格的生成及其在嵌合体技术中的应用

本文应用双曲型方程生成表面网格和空间网格的方法成生Ｃｏｌｌａｒ网格,并将Ｃｏｌｌａｒ网格应用到嵌合体技术中,在保证计算网格的生成方便快捷而且网格质量高的前提下,有效地解决了为物体结合部的内外边界点寻找插值单元的问题。经数值计算验证,用本文方法求解有结合部的复杂组合体绕流可以得到精确可靠的数值结果。     

扑动翼型的低雷诺数气动特性分析

通过求解引入拟压缩项的不可压Navier-Stokes方程,数值模拟了绕扑动翼型的低雷诺数非定常流动。针对厚度在4%-12%之间的NACA对称翼型,分析了翼型厚度等参数对扑动翼型气动特性的影响。在低雷诺数条件下,对于纯俯仰运动,随着翼型厚度的减小,平均阻力系数也变小。而对于纯沉浮运动,发现翼型厚度对气动特性的影响和俯仰运动有很大的差别,平均阻力系数随着翼型厚度的减小而变大。通过对沉浮运动一个周期流线图的分析,认为这是翼型前缘涡的影响造成的。由于前缘涡的影响,翼型厚度增加,平均压差阻力系数变小,甚至会出现负值。雷诺数的影响研究表明,随着雷诺数的增加,扑动翼型的阻力系数减小的趋势越缓慢。     

动静圆柱坐标系下N-S方程的分区算法及子域干涉处理研究

耦合了相对和绝对圆柱坐标系下三维Navier-Stokes方程,统一了动子域和静子域的控制方程;在子域计算中,采用高分辨率的三阶ENN格式及LU-SGS隐式解法;在子域干涉处理中,研究并对比了动静子域干涉面的重合处理与重叠处理;最终发展了一种新的既高效又准确的分区算法。通过对某带进口导叶的三级跨声速轴流压气机叶栅内三维紊流流场的数值模拟,验证了该分区算法的有效性,计算结果与实验数据吻合较好。     

混合网格方法在栅格翼数值模拟中的应用研究

栅格翼是采用一种蜂房结构设计的升力面和控制面,和传统舵面相比具有铰链力矩小、能在较大攻角下保持升力等特点,并得到了越来越广泛的应用。由于栅格尾翼的特殊结构,网格生成具有很大的难度。本文采用了结构、非结构混合网格的办法,在弹体和栅格翼的物面附近区域,生成适合粘性计算的大长宽比的结构网格,弹体网格和栅格翼网格之间采用非结构网格进行填充,满足非结构网格和结构网格交接面的完全对接。本文基于结构/非结构网格体系采用有限体积方法求解NS方程,对不同舵偏角下的栅格翼构形进行了数值模拟,并通过实验结果对数值方法进行了验证。     

三维Euler方程组隐式LU分解的高性能并行计算

许多非定常无粘流体力学问题的数值模拟都需要利用Euler方程组来进行计算,而由于在隐格式下,所选取的时间步长可以比在显格式下时大得多,所以隐格式越来越受到重视,其中隐式LU分解是最常用的方法之一。对三维Euler方程组,采用隐式LU分解进行计算时,网格点所在的各个对角阵面之间存在数据依赖关系,本文分析了采用区域分解且边界上用显格式代替隐格式进行计算的高效性,在长方体建筑物内的爆炸模拟表明,在有112个CPU的某MPP巨型机上,并行计算效率超过60%。本文还分析了计算结果与串行计算时的差异,以及利用区域重叠减小这种差异的方法,同时考虑了对处理器进行合理的逻辑组织,将计算网格映射到处理器网格,以最大限度减少通信开销的方法。文中最后以一个爆炸毁伤的例子实际说明了所述方法的可行性与高效性。     

An adaptive integration surface for predicting transonic rotor noise in hovering and forward flights

In this paper, a new permeable adaptive integration surface is developed in order to evaluate transonic rotor noise in accordance with FW-H_pds equations(Ffowcs Williams-Hawkings equations with penetrable data surface). Firstly, a nonlinear near-field solution is computed on the basis of Navier-Stokes equations, which is developed on moving-embedded grid methodology.The solution calculated through the present CFD method is used as the input for acoustic calculations by FW-H_pds equations. Then, two criteria for constructing integration surfaces are established based on the analysis of the quadrupole source strength and the nonlinear characteristic.A new surface is determined adaptively by the pressure gradient or density in a given flowfield,eschewing the uncertainties associated with determining cylinder-shaped integration surfaces. For varying hover cases, transonic noises are simulated with new integration surfaces for a UH-1 model rotor. Furthermore, numerical results of the new integration surface derived from the density perturbation value conform better to experimental data than results derived from the pressure gradient.Finally, the integration surface given by jrqj being 0.1, which is an applicable criterion obtained from hover cases, is used to predict transonic rotor noise in forward flight. The computational accuracy of the new integration surface method has been validated in predicting transonic rotor noise of an AH-1 model rotor at different advance ratios.