@article{Zhang, author = {Zhang, Yongzheng}, title = {Nonlocal dynamic Kirchhoff plate formulation based on nonlocal operator method}, series = {Engineering with Computers}, volume = {2022}, journal = {Engineering with Computers}, publisher = {Springer}, address = {London}, doi = {10.1007/s00366-021-01587-1}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20220209-45849}, pages = {1 -- 35}, abstract = {In this study, we propose a nonlocal operator method (NOM) for the dynamic analysis of (thin) Kirchhoff plates. The nonlocal Hessian operator is derived based on a second-order Taylor series expansion. The NOM does not require any shape functions and associated derivatives as 'classical' approaches such as FEM, drastically facilitating the implementation. Furthermore, NOM is higher order continuous, which is exploited for thin plate analysis that requires C1 continuity. The nonlocal dynamic governing formulation and operator energy functional for Kirchhoff plates are derived from a variational principle. The Verlet-velocity algorithm is used for the time discretization. After confirming the accuracy of the nonlocal Hessian operator, several numerical examples are simulated by the nonlocal dynamic Kirchhoff plate formulation.}, subject = {Angewandte Mathematik}, language = {en} } @article{RabczukGuoZhuangetal., author = {Rabczuk, Timon and Guo, Hongwei and Zhuang, Xiaoying and Chen, Pengwan and Alajlan, Naif}, title = {Stochastic deep collocation method based on neural architecture search and transfer learning for heterogeneous porous media}, series = {Engineering with Computers}, volume = {2022}, journal = {Engineering with Computers}, publisher = {Springer}, address = {London}, doi = {10.1007/s00366-021-01586-2}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20220209-45835}, pages = {1 -- 26}, abstract = {We present a stochastic deep collocation method (DCM) based on neural architecture search (NAS) and transfer learning for heterogeneous porous media. We first carry out a sensitivity analysis to determine the key hyper-parameters of the network to reduce the search space and subsequently employ hyper-parameter optimization to finally obtain the parameter values. The presented NAS based DCM also saves the weights and biases of the most favorable architectures, which is then used in the fine-tuning process. We also employ transfer learning techniques to drastically reduce the computational cost. The presented DCM is then applied to the stochastic analysis of heterogeneous porous material. Therefore, a three dimensional stochastic flow model is built providing a benchmark to the simulation of groundwater flow in highly heterogeneous aquifers. The performance of the presented NAS based DCM is verified in different dimensions using the method of manufactured solutions. We show that it significantly outperforms finite difference methods in both accuracy and computational cost.}, subject = {Maschinelles Lernen}, language = {en} } @phdthesis{Nouri, author = {Nouri, Hamidreza}, title = {Mechanical Behavior of two dimensional sheets and polymer compounds based on molecular dynamics and continuum mechanics approach}, doi = {10.25643/bauhaus-universitaet.4670}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20220713-46700}, school = {Bauhaus-Universit{\"a}t Weimar}, pages = {152}, abstract = {Compactly, this thesis encompasses two major parts to examine mechanical responses of polymer compounds and two dimensional materials: 1- Molecular dynamics approach is investigated to study transverse impact behavior of polymers, polymer compounds and two dimensional materials. 2- Large deflection of circular and rectangular membranes is examined by employing continuum mechanics approach. Two dimensional materials (2D), including, Graphene and molybdenum disulfide (MoS2), exhibited new and promising physical and chemical properties, opening new opportunities to be utilized alone or to enhance the performance of conventional materials. These 2D materials have attracted tremendous attention owing to their outstanding physical properties, especially concerning transverse impact loading. Polymers, with the backbone of carbon (organic polymers) or do not include carbon atoms in the backbone (inorganic polymers) like polydimethylsiloxane (PDMS), have extraordinary characteristics particularly their flexibility leads to various easy ways of forming and casting. These simple shape processing label polymers as an excellent material often used as a matrix in composites (polymer compounds). In this PhD work, Classical Molecular Dynamics (MD) is implemented to calculate transverse impact loading of 2D materials as well as polymer compounds reinforced with graphene sheets. In particular, MD was adopted to investigate perforation of the target and impact resistance force . By employing MD approach, the minimum velocity of the projectile that could create perforation and passes through the target is obtained. The largest investigation was focused on how graphene could enhance the impact properties of the compound. Also the purpose of this work was to discover the effect of the atomic arrangement of 2D materials on the impact problem. To this aim, the impact properties of two different 2D materials, graphene and MoS2, are studied. The simulation of chemical functionalization was carried out systematically, either with covalently bonded molecules or with non-bonded ones, focusing the following efforts on the covalently bounded species, revealed as the most efficient linkers. To study transverse impact behavior by using classical MD approach , Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) software, that is well-known among most researchers, is employed. The simulation is done through predefined commands in LAMMPS. Generally these commands (atom style, pair style, angle style, dihedral style, improper style, kspace style, read data, fix, run, compute and so on) are used to simulate and run the model for the desired outputs. Depends on the particles and model types, suitable inter-atomic potentials (force fields) are considered. The ensembles, constraints and boundary conditions are applied depends upon the problem definition. To do so, atomic creation is needed. Python codes are developed to generate particles which explain atomic arrangement of each model. Each atomic arrangement introduced separately to LAMMPS for simulation. After applying constraints and boundary conditions, LAMMPS also include integrators like velocity-Verlet integrator or Brownian dynamics or other types of integrator to run the simulation and finally the outputs are emerged. The outputs are inspected carefully to appreciate the natural behavior of the problem. Appreciation of natural properties of the materials assist us to design new applicable materials. In investigation on the large deflection of circular and rectangular membranes, which is related to the second part of this thesis, continuum mechanics approach is implemented. Nonlinear F{\"o}ppl membrane theory, which carefully release nonlinear governing equations of motion, is considered to establish the non-linear partial differential equilibrium equations of the membranes under distributed and centric point loads. The Galerkin and energy methods are utilized to solve non-linear partial differential equilibrium equations of circular and rectangular plates respectively. Maximum deflection as well as stress through the film region, which are kinds of issue in many industrial applications, are obtained.}, subject = {Molekulardynamik}, language = {en} } @phdthesis{Jenabidehkordi, author = {Jenabidehkordi, Ali}, title = {An Efficient Adaptive PD Formulation for Complex Microstructures}, doi = {10.25643/bauhaus-universitaet.4742}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20221124-47422}, school = {Bauhaus-Universit{\"a}t Weimar}, pages = {118}, abstract = {The computational costs of newly developed numerical simulation play a critical role in their acceptance within both academic use and industrial employment. Normally, the refinement of a method in the area of interest reduces the computational cost. This is unfortunately not true for most nonlocal simulation, since refinement typically increases the size of the material point neighborhood. Reducing the discretization size while keep- ing the neighborhood size will often require extra consideration. Peridy- namic (PD) is a newly developed numerical method with nonlocal nature. Its straightforward integral form equation of motion allows simulating dy- namic problems without any extra consideration required. The formation of crack and its propagation is known as natural to peridynamic. This means that discontinuity is a result of the simulation and does not demand any post-processing. As with other nonlocal methods, PD is considered an expensive method. The refinement of the nodal spacing while keeping the neighborhood size (i.e., horizon radius) constant, emerges to several nonphysical phenomena. This research aims to reduce the peridynamic computational and imple- mentation costs. A novel refinement approach is introduced. The pro- posed approach takes advantage of the PD flexibility in choosing the shape of the horizon by introducing multiple domains (with no intersections) to the nodes of the refinement zone. It will be shown that no ghost forces will be created when changing the horizon sizes in both subdomains. The approach is applied to both bond-based and state-based peridynamic and verified for a simple wave propagation refinement problem illustrating the efficiency of the method. Further development of the method for higher dimensions proves to have a direct relationship with the mesh sensitivity of the PD. A method for solving the mesh sensitivity of the PD is intro- duced. The application of the method will be examined by solving a crack propagation problem similar to those reported in the literature. New software architecture is proposed considering both academic and in- dustrial use. The available simulation tools for employing PD will be collected, and their advantages and drawbacks will be addressed. The challenges of implementing any node base nonlocal methods while max- imizing the software flexibility to further development and modification will be discussed and addressed. A software named Relation-Based Sim- ulator (RBS) is developed for examining the proposed architecture. The exceptional capabilities of RBS will be explored by simulating three dis- tinguished models. RBS is available publicly and open to further develop- ment. The industrial acceptance of the RBS will be tested by targeting its performance on one Mac and two Linux distributions.}, subject = {Peridynamik}, language = {en} } @phdthesis{Vogler, author = {Vogler, Verena}, title = {A framework for artificial coral reef design: Integrating computational modelling and high precision monitoring strategies for artificial coral reefs - an Ecosystem-aware design approach in times of climate change}, isbn = {978-3-00-074495-2}, doi = {10.25643/bauhaus-universitaet.4611}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20220322-46115}, school = {Bauhaus-Universit{\"a}t Weimar}, pages = {243}, abstract = {Tropical coral reefs, one of the world's oldest ecosystems which support some of the highest levels of biodiversity on the planet, are currently facing an unprecedented ecological crisis during this massive human-activity-induced period of extinction. Hence, tropical reefs symbolically stand for the destructive effects of human activities on nature [4], [5]. Artificial reefs are excellent examples of how architectural design can be combined with ecosystem regeneration [6], [7], [8]. However, to work at the interface between the artificial and the complex and temporal nature of natural systems presents a challenge, i.a. in respect to the B-rep modelling legacy of computational modelling. The presented doctorate investigates strategies on how to apply digital practice to realise what is an essential bulwark to retain reefs in impossibly challenging times. Beyond the main question of integrating computational modelling and high precision monitoring strategies in artificial coral reef design, this doctorate explores techniques, methods, and linking frameworks to support future research and practice in ecology led design contexts. Considering the many existing approaches for artificial coral reefs design, one finds they often fall short in precisely understanding the relationships between architectural and ecological aspects (e.g. how a surface design and material composition can foster coral larvae settlement, or structural three-dimensionality enhance biodiversity) and lack an integrated underwater (UW) monitoring process. Such a process is necessary in order to gather knowledge about the ecosystem and make it available for design, and to learn whether artificial structures contribute to reef regeneration or rather harm the coral reef ecosystem. For the research, empirical experimental methods were applied: Algorithmic coral reef design, high precision UW monitoring, computational modelling and simulation, and validated through parallel real-world physical experimentation - two Artificial Reef Prototypes (ARPs) in Gili Trawangan, Indonesia (2012-today). Multiple discrete methods and sub techniques were developed in seventeen computational experiments and applied in a way in which many are cross valid and integrated in an overall framework that is offered as a significant contribution to the field. Other main contributions include the Ecosystem-aware design approach, Key Performance Indicators (KPIs) for coral reef design, algorithmic design and fabrication of Biorock cathodes, new high precision UW monitoring strategies, long-term real-world constructed experiments, new digital analysis methods and two new front-end web-based tools for reef design and monitoring reefs. The methodological framework is a finding of the research that has many technical components that were tested and combined in this way for the very first time. In summary, the thesis responds to the urgency and relevance in preserving marine species in tropical reefs during this massive extinction period by offering a differentiated approach towards artificial coral reefs - demonstrating the feasibility of digitally designing such 'living architecture' according to multiple context and performance parameters. It also provides an in-depth critical discussion of computational design and architecture in the context of ecosystem regeneration and Planetary Thinking. In that respect, the thesis functions as both theoretical and practical background for computational design, ecology and marine conservation - not only to foster the design of artificial coral reefs technically but also to provide essential criteria and techniques for conceiving them. Keywords: Artificial coral reefs, computational modelling, high precision underwater monitoring, ecology in design.}, subject = {Korallenriff}, language = {en} } @phdthesis{Legatiuk, author = {Legatiuk, Anastasiia}, title = {Discrete potential and function theories on a rectangular lattice and their applications}, doi = {10.25643/bauhaus-universitaet.4865}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20221220-48654}, school = {Bauhaus-Universit{\"a}t Weimar}, abstract = {The growing complexity of modern engineering problems necessitates development of advanced numerical methods. In particular, methods working directly with discrete structures, and thus, representing exactly some important properties of the solution on a lattice and not just approximating the continuous properties, become more and more popular nowadays. Among others, discrete potential theory and discrete function theory provide a variety of methods, which are discrete counterparts of the classical continuous methods for solving boundary value problems. A lot of results related to the discrete potential and function theories have been presented in recent years. However, these results are related to the discrete theories constructed on square lattices, and, thus, limiting their practical applicability and potentially leading to higher computational costs while discretising realistic domains. This thesis presents an extension of the discrete potential theory and discrete function theory to rectangular lattices. As usual in the discrete theories, construction of discrete operators is strongly influenced by a definition of discrete geometric setting. For providing consistent constructions throughout the whole thesis, a detailed discussion on the discrete geometric setting is presented in the beginning. After that, the discrete fundamental solution of the discrete Laplace operator on a rectangular lattice, which is the core of the discrete potential theory, its numerical analysis, and practical calculations are presented. By using the discrete fundamental solution of the discrete Laplace operator on a rectangular lattice, the discrete potential theory is then constructed for interior and exterior settings. Several discrete interior and exterior boundary value problems are then solved. Moreover, discrete transmission problems are introduced and several numerical examples of these problems are discussed. Finally, a discrete fundamental solution of the discrete Cauchy-Riemann operator on a rectangular lattice is constructed, and basics of the discrete function theory on a rectangular lattice are provided. This work indicates that the discrete theories provide solution methods with very good numerical properties to tackle various boundary value problems, as well as transmission problems coupling interior and exterior problems. The results presented in this thesis provide a basis for further development of discrete theories on irregular lattices.}, subject = {Diskrete Funktionentheorie}, language = {en} } @phdthesis{Tatarin, author = {Tatarin, Ren{\´e}}, title = {Charakterisieren struktureller Ver{\"a}nderungen in zementgebundenen Baustoffen durch akustische zerst{\"o}rungsfreie Pr{\"u}fverfahren}, publisher = {Cuvillier Verlag}, address = {G{\"o}ttingen}, isbn = {978-3-7369-7575-0}, doi = {10.25643/bauhaus-universitaet.4592}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20220215-45920}, school = {Bauhaus-Universit{\"a}t Weimar}, pages = {293}, abstract = {Im Rahmen dieser Arbeit wird das Charakterisieren struktureller Ver{\"a}nderungen zementgebundener Baustoffe durch zwei auf dem Ultraschall-Transmissionsverfahren beruhenden Methoden der zerst{\"o}rungsfreien Pr{\"u}fung (ZfP) mit mechanischen Wellen vorgenommen. Zur kontinuierlichen Charakterisierung der Erstarrung und Erh{\"a}rtung frischer zementgebundener Systeme wird ein auf Ultraschallsensoren f{\"u}r Longitudinal- und Scherwellen basierendes Messsystem in Kombination mit zugeh{\"o}rigen Verfahrensweisen zur Datenauswertung konzipiert, charakterisiert und angewandt. Gegen{\"u}ber der bislang {\"u}blichen alleinigen Bewertung der Verfestigung anhand indirekter Ultraschallparameter wie Ausbreitungsgeschwindigkeit, Signalenergie oder Frequenzgehalt der Longitudinalwelle l{\"a}sst sich damit eine direkte, sensible Erfassung der sich w{\"a}hrend der Strukturbildung entwickelnden dynamischen elastischen Eigenschaften auf der Basis prim{\"a}rer physikalischer Werkstoffparameter erreichen. Insbesondere Scherwellen und der dynamische Schubmodul sind geeignet, den graduellen {\"U}bergang zum Festk{\"o}rper mit {\"U}berschreiten der Perkolationsschwelle sensibel und unabh{\"a}ngig vom Luftgehalt zu erfassen. Die zeitliche Entwicklung der dynamischen elastischen Eigenschaften, die Strukturbildungsraten sowie die daraus extrahierten diskreten Ergebnisparameter erm{\"o}glichen eine vergleichende quantitative Charakterisierung der Strukturbildung zementgebundener Baustoffe aus mechanischer Sicht. Dabei lassen sich typische, oft unvermeidbare Unterschiede in der Zusammensetzung der Versuchsmischungen ber{\"u}cksichtigen. Der Einsatz laserbasierter Methoden zur Anregung und Erfassung von mechanischen Wellen und deren Kombination zu Laser-Ultraschall zielt darauf ab, die mit der Anwendung des konventionellen Ultraschall-Transmissionsverfahrens verbundenen Nachteile zu eliminieren. Diese resultieren aus der Sensorgeometrie, der mechanischen Ankopplung und bei einer Vielzahl von Oberfl{\"a}chenpunkten aus einem hohen pr{\"u}ftechnischen Aufwand. Die laserbasierte, interferometrische Erfassung mechanischer Wellen ist gegen{\"u}ber Ultraschallsensoren rauschbehaftet und vergleichsweise unsensibel. Als wesentliche Voraussetzung der scannenden Anwendung von Laser-Ultraschall auf zementgebundene Baustoffe erfolgen systematische experimentelle Untersuchungen zur laserinduzierten ablativen Anregung. Diese sollen zum Verst{\"a}ndnis des Anregungsmechanismus unmittelbar auf den Oberfl{\"a}chen von zementgebundenen Baustoffen, Gesteinsk{\"o}rnungen und metallischen Werkstoffen beitragen, relevante Einflussfaktoren aus den charakteristischen Materialeigenschaften identifizieren, geeignete Prozessparameter gewinnen und die Verfahrensgrenzen aufzeigen. Unter Einsatz von Longitudinalwellen erfolgt die Anwendung von Laser-Ultraschall zur zeit- und ortsaufgel{\"o}sten Charakterisierung der Strukturbildung und Homogenit{\"a}t frischer sowie erh{\"a}rteter Proben zementgebundener Baustoffe. W{\"a}hrend der Strukturbildung wird erstmals eine simultane ber{\"u}hrungslose Erfassung von Longitudinal- und Scherwellen vorgenommen. Unter Anwendung von tomographischen Methoden (2D-Laufzeit¬tomo¬graphie) werden {\"u}berlagerungsfreie Informationen zur r{\"a}umlichen Verteilung struktureller Gef{\"u}gever{\"a}nderungen anhand der longitudinalen Ausbreitungsgeschwindigkeit bzw. des relativen dynamischen Elastizit{\"a}tsmoduls innerhalb von virtuellen Schnittebenen gesch{\"a}digter Probek{\"o}rper gewonnen. Als beton-sch{\"a}digende Mechanismen werden exemplarisch der kombinierte Frost-Tausalz-Angriff sowie die Alkali-Kiesels{\"a}ure-Reaktion (AKR) herangezogen. Die im Rahmen dieser Arbeit entwickelten Verfahren der zerst{\"o}rungsfreien Pr{\"u}fung bieten erweiterte M{\"o}glichkeiten zur Charakterisierung zementgebundener Baustoffe und deren strukturellen Ver{\"a}nderungen und lassen sich zielgerichtet in der Werkstoffentwicklung, bei der Qualit{\"a}tssicherung sowie zur Analyse von Schadensprozessen und -ursachen einsetzen.}, subject = {Beton}, language = {de} } @article{AlYasiriMutasharGuerlebecketal., author = {Al-Yasiri, Zainab Riyadh Shaker and Mutashar, Hayder Majid and G{\"u}rlebeck, Klaus and Lahmer, Tom}, title = {Damage Sensitive Signals for the Assessment of the Conditions of Wind Turbine Rotor Blades Using Electromagnetic Waves}, series = {Infrastructures}, volume = {2022}, journal = {Infrastructures}, number = {Volume 7, Issue 8 (August 2022), article 104}, editor = {Shafiullah, GM}, publisher = {MDPI}, address = {Basel}, doi = {10.3390/infrastructures7080104}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20220831-47093}, pages = {18}, abstract = {One of the most important renewable energy technologies used nowadays are wind power turbines. In this paper, we are interested in identifying the operating status of wind turbines, especially rotor blades, by means of multiphysical models. It is a state-of-the-art technology to test mechanical structures with ultrasonic-based methods. However, due to the density and the required high resolution, the testing is performed with high-frequency waves, which cannot penetrate the structure in depth. Therefore, there is a need to adopt techniques in the fields of multiphysical model-based inversion schemes or data-driven structural health monitoring. Before investing effort in the development of such approaches, further insights and approaches are necessary to make the techniques applicable to structures such as wind power plants (blades). Among the expected developments, further accelerations of the so-called "forward codes" for a more efficient implementation of the wave equation could be envisaged. Here, we employ electromagnetic waves for the early detection of cracks. Because in many practical situations, it is not possible to apply techniques from tomography (characterized by multiple sources and sensor pairs), we focus here on the question of whether the existence of cracks can be determined by using only one source for the sent waves.}, subject = {Windkraftwerk}, language = {en} } @phdthesis{Jenabidehkordi, author = {Jenabidehkordi, Ali}, title = {An efficient adaptive PD formulation for complex microstructures}, doi = {10.25643/bauhaus-universitaet.4738}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20221116-47389}, school = {Bauhaus-Universit{\"a}t Weimar}, pages = {118}, abstract = {The computational costs of newly developed numerical simulation play a critical role in their acceptance within both academic use and industrial employment. Normally, the refinement of a method in the area of interest reduces the computational cost. This is unfortunately not true for most nonlocal simulation, since refinement typically increases the size of the material point neighborhood. Reducing the discretization size while keep- ing the neighborhood size will often require extra consideration. Peridynamic (PD) is a newly developed numerical method with nonlocal nature. Its straightforward integral form equation of motion allows simulating dynamic problems without any extra consideration required. The formation of crack and its propagation is known as natural to peridynamic. This means that discontinuity is a result of the simulation and does not demand any post-processing. As with other nonlocal methods, PD is considered an expensive method. The refinement of the nodal spacing while keeping the neighborhood size (i.e., horizon radius) constant, emerges to several nonphysical phenomena. This research aims to reduce the peridynamic computational and imple- mentation costs. A novel refinement approach is introduced. The pro- posed approach takes advantage of the PD flexibility in choosing the shape of the horizon by introducing multiple domains (with no intersections) to the nodes of the refinement zone. It will be shown that no ghost forces will be created when changing the horizon sizes in both subdomains. The approach is applied to both bond-based and state-based peridynamic and verified for a simple wave propagation refinement problem illustrating the efficiency of the method. Further development of the method for higher dimensions proves to have a direct relationship with the mesh sensitivity of the PD. A method for solving the mesh sensitivity of the PD is intro- duced. The application of the method will be examined by solving a crack propagation problem similar to those reported in the literature. New software architecture is proposed considering both academic and in- dustrial use. The available simulation tools for employing PD will be collected, and their advantages and drawbacks will be addressed. The challenges of implementing any node base nonlocal methods while max- imizing the software flexibility to further development and modification will be discussed and addressed. A software named Relation-Based Sim- ulator (RBS) is developed for examining the proposed architecture. The exceptional capabilities of RBS will be explored by simulating three distinguished models. RBS is available publicly and open to further develop- ment. The industrial acceptance of the RBS will be tested by targeting its performance on one Mac and two Linux distributions.}, subject = {Peridynamik}, language = {en} } @phdthesis{Partschefeld, author = {Partschefeld, Stephan}, title = {Synthese von Fließmitteln aus St{\"a}rke und Untersuchung der Wechselwirkung mit Portlandzement}, doi = {10.25643/bauhaus-universitaet.4640}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20220505-46402}, school = {Bauhaus-Universit{\"a}t Weimar}, pages = {145}, abstract = {Das Ziel dieser Arbeit war es, neuartige Fließmittel auf Basis von St{\"a}rke als nachwachsenden Rohstoff zu synthetisieren und die Wechselwirkung mit Portlandzement zu charakterisieren. Die Notwendigkeit, Alternativen zu synthetischen Zusatzmittel zu erforschen, ergibt sich aus der ben{\"o}tigten Menge zur Verarbeitung von ca. 4,1 Gt/a, wobei ca. 85 \% der Zusatzmittel auf die Fließmittel entfallen. Um Fließmittel aus St{\"a}rke zu synthetisieren, wurden drei Basisst{\"a}rken unterschiedlicher Herkunft verwendet. Es wurde eine Maniokst{\"a}rke mit einer niedrigen Molekularmasse und eine Weizenst{\"a}rke mit einer hohen Molekularmasse verwendet. Dar{\"u}ber hinaus wurde eine Kartoffelst{\"a}rke mit einer mittleren Molekularmasse, die ein Abfallprodukt der kartoffelverarbeitenden Industrie darstellt, genutzt. Die St{\"a}rkefließmittel wurden durch chemische Modifikation in einem zweistufigen Prozess synthetisiert. Im ersten Schritt wurde das Molekulargewicht der Weizen- und Kartoffelst{\"a}rke durch s{\"a}urehydrolytischen Abbau verringert. F{\"u}r die kurzkettige Maniokst{\"a}rke war eine Degradation der Molekularmasse nicht notwendig. Im zweiten Syntheseschritt wurden anionische Ladungen durch das Versetzen der degradierten St{\"a}rken und Maniokst{\"a}rke mit Natriumvinylsulfonat in die St{\"a}rkemolek{\"u}le eingef{\"u}hrt. Beurteilung der Synthesemethode zur Erzeugung von St{\"a}rkefließmitteln In diesem Zusammenhang sollten molekulare Parameter der St{\"a}rkefließmittel gezielt eingestellt werden, um eine Fließwirkung im Portlandzement zu erhalten. Insbesondere die Molekularmasse und die Menge anionischer Ladungen sollte variiert werden, um Abh{\"a}ngigkeiten mit der Dispergierleistung zu identifizieren. 1. Es konnte durch GPC-Messungen gezeigt werden, dass die Molekularmasse der langkettigen Weizenst{\"a}rke durch die gew{\"a}hlten Modifizierungsbedingungen zum s{\"a}urehydrolytischen Abbau verringert werden konnte. Durch Variation der s{\"a}urehydrolytischen Bedingungen wurden 4 degradierte Weizenst{\"a}rken erzeugt, die eine Reduzierung der Molekularmasse um 27,5 - 43 \% aufwiesen. Die Molekularmasse der Kartoffelst{\"a}rke konnte durch s{\"a}urehydrolytischen Abbau um ca. 26 \% verringert werden. 2. Durch PCD-Messungen wurde gezeigt, dass anionische Ladungen durch Sulfoethylierung der freien Hydroxylgruppen in die degradierten St{\"a}rken eingef{\"u}hrt werden konnten. Durch Variation der Dauer der Sulfoethylierung konnte die Menge der anionischen Ladungen gesteuert und gezielt variiert werden, so dass St{\"a}rkefließmittel mit steigender Ladungsmenge in folgender Reihenfolge synthetisiert wurden: W-3 < W-2 < K-1 < W¬-4 < W¬1 < M-1 Im Ergebnis der chemischen Modifizierung konnten 6 St{\"a}rkefließmittel mit variierten Molekularmassen und anionischen Ladungen erzeugt werden. Es konnte gezeigt werden, dass die Herkunft der St{\"a}rke f{\"u}r die chemische Modifizierung unerheblich ist. Die Fließmittel lagen synthesebedingt als basische, w{\"a}ssrige Suspensionen mit Wirkstoffgehalten im Bereich von 23,5 - 50 \% vor. Beurteilung der Dispergierleistung der synthetisierten St{\"a}rkefließmittel Die Dispergierperformance wurde durch rheologische Experimente mit einem Rotationsviskosimeter erfasst. Dabei wurden der Einfluss auf die Fließkurven und die Viskosit{\"a}tskurven betrachtet. Durch Vergleich der Dispergierleistung mit einem Polykondensat- und einem PCE-Fließmittel konnte eine Einordnung und Bewertung der Fließmittel vorgenommen werden. 3. Die rheologische Experimente haben gezeigt, dass die St{\"a}rkefließmittel eine vergleichbar hohe Dispergierleistung aufweisen, wie das zum Vergleich herangezogen PCE-Fließmittel. Dar{\"u}ber hinaus zeigte sich, dass die Fließwirkung der 6 St{\"a}rkefließmittel gegen{\"u}ber dem Polykondensatfließmittel deutlich h{\"o}her ist. Das aus der Literatur bekannte Einbrechen der Dispergierleistung der Polykondensat-fließmittel bei w/z-Werten < 0,4 konnte best{\"a}tigt werden. 4. Alle 6 St{\"a}rkefließmittel f{\"u}hrten zu einer Verringerung der Fließgrenze und der dynamischen Viskosit{\"a}t des Zementleimes bei einem w/z-Wert von 0,35. 5. Der Vergleich der Dispergierleistung der St{\"a}rkefließmittel untereinander zeigte, dass die anionische Ladungsmenge einen Schl{\"u}sselparameter darstellt. Die St{\"a}rkefließmittel M-1, K-1, W-1 und W-4 mit anionischen Ladungsmengen > 6 C/g zeigten die h{\"o}chste Dispergier¬performance. Die vergleichend herangezogenen klassischen Fließmittel wiesen anionische Ladungsmengen im Bereich von 1,2 C/g (Polycondensat) und 1,6 C/g (PCE) auf. Die Molekularmasse schien f{\"u}r die Dispergierleistung zun{\"a}chst unerheblich zu sein. Aus diesem Grund wurde die Basisweizenst{\"a}rke erneut chemisch modifiziert, indem anionische Ladungen eingef{\"u}hrt wurden, ohne die Molekularmasse jedoch zu verringern. Das St{\"a}rkederivat wies verdickende Eigenschaften im Zementleim auf. Daraus konnte geschlussfolgert werden, dass eine definierte Grenzmolekularmasse (150.000 Da) existiert, die unterschritten werden muss, um Fließmittel aus St{\"a}rke zu erzeugen. Des Weiteren zeigen die Ergebnisse, dass durch die chemische Modifizierung sowohl Fließmittel als auch Verdicker aus St{\"a}rke erzeugt werden k{\"o}nnen. Beurteilung der Beeinflussung der Hydratation und der Porenl{\"o}sung des Portlandzementes Aus der Literatur ist bekannt, dass Fließmittel die Hydratation von Portlandzement maßgeblich beeinflussen k{\"o}nnen. Aus diesem Grund wurden kalorimetrische und konduktometrische Untersuchungen an Zementsuspensionen, die mit den synthetisierten St{\"a}rkefließmitteln versetzt wurden, durchgef{\"u}hrt. Erg{\"a}nzt wurden die Untersuchungen durch Porenl{\"o}sungsanalysen zu verschiedenen Zeitpunkten der Hydratation. 6. Die kalorimetrischen Untersuchungen zur Beeinflussung der Hydratation des Portlandzementes zeigten, dass die dormante Periode durch die Zugabe der St{\"a}rkefließmittel z.T. erheblich verl{\"a}ngert wird. Es konnte gezeigt werden, dass, je h{\"o}her die anionische Ladungsmenge der St{\"a}rkefließmittel ist, desto l{\"a}nger dauert die dormante Periode andauert. Dar{\"u}ber hinaus zeigte sich, dass eine niedrige Molekularmasse der St{\"a}rkefließmittel die Verl{\"a}ngerung der dormanten Periode beg{\"u}nstigt. 7. Durch die konduktometrischen Untersuchungen konnte gezeigt werden, dass alle St{\"a}rkefließmittel die Dauer des freien- und diffusionskontrollierten CSH-Phasenwachstums verlangsamen. Insbesondere die Ausf{\"a}llung des Portlandits, welches mit dem Erstarrungsbeginn korreliert, erfolgt zu deutlich sp{\"a}teren Zeitpunkten. Des Weiteren korrelierten die konduktometrischen Untersuchungen mit der zeitlichen Entwicklung der Calciumkonzentration der Porenl{\"o}sungen. Der Vergleich der St{\"a}rkefließmittel untereinander zeigte, dass die Molekularmasse ein Schl{\"u}sselparameter ist. Das St{\"a}rkefließmittel M-1 mit der geringsten Molekularmasse, welches geringe Mengen kurzkettiger Anhydroglucoseeinheiten aufweist, verz{\"o}gert die Hydratphasenbildung am st{\"a}rksten. Diese Wirkung ist vergleichbar mit der von Zuckern. Dar{\"u}ber hinaus deuteten die Ergebnisse daraufhin, dass die St{\"a}rkefließmittel auf den ersten Hydratationsprodukten adsorbieren, wodurch die Hydratphasenbildung verlangsamt wird. Die kalorimetrischen und konduktometrischen Daten sowie die Ergebnisse der Porenl{\"o}sungsanalytik des Zementes, erforderten eine genauere Betrachtung der Beeinflussung der Hydratation der Klinkerphasen C3A und C3S, durch die St{\"a}rkefließmittel. Demzufolge wurden die Untersuchungen mit den Klinkerphasen C3A und C3S in Analogie zum Portlandzement durchgef{\"u}hrt. Beurteilung der Beeinflussung der Hydratation und der Porenl{\"o}sung des C3A W{\"a}hrend die kalorimetrischen Untersuchungen zur C3A-Hydratation eine Tendenz zur verlangsamten Hydratphasenbildung durch die St{\"a}rkefließmittel aufzeigten, lieferten die konduktometrischen Ergebnisse grundlegende Erkenntnisse zur Beeinflussung der C3A-Hydratation. Das Stadium I der C3A-Hydratation ist durch einen Abfall der elektrischen Leitf{\"a}higkeit gepr{\"a}gt. Dies korreliert mit dem Absinken der Calciumionenkonzentration und dem Anstieg der Aluminiumionenkonzentration in der Porenl{\"o}sung der C3A-Suspensionen. Im Anschluss an das Stadium I bildet sich ein Plateau in den elektrischen Leitf{\"a}higkeitskurven aus. 8. Es konnte gezeigt werden, dass die St{\"a}rkefließmittel das Stadium I der C3A-Hydratation, d.h. die Aufl{\"o}sung und Bildung erster Calciumaluminathydrate verlangsamen. Insbesondere die St{\"a}rkefließmittel mit h{\"o}herer Molekularmasse erh{\"o}hten die Dauer des Stadium I. Das Stadium II wird durch die St{\"a}rkefließmittel in folgender Reihenfolge am st{\"a}rksten verl{\"a}ngert: M-1 > W-3 > K-1 > W-2 ≥ W-4 und verdeutlicht, dass keine Abh{\"a}ngigkeit von der anionischen Ladungsmenge identifiziert werden konnte. Die Ergebnisse zeigten, dass speziell die kurzkettige St{\"a}rke M-1, das Stadium II l{\"a}nger aufrechterhalten. 9. Das Stadium III und IV der C3A-Hydratation wird insbesondere durch die St{\"a}rkefließmittel mit h{\"o}herer Molekularmasse verl{\"a}ngert. Die Ergebnisse der Porenl{\"o}sungsanalytik korrelieren mit den Ergebnissen der elektrischen Leitf{\"a}higkeit. Speziell die zeitlichen Verl{\"a}ufe der Calciumionenkonzentration bildeten die Verl{\"a}ufe der Konduktivit{\"a}tskurven der C3A-Hydratation mit großer {\"U}bereinstimmung ab. Beurteilung der Beeinflussung der Hydratation und der Porenl{\"o}sung des C3S Die Ergebnisse der kalorimetrischen Untersuchungen zur Beeinflussung der C3S-Hydratation durch die St{\"a}rkefließmittel zeigen, dass diese maßgeblich verlangsamt wird. Das Maximum des Haupthydratationspeaks wird zu sp{\"a}teren Zeiten verschoben und auch die H{\"o}he des Maximums wird deutlich verringert. Durch die konduktometrischen Experimente wurde aufgekl{\"a}rt, welche Stadien der C3S-Hydrataion beeinflusst wurden. 10. Es konnte gezeigt werden, dass sowohl die Menge der eingebrachten anionischen Ladungen als auch das Vorhandensein sehr kleiner St{\"a}rkefließmittelmolek{\"u}le (Zucker), Schl{\"u}sselparameter der verz{\"o}gerten Hydratationskinetik des C3S sind. Der grundlegende Mechanismus der Hydratationsverz{\"o}gerung beruht auf einer Kombination aus verminderter CSH-Keimbildung und Adsorptionsprozessen auf den ersten gebildeten CSH-Phasen der C3S-Partikel. Beurteilung des Adsorptionsverhaltens am Zement, C3A und C3S Die Bestimmung des Adsorptionsverhaltens der St{\"a}rkefließmittel erfolgte mit der Phenol-Schwefels{\"a}ure-Methode an Zement,- C3A- und C3S-Suspensionen. Durch den Vergleich der Adsorptionsraten und Adsorptionsmengen in Abh{\"a}ngigkeit von den molekularen Parametern der St{\"a}rkefließmittel wurde ein Wechselwirkungsmodell identifiziert. 11. Die Ursache f{\"u}r die hohe Dispergierleistung der St{\"a}rkefließmittel liegt in Adsorptionsprozessen an den ersten gebildeten Hydratphasen des Zementes begr{\"u}ndet. Die Molekularmasse der St{\"a}rkefließmittel ist ein Schl{\"u}sselparameter der entscheidend f{\"u}r den Mechanismus der Adsorption ist. W{\"a}hrend anionische, langkettige St{\"a}rken an mehreren Zementpartikeln gleichzeitig adsorbieren und f{\"u}r eine Vernetzung der Zementpartikel untereinander sorgen (Verdickerwirkung), adsorbieren kurzkettige anionische St{\"a}rken lediglich an den ersten gebildeten Hydratphasen der einzelnen Zementpartikel und f{\"u}hren zu elektrostatischer Abstoßung (Fließmittelwirkung). 12. Es konnte gezeigt werden, dass die St{\"a}rkefließmittel mit geringerem Molekulargewicht bei h{\"o}heren Konzentrationen an den Hydratphasen des Zementes adsorbieren. Die St{\"a}rkefließmittel mit h{\"o}herer Molekularmasse erreichen bei einer Zugabemenge von 0,7 \% ein Plateau. Daraus wird geschlussfolgert, dass die gr{\"o}ßeren Fließmittelmolek{\"u}le einen erh{\"o}hten Platzbedarf erfordern und zur Abs{\"a}ttigung der hydratisierenden Oberfl{\"a}chen bei geringeren Zugabemengen f{\"u}hren. Dar{\"u}ber hinaus konnte gezeigt werden, dass die St{\"a}rkefließmittel mit h{\"o}herer anionischer Ladungsmenge zu h{\"o}heren Adsorptionsmengen auf den Zement-, C3A- und C3S-Partikeln f{\"u}hren. 13. Die Adsorptionsprozesse finden an den ersten gebildeten Hydratphasen der C3A-Partikel statt, wodurch sowohl die Aufl{\"o}sung des C3A als auch die Bildung der Calciumhydroaluminate verlangsamt wird. Dar{\"u}ber hinaus wurde festgestellt, dass die Verlangsamung des freien- und diffusionskontrollierten Hydratphasenwachstums des C3S, durch die Adsorption der St{\"a}rkefließmittel auf den ersten gebildeten CSH-Phasen hervorgerufen wird. Des Weiteren wurde festgestellt, dass sehr kleine zucker{\"a}hnliche Molek{\"u}le in der kurzkettigen Maniokst{\"a}rke in der Lage sind, die Bildung der ersten CSH-Keime zu unterdr{\"u}cken. Dadurch kann die langanhaltende Plateauphase der elektrischen Leitf{\"a}higkeit der C3S-Hydratation erkl{\"a}rt werden. Beurteilung der Porenstruktur- und Festigkeitsausbildung Die Beurteilung der Qualit{\"a}t der Mikrostruktur erfolgte durch die Bestimmung der Rohdichte und der Porenradienverteilung mit Hilfe der Quecksilberhochdruckporosimetrie. Durch das Versetzen der Zementleime mit den St{\"a}rkefließmitteln konnten bei gleichbleibender Verarbeitbarkeit Zementsteinprobek{\"o}rper mit einem um 17,5 \% geringeren w/z-Wert von 0,35 hergestellt werden. Die Absenkung des w/z-Wertes f{\"u}hrt zu einem Anstieg der Rohdichte des Zementsteins. 14. Durch die Zugabe der St{\"a}rkefließmittel und den verringerten w/z-Wert wird die Porenstruktur der Zementsteinproben im Vergleich zum Referenzzementstein verfeinert, da die Gesamtporosit{\"a}t sinkt. Insbesondere der Kapillarporenanteil wird verringert und der Gelporenanteil erh{\"o}ht. Im Unterschied zu den PCE-Fließmitteln f{\"u}hrt die Zugabe der St{\"a}rkefließmittel zu keinem erh{\"o}hten Eintrag von Luftporen. Dies wiederum hat zur Folge, dass bei der Verwendung der St{\"a}rkefließmittel auf Entsch{\"a}umer verzichtet werden kann. 15. Entsprechend der dichteren Zementsteinmatrix wurden f{\"u}r die Zementsteine mit den St{\"a}rkefließmitteln nach 7 d und 28 d, erh{\"o}hte Biegezug- und Druckfestigkeiten ermittelt. Insbesondere die 28 d Druckfestigkeit wurde durch den verringerten w/z-Wert um die Faktoren 3,5 - 6,6 erh{\"o}ht.}, subject = {Bauchemie}, language = {de} } @masterthesis{Hoinkis, type = {Bachelor Thesis}, author = {Hoinkis, Jule Hannah}, title = {Hitze in der Stadt Jena}, doi = {10.25643/bauhaus-universitaet.4632}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20220414-46323}, school = {Bauhaus-Universit{\"a}t Weimar}, abstract = {Die vorliegende Arbeit befasst sich mit den spezifischen Faktoren und Wechselwirkungen des st{\"a}dtischen Klimas und Strategien zur Pr{\"a}vention und Kompensation lokaler Klimaver{\"a}nderungen. Problematische Merkmale des Stadtklimas werden sich infolge des Klimawandels st{\"a}rker auspr{\"a}gen. Insbesondere die Hitzebelastung wird zunehmen und die Lebensbedingungen in der Stadt negativ beeinflussen. Infolge h{\"o}herer Temperaturen in St{\"a}dten und einer h{\"o}heren Temperaturdifferenz zum Umland ver{\"a}ndern sich Windstr{\"o}me und die Wasserbilanz. Es sind Strategien notwendig, um den Schadstoffausstoß, die Fl{\"a}cheninanspruchnahme, die Abfallproduktion und den Wasser-, Energie- und Ressourcenverbrauch zu verringern, um sowohl langfristig den Klimawandel als auch dessen bereits unvermeidbaren Auswirkungen auf St{\"a}dte zu begrenzen. Beispielhaft untersucht die Arbeit das Stadtklima, dessen zuk{\"u}nftige Ver{\"a}nderungen infolge des Klimawandels, bauliche Maßnahmen und Anpassungsstrategien der Stadt Jena. Jena ist die zweitgr{\"o}ßte Stadt im Bundesland Th{\"u}ringen und geh{\"o}rt heute zu den w{\"a}rmsten und trockensten Großst{\"a}dten Deutschlands. Die Ergebnisse der Arbeit werden anschließend anhand eines st{\"a}dtebaulichen Konzepts und Entwurfs angewendet. Das Bachstraßenareal liegt in der Innenstadt, dem am st{\"a}rksten von Hitze betroffenen Stadtteil. Als ehemaliger Hauptstandort des Jenaer Universit{\"a}tsklinikums, soll es zu einem nachhaltigen Wissenschaftscampus der Lebenswissenschaften umgebaut werden, wobei ein Großteil der denkmalgesch{\"u}tzten, ehemaligen Klinikgeb{\"a}ude erhalten bleibt. Der Fokus liegt dabei auf der Umsetzung der zuvor formulierten, nachhaltigen Strategien zur Verbesserung des lokalen Stadtklimas und einer Abschw{\"a}chung der Auswirkungen des Klimawandels auf den besonders stark betroffenen Innenstadtbereich Jenas.}, subject = {Hitze}, language = {de} } @phdthesis{Zhang, author = {Zhang, Yongzheng}, title = {A Nonlocal Operator Method for Quasi-static and Dynamic Fracture Modeling}, doi = {10.25643/bauhaus-universitaet.4732}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20221026-47321}, school = {Bauhaus-Universit{\"a}t Weimar}, abstract = {Material failure can be tackled by so-called nonlocal models, which introduce an intrinsic length scale into the formulation and, in the case of material failure, restore the well-posedness of the underlying boundary value problem or initial boundary value problem. Among nonlocal models, peridynamics (PD) has attracted a lot of attention as it allows the natural transition from continuum to discontinue and thus allows modeling of discrete cracks without the need to describe and track the crack topology, which has been a major obstacle in traditional discrete crack approaches. This is achieved by replacing the divergence of the Cauchy stress tensor through an integral over so-called bond forces, which account for the interaction of particles. A quasi-continuum approach is then used to calibrate the material parameters of the bond forces, i.e., equating the PD energy with the energy of a continuum. One major issue for the application of PD to general complex problems is that they are limited to fairly simple material behavior and pure mechanical problems based on explicit time integration. PD has been extended to other applications but losing simultaneously its simplicity and ease in modeling material failure. Furthermore, conventional PD suffers from instability and hourglass modes that require stabilization. It also requires the use of constant horizon sizes, which drastically reduces its computational efficiency. The latter issue was resolved by the so-called dual-horizon peridynamics (DH-PD) formulation and the introduction of the duality of horizons. Within the nonlocal operator method (NOM), the concept of nonlocality is further extended and can be considered a generalization of DH-PD. Combined with the energy functionals of various physical models, the nonlocal forms based on the dual-support concept can be derived. In addition, the variation of the energy functional allows implicit formulations of the nonlocal theory. While traditional integral equations are formulated in an integral domain, the dual-support approaches are based on dual integral domains. One prominent feature of NOM is its compatibility with variational and weighted residual methods. The NOM yields a direct numerical implementation based on the weighted residual method for many physical problems without the need for shape functions. Only the definition of the energy or boundary value problem is needed to drastically facilitate the implementation. The nonlocal operator plays an equivalent role to the derivatives of the shape functions in meshless methods and finite element methods (FEM). Based on the variational principle, the residual and the tangent stiffness matrix can be obtained with ease by a series of matrix multiplications. In addition, NOM can be used to derive many nonlocal models in strong form. The principal contributions of this dissertation are the implementation and application of NOM, and also the development of approaches for dealing with fractures within the NOM, mostly for dynamic fractures. The primary coverage and results of the dissertation are as follows: -The first/higher-order implicit NOM and explicit NOM, including a detailed description of the implementation, are presented. The NOM is based on so-called support, dual-support, nonlocal operators, and an operate energy functional ensuring stability. The nonlocal operator is a generalization of the conventional differential operators. Combining with the method of weighted residuals and variational principles, NOM establishes the residual and tangent stiffness matrix of operate energy functional through some simple matrix without the need of shape functions as in other classical computational methods such as FEM. NOM only requires the definition of the energy drastically simplifying its implementation. For the sake of conciseness, the implementation in this chapter is focused on linear elastic solids only, though the NOM can handle more complex nonlinear problems. An explicit nonlocal operator method for the dynamic analysis of elasticity solid problems is also presented. The explicit NOM avoids the calculation of the tangent stiffness matrix as in the implicit NOM model. The explicit scheme comprises the Verlet-velocity algorithm. The NOM can be very flexible and efficient for solving partial differential equations (PDEs). It's also quite easy for readers to use the NOM and extend it to solve other complicated physical phenomena described by one or a set of PDEs. Several numerical examples are presented to show the capabilities of this method. -A nonlocal operator method for the dynamic analysis of (thin) Kirchhoff plates is proposed. The nonlocal Hessian operator is derived from a second-order Taylor series expansion. NOM is higher-order continuous, which is exploited for thin plate analysis that requires \$C^1\$ continuity. The nonlocal dynamic governing formulation and operator energy functional for Kirchhoff plates are derived from a variational principle. The Verlet-velocity algorithm is used for time discretization. After confirming the accuracy of the nonlocal Hessian operator, several numerical examples are simulated by the nonlocal dynamic Kirchhoff plate formulation. -A nonlocal fracture modeling is developed and applied to the simulation of quasi-static and dynamic fractures using the NOM. The phase field's nonlocal weak and associated strong forms are derived from a variational principle. The NOM requires only the definition of energy. We present both a nonlocal implicit phase field model and a nonlocal explicit phase field model for fracture; the first approach is better suited for quasi-static fracture problems, while the key application of the latter one is dynamic fracture. To demonstrate the performance of the underlying approach, several benchmark examples for quasi-static and dynamic fracture are solved.}, subject = {Variationsprinzip}, 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{Zacharias, author = {Zacharias, Christin}, title = {Numerical Simulation Models for Thermoelastic Damping Effects}, doi = {10.25643/bauhaus-universitaet.4735}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20221116-47352}, school = {Bauhaus-Universit{\"a}t Weimar}, pages = {191}, abstract = {Finite Element Simulations of dynamically excited structures are mainly influenced by the mass, stiffness, and damping properties of the system, as well as external loads. The prediction quality of dynamic simulations of vibration-sensitive components depends significantly on the use of appropriate damping models. Damping phenomena have a decisive influence on the vibration amplitude and the frequencies of the vibrating structure. However, developing realistic damping models is challenging due to the multiple sources that cause energy dissipation, such as material damping, different types of friction, or various interactions with the environment. This thesis focuses on thermoelastic damping, which is the main cause of material damping in homogeneous materials. The effect is caused by temperature changes due to mechanical strains. In vibrating structures, temperature gradients arise in adjacent tension and compression areas. Depending on the vibration frequency, they result in heat flows, leading to increased entropy and the irreversible transformation of mechanical energy into thermal energy. The central objective of this thesis is the development of efficient simulation methods to incorporate thermoelastic damping in finite element analyses based on modal superposition. The thermoelastic loss factor is derived from the structure's mechanical mode shapes and eigenfrequencies. In subsequent analyses that are performed in the time and frequency domain, it is applied as modal damping. Two approaches are developed to determine the thermoelastic loss in thin-walled plate structures, as well as three-dimensional solid structures. The realistic representation of the dissipation effects is verified by comparing the simulation results with experimentally determined data. Therefore, an experimental setup is developed to measure material damping, excluding other sources of energy dissipation. The three-dimensional solid approach is based on the determination of the generated entropy and therefore the generated heat per vibration cycle, which is a measure for thermoelastic loss in relation to the total strain energy. For thin plate structures, the amount of bending energy in a modal deformation is calculated and summarized in the so-called Modal Bending Factor (MBF). The highest amount of thermoelastic loss occurs in the state of pure bending. Therefore, the MBF enables a quantitative classification of the mode shapes concerning the thermoelastic damping potential. The results of the developed simulations are in good agreement with the experimental results and are appropriate to predict thermoelastic loss factors. Both approaches are based on modal superposition with the advantage of a high computational efficiency. Overall, the modeling of thermoelastic damping represents an important component in a comprehensive damping model, which is necessary to perform realistic simulations of vibration processes.}, subject = {Werkstoffd{\"a}mpfung}, 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} } @phdthesis{Valizadeh, author = {Valizadeh, Navid}, title = {Developments in Isogeometric Analysis and Application to High-Order Phase-Field Models of Biomembranes}, doi = {10.25643/bauhaus-universitaet.4565}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20220114-45658}, school = {Bauhaus-Universit{\"a}t Weimar}, abstract = {Isogeometric analysis (IGA) is a numerical method for solving partial differential equations (PDEs), which was introduced with the aim of integrating finite element analysis with computer-aided design systems. The main idea of the method is to use the same spline basis functions which describe the geometry in CAD systems for the approximation of solution fields in the finite element method (FEM). Originally, NURBS which is a standard technology employed in CAD systems was adopted as basis functions in IGA but there were several variants of IGA using other technologies such as T-splines, PHT splines, and subdivision surfaces as basis functions. In general, IGA offers two key advantages over classical FEM: (i) by describing the CAD geometry exactly using smooth, high-order spline functions, the mesh generation process is simplified and the interoperability between CAD and FEM is improved, (ii) IGA can be viewed as a high-order finite element method which offers basis functions with high inter-element continuity and therefore can provide a primal variational formulation of high-order PDEs in a straightforward fashion. The main goal of this thesis is to further advance isogeometric analysis by exploiting these major advantages, namely precise geometric modeling and the use of smooth high-order splines as basis functions, and develop robust computational methods for problems with complex geometry and/or complex multi-physics. As the first contribution of this thesis, we leverage the precise geometric modeling of isogeometric analysis and propose a new method for its coupling with meshfree discretizations. We exploit the strengths of both methods by using IGA to provide a smooth, geometrically-exact surface discretization of the problem domain boundary, while the Reproducing Kernel Particle Method (RKPM) discretization is used to provide the volumetric discretization of the domain interior. The coupling strategy is based upon the higher-order consistency or reproducing conditions that are directly imposed in the physical domain. The resulting coupled method enjoys several favorable features: (i) it preserves the geometric exactness of IGA, (ii) it circumvents the need for global volumetric parameterization of the problem domain, (iii) it achieves arbitrary-order approximation accuracy while preserving higher-order smoothness of the discretization. Several numerical examples are solved to show the optimal convergence properties of the coupled IGA-RKPM formulation, and to demonstrate its effectiveness in constructing volumetric discretizations for complex-geometry objects. As for the next contribution, we exploit the use of smooth, high-order spline basis functions in IGA to solve high-order surface PDEs governing the morphological evolution of vesicles. These governing equations are often consisted of geometric PDEs, high-order PDEs on stationary or evolving surfaces, or a combination of them. We propose an isogeometric formulation for solving these PDEs. In the context of geometric PDEs, we consider phase-field approximations of mean curvature flow and Willmore flow problems and numerically study the convergence behavior of isogeometric analysis for these problems. As a model problem for high-order PDEs on stationary surfaces, we consider the Cahn-Hilliard equation on a sphere, where the surface is modeled using a phase-field approach. As for the high-order PDEs on evolving surfaces, a phase-field model of a deforming multi-component vesicle, which consists of two fourth-order nonlinear PDEs, is solved using the isogeometric analysis in a primal variational framework. Through several numerical examples in 2D, 3D and axisymmetric 3D settings, we show the robustness of IGA for solving the considered phase-field models. Finally, we present a monolithic, implicit formulation based on isogeometric analysis and generalized-alpha time integration for simulating hydrodynamics of vesicles according to a phase-field model. Compared to earlier works, the number of equations of the phase-field model which need to be solved is reduced by leveraging high continuity of NURBS functions, and the algorithm is extended to 3D settings. We use residual-based variational multi-scale method (RBVMS) for solving Navier-Stokes equations, while the rest of PDEs in the phase-field model are treated using a standard Galerkin-based IGA. We introduce the resistive immersed surface (RIS) method into the formulation which can be employed for an implicit description of complex geometries using a diffuse-interface approach. The implementation highlights the robustness of the RBVMS method for Navier-Stokes equations of incompressible flows with non-trivial localized forcing terms including bending and tension forces of the vesicle. The potential of the phase-field model and isogeometric analysis for accurate simulation of a variety of fluid-vesicle interaction problems in 2D and 3D is demonstrated.}, subject = {Phasenfeldmodell}, language = {en} } @phdthesis{ShaabanMohamed, author = {Shaaban Mohamed, Ahmed Mostafa}, title = {Isogeometric boundary element analysis and structural shape optimization for Helmholtz acoustic problems}, doi = {10.25643/bauhaus-universitaet.4703}, url = {http://nbn-resolving.de/urn:nbn:de:gbv:wim2-20220816-47030}, school = {Bauhaus-Universit{\"a}t Weimar}, abstract = {In this thesis, a new approach is developed for applications of shape optimization on the time harmonic wave propagation (Helmholtz equation) for acoustic problems. This approach is introduced for different dimensional problems: 2D, 3D axi-symmetric and fully 3D problems. The boundary element method (BEM) is coupled with the isogeometric analysis (IGA) forming the so-called (IGABEM) which speeds up meshing and gives higher accuracy in comparison with standard BEM. BEM is superior for handling unbounded domains by modeling only the inner boundaries and avoiding the truncation error, present in the finite element method (FEM) since BEM solutions satisfy the Sommerfeld radiation condition automatically. Moreover, BEM reduces the space dimension by one from a volumetric three-dimensional problem to a surface two-dimensional problem, or from a surface two-dimensional problem to a perimeter one-dimensional problem. Non-uniform rational B-splines basis functions (NURBS) are used in an isogeometric setting to describe both the CAD geometries and the physical fields. IGABEM is coupled with one of the gradient-free optimization methods, the Particle Swarm Optimization (PSO) for structural shape optimization problems. PSO is a straightforward method since it does not require any sensitivity analysis but it has some trade-offs with regard to the computational cost. Coupling IGA with optimization problems enables the NURBS basis functions to represent the three models: shape design, analysis and optimization models, by a definition of a set of control points to be the control variables and the optimization parameters as well which enables an easy transition between the three models. Acoustic shape optimization for various frequencies in different mediums is performed with PSO and the results are compared with the benchmark solutions from the literature for different dimensional problems proving the efficiency of the proposed approach with the following remarks: - In 2D problems, two BEM methods are used: the conventional isogeometric boundary element method (IGABEM) and the eXtended IGABEM (XIBEM) enriched with the partition-of-unity expansion using a set of plane waves, where the results are generally in good agreement with the linterature with some computation advantage to XIBEM which allows coarser meshes. -In 3D axi-symmetric problems, the three-dimensional problem is simplified in BEM from a surface integral to a combination of two 1D integrals. The first is the line integral similar to a two-dimensional BEM problem. The second integral is performed over the angle of revolution. The discretization is applied only to the former integration. This leads to significant computational savings and, consequently, better treatment for higher frequencies over the full three-dimensional models. - In fully 3D problems, a detailed comparison between two BEM methods: the conventional boundary integral equation (CBIE) and Burton-Miller (BM) is provided including the computational cost. The proposed models are enhanced with a modified collocation scheme with offsets to Greville abscissae to avoid placing collocation points at the corners. Placing collocation points on smooth surface enables accurate evaluation of normals for BM formulation in addition to straightforward prediction of jump-terms and avoids singularities in \$\mathcal{O} (1/r)\$ integrals eliminating the need for polar integration. Furthermore, no additional special treatment is required for the hyper-singular integral while collocating on highly distorted elements, such as those containing sphere poles. The obtained results indicate that, CBIE with PSO is a feasible alternative (except for a small number of fictitious frequencies) which is easier to implement. Furthermore, BM presents an outstanding treatment of the complicated geometry of mufflers with internal extended inlet/outlet tube as an interior 3D Helmholtz acoustic problem instead of using mixed or dual BEM.}, subject = {Randelemente-Methode}, language = {en} }