基于伴随方程和自动微分的雷达散射截面表面灵敏度计算

反隐身技术的发展对军用飞行器的隐身性能提出更高要求，针对经典电磁场伴随方法无法获得表面灵敏度，难以为优化设计提供直观指导，且梯度求解时需反复填充阻抗矩阵，当设计变量、入射角度较多时梯度求解效率下降的问题，提出了基于伴随方程和自动微分的雷达散射截面表面灵敏度计算方法，并结合基函数特点提出了表面灵敏度的稀疏存储和稀疏矩阵-矢量相乘方法，避免了直接计算表面灵敏度时计算量、存储量无法承受的问题。提出的表面灵敏度计算方法可以通过一次阻抗矩阵偏导数求解得到所有网格点的表面灵敏度，避免反复填充阻抗矩阵的问题。当入射角度改变时，求解任意数量设计变量梯度的计算量约为8次矩阵-向量乘，显著提升基于伴随方程的梯度计算效率，为作战飞机气动隐身一体化优化设计提供有力技术支撑。     

桥联重积分方程一种新的数值方法

纤维增强复合材料断裂力学中的桥联问题常常可归结为一个二重积分方程。常见的求解方法有矩阵迭代法和积分迭代法，但这两种方法都产复杂。本文给出一种新的分段插值方法，对此积分方程进行求解。数值计算结果表明，该方法效果令人满意。     

直升机横侧稳定性试飞数据计算法

运用小扰动理论，通过在力矩方程略去姿态项及横侧受扰运动方程组的线化和降价处理，使横侧运动方程组用标准特征矩阵表示，从而转化成求解特征矩阵的特征根问题，并利用Ｚ１１直长操稳试飞数据，用特征矩阵法对Ｚ１１直升机横侧稳定极进行了求解，所得结果与费拉雷法比较，结果是一致的，这表明该方法正确，可靠，完全可用于直升机横侧稳定性的飞行试验研究。     

GMRES算法在悬停旋翼数值模拟中的应用

发展了一种基于广义极小残差(GMRES)算法的悬停旋翼数值模拟方法,并对GMRES算法中矩阵与向量乘积的两种计算方法进行了分析和讨论。应用该方法在旋转坐标系中采用非结构混合网格和格点格式有限体积法对以绝对速度为变量的欧拉方程进行了直接求解,其中对流项的离散应用了基于Roe的Riemann近似解的迎风格式。对Caradonna-Tung旋翼跨声速悬停流场进行了数值模拟,计算结果与相关实验数据吻合较好,并与LU-SGS方法进行了对比,表明GMRES算法可以有效地加速流场的收敛,提高计算效率。     

基于流场/声爆耦合伴随方程的超声速公务机声爆优化

基于自主研发的大规模并行结构化网格CFD求解器PMB3D以及并行化伴随方程求解器PADJ3D，开展了流场/声爆伴随方程的求解研究。首先采用标准算例，对内部CFD代码PMB3D软件和声爆预测代码进行了声爆计算可信度验证，以及声爆强度对近场声压梯度的校核。针对并行环境下多块对接网格的近场声压提取操作的复杂性，提出了"包围盒"的方法实现并行环境下近场声压装配单元编号、网格块编号以及对应的进程编号确定，基于声爆计算坐标将并行传递的数据进行一维排序，为声爆预测、伴随方程以及梯度求解提供输入条件。通过线性插值雅克比矩阵实现均匀坐标系梯度信息向非均匀坐标转换，并进一步根据结构化网格特征提出了插值原则，简化了近场声压转换雅克比矩阵的变分。通过装配单元记录，实现声爆强度对流场守恒变量的变分结果向各个进程装配，将装配结果作为流场伴随方程的右端项实现流场声爆耦合伴随方程的求解。此外，对小型超声速公务机开展了声爆优化，对比分析了设计前后的声压及其频谱特性。     

全隐式TVD格式求解二维Euler方程的对角化算法

