基于因子图的状态估计方法运算实时性研究CSCD

对因子图算法及其理论基础进行研究,使用消元算法分析了一元观测因子图模型的稀疏结构,得到了一种适用于该模型的增量推断算法;并对固定滞后平滑算法进行了研究,实现了滑窗在该增量推断算法中的应用;最后通过仿真验证了基于因子图的状态估计方法的估计能力、该增量推断算法在求解一元观测因子图模型的状态估计问题时运算速度的提升能力,以及滑窗维持该增量推断算法运算速度稳定的能力。     

梁壳组合结构的一种非线性有限元计算方法

提出了梁壳垂直组合结构在有限元计算中存在的问题,并以几何大变形问题为依托提供了解决该问题的计算原理和有限元方法。数值算例证明了该理论和方法的正确性。该法不仅适用于梁壳组合结构的大变形非线性计算,也适用于静动力小变形问题的线弹性分析计算。     

非结构弹性体动网格方法应用策略及改进

针对3种代表性的经典弹性体动网格方法，从网格最大变形能力、质量评价指标和空间分布趋势3个方面综合分析了各方法的网格变形特性，重点总结了各方法的适用范围.结合工程实际，研究了特定问题的弹性体动网格应用策略.最后，整合各经典弹性体动网格方法的优势提出了改进方法，并通过NACA 0006翼型的超声速俯仰振荡加以验证.研究结果表明:基于面积的弹性体法适用于加密型网格，而基于距离的弹性体法适用于均匀型网格，该应用策略可为工程问题提供指导.改进方法效果明显，在大变形时相对基于距离的弹性体法的各项网格质量增幅均在10%以上，相对基于面积的弹性体法的最小网格质量增幅几乎达到4倍，并且在数值模拟中具有较高的精度，可显著增强弹性体法的鲁棒性.      

连续波测速系统精度设计与试验鉴定的改进

本文提出一种更符实际的测速系统小增量最优精度设计法,井探讨了对所设计精度进行相应飞机试验鉴定的可能性。计算表明,这种方法适用于新研制的外测系统工程及旧工程的改造。     

风洞实验数据不确定度评估的偏增量法及其验证

本文从全增量近似方法的角度提出了偏增量法，并对多元函数缀积分理论补充了两个推论，进行相关的数学证明，并利用常用数学函数和风洞试验数据处理公式进行了验证。结果表明，本文提出的偏增量法与全微分法在全增量近似上具有相同的一阶精度，结果具有很好的一致性，偏增量法回避了偏导数，相对于传统方法更易于通用和摊广应用。     

紊流风场下的无人机飞行状态卡尔曼滤波估计

为了模拟紊流风场下的无人机飞行状态响应,减小紊流风场对无人机飞行速度的影响,使用数值法建立了紊流风场有色噪声模型,生成了符合大气紊流相关特性的紊流风场并对其进行数值仿真。在此基础上,通过状态扩增处理将系统噪声白化,设计出一种基于有色噪声的无人机飞行速度卡尔曼滤波器,该滤波器解决了基本卡尔曼滤波仅适用于白噪声的情况。仿真结果表明,使用此滤波器可有效减小紊流风场对无人机飞行速度的影响,进而满足飞行速度控制输入的精度要求。     

电源端口电磁脉冲防护方法研究

为了解决电源端口在不同电磁脉冲条件下的防护问题，提出了一种适用于多种电磁脉冲防护设计方 法。该方法能依据不同电磁脉冲波形精准指导防护电路设计；该方法使用CST 对电源端口电源线及电路建模 进行仿真分析；利用其电磁仿真结果，设计防护电路；并采用数值计算方法对防护电路中的防护器件瞬态电压 抑制器进行精准设计，可保证瞬态电压抑制器本身及后级电路不发生损坏，可为电源端口防护不同电磁脉冲提 供理论依据。     

纤维增强及匀质层层叠板动力分析的解析与半解析方法

本文引入复合板条振型函数系列,提出二种层叠板解法。一种纯解析法适用于各种边界条件的矩形层叠板的动力计算;另一种半解析半数值方法适用于任意形状、任意边界的层叠板动力计算,比现有数值计算方法节省大量内存与时间。文中还给出了一系列常用层叠板计算公式。     

高超声速再入飞行器的多孔流动控制

针对高超声速再入飞行器流动控制问题,提出一种头部进气、尾部排气的多孔流动控制概念。采用Fluent软件针对高超声速再入飞行器无孔模型、排气孔模型在马赫数为6,高度为25km状态下开展了三维雷诺平均N-S方程(RANS)数值模拟研究。结果表明:排气孔模型阻力系数略有增加,下排气孔(SH-DB)模型产生正的升力系数增量和低头力矩增量,上排气孔(SH-UB)模型产生负的升力系数增量和抬头力矩增量,实现了无活动气动面产生控制力矩的目的,为高超声速再入飞行器复合控制提供了参考。      

一种适用于跑步状态的惯性/零速/GPS室内外无缝组合导航定位方法

目前，行人导航定位技术已经深入社会的众多领域，受到诸多学者的广泛关注。针对行人跑步状态，研究了一种惯性/零速/GPS室内外无缝组合导航定位方法。首先提出了可靠的、适用于行人跑步零速检测的方法，有效提高了在行人跑步状态下的零速检测的准确性。针对GPS信号容易受到高楼、高架等环境的干扰及在室内容易完全丢失的特点，提出了基于BP神经网络的GPS可用信号筛选方法，提高了GPS信息的可靠性与精准性。在此基础上，研究了基于可变量测的Kalman滤波器，实现了惯性/零速/GPS信息的有效融合，显著提高了在行人跑步状态下的导航定位精度。试验结果表明，所提出的这种适用于跑步状态的惯性/零速/GPS室内外无缝组合导航定位方法的平均定位误差可减小到行人跑步总里程的1%以内。     

