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

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

增广的Davidson算法

Davidson算法通常用来求解对称特征值问题,并被应用到求解大型线性方程组中,取得了不错的效果。但在求解大型线性方程组的过程中,Davidson算法的重新开始过程将影响残量的收敛速度。由此可以考虑在重新开始时保留一些重要的信息,如把最小特征值对应的近似特征向量添加到Davidson算法的迭代子空间中,这样就可以大大加快其收敛速度。本文将在Davidson算法的基础上给出新的算法,即增广的Davidson算法,通过理论分析比较两算法的收敛速度,并给出数值例子加以说明。     

隐格式并行直接求解方法研究

讨论了由CFD问题得到的大型稀疏线性方程组的并行直接法求解问题.介绍了三对角型方程组的SPP算法,将之推广来求解一般的带宽较窄的带状或者稀疏带状线性方程组.介绍一般稀疏的线性方程组的基于排序后再分解的并行求解方法的步骤及其中的一些关键问题.     

大型稀疏复线性方程组双共轭梯度法

有限元复线性方程组的系数矩阵一般具有稀疏性和对称性的特点,全稀疏存贮方法就是利用这些特点,只存贮对称部分的非零元素,采用链表式管理,即节省存贮空间,又便于动态更改.在一般双共轭梯度法的基础上,本文利用广义变分原理对内积进行了重新定义,使双共轭梯度法求解复线性方程组更为有效.数值算例表明这种双共轭梯度法结合全稀疏存贮方案的求解算法在时间和存贮上都较为占优,可靠高效,能够应用于有限元线性方程组的求解.     

超燃冲压发动机流量匹配的临界流量法

为了使用基于积分方程的超燃冲压发动机性能计算模型模拟跨燃、亚燃工作模态,用临界流量法建立用于表征燃烧室热壅塞程度的残差方程,并以隔离段出口静压为独立变量,通过牛顿迭代法迭代求解残差方程,模拟了燃烧室内存在临界截面（即喉道）但不出现热壅塞的实际跨燃、亚燃工作状态。结果表明：该方法能够确定有喉道的超燃冲压发动机燃烧室一维参数分布,并具有计算精度高、计算速度快、收敛性好等优点,可以为基于积分方程的超燃冲压发动机性能计算模型求解提供一定的参考。     

用于二维无粘跨音速流计算的隐式TVD有限体积法(英文)

本文结合Harten的TVD(Total variation diminishing)概念和有限体公式建立一种求解守恒型Euler方程组的有效方法,采用通量分裂的迎风格式对方程组的隐式算子部分离散化,得到系数矩阵是主对角线占优的线性方程组,采用线Gauss-Seidel迭代推进法求解,提高了数值收敛速度,给出了二维跨音速内流的计算结果。     

并行计算中一种非结构网格分割方法

将递归谱对剖分方法应用于流体力学并行计算中的非结构网格分割,以解决负载平衡和最小切割问题。为使用这种方法,计算了网格伴随图的离散Laplacian矩阵的第二特征矢,然后从该特征矢的分量引入网格的对剖分。特征矢计算中应用了Rayleigh商迭代,并进行了一些修正以使收敛强烈地偏向于第二特征矢及考虑逆迭代步中线性方程组的迭代求解。最后,通过非结构自适应网格上Euler方程分区计算的数值结果验证了所发展的网格分割方法。     

一种制导误差分离的新方法

针对制导误差分离模型中环境矩阵S存在严重病态性,从而影响分离结果精度问题,提出了一种基于动力系统求解的制导误差分离方法。该方法从分析线性迭代求解方法入手,将具有病态特性的线性方程组求解问题转化为对相应刚性动力系统的求解问题。这里给出了该方法收敛性及其他特性的证明。为了验证该方法效果,在遥外测视速度误差分别为0.01m/s、0.02m/s以及0.03 m/s的条件下,选用PB(Primary Bayesian,主成分贝叶斯)估计方法与其进行比较,数值结果表明,该方法可有效地降低环境矩阵病态性对误差分离结果的影响,且分离结果的稳健性和精度都优于PB估计方法得到的结果。     

航空发动机数学模型中方程解法与局限性分析

