三角翼大迎角绕流特性数值模拟的网格处理技术研究

研究网格对三角翼大迎角绕流特性数值模拟的影响.结果表明,网格的生成需要结合流动现象,对网格拓扑结构和网格点分布进行选择与搭配.C-H型网格适宜模拟尖前缘分离涡流态,法向网格在一定范围内应等距增长,沿流向逆压梯度较大的区域内适当增加网格点,尾迹区网格则应做上翘处理.Euler方程具有模拟三角翼旋涡及预测涡破裂特性的能力,但对二次涡等粘性引起的流动细节把握能力不足.利用层流假设的N-S方程,通过合适的网格,也可得到满意的计算结果,但对涡破裂后的强烈非定常湍流流动模拟能力不足.采用旋涡螺旋度可准确反映主涡与二次涡流动,描述旋涡的破裂现象.用轴向速度迅速减小并小于来流速度的点作为涡破裂判据似应更合理.     

带压缩性修正的离散涡方法在结冰数值模拟中的应用

离散涡方法（DVM）是一种无网格的涡运动算法,适用于解决易产生分离流的非定常问题,将其应用于结冰过程中的流场求解,在有效模拟分离流动的同时能避免冰形尖角对网格质量的影响。但离散涡方法基于不可压 N-S 方程,无法应用于预测可压缩流动下的结冰过程。本文在离散涡方法的基础上添加普朗特—格劳尔特压缩性修正,进行基于离散涡方法的可压缩流动下的数值模拟,并将其应用于翼型结冰预测;对流场分布、结冰冰形和结冰模型计算过程等仿真结果进行对比和分析。结果表明：引入压缩性修正后的离散涡方法能较好地模拟可压缩流动,与实验值相比,基于该方法得到的结冰数值模拟结果符合良好,对结冰数值模拟在工程上的应用提供了一定的参考。     

低雷诺数翼型层流分离现象大涡模拟方法

为准确预测翼型低雷诺数条件下出现的层流分离流动现象,发展了基于结构化拼接网格的大涡模拟方法。控制方程为Favre滤波Naver-Stokes方程,并选用多种亚格子模型。空间离散采用AUSM格式以及高阶WENO格式,时间推进采用显式方法和隐式方法。以SD7003翼型在雷诺数60 000及4°迎角下的层流分离流动为研究对象,对比分析了数值格式、亚格子模型、网格尺度对流场预测结果的影响。数值结果表明:划分的计算网格能够有效解析小尺度流动结构,基于隐式亚格子模型、采用AUSM格式和双时间步推进的大涡模拟方法能够准确预测流动分离、转捩、再附等复杂流动现象,计算得到的平均压强系数与雷诺应力分布与文献数据吻合较好。与转捩模式计算结果对比进一步表明:发展的大涡模拟方法能够准确预测翼型低雷诺数层流分离流动非定常演化过程,为下一步的研究工作建立了有效的数值模拟手段。     

圆柱突然起动的快速涡解法

 文中将Chorin的涡方法、格子涡(VIC)方法和快速Fourier变换(FFT)相结合,发展了一套快速有效的求解高Re数下非定常N-S方程的算法(快速涡方法)。通过对二维圆柱突然起动问题的数值计算,得到流动分离及其演化过程(包括二次涡及二次涡对的生成);此外还得到回流区对称轴线上速度分布。这些结果与实验数据及其它数值解相比,是令人满意的。该程序对其它高Re数非定常流动问题也是适用的。     

网格生成和数值模拟的讨论

讨论了多种网格生成方法及其数值模拟。给出一种生成多部件复杂外形飞行器分区网格和Ｅｕｌｅｒ方程分区解。讨论了前缘分离涡为主体的旋涡流场。作了跨声速流分区的并行计算的初步尝试。给出一种二维非结构网格的生成方法。     

基于结构网格的栅格翼绕流数值模拟

