特征向量导数计算的扰动法

评述了有关特征向量导数的计算方法。综合了它们的优缺点后,提出了两种计算特征向量导数的扰动迭代法。方法的概念简单、实施容易。算例表明,迭代的收敛率较高,适于工程应用。     

用Chebyshev多项式加速的子空间迭代法

研究计算大型稀疏对称矩阵的若干个最大或最小特征值的问题，首先引入了求解大型对称特征值问题的子空间迭代法和Chebyshev迭代法，并对后者作了理论分析。为了加速子空间迭代法的收敛速度，作者用Chebyshev多项式来改进原始的子空间迭代法，即讨论Chebyshev迭代法对子空间迭代法的应用，从而给出了Chebyshev-子空间迭代法。最后把原始的方法和改进的方法计算数值例子的结果进行了比较，其结果表明Chebyshev-子空间迭代法比子空间迭代法优越，不仅收敛速度快，并且减少了计算量和计算时间。     

大型结构特征值问题的并行EBE－子空间迭代法

本文利用ＥＢＥ策略和预处理共轭梯度法（ＰＣＧ法），将广义特征值问题子空间迭代法中各步的计算都单元化，从而避免了总刚度和总质量矩阵的组集，大大节省了存储量。由此建立的ＥＢＥ－子空间迭代法尤其适宜于并行计算。在银河－２机上的数值算例结果表明，无论是串行，还是并行计算，该方法都能有效提高计算速度。如对模型问题，若网格取４８０，则在串行计算时，ＥＢＥ计算途径较传统的总体计算途径的速度提高倍数达３．２７，而在挂用４个处理机进行并行计算时的ＥＢＥ－子空间迭代法较串行的总体计算途径的速度提高倍数可达１１．４。总之，该方法为一种有效的大型结构动力分析问题的求解方法。     

适用于大型结构有限元分析的EBE／PCG法

本文在综述一般迭代法基础上,着重介绍采用EBE法作为预处理的共轭梯度法,即EEE/PCG法。其实质在于在结构分析中只需对元素刚阵个别进行分解直接求解,从而避免了通常的组装成结构刚阵并分解求解的有限元步骤,在大型结构分析时,可显著减少所耗容量和机时。文中列举若干例题,由结果可见此法是可行的。在需采用平行解法和内外存交换的大型有限元方程求解时,充分显示其优越性。因此说EBE/PCG法是一种颇有发展前景的有限元方程解法。     

非结构网格上Euler方程有限体积解法的改进

对二维非结构三角形网格上Euler方程有限体积解法的格点格式进行了一些改进，重点在于提高数值解的精度，细致处理人工粘性项的尺度因子以及对该项建立适当的边界条件；发展一种新的基于最长边剖分的三角形网格自适应加密方法。采取多步Runge—Kutta格式推进、当地时间步长、隐式残值光顺等措施加速迭代收敛。文末给出的数值结果非常接近于参考文献中结构网格上的结果，验证了所发展数值方法的精确性。     

特征值重分析的广义伽略金法退化解

特征值的重分析方法在许多技术领域,如结构的动态优化设计分析、利用虚材料结构的模态试验数据提取(辨识)真实材料结构的试验模态参数等工作中都占有重要地位。作者曾为虚材料结构的试验模态处理建立过一种基于动柔度的广义伽略金法和基向量组合法,这两种方法作为特征值重分析技术,其"快速性"还不够理想。为此,本文采用了Kirsch的假设,通过解的结构原理得到广义伽略金法的退化解,其"快速性"非常好。数值结果表明,该方法的精度也满意。     

近似技术在遗传算法和结构优化设计中的应用

在传统结构优化设计中．基于近似技术和灵敏度分析的结构重分析技术可以显降低结构有限元分析的次数。本提出了一套适于遗传算法应用的近似模型——分段近似模型。该分段近似模型可以极大地降低结构有限元分析的次数，并且可以提高遗传算法的收敛性，从而提高遗传算法的计算效率。本通过两个实例验证了该分段近似模型的有效性。     

动力学参数修改的耦合模态子空间摄动法

在结构动力学设计中，如果系统的质量矩阵和刚度矩阵有小修改，则矩阵摄动法是结构动态特性再分析的重要方法，当系统具有重频或近频耦合模态时，常规的矩阵摄动法失效。本文研究了多组耦合模态结构动力学参数修改的子空间摄动法，将系统振动模态分布为耦合模态与弧立模态，由数组相互耦合的模态张成耦合模态子空间，采用初始模态坐标变换后，在维数不大的耦合模态子空间内进行了特征分析，然后作第二次耦合模态子空间坐标交换，最扣     

周期结构和太阳电池翼的模态分析

基于周期结构的振动理论,把折叠状态太阳电池翼当作拟周期结构处理,解释了太阳电池翼在进行模态分析时所出现的“主模态局部化”现象,了解该现象对于在实验模态分析中确定阻尼系数具有重要意义。同时从理论上说明了,简单的基板模型的第1阶特征值是整翼模型的第1阶特征值的近似值。两个算例计算的结果支持了这个观点。     

