一种新PSO混合算法在直升机配平中的应用

直升机配平计算是动力学分析的基础,其实质是求解高维复杂的非线性方程组。针对经典算法与智能算法的特点与不足,提出了一种求解非线性方程组的新粒子群方法。在粒子群(Particle swarm optimization, PSO)算法的基础上,根据模拟退火(Simulated annealing, SA)思想,引入了嵌入式LM (Levenberg-marquardt)优化 算子。该方法充分发挥了3种算法的优势,克服了LM算法初值敏感性,PSO算法易陷入局部极值等问题。通过UH-60A直升机实例配平计算,验证了本文算例模型的准确性。在此基础 上,针对某一前飞状态下的配平算例,在收敛可靠性和计算效率上通过与其他算法进行对比,表明该算法具有可靠的收敛性和较高的计算效率,进一步验证了该算法在配平问题上的可信度与实用性,为直升机飞行动力学问题的处理提供了一种新的有效方法。     

风力机翼型阵风响应计算与分析

基于网格速度方法,对阵风条件下风力机翼型运动的气动特性进行数值计算。中心格式有限体积法、双时间推进法求解预处理后的非定常N-S方程。计算结果表明:在弦向阵风作用下,翼型的运动能有效降低响应波动峰值,即弦向阵风对运动影响较小;而在法向阵风作用下,气动响应波动很大。计算翼型在不同参数阵风作用下的响应过程,阵风的频率变化对响应影响不大。弦向阵风的幅值增大,叶片运动对阵风响应峰值减缓效果增加;法向阵风的幅值增大1倍,响应峰值也增加1倍。     

压力环境下固体推进剂含损伤的统一本构研究

为了建立固体推进剂压力环境下的统一的损伤本构方程,利用改进的万能材料试验机对一种固体推进剂在常温下、不同的压力环境进行了准静态拉伸实验。利用ABAQUS软件的二次开发用户材料子程序UMAT接口,根据实验得到的结果及其应力应变响应特征,确定了本构方程中的系数。计算结果说明本构方程能够很好地描述该推进剂在压力环境下的力学行为。     

高阶精度DG/FV混合方法在二维粘性流动模拟中的推广

DG/FV混合方法因其具有紧致性、易于推广至高阶及相比同阶DGM计算量、存储量小等优点,已成功应用于一维/二维标量方程和Euler方程的求解。在此基础上,将该方法推广于二维三角形/矩形混合网格上的NavierStokes方程数值模拟,将格式形式精度提高至4~5阶。物理量的空间重构及离散使用DG/FV混合重构方法;无粘通量计算采用Roe格式;粘性通量计算采用BR2格式;时间方向离散采用高阶显式R-K方法或隐式方法。利用该方法计算了有解析解的Couette流动问题以验证几种格式的数值精度阶,并计算了层流平板流动和定常、非定常圆柱绕流问题等经典算例。计算结果表明DG/FV混合方法达到了设计的精度阶,在较粗的网格上亦能得到高精度的计算结果;定性分析和数值结果表明相比同阶DG方法单步计算量减少约40%。     

抛物化 Navier-Stokes 方程(PNS)和高氏 PNS 理论及其应用的评述

抛物化 NS 方程得到广泛应用,已经成为工业标准气动计算的基础。现有的八种抛物化 NS 方程有不同的名称,方程中粘性项的形式略有不同,其中的 PNS 和薄层（TL）NS 方程应用最多。但是这些方程都具有类似的数学性质,例如,当流向方向上马赫数大于1时,他们都是抛物型方程,可以采用空间推进算法（SMA）进行求解。与采用时间推进算法求解的 NS 方程或雷诺平均（RA）NS 方程相比,PNS-SMA 计算降低了空间的维数,节省了大量的存储空间和 CPU 计算时间。PNS-SMA 算法也获得了巨大的进展。但是,早期 PNS 研究在理论上是相当模糊的,高智在1990年提出的粘性／无粘干扰剪切流理论（ISF）弥补了这一不足。ISF 理论概括了 PNS 方程所能描述的基本流动,提出了其流动的运动规律及数学定义式,所导出的 ISF 方程组也属于 PNS 方程的一种。为了不增加新的名称,我们将 ISF 方程组也称为高氏 PNS 理论和方程组。这一理论在 NS 方程和RANS 方程的计算中均有重要的应用。例如,计算最优坐标系的选择以减少伪扩散,网格尺度选择及局部网格加密设计以捕捉高超声速流动中物体表面热流等的急剧变化,壁面压力边界条件的选择以及由高 PNS 导出的壁面判据来进行 NS 和 RANS 近壁数值解可信度评估。本文评述了一些初步的应用,进一步的应用和综合 PNS-SMA,RANS-SMA 以及 PSE-SMA 计算值得深入研究,这里 PSE指抛物化稳定性方程。     

