一种高度并行的Schwarz型混乱松弛法及其收敛性证明

文[1—2]把混乱松弛思想引入到Schwarz交替法中,构造了一种Schwarz型混乱松弛法。但这个方法在进行第n+1步迭代时,在拟边界上必须要用到第n步迭代的值,从而影响了算法的并行性,得不到相应的同步或异步MIMD并行算法、为此,本文给出一种高度并行的Schwarz型混乱松弛法,这个方法包括了Schwarz交替法及其相应的同步和异步MIMD并行算法。对于二阶线性与非线性微分方程Dirichlet问题,本文应用微分方程极值原理证明了该方法的收敛性。     

组合QR算法

本文提出了一个求解非对称实矩阵特征值问题的新方法—组合QR算法,给出了一些数值计算实例,并对此方法作了讨论。     

非对称广义特征值问题的拟-Eberlein算法及其并行化

非对称广义特征值问题的并行计算,目前在国内外研究得很少, Ｇ． Ｗ ． Ｓｔｅｗ ａｒｔ 和 Ｐ． Ｊ． Ｅｂｅｒｌｅｉｎ 曾分别研究非 Ｈｅｒｍ ｉｔｅ 矩阵标准特征值的并行拟 Ｊａｃｏｂｉ算法,１９８９ 年 Ｊ． Ｐ． Ｃｈａｒｌｉｅｒ 和 Ｐ． Ｖａｎ Ｄｏｏｒｅｎ 在 Ｇ． Ｗ ． Ｓｔｅｗ ａｒｔ 的工作基础上提出了求解非对称广义特征值问题的拟 Ｊａｃｏｂｉ算法（简称 Ｃ Ｖ 算法）与并行拟 Ｊａｃｏｂｉ算法。文中以 Ｊ． Ｐ． Ｃｈａｒｌｉｅｒ 等人的工作为基础,提出求解大型非对称广义特征值问题的拟 Ｅｂｅｒｌｅｉｎ 算法与并行拟 Ｅｂｅｒｌｅｉｎ 算法, Ｃｈａｌｌｅｎｇｅ Ｌ 并行系统上的数值试验表明,不仅并行效率很高,且敛速远优于 Ｃ Ｖ 算法     

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

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

有限元结构分析并行计算的若干研究进展

有限元结构分析并行处理在大规模科学与工程计算中占有重要地位，本文简要介绍作者在这一领域内的研究工作及其成果。涉及单元刚度矩阵的并行计算、有限元方程组的并行直接解法、有限元方程组的并行迭代解法、结构动力分析并行直接积分法、广义特征值问题的并行算法以及EBE（Element-By-Element）技术在有限元结构分析并行处理中的应用等研究领域。特别提出了“伪单元向量”与“单元分组技术”等概念与新方法；首先考虑了矩阵向量积的多自由度问题；系统、深入而又全面地研究了EBE策略在结构分析并行处理中的应用，填补了国内空白并突破了国际现有成果。     

求解广义特征值问题的并行保域多分法

近年来，随着并行机的发展，提出了代数特征值问题的并行多分法，但国内外的研究工作迄今仅限于对称三对角矩阵的标准特征值问题。在科学与工程众多领域内有着重要应用的广义特征值问题的多分法，因难度大等方面原因尚无人研究。本文提出广义特征值问题的并行保域多分法，该算法适用于大型稀疏实对称矩阵广义特征值问题的求解，它克服了传统的广义特征值问题的对分法（行列式查找法）出现的漏根或迭代不收敛等缺点，并保持其优点。作者在ＹＨ－１向量机上对这一算法进行了数值实验，并与并行保域行列式查找法作了比较。数值结果表明，该算法具有较高的加速比，当系统自由度为２１１４、求解特征对个数为３时，加速比可达７．７；且当问题规模较大时，并行保域多分法优于并行保域行列式查找法。     

瞬态热传导问题的隐式EBE并行算法

结合有限元法,研究了ＥＢＥ策略在热传导问题并行计算方面的应用,给出了瞬态热传导问题的隐式ＥＢＥ并行算法。根据所用ＣＰＵ个数并行计算单元刚度矩阵,而总刚度矩阵Ｋ不需要组装,仅在“单元级”上进行各种计算。就模型实验规模１０２４阶而言,机器存存储量比传统算法降低了４０％左右,在拥有共享内存的并行计算机Ｃｈａｌｌｅｎｇｅ上的数值实验表明,问题的规模越大,并行效率越高,所给算法非常有效。     

不规则结构有限元分析的并行处理及其在YH-1上的实现

在结构分析中,有限元方法是一个重要的数值方法。近些年来,关于单元刚度矩阵计算与总刚度矩阵合成的并行处理问题,A.K.Noor,梁维泰等人已进行了富有成效的工作,但他们的工作都是面向几何形状规则的结构分析问题。本文结合YH—1机的特点,对不规则结构分析问题,在提出单元分组技术的基础上,给出了一个单元刚度矩阵计算的并行算法ESVC,和变带宽存储格式下一个总刚度矩阵合成的并行算法ESVS。通过在YH—1并行机上对实例的计算结果表明:当同时计算的单元数r取120时,加速比s可达9.5,且r愈大,s就愈高。     

改进的行列式查找法及其在阻尼结构振动分析方面的应用

本文改进了K.K.Gupta求解二次特征值问题的行列式查找法及其程序,得出了改进的行列式查找法,指出了原程序出现的一些错误,并研究了该方法在结构系统振动分析方面的应用。根据本文的方法已编制成FORTRAN IV语言程序IDAMP,该程序可以计算阻尼振动及无阻尼振动问题,阻尼又包括粘性阻尼和结构阻尼;还可以计算旋转结构振动问题,旋转结构又包括有阻尼和无阻尼。该程序不仅可以计算全部特征值及相应的特征向量,而且能计算所指定范围内的特征位及相应的特征向量。 本文应用改进的行列式查找法成功地处理了有阻尼扭振系统的固有频率问题。