本文采用有限体积法,Chakravarthy提出的基于近似Riemann解的全隐式TVD差分格式求解二维非定常Euler方程。采用牛顿法对半离散化后的全隐式差分方程进行线性化,然后对线性化后的方程采用近似因式分解法处理,得到了两组三对角块矩阵方程组。在求解该方程组时,采用了矩阵对角化思想使原来的4×4三对角块矩阵方程组转化为4组分离的三对角代数方程组,因而大大节省了求解方程组所耗费的计算时间,提高了计算效率。通过对翼型的跨声速绕流问题及圆柱超声速流动问题的计算,证实了本文算法不但简便可靠,而且还具有较强的通用性。采用TVD格式捕捉到的激波分辨率高,上、下游无任何波动。     

基于LS-DYNA滚动轴承使役可靠度计算方法

利用软件LS-DYNA对滚动轴承使役过程进行了动态有限元计算,基于响应面法（RSM）建立6408滚动轴承输入参数与输出响应的函数关系方程,设定各类随机误差分布,利用矩阵试验法确定抽样点在输入变量抽样空间的位置,进行一系列确定性拟合试验.每次试验借助LS-DYNA进行精确求解,运用求解数据对方程进行最小二乘法回归分析确定方程的组成项及系数.最后利用函数关系方程计算轴承任意时刻的使役可靠度.结果表明:计算得出该轴承在0.020008s时的使役可靠度为97.062%.      

基于遗传算法的多振动台随机振动控制方法

 针对多振动台随机振动控制问题中求解补偿矩阵的广义逆算法,提出利用遗传算法对补偿矩阵的初值进行优化,解决了该类非线性、大范围极值求解问题。双振动台随机振动控制的仿真实验研究表明:将基于遗传算法的广义逆求解补偿矩阵方法引入随机振动控制过程,在同等精度条件下,计算时间大为缩短并改善了控制中共振峰处的奇异性,使广义逆算法在多振动台随机振动控制中的在线应用成为可能。     

基于气动力分段表达的颤振求解方法与模态跟踪技术的应用

在频域颤振方程的求解中,PQI法采用分段气动力的表达形式将颤振方程进行分段直接求解,避免了迭代。由于广义气动力分段插值函数PQI的构建依赖折合频率的选取,将构造的分段插值函数嵌入颤振方程中求解存在模态窜支问题,对此发展了一种基于特征值摄动理论的预测跟踪法来确定各阶模态对应的特征值随参数变化的前后对应关系。文中介绍了PQI颤振分析的基本方法,并采用P-K法进行对比分析来验证算法的计算精度,最后将发展方法应用到某飞机全机轴对称状态颤振分析中来验证算法在实际工程颤振分析中的有效性,数值结果表明该方法计算精度高,模态跟踪效果良好。     

PI型广义预测自校正控制器的简化算法

基于具有PI结构的广义预测自校正控制器,推导了其不需在线求解丢番图方程的算法,从而使计算量大大减小,便于该控制器的在线实现。仿真结果已表明了该算法的快速性。     

圆轴承三维热流体动力润滑的研究

联立求解了广义雷诺方程、三维能量方程、三维固体热传导方程及载荷平衡方程 ,并考虑粘度随温度及压力、密度随粘度及压力的变化 ,在油膜与轴瓦界面使用热流连续性边界条件 ,得到了圆轴承油膜及轴瓦的三维温度场。结果表明最高温度并不出现在油膜与轴瓦交界面 ,不同圆周位置沿膜厚方向的最高温度出现位置不同。证实了考虑三维温度得到的轴承的润滑性能与考虑二维温度得到的轴承润滑性能有很大的差异。采用Newton-Raphson法求解广义雷诺方程和载荷平衡方程 ,采用快速扫描法求解由有限差分离散的能量方程和热传导方程 ,在计算中发现所使用的算法收敛速度快 ,大大节约了计算时间。     

连续广义系统的镇定性分析

本文主要研究奇异系统的稳定性和镇定性问题，通过Lyapunov议程2和Ricaati方程，给出了奇异系统的渐近稳定性和镇定性的判据。     

用于失谐叶盘动力学特性分析的减缩计算方法

