基于CFD-DEM耦合数值模拟的全尺寸直升机沙盲形成机理

为了研究直升机发生沙盲时沙尘云在悬停流场中的状态和分布规律,采用基于雷诺平均Navier-Stokes方程和Menter 剪切应力输运（SST）k-ω湍流模型的数值模型,通过应用程序编程接口耦合基于Hertz-Mindlin （No Slip）碰撞接触模型的离散元模型,并基于"多球法"构建了更加真实的非球形沙尘颗粒,计算了沙尘颗粒在流场中的运动和分布状态。与可得到的试验结果进行对比,验证了数值方法的有效性,并进行了地效状态旋翼拉力系数、桨尖涡位置以及沙尘云形成的宏观轮廓图的计算。应用所建立的方法,对不同悬停高度的直升机地效流场进行了计算,给出了流场的涡量以及速度云图,着重对比了不同悬停高度下沙尘云中沙尘颗粒速度和分布情况,分析了地效流场对沙尘颗粒状态的影响及沙尘云的形成机理,并计算出了沙尘云宏观分布图。结果表明：流场中大部分沙尘颗粒只能在地表随流场扩散而并不能形成沙尘云;沙尘云中外层空间的沙尘浓度比内层高;位于桨盘平面下层区域的沙尘颗粒以径向运动为主,切向速度较小,而位于桨盘平面上层的沙尘颗粒速度方向各异,速度大小接近。     

树形键合剂包覆RDX及其相互作用研究

采用扫描电镜(SEM)，X射线光电子能谱(XPS)和显微红外光谱(MIR)技术，实验研究了树形键合剂(DBA)对RDX晶体的包覆性能及其界面相互作用机理。结果表明，树形键合剂可以在RDX表面形成一层粘流状的吸附层，具有良好的包覆性能，显微红外光谱的数据也表明了键合剂与RDX的NO2基团形成氢键，从而证实了DBA作为RDX键合剂的有效性。     

低雷诺数下二维翼型绕流的流场特性分析

采用高精度有限差分格式,对低雷诺数下二维翼型绕流进行了直接数值模拟,计算了雷诺数为1.0×104,NACA0012翼型0°,4°以及10°攻角下的流场。计算结果表明:在0°和4°攻角条件下,翼型绕流尾迹区的统计特性相似,0°攻角下的统计量值具有很好的对称性;在距翼型尾缘0.3弦长以后的尾迹区,旋涡排列成类似涡街的结构,涡量的极值、压力的极小值都位于旋涡中心,沿着流向,涡量极值的绝对值逐渐减小,压力的极小值逐渐增大。10°攻角下,翼型上表面从前缘开始分离,尾迹区统计分析结果所得图象与0°和4°完全不同,数值上较后者结果大;在翼型尾缘处,涡量的卷吸,高压力区域的形成,是旋涡脱落的条件,正向和反向旋涡的交替脱落,形成了类似涡街的结构。      

纵列式双旋翼悬停状态气动干扰特性参数影响分析

建立了一个纵列式直升机双旋翼气动干扰特性分析的自由尾迹迭代方法。在该方法中,考虑了双旋翼以及旋翼与尾迹之间的相互干扰影响,将旋翼尾迹模型、桨叶气动力模型以及旋翼配平模型进行耦合求解。同时给出了一个适合于双旋翼干扰计算特点的配平方法。计算了悬停时干扰状态下的双旋翼诱导速度分布以及旋翼性能,并与试验数据进行对比,验证了方法的有效性。应用上述方法,对比分析了纵列式双旋翼与单旋翼的性能,结果表明:悬停时纵列式前、后旋翼的性能都比单旋翼时的要差。文中进一步系统地分析了悬停状态双旋翼纵向间距和轴向间距对气动干扰特性的影响,结果显示:当纵向间距为1.3R时,纵列式双旋翼相比两单独旋翼需要额外的附加功率为8.5%,且随纵向间距的增大,附加功率减小并出现负值,当纵向间距为1.85R时,附加功率最小。     

基于自由尾迹方法的共轴式双旋翼流场分析