DLR-F6/FX2B翼身组合体构型高阶精度数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS)和剪切应力输运(SST)两方程湍流模型,开展了DLR-F6和DLR-F6_FX2B 2种翼身组合体构型的高阶精度数值模拟,计算外形来自AIAA第三届阻力预测研讨会。主要目的是确认WCNS模拟跨声速典型运输机构型和预测局部构型变化引起的气动特性变化量的能力。在固定升力系数条件下,采用粗、中、细3套网格开展了网格收敛性研究,从气动力系数、压力系数分布、表面流态等方面研究了网格规模对DLR-F6和DLR-F6_FX2B翼身组合体数值模拟结果的影响;采用中等网格开展了来流迎角对2种翼身组合体气动特性的影响研究。通过与National Transonic Facility(NTF)的试验结果和CFL3D的计算结果对比,表明采用高阶精度计算方法得到了网格收敛的数值模拟结果,较好地模拟了DLR-F6翼身组合体局部修型引起的微小气动特性变化和翼身结合部流动特性的差异。     

Aerodynamic design optimization of nacelle/pylon position on an aircraft

The arbitrary space-shape free form deformation (FFD) method developed in this paper is based on non-uniform rational B-splines (NURBS) basis function and used for the integral parameterization of nacelle-pylon geometry. The multi-block structured grid deformation technique is established by Delaunay graph mapping method. The optimization objects of aerodynamic characteristics are evaluated by solving NavierStokes equations on the basis of multi-block structured grid. The advanced particle swarm optimization (PSO) is utilized as search algorithm, which com-bines the Kriging model as surrogate model during optimization. The optimization system is used for optimizing the nacelle location of DLR-F6 wing-body-pylon-nacelle. The results indicate that the aerodynamic interference between the parts is significantly reduced. The optimization design system established in this paper has extensive applications and engineering value.     

多晶镍基合金循环塑性细观本构关系

基于晶体塑性的本构理论,考虑金属材料晶体细观塑性变形流动的各向异性性质及晶体滑移的非线性运动硬化,以Voronoi多晶集合体作为材料的代表性体积单元(RVE),用晶体塑性模拟描述ZSGH4169多晶镍基合金材料的细观本构关系,对ZSGH4169多晶镍基合金进行了晶粒尺度的对称和非对称循环细观分析.通过对称循环数值模拟表明:该模型适合模拟金属材料循环试验中常见的应变硬化现象和Bauschinger效应;通过非对称循环数值模拟表明该模型具有对棘轮效应的描述能力,计算中发现背应力的演化对滑移切应变率有很大的影响,滑移切应变率的数值很大程度取决于背应力演化主导的非线性硬化过程.      

一种基于共轭方程法求解黏性反问题的简化方法

 对于给定压力分布的黏性气动反问题,考虑到壁面微小扰动造成的压力变化主要由势流作用引起,因此可以简化用以获得目标函数对设计变量敏感性导数的共轭方程。将黏性流场插值到粗网格中作为彻体力模型方程的解,则其相应的共轭方程将以简单的源项取代原方程中复杂的黏性项。由于在粗网格中求解,网格数减少,同时收敛速度加快,简化的共轭方程计算时间可以减少到黏性方程的十分之一。典型的算例结果表明,对于附着的边界层流动简化方法计算得到的敏感性导数具有较高的精度,能够有效完成反问题设计且减少总的计算耗时。     

复合材料厚板双轴非线性刚度特性分析

基于非线性本构关系与厚板理论,研究了一种采用增量步迭代法对复合材料厚层合板渐进失效过程进行双轴刚度特性分析的方法。以三维Hashin失效准则为单层板失效判据,当层合板中有一层或多层被检测到失效时,对失效层进行刚度折减,并在增量过程中不断更新刚度矩阵。采用有限元软件MSC.Patran自带的编程语言PCL,编写了计算程序。应用该程序分别对不同材料和不同铺层的层合板算例进行非线性刚度计算,给出了不同的铺层和材料组合下层合板在不同载荷比下的应力-应变曲线,展示了层合板渐进失效过程中应力-应变的变化关系,为复合材料厚板失效分析过程中的刚度变化提供理论与计算依据。