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)
Die Methode der Finiten Elemente ist ein numerisches Verfahren zur Interpolation vorgegebener Werte und zur numerischen Approximation von Lösungen stationärer oder instationärer partieller Differentialgleichungen bzw. Systemen partieller Differentialgleichungen. Grundlage dieser Verfahren ist die Formulierung geeigneter Finiter Elemente und Finiter Element Zerlegungen. Finite Elemente besitzen in der Regel eine geometrische Basis bestehend aus Strecken im eindimensionalen, Drei- oder Vierecken im zweidimensionalen und Tetra- oder Hexaedern im dreidimensionalen euklidischen Raum, eine Menge von Freiheitsgraden und eine Basis von Funktionen. Die geometrische Basis eines Finiten Elements wird verallgemeinert als geometrische Zelle formuliert. Diese geschlossene geometrische Formulierung führt zu einer geometrieunabhängigen Definition der Basisfunktionen eines Finiten Elements in den Zellkoordinaten der geometrischen Zelle. Finite Elemente auf der Basis geometrischer Zellen werden als Bestandteile Finiter Element Zerlegungen in Finiten Element Interpolationen und Finiten Element Approximationen verwendet. Die Finiten Element Approximationen werden am Beispiel der 2-dimensionalen Diffusionsgleichung über das Standard-Galerkin-Verfahren ermittelt.
Hydro- und morphodynamischen Prozesse in Binnengewässern und im Küstennahbereich erzeugen hochkomplexe Phänomene. Zur Beurteilung der Entwicklung von Küstenzohnen, von Flussbetten sowie von Eingriffen des Menschen in Form von Schutzbauwerken sind geeignete numerische Modellwerkzeuge notwendig. Es wird ein holistischer Modellansatz zur Approximation gekoppelter Seegangs-, Strömungs- und Morphodynamischer Prozesse auf der Basis stabilisierter Finiter Elemente vorgestellt. Der Großteil der Modellgleichungen der Hydro- und Morphodynamik sind Transportgleichungen. Dem Transportcharakter dieser Gleichungen entsprechend wird ein stabilisiertes Finites Element Verfahren auf Dreiecken vorgestellt. Die vorgestellte Approximation entspricht einem streamline upwinding Petrov-Galerkin-Verfahrens für vektorwertige mehrdimensionale Probleme, bei dem der Fehler eines Standard-Galerkin-Verfahrens mit Hilfe eines Upwinding-Koeffizienten minimiert wird. Die Wahl des Upwinding-Koeffizienten ist übertragbar auf andere Problemklassen und basiert ausschließlich auf dem Charakter der zugrundeliegene Das Modell wurde für Seegangs- und Strömungs-Untersuchungen im Jade-Weser-Ästuar an der deutschen Nordseeküste eingesetzt.
A geometrical inclusion-matrix model for the finite element analysis of concrete at multiple scales
(2003)
This paper introduces a method to generate adequate inclusion-matrix geometries of concrete in two and three dimensions, which are independent of any specific numerical discretization. The article starts with an analysis on shapes of natural aggregates and discusses corresponding mathematical realizations. As a first prototype a two-dimensional generation of a mesoscale model is introduced. Particle size distribution functions are analysed and prepared for simulating an adequate three-dimensional representation of the aggregates within a concrete structure. A sample geometry of a three-dimensional test cube is generated and the finite element analysis of its heterogeneous geometry by a uniform mesh is presented. Concluding, aspects of a multiscale analysis are discussed and possible enhancements are proposed.
Iso-parametric finite elements with linear shape functions show in general a too stiff element behavior, called locking. By the investigation of structural parts under bending loading the so-called shear locking appears, because these elements can not reproduce pure bending modes. Many studies dealt with the locking problem and a number of methods to avoid the undesirable effects have been developed. Two well known methods are the >Assumed Natural Strain< (ANS) method and the >Enhanced Assumed Strain< (EAS) method. In this study the EAS method is applied to a four-node plane element with four EAS-parameters. The paper will describe the well-known linear formulation, its extension to nonlinear materials and the modeling of material uncertainties with random fields. For nonlinear material behavior the EAS parameters can not be determined directly. Here the problem is solved by using an internal iteration at the element level, which is much more efficient and stable than the determination via a global iteration. To verify the deterministic element behavior the results of common test examples are presented for linear and nonlinear materials. The modeling of material uncertainties is done by point-discretized random fields. To show the applicability of the element for stochastic finite element calculations Latin Hypercube Sampling was applied to investigate the stochastic hardening behavior of a cantilever beam with nonlinear material. The enhanced linear element can be applied as an alternative to higher-order finite elements where more nodes are necessary. The presented element formulation can be used in a similar manner to improve stochastic linear solid elements.
A realistic and reliable model is an important precondition for the simulation of revitalization tasks and the estimation of system properties of existing buildings. Thereby, the main focus lies on the parameter identification, the optimization strategies and the preparation of experiments. As usual structures are modeled by the finite element method. This as well as other techniques are based on idealizations and empiric material properties. Within one theory the parameters of the model should be approximated by gradually performed experiments and their analysis. This approximation method is performed by solving an optimization problem, which is usually non-convex, of high dimension and possesses a non-differentiable objective function. Therefore we use an optimization procedure based on genetic algorithms which was implemented by using the program package SLang...
A large-scale computer modeling and simulation method is presented for environmental flows in urban area. Several GIS and CAD data were used for the preparation of shape model and an automatic mesh generation method based on Delaunay method was developed. Parallel finite element method based on domain decomposition method was employed for the numerical simulation of natural phenomena. The present method was applied to the simulation of flood flow and wind flow in urban area. The present method is shown to be a useful planning and design tool for the natural disasters and the change of environments.
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.
Transport problems, as, for instance, the transport of sediment in hydraulic engineering and the transport of harmful substances through porous media, play an important role in many fields of civil engineering. Other examples include the dissipation of heat or sound as well as the simulation of traffic with macroscopic models. The contribution explains the analysis of the applicability of Voronoi-based finite volume methods for the approximation of solutions of transport problems. A special concern is the discretisation of the transport equation. Current limitations of the method as well as ideas for stabilisation are explained with examples.
Analysis System for Bridge Test (Chinese name abbr.: QLJC) is an application software specially designed for bridge test to analyze the static and dynamic character of bridge structures, calculate efficiency ratio of load test, pick up the results of observation points and so on. In this paper, research content, system design, calculation theory, characteristics and practical application of QLJC is introduced in detail.
The influence of vortex-induces vibrations on vertical tie rods has been proved as a determinant load factor in the lifetime-oriented dimensioning of arched steel bridges. Particularly, the welded connection plates between the suspenders and the arches often exhibit cracks induced primarily rods. In this context, the synchronization of the vortex-shedding to the rod motion in a critical wind velocity range, the so-called lock-in effect, is of essential interest.
The method of the finite elements is an adaptable numerical procedure for interpolation as well as for the numerical approximation of solutions of partial differential equations. The basis of these procedure is the formulation of suitable finite elements and element decompositions of the solution space. Classical finite elements are based on triangles or quadrangles in the two-dimensional space and tetrahedron or hexahedron in the threedimensional space. The use of arbitrary-dimensional convex and non-convex polyhedrons as the geometrical basis of finite elements increases the flexibility of generating finite element decompositions substantially and is sometimes the only way to get a clear decomposition...
In this paper, systematic analyses for the shoring systems installed to support the applied loads during construction are performed on the basis of the numerical approach. On the basis of a rigorous time-dependent analysis, structural behaviors of reinforced concrete (RC) frame structures according to the changes in design variables such as the types of shoring systems, shore stiffness and shore spacing are analyzed and discussed. The time-dependent deformations of concrete such as creep and shrinkage and construction sequences of frame structures are also taken into account to minimize the structural instability and to reach to an improved design of shoring system because these effects may increase the axial forces delivered to the shores. In advance, the influence of the column shortening effect, generally mentioned in a tall building structure, is analyzed. From many parametric studies, it has been finally concluded that the most effective shoring system in RC frame structures is 2S1R (two shores and one reshore) regardless of the changes in design variables.
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...
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.
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.
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 displacements and stresses in arch dams and their abutments are frequently determined with 20-node brick elements. The elements are distorted near the contact plane between the wall and the abutment. A cantilever beam testbed has been developed to investigate the consequences of this distortion. It is shown that the deterioration of the accuracy in the computed stresses is significant. A compatible 18-node wedge element with linear stress variation is developed as an alternative to the brick element. The shape of this element type is readily adapted to the shape of the contact plane. It is shown that the accuracy of the computed stresses in the vicinity of the contact plane is improved significantly by the use of wedge elements.
The primary objective of initial shape analysis of a cable stayed bridge is to calculate initial installation cable tension forces and to evaluate fabrication camber of main span and pylon providing the final longitudinal profile of the bridge at the end of construction. In addition, the initial cable forces depending on the alternation of the bridge’s shape can be obtained from the analysis, and will be used to provide construction safety during construction. In this research, we conducted numerical experiments for initial shape of Ko-ha bridge, which will be constructed in the near future, using three different typical methods such as continuous beam method, linear truss method, and IIMF (Introducing Initial Member Force) method
The worldwide growth of communication networks and associated technologies provide the basic infrastructure for new ways of executing the engineering process. Collaboration amongst team members seperated in time and location is of particular importance. Two broad themes can be recognized in research pertaining to distributed collaboration. One theme focusses on the technical and technological aspects of distributed work, while the other emphasises human aspects thereof. The case of finite element structural analysis in a distributed collaboratory is examined in this paper. An approach is taken which has its roots in human aspects of the structural analysis task. Based on experience of how structural engineers currently approach and execute this task while utilising standard software designed for use on local workstations only, criteria are stated for a software architechture that could support collaborative structural analysis. Aspects of a pilot application and the results of qualitative performance measurements are discussed.
The optimization of continuous structures requires careful attention to discretization errors. Compared to ordinary low order formulation (h-elements) in conjunction with an adaptive mesh refinement in each optimization step, the use of high order finite elements (so called p-elements) has several advantages. However, compared to the h-method a higher order finite element analysis program poses higher demands from a software engineering point of view. In this article the basics of an object oriented higher order finite element system especially tailored to the use in structural optimization is presented. Besides the design of the system, aspects related to the employed implementation language Java are discussed.
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.