提出了一种用于失谐叶盘动力学特性分析的减缩计算方法.该方法采用基于周期对称的谐波平衡法,建立非线性的迭代方程,能够基于单扇区的矩阵求解整个失谐叶盘的强迫响应,且不存在任何截断误差,因而效率高于部件模态综合法但精度并不降低.从周期对称结构理论出发给出了方法原理的数学推导,在此基础上通过一个失谐叶盘有限元模型算例对该方法的效率与精度进行了讨论.结果表明:相较部件模态综合法,周期对称谐波平衡法矩阵维度减小59.5%,计算复杂度减少93.4%,计算时间节省57.4%.     

薄板弯曲问题的Trefftz方法

 本文从弹性力学广义变分原理出发,提出了一种求解薄板弯曲问题的Trefftz方法。为了让近似函数满足平衡微分方程,利用复变函数理论研究了试函数空间的构造问题。最后对矩形板给出了3个算例,计算结果表明该方法具有较好的收敛性和稳定性。     

柔性臂空间机器人的一类运动学方程及轨迹规划

用柔性机械臂连杆末端的弹性变形以及变形角度来表示空间机器人柔性臂的弹性运动变量,克服了用无穷维振动模态变量来表示弹性变形给系统运动学建模带来的困难;基于广义雅可比矩阵的思想,建立了柔性臂空间机器人"双广义雅可比矩阵"形式的运动学模型,该运动学模型描述了柔性臂弹性变形对空间机器人的运动影响;以运动学方程为基础,设计了柔性臂空间机器人的惯性空间内连续轨迹规划算法。仿真表明,规划的机械臂关节运动规律可以补偿柔性连杆振动给机械臂末端位置带来的影响,使机械臂末端位置准确沿着期望的轨迹运动。     

圆板混合边值问题的一种数值方法

 利用弹性薄板的广义变分原理和双调和方程的特解序列,提出了一种求解圆板混合边值问题的数值方法。最后通过两个实例计算说明该方法的收敛性和可靠性。     

线性时不变系统的完全跟踪控制器

从线性时不变系统的状态空间模型出发,借助于广义逆矩阵概念和相容线性方程组的最小范数解,从理论上推导了线性时不变系统的完全跟踪控制器方案。     

压电材料壳的广义变分原理和状态向量方程

将三维弹性材料的H-R变分原理引入到具有机-电耦合效应的三维压电弹性材料圆柱壳问题中,建立了对应于正交各向异性压电材料圆柱壳的Hamilton型广义变分原理。通过变分运算得到了柱坐标系下的状态向量方程,该状态向量方程也是Hamilton正则方程。最后指出纯正交各向异性弹性材料壳Hamilton型广义变分原理是变分公式的特例。变分原理的建立将有利于压电材料圆柱壳静力学和动力学问题的有限元方法或半解析法的推导。     

弱界面黏结通用单胞模型理论及应用

将通用单胞模型与弱界面分离模型结合.首先推导了增量型通用单胞模型计算公式, 然后采用界面位移公式对平均位移连续条件进行改造, 建立子胞平均应变增量和宏观应变增量的函数关系.应用子胞增量本构方程得到子胞平均应力增量表达式, 采用均匀化方法获得宏观应力增量表达式.最后采用渐进适应界面模型模拟了长纤维复合材料的横向拉伸过程.算例分析表明, 将界面分离模型和通用单胞模型结合能较好地模拟复合材料在拉伸载荷下的宏观响应.      

DOA estimation and mutual coupling calibration with the SAGE algorithm

In this paper, a novel algorithm is presented for direction of arrival（DOA） estimation and array self-calibration in the presence of unknown mutual coupling. In order to highlight the relationship between the array output and mutual coupling coefficients, we present a novel model of the array output with the unknown mutual coupling coefficients. Based on this model, we use the space alternating generalized expectation-maximization（SAGE） algorithm to jointly estimate the DOA parameters and the mutual coupling coefficients. Unlike many existing counterparts, our method requires neither calibration sources nor initial calibration information. At the same time,our proposed method inherits the characteristics of good convergence and high estimation precision of the SAGE algorithm. By numerical experiments we demonstrate that our proposed method outperforms the existing method for DOA estimation and mutual coupling calibration.