针对栅格翼这种复杂的构造形式,采取有效的多块结构网格生成策略,对栅格翼导弹生成高质量贴体的结构网格。在此基础上,通过求解N-S方程,对栅格翼导弹亚、跨、超声速绕流流场进行了数值模拟。计算结果与风洞试验结果吻合较好,能够准确地捕捉到栅格框内复杂流场的流动特性,从而结构网格生成策略与计算方法得到验证,为栅格翼导弹绕流流场的精确数值模拟提供了技术途径。     

尖拱弹身横向喷流数值模拟

本文通过数值求解层流Navier—Stakes方程来模拟超声速流场中横向超声速喷流的干扰流场，采用分块对接网格和“O”型网格技术，精确模拟喷口截面形状并生成高质量的贴体计算网格，对尖拱弹身的横向超声速喷流的干扰流场进行了数值模拟，研究了喷口附近的流场结构、涡系结构、波系结构及流动分离等流场特性，得到喷流流场的干扰效应，计算结果与实验值吻合较好。     

叶轮机械中间隙流与通道二次流相互作用的数值研究

采用数值解法对一压气机静叶栅中叶顶泄漏流与通道二次流的相互作用进行研究 ,详细揭示了在顶部区通道涡、泄漏涡及分离涡的生成和发展。结果表明 ,由于泄漏流动与通道二次流动的方向正好相反 ,适量的泄漏流动可以在一定程度上起到削弱通道涡强度、改善顶部区流动状况的作用。     

分离涡模拟方法在内埋弹舱机弹分离模拟中的应用

为了研究内埋弹舱非定常流场对武器投放分离特性的影响,基于分离涡模拟（DES）方法和非结构重叠网格技术,建立了内埋武器机弹分离高精度数值模拟方法。在此基础上,首先通过典型空腔流动和机弹分离算例验证了数值方法的有效性,然后采用平均雷诺应力模拟（RANS）和DES两种方法探究了机弹分离过程中弹舱腔内涡结构的演化历程,分析了两种方法在机弹分离流场结构模拟中的差异,研究了内埋弹舱流场结构对导弹分离特性的影响。计算结果表明：机弹分离过程中,高强度涡结构逐渐向弹舱腔体后缘聚集,同时弹舱剪切层流动结构被导弹破坏,导弹出弹舱阶段俯仰方向所受气动力矩急剧变化。相比RANS方法,DES方法捕捉到的高强度小尺度涡结构对导弹压力分布影响较大,使得分离后期导弹俯仰角、俯仰角速度产生明显差异,DES方法更适合内埋武器舱分离特性精细化研究。该研究结果可为内埋武器机弹分离研究的计算方法选择和安全分离策略提供一定的参考依据。     

基于混合笛卡儿网格方法的可压流动数值模拟

以可压缩黏性流动的数值模拟为研究背景,发展了一套自适应混合笛卡儿网格(AHCG)方法以及基于有限体积方法的雷诺数平均Navier-Stokes(RANS)的数值求解方法.为更好地模拟边界层的黏性流动在近壁面处采用贴体结构网格,剩余计算区域自动生成与之相重叠的笛卡儿网格,并同时发展了基于流场特征的笛卡儿网格自适应技术.结合ADT(alternating digital tree)算法显著减少了网格生成中“挖洞”和“贡献单元”搜索的消耗机时,50万左右的网格数目下,搜索耗时为0.062s,仅为普通遍历方法的1/1847.通过二维圆柱与两段翼型绕流的数值算例显示,定常AHCG方法能够准确地预测物面压力分布与升阻力系数并且具备处理复杂外形的能力;同时通过二维非定常圆柱绕流问题与NACA0015矩形机翼翼尖尾涡的捕捉算例显示,结合了动态自适应网格加密的非定常AHCG方法尤其适用于旋涡主导流动.      

斜爆轰波并行数值模拟的化学加速算法性能研究

