Improving the robustness of the finite cell method for problems with large deformations

Produktbeschreibung

Zahlreiche Ingenieuranwendungen mit komplexen Strukturen und mikrostrukturierten Materialien wie Schaumstoffen oder 3D-gedruckten Objekten stellen eine Herausforderung für die traditionelle Finite-Elemente-Methode (FEM) dar. Die Finite-Cell-Methode (FCM) kombiniert den fiktiven Gebietsansatz mit hierarchischen Ansatzfunktionen, um den Netzgenerierungsaufwand zu reduzieren. Eine nicht-negative Moment-Fitting-Quadratur-Regel wird entwickelt, um die Integration von gebrochenen Zellen zu verbessern und die NewtonRaphson-Methode in der nichtlinearen Berechnung zu stabilisieren. Zusätzlich wird eine Eigenwert-Stabilisierungs-Technik entwickelt, um die Konditionszahl der Schnittzellen zu reduzieren. Ein Remeshing-Verfahren wird entwickelt, um die FCM-Robustheit zu erhöhen und größere Lastschritte zu ermöglichen. Numerische Beispiele veranschaulichen die Wirksamkeit dieser Ansätze.

Engineering applications often involve complex structures and microstructured materials, such as foams or 3D printed objects, posing challenges for traditional finite element methods due to the complex geometry. The finite cell method (FCM) is employed to eliminate mesh creation overhead, combining the fictitious domain approach with high-order basis functions. Special quadrature schemes are implemented to integrate broken cells effectively, addressing issues with negative weights through a non-negative moment fitting quadrature rule. An eigenvalue stabilization technique is developed to avoid ill-conditioning caused by broken cells, while a remeshing technique is applied to enhance the FCM robustness for large deformation analysis. Different numerical examples showcase the effectiveness of these approaches.

