Delaunay graph-based moving mesh method with damping functions

The fluid–structure interaction and aerodynamic shape optimization usually involve the moving or deforming boundaries, thus the dynamic mesh techniques are the key techniques to cope with such deformation. A novel dynamic mesh method was developed based on the Delaunay graph in this paper. According to the Delaunay graph, the mesh points were divided into groups. In each group, a factor ranging from 0 to 1 was calculated based on the area/volume ratio. By introducing a proper function for this factor, this method can control the mesh quality with high efficiency. Several test cases were compared with other dynamic mesh methods regarding mesh quality and CPU time, such as radial basis function method and Delaunay graph mapping method.     

基于动态边界控制的优化方法在高升阻比气动外形优化中的初步应用

本文是利用动网格技术和有限体积法来研究基于动态边界控制的气动优化的方法.该优化方法是通过在求解非定常流动的过程中计算优化的目标函数随设计变量变化的梯度.而利用参数化的结构动网格技术能有效、快速地生成与优化过程出现的外形对应的网格并进行计算,具有较高的效率.本文以典型的高升阻比外形优化为例,用非定常NS方程为流动控制方程进行优化,取得了初步的结果.     

基于网格重生成的涡轮叶片变形传递方法

在涡轮叶片多学科优化中需要在气动模型和结构模型之间进行耦合信息(气压、温度和变形等数据)传递,来考虑学科耦合的性质.采用修正几何外形和网格重生成的方法实现了叶片结构变形向气动模型传递.传递时在叶片几何模型上增加若干外形控制线,叶片的外形由这些控制线的形状来决定.根据叶片结构变形,调整控制线节点位置坐标,以修改叶片外形控制线的形状,可以实现变形向叶片几何模型的传递.进一步,重新生成气动网格,可以得到结构变形后叶片的气动网格模型.      

快速实现二次曲面片三角网格剖分的新方法

三维表面的三角网格剖分,对图形学和数值分析来说,都是极其重要的一项预处理技术.利用二维动态、带约束的Delaunay三角网格剖分算法,可以有效地实现二维空间中带边界约束的点集的三角网格剖分.提出了对二次曲面进行快速三角网格剖分的一种新方法.     

复合固体推进剂颗粒的三维网格自动生成

为了分析圆柱颗粒对复合固体推进剂性能的影响,必须对复合固体推进剂中三维圆柱颗粒进行网格划分。根据"力平衡"网格划分原理设计了一种适合于三维圆柱颗粒的网格自动剖分方法。编程运行结果表明,方法能自动完成划分,网格质量较高,满足要求。     

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.     

基于离散等价方程的非结构网格有限差分法

激波装配法把解析关系式嵌入流场避免了间断引起的理论问题，使全场一致高精度的实现成为可能，有望解决超声速流动转捩研究中的感受性模拟难题。但装配复杂激波结构以后，被分割的流场空间常出现不规则的几何形状，这给常规的基于结构网格的有限差分法应用带来困难。基于离散等价方程理论，提出了一种新的基于非结构网格的有限差分法，在空间二维离散点附近仅用3条网格线就可以构造出一阶迎风格式。数值算例表明收敛过程对网格质量不敏感，解决了激波装配法模拟激波相交出现小夹角后使用结构网格进行计算存在的难题。根据这种方法的特点展望了未来的应用前景。     

动载下的圆柱齿轮接触应力分析

本文研究了一般三维弹性接触问题的有限元局部网格密化技术 ,并在此基础上结合有限元—线性规划法首次提出了一种计算直齿和斜齿圆柱齿轮啮合过程接触应力分布的方法。文中以一对航空斜齿圆柱齿轮为例 ,计算了在系统动载和静载作用下该对斜齿轮不同接触线上的接触应力 ,并与航标进行了对比分析。计算结果表明本文方法及其微机用程序系统是快速可靠的。     

贴体与边界正交网格及大迎角涡流计算