考虑基元反应条件下的爆轰波精细结构的数值模拟计算量巨大,发展高精度和高效率的计算方法十分必要。本文以Ma马赫数为7的H2/O2/N2预混气来流形成的斜爆轰波为数值模拟对象,研究了一种基于并行计算架构的用于加速化学反应计算的建表算法的计算性能,考察了不同建表策略和不同并行分区数量对算法性能的影响。研究结果显示,本文采用的建表算法能够很好地再现斜爆轰结构,其计算精度不受建表策略和并行计算分区数量的影响；而算法的计算效率则取决于不同分区对应的数据表之间操作的同步性,其中,数据表中节点数据的取回率和数据表设定的尺寸上限都会影响数据表操作的同步性。本文采用的TP和PLP两种建表策略的计算结果表明,TP策略数据表中节点取回率高于PLP策略,故计算效率更高；而计算分区数量越少,则分区对应的数据表尺寸上限越大,数据表的同步性就越好,计算效率也越高。上述研究为发展适于气相斜爆轰波数值模拟的高保真和高效率的计算方法提供了有益参考。     

天津空港国际物流区的功能和发展模式研究

天津机场发展的关键在于天津空港国际物流区的建立，它是天然物业发展的引擎，是提高天津形象的窗口，有利于推动天津产业发展和产业的转型，空港国际物流区发展的模式是自由区运行模式，空港和海港齐飞的空间模式与首都国际机场联盟发展的货客分流模式。     

基于复合地标导航的动平台四旋翼飞行器自主优化降落技术

针对四旋翼飞行器在依靠地标导航完成动平台自主降落任务中存在的目标易丢失、地标定位死区大和降落可靠性差等问题，设计了一种基于复合地标导航的动平台四旋翼飞行器自主优化降落系统。该系统以圆环地标和二维码构成复合地标来解决仅用单一地标定位时存在的定位死区大和定位范围小等问题。针对地标识别丢失、动平台车轮打滑和码盘标定不精确等问题，建立动平台的精确运动模型同时考虑码盘包含未知测量偏差，基于扩展卡尔曼滤波器实现了对动平台连续位姿的在线估计。最终，基于动平台位姿估计结果以最小Jerk指标设计降落轨迹和降落策略，实现了四旋翼飞行器在动平台上高效平稳的降落。为验证所提系统的有效性，设计了仿真和实际降落实验，验证了所提复合地标实现摄像头距离在0.5~6.0 m内的综合完整定位；所设计的动平台状态估计器能在码盘存在未知测量偏差的情况下准确估计出平台的实时位姿，同时所提自主优化降落策略和轨迹规划算法保证了可靠的动平台降落。     

基于效能的先进战斗机航电系统动态重构方法

为满足未来先进战斗机全战斗过程的对敌压制能力需求，分析了作战任务与航电系统支持能力间的关系，建立了"作战任务-航电能力-资源需求"间的关系矩阵和航电系统效能模型；以最大化全飞行阶段航电功能整体效能和飞行安全为目标设计了针对不同作战场景的航电系统动态重构策略及重构流程；通过数值分析，对比了动态重构航电系统与静态配置航电系统在不同作战区边界的效能，表明动态重构特性能有效提高战斗机各阶段的作战效能，提高有限资源条件下的阶段优势。     

面向未知区域深度测量的序列图像稠密点特征生成算法

对未知着降区平坦度测量是无人机在复杂地形下安全着陆的关键问题。首先,根据小孔成像原理推导出基于单目序列图像的未知区域深度计算方程;其次,针对稀疏匹配存在深度信息重构误差大而稠密匹配在平滑区域误匹配率高的问题,提出一种基于Delaunay三角剖分的稠密点特征生成算法;然后,分别对序列图像中的2帧图像提取亚像素级Harris角点和尺度不变特征变换(SIFT)特征点,并分别进行特征点匹配;再以2种特征点间的欧氏距离作为约束条件将2种特征点进行融合,生成准稠密特征点;最后,将准稠密特征点进行Delaunay三角剖分,并根据每个剖分三角形上3个顶点像素偏差的方差值制定稠密特征点的生成策略,并结合所提出的深度计算方程计算整个未知区域各点的深度信息。通过Vega Prime(VP)搭建仿真演示验证系统,实验结果表明在机载相机距地面400m处计算高度分别为90m和55m的物体深度信息时,其深度测量相对误差不超过0.89%,具有较高的精度。     

