Neural networks in a multiscale approach for concrete

From a macroscopic point of view, failure within concrete structures is characterized by the initiation and propagation of cracks. In the first part of the thesis, a methodology for macroscopic crack growth simulations for concrete structures using a cohesive discrete crack approach based on the extended finite element method is introduced. Particular attention is turned to the investigation of criteria for crack initiation and crack growth. A drawback of the macroscopic simulation is that the real physical phenomena leading to the nonlinear behavior are only modeled phenomenologically. For concrete, the nonlinear behavior is characterized by the initiation of microcracks which coalesce into macroscopic cracks. In order to obtain a higher resolution of this failure zones, a mesoscale model for concrete is developed that models particles, mortar matrix and the interfacial transition zone (ITZ) explicitly. The essential features are a representation of particles using a prescribed grading curve, a material formulation based on a cohesive approach for the ITZ and a combined model with damage and plasticity for the mortar matrix. Compared to numerical simulations, the response of real structures exhibits a stochastic scatter. This is e.g. due to the intrinsic heterogeneities of the structure. For mesoscale models, these intrinsic heterogeneities are simulated by using a random distribution of particles and by a simulation of spatially variable material parameters using random fields. There are two major problems related to numerical simulations on the mesoscale. First of all, the material parameters for the constitutive description of the materials are often difficult to measure directly. In order to estimate material parameters from macroscopic experiments, a parameter identification procedure based on Bayesian neural networks is developed which is universally applicable to any parameter identification problem in numerical simulations based on experimental results. This approach offers information about the most probable set of material parameters based on experimental data and information about the accuracy of the estimate. Consequently, this approach can be used a priori to determine a set of experiments to be carried out in order to fit the parameters of a numerical model to experimental data. The second problem is the computational effort required for mesoscale simulations of a full macroscopic structure. For this purpose, a coupling between mesoscale and macroscale model is developed. Representative mesoscale simulations are used to train a metamodel that is finally used as a constitutive model in a macroscopic simulation. Special focus is placed on the ability of appropriately simulating unloading.Makroskopisch betrachtet kann das Versagen von Beton durch die Entstehung und das Wachstum von Rissen beschrieben werden. Im ersten Teil der Arbeit wird eine Methode zur Simulation der makroskopischen Rissentwicklung von Beton unter Verwendung von kohäsiven diskreten Rissen basierend auf der erweiterten Finiten Elemente Methode vorgestellt. Besondere Bedeutung liegt dabei auf der Untersuchung von Kriterien zur Rissentstehung und zum Risswachstum. Ein Nachteil von makroskopischen Simulationen liegt in der nur phänomenologischen Berücksichtigung der tatsächlichen Vorgänge. Nichtlineares Verhalten von Beton ist durch die Entstehung von Mikrorissen gekennzeichnet, die bei weiterer Belastung zu makroskopischen Rissen zusammenwachsen. Um die Versagenszone realitätsnah abbilden zu können, wurde ein Mesoskalenmodell von Beton entwickelt, welches Zuschläge, Matrix und Übergangszone zwischen beiden Materialien (ITZ) direkt abbildet. Hauptmerkmal sind die Simulation der Zuschläge nach einer Sieblinie, eine kohäsive Materialformulierung der ITZ und ein kombiniertes Model aus Schädigung und Plastizität für das Matrixmaterial. Im Gegensatz zu numerischen Simulationen ist die Systemantwort reeller Strukturen eine unscharfe Größe. Dies liegt u.a. an Heterogenitäten innerhalb der Struktur, die im Rahmen der Arbeit durch eine zufällige Verteilung der Zuschläge und über räumlich variierende Materialparameter unter Verwendung von Zufallsfeldern simuliert werden. Zwei Hauptprobleme sind bei den Mesoskalensimulationen aufgetreten. Einerseits sind Materialparameter auf der Mesoskala oft schwer zu bestimmen. Deswegen wurde eine Methode basierend auf Bayes neuronalen Netzen entwickelt, die eine Parameteridentifikation unter Verwendung von makroskopischen Versuchen erlaubt. Diese Methode ist aber universell anwendbar auf alle Parameteridentifikationsprobleme in numerischen Simulationen basierend auf experimentellen Daten. Der Ansatz liefert sowohl Informationen über den wahrscheinlichsten Parametersatz des Models zur numerischen Simulation eines Experiments als auch eine Einschätzung der Genauigkeit dieses Schätzers. Die Methode kann auch verwendet werden, um a priori einen Satz von Experimenten auszuwählen der notwendig ist, um die Parameter eines numerischen Modells zu bestimmen. Ein zweites Problem ist der numerische Aufwand von Mesoskalensimulationen für makroskopische Strukturen. Aus diesem Grund wurde eine Kopplungsstrategie zwischen Meso- und Makromodell entwickelt, bei dem repräsentative Simulationen auf der Mesoebene verwendet werden, um ein Metamodell zu generieren, welches dann die Materialformulierung in einer makroskopischen Simulation darstellt. Ein Fokus liegt dabei auf der korrekten Abbildung von Entlastungen

