TY - THES A1 - Häfner, Stefan T1 - Grid-based procedures for the mechanical analysis of heterogeneous solids N2 - The importance of modern simulation methods in the mechanical analysis of heterogeneous solids is presented in detail. Thereby the problem is noted that even for small bodies the required high-resolution analysis reaches the limits of today's computational power, in terms of memory demand as well as acceptable computational effort. A further problem is that frequently the accuracy of geometrical modelling of heterogeneous bodies is inadequate. The present work introduces a systematic combination and adaption of grid-based methods for achieving an essentially higher resolution in the numerical analysis of heterogeneous solids. Grid-based methods are as well primely suited for developing efficient and numerically stable algorithms for flexible geometrical modeling. A key aspect is the uniform data management for a grid, which can be utilized to reduce the effort and complexity of almost all concerned methods. A new finite element program, called Mulgrido, was just developed to realize this concept consistently and to test the proposed methods. Several disadvantages which generally result from grid discretizations are selectively corrected by modified methods. The present work is structured into a geometrical model, a mechanical model and a numerical model. The geometrical model includes digital image-based modeling and in particular several methods for the theory-based generation of inclusion-matrix models. Essential contributions refer to variable shape, size distribution, separation checks and placement procedures of inclusions. The mechanical model prepares the fundamentals of continuum mechanics, homogenization and damage modeling for the following numerical methods. The first topic of the numerical model introduces to a special version of B-spline finite elements. These finite elements are entirely variable in the order k of B-splines. For homogeneous bodies this means that the approximation quality can arbitrarily be scaled. In addition, the multiphase finite element concept in combination with transition zones along material interfaces yields a valuable solution for heterogeneous bodies. As the formulation is element-based, the storage of a global stiffness matrix is superseded such that the memory demand can essentially be reduced. This is possible in combination with iterative solver methods which represent the second topic of the numerical model. Here, the focus lies on multigrid methods where the number of required operations to solve a linear equation system only increases linearly with problem size. Moreover, for badly conditioned problems quite an essential improvement is achieved by preconditioning. The third part of the numerical model discusses certain aspects of damage simulation which are closely related to the proposed grid discretization. The strong efficiency of the linear analysis can be maintained for damage simulation. This is achieved by a damage-controlled sequentially linear iteration scheme. Finally a study on the effective material behavior of heterogeneous bodies is presented. Especially the influence of inclusion shapes is examined. By means of altogether more than one hundred thousand random geometrical arrangements, the effective material behavior is statistically analyzed and assessed. N2 - Die wichtige Bedeutung moderner Simulationsverfahren in der mechanischen Analyse heterogener Festkörper wird eingangs ausführlich dargestellt. Dabei wird als Problem festgestellt, dass die erforderliche hochauflösende Analyse bereits für relativ kleine Körper an die Grenzen heutiger Rechenleistung stößt, sowohl bezüglich Speicherbedarf als auch akzeptablen Rechenaufwands. Ein weiteres Problem stellt die häufig unzureichend genaue geometrische Modellierung der Zusammensetzung heterogener Körper dar. Die vorliegende Arbeit führt eine systematische Kombination und Anpassung von gitterbasierten Methoden ein, um dadurch eine wesentlich höhere Auflösung in der numerischen Analyse heterogener Körper zu erzielen. Gitterverfahren eignen sich ebenfalls ausgezeichnet, um effiziente und numerisch stabile Algorithmen zur flexiblen geometrischen Modellierung zu entwickeln. Ein Schlüsselaspekt stellt ein gleichmäßiges Datenmanagement für Gitter dar, welches dafür eingesetzt werden kann, um den Aufwand und die Komplexität von nahezu allen beteiligten Methoden zu reduzieren. Ein neues Finite-Elemente Programm, namens Mulgrido, wurde eigens dafür entwickelt, um das vorgeschlagene Konzept konsistent zu realisieren und zu untersuchen. Einige Nachteile, die sich klassischerweise aus Gitterdiskretisierungen ergeben, werden gezielt durch modifizierte Verfahren korrigiert. Die gegenwärtige Arbeit gliedert sich in ein geometrisches Modell, ein mechanisches Modell und ein numerisches Modell. Das geometrische Modell beinhaltet neben Methoden der digitalen Bildverarbeitung, insbesondere sämtliche Verfahren zur künstlichen Generierung von Einschluss-Matrix Geometrien. Wesentliche Beiträge werden bezüglich variabler Form, Größenverteilung, Überschneidungsabfragen und Platzierung von Einschlüssen geleistet. Das mechanische Modell bereitet durch Grundlagen der Kontinuumsmechanik, der Homogenisierung und der Schädigungsmodellierung auf eine numerische Umsetzung vor. Als erstes Thema des numerischen Modells wird eine besondere Umsetzung von B-Spline Finiten Elementen vorgestellt. Diese Finite Elemente können generisch für eine beliebige Ordnung k der B-Splines erzeugt werden. Für homogene Körper verfügen diese somit über beliebig skalierbare Approximationseigenschaften. Mittels des Konzepts mehrphasiger Finite Elemente in Kombination mit Übergangszonen entlang von Materialgrenzen gelingt eine hochwertige Erweiterung für heterogene Körper. Durch die Formulierung auf Elementebene, kann die Speicherung der globalen Steifigkeitsmatrix und somit wesentlicher Speicherplatz eingespart werden. Dies ist möglich in Kombination mit iterativen Lösungsverfahren, die das zweite Thema des numerischen Modells darstellen. Dabei liegt der Fokus auf Mehrgitterverfahren. Diese zeichnen sich dadurch aus, dass die Anzahl der erforderlichen Operationen um ein lineares Gleichungssystem zu lösen, nur linear mit der Problemgröße ansteigt. Durch Vorkonditionierung wird für schlecht konditionierte Probleme eine ganz wesentliche Verbesserung erreicht. Als drittes Thema des numerischen Modells werden Aspekte der Schädigungssimulation diskutiert, die in engem Zusammenhang mit der Gitterdiskretisierung stehen. Die hohe Effizienz der linearen Analyse kann durch ein schädigungskontrolliertes, schrittweise lineares Iterationsschema für die Schädigungsanalyse aufrecht erhalten werden. Abschließend wird eine Studie über das effektive Materialverhalten heterogener Körper vorgestellt. Insbesondere wird der Einfluss der Form von Einschlüssen untersucht. Mittels insgesamt weit über hunderttausend zufälliger geometrischer Anordnungen wird das effektive Materialverhalten statistisch analysiert und bewertet. T2 - Gitterbasierte Verfahren zur mechanischen Analyse heterogener Festkörper KW - B-Spline KW - Finite-Elemente-Methode KW - Mehrgitterverfahren KW - Homogenisieren KW - Schädigung KW - Festkörpermechanik KW - Numerische Mathematik KW - B-Spline Finite Elemente KW - Homogenisierung KW - mehrphasig KW - Lösungsverfahren KW - Modellierung KW - B-spline KW - finite element KW - multigrid KW - multiphase KW - effective properties Y1 - 2006 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:gbv:wim2-20070830-9185 ER - TY - THES A1 - López Zermeño, Jorge Alberto T1 - Isogeometric and CAD-based methods for shape and topology optimization: Sensitivity analysis, Bézier elements and phase-field approaches N2 - The Finite Element Method (FEM) is widely used in engineering for solving Partial Differential Equations (PDEs) over complex geometries. To this end, it is required to provide the FEM software with a geometric model that is typically constructed in a Computer-Aided Design (CAD) software. However, FEM and CAD use different approaches for the mathematical description of the geometry. Thus, it is required to generate a mesh, which is suitable for FEM, based on the CAD model. Nonetheless, this procedure is not a trivial task and it can be time consuming. This issue becomes more significant for solving shape and topology optimization problems, which consist in evolving the geometry iteratively. Therefore, the computational cost associated to the mesh generation process is increased exponentially for this type of applications. The main goal of this work is to investigate the integration of CAD and CAE in shape and topology optimization. To this end, numerical tools that close the gap between design and analysis are presented. The specific objectives of this work are listed below: • Automatize the sensitivity analysis in an isogeometric framework for applications in shape optimization. Applications for linear elasticity are considered. • A methodology is developed for providing a direct link between the CAD model and the analysis mesh. In consequence, the sensitivity analysis can be performed in terms of the design variables located in the design model. • The last objective is to develop an isogeometric method for shape and topological optimization. This method should take advantage of using Non-Uniform Rational B-Splines (NURBS) with higher continuity as basis functions. Isogeometric Analysis (IGA) is a framework designed to integrate the design and analysis in engineering problems. The fundamental idea of IGA is to use the same basis functions for modeling the geometry, usually NURBS, for the approximation of the solution fields. The advantage of integrating design and analysis is two-fold. First, the analysis stage is more accurate since the system of PDEs is not solved using an approximated geometry, but the exact CAD model. Moreover, providing a direct link between the design and analysis discretizations makes possible the implementation of efficient sensitivity analysis methods. Second, the computational time is significantly reduced because the mesh generation process can be avoided. Sensitivity analysis is essential for solving optimization problems when gradient-based optimization algorithms are employed. Automatic differentiation can compute exact gradients, automatically by tracking the algebraic operations performed on the design variables. For the automation of the sensitivity analysis, an isogeometric framework is used. Here, the analysis mesh is obtained after carrying out successive refinements, while retaining the coarse geometry for the domain design. An automatic differentiation (AD) toolbox is used to perform the sensitivity analysis. The AD toolbox takes the code for computing the objective and constraint functions as input. Then, using a source code transformation approach, it outputs a code for computing the objective and constraint functions, and their sensitivities as well. The sensitivities obtained from the sensitivity propagation method are compared with analytical sensitivities, which are computed using a full isogeometric approach. The computational efficiency of AD is comparable to that of analytical sensitivities. However, the memory requirements are larger for AD. Therefore, AD is preferable if the memory requirements are satisfied. Automatic sensitivity analysis demonstrates its practicality since it simplifies the work of engineers and designers. Complex geometries with sharp edges and/or holes cannot easily be described with NURBS. One solution is the use of unstructured meshes. Simplex-elements (triangles and tetrahedra for two and three dimensions respectively) are particularly useful since they can automatically parameterize a wide variety of domains. In this regard, unstructured Bézier elements, commonly used in CAD, can be employed for the exact modelling of CAD boundary representations. In two dimensions, the domain enclosed by NURBS curves is parameterized with Bézier triangles. To describe exactly the boundary of a two-dimensional CAD model, the continuity of a NURBS boundary representation is reduced to C^0. Then, the control points are used to generate a triangulation such that the boundary of the domain is identical to the initial CAD boundary representation. Thus, a direct link between the design and analysis discretizations is provided and the sensitivities can be propagated to the design domain. In three dimensions, the initial CAD boundary representation is given as a collection of NURBS surfaces that enclose a volume. Using a mesh generator (Gmsh), a tetrahedral mesh is obtained. The original surface is reconstructed by modifying the location of the control points of the tetrahedral mesh using Bézier tetrahedral elements and a point inversion algorithm. This method offers the possibility of computing the sensitivity analysis using the analysis mesh. Then, the sensitivities can be propagated into the design discretization. To reuse the mesh originally generated, a moving Bézier tetrahedral mesh approach was implemented. A gradient-based optimization algorithm is employed together with a sensitivity propagation procedure for the shape optimization cases. The proposed shape optimization approaches are used to solve some standard benchmark problems in structural mechanics. The results obtained show that the proposed approach can compute accurate gradients and evolve the geometry towards optimal solutions. In three dimensions, the moving mesh approach results in faster convergence in terms of computational time and avoids remeshing at each optimization step. For considering topological changes in a CAD-based framework, an isogeometric phase-field based shape and topology optimization is developed. In this case, the diffuse interface of a phase-field variable over a design domain implicitly describes the boundaries of the geometry. The design variables are the local values of the phase-field variable. The descent direction to minimize the objective function is found by using the sensitivities of the objective function with respect to the design variables. The evolution of the phase-field is determined by solving the time dependent Allen-Cahn equation. Especially for topology optimization problems that require C^1 continuity, such as for flexoelectric structures, the isogeometric phase field method is of great advantage. NURBS can achieve the desired continuity more efficiently than the traditional employed functions. The robustness of the method is demonstrated when applied to different geometries, boundary conditions, and material configurations. The applications illustrate that compared to piezoelectricity, the electrical performance of flexoelectric microbeams is larger under bending. In contrast, the electrical power for a structure under compression becomes larger with piezoelectricity. T3 - ISM-Bericht // Institut für Strukturmechanik, Bauhaus-Universität Weimar - 2022,4 KW - CAD KW - Gestaltoptimierung KW - Topologieoptimierung KW - Isogeometrische Analyse KW - Finite-Elemente-Methode KW - Computer-Aided Design KW - Shape Optimization KW - Topology Optimization KW - Isogeometric Analysis KW - Finite Element Method Y1 - 2022 U6 - http://nbn-resolving.de/urn/resolver.pl?urn:nbn:de:gbv:wim2-20220831-47102 ER -