@inproceedings{KirichukKoeppler1997, author = {Kirichuk, A. and K{\"o}ppler, H.}, title = {Numerical Algorithms and Computer Modeling for nonlinear Analysis of Shell Structures}, doi = {10.25643/bauhaus-universitaet.438}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20111215-4382}, year = {1997}, abstract = {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.}, subject = {Schale}, language = {en} } @inproceedings{vanRooyenOlivier2004, author = {van Rooyen, G.C. and Olivier, A. H.}, title = {Notes on structural analysis in a distributed collaboratory}, doi = {10.25643/bauhaus-universitaet.145}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20111215-1451}, year = {2004}, abstract = {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.}, subject = {Ingenieurbau}, language = {en} } @inproceedings{EbertBucher2000, author = {Ebert, Matthias and Bucher, Christian}, title = {Modelling of changing of dynamic and static parameters of damaged R/C}, doi = {10.25643/bauhaus-universitaet.582}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20111215-5825}, year = {2000}, abstract = {Dynamic testing for damage assessment as non-destructive method has attracted growing in-terest for systematic inspections and maintenance of civil engineering structures. In this con-text the paper presents the Stochastic Finite Element (SFE) Modeling of the static and dy-namic results of own four point bending experiments with R/C beams. The beams are dam-aged by an increasing load. Between the load levels the dynamic properties are determined. Calculated stiffness loss factors for the displacements and the natural frequencies show differ-ent histories. A FE Model for the beams is developed with a discrete crack formulation. Cor-related random fields are used for structural parameters stiffness and tension strength. The idea is to simulate different crack evolutions. The beams have the same design parameters, but because of the stochastic material properties their undamaged state isn't yet the same. As the structure is loaded a stochastic first crack occurs on the weakest place of the structure. The further crack evolution is also stochastic. These is a great advantage compared with de-terministic formulations. To reduce the computational effort of the Monte Carlo simulation of this nonlinear problem the Latin-Hypercube sampling technique is applied. From the results functions of mean value and standard deviation of displacements and frequencies are calcu-lated. Compared with the experimental results some qualitative phenomena are good de-scribed by the model. Differences occurs especially in the dynamic behavior of the higher load levels. Aim of the investigations is to assess the possibilities of dynamic testing under consideration of effects from stochastic material properties}, subject = {Stahlbetonbauteil}, language = {en} } @article{MortazaviPereiraJiangetal., author = {Mortazavi, Bohayra and Pereira, Luiz Felipe C. and Jiang, Jin-Wu and Rabczuk, Timon}, title = {Modelling heat conduction in polycrystalline hexagonal boron-nitride films}, series = {Scientific Reports}, journal = {Scientific Reports}, doi = {10.1038/srep13228}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20170425-31534}, abstract = {We conducted extensive molecular dynamics simulations to investigate the thermal conductivity of polycrystalline hexagonal boron-nitride (h-BN) films. To this aim, we constructed large atomistic models of polycrystalline h-BN sheets with random and uniform grain configuration. By performing equilibrium molecular dynamics (EMD) simulations, we investigated the influence of the average grain size on the thermal conductivity of polycrystalline h-BN films at various temperatures. Using the EMD results, we constructed finite element models of polycrystalline h-BN sheets to probe the thermal conductivity of samples with larger grain sizes. Our multiscale investigations not only provide a general viewpoint regarding the heat conduction in h-BN films but also propose that polycrystalline h-BN sheets present high thermal conductivity comparable to monocrystalline sheets.}, subject = {W{\"a}rmeleitf{\"a}higkeit}, language = {en} } @phdthesis{Jia, author = {Jia, Yue}, title = {Methods based on B-splines for model representation, numerical analysis and image registration}, doi = {10.25643/bauhaus-universitaet.2484}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20151210-24849}, school = {Bauhaus-Universit{\"a}t Weimar}, pages = {200}, abstract = {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.}, subject = {Finite-Elemente-Methode}, language = {en} } @inproceedings{PopovaDatchevaIankov2003, author = {Popova, E. D. and Datcheva, Maria and Iankov, Roumen}, title = {Mechanical Models with Interval Parameters}, doi = {10.25643/bauhaus-universitaet.348}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20111215-3482}, year = {2003}, abstract = {In this paper we consider modelling of composite material with inclusions where the elastic material properties of both matrix and inclusions are uncertain and vary within prescribed bounds. Such mechanical systems, involving interval uncertainties and modelled by finite element method, can be described by parameter dependent systems of linear interval equations and process variables depending on the system solution. A newly developed hybrid interval approach for solving parametric interval linear systems is applied to the considered model and the results are compared to other interval methods. The hybrid approach provides very sharp bounds for the process variables - element strains and stresses. The sources for overestimation when dealing with interval computations are demonstrated. Based on the element strains and stresses, we introduce a definition for the values of nodal strains and stresses by using a set-theoretic approach.}, subject = {Verbundwerkstoff}, language = {en} } @phdthesis{Kessler2018, author = {Keßler, Andrea}, title = {Matrix-free voxel-based finite element method for materials with heterogeneous microstructures}, doi = {10.25643/bauhaus-universitaet.3844}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20190116-38448}, school = {Bauhaus-Universit{\"a}t Weimar}, pages = {113}, year = {2018}, abstract = {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.}, subject = {Dissertation}, language = {en} } @article{GalffyWellmannJelicHartmann2004, author = {Galffy, Mozes and Wellmann Jelic, Andres and Hartmann, Dietrich}, title = {Lifetime-oriented modelling of vortex-induced across-wind vibrations on bridge tie rods}, doi = {10.25643/bauhaus-universitaet.253}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20111215-2536}, year = {2004}, abstract = {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.}, subject = {Finite-Elemente-Methode}, language = {en} } @article{KashiyamaHamadaTaniguchi2004, author = {Kashiyama, Kazuo and Hamada, Hidetaka and Taniguchi, Takeo}, title = {Large Scale Finite Element Simulation and Modeling Using GIS/CAD for Environmental Flows in Urban Area}, doi = {10.25643/bauhaus-universitaet.267}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20111215-2675}, year = {2004}, abstract = {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.}, subject = {Geoinformationssystem}, language = {en} } @phdthesis{LopezZermeno, author = {L{\´o}pez Zerme{\~n}o, Jorge Alberto}, title = {Isogeometric and CAD-based methods for shape and topology optimization: Sensitivity analysis, B{\´e}zier elements and phase-field approaches}, doi = {10.25643/bauhaus-universitaet.4710}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20220831-47102}, school = {Bauhaus-Universit{\"a}t Weimar}, abstract = {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{\´e}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{\´e}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{\´e}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{\´e}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.}, subject = {CAD}, language = {en} } @phdthesis{Hossain, author = {Hossain, Md Naim}, title = {Isogeometric analysis based on Geometry Independent Field approximaTion (GIFT) and Polynomial Splines over Hierarchical T-meshes}, doi = {10.25643/bauhaus-universitaet.4037}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20191129-40376}, school = {Bauhaus-Universit{\"a}t Weimar}, pages = {157}, abstract = {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.}, subject = {Finite-Elemente-Methode}, language = {en} } @phdthesis{Nanthakumar, author = {Nanthakumar, S.S.}, title = {Inverse and optimization problems in piezoelectric materials using Extended Finite Element Method and Level sets}, doi = {10.25643/bauhaus-universitaet.2709}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20161128-27095}, school = {Bauhaus-Universit{\"a}t Weimar}, abstract = {Piezoelectric materials are used in several applications as sensors and actuators where they experience high stress and electric field concentrations as a result of which they may fail due to fracture. Though there are many analytical and experimental works on piezoelectric fracture mechanics. There are very few studies about damage detection, which is an interesting way to prevent the failure of these ceramics. An iterative method to treat the inverse problem of detecting cracks and voids in piezoelectric structures is proposed. Extended finite element method (XFEM) is employed for solving the inverse problem as it allows the use of a single regular mesh for large number of iterations with different flaw geometries. Firstly, minimization of cost function is performed by Multilevel Coordinate Search (MCS) method. The XFEM-MCS methodology is applied to two dimensional electromechanical problems where flaws considered are straight cracks and elliptical voids. Then a numerical method based on combination of classical shape derivative and level set method for front propagation used in structural optimization is utilized to minimize the cost function. The results obtained show that the XFEM-level set methodology is effectively able to determine the number of voids in a piezoelectric structure and its corresponding locations. The XFEM-level set methodology is improved to solve the inverse problem of detecting inclusion interfaces in a piezoelectric structure. The material interfaces are implicitly represented by level sets which are identified by applying regularisation using total variation penalty terms. The formulation is presented for three dimensional structures and inclusions made of different materials are detected by using multiple level sets. The results obtained prove that the iterative procedure proposed can determine the location and approximate shape of material subdomains in the presence of higher noise levels. Piezoelectric nanostructures exhibit size dependent properties because of surface elasticity and surface piezoelectricity. Initially a study to understand the influence of surface elasticity on optimization of nano elastic beams is performed. The boundary of the nano structure is implicitly represented by a level set function, which is considered as the design variable in the optimization process. Two objective functions, minimizing the total potential energy of a nanostructure subjected to a material volume constraint and minimizing the least square error compared to a target displacement, are chosen for the numerical examples. The numerical examples demonstrate the importance of size and aspect ratio in determining how surface effects impact the optimized topology of nanobeams. Finally a conventional cantilever energy harvester with a piezoelectric nano layer is analysed. The presence of surface piezoelectricity in nano beams and nano plates leads to increase in electromechanical coupling coefficient. Topology optimization of these piezoelectric structures in an energy harvesting device to further increase energy conversion using appropriately modified XFEM-level set algorithm is performed .}, subject = {Finite-Elemente-Methode}, language = {de} } @phdthesis{Schwedler, author = {Schwedler, Michael}, title = {Integrated structural analysis using isogeometric finite element methods}, doi = {10.25643/bauhaus-universitaet.2737}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20170130-27372}, school = {Bauhaus-Universit{\"a}t Weimar}, pages = {209}, abstract = {The gradual digitization in the architecture, engineering, and construction industry over the past fifty years led to an extremely heterogeneous software environment, which today is embodied by the multitude of different digital tools and proprietary data formats used by the many specialists contributing to the design process in a construction project. Though these projects become increasingly complex, the demands on financial efficiency and the completion within a tight schedule grow at the same time. The digital collaboration of project partners has been identified as one key issue in successfully dealing with these challenges. Yet currently, the numerous software applications and their respective individual views on the design process severely impede that collaboration. An approach to establish a unified basis for the digital collaboration, regardless of the existing software heterogeneity, is a comprehensive digital building model contributed to by all projects partners. This type of data management known as building information modeling (BIM) has many benefits, yet its adoption is associated with many difficulties and thus, proceeds only slowly. One aspect in the field of conflicting requirements on such a digital model is the cooperation of architects and structural engineers. Traditionally, these two disciplines use different abstractions of reality for their models that in consequence lead to incompatible digital representations thereof. The onset of isogeometric analysis (IGA) promised to ease the discrepancy in design and analysis model representations. Yet, that initial focus quickly shifted towards using these methods as a more powerful basis for numerical simulations. Furthermore, the isogeometric representation alone is not capable of solving the model abstraction problem. It is thus the intention of this work to contribute to an improved digital collaboration of architects and engineers by exploring an integrated analysis approach on the basis of an unified digital model and solid geometry expressed by splines. In the course of this work, an analysis framework is developed that utilizes such models to automatically conduct numerical simulations commonly required in construction projects. In essence, this allows to retrieve structural analysis results from BIM models in a fast and simple manner, thereby facilitating rapid design iterations and profound design feedback. The BIM implementation Industry Foundation Classes (IFC) is reviewed with regard to its capabilities of representing the unified model. The current IFC schema strongly supports the use of redundant model data, a major pitfall in digital collaboration. Additionally, it does not allow to describe the geometry by volumetric splines. As the pursued approach builds upon a unique model for both, architectural and structural design, and furthermore requires solid geometry, necessary schema modifications are suggested. Structural entities are modeled by volumetric NURBS patches, each of which constitutes an individual subdomain that, with regard to the analysis, is incompatible with the remaining full model. The resulting consequences for numerical simulation are elaborated in this work. The individual subdomains have to be weakly coupled, for which the mortar method is used. Different approaches to discretize the interface traction fields are implemented and their respective impact on the analysis results is evaluated. All necessary coupling conditions are automatically derived from the related geometry model. The weak coupling procedure leads to a linear system of equations in saddle point form, which, owed to the volumetric modeling, is large in size and, the associated coefficient matrix has, due to the use of higher degree basis functions, a high bandwidth. The peculiarities of the system require adapted solution methods that generally cause higher numerical costs than the standard procedures for symmetric, positive-definite systems do. Different methods to solve the specific system are investigated and an efficient parallel algorithm is finally proposed. When the structural analysis model is derived from the unified model in the BIM data, it does in general initially not meet the requirements on the discretization that are necessary to obtain sufficiently accurate analysis results. The consequently necessary patch refinements must be controlled automatically to allowfor an entirely automatic analysis procedure. For that purpose, an empirical refinement scheme based on the geometrical and possibly mechanical properties of the specific entities is proposed. The level of refinement may be selectively manipulated by the structural engineer in charge. Furthermore, a Zienkiewicz-Zhu type error estimator is adapted for the use with isogeometric analysis results. It is shown that also this estimator can be used to steer an adaptive refinement procedure.}, subject = {Finite-Elemente-Methode}, language = {en} } @misc{Kurukuri2005, type = {Master Thesis}, author = {Kurukuri, Srihari}, title = {Homogenization of Damaged Concrete Mesostructures using Representative Volume Elements - Implementation and Application to SLang}, doi = {10.25643/bauhaus-universitaet.667}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20111215-6670}, school = {Bauhaus-Universit{\"a}t Weimar}, year = {2005}, abstract = {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.}, subject = {Schadensmechanik}, language = {en} } @inproceedings{GurtovyPiskunov2000, author = {Gurtovy, O. G. and Piskunov, V. G.}, title = {HIGH-PRECISION MODELING AND FINITE-ELEMENT INVESTIGATION OF ELASTOPLASTIC DEFORMATION OF NON-ISOTROPIC THICK SANDWICH PLATES AND SHELLS}, doi = {10.25643/bauhaus-universitaet.584}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20111215-5848}, year = {2000}, abstract = {There was suggested a phenomenological modified quadratic condition of the beginning of plasticity for plastic and quasifragile orthotropic materials. Limiting surface in the shape of a paraboloid with an axis bend over hydrostatic axis corresponds to the condition. The equations of theory of current with the isotropic and anisotropic hardenings, associated with the suggested yield condition, modified into the version of determining equations of strain theory of plasticity are received. These defining equations formed the basis of highlyprecise non-classic continual (along thickness) theory of non-linear deformation of thick sandwich plates and sloping shells. In the approximations along the cross coordinate the specificity of flexural and non-flexural deformations is taken into account. The necessity of introducing the approximations of higher order, as well as accounting for the cross compression while decreasing of the relatively cross normal and shear layer rigidness is shown. The specifications, obtained in comparison with the known physically nonlinear specified model of the bending of plates with orthotropic layers are distinguished. An effective procedure of linearization of the solving equations and getting the solutions in frames of the discrete-continual scheme of the finite-element method is suggested. The approximations of higher order let to model the appearance of the cracs of layers being split by the introducing of slightly hard thin layers into the finite element, not violating the idea of continuality of theory. Calculation of a threelayer plate with rigid face diaphragms on the contour is considered}, subject = {Platte}, language = {en} } @phdthesis{Haefner2006, author = {H{\"a}fner, Stefan}, title = {Grid-based procedures for the mechanical analysis of heterogeneous solids}, doi = {10.25643/bauhaus-universitaet.858}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20070830-9185}, school = {Bauhaus-Universit{\"a}t Weimar}, year = {2006}, abstract = {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.}, subject = {B-Spline}, language = {en} } @phdthesis{Habtemariam, author = {Habtemariam, Abinet Kifle}, title = {Generalized Beam Theory for the analysis of thin-walled circular pipe members}, doi = {10.25643/bauhaus-universitaet.4572}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20220127-45723}, school = {Bauhaus-Universit{\"a}t Weimar}, pages = {188}, abstract = {The detailed structural analysis of thin-walled circular pipe members often requires the use of a shell or solid-based finite element method. Although these methods provide a very good approximation of the deformations, they require a higher degree of discretization which causes high computational costs. On the other hand, the analysis of thin-walled circular pipe members based on classical beam theories is easy to implement and needs much less computation time, however, they are limited in their ability to approximate the deformations as they cannot consider the deformation of the cross-section. This dissertation focuses on the study of the Generalized Beam Theory (GBT) which is both accurate and efficient in analyzing thin-walled members. This theory is based on the separation of variables in which the displacement field is expressed as a combination of predetermined deformation modes related to the cross-section, and unknown amplitude functions defined on the beam's longitudinal axis. Although the GBT was initially developed for long straight members, through the consideration of complementary deformation modes, which amend the null transverse and shear membrane strain assumptions of the classical GBT, problems involving short members, pipe bends, and geometrical nonlinearity can also be analyzed using GBT. In this dissertation, the GBT formulation for the analysis of these problems is developed and the application and capabilities of the method are illustrated using several numerical examples. Furthermore, the displacement and stress field results of these examples are verified using an equivalent refined shell-based finite element model. The developed static and dynamic GBT formulations for curved thin-walled circular pipes are based on the linear kinematic description of the curved shell theory. In these formulations, the complex problem in pipe bends due to the strong coupling effect of the longitudinal bending, warping and the cross-sectional ovalization is handled precisely through the derivation of the coupling tensors between the considered GBT deformation modes. Similarly, the geometrically nonlinear GBT analysis is formulated for thin-walled circular pipes based on the nonlinear membrane kinematic equations. Here, the initial linear and quadratic stress and displacement tangent stiffness matrices are built using the third and fourth-order GBT deformation mode coupling tensors. Longitudinally, the formulation of the coupled GBT element stiffness and mass matrices are presented using a beam-based finite element formulation. Furthermore, the formulated GBT elements are tested for shear and membrane locking problems and the limitations of the formulations regarding the membrane locking problem are discussed.}, subject = {Finite-Elemente-Methode}, language = {en} } @inproceedings{MilbradtSchwoeppe2003, author = {Milbradt, Peter and Schw{\"o}ppe, Axel}, title = {Finite Element Approximation auf der Basis geometrischer Zellen}, doi = {10.25643/bauhaus-universitaet.333}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20111215-3333}, year = {2003}, abstract = {Die Methode der Finiten Elemente ist ein numerisches Verfahren zur Interpolation vorgegebener Werte und zur numerischen Approximation von L{\"o}sungen station{\"a}rer oder instation{\"a}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{\"u}hrt zu einer geometrieunabh{\"a}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 {\"u}ber das Standard-Galerkin-Verfahren ermittelt.}, subject = {Finite-Elemente-Methode}, language = {de} } @article{MilbradtSchierbaumSchwoeppe2004, author = {Milbradt, Peter and Schierbaum, Jochen and Schw{\"o}ppe, Axel}, title = {Finite Cell-Elements of Higher Order}, doi = {10.25643/bauhaus-universitaet.252}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20111215-2524}, year = {2004}, abstract = {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...}, subject = {Finite-Elemente-Methode}, language = {en} } @inproceedings{BergerGraeffWeinberg1997, author = {Berger, Hans and Graeff-Weinberg, K.}, title = {FEM-Detailuntersuchungen an Tragwerken unter Einsatz von pNh-{\"U}bergangselementen}, doi = {10.25643/bauhaus-universitaet.426}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20111215-4267}, year = {1997}, abstract = {Detailuntersuchungen an Tragwerken f{\"u}hren bei FE-Berechnungen immer wieder auf das Problem einer geeigneten Netzgestaltung. W{\"a}hrend in weiten Bereichen ein grobes Netz ausreicht, muß an kritischen Stellen ein sehr feines Netz gew{\"a}hlt werden, um gerade dort hinreichend genaue Ergebnisse zu erhalten. Bei der Realisierung lokaler Netzverdichtungen stellt die Gestaltung des {\"U}bergangs vom groben zum feinen Netz das Hauptproblem dar. Im Beitrag wird hierzu eine Familie von FE-{\"U}bergangselementen vorgestellt, mit denen sich eine voll-kompatible Kopplung von wenigen großen Elementen mit vielen kleinen Elementen bereits {\"u}ber nur eine Stufe erzielen l{\"a}ßt. Diese neu entwickelten sogenannten pNh-Elemente erm{\"o}glichen an einer oder mehreren Seiten den Anschluß von N kleineren Elementen (Elementseiten f{\"u}r h-Verfeinerung). Das wird durch N st{\"u}ckweise definierte Ansatzfunktionen an den entsprechenden Seiten erreicht, wobei die Teilung nicht {\"a}quidistant sein braucht. Dar{\"u}ber hinaus ist es m{\"o}glich, Elemente unterschiedlichen Polynomgrades p an den Standardseiten und den Verfeinerungsseiten anzuschließen. Der praktische Einsatz der {\"U}bergangselemente setzt geeignete automatische oder halbautomatische Netzgeneratoren voraus, die diese Elemente einbeziehen. Im Rahmen einer substrukturorientierten Modellierung l{\"a}ßt sich dies besonders g{\"u}nstig realisieren. Im Beitrag wird gezeigt, wie durch Zerlegung des Gesamtmodells in Bereiche mit grobem Netz, mit {\"U}bergangsnetz und mit feinem Netz, eine effektive Generierung der Netzverdichtungen zu erreichen ist. An einem praktischen Beispiel aus dem Bauingenieurwesen werden die Vorteile des vorgestellten {\"U}bergangselementkonzeptes umfassend demonstriert.}, subject = {Tragwerk}, language = {de} }