Implicit–explicit integration of gradient-enhanced damage models

This material may be downloaded for personal use only. Any other use requires prior permission of the American Society of Civil Engineers. This material may be found at https://ascelibrary.org/doi/10.1061/%28ASCE%29EM.1943-7889.0001608.Quasi-brittle materials exhibit strain softening. Their modeling requires regularized constitutive formulations to avoid instabilities on the material level. A commonly used model is the implicit gradient-enhanced damage model. For complex geometries, it still shows structural instabilities when integrated with classical backward Euler schemes. An alternative is the implicit–explicit (IMPL-EX) integration scheme. It consists of the extrapolation of internal variables followed by an implicit calculation of the solution fields. The solution procedure for the nonlinear gradient-enhanced damage model is thus transformed into a sequence of problems that are algorithmically linear in every time step. Therefore, they require one single Newton–Raphson iteration per time step to converge. This provides both additional robustness and computational acceleration. The introduced extrapolation error is controlled by adaptive time-stepping schemes. This paper introduced and assessed two novel classes of error control schemes that provide further performance improvements. In a three-dimensional compression test for a mesoscale model of concrete, the presented scheme was about 40 times faster than an adaptive backward Euler time integration.The research was supported by the Federal Institute for Materials Research and Testing, Berlin, Germany and by the German Research Foundation (DFG) under project Un224/7-1. Additionally, the research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Programme (FP/2007-2013) / ERC Grant Agreement n. 320815 (ERC Advanced Grant Project "Advanced tools for computational design of engineering materials" COMP-DES-MAT).Peer ReviewedPostprint (author's final draft

DISCRETE CRACK SIMULATION OF CONCRETE USING THE EXTENDED FINITE ELEMENTMETHOD

The extended finite element method (XFEM) offers an elegant tool to model material discontinuities and cracks within a regular mesh, so that the element edges do not necessarily coincide with the discontinuities. This allows the modeling of propagating cracks without the requirement to adapt the mesh incrementally. Using a regular mesh offers the advantage, that simple refinement strategies based on the quadtree data structure can be used to refine the mesh in regions, that require a high mesh density. An additional benefit of the XFEM is, that the transmission of cohesive forces through a crack can be modeled in a straightforward way without introducing additional interface elements. Finally different criteria for the determination of the crack propagation angle are investigated and applied to numerical tests of cracked concrete specimens, which are compared with experimental results

PARAMETER IDENTIFICATION OF MESOSCALE MODELS FROM MACROSCOPIC TESTS USING BAYESIAN NEURAL NETWORKS

In this paper, a parameter identification procedure using Bayesian neural networks is proposed. Based on a training set of numerical simulations, where the material parameters are simulated in a predefined range using Latin Hypercube sampling, a Bayesian neural network, which has been extended to describe the noise of multiple outputs using a full covariance matrix, is trained to approximate the inverse relation from the experiment (displacements, forces etc.) to the material parameters. The method offers not only the possibility to determine the parameters itself, but also the accuracy of the estimate and the correlation between these parameters. As a result, a set of experiments can be designed to calibrate a numerical model

Model bias identification for Bayesian calibration of stochastic digital twins of bridges

Simulation-based digital twins must provide accurate, robust and reliable digital representations of their physical counterparts. Quantifying the uncertainty in their predictions plays, therefore, a key role in making better-informed decisions that impact the actual system. The update of the simulation model based on data must be then carefully implemented. When applied to complex standing structures such as bridges, discrepancies between the computational model and the real system appear as model bias, which hinders the trustworthiness of the digital twin and increases its uncertainty. Classical Bayesian updating approaches aiming to infer the model parameters often fail at compensating for such model bias, leading to overconfident and unreliable predictions. In this paper, two alternative model bias identification approaches are evaluated in the context of their applicability to digital twins of bridges. A modularized version of Kennedy and O'Hagan's approach and another one based on Orthogonal Gaussian Processes are compared with the classical Bayesian inference framework in a set of representative benchmarks. Additionally, two novel extensions are proposed for such models: the inclusion of noise-aware kernels and the introduction of additional variables not present in the computational model through the bias term. The integration of such approaches in the digital twin corrects the predictions, quantifies their uncertainty, estimates noise from unknown physical sources of error and provides further insight into the system by including additional pre-existing information without modifying the computational model.Comment: 31 pages, 21 figures, 5 tables. Submitted for consideration to Applied Stochastic Models in Business and Industr

