一种基于有限元/边界元耦合方法的管道进口声传播及声辐射模型

本文发展了一种基于有限元和边界元耦合方法的管道进口声传播及声辐射计算模型。该模型将整个声场分为内部有限域和外部无界域,分别用有限元和边界元方法求解控制方程,在二者之间的界面上使用具有物理意义的声阻抗参数进行匹配,并通过一种快速迭代方法实现全声场求解。这种迭代方法可以保证有限元刚度矩阵等带宽以及对称的特性不被破坏,有助于提高计算效率。该模型先得到了Levine and Schwinger标准解的检验,进而在无流动情况下对于简化的航发短舱进口管道模型进行了噪声辐射现象的数值模拟,最后基于计算结果分析了声衬对远场声辐射的影响。     

基于有限元/边界元耦合方法的管道进口声传播及声辐射模型

发展了一种基于有限元和边界元耦合方法的管道进口声传播及声辐射计算模型.该模型将整个声场分为内部有限域和外部无界域,分别用有限元和边界元方法求解控制方程,在两者之间的界面上使用具有物理意义的声阻抗参数进行匹配,并通过一种快速迭代方法实现全声场求解.这种迭代方法可以保证有限元刚度矩阵等带宽以及对称的特性不被破坏,有助于提高计算效率.该模型先得到了Levine-Schwinger标准解的检验,进而在无流动情况下对于简化的航发短舱进口管道模型进行了噪声辐射现象的数值模拟,最后基于计算结果分析了声衬对远场声辐射的影响.      

基于降阶模型和梯度优化的流场加速收敛方法

CFD是解决流体相关复杂工程问题的重要手段，随着流体力学研究和工程应用问题的日益复杂，更高效、精确的CFD算法研究至关重要。提出了一种基于降阶模型（ROM）和梯度优化用于加速稳态流场的收敛方法。该方法首先收集CFD迭代求解过程中的流场快照构造基于正交分解（POD）的降阶模型，然后使用梯度优化方法求解降阶模型得到残差更低的流场，最后以该流场为初值继续进行CFD迭代求解从而加速CFD的收敛过程。将该方法用于基于有限体积法的CFD求解器中，针对无黏、湍流算例验证了所提加速收敛方法的有效性。结果表明，加速收敛方法相比于初始CFD方法迭代步数能显著减少2～3倍，是一种高效且鲁棒的加速收敛方法，且在复杂工程问题中具有广泛的应用前景。     

涡桨发动机螺旋桨实时建模技术

基于螺旋桨片条理论对叶素进行受力分析,推导了螺旋桨拉力和功率等参数的计算公式,建立了螺旋桨实时数学模型,将模型求解归结于干涉角的迭代,并指出模型保证实时性的关键在于迭代算法的收敛速度.通过分析迭代函数及其导数关系,提出一种干涉角初值设置方法,并提出采用割线法代替导数法能加快迭代运算.仿真结果与实验数据对比分析表明:基于叶素受力分析得到螺旋桨拉力和功率的计算精度满足要求,干涉角初值设置以及基于割线法的迭代收敛速度能满足涡桨发动机控制系统实时仿真的需要.      

崎岖地形中关节式月球车姿态估计数值求解方法

 以被动关节式地形自适应月球车为研究对象,融合关节机器人D-H坐标建模方法构建了月球车悬架运动学模型和以侧倾、俯仰和偏转角表示姿态的欧拉角方式建立了车轮到世界坐标系的悬架完整表达模型。利用悬架完整表达模型,用带有迭代因子的连续迭代和离散迭代数值求解方法,建立了光滑和离散崎岖地形姿态估计算法。最后以8轮扭杆双摇杆摇臂月球车原理样机为例来验证此模型求解方法。仿真结果验证了所建立的基于数值方法的崎岖地形被动关节式月球车姿态估计模型的正确性,其求解模型精度能够满足仿真要求。     

基于改进非支配排序遗传算法 的维修资源优化配置

针对在武器装备维修保障过程中出现的资源短缺、资源冲突等问题,构建了多目标的维修资源优化配置模型,并且基于加权思想将多目标转化为单目标以方便求解。针对遗传算法求解多目标优化问题存在的解空间过大、收敛速度慢、计算效率低等问题,提出了基于约束的改进非支配排序遗传算法对资源优化配置模型进行求解。实例分析结果验证了模型及算法的可行性和有效性。     

涡轴发动机高精度实时部件级模型研究

提出了基于神经网络的涡轴发动机共同工作方程求解方法。在基于牛顿-拉夫逊迭代法求解共同工作方程的模型上采集离线训练数据,以共同工作方程迭代求解前的残差为输入,迭代收敛后的共同工作方程猜值修正量为输出,训练BP神经网络,对共同工作方程进行求解。采用变缩放因子的萤火虫算法优化神经网络参数,提高了猜值修正量的预测精度。在飞行包线的某一区域内,采集额定发动机在直升机前飞过程的数据进行神经网络离线训练,并将网络参数代入部件级模型对共同工作方程进行求解,在训练数据采集区域附近的爬升状态、远离训练数据采集区域的前飞状态下进行测试,计算模型输出与牛顿-拉夫逊迭代算法模型输出的偏差,与一次通过算法相比,本文提出方法模型输出最大偏差约为一次通过算法的1/34到1/4,模型运行耗时约为一次通过算法的2/5,验证了算法的有效性。     

