一类广义特征值问题的扰动研究

研究了实正则对的广义特征值扰动问题。对于求解结构动力学、弹性动力学中的振动模态及分析受外界扰动后引起的谐振频率的变化都有一定的意义。     

混合移频静刚度法和直接扰动迭代法

本文对混合移频动柔度法的失代式进行等价变换，从而得到一个“变种”，该变种不包含混合移频系统的动柔度式F^△（λ*），这便消除了F^△（λ*）带来的近似性，另外，对直接扰动法也进行等价变换，引出可以省时地用于许多特征向量导数计算的直接扰动迭代法。最后，为了供读者参数，本文还导出了三种混合移频动柔度法的非迭代式，其中有两种恰好是另外两种动柔度法经混合移频后的非迭代式，这实际上已将另两种动柔度法也扩展到密集根情况了。     

密集根之特征向量导数的改进高精度动柔度法

作者曾为许多特征向量导数计算，提出过一种介于直接法（指直接求解线代方程组的，以Nelson为代表的一类方法）和模态法之间的高精度动柔度法。这种方法与作者建立的一般动柔度法一样，不能用于密集根之特征向量导数的计算。为使该方法扩展到密集根状态，本文将作者发展的混合移频技术应用于原高精度动柔度法，并重新推导了混合移频系统 的高精度动柔度式，从而形成能适用于许多密集根之特征向量导数计算的改进高精度动柔度法。     

关于Davidson—Lanczos方法的收敛率

本文对文[1]提出的求解大型对称矩阵A的极端(几个最大或最小)特征值及相应特征向量的Davidson—Lanczos方法,用Rayleigh—Ritz逼近理论,研究了该方法的收敛率。证明了由该方法产生的规范正交向量{v_i}_i~m=1是Krylov子空间K_m≡Span(v_1,Av_1,…,A~(m-1)v_1)的一组基。设A的k个最大特征值为又,λ_1>λ_2>…>λ_k,相应的近似特征值为λ_i~(m)(i=1,…,k),得到 这里γ_i(γ_i>1),W_i和W_i~(m)是常数。     

求解大型矩阵特征值问题的并行块Davidson方法

针对拥有共享内存的并行计算环境和微机网络并行计算环境,给出了求解大型稀疏对称矩阵部分极端特征对的并行块Davidson方法。该方法将矩阵A按行块分配到各处理器上,各处理器利用矩阵A的行块和投影子空间的正交基所组成矩阵V的行块进行运算,减少了处理机之间的通讯次数,实现了算法的并行计算。在微机网络并行计算环境和拥有共享内存并行计算环境IBMP650上的数值试验表明,该算法非常有效。     

一种基于数值计算的轴流压气机旋转失速稳定性模型——之一:理论模型及算法分析

本文提出了一个基于直接数值计算结果基础上的旋转失速线性稳定性模型,描述了如何从基本的Navier－Stokes方程出发,把轴流压气机流动稳定性问题转化为一个特征值问题的基本理论框架。对压气机不稳定边界的判断通过计算和检验特征值的方式来实现。本文给出了有限差分法和基于Gauss－Lobatto节点的谱配置方法两种数值离散的方式,并分析了这两种离散方式的收敛特性。     

模态综合的直接变换法

基于固定界面模态综合的Craig-Bampton方法,提出了模态综合的直接变换方法,该法在形式上更简洁,适用于大规模动力学系统的模态分析和动力学模型减缩。对采用直接变换法模态综合过程的计算误差做了分析,给出了计算误差界,阐明了模态综合后总体结构的固有频率与子结构保留主模态个数之间的关系。针对子结构界面自由度太大的情况,提出了直接变换法的分部算法,不需组集系统总刚度、总质量矩阵,仅在子结构的基础上完成系统的减缩。计算实例表明该算法的有效性及计算的稳定性。     

民用飞机壁板蒙皮及长桁布置结构优化设计

机身壁板是飞机结构设计的重要承载组件,轻量化、高效率、共通性设计及优化是民机设计关注的重点。首先提出一种耦合ABAQUS的Buckle分析及ISIGHT优化的设计方法,利用自编子程序获取ABAQUS屈曲特征值,将特征值输入ISIGHT中计算临界屈曲载荷,同步更新变量参数及ABAQUS文件并提交计算,迭代分析直至优化流程结束。采用上述方法考虑轴向压缩载荷情况,以壁板整体重量最小为优化目标,疲劳应力值为约束条件,对单曲度金属机身壁板的蒙皮厚度,长桁数量及长桁截面厚度等几何参数进行优化。在满足壁板结构承载能力及总重量最小条件下,综合考虑结构载重比,临界应力及壁板加筋比,对比分析出一组最优参数,并与工程算法结果对比吻合程度较好,两者相对误差为3.73%。该优化思路实现FEA平台与优化工作一体化,可用于复合材料壁板设计及结构件减重优化工作,一定程度上可缩短零组件设计周期。     

扩展酉矩阵束算法实现稀疏可重构天线阵的优化设计

针对方向图可重构天线阵列的稀疏布阵问题,提出一种基于扩展酉矩阵束算法的优化设计算法。首先建立以阵元位置和多组激励为变量的多方向图联合稀疏优化模型,并利用期望方向图采样数据构建Hankel块矩阵。然后通过centro-Hermitian化处理和酉变换将采样矩阵从复数域转换到实数域,舍弃实矩阵中较小的奇异值对可重构线阵进行稀疏。最后通过对等价矩阵束的广义特征值分解估计稀布阵元位置,进而得到每个方向图对应的激励。仿真验证了该方法能够以阵元非均匀稀疏分布的阵列形式有效实现多个方向图的精确重构。     

求解广义特征值反问题的数值方法

讨论一类广义特征值反问题的数值解法,这类问题包括加法、乘法和经典特征值反问题作为其特殊情况。基于行列式和最小奇异值的计算,文中给出了求解这类问题的两个二次收敛的数值方法,描述了在出现重特征值的情况下如何改进其中的一个方法以保持二次收敛性,并且给出了两个数值例子以解释收敛性结果