建立了一个适合于共轴式直升机双旋翼流场特性分析的自由尾迹计算方法.该方法中,桨叶模型采用了升力面/涡格法,尾迹则使用畸变的自由尾迹模型,以更好地模拟流场特性及计入桨尖三维效应.通过旋翼下洗速度的计算值与可得到的实验值对比,表明了该计算方法的有效性.然后,着重计算分析了共轴式直升机的旋翼尾迹结构和诱导速度分布,并与单旋翼情况进行了对比.结果表明:悬停时,与单旋翼相比,共轴式双旋翼在旋翼下方的诱导速度更大,但并非两幅单旋翼诱导速度的迭加;前飞时,随前飞速度的增加,上下旋翼间的相互干扰逐渐减弱.     

某高层建筑实测风场和风压的相关性研究

通过对台风“天兔”登陆时厦门沿海某高层建筑的风场及建筑物迎风墙面风压的现场实测数据进行分析,研究了高层建筑风场和其迎风墙面风压的状况及它们的相关性。对台风登陆前后的风场和风压数据进行的分析表明:城市上空强风“天兔”的湍流度不是很大,且随着风速的增加变化比较平稳;迎风墙面各测点的风压系数呈现出中间大、两头小的特征,但差别不是很大。在不同平均风时距和不同来流风向角条件下,对迎风墙面风场和风压的相关性进行了分析。结果表明:当基本时距取为3s 时风场同风压的相关性较高,当来流风向垂直于结构迎风墙面时,风场同风压场的相关性比较大。     

基于CFD方法的主动襟翼控制旋翼翼型涡特性研究

基于计算流体力学(Computational fluid dynamics,CFD)技术,建立了一个适用于旋翼二维主动襟翼控制(Active flap control,AFC)数值模拟的方法。在满足对翼型参数化分析的前提下,使用Euler方程求解以提高计算速度,并采用嵌套网格方法对AFC旋翼后缘襟翼进行运动控制。应用所建立的方法,首先进行了算例验证计算,然后着重对AFC旋翼翼型进行了数值模拟。在此基础上,进一步开展了AFC旋翼翼型主要参数对后缘涡影响的计算分析。结果表明:提高桨尖马赫数、增加后缘小翼摆动频率能加快涡产生速度;而提高桨尖马赫数、增大后缘小翼摆动幅度和后缘小翼长度能增大涡的强度;但增大后缘小翼与主桨叶缝隙间距仅在一定范围内能够增加涡强度。     

基于CFD的直升机旋翼流场及气动力计算

建立了一个基于结构运动嵌套网格的流场求解器,用来精确模拟复杂的旋翼流场,为更好地预测旋翼气动载荷提供一套计算方法.在该求解器中,控制方程为惯性坐标系下的三维非定常Navier-Stokes方程,空间方向上采用低数值耗散的Roe格式结合三阶逆风格式(Monotonic upwind scheme for conservation law,MUSCL),湍流模式采用了一方程的Spalart-Allmaras模型.应用所建立的方法,分别针对Caradonna模型旋翼、UH-60A直升机旋翼的悬停流场和7A旋翼、SA349/2直升机旋翼的前飞流场以及旋翼气动载荷进行了数值模拟.计算结果表明,本文的方法在一定程度上提高了对旋翼流场的模拟能力,进而提高了旋翼气动载荷的计算精度.     

基于生产成本的作业调度

提出了一个基于净现值的调度指标，该指标充分考虑到影响调度决策的各项费用．如在制品库存费用、机床工时费、直接工人的工资、工件提前或拖期完工造成的损失等。然后以遗传算法为工具，研究了作业调度问题，设计了一种很实用的算法。该算法采用基于工序的编码方法把调度编码成染色体，并通过遗传操作来搜索最佳染色体。解码算法根据工序在染色体中的次序和最短加工时间一紧迫度规则得到调度。仿真结果表明该方法是可行的。并有一定的优越性。     

基于非定常气动模型的翼型动态响应计算

基于附着流状态下的旋翼非定常翼型气动模型,给出了计算该状态下的翼型非定常升力、阻力和俯仰力矩的数值计算方法。该方法采用半经验指数响应公式,将响应分为环量和非环量载荷部分,利用数值离散来得到非定常气动响应。应用该方法,分别进行了法向力、弦向力和俯仰力矩等的算例计算,并与可得到的试验数据进行了对比,验证了方法的有效性。同时,还对缩减频率和马赫数对法向力和俯仰力矩的影响进行了计算和分析,得出了一些结论。