本文应用空间理想非均匀线性弹性力学平衡关系生成三维贴体与边界正交网格。通过弹线性的刚度来控制网格的密度，并且自动调整边界节点上外力分布使网格线与边界正交。文中还用欧拉方程有限体积法三步显式时间推进及焓修正方法求解了绕三角翼的大迎角流动。计算结果表明这种方法能够较好模拟跨声速胶体涡流动。     

高亚声速压气机叶片优化设计

为实现压气机叶片的优化设计,采用Hicks Henne函数进行叶型参数化,N S方程流场计算与混合遗传算法结合构成设计软件。以给定叶片表面压力分布为目标,以损失小而扩压度大和给定压比损失最小作为目标,所得优化叶片吸力面等熵马赫数分布合理、符合控制扩散规律,具有较好的压比和损失指标。采用几何方法与椭圆型方程方法结合生成壁面正交H型网格,可提高计算精度和便于采用代数紊流模型的流场计算。     

基于网格参数化变形的单级涡轮多学科可靠性设计优化

基于自由变形技术发展了一种可以协调实现流场、结构分析网格参数化变形的方法,并实现单级涡轮的多学科可靠性设计优化。针对流场分析网格展向积叠的特征,将三维网格变形转化为不同叶高二维截面的网格变形;根据二维流场分析网格的拓扑结构,设计二维控制面,基于叶型设计参数建立叶型控制点、其他控制点伴随叶型控制点移动,实现了网格的参数化变形,避免了网格畸变、提高了网格变形的质量;不同叶高的二维控制面组成控制体,实现三维流场分析网格的参数化变形。采用相同的控制体进行结构分析网格的变形,实现了耦合界面网格协调变形。基于Kriging代理模型进行单级涡轮的多学科可靠性优化,在满足可靠性约束的情况下效率提高4.97%、寿命提高40.86%,证明了方法的有效性。      

涡轮叶片叶面B-spline曲面层次化拟合方法

航空发动机涡轮叶片表面形状复杂,其三维重构是叶片成型精度分析以及质量控制的重要手段。涡轮叶片叶面是一个复杂的自由曲面,其拟合是重构涡轮叶片的关键步骤之一。本文以涡轮叶片叶面拟合为研究对象,提出一种B-spline曲面层次化拟合方法,该方法首先构建网格曲面的多分辨率模型,然后完成该多分辨率三角网格模型的参数化,最后在多分辨率三角网格的基础上通过粗略拟合、精确拟合以及最终拟合3个主要步骤完成曲面拟合。该方法可以完成高阶连续自由曲面的拟合,适用于海量测量数据的快速拟合,具有较高效率并且拟合结果精度高。     

具有孔洞和不规则边界的网格参数化技术

针对凸组合参数化方法难以处理具有孔洞的网格模型,以及对不规则边界的网格模型参数化的边界处变形较大的问题,提出了一种新的参数化方法,首先通过拟合变分隐式曲面修补网格上的孔洞,并构造出较规则的附加边界,最后采用凸组合参数化方法进行参数化.通过曲面拟合实验表明,该方法引入的附加数据作为曲面拟合数据,最终的拟合曲面在孔洞以及边界处具有较好的光顺性,且该方法已应用于复合材料成型模具快速设计中的参数化曲面构建.     

一种适用于飞行器复杂外形的表面与空间网格生成方法

本文通过建立反映机身横截面外形特征的控制点概念，及在控制点间分配网格点数方式，给出了一种非常实用、交互式的表面网格生成方法，所生成的网格失真度小，物体保形好，生成适用于CFD计算需要的表面网络花费的时间少，适用性强。基于无限插值理论，引入本文提出的网格正交控制、物面法向量控制、加权平均光顺措施，给出了一种改进的空间网格生成方法，有效克服了传统无限插值网格生成方法在复杂外形网格生成方面的缺陷，有效改善了网格的法向疏密性、贴体性及周向网格的均匀性，对真实飞行器部件网格生成非常迅速。以某飞行器复杂机身表面与空间网格生成，及无人机翼身组合体空间网格生成为例，检验了本文网格生成方法的实用性和有效性。     

