Refine
Document Type
- Article (27)
- Conference Proceeding (24)
- Doctoral Thesis (17)
- Master's Thesis (4)
- Diploma Thesis (1)
Institute
- Professur Informatik im Bauwesen (41)
- Institut für Strukturmechanik (ISM) (17)
- Professur Baumechanik (3)
- Professur Stahlbau (3)
- Bauhaus-Institut für zukunftsweisende Infrastruktursysteme (b.is) (2)
- Graduiertenkolleg 1462 (2)
- Professur Stahl- und Hybridbau (2)
- Professur Angewandte Mathematik (1)
- Professur Baustatik und Bauteilfestigkeit (1)
- Professur Modellierung und Simulation - Mechanik (1)
Keywords
- Finite-Elemente-Methode (73) (remove)
The nonlinear behavior of concrete can be attributed to the propagation of microcracks within the heterogeneous internal material structure. In this thesis, a mesoscale model is developed which allows for the explicit simulation of these microcracks. Consequently, the actual physical phenomena causing the complex nonlinear macroscopic behavior of concrete can be represented using rather simple material formulations. On the mesoscale, the numerical model explicitly resolves the components of the internal material structure. For concrete, a three-phase model consisting of aggregates, mortar matrix and interfacial transition zone is proposed. Based on prescribed grading curves, an efficient algorithm for the generation of three-dimensional aggregate distributions using ellipsoids is presented. In the numerical model, tensile failure of the mortar matrix is described using a continuum damage approach. In order to reduce spurious mesh sensitivities, introduced by the softening behavior of the matrix material, nonlocal integral-type material formulations are applied. The propagation of cracks at the interface between aggregates and mortar matrix is represented in a discrete way using a cohesive crack approach. The iterative solution procedure is stabilized using a new path following constraint within the framework of load-displacement-constraint methods which allows for an efficient representation of snap-back phenomena. In several examples, the influence of the randomly generated heterogeneous material structure on the stochastic scatter of the results is analyzed. Furthermore, the ability of mesoscale models to represent size effects is investigated. Mesoscale simulations require the discretization of the internal material structure. Compared to simulations on the macroscale, the numerical effort and the memory demand increases dramatically. Due to the complexity of the numerical model, mesoscale simulations are, in general, limited to small specimens. In this thesis, an adaptive heterogeneous multiscale approach is presented which allows for the incorporation of mesoscale models within nonlinear simulations of concrete structures. In heterogeneous multiscale models, only critical regions, i.e. regions in which damage develops, are resolved on the mesoscale, whereas undamaged or sparsely damage regions are modeled on the macroscale. A crucial point in simulations with heterogeneous multiscale models is the coupling of sub-domains discretized on different length scales. The sub-domains differ not only in the size of the finite elements but also in the constitutive description. In this thesis, different methods for the coupling of non-matching discretizations - constraint equations, the mortar method and the arlequin method - are investigated and the application to heterogeneous multiscale models is presented. Another important point is the detection of critical regions. An adaptive solution procedure allowing the transfer of macroscale sub-domains to the mesoscale is proposed. In this context, several indicators which trigger the model adaptation are introduced. Finally, the application of the proposed adaptive heterogeneous multiscale approach in nonlinear simulations of concrete structures is presented.
Framed-tube system with multiple internal tubes is analysed using an orthotropic box beam analogy approach in which each tube is individually modelled by a box beam that accounts for the flexural and shear deformations, as well as the shear-lag effects. A simple numerical modeling technique is proposed for estimating the shear-lag phenomenon in tube structures with multiple internal tubes. The proposed method idealizes the framed-tube structures with multiple internal tubes as equivalent multiple tubes, each composed of four equivalent orthotropic plate panels. The numerical analysis is based on the minimum potential energy principle in conjunction with the variational approach. The shear-lag phenomenon of such structures is studied taking into account the additional bending moments in the tubes. A detailed work is carried out through the numerical analysis of the additional bending moment. The moment factor is further introduced to identify the shear lag phenomenon along with the additional moment.
The steel structure design codes require to check up the member strength when evaluating plastic deformations. The model of perfectly plastic material is accepted. The strength criteria for simple cross-sections (I section, etc.) of steel members are given in design codes. The analytical strength criteria for steel cross-sections and numerical approaches based on stepwise procedure are investigated in many articles. Another way for checking the carrying capacity of cross-sections is the use of methods that are applied for defining strain-deformed state of elastic perfectly plastic systems. In this paper non-iterative methods are suggested for checking strength of cross-sections. Carrying capacity of cross section is verified according to extremum principle of plastic fail under monotonically loading and the strain-deformed state of cross-section is defined according to extremum energy principals of elastic potential of residual stresses and complementary work of residual displacements. The mathematical expressions of these principals for discrete cross-section are formulated as problems of convex mathematical programming. The cross-section of steel member using finite element method is divided into free form plane elements. The constant distribution of stresses along the finite element is accepted. The relationships of finite elements for static formulation of the problem are formed so, that kinematics formulation relationships could be obtained in a formal way using the theory of duality. Numerical examples of determination of cross-section strength, composition of interactive curves and composition of moment-curvature curves for different axial force levels are presented.
This paper presents the combination of two different parallelization environments, OpenMP and MPI, in one numerical simulation tool. The computation of the system matrices and vectors is parallelized with OpenMP and the solution of the system of equations is done with the MPIbased solver MUMPS. The efficiency of both algorithms is shown on several linear and nonlinear examples using the Finite Element Method and a meshless discretization technique.
Creation of hierarchical sequence of the plastic and viscoplastic models according to different levels of structure approximations is considered. Developed strategy of multimodel analysis, which consists of creation of the inelastic models library, determination of selection criteria system and caring out of multivariant sequential clarifying computations, is described. Application of the multimodel approach in numerical computations has demonstrated possibility of reliable prediction of stress-strain response under wide variety of combined nonproportional loading.
Die vorliegende Arbeit beschäftigt sich mit der Berechnung der Sicherheit von Strukturen mit sowohl geometrisch als auch physikalisch nichtlinearem Verhalten. Die Berechnung der Versagenswahrscheinlichkeit einer Struktur mit Hilfe von Monte-Carlo-Simulationsmethoden erfordert, dass die Funktion der Strukturantwort implizit berechnet wird, zum Beispiel durch nichtlineare Strukturanalysen für jede Realisation der Zufallsvariablen. Die Strukturanalysen bilden jedoch den Hauptanteil am Berechnungsaufwand der Zuverlässigkeitsanalyse, so dass die Analyse von realistischen Strukturen mit nichtlinearem Verhalten durch die begrenzten Computer-Ressourcen stark eingeschränkt ist. Die klassischen Antwortflächenverfahren approximieren die Funktion der Strukturantwort oder aber die Grenzzustandsfunktion durch Polynome niedriger Ordnung. Dadurch ist für die Auswertung des Versagens-Kriteriums nur noch von Interesse, ob eine Realisation der Basisvariablen innerhalb oder außerhalb des von der Antwortflächenfunktion gebildeten Raumes liegt - die Strukturanalyse kann dann entfallen. Bei stark nichtlinearen Grenzzustandsfunktionen versagt die polynomiale Approximation. Das directional sampling neigt bei Problemen mit vielen Zufallsvariablen zu einem systematischen Fehler. Das adaptive importance directional sampling dagegen beseitigt diesen Fehler, verschenkt jedoch Informationen über den Verlauf der Grenzzustandsfunktion, da die aufgefundenen Stützstellen aus den vorangegangenen Simulationsläufen nicht berücksichtigt werden können. Aus diesem Grund erscheint eine Kombination beider Simulationsverfahren und eine Interpolation mittels einer Antwortfläche geeignet, diese Probleme zu lösen. Dies war die Motivation für die Entwicklung eines Verfahren der adaptiven Simulation der Einheitsvektoren und anschließender Interpolation der Grenzzustandsfunktion durch eine Antwortflächenfunktion. Dieses Vorgehen stellt besondere Anforderungen an die Antwortflächenfunktion. Diese muss flexibel genug sein, um stark nichtlineare Grenzzustandsfunktionen beliebig genau annähern zu können. Außerdem sollte die Anzahl der verarbeitbaren Stützstellen nicht begrenzt sein. Auch ist zu berücksichtigen, dass die Ermittlung der Stützstellen auf der Grenzzustandsfunktion nicht regelmäßig erfolgt. Die in dieser Arbeit entwickelten Methoden der lokalen Interpolation der Grenzzustandsfunktion durch Normalen-Hyperebenen bzw. sekantialen Hyperebenen und der sowohl lokalen als auch globalen Interpolation durch gewichtete Radien erfüllen diese Anforderungen. ungen. dieser Arbeit entwickelten Methoden der lokalen Interpolation der Grenzzustandsfunktion durch Normalen-Hyperebenen bzw. sekantialen Hyperebenen und der sowohl lokalen als auch globalen Interpolation durch gewichtete Radien erfüllen diese Anforderungen.
Für geometrisch imperfekte Strukturen wird die Versagenswahrscheinlichkeit bezüglich Stabilitätskriterien bestimmt. Eine probabilistische Beschreibung der geometrischen Imperfektionen erfolgt mit skalaren ortsdiskretisierten Zufallsfeldern. Die Stabilitätsberechnungen werden mit der Finite Elemente Methode durchgeführt. Ausgangspunkt der Berechnung ist eine systematische Formulierung probabilistisch gewichteter Imperfektionsformen durch eine Eigenwertzerlegung der Kovarianzmatrix. Wenn mit einer strukturmechanisch orientierten Sensitivitätsanalyse ein Unterraum zur näherungsweisen Beschreibung des probabilistischen Strukturverhaltens gefunden wird, kann die Versagenswahrscheinlichkeit numerisch sehr effizient durch ein Interaktionsmodell bestimmt werden. Es zeigte sich, daß dies genau dann möglich ist, wenn die Beulform merklich im Imperfektionsfeld enthalten ist. Die Imperfektionsform am Bemessungspunkt entspricht dann, unabhängig vom Lastniveau, gerade der Beulform. Wenn die Beulform im Imperfektionsfeld einen untergeordneten Beitrag liefert, erscheint eine Reduktion des stochastischen Problems auf wenige Zufallsvariablen dagegen nicht möglich.
When it comes to monitoring of huge structures, main issues are limited time, high costs and how to deal with the big amount of data. In order to reduce and manage them, respectively, methods from the field of optimal design of experiments are useful and supportive. Having optimal experimental designs at hand before conducting any measurements is leading to a highly informative measurement concept, where the sensor positions are optimized according to minimal errors in the structures’ models. For the reduction of computational time a combined approach using Fisher Information Matrix and mean-squared error in a two-step procedure is proposed under the consideration of different error types. The error descriptions contain random/aleatoric and systematic/epistemic portions. Applying this combined approach on a finite element model using artificial acceleration time measurement data with artificially added errors leads to the optimized sensor positions. These findings are compared to results from laboratory experiments on the modeled structure, which is a tower-like structure represented by a hollow pipe as the cantilever beam. Conclusively, the combined approach is leading to a sound experimental design that leads to a good estimate of the structure’s behavior and model parameters without the need of preliminary measurements for model updating.
This study is focused on finite element analysis of a model comprising femur into which a femoral component of a total hip replacement was implanted. The considered prosthesis is fabricated from a functionally graded material (FGM) comprising a layer of a titanium alloy bonded to a layer of hydroxyapatite. The elastic modulus of the FGM was adjusted in the radial, longitudinal, and longitudinal-radial directions by altering the volume fraction gradient exponent. Four cases were studied, involving two different methods of anchoring the prosthesis to the spongy bone and two cases of applied loading. The results revealed that the FG prostheses provoked more SED to the bone. The FG prostheses carried less stress, while more stress was induced to the bone and cement. Meanwhile, less shear interface stress was stimulated to the prosthesis-bone interface in the noncemented FG prostheses. The cement-bone interface carried more stress compared to the prosthesis-cement interface. Stair climbing induced more harmful effects to the implanted femur components compared to the normal walking by causing more stress. Therefore, stress shielding, developed stresses, and interface stresses in the THR components could be adjusted through the controlling stiffness of the FG prosthesis by managing volume fraction gradient exponent.
Encapsulation-based self-healing concrete (SHC) is the most promising technique for providing a self-healing mechanism to concrete. This is due to its capacity to heal fractures effectively without human interventions, extending the operational life and lowering maintenance costs. The healing mechanism is created by embedding capsules containing the healing agent inside the concrete. The healing agent will be released once the capsules are fractured and the healing occurs in the vicinity of the damaged part. The healing efficiency of the SHC is still not clear and depends on several factors; in the case of microcapsules SHC the fracture of microcapsules is the most important aspect to release the healing agents and hence heal the cracks. This study contributes to verifying the healing efficiency of SHC and the fracture mechanism of the microcapsules. Extended finite element method (XFEM) is a flexible, and powerful discrete crack method that allows crack propagation without the requirement for re-meshing and has been shown high accuracy for modeling fracture in concrete. In this thesis, a computational fracture modeling approach of Encapsulation-based SHC is proposed based on the XFEM and cohesive surface technique (CS) to study the healing efficiency and the potential of fracture and debonding of the microcapsules or the solidified healing agents from the concrete matrix as well. The concrete matrix and a microcapsule shell both are modeled by the XFEM and combined together by CS. The effects of the healed-crack length, the interfacial fracture properties, and microcapsule size on the load carrying capability and fracture pattern of the SHC have been studied. The obtained results are compared to those obtained from the zero thickness cohesive element approach to demonstrate the significant accuracy and the validity of the proposed simulation. The present fracture simulation is developed to study the influence of the capsular clustering on the fracture mechanism by varying the contact surface area of the CS between the microcapsule shell and the concrete matrix. The proposed fracture simulation is expanded to 3D simulations to validate the 2D computational simulations and to estimate the accuracy difference ratio between 2D and 3D simulations. In addition, a proposed design method is developed to design the size of the microcapsules consideration of a sufficient volume of healing agent to heal the expected crack width. This method is based on the configuration of the unit cell (UC), Representative Volume Element (RVE), Periodic Boundary Conditions (PBC), and associated them to the volume fraction (Vf) and the crack width as variables. The proposed microcapsule design is verified through computational fracture simulations.
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 paper investigates accuracy of deflection predictions made by the finite element package ATENA and design code methods ACI and EC2. Deflections have been calculated for a large number of experimental reinforced concrete beams reported by three investigators. Statistical parameters have been established for each of the technique at different load levels, separately for the beams with small and moderate reinforcement ratio.
For modeling of singular fields of stresses and deformations in elasters with a crack is offered to use of three-dimesional a special finite element. Weak compessible of elasters is taken into account on the basis of threefold approximation of fields of displacements, deformations and function of volume change. At intensive cyclic loading of the elastomer constructions with a crack it is necessary to take into account warming and large deformations at the crack top. The stress-deformed state elasters with a crack is determined from the decision of a nonlinear problem by a modified method Newton-Kantorovich. Account stress intensity factors for a rectangular plate with a various arrangement of a through crack is executed. Process of development of a surface crack and dissipative warming in prismatic a element of shift is investigated.
Development and Analysis of Sparse Matrix Concepts for Finite Element Approximation on general Cells
(2004)
In engineering and computing, the finite element approximation is one of the most well-known computational solution techniques. It is a great tool to find solutions for mechanic, fluid mechanic and ecological problems. Whoever works with the finite element method will need to solve a large system of linear equations. There are different ways to find a solution. One way is to use a matrix decomposition technique such as LU or QR. The other possibility is to use an iterative solution algorithm like Conjugate Gradients, Gauß-Seidel, Multigrid Methods, etc. This paper will focus on iterative solvers and the needed storage techniques...
Die Eisenbahnbrücken des Lehrter Bahnhofs in Berlin - Ein ganzheitliches FE-Berechnungskonzept
(2000)
Der Komplexität moderner Brückenbauwerke scheinen die verwendeten Berechungsmodelle oft nicht angemessen. Tragwerksberechnungen basieren in vielen Fällen noch auf der Vorgehensweise, das Brückenbauwerk in Einzelbauteile zu zerlegen und mit unterschiedlichen Teilmodellen zu behandeln. Das erscheint, auch vor dem Hintergrund ständig wachsender Rechnerleistung, nicht mehr zeitgemäß. Dies gilt zum Beispiel auch für die gängige Praxis, flächenhafte Brückenüberbauten mit Balkenmodellen zu berechnen. Der vorliegende Beitrag stellt ein ganzheitliches Berech-nungskonzept vor, welches auf der Basis eines einzigen FE-Modells die Berechnung des Gesamtbauwerks erlaubt. Damit wird für alle Bauteile neben der Zustandsgrößenberechnung auch die Bemessung von Stahl- und Spannbetonbauteilen bis hin zu Nachweisen wie zur Beschränkung der Rissbreite geführt. Die Anwendung dieses Berechnungskonzeptes wird am Beispiel der Eisenbahnüberführung des neuen Lehrter Bahn-hofs in Berlin gezeigt. Das verwendete FE-Modell umfasst Baugrund, Fundamente, Stahl- bzw. Gußstahlunterkonstruktion sowie den Stahl- bzw. Spannbetonüberbau. Besonderheiten sind unter anderem die Modellierung des plattenbalkenartigen Überbaus durch exzentrische, vorspannbare Schalenelemente und das getrennte Vorhalten von tragwerks- und lastbezogenen Eingabefiles. Damit gelingt die sequentielle Erfassung unterschiedlicher Bettungsmoduli zur Simulation statischer und dynamischer Beanspruchungen, die Berücksichtigung des Anspannens und der Interaktion zwischen vorgespannten Stahlverbänden zur Aufnahme von Horizontallasten sowie die Berücksichtigung unterschiedlicher statischer Systeme bei der Herstellung des Spannbetonüberbaus.
Discrete-continual Finite Element Method of Analysis for Three-dimensional Curvilinear Structures
(2003)
This paper is devoted to discrete-continual finite element method (DCFEM) of analysis for three-dimensional curvilinear structures. Operational and variational formulations of the problem in the ring coordinate system are presented. The discrete-continual design model for structures with constant physical and geometrical parameters in longitudinal direction is offered on the basis of so-called curvilinear discrete-continual finite elements. Element coordinate system, approximation of nodal unknowns, construction of element nodal load vector are under consideration. Element system of differential equations is formulated with use of special generalized block-structured stiffness matrix of discrete-continual finite element. Local differential relations are formulated. Resultant multipoint boundary problem for system of ordinary differential equations is given. Method of analytical solution of multipoint boundary problems in structural analysis is offered as well. Its major peculiarities include universality, computer-oriented algorithm involving theory of distributions, computational stability, optimal conditionality of resultant systems, partial Jordan decomposition of matrix of coefficients, eliminating necessity of calculation of root vectors. Brief information concerning developed software is provided.
Bei komplexen Gründungskonstruktionen sind Planungsfehler durch eine konsistente Modellierung vermeidbar. Manuelle Berechnungsmethoden ermöglichen im allgemeinen ein dreidimensionales Vorgehen nicht. Numerische Berechnungsmethoden, wie z.B. die Finite-Element-Methode, sind ein optimales Werkzeug zur ganzheitlichen Simulation des Problems. Die für die Finite-Element-Analyse notwendige Diskretisierung komplexer Bau- grundstrukturen ist manuell nicht zu bewältigen. Der vorliegende Beitrag zeigt wie ein Finite-Element-Modell automatisch aus einem geotechnischen Modell unter Berücksichtigung der spezifischen Anforderungen der Baugrund-Tragwerk-Struktur und des Bauablaufes erzeugt werden kann. Hierbei wird die Berücksichtigung der geometrischen und der mechanischen Besonderheiten bei der Netzgenerierung dargestellt.
The paper deals with the simulation of the non-linear and time dependent behaviour of complex structures in engineering. Such simulations have to provide high accuracy in the prediction of deformations and stability, by taking into account the long term influences of the non-linear behaviour of the material as well as the large deformation and contact conditions. The limiting factors of the computer simulation are the computer run time and the memory requirement during solving large scale problems. To overcome these problems we use a dynamic-explicit time integration procedure for the solution of the semi-discrete equations of motion, which is very suited for parallel processing. In the paper at first we give a brief review of the theoretical background of the mechanical modelling and the dynamic-explicit technique for the solution of the semi-discrete equations of motion. Then the concept of parallel processing will be discussed . A test example concludes the paper.
Die vorliegende Arbeit beschäftigt sich mit der dynamischen Analyse der Sprottetalbrücke infolge aufgetretener Asphaltschäden. Sie beinhaltet die Erstellung eines FE-Modells, der Darstellung der theoretischen Grundlagen der Dynamik sowie die Auswertung von berechneten Eigenformen und Asphaltspannungen unter Berücksichtigung der derzeit gültigen Normen.