基于PVM的网络并行有限元面向对象的设计

子结构是有限元并行计算常用的一种方法,本文采用面向对象的方法,首先对子结构进行了面向对象的设计,得到了其类层次结构图;然后针对工作站网络有限元并行计算环境。提出了基于PVM消息传递平台上的Shadow—Mirror数据传输模型,该模型在有限元并行计算数据传输时,充分发挥数据面向对象的特性,采用设置数据缓冲区、短消息合并等方法以缩短数据通信时间,并据此编制了相应的程序。计算结果表明,使用文中提出的面向对象的Shadow—Mirror数据传输模型可以得到较为理想的并行加速比,而且随着问题规模增大,并行加速比增高。本文研究内容为进一步开展基于工作站网络的并行有限元研究提供了一个可参考的基础。     

用多色染色理论实现有限元刚度矩阵并行组集

根据多色染色理论，在同一节点有任意多个单元邻接的情况下，对有限元的单元进行了分类；在刚度矩阵组集时，同类单元可以并行计算，从而提高了组集效率．该并行算法在PVM(并行虚拟机器)并行平台上进行了具体实现，取得了较好的并行效率．     

????????????????????????м???

面板堆石坝的三维非线性有限元分析,对并行计算提出了较高的要求. 在建模阶段挖掘出模 型内在的并行性,由相邻坝段构成子结构,既方便灵活又能满足区域分解法的要求. 编制了 适用于机群的并行多波前法计算程序,并应用于西龙池下库沥青混凝土面板堆石坝的有限元 数值分析,其计算结果和串行完全相同,通过对加速比和并行效率进行分析,结果表明采用 并行多波前法求解面板堆石坝的三维非线性有限元问题有效且算法简单,易于实现,无需辅 助的任务分配程序.     

大型结构特征值问题的混合粒度并行算法

本文提出一种求解大形结构特征值问题的粗细粒度混合并行算法：在子结构模态综合粗粒度并行算基础上,综合系统的特性值问题采用细粒度并行方式求解。细粒度并行包括子空间迭代法的子结构并行算法、雅可比分块并行计算的方法和一种Ｎｅｗｔｏｎ－Ｒａｐｈｏｎ迭代法在多处理器上任力均衡分配的有效策略。子空间迭代法的子结构并行计算的实施是利用子结构的刚度阵和质量阵而不必完全组集系统刚度阵和国求综合系统的特征值问题。利用雅     

无单元伽辽金法的并行计算

对无单元伽辽金法的并行计算进行了详细研究,并将其应用于弹性动力学问题。使用并行桶搜索算法进行节点搜索,使用并行几何搜索算法进行样点搜索,讨论了移动最小二乘MLS(Moving Least Squares)形函数及其导数的并行计算和方程组的并行求解,并利用多层图形划分实现负载平衡。最后给出了并行无单元伽辽金法应用于弹性动力学的计算流程和实例。计算结果表明无单元伽辽金法具有很高的并行性和很好的并行效率,对其进行并行计算具有非常重要的意义。     

结构的多机并行分析I——结构的并行有限元方法

本文详细分析和讨论了结构分析的并行有限元方法—并行预处理共轭梯度法(以下简称PPCG法)。着重讨论了基于自带存储器的多处理机系统的并行预处理算法问题,并由此提出了两种PPCG法:PPCG1和PPCG2法。这两种方法适用于以单道剖分(one-way dissection)的子结构法为基础的并行分析。由于这种剖分法产生的结构刚度矩阵具有箭头形状,可独立地消除各子结构的内部自由度,并且不会在刚度矩阵中产生新的非零元素,因此很适合具有较多处理机的并行机系统对复杂结构进行的并行分析。     

裂纹扩展独立于网格的一种有限元单元子划分方法

提出了一种有限元模拟裂纹扩展的单元子划分结合子结构的方法。本方法中,裂纹可以进入或穿过一个单元,或沿单元的边界扩展,因此裂纹可以沿任意路径扩展而不受初始网格的限制。对上述几类包含裂纹的单元按照裂纹的路径进行子划分,覆盖一条裂纹的所有子划分单元就组成了一个子结构,子结构规模随裂纹的扩展而增大。子结构中因单元子划分而新增的结点自由度,通过自由度的凝聚用初始网格结点的自由度表示,因此结构整体分析的总自由度不变。以上述方法为基础建立了裂纹萌生和扩展的准则。用本文的方法分析了单（双）材料无限大平面中心（界面）裂纹的裂尖场,验证了本文方法的精度,并模拟了颗粒复合材料中微裂纹在颗粒、基体和界面中逐步扩展的过程,考核了本文方法对复杂裂纹扩展问题模拟的适用性。     

并行元胞单元法

在对元胞单元法的基本原理与特点进行简要介绍之后,探讨了元胞单元法的并行计算方法,给出了算法描述和解题步骤。以平面应力问题分析为例,进行了数值试验,给出了详细的计算步骤和结果。结果表明元胞单元法可形成一种高度并行的算法,适用于未来并行机的要求,有望在大型结构和纳米力学的大规模摸拟方面得到应用。     