基于部件法建立的发动机数学模型中含有大量非线性方程求解过程,目前所采用的插值法在理论上存在一定的局限,使得发动机性能仿真时常引起模型不收敛.本文首先以发动机部件非定常流计算过程为例,归结给出模型中非线性方程的数学形式及其特点,然后采用实例对解法的局限性进行分析,并得出结论:建议使用插值法求解时应注意两点:迭代初值选取直接影响算法的收敛性;迭代过程采用最近三个近似解来拟合,不一定能真实反映原方程轨迹,会导致迭代发散.     

一种改进的流场计算加速收敛方法

探讨了二维非结构网格无粘流场计算的加速方法。在对残差平均和多重网格这两种主要加速方法进行分析后,提出了一种改进的兼具多重网格思想的残差平均方法。经对绕NACA0012翼型无粘流动计算的验证,收效明显。在几乎不增加计算量的情况下,比原有的残差平均方法可减少迭代步25%～30%.该思想可推广用于三维、有粘流动计算和结构化网格中。     

航空发动机气路故障诊断的平方根UKF方法研究

设计了适用于双轴涡扇发动机健康参数估计的平方根UKF滤波算法,解决了线性卡尔曼滤波器估计结果准确性依赖于线性模型精度;常规UKF算法中由于计算误差及噪声信号影响引起误差协方差矩阵负定而导致滤波结果发散等问题.提出了根据测量残差变化改进滤波收敛速度与稳定性的方法.发动机渐变与突变故障模式下仿真结果表明,平方根UKF估计算法收敛速度快,稳定性强,精度高,是一种有效的发动机气路部件健康参数估计与故障诊断方法.      

薄壁结构双侧异步激光喷丸强化试验研究

航空发动机压气机叶片是典型薄壁结构,通过激光喷丸强化在叶缘引入残余压应力,是提高其抗异物撞击能力,延长疲劳寿命的有效途径。为解决双侧同步强化方法存在的材料层裂损伤风险的问题,提出薄壁结构双侧异步激光喷丸强化方法。试验研究发现:单侧薄壁激光喷丸试验中,光斑功率密度的改变会使薄壁结构相对于激光入射方向呈现"∧"或"∨"两种扭曲变形趋势;另一面采用同样参数进行激光喷丸后,扭曲变形会恢复;合理采用激光喷丸参数,双侧异步激光喷丸强化与同步强化相似,可以在两侧表面均得到残余压应力场,并且扭曲变形能满足形状精度要求。     

起落架缓冲支柱柱塞组件故障分析及改进研究

某起落架缓冲支柱修理时发现柱塞组件的螺母螺纹发生咬合变形,为了探究该故障发生的原因,根据柱塞组件的结构特点与工作原理,构建螺纹咬合变形故障的故障树,基于历史试验数据、螺纹剩余强度分析以及产品状态检查等方法对故障树底事件逐一分析并进行改进研究。结果表明：螺母设计缺陷和螺母未装配到位是导致螺纹发生咬合变形的原因;将螺距改为...     

GH3536电子束焊接残余应力分布分析

采用应力释放法测量了10mm厚GH3536高温合金电子束焊接纵向残余应力分布.结果表明:高温合金电子束焊接残余应力分布范围较小,约为焊缝宽度的4-7倍.残余应力的最大值可达到材料的屈服极限.研究成果为优化电子束焊接工艺和应力与变形分析提供了理论依据和指导.     

存在外部干扰的网络控制系统故障检测和优化

针对存在时延和丢包的网络控制系统,研究了系统的故障检测滤波器设计和优化设计问题.通过泰勒展开将未知时延处理成外界干扰,采用伯努利随机过程来描述丢包现象.设计观测器产生残差信号,将故障检测问题转化为H∞滤波问题.通过Lyapunov函数方法分析了系统的稳定性.以线性矩阵不等式的形式给出了滤波器存在的充分条件和求解方法.为进一步提高系统对故障信号的灵敏度,采用了一种优化设计方法对故障检测系统进行了优化.最后通过仿真验证了算法的有效性.     

基于半随机滤波-期望最大化算法的剩余寿命在线预测

