We develop a density based topology optimization method for linear elasticity based on the cut finite element method. More precisely, the design domain is discretized using cut finite elements which allow complicated geometry to be represented on a structured fixed background mesh. The geometry of the design domain is allowed to cut through the background mesh in an arbitrary way and certain stabilization terms are added in the vicinity of the cut boundary, which guarantee stability of the method. Furthermore, in addition to standard Dirichlet and Neumann conditions we consider interface conditions enabling coupling of the design domain to parts of the structure for which the design is already given. These given parts of the structure, called the nondesign domain regions, typically represents parts of the geometry provided by the designer. The nondesign domain regions may be discretized independently from the design domains using for example parametric meshed finite elements or isogeometric analysis. The interface and Dirichlet conditions are based on Nitsche's method and are stable for the full range of density parameters. In particular we obtain a traction-free Neumann condition in the limit when the density tends to zero.

We design and analyze a hybridized cut finite element method for elliptic interface problems. In this method very general meshes can be coupled over internal unfitted interfaces, through a skeletal variable, using a Nitsche type approach. We discuss how optimal error estimates for the method are obtained using the tools of cut finite element methods and prove a condition number estimate for the Schur complement. Finally, we present illustrating numerical examples.

We present a cut finite element method for shape optimization in the case of linear elasticity. The elastic domain is defined by a level-set function, and the evolution of the domain is obtained by moving the level-set along a velocity field using a transport equation. The velocity field is the largest decreasing direction of the shape derivative that satisfies a certain regularity requirement and the computation of the shape derivative is based on a volume formulation. Using the cut finite element method no re-meshing is required when updating the domain and we may also use higher order finite element approximations. To obtain a stable method, stabilization terms are added in the vicinity of the cut elements at the boundary, which provides control of the variation of the solution in the vicinity of the boundary. We implement and illustrate the performance of the method in the two-dimensional case, considering both triangular and quadrilateral meshes as well as finite element spaces of different order.

In this paper, we propose a stabilized finite element method for the numerical solution of contact between a small deformation elastic membrane and a rigid obstacle. We limit ourselves to friction-free contact, but the formulation is readily extendable to more complex situations.

We propose a boundary value correction approach for cases when curved boundaries are approximated by straight lines (planes) and Lagrange multipliers are used to enforce Dirichlet boundary conditions. The approach allows for optimal order convergence for polynomial order up to 3. We show the relation to a Taylor series expansion approach previously used in the context of Nitsche's method and, in the case of inf-sup stable multiplier methods, prove a priori error estimates with explicit dependence on the meshsize and distance between the exact and approximate boundary.

We develop a finite element method for the Laplace–Beltrami operator on a surface with boundary and nonhomogeneous Dirichlet boundary conditions. The method is based on a triangulation of the surface and the boundary conditions are enforced weakly using Nitsche's method. We prove optimal order a priori error estimates for piecewise continuous polynomials of order k ≥ 1 in the energy and L^{2} norms that take the approximation of the surface and the boundary into account.

We develop a cut finite element method (CutFEM) for the convection problem in a so called fractured domain which is a union of manifolds of different dimensions such that a d dimensional component always resides on the boundary of a d+1 dimensional component. This type of domain can for instance be used to model porous media with embedded fractures that may intersect. The convection problem is formulated in a compact form suitable for analysis using natural abstract directional derivative and divergence operators. The cut finite element method is posed on a fixed background mesh that covers the domain and the manifolds are allowed to cut through a fixed background mesh in an arbitrary way. We consider a simple method based on continuous piecewise linear elements together with weak enforcement of the coupling conditions and stabilization. We prove a priori error estimates and present illustrating numerical examples.

8. Cut Finite Element Methods for Linear Elasticity Problems

Hansbo, Peter

et al.

Jönköping University, School of Engineering, JTH, Product Development.

Larson, Mats G.

Umeå University, Faculty of Science and Technology, Department of Mathematics and Mathematical Statistics.

Larsson, Karl

Umeå University, Faculty of Science and Technology, Department of Mathematics and Mathematical Statistics.

Cut Finite Element Methods for Linear Elasticity Problems2018In: Geometrically Unfitted Finite Element Methods and Applications: Proceedings of the UCL Workshop 2016 / [ed] Stéphane P. A. Bordas; Erik Burman; Mats G. Larson; Maxim A. Olshanskii, Springer, 2018, p. 25-63Chapter in book (Refereed)

Abstract [en]

We formulate a cut finite element method for linear elasticity based on higher order elements on a fixed background mesh. Key to the method is a stabilization term which provides control of the jumps in the derivatives of the finite element functions across faces in the vicinity of the boundary. We then develop the basic theoretical results including error estimates and estimates of the condition number of the mass and stiffness matrices. We apply the method to the standard displacement problem, the frequency response problem, and the eigenvalue problem. We present several numerical examples including studies of thin bending dominated structures relevant for engineering applications. Finally, we develop a cut finite element method for fibre reinforced materials where the fibres are modeled as a superposition of a truss and a Euler-Bernoulli beam. The beam model leads to a fourth order problem which we discretize using the restriction of the bulk finite element space to the fibre together with a continuous/discontinuous finite element formulation. Here the bulk material stabilizes the problem and it is not necessary to add additional stabilization terms.