Multi-UAV coordination control by chaotic grey wolf optimization based distributed MPC with event-triggered strategy

The paper proposes a new swarm intelligence-based distributed Model Predictive Control (MPC) approach for coordination control of multiple Unmanned Aerial Vehicles (UAVs). First, a distributed MPC framework is designed and each member only shares the information with neighbors. The Chaotic Grey Wolf Optimization (CGWO) method is developed on the basis of chaotic initialization and chaotic search to solve the local Finite Horizon Optimal Control Problem (FHOCP). Then, the distributed cost function is designed and integrated into each FHOCP to achieve multi-UAV formation control and trajectory tracking with no-fly zone constraint. Further, an event-triggered strategy is proposed to reduce the computational burden for the distributed MPC approach, which considers the predicted state errors and the convergence of cost function. Simulation results show that the CGWO-based distributed MPC approach is more computationally efficient to achieve multi-UAV coordination control than traditional method.     

断裂问题研究中裂尖塑性区的影响及处理原则

通过对K判据的适用区与裂纹尖端塑性区的比较,在塑性区一定要小于K判据适用区的条件下,给出了K判据的适用边界线和判断脆性断裂的判断式;确定了K判据在结构静强度计算中塑性修正的基本原则,并在此基础上讨论得出了塑性区尺寸小于裂纹长度,裂纹扩展寿命的估算方法仍适用,塑性对裂纹扩展寿命估算影响不大,而选用不同的经验公式估算裂纹扩展寿命所出现的误差却很大的结论。     

基于ARLEQUIN方法的复合材料螺栓连接失效模拟

将ARLEQUIN方法应用于螺栓连接件的挤压失效模拟中，将连接件层合板分为孔边局部损伤区、无损伤区和过渡区。局部损伤区精细化划分网格，无损伤区采用粗网格划分，过渡区通用ARLEQUIN方法将损伤区和无损伤区进行耦合。采用用户自定义单元在ABAQUS中引入ARLEQUIN方法，实现了基于ARLEQUIN方法的复合材料螺栓连接的三维有限元失效模拟。建立了传统有限元模型和ARLEQUIN模型，并与试验结果进行了对比，表明在保证与传统有限元模型相同计算精度的情况下，采用ARLEQUIN方法进行复合材料螺栓连接失效模拟大大降低了计算资源，具有更高的效率。     

尖楔结构低速高/中温双路气流组合热试验方法

在典型高超声速飞行工况下,数值模拟分析了高温合金尖楔前缘结构沿气流方向大温度分布梯度将带来严重的热强度问题,产生大温度分布梯度的根本原因是尖楔结构头部区域平均热流密度与后段平板区域平均热流密度之差,而受头部区域热流密度具体分布的影响不大;进而提出了一种低速高/中温双路气流组合热试验方案,并通过数值模拟方法证明了该方案具有两股气流参数可以独立调节分别满足尖楔结构头部驻点区域及后段平板区域大、小两种热流密度的优点,进而解决单喷口低速高温燃气流热模拟试验难题,满足尖楔结构高超声速飞行工况下大温度梯度模拟要求.同时,该方案通过高/中温气源的合理组合搭配可以大大降低尖楔结构热试验所需高温气源发生功率,推广应用于电弧风洞可拓展其热试验范围.      

AF1410钢电子束焊接接头组织及性能研究

采用电子束扫描焊和修饰焊工艺对5mm AF1410钢板进行电子束焊接。针对热处理前、后的接头,采用光学显微镜、扫描电镜、硬度仪等对其组织、硬度及拉伸断口进行分析。结果表明:接头热影响区分为浅腐蚀区和深腐蚀区,分别为单一马氏体(M)、M+少量逆转奥氏体(Ar)组织;焊缝区柱状晶分为重熔区、正火区和回火区,晶内为马氏体,晶界附近为残留奥氏体。热处理后,接头组织、显微硬度与母材趋于一致,抗拉强度也达到了母材98%以上,失效裂纹始于热影响区,呈约60°方向扩展并贯穿整个接头。