用边界元-数学规划法解弹塑性接触问题

 本文基于塑性增量理论导出一种不显含应变的矩阵增量物理关系及平面应变问题σ_z的表达式,提出解决无摩擦弹塑性接触问题的边界元-数学规划法。据此编写了用边界元法解二维问题的通用程序,并用之解算了两圆柱体的弹性和弹塑性接触问题,将数值解与现有结果进行了比较。     

三维弹塑性接触问题的边界元法

将边界元法发展用于解决三维弹塑性接触问题,推导了无摩擦的三维弹塑性接触问题的积分方程,讨论了其数值解技术。数值结果与有限无法的数值结果对比表明,所发展的方法是成功的。     

干摩擦面上接触应力分布的混合分析法

为准确而方便地获得诸如航空发动机风扇（转子）叶片凸肩（叶冠）之间的干摩擦接触面上的应用力分布，将光测弹性实验方法与有限元数值计算方法结合，建立一种求解接触应力的滋合法。用云纹干涉法提取光弹模型接触面的位称作为计算的边界条件，获得接触面上的接触应力。该方法提高了计算方法的精度和效率，也克服了单纯用实验方法求解接触应力的困难和花费大的缺点。应用混合法分析了航空发动机带凸肩风扇叶片在高速旋转时凸肩间接触面上的接触应力，比较混合法一光弹性法所得接触应力分布，两的结果比较一致。     

榫头/榫槽接触问题边界元分析

本文用增量理论建立了二维弹性接触问题的边界元法方程,讨论了求解接触问题过程中的一些具体做法。文中分别对盘/片组件的燕尾型与枞树型榫头与榫槽的接触应力作了数值分析。与有限元及光弹试验结果对比表明,边界元法的结果是令人满意的。      

复合材料机械连接结构分析中的接触算法适用性研究

采用有限元方法计算复合材料机械连接结构孔边接触应力分布时,需要保证其结果的可靠度。分别以Nastran、Marc和Abaqus为平台,研究小滑移接触检测技术和有限滑移接触检测技术以及拉格朗日乘子法、直接约束法和罚函数法等接触约束施加技术的基本原理与特点,结合复合材料销钉连接结构应力分析实验与数值计算,验证各接触算法的准确性和适用性。结果表明:在复合材料机械连接结构有限元分析中,应当采用有限滑移接触检测技术和罚函数接触约束施加方法以提高计算精度。     

用有限元法分析接触问题的一种新的力学模型及计算方法

阐明了用有限元法计算物体接触问题的一种新的力学模型:主动接触体网格中一个节点与被动接触体单元中任意点相接触的模型,放弃了过去常用的节点对模型,为网格划分带来了便利。在这种模型下,给出了粘式、摩擦滑移、开式等接触边界条件。按照修正势能原理导出了对称的接触控制方程,可求解接触体具有非线性的静态/动态的接触问题。本文给出了在一次求解迭代后判断各节点新的接触状态的方法,可进一步决定是否和如何修改接触控制方程。      

Stability analysis of an aeroelastic system with friction

In this paper, harmonic balance method, exact formulation and numerical simulation method are adopted to study the effects of different friction stiffness on the stability of 1.5 degrees of freedom aeroelastic system. On this basis, the expressions of input energy and dissipated energy are deduced, and the energy method is used to reveal the mechanisms of the stable boundary and unstable boundary existing in the system and the effects of different friction stiffness on the stability of the system. Studies have shown that the stability region and the critical aerodynamic damping ratio of the system rise with the increase of the friction stiffness, while the friction stiffness has little effect on the stability boundary. In the analysis of the stability of system, the results of harmonic balance method, exact formulation and Newmark of numerical simulation method are in good agreement. Compared with exact formulation and numerical simulation method, the concept and conclusion of harmonic balance method are simple in the system stability analysis.     

燕尾形榫联接修形的数值研究

将边界元法用于榫联接修形的数值分析研究。对分析具有摩擦接触问题的边界元法的理论及实现方法作了讨论;将边界元分析应用于航空发动机轮盘榫槽与叶片样头联接的计算,并与激光散斑法的实测结果进行比较;对修形的榫联接进行了计算分析,初步揭示了修形榫联接接触边界上力学参数的变化规律,与试验结果相符合。     

接触问题形状优化的边界元法及其工程应用

提出一种用边界元法对平面弹性接触问题进行形状优化的设计分析方法。讨论了平面接触问题形状优化中设计变量的选取、形状描述的方法以及目标函数和约束条件的确定,进行了灵敏度和形状优化的算例考核,并对航空发动机枞树形榫头／榫槽进行了形状优化设计的应用分析。     

Topology optimization in lightweight design of a 3D-printed flapping-wing micro aerial vehicle

Topology optimization is an effective method to obtain a lightweight structure that meets the requirements of structural strength. Whether the optimization results meet the actual needs mainly depends on the accuracy of the material properties and the boundary conditions, especially for a tiny Flapping-wing Micro Aerial Vehicle (FMAV) transmission system manufactured by 3D printing. In this paper, experimental and numerical computation efforts were undertaken to gain a reliable topology optimization method for the bottom of the transmission system. First, the constitutive behavior of the ultraviolet (UV) curable resin used in fabrication was evaluated. Second, a numerical computation model describing further verified via experiments. Topology optimization modeling considering nonlinear factors, e.g. contact, friction and collision, was presented, and the optimization results were verified by both dynamic simulation and experiments. Finally, detailed discussions on different load cases and constraints were presented to clarify their effect on the optimization. Our methods and results presented in this paper may shed light on the lightweight design of a FMAV.