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)
This thesis addresses an adaptive higher-order method based on a Geometry Independent Field approximatTion(GIFT) of polynomial/rationals plines over hierarchical T-meshes(PHT/RHT-splines).
In isogeometric analysis, basis functions used for constructing geometric models in computer-aided design(CAD) are also employed to discretize the partial differential equations(PDEs) for numerical analysis. Non-uniform rational B-Splines(NURBS) are the most commonly used basis functions in CAD. However, they may not be ideal for numerical analysis where local refinement is required.
The alternative method GIFT deploys different splines for geometry and numerical analysis. NURBS are utilized for the geometry representation, while for the field solution, PHT/RHT-splines are used. PHT-splines not only inherit the useful properties of B-splines and NURBS, but also possess the capabilities of local refinement and hierarchical structure. The smooth basis function properties of PHT-splines make them suitable for analysis purposes. While most problems considered in isogeometric analysis can be solved efficiently when the solution is smooth, many non-trivial problems have rough solutions. For example, this can be caused by the presence of re-entrant corners in the domain. For such problems, a tensor-product basis (as in the case of NURBS) is less suitable for resolving the singularities that appear since refinement propagates throughout the computational domain. Hierarchical bases and local refinement (as in the case of PHT-splines) allow for a more efficient way to resolve these singularities by adding more degrees of freedom where they are necessary. In order to drive the adaptive refinement, an efficient recovery-based error estimator is proposed in this thesis. The estimator produces a recovery solution which is a more accurate approximation than the computed numerical solution. Several two- and three-dimensional numerical investigations with PHT-splines of higher order and continuity prove that the proposed method is capable of obtaining results with higher accuracy, better convergence, fewer degrees of freedom and less computational cost than NURBS for smooth solution problems. The adaptive GIFT method utilizing PHT-splines with the recovery-based error estimator is used for solutions with discontinuities or singularities where adaptive local refinement in particular domains of interest achieves higher accuracy with fewer degrees of freedom. This method also proves that it can handle complicated multi-patch domains for two- and three-dimensional problems outperforming uniform refinement in terms of degrees of freedom and computational cost.
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.
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.
Global structural analyses in civil engineering are usually performed considering linear-elastic material behavior. However, for steel structures, a certain degree of plasticization depending on the member classification may be considered. Corresponding plastic analyses taking material nonlinearities into account are effectively realized using numerical methods. Frequently applied finite elements of two and three-dimensional models evaluate the plasticity at defined nodes using a yield surface, i.e. by a yield condition, hardening rule, and flow rule. Corresponding calculations are connected to a large numerical as well as time-consuming effort and they do not rely on the theoretical background of beam theory, to which the regulations of standards mainly correspond. For that reason, methods using beam elements (one-dimensional) combined with cross-sectional analyses are commonly applied for steel members in terms of plastic zones theories. In these approaches, plasticization is in general assessed by means of axial stress only. In this paper, more precise numerical representation of the combined stress states, i.e. axial and shear stresses, is presented and results of the proposed approach are validated and discussed.
Methods based on B-splines for model representation, numerical analysis and image registration
(2015)
The thesis consists of inter-connected parts for modeling and analysis using newly developed isogeometric methods. The main parts are reproducing kernel triangular B-splines, extended isogeometric analysis for solving weakly discontinuous problems, collocation methods using superconvergent points, and B-spline basis in image registration applications.
Each topic is oriented towards application of isogeometric analysis basis functions to ease the process of integrating the modeling and analysis phases of simulation.
First, we develop reproducing a kernel triangular B-spline-based FEM for solving PDEs. We review the triangular B-splines and their properties. By definition, the triangular basis function is very flexible in modeling complicated domains. However, instability results when it is applied for analysis. We modify the triangular B-spline by a reproducing kernel technique, calculating a correction term for the triangular kernel function from the chosen surrounding basis. The improved triangular basis is capable to obtain the results with higher accuracy and almost optimal convergence rates.
Second, we propose an extended isogeometric analysis for dealing with weakly discontinuous problems such as material interfaces. The original IGA is combined with XFEM-like enrichments which are continuous functions themselves but with discontinuous derivatives. Consequently, the resulting solution space can approximate solutions with weak discontinuities. The method is also applied to curved material interfaces, where the inverse mapping and the curved triangular elements are considered.
Third, we develop an IGA collocation method using superconvergent points. The collocation methods are efficient because no numerical integration is needed. In particular when higher polynomial basis applied, the method has a lower computational cost than Galerkin methods. However, the positions of the collocation points are crucial for the accuracy of the method, as they affect the convergent rate significantly. The proposed IGA collocation method uses superconvergent points instead of the traditional Greville abscissae points. The numerical results show the proposed method can have better accuracy and optimal convergence rates, while the traditional IGA collocation has optimal convergence only for even polynomial degrees.
Lastly, we propose a novel dynamic multilevel technique for handling image registration. It is application of the B-spline functions in image processing. The procedure considered aims to align a target image from a reference image by a spatial transformation. The method starts with an energy function which is the same as a FEM-based image registration. However, we simplify the solving procedure, working on the energy function directly. We dynamically solve for control points which are coefficients of B-spline basis functions. The new approach is more simple and fast. Moreover, it is also enhanced by a multilevel technique in order to prevent instabilities. The numerical testing consists of two artificial images, four real bio-medical MRI brain and CT heart images, and they show our registration method is accurate, fast and efficient, especially for large deformation problems.
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.
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.
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.
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.
Matrix-free voxel-based finite element method for materials with heterogeneous microstructures
(2019)
Modern image detection techniques such as micro computer tomography
(μCT), magnetic resonance imaging (MRI) and scanning electron microscopy (SEM) provide us with high resolution images of the microstructure of materials in a non-invasive and convenient way. They form the basis for the geometrical models of high-resolution analysis, so called image-based analysis.
However especially in 3D, discretizations of these models reach easily the size of 100 Mill. degrees of freedoms and require extensive hardware resources in terms of main memory and computing power to solve the numerical model. Consequently, the focus of this work is to combine and adapt numerical solution methods to reduce the memory demand first and then the computation time and therewith enable an execution of the image-based analysis on modern computer desktops. Hence, the numerical model is a straightforward grid discretization of the voxel-based (pixels with a third dimension) geometry which omits the boundary detection algorithms and allows reduced storage of the finite element data structure and a matrix-free solution algorithm.
This in turn reduce the effort of almost all applied grid-based solution techniques and results in memory efficient and numerically stable algorithms for the microstructural models. Two variants of the matrix-free algorithm are presented. The efficient iterative solution method of conjugate gradients is used with matrix-free applicable preconditioners such as the Jacobi and the especially suited multigrid method. The jagged material boundaries of the voxel-based mesh are smoothed through embedded boundary elements which contain different material information at the integration point and are integrated sub-cell wise though without additional boundary detection. The efficiency of the matrix-free methods can be retained.
The dynamic behaviour of shells, which are widely used in construction and mechanical engineering as critical components of machinery and 3-D structures, under static and dynamic loadings is described by system of deep nonlinear differential equations. Solution of these equations can be received with assistance of technique basing on a modern numerical algorithms and computer modeling.. The system of nonlinear differential equations of vibration of the shells is proposed taking into account the inertia forces in the tangential and normal directions. Its solution is based on combination of parameter prolongation method, finite-difference method and the Newton-Kantorovich iterative algorithm that allows plotting the loading trajectories and determination of bifurcation points on them. Package of Applied Programs >SEVSOR< is a computation means to be used in research of deformation, stability and vibration in thin axically-symmetric shells of complicated shape Input data include information on shell geometry, physical and mechanical properties, bearing conditions, types of loadings and load application. Frame output of motion forms in real time or either in decelerated or accelerated time scales for creating cartoons or video films is used for analysis of the compound dynamic processes in shell-type structures.
According to Eurocode, the computation of bending strength for steel cantilever beams is a straightforward process. The approach is based on an Ayrton-Perry formula adaptation of buckling curves for steel members in compression, which involves the computation of an elastic critical buckling load for considering the instability. NCCI documents offer a simplified formula to determine the critical bending moment for cantilevers beams with symmetric cross-section. Besides the NCCI recommendations, other approaches, e.g. research literature or Finite-Element-Analysis, may be employed to determine critical buckling loads. However, in certain cases they render different results. Present paper summarizes and compares the abovementioned analytical and numerical approaches for determining critical loads and it exemplarily analyses corresponding cantilever beam capacities using numerical approaches based on plastic zones theory (GMNIA).
This master thesis explores an important and under-researched topic on the so-called bridging of length scales (from >meso< to >macro<), with the concept of homogenization in which the careful characterization of mechanical response requires that the developed material model >bridge< the representations of events that occur at two different scales. The underlying objective here is to efficiently incorporate material length scales in the classical continuum plasticity/damage theories through the concept of homogenization theory. The present thesis is devoted to computational modeling of heterogeneous materials, primarily to matrix-inclusion type of materials. Considerations are focused predominantly on the elastic and damage behavior as a response to quasistatic mechanical loading. Mainly this thesis focuses to elaborate a sound numerical homogenization model which accounts for the prediction of overall properties with the application of different types of boundary conditions namely: periodic, homogeneous and mixed type of boundary conditions over two-dimensional periodic and non-periodic RVEs and three-dimensional non-periodic RVEs. Identification of the governing mechanisms and assessing their effect on the material behavior leads one step further. Bringing together this knowledge with service requirements allows for functional oriented materials design. First, this thesis gives attention on providing the theoretical basic mechanisms involved in homogenization techniques and a survey will be made on existing analytical methods available in literature. Second, the proposed frameworks are implemented in the well known finite element software programs ANSYS and SLang. Simple and efficient algorithms in FORTRAN are developed for automated microstructure generation using RSA algorithm in order to perform a systematic numerical testing of microstructures of composites. Algorithms are developed to generate constraint equations in periodic boundary conditions and different displacements applied spatially over the boundaries of the RVE in homogeneous boundary conditions. Finally, nonlinear simulations are performed at mesolevel, by considering continuum scalar damage behavior of matrix material with the linear elastic behavior of aggregates with the assumption of rigid bond between constituents.
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.
Wirklichkeitsnahe Erfassung und Beschreibung des Trag- und Verformungsverhaltens von Strukturen baulicher Anlagen hat in den letzten Jahrzehnten ständig an Bedeutung gewonnen. Konstruktionen im Hoch- und Industriebau werden zunehmend multifunktional genutzt - die >Grenzen< zwischen Bauwerk und Tragwerk, zwischen Hüll- und Tragkonstruktion lösen sich auf. Werden raumabschließende Elemente (Wände, Decken, Dächer) gleichzeitig als Tragelemente und wärme- und schalldämmende Konstruktionen ausgeführt, so entstehen beispielsweise Sandwichplatten, deren Schichten sehr stark differierende Materialeigenschaften aufweisen. Beim Aufbau des FEM-Modells für vielschichtige Schalen können die Formänderungshypothesen für jede Schicht einzeln als auch für die Schale insgesamt gegeben werden. Im ersten Fall ist der Knotenfreiheitsgrad von der Schichtenzahl abhängig, im zweiten Fall nicht. Im weiteren wird eine Formänderungshypothese für das Schichtenpaket angenommen. Ausgegangen wird von den Gleichungen der 3D-Elastizitätstheorie. Die Berücksichtigung der Querkraftschubverformungen ergibt die Möglichkeit einer adäquaten Beschreibung der Verformungen sowohl dünner Schalen als auch von Schalen mittlerer Dicke; die Berechnung der Krümmungen und der LAMEschen Parameter der Bezugsfläche zu umgehen, was für komplizierte Schalenformen eine selbständige Aufgabe ist; eines natürlichen Übergangs von homogenen zu geschichteten Schalen. Das vielschichtige isoparametrische Schalen-FE wird vorgestellt, seine Implementierung in das in Entwicklung befindliche Programmsystem SLANG wird vorbereitet.
We propose an enhanced iterative scheme for the precise reconstruction of piezoelectric material parameters from electric impedance and mechanical displacement measurements. It is based on finite-element simulations of the full three-dimensional piezoelectric equations, combined with an inexact Newton or nonlinear Landweber iterative inversion scheme. We apply our method to two piezoelectric materials and test its performance. For the first material, the manufacturer provides a full data set; for the second one, no material data set is available. For both cases, our inverse scheme, using electric impedance measurements as input data, performs well.
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.
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.
In der täglichen Ingenieurpraxis werden in zunehmenden Maße numerische Analysen im Rahmen der Finite-Elemente-Methode auch zur Untersuchung stabilitätsgefährdeter Strukturen eingesetzt. Für die aktuelle Praxis, insbesondere im konstruktiven Stahlbau, ist jedoch festzustellen, dass zwischen der fortgeschrittenen Theorie und dem Niveau der praktischen Anwendung numerischer Stabilitätsanalysen eine große Kluft besteht. Aus praktischer Sicht erscheint es unumgänglich, die weiter wachsende Diskrepanz zwischen den umfangreichen theoretischen Möglichkeiten und der gegenwärtigen Praxis abzubauen. Damit steht der praktisch tätige Ingenieur vor der Aufgabe, sein Wissen auf dem Gebiet numerischer Stabilitätsanalysen zu vertiefen und bereits vorhandene FE-Programme um Berechnungsalgorithmen für umfassende numerische Stabilitätsanalysen zu erweitern. Dafür werden in der Arbeit die Grundlagen einer FEM- orientierten modernen Stabilitätstheorie einheitlich und aus Sicht einer praktischen Anwendung aufbereitet. Die Darstellung von realisierten programmtechnischen Umsetzungen für erweiterte Analysenmethoden wie Nachbeulanalysen, Pfadwechsel und Approximationen imperfekter Pfade ermöglicht eine Erweiterung des Methodenvorrates. Die innerhalb der Arbeit untersuchten Beispiele zeigen, dass durch die Anwendung der behandelten Verfahren das Tragverhalten einer stabilitätsgefährdeten Struktur wesentlich besser eingeschätzt werden kann als bei Beschränkung auf die herkömmlichen Analysemethoden.
Auf der Basis der Literaturrecherche wird in dieser Arbeit eine 5-lagige MAG-geschweißte Stumpfnaht an austenitisch-ferritischen Stahl X2CrNiMoN22-5-3 (Duplex-Stahl 1.4462) mit dem FE-Programm „SYSWELD®“ simuliert. Die Berech-nungen der Temperaturfelder werden unter der Berücksichtigung sowohl von tempe-raturunabhängigen als auch temperaturabhängigen thermophysikalischen Material-eigenschaften am drei-dimensionalen und zwei-dimensionalen Modell durchgeführt. Die berechneten Temperatur-Zeit-Verläufe und Gefügeumwandlungen beim MAG-Schweißen der Stumpfnaht werden hinsichtlich der Einflüsse und Veränderun-gen analysiert und die ermittelten Abkühlzeiten t12/8 werden für jede Schweißlage bewertet. Anschließend werden die Berechnungen des Eigenspannungszustandes für einzelne Schweißlagen untersucht.