多波前并行处理的弹塑性子结构并行有限元

研究一种基于ＰＶＭ的弹塑性子结构并行有限元法。利用多个波前,对各子结构并行地进行静凝聚。再采用预条件共轭梯度法（ＰＣＧ）并行求解界面方程。算例表明该方法能获得较好的并行加速比,同时也能有癣地节省内存量。     

结构模态分析的有限元并行模拟

以子结构模态综合分析为基础,提出一种求解大型结构特征值问题的并行解法．采用子结构模态综合算法,结构特征模态采用子空间迭代方式并行求解．这种子空间迭代法的子结构并行计算的实施是利用子结构的刚度阵和质量阵而不必完全组集系统刚度阵和质量阵求解综合系统的特征值问题．数值结果表明这种求解大型结构特征值问题的并行算法是可行有效的．     

大型结构特征值问题的Lanczos子结构并行算法

本文利用子结构和Lanczos方法,提出了大型结构固有频率与模态的并行解法。该方法在Lanczos方法的求解过程中,仅利用子结构刚度阵和质量阵并行进行凝聚,进而求得新的迭代矢量,最终求得三对角阵对应的特征值和特征向量。该算法在西安交通大学ELXSI-6400并行计算机上程序实现,计算结果表明能有效地节省计算时间和计算机的内存,为一种有效的大型工程结构动力问题的求解方法。     

Modelling of the micromechanical behaviour of unidirectional fibre-reinforced ceramics by the example of SiC/SiC

Summary The modelling of the mechanical properties of a substructure, based on the properties of the individual phases and their interactions, is presented on the example of the unidirectional SiC fibre-reinforced SiC composite. The substructure is selected in such a way that the total structure can be modelled from a wide number of substructures. The numerical evaluation of the model is accomplished by means of the finite element method (FEM). Finally, in numerical simulations of a particular example, the statistically verified events of damage in the substructure are described.     

有限元子结构并行算法的性能分析

 对有限元子结构并行算法的主要步骤进行了归纳，比较准确地 计算了算法的计算量和通信量. 并在此基础上分析了算法的性能和可 扩展性，最后指出了影响算法并行效率的关键因素，并给出了算法优 化方向和改进建议.     

求解弱不连续问题的p型自适应有限元方法

在实际工程计算中,存在大量的弱不连续问题,如含夹杂问题。利用通常的有限元方法,为确保界面上各点满足给定高精度,往往需要采用全域网格加密或全域提高单元阶次的方法,这将会导致计算机的物理内存和CPU时间的剧烈增长。p-型自适应有限元方法是一种能通过自适应分析逐步增加单元阶次以改善计算精度的数值方法。本文,我们针对弱不连续问题设计了相应的p-型自适应有限元方法,重点讨论了容许误差控制标准对界面上各点计算结果的影响,并对几类典型的弱不连续问题进行了数值计算与模拟。数值结果表明,本文设计的p-型自适应有限元方法对求解弱不连续问题是非常有效的,用较少的单元得到精度可靠的数值结果,可大大提高其有限元分析效率。     

Cyclo-symmetric modal synthesis of rotating bladed disc system

The natural mode analysis of bladed disc system at an arbitrary constant rotational speed is investigated in this paper and the large deformation effect caused by the pretorsion of the blades is considered. The computational method of group theory is employed in substructure technique and equilibrium position and stress distribution of the bladed disc system in centrifugal field are obtained through simultaneous iteration, furthermore, the initial stress stiffness matrix, the large deformation stiffness matrix and the corresponding rotating natural frequencies and modal shapes are determined. Because the large deformation effect of the blades is considered and all the interface degrees of freedom between the substructures can be eliminated in the present method, the computation work can be cut down remarkably under the precondition that the computational accuracy is ensured. This shows that, in general, the substructure simultaneous iteration method presented in this paper has a good convergency, thus greatly reducing the computation work.     

基于界面主从位移控制的频响子结构方法

在复杂工程结构频响子结构分析中,各个子结构的频响函数往往是由不同的部门获得,各个子结构的频响函数可能是在不同的坐标系下得到,亦或有些子结构具有某种对称性,这需要对相应的子结构频响函数进行旋转和镜像变换后才能进行整体拼接。此外,子结构拼接时交界面的连接点有时可能存在不一致,为了满足位移协调和力平衡条件,本文基于界面主从位移控制原理,提出一种改进的频响子结构方法,该方法既能考虑子结构的旋转和镜像变换,也能考虑子结构界面间的柔性连接,并可大幅提高计算效率。当界面位移控制阵为单位阵时,改进的方法退化为传统的频响子结构方法。数值算例验证了所提方法的正确性和有效性。     

有限元子结构并行算法的效能分析

引入了有限元子结构的并行算法,该方法通过划分子结构,由多个CPU并行进行静凝聚,然后求解界面位移,最后返回各小结构并行求解内点位移和应力.该方法在ELXSI 6 400并行机上程序实现,表明能大幅度节省机时.本文同时也对该方法的效能进行了分析.