Fluid structure interaction simulation of supersonic parachute inflation by an interface tracking method

An Arbitrary Lagrangian–Eulerian (ALE) approach with interface tracking is developed in this paper to simulate the supersonic parachute inflation. A two-way interaction between a nonlinear finite element method and a finite volume method is accomplished. In order to apply this interface tracking method to problems with instantaneous large deformation and self-contact, a new virtual structure contact method is proposed to leave room for the body-fitted mesh between the contact structural surfaces. In addition, the breakpoint due to the fluid mesh with negative volume is losslessly restarted by the conservative interpolation method. Based on this method, fluid and structural dynamic behaviors of a highly folded disk-gap-band parachute are obtained. Numerical results such as maximum Root Mean Square (RMS) drag, general canopy shape and the smallest canopy projected areas in the terminal descent state are in accordance with the wind tunnel test results. This analysis reveals the inflation law of the disk-gap-band parachute and provides a new numerical method for supersonic parachute design.     

Sensitivity analysis of flowfield modeling parameters upon the flow structure and aerodynamics of an opposing jet over a hypersonic blunt body

Implementation of an opposing jet in design of a hypersonic blunt body significantly modifies the external flowfield and yields a considerable reduction in the aerodynamic drag. This study aims to investigate the effects of flowfield modeling parameters of injection and freestream on the flow structure and aerodynamics of a blunt body with an opposing jet in hypersonic flow. Reynolds-Averaged Navier-Stokes (RANS) equations with a Shear Stress Transport (SST) turbulence model are employed to simulate the intricate jet flow interaction. Through utilizing a Non-Intrusive Polynomial Chaos (NIPC) method to construct surrogates, a functional relation is established between input modeling parameters and output flowfield and aerodynamic quantities in concern. Sobol indices in sensitivity analysis are introduced to represent the relative contribution of each parameter. It is found that variations in modeling parameters produce large variations in the flow structure and aerodynamics. The jet-to-freestream total-pressure ratio, jet Mach number, and freestream Mach number are the major contributors to variation in surface pressure, demonstrating an evident location-dependent behavior. The penetration length of injection, reattachment angle of the shear layer, and aerodynamic drag are also most sensitive to the three crucial parameters above. In comparison, the contributions of freestream temperature, freestream density, and jet total temperature are nearly negligible.     

整体弹性结构法及在叶片流固耦合分析中的应用

针对三维叶片时域流固耦合振动响应计算普遍耗时的问题,采用一种假设整体结构的模态位移,求解固体时域响应和实现高效网格变形,发展了一种3维时域流固耦合分析的整体弹性结构方法,并应用于压气机叶片0°和180°相位角的气动稳定性分析中。结果表明,所发展方法的计算结果与传统双向时域算法和文献的结果较为吻合,而计算效率相比于传统算法显著提升;在所分析的两个相位角下,叶片振动的气动阻尼均随流量减小先增大后降低,相比于0°相位角,180°下叶片的气动稳定性大幅提高,表明该方法能有效应用于叶片的工程流固耦合研究。     

基于拓扑优化和形状优化的桨叶结构设计

直升机旋翼桨叶剖面结构设计是研究旋翼动力学设计的基础。提出一种联合拓扑优化和形状优化的 两级优化设计方法,可用于桨叶剖面结构的设计。第一级优化以变密度法(SIMP)作为拓扑优化的材料插值方 法,以桨叶体积最小化为设计目标,约束桨叶节点位移和应力,建立桨叶的拓扑优化求解模型;第二级优化以重 构的桨叶模型为基础开展形状优化,降低局部应力集中以及找到合理的边界节点位置。对优化后的模型进行 有限元分析,结果表明:通过拓扑优化和形状优化的两级优化,能够得到满足强度和稳定性要求的结构布局,为 桨叶结构设计提供指导方案。     

MTPS蜂窝夹芯结构传热性能及热应力分析