基于闭路制导的弹道导弹发射诸元的仿真算法

提出了一种基于闭路制导的弹道导弹的发射诸元算法.给定发射点与瞄准点的地理位置,以压低弹道为目标,通过一种迭代方法求解关机点弹道倾角,并基于地球引力级数展开模型提出了一种虚拟目标的快速计算方法,以主动段弹道平面和被动段弹道平面二面角为指标完成了发射方位角的迭代计算,最后通过仿真算例论证了该算法的可行性.     

机动发射下空间飞行器上升段弹道快速计算方法

机动发射条件下,为提高空间飞行器变射面横向机动模型的诸元计算速度,设计了基于牛顿迭代算法的诸元快速计算模型。以基准弹道诸元参数为迭代初值,选取了6个关键参数,在同一弹道模型基础上,采用牛顿迭代法对诸元进行快速求解。仿真结果表明,在一定范围内,迭代后的变射面横向机动弹道能够与飞行中段弹道高精度交班,且计算耗时不超过2 s,能够满足机动发射需求。     

求解涡扇发动机数学模型的有限域搜索方法

提出一种基于N-R迭代方法的有限域优化搜索法,用于求解涡扇发动机数学模型。该方法不仅具有N-R优化迭代求解特色,还具有将初猜解附近划定域进行泛搜能力,使迭代运算具有强收敛性和平稳性。该方法所定义的迭代因子可通过编程自动调整,进而可以确定最小迭代次数及此意义下的最优迭代因子。      

流体网络节点残量修正算法的改进

为了提高流体网络求解的收敛性与稳定性,对网络节点残量修正算法进行改进。基于理论分析,提出使残量修正算法大范围收敛的初值赋值规律,并通过简单模型的数值分析加以验证;在节点残量修正显式迭代的基础上,提出隐式梯形迭代格式,扩大最大迭代步长。基于改进方法,以某典型航空发动机空气系统为例进行计算,获得空气系统各节点的压力、温度以及各支路的质量流量。该计算过程显示,相对于残量修正显式迭代,改进计算振荡大幅度减弱,收敛速度提高了80%以上。结果表明:改进算法有效地提高了流体网络求解的稳定性和收敛性。     

基于分布式模型预测控制的多UAV协同区域搜索

针对多无人机(UAV)协同区域搜索问题展开研究。提出了一种基于分布式模型预测控制(DMPC)的多UAV分布式优化搜索方法。首先基于传统的搜索图模型,建立了多UAV协同搜索的问题描述和状态空间模型,然后在DMPC框架下,将集中式多UAV在线优化决策问题转化为各架UAV的小规模分布式优化问题,采用基于纳什最优和粒子群优化(PSO)相结合的算法实现对每个子系统优化问题的迭代求解。仿真结果表明:DMPC方法能够有效地降低多UAV协同搜索决策问题的求解规模,是一种可行的方法。     

考虑变量相关性的非概率可靠性灵敏度分析方法

基于不确定性的凸模型描述,研究考虑变量相关情况下非概率可靠性指标及其参数灵敏度求解方法.从非概率可靠性指标的物理意义出发,建立了相关变量和相互独立变量情况下的非概率可靠性统一模型,在此基础上基于Monte-Carlo法,在变量的标准化扩展空间内通过全局抽样得到结构失效域内距离原点最近的设计点,然后通过迭代策略求解得到非概率可靠性指标的收敛解.在此基础上基于差分理论,提出了非概率可靠性参数灵敏度分析方法.数值算例验证了所提出方法的可行性和有效性.     

带非线性支撑的转子有限元模型求解方法

用数值方法研究了非线性支撑的柔性转子系统的动学行为,提出了一种将有限元与非线性支撑结合的模型和求解方法。利用有限元法(FEM)构建转轴和转盘部分的模型,通过矩阵进行组合;利用离散元方法对包含滚动轴承和挤压油膜阻尼器(SFD)的支撑部分进行建模,此部分包含4个单元,分别为轴承内圈、外圈、SFD内圈和支撑鼠笼。有限元部分和离散元部分通过轴端节点相连,仿真过程中轴端位移传递给非线性支撑部分,支撑部分通过位移计算得到的非线性力反过来作用于有限元转子轴端部分。为了耦合求解有限元转子和非线性支撑组成的数学模型,提出了一种综合的迭代求解方法,克服传统的有限元求解方法对轴端隐性非线性支撑的求解局限性。由于转轴部分采用了Timoshenko梁单元建模,对比与简单转子模型,可以考虑陀螺力矩和轴的柔性特征,更能体现非线性支撑对振动真实影响。在建立的20个轴单元的有限元转子模型中,非线性响应更多体现在靠近非线性支撑的节点1和节点21处,响应频谱中靠近轴端的节点能体现出滚动轴承的2倍和3倍变柔振动频率。     

