ADAPTIVE DYNAMIC RELAXATION ALGORITHM FOR NONLINEAR HYPERELASTIC STRUCTURES .2. SINGLE-PROCESSOR IMPLEMENTATION

An efficient single-processor implementation of the Adaptive Dynamic R elaxation (ADR) algorithm is developed. It is designed to exploit data locality and can exploit vectorization of the finite element computat ions. The ADR algorithm is used to solve for the non-linear static res ponse of two- and three-dimensional hyperelastic systems involving fri ctionless contact. Performance is compared with an existing finite ele ment code which utilizes a direct solution method. ADR is found to be reliable and highly vectorizable, and it outperforms the direct soluti on method for the highly non-linear problems considered. In addition, it permits the use of a very simple and efficient contact algorithm. I n contrast to direct solution methods, ADR has minimal memory requirem ents and is easily parallelizable and scalable to more processors. For the class of problems addressed, it represents a very promising appro ach for parallel-vector processing.