对于金属防热结构的蜂窝夹芯结构的稳态情况,基于材料的全灰体假设,同时考虑热传导和热辐射两种传热形式对温度场的耦合作用,利用热流量守恒建立了蜂窝芯层温度场的非线性积分方程。离散化后利用数值方法得到方程组的数值解。对于美国兰利研究中心的实验结果,与本文方法的对比计算结果基本吻合。进一步,利用计算结果讨论了给定面板温度边界情况下,下面板、柱体层的灰度、蜂窝结构长径比对夹芯温度场的影响。并根据温度场和近似的应力分析模型,用半解析结果讨论了稳态情况下蜂窝芯层上的热应力。     

超临界航空煤油喷射到大气环境的喷射特性

用纹影系统对超临界RP-3航空煤油喷射到大气环境中的喷射特性进行了试验研究.在纹影照片中可以看出喷嘴出口下游产生了桶形激波、马赫盘和斜激波,射流近场结构与高度欠膨胀空气射流非常相似.测试结果表明马赫盘的位置随着对比压力的增大而增大,对比温度对马赫盘位置的影响不大;马赫盘的直径随着对比温度的升高而减小,随着对比压力的增大而增大;射流的喷射长度随着对比压力的增大而增大,随着对比温度的升高而减小.      

基于有限元模拟的三维型材拉弯轨迹设计

针对三维型材拉弯,提出一种基于变形控制有限元模拟的轨迹设计方法。首先提取模具长度方向特征线并将其离散成诸多线段单元,这样毛料逐步包覆模具的过程就变成毛料包覆这些线段单元的过程。给定毛料依次包覆线段单元发生的变形,根据切线接触条件（拉弯成形过程中毛料离开模具的位置处两者相切）计算出每步毛料末端位移。将这些位移作为边界条件输入有限元模型中计算毛料应力、应变和回弹。根据成形极限图和规定最大回弹超差量来调整变形模式,找到优化的变形模式和拉弯轨迹。以中空矩形截面型材三维拉弯为例,给出了轨迹设计的详细流程。     

微细电解铣削加工模型及实验研究

 对微细电解铣削加工技术进行了深入研究。将分层加工技术应用到微细电解加工过程中,显著改善了加工稳定性;建立了微细电解铣削加工的数学模型;基于电化学刻蚀原理,在线制得直径小至10 μm的圆柱电极;分组实验并验证了加工模型中各参数如：电极直径、加工电压、电解液浓度、铣削层厚度等对微细电解铣削加工精度的影响。通过优化加工参数,成功加工出了深三角结构和四棱台微型腔,形状精度高,加工稳定性好。     

C/E复合材料网格缠绕结构三维建模与模具设计

通过对圆柱螺旋线和圆锥螺旋线参数方程的推导,实现了网格结构及模具的三维建模。三维模型的建立,可实现结构空间协调,缩短设计周期,节省研制经费,提高生产效率。在此基础上提出了可行的金属模具设计方案,找到了切实可行的模具加工方法,并在工程中加以应用,取得了较好效果。     

MA 及原位-热压合成TiC/ Ti5 Si3 复合材料

以Ti粉、Si粉和C粉为原料,利用高能球磨及热压工艺合成了TiC/Ti5Si3陶瓷复合材料。研究了工艺条件尤其是热压温度对合成产物相组成及微观结构的影响,并结合DSC、XRD和SEM对反应合成机理进行探讨。结果表明:通过优化合成工艺,高能球磨12 h,热压温度1 400℃时,烧结6 h得到了高纯度的TiC/Ti5Si3陶瓷复合材料;合成过程为:反应开始时发生Ti+C■TiC,反应ΔG=-167.72 kJ/mol。2 h时发生5TiC+8Si■Ti5Si3+5SiC,反应ΔG=-62.12 kJ/mol,当6 h时发生3SiC+8Ti■Ti5Si3+3TiC,反应ΔG=-697.8 kJ/mol。显微结构表明:TiC/Ti5Si3复合材料的合成过程伴随Si熔融,该材料以TiC-Si-Ti5Si3形式相结合,其中Si为黏结剂。