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

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

CFD问题显式有限差分方法的并行化

基于ＭＩＭＤ并行计算机模型，本文讨论了ＣＦＤ问题应用显式有限差分方法的并行化问题。利用区域分解法将计算问题分解为多个子问题，每个子问题由不同的处理器分别处理。针对ＳＣＢ计算格式的特点，尽量降低各处理器间的数据通讯，提高了并行计算效率。最后给出二维Ｅｕｌｅｒ方程组计算实例，计算结果令人满意，且算法有较好的可扩放性。     

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

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

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

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

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

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

线性矩阵方程AXB－CXD＝E的误差界

首先讨论线性矩阵方程ＡＸＢ－ＣＸＤ＝Ｅ在有唯一解的条件下方程之解Ｘ的一个上界，再考虑系数矩阵Ａ，Ｂ，Ｃ，Ｄ，Ｅ发生小扰动时，线性矩阵方程（Ａ＋δＡ）Ｘ（Ｂ十δＢ）－（Ｃ＋δＣ）Ｘ（Ｄ＋δＤ）＝Ｅ＋δＥ之解Ｘ作为Ｘ的近似值的相对误差的一个上界。     

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

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

DFP法在消息传递型MIMD上的并行实现

研究了ＤＦＰ算法在消息传递型ＭＩＭＤ并行机上的并行实现问题，首先讨论并行一维寻优问题，然后利用矩阵和向量分解，并充分考虑到Ｔｒａｎｓｐｕｔｅｒ并行机的特点，给出了一个负载平衡度高，通讯最小的并行化方案。理论分析和实际算例表明，当问题的规模充分大时，其并行效率趋于１。     

卷帘行存储下的一种并行Cholesky分解及其在PAR95上的实现

Ｃｈｏｌｅｓｋｙ 分解在科学与工程计算中占有重要的地位，串行的 Ｃｈｏｌｅｓｋｙ 分解已有成熟的方法，但并行的 Ｃｈｏｌｅｓｋｙ 分解方法要充分考虑机器体系的结构，在拥有共享内存的 Ｍ Ｉ Ｍ Ｄ 型多处理机、 Ｓ Ｉ Ｍ Ｄ 型向量机系统上前人已有较好的工作，本文给出适用于 Ｍ Ｐ Ｐ 大规模并行计算机的卷帘行存储行格式的并行 Ｃｈｏｌｅｓｋｙ 分解算法，该算法使用了优先计算优先发送的策略，减少了结点机之间相互等待的时间，建立了结点机之间用于通讯的通讯数组，避免了使用撒播这一通讯模式，减少了通讯时间，通过在 Ｐ Ａ Ｒ９５ 上的数值试验表明，随着问题规模的扩大，并行效率越来越高，并且该法容易推广到多行卷帘存储的形式。     

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

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

A PARALLEL IMAGINARY EBE METHOD FOR SOLVING POSITIVE DEFINITE LINEAR SYSTEMS

ＡＰＡＲＡＬＬＥＬＩＭＡＧＩＮＡＲＹＥＢＥＭＥＴＨＯＤＦＯＲＳＯＬＶＩＮＧＰＯＳＩＴＩＶＥＤＥＦＩＮＩＴＥＬＩＮＥＡＲＳＹＳＴＥＭＳＰａｎＸｉａｏｓｕ（ＤｅｐａｒｔｍｅｎｔｏｆＡｐｐｌｉｅｄＭａｔｈｅｍａｔｉｃｓ，ＰｈｙｓｉｃｓａｎｄＭｅｃｈａｎｉｓ...     

三维非结构混合网格高超声速流场并行计算

在非结构混合网格上对三维高超声速化学非平衡粘性绕流进行了基于PC-Cluster的分布式并行数值模拟.本采用区域分裂思想,研究了三维非结构混合网格区域自动分解技术,并以此为基础对高超声速化学非平衡绕流进行了并行数值计算.控制方程为多组分N-S方程,空间离散采用有限体积格心格式,时间推进为显式Runge-Kutta格式.化学非平衡动力学模型为七组元带电离反应模型,对化学反应源项进行了点隐式处理,温度场的计算采用牛顿迭代法.在PC-Cluster上对三维双椭球模型的高超声速绕流流场进行了基于区域分解技术的并行数值模拟,所得数值结果与参考献中的结果作了对比验证。     

Increment-Dimensional Scaled Boundary Finite Element Method for Solving Transient Heat Conduction Problem

An increment-dimensional scaled boundary finite element method (ID-SBFEM) is developed to solve the transient temperature field.To improve the accuracy of SBFEM,the effect of high frequency factor on dynamic stiffness is considered,and the first-order continued fraction technique is used.After the derivation,the SBFE equations are obtained,and the dimensions of thermal conduction,the thermal capacity matrix and the vector of the right side term in the equations are doubled.An example is presented to illustrate the feasibility and good accuracy of the proposed method.     

Identification of Magnetic Bearing Stiffness and Damping Based on Hybrid Genetic Algorithm

Identifying the stiffness and damping of active magnetic bearings(AMBs)is necessary since those parameters can affect the stability and performance of the high-speed rotor AMBs system.A new identification method is proposed to identify the stiffness and damping coefficients of a rotor AMB system.This method combines the global optimization capability of the genetic algorithm(GA)and the local search ability of Nelder-Mead simplex method.The supporting parameters are obtained using the hybrid GA based on the experimental unbalance response calculated through the transfer matrix method.To verify the identified results,the experimental stiffness and damping coefficients are employed to simulate the unbalance responses for the rotor AMBs system using the finite element method.The close agreement between the simulation and experimental data indicates that the proposed identified algorithm can effectively identify the AMBs supporting parameters.     

挠性膜片联轴器优化设计

本文改进了Ｐｏｗｅｌｌ的直接优化方法，提出了有效的约束函数法，并用此法研究了膜片联轴器在强度刚度振动条件下，以重量轻，刚度最小为目标的优化设计计算，基于膜片的强度，刚度是在有限元分析基础上的，必须要很好地解决有限元前置处理（自动分单元、分载荷）和后置处理（绘等值线）。开发了一套比较完善的膜片外形优化设计软件系统，简称ＯＤＦＣ系统。     

三维有限元网格的自动划分

本文作者从研究结构的数学模型入手,讨论了建立三维结构有限元计算模型的方法。成功地获得了一些具有应用代表性的三维板壳和块体结构有限元计算网格的自动划分原理和算法。     

面向网络机群系统的结构有限元并行分析

提出了基于网络机群的有限元并行分析方法,结合网络机群体系的特点,分析了这一方法所应具有的基本特点:分布式并行建模;子区域间计算任务弱相关;节点负载平衡。在自建的网络并行平台上,实现了无重叠区域分解直接算法的并行处理,并进行了算例验证。