Umeå universitets logga

Thin and slender structures are widely occurring both in nature and in human creations. Clever geometries of thin structures can produce strong constructions while requiring a minimal amount of material. Computer modeling and analysis of thin and slender structures have their own set of problems, stemming from assumptions made when deriving the governing equations. This thesis deals with the derivation of numerical methods suitable for approximating solutions to problems on thin geometries. It consists of an introduction and four papers.

In the first paper we introduce a thread model for use in interactive simulation. Based on a three-dimensional beam model, a corotational approach is used for interactive simulation speeds in combination with adaptive mesh resolution to maintain accuracy.

In the second paper we present a family of continuous piecewise linear finite elements for thin plate problems. Patchwise reconstruction of a discontinuous piecewise quadratic deflection field allows us touse a discontinuous Galerkin method for the plate problem. Assuming a criterion on the reconstructions is fulfilled we prove a priori error estimates in energy norm and L²-norm and provide numerical results to support our findings.

The third paper deals with the biharmonic equation on a surface embedded in R³. We extend theory and formalism, developed for the approximation of solutions to the Laplace-Beltrami problem on an implicitly defined surface, to also cover the biharmonic problem. A priori error estimates for a continuous/discontinuous Galerkin method is proven in energy norm and L²-norm, and we support the theoretical results by numerical convergence studies for problems on a sphere and on a torus.

In the fourth paper we consider finite element modeling of curved beams in R³. We let the geometry of the beam be implicitly defined by a vector distance function. Starting from the three-dimensional equations of linear elasticity, we derive a weak formulation for a linear curved beam expressed in global coordinates. Numerical results from a finite element implementation based on these equations are compared with classical results.