This project is studying a new class of very promising, highly parallel iterative methods for the solution of large-scale systems which arise from high-order elements (the p-version finite element method) in three-dimensional elasticity for structures with plates, shells, and junctures. The size of problems of interest now rules out direct solvers, but high order elements, particularly thin elements, give rise to very ill-conditioned systems of equations, which necessitate new specific preconditioners.Researchers are investigating new methods for preconditioning thin, solid, high-order elements and studying their performance on three-dimensional problems with up to over one million degrees of freedom. The design and development is guided by theoretical analyses as well as extensive experimentation. Research will result in a coarse grained parallel prototype code running on CRAY-YMP. The project will further knowledge in the areas of multilevel and domain decomposition computational methods on parallel supercomputers. The availability of the methods will have a significant impact on mechanical engineering practice.Commercial Applications:This project will develop a fast and robust solver for large-scale system arising in the analysis of structures including shells and plates by the p-version finite element method on supercomputers. This software will be sold to developers of p-version Commercial packages. The method will be of interest in large-scale finite element computations where the quality of the solution is essential.
