50.31 Technische Mechanik
In this paper we present a theoretical background for a coupled analytical–numerical approach to model a crack propagation process in two-dimensional bounded domains. The goal of the coupled analytical–numerical approach is to obtain the correct solution behaviour near the crack tip by help of the analytical solution constructed by using tools of complex function theory and couple it continuously with the finite element solution in the region far from the singularity. In this way, crack propagation could be modelled without using remeshing. Possible directions of crack growth can be calculated through the minimization of the total energy composed of the potential energy and the dissipated energy based on the energy release rate. Within this setting, an analytical solution of a mixed boundary value problem based on complex analysis and conformal mapping techniques is presented in a circular region containing an arbitrary crack path. More precisely, the linear elastic problem is transformed into a Riemann–Hilbert problem in the unit disk for holomorphic functions. Utilising advantages of the analytical solution in the region near the crack tip, the total energy could be evaluated within short computation times for various crack kink angles and lengths leading to a potentially efficient way of computing the minimization procedure. To this end, the paper presents a general strategy of the new coupled approach for crack propagation modelling. Additionally, we also discuss obstacles in the way of practical realisation of this strategy.
The purpose of this study is to develop self-contained methods for obtaining smooth meshes which are compatible with isogeometric analysis (IGA). The study contains three main parts. We start by developing a better understanding of shapes and splines through the study of an image-related problem. Then we proceed towards obtaining smooth volumetric meshes of the given voxel-based images. Finally, we treat the smoothness issue on the multi-patch domains with C1 coupling. Following are the highlights of each part.
First, we present a B-spline convolution method for boundary representation of voxel-based images. We adopt the filtering technique to compute the B-spline coefficients and gradients of the images effectively. We then implement the B-spline convolution for developing a non-rigid images registration method. The proposed method is in some sense of “isoparametric”, for which all the computation is done within the B-splines framework. Particularly, updating the images by using B-spline composition promote smooth transformation map between the images. We show the possible medical applications of our method by applying it for registration of brain images.
Secondly, we develop a self-contained volumetric parametrization method based on the B-splines boundary representation. We aim to convert a given voxel-based data to a matching C1 representation with hierarchical cubic splines. The concept of the osculating circle is employed to enhance the geometric approximation, where it is done by a single template and linear transformations (scaling, translations, and rotations) without the need for solving an optimization problem. Moreover, we use the Laplacian smoothing and refinement techniques to avoid irregular meshes and to improve mesh quality. We show with several examples that the method is capable of handling complex 2D and 3D configurations. In particular, we parametrize the 3D Stanford bunny which contains irregular shapes and voids.
Finally, we propose the B´ezier ordinates approach and splines approach for C1 coupling. In the first approach, the new basis functions are defined in terms of the B´ezier Bernstein polynomials. For the second approach, the new basis is defined as a linear combination of C0 basis functions. The methods are not limited to planar or bilinear mappings. They allow the modeling of solutions to fourth order partial differential equations (PDEs) on complex geometric domains, provided that the given patches are G1
continuous. Both methods have their advantages. In particular, the B´ezier approach offer more degree of freedoms, while the spline approach is more computationally efficient. In addition, we proposed partial degree elevation to overcome the C1-locking issue caused by the over constraining of the solution space. We demonstrate the potential of the resulting C1 basis functions for application in IGA which involve fourth order PDEs such as those appearing in Kirchhoff-Love shell models, Cahn-Hilliard phase field application, and biharmonic problems.
Numerische Berechnung von Mauerwerkstrukturen in homogenen und diskreten Modellierungsstrategien
(2004)
Im Zentrum der Arbeit stehen die Entwicklung, Verifikation, Implementierung und Leistungsfähigkeit numerischer Berechnungsmodelle für Mauerwerk im Rahmen der Kontinuums- und Diskontinuumsmechanik. Makromodelle beschreiben das Mauerwerk als verschmiertes Ersatzkontinuum. Mikromodelle berücksichtigen durch die Modellierung der einzelnen Steine und Fugen die Struktur des Mauerwerkverbandes. Soll darüber hinaus der durch die Querdehnungsinteraktion zwischen Stein und Mörtel hervorgerufene heterogene Spannungszustand im Mauerwerk abgebildet werden, so ist ein detailliertes Mikromodell, welches Steine und Fugen in ihren exakten geometrischen Dimensionen berücksichtigt, erforderlich. Demgegenüber steht die vereinfachte Mikromodellierung, bei der die Fugen mit Hilfe von Kontaktalgorithmen beschrieben werden. Im Rahmen der Makromodellierung werden neue räumliche Materialmodelle für verschiedene ein- und mehrschalige Mauerwerkarten hergeleitet. Die vorgestellten Modelle berücksichtigen die Anisotropie der Steifigkeiten, der Festigkeiten sowie des Ver- und Entfestigungsverhaltens. Die numerische Implementation erfolgt mit Hilfe moderner elastoplastischer Algorithmen im Rahmen der impliziten Finite Element Methode in das Programm ANSYS. Innerhalb der detaillierten Mikromodellierung wird ein neues, aus Materialbeschreibungen für Stein, Mörtel sowie deren Verbund bestehendes nichtlineares Berechnungsmodell entwickelt und in das Programm ANSYS implementiert. Die diskontinuumsmechanische Beschreibung von Mauerwerk im Rahmen der vereinfachten Mikromodellierung erfolgt unter Verwendung der expliziten Distinkt Element Methode mit Hilfe der Programme UDEC und 3DEC. An praktischen Beispielen werden Probleme der Tragfähigkeitsbewertung gemauerter Bogenbrücken, Möglichkeiten zur Bewertung vorhandener Rissbildungen und Schädigungen an historischen Mauerwerkstrukturen und Traglastberechnungen an gemauerten Stützen ausgewertet und analysiert.