复杂构件MDO六面体网格重构方法

从参数化模型和模型拓扑特征出发,提出了一套六面体网格重构的有限元模型建模方法.同时,使用随动点来约束网格的拓扑边界,保证网格的正交性,确保了网格重构后的网格质量;结合现有的商用软件ICEM的网格生成算法,开发了基于特征的参数化有限元(FE)建模程序.最终,针对航空发动机低压涡轮带冠叶片的叶冠参数进行的多学科优化 (MDO),实例验证了网格重构方法的有效性和可行性.      

非结构混合网格鲁棒自适应技术

数值格式、湍流模型和计算网格是影响CFD数值模拟精度的3个主要因素。结合流场信息的网格自适应技术具备动态优化计算网格的能力，被NASA列为未来CFD发展的一项关键技术。本文针对非结构混合网格，发展了网格单元分布优化、表面网格几何投影和空间网格协调匹配3项关键技术，建立了高鲁棒性几何保真的网格自适应系统。首先，为了提高自适应方法的鲁棒性和通用性，发展了基于标准面网格的多面体网格单元分布优化方法。其次，发展了仅依赖表面网格信息的局部曲面重构技术，采用参数点映射方法实现了新增表面网格点的几何投影，消除了自适应系统对几何CAD系统的依赖。再次，采用改进的距离函数方法实现了空间网格与投影后表面网格的快速匹配。最后，结合基于流场特征的自适应探测器，采用二阶格式的有限体积方法，开展了30P30N三段翼绕流和三角翼大迎角绕流的网格自适应数值模拟。结果表明，通过网格自适应对网格单元的分布进行优化后，流场求解的收敛性和模拟精度都得到了显著提高。     

多块协调变形的网格变形技术及其应用

基于多块网格协调变形的理论,提出了适用于叶轮机内流多物边界的更具一般性的网格变形技术,先获得块边界处的变形量,然后插值得到块内网格各点的变形量。通过二维振荡叶栅和带叶尖间隙的三维跨音风扇叶片的非定常绕流算例,证明网格变形技术在叶轮机气弹稳定性数值模拟中的应用是有效的,能够在保证网格高精度的条件下,提高网格变形速度。     

气动弹性计算中网格变形方法研究进展

网格变形是气动弹性计算中实现计算网格随运动边界变形的主要方法。在总结目前网格变形方法发展现状的基础上,对近几年常用的网格变形方法,即弹簧法、弹性体法、超限插值法、Delaunay背景网格法、径向基函数插值法和温度体法等做了简要的总结。根据各方法构建模型的不同,将它们分成物理模型法、数学插值法和混合方法3类,简要介绍了各方法的基本思想和研究进展,重点比较了各方法的网格变形特性（变形能力、变形质量和变形效率）和优缺点。总结了当前流场边界与结构边界的数据传递方法,对当前气动弹性计算中遇到的网格变形的难点问题作了简要评述并对未来网格变形方法的发展方向进行了探讨。     

过驱动飞行器积分滑模容错控制方法

针对过驱动飞行控制系统中执行器故障时的协调分配问题,提出了一种基于积分滑模的容错飞行控制器设计方法.引入虚拟控制思想,构造了多执行器故障条件下飞行器层级结构控制的数学模型.运用小增益定理,推导了系统闭环稳定的充分条件,并基于线性矩阵不等式,建立了最优状态反馈的凸优化模型.以积分滑模面切换函数为变量,选取李雅普诺夫能量函数,设计了渐近稳定的积分滑模控制器.仿真结果表明,该方法可综合考虑执行器完好和故障时的控制效能,能够实现多操纵面容错飞行控制,具有较好的鲁棒性.