结构的模糊随机可靠性分析

根据工程设计中存在的大量的模糊性因素，将模糊理论与结构可靠性理论相结合，分别建立了基本空间为连续集和离散集时结构的模糊随机可靠度公式及一些特殊可靠度公式；运用模糊集合基本运算法则，获得结构安全这一模糊随机事件的隶属度表达式；通过算例形象地阐述了计算方法，并分别给邮由随机模糊性，随机性和模糊性所引起的结构可靠度值。     

大型捆绑火箭振动试验液体取代与支承影响修正方法

在大型捆绑火箭振动中，液体取代和支承边界影响作为结构质量和刚度的变化，对火箭振动特性有重要影响。本文采用有限元再分析法，矩阵摄动法和拓广灵敏度法，对模态试验数据进行了合理修正，并获得了满意的修正结果。研究表明，在秒状态足够多的情况下，用最小二乘曲线拟合技术修正振型斜率也是可行的。     

求解非对称线性方程组的松弛混合算法

求解大型稀疏非对称线性方程组的混合迭代算法通常会由于系数矩阵的谱分布较广而导致收敛失败。本文通过在迭代多项式中加入变化的松驰因子定义了一类松驰混合算法。选择适当的松驰因子可以显著地改善算法的收敛效果。     

复合材料层压板热屈曲优化设计的自适应免疫遗传算法

复合材料层压壁板的热屈曲优化问题是高速飞行器结构设计的重点考虑内容。通过对免疫遗传算法引入自适应交叉和变异,构造了一种自适应免疫遗传算法(AIGA),并将该算法应用于考虑强度约束的层压板热屈曲铺层顺序优化设计。并将算法的优化结果与简单遗传算法(SGA)、免疫遗传算法(IGA)的优化结果进行了比较,结果表明该算法收敛速度快,优化解的质量最好,并有效的克服了SGA易于早熟收敛,IGA收敛缓慢的缺点。同时研究了抗体调节系数对AIGA算法性能的影响。     

基于非确定性的翼型鲁棒反设计方法

实际的翼型是在一个非确定性的状态范围内工作的,为此应用鲁棒控制理论建立了翼型的鲁棒气动反设计方法及数值求解过程,基本的优化器为基于伴随方程的位流优化方法。对守恒型位流方程应用高效AF 2格式求解,对非守恒型的伴随方程采用了旋转差分和高效AF 3格式,提高了梯度计算的精度和效率。优化算例表明:本文方法所得到的翼型在一个较宽的工作状态范围内都具有良好的性能。     

三角面元灰度图象三维重建算法研究

三角面元SFS算法是结合三角面元模型和线性反射图来阐述SFS问题的,该算法的基本思想是用三角面元素来近似光滑表面,并用一系列基于节点的线性方程组合来表达该表面。本文在此算法的基础上,求解时通过加权平均迭代,提高了恢复的准确性和迭代的收敛速度。     

非光滑动力系统周期响应的数值解法

工程中有重要意义的非线性动力系统周期响应分析,一般因系统维数高、特性复杂而需要用数值方法进行。但目前有影响的打靶法、增量谐波平衡法均不能处理含弹塑性、间隙、干摩擦等非光滑因素的工程系统,计算光滑动力系统的强非线性周期响应时也常不收敛。本文分析了上述两种方法的欠缺,指出其原因在于方法中直接或隐含的Newton迭代格式。文中提出了用拟Newton迭代格式和无约束优化方法改进打靶法和增量谐波平衡法的8种方案,并对这些方案就收敛性、精度和效率进行了考核比较,给出了不同方案的适用对象,精选出了几种方案推荐工程界应用。     

特征值重分析的递推算法和若干技术比较

对于利用虚材料结构的模态试验数据提取（辨识）真实材料结构的试验模态参数（简称虚材料结构试验模态）的问题,作者曾建立过一种基于动柔度的广义伽略金法（GG）和基向量组合法（DF）。为了特征值重分析的“快速性”要求,作者又得到广义伽略金法的一种退化解（DS）,其“快速性”非常好,但精度比广义伽略金法差。为了全面了解,还对若干相关方法,如GG、DF、DS和CA（combinedapproximation）等,在满足修改结构特征方程的近似性方面作了探讨,同时还就它们的数值结果进行了比较。其探讨表明,GG法的精度最好,而其余方法的精度是相同或相当的、且基本满意。这些方法在修改结构之特征方程Ku=2Mu的右边惯性项中都不同程度引入了Kirsch的近似假设,故而会引入不同程度的误差。现在要问：右边惯性载荷项的误差对方程解之精度的影响倒底有多大？为此,首次提出了一种递推算法。该递推算法的第一级递推解恰好是广义伽略金法的退化解（DS法）。递推算法的概念适用于上述各种方法。为此,作者又首次建立了CA法的递推公式。最后,该递推算法从数值上得出一个重要结论：导致右边惯性载荷项误差的“Kirsch近似假设”对特征值重分析结果的影响不严重。同时,该递推算法还能逐级减小由Kirsch的近似假设给重分析结果带来的误差。