Subdivision-based multilevel methods for large scale engineering simulation of thin shells

This paper presents a multilevel algorithm to accelerate the numerical solution of thin shell finite element problems de-scribed by subdivision surfaces. Subdivision surfaces have become a widely used geometric representation for general curved three dimensional boundary models and thin shells as they provide a compact and robust framework for mod-eling 3D geometry. More recently, the shape functions used in the subdivision surfaces framework have been proposed as candidates for use as finite element basis functions in the analysis and simulation of the mechanical deformation of thin shell structures. When coupled with standard solvers, however, such simulations do not scale well. Run time costs associated with high-resolution simulations (105 degrees of freedom or more) become prohibitive. The main contribution of the paper is to show that the subdivision framework can be used for accelerating such sim-ulations. Specifically the subdivision matrix is used as the intergrid information transfer operator in a multilevel pre-conditioner. The method described in the paper allows the practical simulation or a broad range of problems. Included examples show that the run time of the algorithm presented scales nearly linearly in time with problem size.

[1]  J. Hall,et al.  The multigrid method in solid mechanics: Part I—Algorithm description and behaviour , 1990 .

[2]  J. C. Simo,et al.  On stress resultant geometrically exact shell model. Part I: formulation and optimal parametrization , 1989 .

[3]  Tony DeRose,et al.  Piecewise smooth surface reconstruction , 1994, SIGGRAPH.

[4]  Leonidas J. Guibas,et al.  Primitives for the manipulation of general subdivisions and the computation of Voronoi diagrams , 1983, STOC.

[5]  Michael Ortiz,et al.  Fully C1‐conforming subdivision elements for finite deformation thin‐shell analysis , 2001, International Journal for Numerical Methods in Engineering.

[6]  James Demmel,et al.  A Supernodal Approach to Sparse Partial Pivoting , 1999, SIAM J. Matrix Anal. Appl..

[7]  S. Ashby,et al.  A parallel multigrid preconditioned conjugate gradient algorithm for groundwater flow simulations , 1996 .

[8]  Jos Stam,et al.  Exact evaluation of Catmull-Clark subdivision surfaces at arbitrary parameter values , 1998, SIGGRAPH.

[9]  Hong Qin,et al.  A novel FEM-based dynamic framework for subdivision surfaces , 1999, SMA '99.

[10]  D. Zorin Stationary Subdivision and Multiresolution Surface Representations , 1997 .

[11]  Charles T. Loop,et al.  Smooth Subdivision Surfaces Based on Triangles , 1987 .

[12]  K. Bathe Finite Element Procedures , 1995 .

[13]  Leonidas J. Guibas,et al.  Primitives for the manipulation of general subdivisions and the computation of Voronoi diagrams , 1983, STOC.

[14]  Hong Qin,et al.  Dynamic Catmull-Clark Subdivision Surfaces , 1998, IEEE Trans. Vis. Comput. Graph..

[15]  Juan C. Meza,et al.  A Multigrid Preconditioner for the Semiconductor Equations , 1996, SIAM J. Sci. Comput..

[16]  Tony DeRose,et al.  Subdivision surfaces in character animation , 1998, SIGGRAPH.

[17]  Ulrich Reif,et al.  Curvature integrability of subdivision surfaces , 2001, Adv. Comput. Math..

[18]  Peter Schröder,et al.  Integrated modeling, finite-element analysis, and engineering design for thin-shell structures using subdivision , 2002, Comput. Aided Des..

[19]  Gabriel Taubin Is this a quadrisected mesh? , 2001, SMA '01.

[20]  George Turkiyyah,et al.  MULTIGRID PRECONDITIONER FOR UNSTRUCTURED NONLINEAR 3D FE MODELS , 1999 .

[21]  M. Ortiz,et al.  Subdivision surfaces: a new paradigm for thin‐shell finite‐element analysis , 2000 .

[22]  Yousef Saad,et al.  Iterative methods for sparse linear systems , 2003 .