剩余寿命(RL)预测是设备预测维护的关键环节。准确在线预测能够为维护策略的实时安排提供更加精确的技术支持,有效避免失效的发生。工程实际中,反映设备退化过程的性能指标往往不能直接监测,为解决隐含退化过程的剩余寿命在线预测问题,提出一种基于半随机滤波-期望最大化(EM)算法的预测方法。首先以剩余寿命为隐含状态,构建状态空间模型描述直接监测数据与设备剩余寿命间的随机关系。为实现单个设备剩余寿命的在线预测,依据到当前时刻为止的监测数据,采用扩展卡尔曼滤波(EKF)与期望最大化算法相互协作的方法实时估计与更新模型未知参数和剩余寿命分布。最后,将该方法用于惯性测量组合(IMU)剩余寿命在线预测问题,实验结果表明该方法能够提高预测的准确性并减少预测的不确定性。     

圆锥面等测地曲率曲线的轨迹规划方法

为解决复合材料锥壳铺放过程中的纤维残余应变过大的问题,提出了等测地曲率曲线的轨迹规划方法。通过分析发现轨迹的测地曲率决定了纤维的残余应变,于是通过计算圆锥上任意曲线的测地曲率推导出了等测地曲率曲线的表达式,并证明在圆锥展开面上等测地曲率曲线是圆弧。最后,通过计算等测地曲率轨迹的残余应变及其与圆锥母线的夹角分析该轨迹铺放工艺性。结果表明:等测地曲率轨迹相较于固定角轨迹与测地线轨迹具有良好的可设计性和铺放工艺性。     

Generation of Accelerated Stability Experiment Profile of Inertial Platform Based on Finite Element

The residual stress generated in the manufacturing process of inertial platform causes the drift of inertial platform parameters in long-term storage condition.However,the existing temperature cycling experiment could not meet the increased repeatability technical requirements of inertial platform parameters.In order to solve this problem,in this paper,firstly the Unigraphics(UG) software and the interface compatibility of ANSYS software are used to establish the inertial platform finite element model.Secondly,the residual stress is loaded into finite element model by ANSYS function editor in the form of surface loads to analyze the efficiency.And then,the generation based on ANSYS simulation inertial platform to accelerate the stability of experiment profile is achieved by the application of the analysis method of orthogonal experimental design and ANSYS thermal-structural coupling.The optimum accelerated stability experiment profile is determined finally,which realizes the rapid,effective release of inertial platform residual stress.The research methodology and conclusion of this paper have great theoretical and practical significance to the production technology of inertial platform.     

基于Wiener过程的民用航空发动机性能可靠性预测

通过对民用航空发动机性能退化数据的分析,提出了一种有效融合先验退化数据和现场退化数据的性能可靠性评估和剩余寿命预测方法.首先在先验退化数据的基础上确定Wiener过程参数的先验分布,然后利用贝叶斯方法融合新增的现场数据,对Wiener过程参数进行更新,并在此基础上对单台发动机进行性能可靠性评估和剩余寿命预测.该方法能根据现场退化数据不断地对可靠性和剩余寿命进行更新.最后通过某航空公司发动机性能退化数据验证文中提出的方法,结果表明41号发动机在2000循环和3000循环时预测的剩余寿命相对误差分别为0.060和0.018,可以满足航空公司发动机下发计划制定的实际需要.      

典型机身框地板梁缘条裂纹修理损伤容限分析

采用简单有效的方法对复杂的飞机结构进行损伤容限评定具有重要的意义,提出一种简单有效的应力强度因子获取方法,并结合损伤容限分析的一般流程,分析某机身框地板梁缘条含裂纹修补结构的疲劳寿命及使用寿命期内结构的剩余强度。根据机身框地板梁结构受载特点建立简化的分析模型,计算单位载荷时不同长度下裂纹尖端应力强度因子,再由结构边界载荷与应力强度因子的关系确定无量纲应力强度因子;根据损伤容限分析方法编制程序,计算结构在飞行载荷谱下从初始裂纹扩展到临界长度的寿命及各裂纹长度下结构的剩余强度,给出结构检查间隔。结果表明:结构修补后的疲劳寿命及剩余强度均满足损伤容限设计要求。本文给出的损伤容限分析过程及方法可应用于工程中类似结构的损伤容限评定。