Finite element analysis & design optimization in a distributed computing environment

As product designs become more sophisticated, both the finite element (FE) and the design optimization (DO) models have grown bigger. At the same time there is increasing evidence that computing clusters created with commodity chips are capable of outperforming traditional supercomputers. In this paper, the HYI-3D design optimization software system is discussed. HYI-3D has been developed to work in sequential and distributed processing environments. The two major objectives are as follows (a) implement the FE methodology using the well-known domain decomposition technique where the original FE model is split into a number of smaller subdomains, and (b) implement a design optimization methodology for sizing, shape and topology optimization using coarse and fine-grain parallelism. For the FE engine, a direct sparse solver is used at the subdomain level and preconditioned conjugate gradient (PCG) is used at the interface system level. The finite element equations are then generated and assembled at the individual domain level in parallel. Matrix and vector operations involving sparse matrices form the bulk of the computations in this step. Once these equations are assembled, the condensed system level equations are formed. These condensed system level equations are usually much smaller (but denser) than the original system equations, and hence can be computationally expensive. With respect to design optimization, multi-level parallelism is employed. Not only can the finite element analysis be carried out in parallel but also other steps in the design optimization algorithm can be computed in parallel – gradients, line search and direction-finding problem. Numerical examples show the gains obtained from coarse and fine grain parallelism.

[1]  Manfred Grauer,et al.  Grid -based Computing for Multidisciplinary Analysis and Optimization , 2004 .

[2]  Gerhard Venter,et al.  EFFICIENT OPTIMIZATION ALGORITHMS FOR PARALLEL APPLICATIONS , 2000 .

[3]  Duc T. Nguyen,et al.  An algorithm for domain decomposition in finite element analysis , 1991 .

[4]  Willie R. Watson,et al.  Parallel Finite Element Domain Decomposition for Structural/Acoustic Analysis , 2005 .

[5]  A. Belegundu,et al.  A shape optimization approach based on natural design variables and shape functions , 1988 .

[6]  Ashok D. Belegundu,et al.  Parallel Line Search in Method of Feasible Directions , 2004 .

[7]  Ole Sigmund,et al.  A 99 line topology optimization code written in Matlab , 2001 .

[8]  Jack Dongarra,et al.  High performance computing: Clusters, constellations, MPPs, and future directions , 2003 .

[9]  Barry Hilary Valentine Topping,et al.  Parallel simulated annealing for structural optimization , 1999 .

[10]  Subramaniam Rajan,et al.  Toward a practical design optimization tool , 1996 .

[11]  M. Zhou,et al.  The COC algorithm, Part II: Topological, geometrical and generalized shape optimization , 1991 .

[12]  Ole Sigmund,et al.  On the Design of Compliant Mechanisms Using Topology Optimization , 1997 .

[13]  Bernard Grossman,et al.  A Coarse-Grained Parallel Variable-Complexity Multidisciplinary Optimization Paradigm , 1996, Int. J. High Perform. Comput. Appl..

[14]  M. Bendsøe Optimal shape design as a material distribution problem , 1989 .

[15]  M. Bendsøe,et al.  Material interpolation schemes in topology optimization , 1999 .

[16]  S. L. Padula,et al.  Parallel implementation of large-scale structural optimization , 1998 .

[17]  Duc T. Nguyen,et al.  Nonlinear constrained optimizer and parallel processing for golden block line search , 1999 .