This work presents a novel iterative approach for mesh partitioning optimization to promote the efficiency of parallel nonlinear dynamic finite element analysis with the direct substructure method, which involves static condensation of substructures' internal degrees of freedom. The proposed approach includes four major phasesinitial partitioning, substructure workload prediction, element weights tuning, and partitioning results adjustment. The final three phases are performed iteratively until the workloads among the substructures are balanced reasonably. A substructure workload predictor that considers the sparsity and ordering of the substructure matrix is used in the proposed approach. Several numerical experiments conducted herein reveal that the proposed iterative mesh partitioning optimization often results in a superior workload balance among substructures and reduces the total elapsed time of the corresponding parallel nonlinear dynamic finite element analysis.
IntroductionParallel direct substructure method, which involves static condensation of substructures' internal degrees of freedom, is one of the most popular domain decomposition methods (Smith et al., 1996) for parallel finite element analysis. However, the parallel efficiency of the method remains a topic of concern. For example, computational workload imbalance among substructures inhibits good efficiency. For an efficient parallel finite element analysis, it is essential that the finite element mesh be partitioned such that among processors (or substructures) computational workloads are well balanced and inter-process communication is minimized. Although many mesh partitioning algorithms have been proposed (Fox, 1988;Farhat and Simon, 1995;Pothen, 1997;Hsieh et al., 1997), most employ simple and general assumptions in their optimization process for load balancing without considering the characteristics of a parallel solution algorithm. For example, these algorithms generally assume that a balanced number of elements among substructures implies a balance of workloads among processors. Therefore, the mesh partitions produced by these algorithms normally do not result in satisfactory parallel efficiency for a specific parallel solution algorithm, such as the parallel direct substructure method (Vanderstraeten et al., 1996;Yang and Hsieh, 1997). Hsieh et al. (1999) proposed an iterative mesh partitioning approach to improve workload balance among processors in parallel substructure finite element computations. This approach attempts to balance the computational workloads among substructures by iteratively adjusting the mesh partitions through element weight tuning. The element weights within a substructure are tuned in each iterative step according to the substructure workloads predicted by a set of empirical rules, which are derived from numerical experiments on numerous finite element meshes.Based on the basic idea of iterative mesh partitioning approach of Hsieh et al. (1999), a new iterative mesh partitioning approach with a more accurate work...