基于ADE-ELM的涡轴发动机建模方法

提出了基于自适应微分进化-极端学习机(ADE-ELM)求解平衡方程的高精度涡轴发动机实时部件级模型建立方法.基于牛顿-拉夫逊(N-R)迭代模型,以迭代计算前模型平衡方程残差为输入,迭代收敛后平衡方程猜值修正量为输出,训练极端学习机,并采用自适应微分进化(ADE)算法优化极端学习机(ELM)参数,提高猜值修正量映射精度.ADE算法中采用sigmoid型自适应缩放因子,提高了微分进化算法的寻优能力.在涡轴发动机不同飞行状态下的测试结果表明,以N-R迭代算法模型为基准,基于ADE-ELM的发动机模型,最大建模误差约为一次通过算法的1/3,运算耗时约为一次通过算法的1/3,验证了算法的有效性.      

复合材料盒段结构屈曲稳定性分析及优化技术

随着计算机科学和有限元技术的发展,屈曲稳定性问题有限元数值求解技术已经比较成熟,但是在大型飞机结构工程应用中还是由于计算量大、收敛困难而受到限制.特别是在需要反复迭代计算的优化过程中,更是受到该问题的困扰.针对飞机机翼结构中的典型盒段结构,本文研究了基于ANSYS子模型的考虑线性屈曲稳定性约束的优化方法,以及非线性数值计算技术.研究表明,在相同精度要求的条件下,位移求解模型所需网格尺寸可远大于屈曲稳定性模型.本文利用验证结论,采用网格尺寸100mm建立粗模型进行位移求解并确定危险部位,采用网格尺寸10mm建立子模型进行屈曲稳定性求解,从而完成全结构的优化设计.此方法计算效率高,在非线性屈曲稳定性求解以及优化迭代计算中优势明显.     

一种扑翼飞行器气动参数的辨识方法

针对实际飞行中无法直接测量的扑翼飞行器气动参数辨识问题,结合刚体六自由度模型,提出一种基于迭代学习和人工鱼群法的扑翼飞行器隐式气动参数辨识方法。鉴于扑翼飞行器飞行试验中待辨识气动参数值难以直接测量、导致一般辨识算法中梯度难以求解的问题,提出基于摄动法的梯度方向寻优方法。考虑到待辨识参数数量及辨识结果对参数初始值的敏感性,该方法采用人工鱼群算法优化计算待辨识参数初始值。且针对迭代过程中损失函数易陷入局部最优和优化速度受限问题,采用变学习因子迭代学习策略。试验结果表明,所提出的算法能有效估算出扑翼飞行器气动参数。     

喷管中泡状气液两相流影响因素分析

为了获得边界条件对缩放喷管中泡状流的影响,采用包含气泡动力学的连续均质混合流模型,基于求解Rayleigh-Plesset方程的变步长龙格-库塔方法,提出一种稳态泡状流求解的迭代修正计算方法,实现了喷管中流场计算。获得了气泡半径、混合密度在喷管中变化情况。将不同边界条件下的计算结果进行对比分析,得到了气泡平衡半径、喷管入口含气率及入口速度对流动的影响。     

基于SOA-Newton迭代的六自由度平台正解算法

运动学正解是实现六自由度平台高精度控制和三维可视化仿真的基础,但是传统方法存在着求解难度大、精度低的缺陷。针对Newton-Raphson迭代法求解时依赖迭代初值的问题进行了研究,提出了一种基于SOA-Newton迭代的六自由度平台正解混合算法。该混合算法充分发挥了SOA算法的群体搜索性和拟Newton法的局部细致搜索性,同时也克服了粒子群算法后期搜索效率降低和Newton法对初始点敏感的缺陷。以研制的车载位姿平台为例,在Simulink软件中建立了混合求解算法的仿真模型,给出了计算实例,并与单独的Newton-Raphson迭代法和SOA算法进行了对比。对比结果表明,SOA-Newton混合算法具有极好的稳定性和较高的收敛速度及精度,更能满足工程实际应用。     

RANS方程和附面层方程耦合求解转捩位置的方法

研究了基于线性稳定性分析的RANS（Reynolds-Averaged Navier-Stokes）方程和附面层方程耦合求解转捩位置的方法。RANS方程计算得附面层外边界的速度分布作为附面层方程计算的输入,求解附面层方程得到稳定性分析需要的附面层参数。由于附面层的解在分离点处的奇异性,用交互式附面层方程求解方法。用基于线性稳定性分析的e^N方法计算转捩位置,并考虑了Tollmien-Schliching波和层流附面层分离造成的转捩。迭代求解RANS方程和附面层方程,直到转捩位置收敛。RANS求解中使用转捩过渡区模型,避免了点转捩模型引入的数值扰动。通过对NA-CA0012翼型的计算并和实验结果及Xfoil计算的结果进行比较,吻合较好,有较好的工程实用价值。