Contents
Abstract VIII
1 Introduction 1
1.1 Motivation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1
1.2 Aims and scope of this work . . . . . . . . . . . . . . . . . . . . . . . . . . 3
1.3 Outline of the thesis . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2 Basic equations of continuum mechanics 6
2.1 Kinematics . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2.1.1 Motion and deformation gradient . . . . . . . . . . . . . . . . . . . 6
2.1.2 Strain measures . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
2.2 Equilibrium and stress measures . . . . . . . . . . . . . . . . . . . . . . . . 8
2.2.1 Equilibrium . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
2.2.2 Stress measures . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
2.3 Constitutive equations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.3.1 Linear elasticity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 10
2.3.2 Hyperelasticity . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.3.3 Finite J2 elastoplasticity . . . . . . . . . . . . . . . . . . . . . . . . 11
2.4 Weak form of equilibrium . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.5 Weak form of equilibrium in the initial configuration . . . . . . . . . . . . 14
2.6 Weak form of equilibrium in the current configuration . . . . . . . . . . . . 14
2.7 Linearization of the weak form . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.7.1 Linearization of the weak form in the initial configuration . . . . . . 15
2.7.2 Linearization of the weak form in the current configuration . . . . . 16
3 The finite cell method 17
3.1 The weak forms in the extended domain . . . . . . . . . . . . . . . . . . . 17
3.2 Spatial discretization . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
3.2.1 Displacement field and mapping functions . . . . . . . . . . . . . . 19
3.2.2 Discretization of the weak forms . . . . . . . . . . . . . . . . . . . . 21
3.2.3 Boundary conditions . . . . . . . . . . . . . . . . . . . . . . . . . . 23
3.3 Numerical integration . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24
3.3.1 Gaussian quadrature . . . . . . . . . . . . . . . . . . . . . . . . . . 24
3.3.2 Adaptive Gaussian quadratures . . . . . . . . . . . . . . . . . . . . 25
4 Non-negative moment fitting quadrature for cut finite elements and cells 27
4.1 Moment fitting . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
4.2 Adaptive moment fitting . . . . . . . . . . . . . . . . . . . . . . . . . . . . 31
4.3 Reduced adaptive moment fitting . . . . . . . . . . . . . . . . . . . . . . . 31
4.4 Non-negative moment fitting . . . . . . . . . . . . . . . . . . . . . . . . . . 33
4.5 Non-negative least square solver . . . . . . . . . . . . . . . . . . . . . . . . 35
4.6 Preliminary investigations . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
4.6.1 Recovery of the Gauss-Legendre quadrature . . . . . . . . . . . . . 37
4.6.2 Cell cut by a sphere with radius r = 1.55 mm . . . . . . . . . . . . 37
4.6.3 Cell cut by a sphere with radius r = 0.3 mm . . . . . . . . . . . . . 41
4.7 Stabilization of the solution in the fictitious domain . . . . . . . . . . . . . 43
4.8 Application to problems in hyperelasticity . . . . . . . . . . . . . . . . . . 44
4.8.1 Plate with a cylindrical hole . . . . . . . . . . . . . . . . . . . . . . 44
4.8.2 Single pore of a foam . . . . . . . . . . . . . . . . . . . . . . . . . . 49
4.9 Application to problems in elastoplasticity . . . . . . . . . . . . . . . . . . 52
4.9.1 Plate with a cylindrical hole . . . . . . . . . . . . . . . . . . . . . . 53
4.9.2 Single pore of a foam . . . . . . . . . . . . . . . . . . . . . . . . . . 55
5 Eigenvalue stabilization technique to increase the robustness of the FCM 59
5.1 Material stabilization (α-method) . . . . . . . . . . . . . . . . . . . . . . . 61
5.2 Eigenvalue stabilization technique . . . . . . . . . . . . . . . . . . . . . . . 62
5.2.1 Excluding the rigid body modes from the stabilization . . . . . . . 66
5.2.2 Compute the rigid body modes analytically . . . . . . . . . . . . . 67
5.3 Preliminary investigations . . . . . . . . . . . . . . . . . . . . . . . . . . . 69
5.3.1 Plate with a circular inclusion . . . . . . . . . . . . . . . . . . . . . 69
5.3.2 Plate with a cylindrical hole . . . . . . . . . . . . . . . . . . . . . . 73
5.4 Nonlinear applications in finite strains . . . . . . . . . . . . . . . . . . . . 78
5.4.1 Combining the eigenvalue stabilization with the α-method . . . . . 78
5.4.2 Single cube connector . . . . . . . . . . . . . . . . . . . . . . . . . . 78
5.4.3 Complex cube connector . . . . . . . . . . . . . . . . . . . . . . . . 83
5.4.4 Single pore of a foam . . . . . . . . . . . . . . . . . . . . . . . . . . 86
6 Remeshing strategy for large deformations in the FCM 89
6.1 Kinematic relations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
6.2 A procedure for performing the remeshing technique . . . . . . . . . . . . . 93
6.3 Mesh generation and geometry description on any configuration state . . . 96
6.4 Remeshing criteria . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
6.4.1 Ratio of Jacobians . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
6.4.2 Orthogonality . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
6.4.3 Aspect ratio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
6.4.4 Investigation of the remeshing criteria on a single element . . . . . 98
6.5 Data transfer . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 100
6.5.1 Radial basis function interpolation . . . . . . . . . . . . . . . . . . 101
6.5.2 Inverse distance weighting interpolation . . . . . . . . . . . . . . . . 103
6.6 Numerical investigations in finite strain hyperelasticity . . . . . . . . . . . 104
6.6.1 Plate with a cylindrical hole . . . . . . . . . . . . . . . . . . . . . . 104
6.6.2 Single cube connector . . . . . . . . . . . . . . . . . . . . . . . . . . 107
6.6.3 Complex cube connector . . . . . . . . . . . . . . . . . . . . . . . . 110
6.6.4 Foam . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 113
6.7 Numerical investigations in finite strain elastoplasticity . . . . . . . . . . . 116
6.7.1 Plate with a cylindrical hole . . . . . . . . . . . . . . . . . . . . . . 116
6.7.2 Single cube connector . . . . . . . . . . . . . . . . . . . . . . . . . . 120
6.7.3 Complex cube connector . . . . . . . . . . . . . . . . . . . . . . . . 124
6.7.4 Single pore of a foam . . . . . . . . . . . . . . . . . . . . . . . . . . 127
7 Summary and outlook 130
Publications resulting from this project 134
Bibliography 135

Keywords: Finite-Cell Methode, Moment-Fitting-Quadratur, Eigenwert-Stabilisierung, Remeshing, Große Verformung, Finite cell method, Moment fitting quadrature, Eigenvalue stabilization, Remeshing, Large deformation