The Shapes of Cooperatively Rearranging Regions in Glass Forming Liquids

The shapes of cooperatively rearranging regions in glassy liquids change from being compact at low temperatures to fractal or ``stringy'' as the dynamical crossover temperature from activated to collisional transport is approached from below. We present a quantitative microscopic treatment of this change of morphology within the framework of the random first order transition theory of glasses. We predict a correlation of the ratio of the dynamical crossover temperature to the laboratory glass transition temperature, and the heat capacity discontinuity at the glass transition, Delta C_p. The predicted correlation agrees with experimental results for the 21 materials compiled by Novikov and Sokolov.Comment: 9 pages, 6 figure

Measurement of the cosmic ray spectrum above 4×10184{\times}10^{18} eV using inclined events detected with the Pierre Auger Observatory

A measurement of the cosmic-ray spectrum for energies exceeding $4{\times}10^{18}$ eV is presented, which is based on the analysis of showers with zenith angles greater than $60^{\circ}$ detected with the Pierre Auger Observatory between 1 January 2004 and 31 December 2013. The measured spectrum confirms a flux suppression at the highest energies. Above $5.3{\times}10^{18}$ eV, the "ankle", the flux can be described by a power law $E^{-\gamma}$ with index $\gamma=2.70 \pm 0.02 \,\text{(stat)} \pm 0.1\,\text{(sys)}$ followed by a smooth suppression region. For the energy ($E_\text{s}$) at which the spectral flux has fallen to one-half of its extrapolated value in the absence of suppression, we find $E_\text{s}=(5.12\pm0.25\,\text{(stat)}^{+1.0}_{-1.2}\,\text{(sys)}){\times}10^{19}$ eV.Comment: Replaced with published version. Added journal reference and DO

Energy Estimation of Cosmic Rays with the Engineering Radio Array of the Pierre Auger Observatory

The Auger Engineering Radio Array (AERA) is part of the Pierre Auger Observatory and is used to detect the radio emission of cosmic-ray air showers. These observations are compared to the data of the surface detector stations of the Observatory, which provide well-calibrated information on the cosmic-ray energies and arrival directions. The response of the radio stations in the 30 to 80 MHz regime has been thoroughly calibrated to enable the reconstruction of the incoming electric field. For the latter, the energy deposit per area is determined from the radio pulses at each observer position and is interpolated using a two-dimensional function that takes into account signal asymmetries due to interference between the geomagnetic and charge-excess emission components. The spatial integral over the signal distribution gives a direct measurement of the energy transferred from the primary cosmic ray into radio emission in the AERA frequency range. We measure 15.8 MeV of radiation energy for a 1 EeV air shower arriving perpendicularly to the geomagnetic field. This radiation energy -- corrected for geometrical effects -- is used as a cosmic-ray energy estimator. Performing an absolute energy calibration against the surface-detector information, we observe that this radio-energy estimator scales quadratically with the cosmic-ray energy as expected for coherent emission. We find an energy resolution of the radio reconstruction of 22% for the data set and 17% for a high-quality subset containing only events with at least five radio stations with signal.Comment: Replaced with published version. Added journal reference and DO

Measurement of the Radiation Energy in the Radio Signal of Extensive Air Showers as a Universal Estimator of Cosmic-Ray Energy

We measure the energy emitted by extensive air showers in the form of radio emission in the frequency range from 30 to 80 MHz. Exploiting the accurate energy scale of the Pierre Auger Observatory, we obtain a radiation energy of 15.8 \pm 0.7 (stat) \pm 6.7 (sys) MeV for cosmic rays with an energy of 1 EeV arriving perpendicularly to a geomagnetic field of 0.24 G, scaling quadratically with the cosmic-ray energy. A comparison with predictions from state-of-the-art first-principle calculations shows agreement with our measurement. The radiation energy provides direct access to the calorimetric energy in the electromagnetic cascade of extensive air showers. Comparison with our result thus allows the direct calibration of any cosmic-ray radio detector against the well-established energy scale of the Pierre Auger Observatory.Comment: Replaced with published version. Added journal reference and DOI. Supplemental material in the ancillary file