5.5 Materialmodellierung
Filtern
Dokumenttyp
- Vortrag (32)
- Zeitschriftenartikel (29)
- Beitrag zu einem Tagungsband (9)
- Posterpräsentation (3)
- Buchkapitel (2)
- Forschungsdatensatz (2)
- Dissertation (1)
- Forschungsbericht (1)
Sprache
- Englisch (67)
- Deutsch (11)
- Französisch (1)
Schlagworte
- Creep (9)
- Nickel-base superalloys (8)
- Fatigue (7)
- LCF (7)
- Thermomechanical fatigue (7)
- Additive manufacturing (5)
- Anisotropy (5)
- Crack growth (5)
- Crystal plasticity (5)
- HIP (5)
Organisationseinheit der BAM
- 5.2 Metallische Hochtemperaturwerkstoffe (79) (entfernen)
Eingeladener Vortrag
- nein (32)
The current competitive situation on electricity markets forces conventional power plants into cyclic operation regimes with frequent load shifts and starts/shutdowns. In the present work, the cyclic mechanical behavior of ferritic-martensitic 9-12 % Cr steels under isothermal and thermomechanical loading was investigated for the example of grade P92 material. A continuous softening was observed under all loading conditions. The introduction of hold periods to the applied cycles reduced material lifetime, with most prominent effects at technologically relevant small strain levels. The microstructural characterization reveals a coarsening of the original “martensitic” lath-type microstructure to a structure with polygonal subgrains and reduced dislocation density. The microstructural data forms the input for a physically-based modelling approach.
A numerical pre-study has shown that cracks in a flat sample featuring a drilled hole can be classified into one of three crack shape classes based on the combined evaluation of various types of test data.
Es wird eine neu entwickelte Methode zur Thermographiebasierten Rissmessung vorgestellt. Darüber hinaus wird eine numerische Vorarbeit präsentiert, die zeigt, dass anhand der gemeisamen Auswertung der Versuchsdaten aus unterschiedlicher Sensorik die Möglichkeit besteht, die unter Ermüdungsbelastung in Bohrlochproben auftretenden Risse in Geometriekategorien zu unterteilen.
Es konnte anhand einer numerischen Voruntersuchung gezeigt werden, dass anhand der kombinierten Auswertung der im Versuch verwendeten Sensorik eine Einteilung der unter Ermüdung in Bohrlochproben auftretenden Rissformen in verschiedene Hauptkategorien (Eckriss, Oberflächenriss, Durchgangsriss) möglich ist.
Additive manufacturing (AM) offers significantly greater freedom of design compared to conventional manufacturing processes since the final parts are built layer by layer. This enables metal AM, also known as metal 3D printing, to be utilized for improving efficiency and functionality, for the production of parts with very complex geometries, and rapid prototyping. However, despite many technological advancements made in recent years, several challenges hinder the mass adoption of metal AM. One of these challenges is mechanical anisotropy which describes the dependency of material properties on the material orientation. Therefore, in this work, stainless steel 316L parts produced by laser-based powder bed fusion are used to isolate and understand the root cause of anisotropy in AM parts. Furthermore, an efficient and accurate multiscale numerical framework is presented for predicting the deformation behavior of actual AM parts on the macroscale undergoing large plastic deformations. Finally, a novel constitutive model for the plastic spin is formulated to capture the influence of the microstructure evolution on the material behavior on the macroscale.
Virtual-lab-based determination of a macroscopic yield function for additively manufactured parts
(2019)
This work aims for a yield function description of additively manufactured parts of S316L steel at the continuum-mechanical macro-scale by means of so-called virtual experiments using a crystal plasticity (CP) model at meso-scale. Additively manufactured parts require the consideration of the specific process-related microstructure, which prevents this material to be macroscopically treated as isotropic, because of crystallographic as well as topological textures. From virtual experiments, yield loci under various loading conditions are simulated. The scale bridging from meso- to macro-scale is realised by the identification of the simulated yield loci as a modified anisotropic Barlat-type yield model representation.
Virtual-lab-based determination of a macroscopic yield function for additively manufactured parts
(2019)
This work aims for an yield function description of additively manufactured (AM) parts of S316L steel at the continuum-mechanical macro-scale by means of so-called virtual experiments using a crystal plasticity (CP) model at meso-scale. Additively manufactured parts require the consideration of the specific process-related microstructure, which prevents this material to be macroscopically treated as isotropic, because of crystallographic as well as topological textures.
EBSD/CT-Scans from in-house additively manufactured specimen extract the unique microstructural topology which is converted to a representative volume element (RVE) with grain structure and crystal orientations.
Crystal plasticity model parameters on this RVE are calibrated and validated by means of mechanical testing under different texture angles. From virtual experiments on this RVE, yield loci under various loading conditions are simulated. The scale bridging from meso- to macro-scale is realised by the identification of the simulated yield loci as a modified anisotropic Barlat-type yield model representation.
The underlying cause of mechanical anisotropy in additively manufactured (AM) parts is not yet fully understood and has been attributed to several different factors like microstructural defects, residual stresses, melt pool boundaries, crystallographic and morphological textures. To better understand the main contributing factor to the mechanical anisotropy of AM stainless steel 316L, bulk specimens were fabricated via laser powder bed fusion (LPBF). Tensile specimens were machined from these AM bulk materials for three different inclinations: 0◦, 45◦, and 90◦ relative to the build plate. Dynamic Young’s modulus measurements and tensile tests were used to determine the mechanical anisotropy. Some tensile specimens were also subjected to residual stress measurement via neutron diffraction, porosity determination with X-ray micro-computed tomography (μCT), and texture analysis with electron backscatter diffraction (EBSD). These investigations revealed that the specimens exhibited near full density and the detected defects were spherical. Furthermore, the residual stresses in the loading direction were between −74 ± 24 MPa and 137 ± 20 MPa, and the EBSD measurements showed a preferential ⟨110⟩ orientation parallel to the build direction. A crystal plasticity model was used to analyze the elastic anisotropy and the anisotropic yield behavior of the AM specimens, and it was able to capture and predict the experimental behavior accurately. Overall, it was shown that the mechanical anisotropy of the tested specimens was mainly influenced by the crystallographic texture.
A constitutive model for the mechanical behaviour of single crystal Ni-base superalloys under high temperature conditions has been developed in the framework of a Cooretec project in cooperation with Siemens AG, MTU Aero Engines AG and University Bayreuth. In addition to the conventional material properties e.g. elastic constants, the model requires the parameters of the initial microstructure as an input. Thus, the γ’-precipitate size and the channel width of the γ-matrix were obtained from SEM micrographs. The model uses the slip system theory and describes the movement, multiplication and annihilation of dislocations in the channels. Furthermore, the cutting of precipitates is another mechanism contributing to the plastic flow. The evolution of the morphology due to rafting and its effects on the deformation have been implemented according to. The kinematic hardening is introduced as a stress tensor to realistically represent the strain hardening of arbitrary oriented single crystals.
The mechanical behaviour of single crystal specimens has been experimentally investigated in tension tests at different strain rates and in creep tests under various loads. The constitutive model has been calibrated based on the experimental data for temperatures of 950°C and 850°C and the [001] and [111] crystallographic orientations.
Finally, a micromechanical model was created to simulate the creep response of additive manufactured polycrystalline structures. An EBSD image is taken to obtain the grain geometry and their respective orientation. The grain boundaries are discretised using cohesive elements, whereas the single crystal model was applied to each grain in the representative volume. The polycrystal model is generated using Dream3D, NetGen and other software previously developed at the BAM.
Systematic microstructure design requires reliable thermodynamic descriptions of each and all microstructure elements. While such descriptions are well established for most bulk phases, thermodynamic assessment of microstructure defects is challenging because of their individualistic nature. In this paper, a model is devised for assessing grain boundary thermodynamics based on available bulk thermodynamic data. We propose a continuous relative atomic density field and its spatial gradients to describe the grain boundary region with reference to the homogeneous bulk and derive the grain boundary Gibbs free energy functional. The grain boundary segregation isotherm and phase diagram are computed for a regular binary solid solution, and qualitatively benchmarked for the Pt–Au system. The relationships between the grain boundary's atomic density, excess free volume, and misorientation angle are discussed. Combining the current density-based model with available bulk thermodynamic databases enables constructing databases, phase diagrams, and segregation isotherms for grain boundaries, opening possibilities for studying and designing heterogeneous microstructures.
Interfacial Spinodals
(2021)
Despite their finite spatial extent, interfaces can have profound impacts on microstructure properties. This is because of their distinct phase-like properties distinguishing them from the adjacent homogeneous bulk structure. When noticed by solute atoms, interfaces can experience their own chemical phase changes. In this talk, we investigate the constrained states of chemically decomposed phases at interfaces. A density-based theory of interfaces is proposed to describe the confined chemical decomposition at general grain boundaries. Here the grain boundary is viewed as a lesser dense, defected structure with reference to the corresponding bulk structure. Using this picture, the thermodynamic origins of interfacial spinodal phenomena are revealed. We also show that transient interfacial spinodals can be activated over a large alloy composition range, enabling kinetic engineering of interfacial chemistry.
Segregation to grain boundaries affects their cohesion, corrosion, and embrittlement and plays a critical role in heterogeneous nucleation. In order to quantitatively study segregation and low-dimensional phase separation at grain boundaries, here, we apply a density-based phase-field model. The current model describes the grain-boundary thermodynamic properties based on available bulk thermodynamic data, while the grain-boundary-density profile is obtained using atomistic simulations. To benchmark the performance of the model, Mn grain-boundary segregation in the Fe–Mn system is studied. 3D simulation results are compared against atom probe tomography measurements conducted for three alloy compositions. We show that a continuous increase in the alloy composition results in a discontinuous jump in the segregation isotherm. The jump corresponds to a spinodal Phase separation at grain boundary. For alloy compositions above the jump, we reveal an interfacial transient spinodal phase separation.
The transient spinodal phenomenon opens opportunities for knowledge-based microstructure design through the chemical manipulation of grain boundaries. The proposed density-based model provides a powerful tool to study thermodynamics and kinetics of segregation and phase changes at grain boundaries.
The elastic energy of mixing for multi-component solid solutions is derived by generalizing Eshelby's sphere-in-hole model. By surveying the dependence of the elastic energy on the chemical composition and lattice misfit, we derive a lattice strain coefficient λ*. Studying several high-entropy alloys and superalloys, we propose that most solid solution multi-component alloys are stable when λ*<0.16, generalizing the Hume-Rothery atomic-size rule for binary alloys. We also reveal that the polydispersity index δ, frequently used for describing strain in multi-component alloys, directly represents the elastic energy e with e=qδ², q being an elastic constant. Furthermore, the effects of (i) the number and (ii) the atomic-size distribution of constituting elements on the phase stability of high-entropy alloys were quantified. The present derivations and discussions open for richer considerations of elastic effects in high-entropy alloys, offering immediate support for quantitative assessments of their thermodynamic properties and studying related strengthening mechanisms.
The Gibbs free energy of a grain boundary is a complex thermodynamic function of temperature, pressure, and composition. These complexities add to the intrinsic crystallographic and chemical constraints imposed by the adjacent bulk phase. Recently we have proposed a density-based model for assessing grain boundary thermodynamics that enables a CALPHAD-informed description of the grain boundary. As such, the Gibbs free energy of the grain boundary is directly linked with available CALPHAD thermodynamic data. In this talk, new aspects of interfacial segregation and phase transformation are revealed by benchmarking the current model for various experimental cases, including several steels, high-entropy alloys and aluminum alloys. The effects of elastic interactions on the grain boundary segregation and the application of the model to a nanocrystalline Pt-Au alloy, with numerous grain boundaries of various characters, will be discussed.
The onset of plasticity in a single crystal C60 fullerite was investigated by nanoindentation on the (111) crystallographic plane. The transition from elastic to plastic deformation in a contact was observed as pop-in events on loading curves. The respective resolved shear stresses were computed for the octahedral slip systems ⟨011¯¯¯⟩{111}, supposing that their activation resulted in the onset of plasticity. A finite element analysis was applied, which reproduced the elastic loading until the first pop-in, using a realistic geometry of the Berkovich indenter blunt tip. The obtained estimate of the C60 theoretical shear strength was about 1/11 of the shear modulus on {111} planes.
An improved diffusion model is proposed for pore annihilation during HIP of single-crystal nickel-base superalloys. The model assumes the pore dissolution by emission of vacancies and their sink to the low angle boundaries. Calculation, considering distribution of the pore sizes, predicts the kinetics of pore annihilation similar to the experimental one.
Initially, as-cast and homogenized single crystals of nickel-base superalloy CMSX-4 are subjected to hot isostatic pressing at 1288 °C. Two series of experiments are conducted: under the same pressure of 103 MPa but with different durations, between 0.5 and 6 h, and under different pressures, between 15 and 150 MPa, but for the same time of 0.5 h. The porosity annihilation is investigated metallographically and by high-resolution synchrotron X-ray tomography. The obtained experimental results are compared with the predictions of the vacancy model proposed recently in the group. Herein, the model is further refined by coupling with X-ray tomography. The model describes the evolution of the pore arrays enclosed in the 3D synchrotron tomograms during hot isostatic pressing and properly predicts the time and stress dependences of the pore annihilation kinetics. The validated model and the obtained experimental results are used for selecting the optimal technological parameters such as applied pressure and processing time
The elastic properties of the single-crystal nickel-base superalloy CMSX-4 used as a blade material in gas turbines were investigated by the sonic resonance method in the temperature interval between room temperature and 1300 °C. Elastic constants at such high temperatures are needed to model the mechanical behavior of blade material during manufacturing (hot isostatic pressing) as well as during technical accidents which may happen in service (overheating). High reliability of the results was achieved using specimens of different crystallographic orientations, exciting various vibration modes as well as precise measurement of the material density and thermal Expansion required for modeling the resonance frequencies by finite element method. Combining the results measured in this work and literature data the elastic constants of the gamma and gamma' phases were predicted.
This prediction was supported by measurement of the temperature dependence of the gamma'fraction.
All data obtained in this work are given in numerical or analytical forms and can be easily used for different scientific and engineering calculations.
The creep behavior of single crystals of the nickel-based superalloy CMSX-4 was investigated at 1288 °C, which is the temperature of the hot isostatic pressing treatment applied to this superalloy in the industry. It was found that at this super-solvus temperature, where no gammaPrime-strengthening occurs, the superalloy is very soft and rapidly deforms under stresses between 4 and 16 MPa. The creep resistance was found to be very anisotropic, e.g., the creep rate of [001] crystals was about 11 times higher than that of a [111] crystal. The specimens of different orientations also showed a very different necking behavior. The reduction of the cross-sectional area psi of [001] crystals reached nearly 100 pct, while for a [111] crystal psi = 62 pct. The EBSD analysis of deformed specimens showed that despite such a large local strain the [001] crystals did not recrystallize, while a less deformed [111] crystal totally recrystallized within the necking zone. The recrystallization degree was found to be correlated with deformation behavior as well as with dwell time at high temperature. From the analysis of the obtained results (creep anisotropy, stress dependence of the creep rate, traces of shear deformation, and TEM observations), it was concluded that the main strain contribution resulted from <01-1>{111} octahedral slip.
A γ-analogue of the superalloy CMSX-4 that does not contain the strengthening γ′ -phase and only consists of the γ-solid solution of nickel has been designed, solidified as single-crystals of different orientations, and tested under creep conditions in the temperature range between 1150 and 1288 °C. The tests have revealed a very high creep anisotropy of this alloy, as was previously found for CMSX-4 at supersolvus temperature of 1288 °C. This creep anisotropy could be explained by the dominance of 〈011〉{111} octahedral slip. Furthermore, the analysis of the creep data has yielded a high value of the creep activation energy, Qc≈442 kJ/mol, which correlates with the high activation energy of Re diffusion in Ni. This supports the hypothesis that dislocation motion in the γ-matrix of Re-containing superalloys is controlled by the diffusion of the Re atoms segregating at the dislocation core. The Norton stress exponent n is close to 5, which is a typical value for pure metals and their alloys. The absence of γ′ -reprecipitation after high-temperature creep tests facilitates microstructural investigations. It has been shown by EBSD that creep deformation results in an increasing misorientation of the existing low angle boundaries. In addition, according to TEM, new low angle boundaries appear due to reactions of the a/2 〈011〉 mobile dislocations and knitting of new networks.
Entwicklung der Mikrostruktur der mechanischen Eigenschaften und der Eigenspannungen in L-PBF 316L
(2022)
Die additive Fertigung (AM) metallischer Werkstoffe mittels Laser Powder Bed Fusion (L-PBF) ermöglicht einzigartige hierarchische Mikrostrukturen, die zu Verbesserungen bestimmter mechanischer Eigenschaften gegenüber konventionell hergestellten Varianten derselben Legierung führen können. Allerdings ist das L-PBF-Verfahren häufig durch das Vorhandensein hoher Eigenspannungen gekennzeichnet, die es zu verstehen und zu mindern gilt. Daher ist das Verständnis der Mikrostrukturen, der Eigenspannungen und der daraus resultierenden mechanischen Eigenschaften entscheidend für eine breite Akzeptanz bei sicherheitskritischen Anwendungen. Die BAM hat ein multidisziplinäres Forschungsprogramm gestartet, um diese Aspekte bei LPBF 316L zu untersuchen. Der vorliegende Beitrag stellt einige der wichtigsten Ergebnisse vor: der Einfluss von Prozessparametern auf die Mikrostruktur, der Einfluss von Mikrostruktur und Textur auf die Festigkeit, Kriechverhalten und Schädigung und die Stabilität von Eigenspannungen und Mikrostruktur unter Wärmebehandlungsbedingungen.
Laser powder bed fusion (LPBF) is an additive manufacturing process for materials which inherently tends to yield various degrees of metastable hierarchical microstructures, defects and high residual stresses in the as-built condition depending on the process parameters. The understanding of the evolution of these typical features during heat treatment and subsequent thermal and mechanical ageing is crucial for the wider acceptance for safety critical structures. A multi-disciplinary research project at BAM studying the development of the microstructure, defects, residual stresses typical of LPBF 316L and their evolution during thermal and mechanical ageing has led to insights into the stability of these inherent features. This presentation aims to give a broad overview of the project with a few specific cases of investigation. Firstly, the formation of residual stresses, the nature of the initial microstructure, the tensile properties and a modelling approach to understand the anisotropy will be presented. This will be followed by examples of studies of their evolution during heat treatment, long term thermal exposure, and room temperature and high temperature mechanical testing compared to a baseline of conventional wrought variant of the same alloy.
An overview of the BAM funed Focus Area Materials Project "AGIL" will be presented. AGIL focussed on the stdiy of the ageing characteristics of additively manufactured austenitic stainless steel with a "powder to mechanical failure" Approach. Recent Highlights are presented and a perspective for future studies.
With the introduction of a mass transport mechanism the entire problem is subjected to a time frame that dictates the time-dependent action of soluted species on mechanical properties. A numerical framework within the phase-field approach is presented with an embrittlement-based coupling mechanism. The underlying functionals are expressed in terms of the displacement, mass concentration and crack phase-field. Within the phase-field approach the modelling of sharp crack discontinuities is replaced by a diffusive crack model facilitating crack initiation and complex crack topologies without the requirement of a predefined crack path. The isotropic hardening of the elasto-plastic deformation model and the local fracture criterion are affected by the species concentration. This allows for embrittlement and leads to an accelerated crack propagation. An extended mass transport equation for hydrogen embrittlement, accounting for mechanical stresses and deformations, is implemented. For stabilisation purposes a staggered scheme is applied to solve the system of partial differential equations by a multi-field finite-element method. A thermodynamically consistent coupling relation that accommodates the required mechanisms is presented.
Modelling of environmentally assisted material degradation in the crack phase-field framework
(2019)
The simulation of crack propagation was conducted with a diffusive crack model in a variational framework. Moreover, the physically sound introduction of mass transport and coupling mechanisms due to environmentally assisted effects could be realised in this framework. The objective consists of the application of the phase-field Approach towards the simulation of environmentally assisted material degradation with the advantage of a non-required predefined crack path and a mesh-independent non-local formulation that facilitates the damage evolution with respect to material softening. The sharp crack is regularised by the introduction of a phase-field order parameter leading to a diffusive crack formulation.
Besides the equations originating from the linear momentum balance an additional evolution equation for the crack phase-field is introduced. Furthermore, mass transport is simulated by a Diffusion equation. The description delivered by the variational phase-field framework is able to simulate crack propagation according to published numerical test cases.
Additionally, the calculation of stress intensity factors is possible as well as crack resistance curves that describe stable crack propagation.
Virtual-lab-based determination of a macroscopic yield function for additively manufactured parts
(2018)
This work presents a method for the yield function determination of additively manufactured parts of S316L steel. A crystal plasticity model is calibrated with test results and used afterwards to perform so-called virtual experiments, that account for the specific process-related microstructure including crystallographic and morphological textures. These simulations are undertaken on a representative volume element (RVE), that is generated from EBSD/CT-Scans on in-house additively manufactured specimen, considering grain structure and crystal orientations. The results of the virtual experiments are used to determine an anisotropic Barlat yield function, that can be used in a macroscopical continuum-sense afterwards. This scale-bridging approach enables the calculation of large-scale parts, that would be numerically too expensive to be simulated by a crystal plasticity model.
Virtual-lab-based determination of a macroscopic yield function for additively manufactured parts
(2018)
This work presents a method for the yield function determination of additively manufactured parts of S316L steel. A crystal plasticity model is calibrated with test results and used afterwards to perform so-called virtual experiments, that account for the specific process-related microstructure including crystallographic and morphological textures. These simulations are undertaken on a representative volume element (RVE), that is generated from EBSD/CT-Scans on in-house additively manufactured specimen, considering grain structure and crystal orientations. The results of the virtual experiments are used to determine an anisotropic Barlat yield function, that can be used in a macroscopical continuum-sense afterwards. This scale-bridging approach enables the calculation of large-scale parts, that would be numerically too expensive to be simulated by a crystal plasticity model.
The elastic constants are the most basic mechanical properties of a material and are needed for any structural analysis of a component. For example, they have a major influence on the eigenfrequencies of vibrating parts. Single crystals of Ni-base superalloys are strongly anisotropic, which means that the observed properties are orientation dependent. Tensor algebra is then required to mathematically formulate the elastic properties and their relations to the crystal orientation. Hence, this chapter first summarizes some basic definitions and calculation rules for Rotation matrices, including the definition of the Euler angles, which are most commonly used to define the relative orientations of the crystal and the component. Parts of this chapter closely follow the lines of the excellent exposition of the topic by Olschewski.
The large number of TEM investigations and the regular microstructure
of single-crystal nickel-base superalloys has boosted the development of a number of physically motivated constitutive laws. In contrast to the more phenomenological models discussed in the next chapter, these models use dislocation densities as internal variables. Obvious advantages are that the computed densities can be compared to TEM observations and the Deformation mechanisms can be easier translated into mathematical equations.
Mit Hilfe des HIP-Verfahrens („Hot Isostatic Pressing“) werden Poren in der einkristallinen Nickel-Basis Superlegierung CMSX-4 kontinuierlich geschrumpft und dadurch die nach der Erstarrung und der Wärmebehandlung vorhandene Porosität stark reduziert. In diesem Beitrag werden experimentelle und numerische Untersuchungen zu den Mechanismen der Porenschrumpfung zusammengefasst. Es zeigt sich, dass das Verformungsverhalten während Kriechversuchen bei der HIP-Temperatur durch Versetzungsgleitung auf oktaedrischen Ebenen dominiert wird.
Dagegen zeigen Messungen der Porositätsabnahme und Simulationen des Porenschließens, dass die Kinetik der Porenschrumpfung durch das Phänomen der Leerstellendiffusion zwischen Poren und Kleinwinkelkorngrenzen („Low Angle Boundary“, LAB) bestimmt wird. Im Gegensatz führt die klassische Kristallviskoplastizität zu einer systematischen Überschätzung dieser Kinetik. Der scheinbare Widerspruch lässt sich auflösen, wenn man bedenkt, dass auf der Skala der Poren Versetzungsquellen nicht gleichmäßig verteilt sind, wie in der konventionellen Kristallplastizität implizit angenommen wird. Stattdessen wird in einem weiterführenden Modell davon ausgegangen, das Kleinwinkelkorngrenzen (LABs) als Versetzungsquellen fungieren, während die Scherspannungen sehr stark in der Nähe der Poren lokalisiert sind, was die Emission von Versetzungen deutlich reduziert.
The creep behavior of single-crystals of the nickel-base superalloy CMSX-4 was investigated at 1288°C, which is the temperature of the hot isostatic pressing (HIP) treatment applied to this superalloy in the industry. It was found that at this super-solvus temperature, where no Gamma’-strengthening occurs, the superalloy is very soft and rapidly deforms under stresses between 4 and 16 MPa. The creep resistance was found to be very anisotropic, e.g. the creep rate of [001] crystals was about 11 times higher than that of a [111] crystal. The specimens of different orientations also showed a very different necking behavior. The reduction of the cross-section area psi of [001] crystals reached nearly 100%, while for a [111] crystal psi=62%. The EBSD analysis of deformed specimens showed that despite such a large local strain the [001] crystals didn’t not recrystallize, while a less deformed [111] crystal totally recrystallized within the necking zone. From the shape of deformed specimens and TEM investigations it was concluded that the main strain contribution resulted from <011> {111} octahedral slip.
Turbine blades often contain cylindric holes used to generate an air film that protects the blade alloy from the hot gases. These cooling holes of diameter around one mm are drilled by laser through the thickness of the blades. Unfortunately, the resulting stress concentration and the drilling-induced damage are known to favor crack initiation from the holes. It is thus necessary to assess the impact of these cooling holes on the structural integrity of the blades. Since cracks initiate very readily, the fatigue life of the components is mainly controlled by the propagation of the cracks in the stress gradient induced by the holes.
For this purpose, displacement controlled high-temperature LCF (Low-Cycle-Fatigue) tests were performed with center hole specimens of a coarse-grained Nickel base Superalloy. The tests were stopped after a defined load drop. In addition, crack propagation tests with Double Edge Notch specimens were performed. Moreover, specimens with different hole surface finishes were investigated, which showed a detrimental effect of the hole surface roughness. In parallel, an evaluation of the LCF tests based on a fracture mechanics-based model (Madia et al., Eng. Fract. Mech., 2018) has been applied. Thereby, the specimen life is controlled by the crack propagation time until failure. Crack growth is controlled by a modified NASGRO equation accounting for large-scale yielding and a progressive build-up of crack closure. The initial crack size has been derived from the measurements of defects around the borehole. A reasonable agreement between predicted and measured lifetimes is observed if one keeps in mind the large uncertainty regarding the effective shape of the cracks.
The scatter of fatigue crack growth data can become significant for coarse grained materials. By using a probabilistic description of crack propagation as the foundation of a lifetime prediction model, lifetime scatter of laboratory specimens can be reproduced. However, the lifetime of real components is subjected to additional scattering factors such as surface condition or uncertainty regarding direction and shape of emerging cracks. These factors need to be addressed in order to exploit the advantages of probabilistic description, i.e. the reduction of unnecessary conservatisms.
High temperature LCF (Low-Cycle-Fatigue) tests were performed with center hole specimens of a coarse-grained Nickel base Superalloy. In addition, crack propagation tests with Double Edge Notch specimens were performed. A procedure to detect the shape of the starting crack that combines the potential drop method and induction thermography was developed. The geometry and the number of notches were varied. Moreover, specimen with different hole surface finishes were investigated, which showed a detrimental effect of the roughness of the hole surface.
The results have been compared to predictions of a probabilistic tool for the estimation of lcf lifetime that has been calibrated beforehand on laboratory specimens with and without notches. In parallel, a fracture mechanics-based lifetime model was developed, which includes the initial crack size as a critical parameter. Thereby, the influence of the large grain size (>1mm), and the shape of the starting crack at the notch were considered.
Acknowledgements
The investigations are conducted as part of the joint research program COOREFLEX-Turbo in the frame of AG Turbo. The work is supported by Siemens AG and the Bundesministerium für Wirtschaft und Technologie (BMWi) as per resolution of the German Federal Parliament under grant number 03ET7071E.
Single crystal superalloys usually contain pores of sizes 5-10 micro-m after casting and heat treatment. These pores can be reduced under compression by combined creep and diffusion in a subsequent treatment called Hot Isostatic Pressing (HIP). The paper presents a methodology to simulate pore shrinkage under HIP conditions in two dimensions (2D).
At the scale of the pores, which is also the scale of the sub-grains (<50 micro-m) the dislocation sources cannot be assumed to be homogeneously distributed. Thus, the applicability of classical crystal plasticity is questionable. In this case, the transport of dislocations under an applied stress from the location where they are nucleated must be explicitly modelled. This is done by solving the transport equations for the dislocation densities and the elasticity equations in 2D. The dislocations are assumed to be nucleated at Low Angle Boundaries. They glide or climb through the sub-grains with a stress dependent velocity.
The transport equations are solved by the Flux-Corrected Transport method, which belongs to the predictor-corrector class of algorithms. In the first step, an artificial diffusion is introduced, which suppresses spurious oscillations of the solution. In a second step, the solution is corrected in such a way that no additional extremes appear and that the extremes do not grow. The algorithm is validated by simulating the transport of simple distributions with a constant velocity field.
With the dislocation velocities and the computed dislocation densities, the inelastic shear rate at the slip system level is computed by integrating the Orowan equation. In the 2D-setting, three slip systems are considered. The contributions of these slip systems are summed up to obtain the total inelastic strain rate. Dislocation glide and climb and the coupling of climb with vacancies diffusion are considered.
The resolution of the equilibrium equations from the inelastic strains turned out to be prone to numerical instabilities. As an alternative, the stresses are directly computed from the distribution of geometrically necessary dislocations following the method presented in. The resulting boundary value problem is solved by the Least-Square Finite Element method.
Examples of simulations are presented for a representative region under creep tension and for a pore shrinking under external pressure.
Mit Hilfe des HIP-Verfahrens („Hot Isostatic Pressing“) werden Poren in der einkristallinen Nickel-Basis Superlegierung CMSX-4 kontinuierlich geschrumpft und dadurch die nach der Erstarrung und der Wärmebehandlung vorhandene Porosität stark reduziert. In diesem Beitrag werden experimentelle und numerische Untersuchungen zu den Mechanismen der Porenschrumpfung zusammengefasst. Es zeigt sich, dass das Verformungsverhalten während Kriechversuchen bei der HIP-Temperatur durch Versetzungsgleitung auf oktaedrischen Ebenen dominiert wird.
Dagegen zeigen Messungen der Porositätsabnahme und Simulationen des Porenschließens, dass die Kinetik der Porenschrumpfung durch das Phänomen der Leerstellendiffusion zwischen Poren und Kleinwinkelkorngrenzen („Low Angle Boundary“, LAB) bestimmt wird. Im Gegensatz führt die klassische Kristallviskoplastizität zu einer systematischen Überschätzung dieser Kinetik. Der scheinbare Widerspruch lässt sich auflösen, wenn man bedenkt, dass auf der Skala der Poren Versetzungsquellen nicht gleichmäßig verteilt sind, wie in der konventionellen Kristallplastizität implizit angenommen wird. Stattdessen wird in einem weiterführenden Modell davon ausgegangen, das Kleinwinkelkorngrenzen (LABs) als Versetzungsquellen fungieren, während die Scherspannungen sehr stark in der Nähe der Poren lokalisiert sind, was die Emission von Versetzungen deutlich reduziert.
Single-crystal components made of nickel base superalloys contain pores after casting and homogenization heat treatment. Hot isostatic pressing (HIP), which is carried above the γ' -solvus temperature of the alloy, is industrially applied to reduce porosity. A modeling of HIP based on continuously distributed dislocations is developed in a 2D setting. Glide and climb of straight-edge dislocations, as well as vacancy diffusion, are the deformation mechanisms taken into account.
Thereby, dislocation glide is controlled by dragging a cloud of large atoms, and climb is controlled by vacancy diffusion. Relying on previous investigations of the creep behavior at HIP temperatures, it is assumed that new dislocations are nucleated at low-angle boundaries (LAB) and move through subgrains until they either reach the opposite LABs or react with other dislocations and annihilate.
Vacancies are created at the pore surface and diffuse through the alloy until they are either consumed by climbing dislocations or disappear at the LABs. The field equations are solved by finite elements. It is shown that pore shrinking is mostly controlled by vacancy diffusion as the shear stresses at the LABs are too low to
nucleate a sufficient amount of dislocations.
The finite element software ProCrackPlast is developed for the automated simulation of fatigue crack growth in arbitrarily loaded three-dimensional components with large scale plastic deformations, in particular under cyclic thermomechanical (TMF) loading.
ProCrackPlast is developed by extending the software Procrack, created at TU Freiberg for linear-elastic crack growth simulations. The pre-processing, FEM analysis, and the post-processing in ProCrackPlast are done by the commercial software ABAQUS. ProCrackPlast resorts to a crack growth procedure which adaptively updates the crack in finite increments based on the fracture-mechanical parameter, cyclic crack tip opening displacement ΔCTOD . Features of this software along with two application examples of fatigue crack growth in a typical cast steel, Ni-Resist, are presented in this paper to show its capability and performance.
The cyclic crack tip opening displacement ΔCTOD is a promising loading Parameter to quantify the crack propagation under thermomechanical fatigue (TMF). In this work, suitable techniques are investigated and compared for an accurate calculation of ΔCTOD under TMF loading using a viscoplastic temperature dependent material model.
It turned out that collapsed special crack tip elements give the best results. An efficient FEM-technique is developed to simulate the incremental crack growth by successive remeshing, whereby the deformations and internal state variables are mapped from the old mesh onto the new one. The developed techniques are demonstrated and discussed for two-dimensional examples like TMF-specimens. Recommendations are made regarding important numerical control parameters like optimal size of crack tip elements, length of crack growth increment in relation to plastic zone size and ΔCTOD value.
The cyclic crack tip opening displacement ΔCTOD is a promising loading Parameter to quantify the crack propagation under thermomechanical fatigue (TMF). In this work, suitable techniques are investigated and compared for an accurate calculation of ΔCTOD under TMF loading using a viscoplastic temperature dependent material model.
It turned out that collapsed special crack tip elements give the best results. An efficient FEM-technique is developed to simulate the incremental crack growth by successive remeshing, whereby the deformations and internal state variables are mapped from the old mesh onto the new one. The developed techniques are demonstrated and discussed for two-dimensional examples like TMF-specimens. Recommendations are made regarding important numerical control parameters like optimal size of crack tip elements, length of crack growth increment in relation to plastic zone size and ΔCTOD value.
The cyclic crack tip opening displacement ΔCTOD is a promising loading Parameter to quantify the crack propagation under thermomechanical fatigue (TMF). In this work, suitable techniques are investigated and compared for an accurate calculation of ΔCTOD under TMF loading using a viscoplastic temperature dependent material model.
It turned out that collapsed special crack tip elements give the best results. An efficient FEM-technique is developed to simulate the incremental crack growth by successive remeshing, whereby the deformations and internal state variables are mapped from the old mesh onto the new one. The developed techniques are demonstrated and discussed for two-dimensional examples like TMF-specimens. Recommendations are made regarding important numerical control parameters like optimal size of crack tip elements, length of crack growth increment in relation to plastic zone size and ΔCTOD value.
Die Steigerung der Leistung und des thermodynamischen Wirkungsgrades von Verbrennungsmotoren führt zu erhöhten Anforderungen an die Festigkeit abgasführender Komponenten wie z. Bsp. Abgasturbolader. Als Folge erhöhter thermomechanischer Wechselbeanspruchungen (TMF) im Betrieb kommt es an den mechanisch und/oder thermisch höchst beanspruchten Stellen der Bauteile zur Bildung von Rissen, wodurch die Lebensdauer der Komponenten begrenzt wird. Derzeit werden bei Turboladern heißgehende Bauteile mit detektierten Rissen zumeist prophylaktisch ersetzt, da die weitere Ausbreitung der Risse während des Betriebs nicht vorhergesagt werden kann. Um diese aufwändige und un- ökonomische Praxis zu vermeiden, wurde im vorliegenden Forschungsvorhaben eine rechnerische Bewertungsmethode auf Basis der experimentellen und numerischen Bruchmechanik erarbeitet, mit deren Hilfe bereits in der Auslegungsphase oder während des Betriebs die (restliche) Lebensdauer der abgasführenden Komponenten vorhergesagt werden kann. Damit wird erstmalig die quantitative Vorhersage der Rissentwicklung bei TMF-Beanspruchungsbedingungen unter Berücksichtigung großer zyklischer viskoplastischer Verformungen ermöglicht.
Zentrales Ergebnis des Vorhabens ist eine automatisierte Berechnungsprozedur auf der Basis spezieller Finite-Elemente-Techniken (FEM), womit sowohl der Pfad als auch die Größe eines Risses als Funktion der Anzahl der Lastwechsel in Bauteilen unter TMF-Bedingungen berechnet werden kann. Als geeigneter Beanspruchungsparameter zur Bewertung des Rissfortschritts unter TMF wurde die zyklische Rissöffnungsverschiebung ΔCTOD verwendet. Das Werkstoffverhalten des betrachteten austenitischen Gusseisens Ni-Resist D-5S wurde mit einem validierten viskoplastischen, temperaturabhängigen Materialmodell modelliert, das zur Berücksichtigung große Verzerrungen und Rotationen am Riss erweitert wurde. Für die genaue Berechnung des ΔCTOD bei TMF wurden effiziente FEM-Techniken erarbeitet. Zur Simulation der Rissausbreitung wurde ein automatischer FEM-Algorithmus mit inkrementeller adaptiver Neuvernetzung entwickelt, bei dem die Verformungen und inelastischen Zustandsvariablen jeweils vom alten auf das neue Netz übertragen werden. Dieser Algorithmus wurde im Software-Paket ProCrackPlast implementiert, das in Verbindung mit dem kommerziellen FEM-Code Abaqus zur Lösung dreidimensionaler Rissprobleme zur Verfügung steht. Ziel der umfangreichen experimentellen Arbeiten war es, an isothermen LCF und anisothermen TMF-Versuchen mit gekerbten Flachzugproben (SENT) das Risswachstum im Temperaturbereich von 20 °C bis 700 °C zu ermitteln. Mit Hilfe begleitender 2D FEM Simulationen wurden anhand dieser Datenbasis die Rissfortschrittskurven des Werkstoffs unter Anwendung des ΔCTOD-Konzepts bestimmt und in geeigneter, parametrisierter Form den Nutzern zur Verfügung gestellt. Die Versuche an SENT-Proben wurden mit der entwickelten Software ProCrackPlast als 3D Modell simuliert. Der Vergleich der 2D und 3D Simulationen ergab einen systematischen Unterschied im CTOD und CTOD, der mit Hilfe eines Übertragungsfaktors korrigiert wurde. Der Vergleich der 3D Berechnungen mit den Experimenten zeigte eine zufriedenstellende Übereinstimmung der er- reichten Risslänge mit der Zahl der Lastzyklen im gesamten Temperaturbereich, wobei die numerische Prognose meist auf der konservativen / sicheren Seite lag. Die Übertragbarkeit der Ergebnisse der 2D Parameteridentifikation auf 3D Risskonfigurationen mit Mixed-Mode Beanspruchung ist mit zusätzlichen Versagenshypothesen verbunden, die aufgrund fehlender Versuchsdaten im Vorhaben nicht endgültig geklärt werden konnten. Zur Validierung des Gesamtkonzeptes wurden LCF-Proben mit einem bauteil- typischen Oberflächenriss experimentell und numerisch untersucht. In der Simulation konnte die komplexe Form und Größe der Rissentwicklung zufriedenstellend (richtig) vorhergesagt werden. Die Leis- tungsfähigkeit der erarbeiteten rechnerische Bewertungsmethode wurde an weiteren TMF-Beispielen vorgestellt und diskutiert. Die Software ProCrackPlast und die viskoplastische Materialroutine wurden dem Anwenderkreis des Vorhabens zusammen mit einem Nutzer-Handbuch und Verifikationsbeispielen zur Verfügung gestellt.
Das Ziel des Forschungsvorhabens ist erreicht worden.
The sound numerical prediction of welding-induced thermal stresses, residual stresses, and distortions strongly depends on the accurate description of a welded material’s thermomechanical deformation behaviour. In this work, we provide experimental data on the viscoplastic deformation behaviour of a grade-s960ql steel up to a temperature of 1000 ◦C. In addition, a multi-phase viscoplastic material model is proposed, which accounts for the experimentally observed isothermal deformation behaviour of grade-s960ql steel base and austenitised material, as well as for athermal contributions that originate from solid-state phase transformations. The multi-phase viscoplastic and a classic rateindependent isotropic hardening material model were applied in the numerical simulations of both-ends-fixed bar Satoh tests and a single-pass gas metal arc weld. The influence of material modelling choices on the agreement between numerical simulation and experimental results is discussed, and recommendations for further work are given.
Due to their excellent creep resistance and good oxidation resistance, 9–12% Cr ferritic–martensitic stainless steels are widely used as high temperature construction materials in power plants. However, the mutual combination of different loadings (e.g., creep and fatigue), due to a “flexible” operation of power plants, may seriously reduce the lifetimes of the respective components. In the present study, low cycle fatigue (LCF) and relaxation fatigue (RF) tests performed on grade P92 helped to understand the behavior of ferritic–martensitic steels under a combined loading. The softening and lifetime behavior strongly depend on the temperature and total strain range. Especially at small strain amplitudes, the lifetime is seriously reduced when adding a hold time which indicates the importance of considering technically relevant small strains.
Results of an extended TMF test program on grade P92 steel in the temperature range of 620 °C - 300 °C, comprising in-phase (IP) and out-of-phase (OP) tests, partly performed with symmetric dwells at Tmax/Tmin, are presented. In contrast to previous studies, the low-strain regime is also illuminated, which approaches flexible operation in a power plant with start/stop cycles. At all strain amplitudes, the material performance is characterized by continuous cyclic softening, which is retarded in tests at lower strains but reaches similar magnitudes in the course of testing. In the investigated temperature range, the phase angle does not affect fatigue life in continuous experiments, whereas the IP condition is more detrimental in tests with dwells. Fractographic analyses indicate creep-dominated and fatigue-dominated damage for IP and OP, respectively. Analyses of the (micro)hardness distribution in the tested specimens suggest an enhanced microstructural softening in tests with dwell times for the low- but not for the high-strain regime. To rationalize the obtained fatigue data, the fracture-mechanics-based D_TMF concept, which was developed for TMF life assessment of ductile alloys, was applied. It is found that the D_TMF parameter correlates well with the measured fatigue lives, suggesting that subcritical growth of cracks (with sizes from a few microns to a few millimeters) governs failure in the investigated range of strain amplitudes.
The methods of computational damage mechanics are well-established for the description of degradation of materials under monotone loading. An extension to structural damage induced by cyclic loading is however significantly limited. This is due to enormous computational costs required to resolve each load cycle by conventional temporal incremental integration schemes while a typical fatigue loading history comprises between thousands and millions of cycles. Despite the permanent increase of computational resources and algorithmic performance, a successful approach is rather based on the development of novel multiscale in time integration schemes.
A Fourier transformation-based temporal integration (FTTI) is represented, which takes advantage of temporal scale separation incorporated into the cycle jump method. The response fields are approximated by a Fourier series whose coefficients undergo the evolution on a long-time scale. This is correlated with the evolution of the history variables, including damage, by means of the adaptive cycle jump method of various orders. The necessary extrapolation rates are obtained from the underlying solution of a short-time scale problem, which results from the oscillatory boundary condition and fulfills the global equilibrium of the Fourier coefficients. In this way, a remarkable speedup is achieved because the number of cycles to be fully integrated dramatically decreases.
The key idea behind the FTTI method is that the global in space equilibrium problem is linear since it is decoupled from the evolution equations. The latter are solved in the quadrature points under response fields prescribed throughout the whole load cycle. Consequently, integration of a single load cycle is much more efficient than the conventional single scale integration where the global equilibrium iteration and the local iteration of the evolution equations are coupled. This results in an additional speedup of the FTTI method.
The performance of the FTTI technique is demonstrated for two different constitutive behaviors: a viscoplastic model with a damage variable governed by the local equivalent viscoplastic strain; a quasi-brittle response where the damage variable is driven by a non-local equivalent strain. The latter is implicitly introduced as proposed by Peerlings. Both, the explicit and implicit extrapolation schemes are validated. The FTTI solutions agree very well with the reference cycle-by -cycle solutions, while significantly reducing the computational costs. The adaptive determination of the jump length can properly recognize the particular responses throughout the fatigue loading history (stationary fatigue, acceleration of fatigue damage when approaching failure) as well as stress redistribution phenomena.
A continuum damage model for concrete is developed with a focus on fatigue under compressive stresses. This includes the possibility to model stress redistributions and capture size effects. In contrast to cycle based approaches, where damage is accumulated based on the number of full stress cycles, a strain based approach is developed that can capture cyclic degradation under variable loading cycles including different amplitudes and loading frequencies. The model is designed to represent failure under static loading as a particular case of fatigue failure after a single loading cycle. As a consequence, most of the material parameters can be deduced from statictests. Only a limit set of additional constitutive parameters is required to accurately describe the evolution under fatigue loading. Another advantage of the proposed model is the possibility to directly incorporate other multi-physics effects such as creep and shrinkage or thermal loading on the constitutive level. A multiscale approach in time is presented to enable structural computations of fatigue failure with a reduced computational effort. The damage rate within the short time scale corresponding to a single cycle is computed based on a Fourier based approach. This evolution equation is then solved on the long time scale using different implicit and explicit time integration schemes. Their performance and some limitations for specific loading regimes is discussed.
A continuum damage model for concrete is developed with a focus on fatigue under compressive stresses. This includes the possibility to model stress redistributions and capture size effects. In contrast to cycle based approaches, where damage is accumulated based on the number of full stress cycles, a strain based approach is developed that can capture cyclic degradation under variable loading cycles including different amplitudes and loading frequencies. The model is designed to represent failure under static loading as a particular case of fatigue failure after a single loading cycle. As a consequence, most of the material parameters can be deduced from statictests. Only a limit set of additional constitutive parameters is required to accurately describe the evolution under fatigue loading. Another advantage of the proposed model is the possibility to directly incorporate other multi-physics effects such as creep and shrinkage or thermal loading on the constitutive level. A multiscale approach in time is presented to enable structural computations of fatigue failure with a reduced computational effort. The damage rate within the short time scale corresponding to a single cycle is computed based on a Fourier based approach. This evolution equation is then solved on the long time scale using different implicit and explicit time integration schemes. Their performance and some limitations for specific loading regimes is discussed.
Computational methods for lifetime prediction of metallic components under high-temperature fatigue
(2019)
The issue of service life prediction of hot metallic components subjected to cyclic loadings is addressed. Two classes of lifetime models are considered, namely, the incremental lifetime rules and the parametric models governed by the fracture mechanics concept. Examples of application to an austenitic cast iron are presented. In addition, computational techniques to accelerate the time integration of the incremental models throughout the fatigue loading history are discussed. They efficiently solve problems where a stabilized response of a component is not observed, for example due to the plastic strain which is no longer completely reversed and accumulates throughout the fatigue history. The performance of such an accelerated Integration technique is demonstrated for a finite element simulation of a viscoplastic solid under repeating loading–unloading cycles.
A phase field model of brittle fracture has been developed to simulate the Hertzian crack induced by penetration of a rigid sphere to an isotropic linear-elastic half-space. The fracture formation is regarded as a diffusive field variable, which is zero for the intact material and unity if there is a crack. Crack growth is assumed to be driven by a strain invariant. The numerical implementation is performed with the finite element method and an implicit time integration scheme. The mechanical equilibrium and the phase field equations are solved in a staggered manner, sequentially updating the displacement field and the phase field variable. Numerical examples demonstrate the capability of the model to reproduce the nucleation and growth of the Hertzian cone crack.
Lifetime aspects including fatigue failure of concrete structures were traditionally only of minor importance. Because of the growing interest in maxing out the capacities of concrete, its fatigue failure under compression has become an issue. A variety of interacting phenomena such as e.g. loss of prestress, degradation due to chemical reactions or creep and shrinkage influence the fatigue resistance. Failure due to cyclic loads is generally not instantaneous, but characterized by a steady damage accumulation. Therefore, a reliable numerical model to predict the performance of concrete over its lifetime is required, which accurately captures order effects and full three-dimensional stress states.
Many constitutive models for concrete are currently available, which are applicable for specific loading regimes, different time scales and different resolution scales.
However, a key limitation of those models is that they generally do not address issues related to fatigue on a structural level. Very few models can be found in the literature that reproduce deterioration of concrete under repeated loading-unloading cycles. This is due to the computational effort necessary to explicitly resolve every cycle which exceeds the currently available computational resources. The limitation can only be overcome by the application of multiscale methods in time.
The objective of the paper is the development of numerical methods for the simulation of concrete under fatigue loading using temporal multiscale methods.
First, a continuum damage model for concrete is developed with a focus on fatigue under compressive stresses [1]. This includes the possibility to model stress redistributions and capture size effects. In contrast to cycle based approaches, where damage is accumulated based on the number of full stress cycles, a strain based approach is developed that can capture cyclic degradation under variable loading cycles including different amplitudes and loading frequencies. The model is designed to represent failure under static loading as a particular case of fatigue failure after a single loading cycle. As a consequence, most of the material parameters can be deduced from static tests. Only a limit set of additional constitutive parameters is required to accurately describe the evolution under fatigue loading. Another advantage of the proposed model is the possibility to directly incorporate other multi-physics effects such as creep and shrinkage or thermal loading on the constitutive level.
Second, a multiscale approach in time is presented to enable structural computations of fatigue failure with a reduced computational effort. The damage rate within the short time scale corresponding to a single cycle is computed based on a Fourier based approach [2]. This evolution equation is then solved on the long time scale using different implicit and explicit time integration schemes. Their performance and some limitations for specific loading regimes is discussed.
Finally, the developed methods will be validated and compared to experimental data.
[1] Vitaliy Kindrachuk, Marc Thiele, Jörg F. Unger. Constitutive modeling of creep-fatigue interaction for normal strength concrete under compression, International Journal of Fatigue, 78:81-94, 2015
[2] Vitaliy Kindrachuk, Jörg F. Unger. A Fourier transformation-based temporal integration scheme for viscoplastic solids subjected to fatigue deterioration, International Journal of Fatigue, 100:215-228, 2017
A key limitation of the most constitutive models that reproduce a Degradation of quasi-brittle materials is that they generally do not address issues related to fatigue. One reason is the huge computational costs to resolve each load cycle on the structural level. The goal of this paper is the development of a temporal Integration scheme, which significantly increases the computational efficiency of the finite element method in comparison to conventional temporal integrations.
The essential constituent of the fatigue model is an implicit gradient-enhanced formulation of the damage rate. The evolution of the field variables is computed as amultiscale Fourier series in time.On a microchronological scale attributed to single cycles, the initial boundary value problem is approximated by linear BVPs with respect to the Fourier coefficients. Using the adaptive cycle jump concept, the obtained damage rates are transferred to a coarsermacrochronological scale associated with the duration of material deterioration. The performance of the developedmethod is hence improved due to an efficient numerical treatment of the microchronological problem in combination with the cycle jump technique on the macrochronological scale. Validation examples demonstrate the convergence of the obtained solutions to the reference simulations while significantly reducing the computational costs.
Grain boundary energy effect on grain boundary segregation in an equiatomic high-entropy alloy
(2020)
Grain boundary (GB) Segregation has a substantial effect on the microstructure evolution and properties of polycrystalline alloys. The mechanism of nanoscale segregation at the various GBs in multicomponent alloys is of great challenge to reveal and remains elusive so far. To address this issue, we studied the GB segregation in a representative equiatomic FeMnNiCoCr high-entropy alloy (HEA) aged at 450 °C. By combining transmission Kikuchi diffraction, atom probe tomography analysis and a density-based thermodynamics modeling, we uncover the nanoscale segregation behavior at a series of well-characterized GBs of different characters. No segregation occurs at coherent twin boundaries; only slight nanoscale segregation of Ni takes place at the low-angle GBs and vicinal \Sigma 29b coincidence site lattice GBs. Ni and Mn show cosegregation of high levels at the general high-angle GBs with a strong depletion in Fe, Cr, and Co. Our density-based thermodynamic model reveals that the highly negative energy of mixing Ni and Mn is the main driving force for nanoscale cosegregation to the GBs. This is further assisted by the opposite segregation of Ni and Cr atoms with a positive enthalpy of mixing. It is also found that GBs of higher interfacial energy, possessing lower atomic densities (higher disorder and free volume), show higher segregation levels. By clarifying the origins of GB segregations in the FeMnNiCoCr HEA, the current work provides fundamental ideas on nanoscale segregation at crystal defects in multicomponent alloys.
The effectiveness of the mechanism of precipitation strengthening in metallic alloys depends on the shapes of the precipitates. Two different material systems are considered: tetragonal γ′′ precipitates in Ni-based alloys and tetragonal θ′ precipitates in Al-Cu-alloys. The shape formation and evolution of the tetragonally misfitting precipitates was investigated by means of experiments and phase-field simulations. We employed the method of invariant moments for the consistent shape quantification of precipitates obtained from the simulation as well as those obtained from the experiment. Two well-defined shape-quantities are proposed: (i) a generalized measure for the particles aspect ratio and (ii) the normalized λ2, as a measure for shape deviations from an ideal ellipse of the given aspect ratio. Considering the size dependence of the aspect ratio of γ′′ precipitates, we find good agreement between the simulation results and the experiment. Further, the precipitates’ in-plane shape is defined as the central 2D cut through the 3D particle in a plane normal to the tetragonal c-axes of the precipitate. The experimentally observed in-plane shapes of γ′′-precipitates can be quantitatively reproduced by the phase-field model.
Results of an extended TMF test program on grade P92 steel in the temperature range of 620 ◦C–300 ◦C, comprising in-phase (IP) and out-of-phase (OP) tests, partly performed with symmetric dwells at Tmax/Tmin, are presented. In contrast to previous studies, the low-strain regime is also illuminated, which approaches flexible operation in a power plant with start/stop cycles. At all strain amplitudes, the material performance is characterized by continuous cyclic softening, which is retarded in tests at lower strains but reaches similar magnitudes in the course of testing. In the investigated temperature range, the phase angle does not affect fatigue life in continuous experiments, whereas the IP condition is more detrimental in tests with dwells. Fractographic analyses indicate creep-dominated and fatigue-dominated damage for IP and OP, respectively. Analyses of the (micro)hardness distribution in the tested specimens suggest an enhanced microstructural softening in tests with dwell times for the low- but not for the high-strain regime. To rationalize the obtained fatigue data, the fracturemechanics-based DTMF concept, which was developed for TMF life assessment of ductile alloys, was applied. It is found that the DTMF parameter correlates well with the measured fatigue lives, suggesting that subcritical growth of cracks (with sizes from a few microns to a few millimeters) governs failure in the investigated range of strain amplitudes.
Cyclic mechanical performance and microstructure evolution of P92 under LCF and TMF conditions
(2019)
9-12% Cr ferritic-martensitic stainless steels are widely used as high temperature construction materials in fossil fueled power plants due to their excellent creep and oxidation resistance, but changes in electricity markets during the last two decades have considerably changed the typical working conditions of these facilities. The growing contribution of renewable energy sources in power generation forces most of these plants into flexible operation with frequent load shifts or shutdowns. These cyclic operation profiles constitute a major lifetime issue, raising the question which fundamental processes govern the reaction of ferritic-martensitic steels to cyclic load and temperature variations.
The present contribution reports on current findings obtained in a multidisciplinary project funded by German Ministry of Education and Research (BMBF) which combines cyclic mechanical and cyclic oxidation testing of different 9-12% Cr grades with detailed microstructural analyses and related micromechanical modeling.
In this contribution, an overview will be given on the results obtained in the mechanical testing programme of the project. Mechanical analyses were carried out on P91 and (mainly) P92 steel grades, particularly looking at softening phenomena and lifetimes obtained in isothermal cyclic loading (low cycle fatigue, LCF), non-isothermal cyclic loading (thermo-mechanical fatigue, TMF), and service-like combinations of fatigue and creep/relaxation periods. For this purpose, cylindrical specimens were extracted from thick-walled steam pipes, orthogonal to the pipe axis, and subjected to strain controlled cyclic loading (± 0.2 to ±0.5 % mechanical strain). Temperature intervals of TMF tests were chosen as either 300-620°C or 500-620°C, resembling so-called warm or hot start conditions of a power plant.
The test results will be presented and discussed with a focus on the impact of hold periods during testing (combined creep/relaxation-fatigue conditions) on mechanical softening, lifetime and formation of cracks. The findings will be complemented by results on the modification of the hierarchical ferritic-martensitic microstructure under different loading scenarios.
9-12% Cr ferritic-martensitic stainless steels are widely used as high temperature construction materials in power plants due to their excellent creep and oxidation resistance. The growing share of renewable energy sources in power generation forces many of these plants into more flexible operation with frequent load shifts or shutdowns. These cyclic operation profiles constitute a major lifetime issue.
The present contribution reports on current findings obtained in a multidisciplinary project which combines cyclic mechanical and cyclic oxidation testing with detailed microstructural analyses. Mechanical analyses are carried out on P92 and P91 steel grades to give an overview of softening phenomena and lifetimes obtained in isothermal cyclic loading (low cycle fatigue, LCF), non-isothermal cyclic loading (thermo-mechanical fatigue, TMF), and service-like combinations of creep and fatigue periods (creep-fatigue interaction).
Oxidation testing focuses on the grades P92 and VM12 with the intention of clarifying the impact of frequent passes through intermediate temperature levels on the kinetics of steam-side oxidation and the characteristics of the evolving oxide scales. An attempt is made to evaluate their composition, strength, integrity and adhesion after up to 250 temperature cycles. Flat coupons as well as curved tube sections are tested to assess the mutual influence of geometry on oxide scale integrity.
Complementary microstructural investigations by scanning and transmission electron microscopy plus EBSD are used for phase identification and substrate/oxide interface characterisation. The evolutions of grain size and dislocation density under different test conditions are quantified.
9-12% Cr ferritic-martensitic stainless steels are widely used as high temperature construction materials in power plants due to their excellent creep and oxidation resistance. The growing share of renewable energy sources in power generation forces many of these plants into more flexible operation with frequent load shifts or shutdowns. These cyclic operation profiles constitute a major lifetime issue. The present contribution reports on current findings obtained in a multidisciplinary project which combines cyclic mechanical and cyclic oxidation testing of different 9-12% Cr grades with detailed microstructural analyses.
Mechanical analyses are carried out on P92 and P91 steel grades to give an overview of softening phenomena and lifetimes obtained in isothermal cyclic loading (low cycle fatigue, LCF), non-isothermal cyclic loading (thermo-mechanical fatigue, TMF), and service-like combinations of creep and fatigue periods.
Complementary microstructural investigations by scanning and transmission electron microscopy plus EBSD are used for phase identification, substrate/oxide interface characterization and quantification of the microstructure evolution under cyclic conditions.
The current competitive situation on electricity markets forces power plants into cyclic operation regimes with frequent load shifts and starts/shutdowns. In the present work, the cyclic mechanical behavior of ferritic-martensitic 9-12 % Cr steels under isothermal and thermomechanical loading was investigated for the example of grade P92 material. A continuous softening was observed under all loading conditions. The introduction of hold periods to the applied cycles reduced material lifetime, with most prominent effects at technologically relevant small strain levels. The microstructural characterization reveals a coarsening of the original “martensitic” lath-type microstructure to a structure with polygonal subgrains and reduced dislocation density. The microstructural data forms the input for a physically-based modelling approach.
The aluminum alloy 2618A is applied for engine components such as radial compressor wheels which operate for long time at elevated temperatures. This results in coarsening of the hardening precipitates and degradation in mechanical properties during the long-term operation, which is not taken into account in the current lifetime prediction models due to the lack of quantitative microstructural and mechanical data. To address this issue, a quantitative investigation on the evolution of precipitates during long-term aging at 190 °C for up to 25,000 h was conducted. Detailed transmission electron microscopy (TEM) was combined with Brinell hardness measurements and thorough differential scanning calorimetry (DSC) experiments. The results showthat GPB zones and S-phase Al2CuMg grow up to < 1,000 h during which the GPB zones dissolve and S-phase precipitates form. For longer aging times, only S-phase precipitates coarsen, which can be well described using the Lifshitz–Slyozov Wagner theory of ripening. A thorough understanding of the underlying microstructural processes is a prerequisite to enable the integration of aging behavior into the established lifetime models for components manufactured from alloy 2618A.
The employment of renewable energy resources, which are naturally intermittent, for electricity generation has altered the working conditions of conventional power plants from continuous (baseload) to cyclic or flexible operation. For a long time, 9-12% Cr ferritic-martensitic stainless steels have been widely used in power plants due to their favorable characteristics such as high creep strength at high temperatures and oxidation and corrosion resistance. The components of power plants are subjected to long term cyclic loadings including fatigue and creep-fatigue at high temperatures. As ferritic-martensitic steels are known to exhibit cyclic softening when subjected to such loading scenarios, it is crucial to study the material response in such conditions. Since it is impossible to test the material behavior exactly as the operation conditions of power plants, due to technical difficulties and cost issues, it is necessary to develop physically-based material models that can predict the material behavior in more realistic situations.
In recent years, many material models have been proposed to describe the behavior of 9-12% Cr ferritic-martensitic stainless steels, which follow phenomenological or physically-based approaches. Phenomenological models provide a stress-strain relation based on empirical observations although they usually lack physical background. To alleviate this drawback and to allow for more flexibility and wider ranges of strain-rate and temperature, physically-based models are suggested. In this approach, microstructural evolution, dislocation movement and/or kinetics of plastic deformation processes are included. The physically-based models allow for a better extrapolation from the experiment results to other operation conditions and their material constants can be interpreted physically. Compared to the former approach, the material behavior can be described more accurately and flexibly and the number of material constants is less in general. In the presented work, the existing micromechanical models developed for P92 steel are compared and extended to allow for new dislocation-based strengthening/cyclic mechanisms. Their performance is assessed in the light of mechanical test data from creep-fatigue and thermo-mechanical fatigue experiments and detailed characterizations of the microstructure evolution in the fatigued material.
The formation of pores in CMSX-4 nickel based superalloys is detrimental to the service life of the material. A way to avoid the problem is to treat the superalloys under Hot Isostatic Pressing (HIP), which enables a large volume fraction of pores to be annihilated. This paper aims to understand the contribution of plastic activity related to the gliding of dislocations on the pore annihilation. Simulations based on a phase-field model of dislocation are performed and make it possible to consider the strong anisotropy of the CMSX-4 under HIP conditions in conjunction to the strong elastic heterogeneity introduced by the pore. For pores with a radius of few micrometers, it is shown that edge parts of dislocation lines that present an extra half atomic plane oriented towards the pore are stacked above and under it in the direction which is perpendicular to their slip-planes, causing an increase of the number of dislocation along the four octahedral directions of the FCC single crystal which intersect the pore center. Results are streamlined within the isotropic elastic theory of dislocations. Effects of elastic anisotropy and dislocation reactions are also investigated in order to specify what would be the dislocation configuration around a pore in CMSX-4 under HIP conditions. Notably, the elastic anisotropy is shown to significantly modify the arrangement of dislocations close to the pore equator. Simulations also allow for the characterization of pore/dislocation interactions when dislocations are involved in Low Angle Boundaries as experimentally observed.
This is the stable version of the full-notch creep test ontology (OntoFNCT) that ontologically represents the full-notch creep test. OntoFNCT has been developed in accordance with the corresponding test standard ISO 16770:2019-09 Plastics - Determination of environmental stress cracking (ESC) of polyethylene - Full-notch creep test (FNCT).
The OntoFNCT provides conceptualizations that are supposed to be valid for the description of full-notch creep tests and associated data in accordance with the corresponding test standard. By using OntoFNCT for storing full-notch creep test data, all data will be well structured and based on a common vocabulary agreed on by an expert group (generation of FAIR data) which is meant to lead to enhanced data interoperability. This comprises several data categories such as primary data, secondary data and metadata. Data will be human and machine readable. The usage of OntoFNCT facilitates data retrieval and downstream usage. Due to a close connection to the mid-level PMD core ontology (PMDco), the interoperability of full-notch creep test data is enhanced and querying in combination with other aspects and data within the broad field of materials science and engineering (MSE) is facilitated.
The class structure of OntoFNCT forms a comprehensible and semantic layer for unified storage of data generated in a full-notch creep test including the possibility to record data from analysis and re-evaluation. Furthermore, extensive metadata allows to assess data quality and reliability. Following the open world assumption, object properties are deliberately low restrictive and sparse.
Mit der Digitalisierung von Materialien und Prozessen ist ein Paradigmenwechsel in der Materialentwicklung, -gestaltung und -optimierung verbunden, welcher gleichermaßen vielseitige Möglichkeiten eröffnet und eine große Herausforderung darstellt. Insbesondere sind bei den Digitalisierungsbestrebungen die Qualitätssicherung von Prozessen und Ausgabedaten sowie die Interoperabilität zwischen Anwendungen nach FAIR-Prinzipien (Findability, Accessibility, Interoperability, Reusability) sicherzustellen. Dies umfasst die Speicherung, Verarbeitung und Abfrage von Daten in möglichst standardisierter Form, wobei entsprechend auch Normungs- und Standardisierungsgremien beteiligt werden müssen. Um der Herausforderung gerecht zu werden, Materialdaten für alle Beteiligten konsistent zu kontextualisieren, müssen alle notwendigen Informationen zum Zustand des Materials inklusive produktions- und anwendungsbedingter Veränderungen über eine einheitliche, maschinenlesbare Beschreibung verfügbar gemacht werden. Hierfür sollen Ontologien genutzt werden, da sie maschinenverständliche und -interpretierbare Wissensrepräsentationen durch semantische Konzeptualisierungen ermöglichen, die für das Datenmanagement und die Digitalisierung im Bereich der Materialwissenschaften benötigt werden.
Dieses hochaktuelle Thema der Integration und Wiederverwendung von Wissen und Daten aus Herstellung, Bearbeitung und Charakterisierung von Materialien wird in den Projekten Innovationsplattform MaterialDigital (PMD, materialdigital.de) und Materials-open-Lab (Mat-o-Lab, matolab.org) adressiert. Diese beiden unter der Beteiligung der Bundeanstalt für Materialforschung und -prüfung (BAM) durchgeführten Projekte wurden auf der 39. Vortrags- und Diskussionstagung „Werkstoffprüfung“ (2021) grundlegend vorgestellt. In dieser Präsentation sollen die Weiterentwicklungen hinsichtlich der Speicherung von Zugversuchsdaten gemäß einer normenkonformen ontologischen Repräsentation vorgestellt werden. Das umfasst den Weg von der Entwicklung einer Ontologie nach Norm, der Konvertierung von Daten aus Standardtests in das interoperable RDF-Format bis hin zur Verknüpfung von Ontologie und Daten. Letztendlich können die entsprechenden Daten in einem Triple Store abgelegt und abgefragt werden.
Die übliche Standardisierung von mechanischen Prüfverfahren im Bereich der Materialwissenschaft und Werkstofftechnik dient als solide Grundlage für die Ontologieentwicklung. Daher wurde der Zugversuch von Metallen bei Raumtemperatur nach DIN EN ISO Norm 6892-1:2019-11 als einer der ersten Anwendungsfälle in den genannten Projekten ausgewählt. Die Betrachtung und Beschreibung des Zugversuchs beinhaltet sowohl dessen ontologische Darstellung als auch eine exemplarische Datengenerierung. Die semantische Verbindung von Ontologie und Daten führt zu Interoperabilität und einer verbesserten Abfragefähigkeit.
Following the new paradigm of materials development, design and optimization, the digitalization of materials and processes is the main goal which imposes a huge challenge. In this respect, the quality assurance of processes and output data as well as the interoperability between applications following FAIR (findability, accessibility, interoperability, reusability) principles are to be ensured. This includes storage, processing and querying of data in a preferably standardized form, also addressing the incorporation of standardization bodies. To meet the challenge to contextualize material data in a way that is consistent with all stakeholders, all necessary information on the condition of the material including production and application-related changes have to be made available via a uniform, machine-readable description. For this purpose, ontologies are to be used since they allow for machine-understandable knowledge representations and semantic conceptualizations that are needed for data management and the digitalization in the field of materials science.
With respect to this currently ever-growing topic of integration and reuse of data and knowledge from synthesis, production and characterization of materials, this presentation shows the efforts taken within the project Platform MaterialDigital (PMD, materialdigital.de) to store tensile test data in accordance with a standard-compliant ontological representation. The includes the path from developing an ontology in accordance with the respective standard, converting ordinary and arbitrarily selected data gained from standard tests into the interoperable RDF format, up to connecting the ontology and data, respectively. Finally, such data can be queried from a triple store.
In the field of material science and engineering (MSE), most mechanical test methods are standardized which serves as a valid basis for ontology development. Therefore, the well-known tensile test of metals at room temperature (DIN EN ISO standard 6892-1:2019-11) was selected to be considered as one of the first use cases in PMD. This consideration within the PMD features both, the ontological representation of such a tensile test in accordance with the standard as well as exemplary data generation. The semantic connection of the ontology and data leads to interoperability and an enhanced ability of querying.
The dataset provided in this repository comprises data obtained from a series of full-notch creep tests (FNCT) performed on selected high-density polyethylene (PE-HD) materials (for further details, see section 1 Materials in this document) in accordance with the corresponding standard ISO 16770.
The FNCT is one of the mechanical testing procedures used to characterize polymer materials with respect to their environmental stress cracking (ESC) behavior. It is widely applied for PE-HD materials, that are predominantly used for pipe and container applications. It is based on the determination of the time to failure for a test specimen under constant mechanical load in a well-defined and temperature controlled liquid environment. The test device used here also allows for continuous monitoring of applied force, specimen elongation and temperature.
Relation of craze to crack length during slow crack growth phenomena in high‐density polyethylene
(2024)
The craze‐crack mechanism occurring in high‐density polyethylene (HDPE) causing slow crack growth and environmental stress cracking is investigated in detail with respect to the relation of crack length and the related craze zone. This is essential for the understanding of the resulting features of the formed fracture surface and their interpretation in the context of the transition from crack propagation to ductile shear deformation. It turns out that an already formed craze zone does not inevitably result in formation of a propagating crack, but could also undergo ductile failure. For the examination, the full notch creep test (FNCT) was employed with a subsequent advanced fracture surface analysis that was performed using various imaging techniques: light microscopy, laser scanning microscopy, scanning electron microscopy, and X‐ray micro computed tomography scan. FNCT specimens were progressively damaged for increasing durations under standard test conditions applying Arkopal, the standard surfactant solution, and biodiesel as test media were used to analyze the stepwise growth of cracks and crazes. From considerations based on well‐established fracture mechanics approaches, a theoretical correlation between the length of the actual crack and the length of the preceding craze zone was established that could be evidenced and affirmed by FNCT fracture surface analysis. Moreover, the yield strength of a HDPE material exposed to a certain medium as detected by a classic tensile test was found to be the crucial value of true stress to induce the transition from crack propagation due to the craze‐crack mechanism to shear deformation during FNCT measurements.
Highlights
- Progress of crack formation in high‐density polyethylene is analyzed by different imaging techniques
- Determined growth rates depend on distinction between craze zone and crack
- The ratio of the present crack to the anteceding craze zone is validated theoretically
- The transition from crack propagation to ductile shear deformation is identified
- An already formed craze zone may still fail by ductile mechanisms
The heat-resistant cast iron EN-GJSA-XNiSiCr35-5-2 (Ni-Resist D-5S) was investigated for its fatigue crack growth behavior at room and high temperatures. Force-controlled tests were carried out at constant temperatures (20 °C, 500 °C, 700 °C) without and with hold time and different load ratios. The crack growth behavior was also characterized under TMF loading (Tmin = 400 °C, Tmax = 700 °C) by applying IP and OP conditions and different load ratios. Three different techniques were combined to monitor crack growth: potential drop, thermography, and compliance method. The effect of the different loading conditions on the fatigue crack growth behavior will be presented and discussed.
Simulation of the θ′ Precipitation Process with Interfacial Anisotropy Effects in Al-Cu Alloys
(2021)
The effects of anisotropic interfacial properties and heterogeneous elasticity on the growth and ripening of plate-like θ′-phase (Al2Cu) in Al-1.69 at.% Cu alloy are studied. Multi-phase-field simulations are conducted and discussed in comparison with aging experiments. The precipitate/matrix interface is considered to be anisotropic in terms of its energy and mobility. We find that the additional incorporation of an anisotropic interfacial mobility in conjunction with the elastic anisotropy result in substantially larger aspect ratios of the precipitates closer to the experimental observations. The anisotropy of the interfacial energy shows comparably small effect on the precipitate’s aspect ratio but changes the interface’s shape at the rim. The effect of the chemo-mechanical coupling, i.e., the composition dependence of the elastic constants, is studied as well. We show that the inverse ripening phenomenon, recently evidenced for δ’ precipitates in Al-Li alloys (Park et al. Sci. Rep. 2019, 9, 3981), does not establish for the θ′ precipitates. This is because of the anisotropic stress fields built around the θ′ precipitates, stemming from the precipitate’s shape and the interaction among different variants of the θ′ precipitate, that disturb the chemo-mechanical effects. These results show that the chemo-mechanical effects on the precipitation ripening strongly depend on the degree of sphericity and elastic isotropy of the precipitate and matrix phases.
Lifetime aspects including fatigue failure of concrete structures were traditionally only of minor importance. Because of the growing interest in maxing out the capacities of concrete, its fatigue failure under compression has become an issue. A variety of interacting phenomena such as e.g. loss of prestress, degradation due to chemical reactions or creep and shrinkage influence the fatigue resistance. Failure due to cyclic loads is generally not instantaneous, but characterized by a steady damage accumulation. Therefore, a reliable numerical model to predict the performance of concrete over its lifetime is required, which accurately captures order effects and full three-dimensional stress states.
Many constitutive models for concrete are currently available, which are applicable for specific loading regimes, different time scales and different resolution scales.
However, a key limitation of those models is that they generally do not address issues related to fatigue on a structural level. Very few models can be found in the literature that reproduce deterioration of concrete under repeated loading-unloading cycles. This is due to the computational effort necessary to explicitly resolve every cycle which exceeds the currently available computational resources. The limitation can only be overcome by the application of multiscale methods in time.
The objective of the paper is the development of numerical methods for the simulation of concrete under fatigue loading using temporal multiscale methods.
First, a continuum damage model for concrete is developed with a focus on fatigue under compressive stresses. This includes the possibility to model stress redistributions and capture size effects. In contrast to cycle based approaches, where damage is accumulated based on the number of full stress cycles, a strain based approach is developed that can capture cyclic degradation under variable loading cycles including different amplitudes and loading frequencies. The model is designed to represent failure under static loading as a particular case of fatigue failure after a single loading cycle. As a consequence, most of the material parameters can be deduced from static tests. Only a limit set of additional constitutive parameters is required to accurately describe the evolution under fatigue loading. Another advantage of the proposed model is the possibility to directly incorporate other multi-physics effects such as creep and shrinkage or thermal loading on the constitutive level.
Second, a multiscale approach in time is presented to enable structural computations of fatigue failure with a reduced computational effort. The damage rate within the short time scale corresponding to a single cycle is computed based on a Fourier based approach. This evolution equation is then solved on the long time scale using different implicit and explicit time integration schemes. Their performance and some limitations for specific loading regimes is discussed.
Finally, the developed methods will be validated and compared to experimental data.
Les aubes de turbines à gaz utilisées en particulier pour les turboréacteurs de l’aéronautique sont élaborées par fonderie en superalliage monocristallin à base de nickel. Le procédé de fonderie, ainsi que les traitements thermique d’homogénéisation réalisés à très haute température, induisent la présence de pores au sein des pièces qui affectent les propriétés mécaniques et la durée de vie des aubes. Afin de réduire cette porosité les motoristes effectuent un traitement de compression isostatique à chaud (CIC) au cours duquel la porosité diminue par fermeture des pores. Afin de mieux comprendre les mécanismes impliqués au cours du traitement de CIC, nous avons lancé un programme de recherche dans le cadre du projet ERA-Net MICROPORE. La modélisation par champ de phase des mécanismes en jeu est présentée au cours de ce colloque. Nous présentons dans cette affiche un des volets de la caractérisation expérimentale du projet.
Des échantillons de superalliage CMSX4 sont observés après traitement de mise en solution et CIC sous 103 MPa à 1288°C pour différentes durées. Les pores présents sont caractérisés par microscopie électronique à balayage (MEB) afin de suivre l’évolution du taux de porosité au cours du traitement. Une caractérisation plus détaillée de pores partiellement refermés est menée par MEB et grâce à la diffraction des électrons rétrodiffusés (EBSD). Une vision tridimensionnelle de ces défauts est obtenue par des coupes métallographiques effectuées par découpe ionique (FIB).
Le projet ERA – Net MICROPORE est financé en Allemagne par la DFG (projects EP 136/1-1 and FE933/2-1) et en France par l’ANR (projects ANR15-MERA-000-03 and ANR15-MERA-0003-04).
Porosity in single-crystal nickel-base superalloys is removed by hot isostatic pressing (HIP) at temperatures above gamma’-solvus where the material is very soft and ductile. For example, single-crystal nickel-base superalloy CMSX-4 is HIPed at temperature 1288 °C, which is slightly higher than the gamma’-solvus temperature of this alloy equal to about 1280 °C. It is assumed that pore shrinking during HIP is mostly due to dislocation creep. Such a modelling of HIP of CMSX-4 was started in our group on the base of results of creep tests of [001] single-crystals at 1288 °C [1]. However, it was found later [2] that the alloy CMSX-4 shows very strong creep anisotropy at 1288 °C. Therefore, for calibration of the creep law, creep tests of different orientations under different stress levels are required at the HIP temperature. This was the main task of present work. Single-crystals of CMSX-4 of axial orientations [001], [011], [123] and [111] were cast by VIAM Moscow and tested by BAM Berlin under creep conditions at 1288 °C and stress levels between 4 MPa and 16 MPa. At all stress levels, the creep rate increases by an order of magnitude when changing the orientation from [001] to [111] with [011] and [123] orientations in between. Such a character of creep anisotropy corresponds to the orientation dependence of the Schmid factor for octahedral glide.
The crystal viscoplasticity model developed in [1] was improved to better represent the time induced softening observed during creep. The creep tests for different stresses and orientations as well as pore closure were simulated. The results of pore closure simulation are compared with measurements of porosity decrease during Hiping.
The phase-like behavior of grain boundaries (GBs), recently evidenced in several materials, is opening up new possibilities in the design of alloy microstructures. In this context, GB phase diagrams are contributing to a predictive description of GB segregation and (interfacial) phase changes. The influence of chemo-mechanical solute-GB interactions on the GB phase diagram remains elusive so far. This is particularly important for multi-component alloys where the elastic interactions among solute atoms, of various sizes and bonding energies, can prevail, governing a complex co-segregation phenomenon. Recently, we developed a density-based model for GB thermodynamics that intrinsically accounts for GB elasticity in pure elements. In this work, we incorporate the homogeneous and heterogeneous elastic energies associated with the solutes into the density-based framework. We derive the multi-component homogeneous elastic energy by generalizing the continuum misfitting sphere model and extend it for GBs. The density-based free energy functional directly uses bulk CALPHAD thermodynamic data. The model is applied to binary and ternary Al alloys. We reveal that the elastic energy can profoundly affect the GB solubility and segregation behavior, leading to Cu segregation in otherwise Cu-depleted Al GBs. Consequently, GB segregation transition, i.e., a jump in the GB segregation as a function of alloy composition, is revealed in Al-Cu and Al-Cu-Mg alloy systems with implications for subsequent GB precipitation in these alloys. CALPHAD-informed elasticity-incorporated GB phase diagrams enable addressing a broader range of GB phenomena in engineering multi-component alloys.
Systematic microstructure design requires reliable thermodynamic descriptions and phase diagrams of each and all microstructure elements. While such descriptions are well established for most bulk phases, thermodynamic assessment of crystal defects is greatly challenged by their individualistic aspects. In this talk, we present a density-based thermodynamic concept to describe defects based on available bulk thermodynamic data. Here dealing with grain boundaries (GBs), we apply this concept to compute GB (phase) diagram. Applications to segregation engineering of GBs in bulk and nanocrystalline alloys will be presented. We further develop this model to include the effect of elastic interactions due to atom size mismatch and obtain the corresponding GB (phase) diagram for the ternary Al-Cu-Li system.
Phase diagrams are the roadmaps for designing bulk phases. Similar to bulk, grain boundaries can possess various phases, but their phase diagrams remain largely unknown. Using a recently introduced density-based model, here we devise a strategy for computing multi-component grain boundary phase diagrams based on available bulk (CALPHAD) thermodynamic data. Fe-Mn-Cr, Fe-Mn-Ni, Fe-Mn-Co, Fe-Cr-Ni and Fe-Cr-Co alloy systems, as important ternary bases for several trending steels and high-entropy alloys, are studied. We found that despite its solute segregation enrichment, a grain boundary can have lower solubility limit than its corresponding bulk, promoting an interfacial chemical decomposition upon solute segregation. This is revealed here for the Fe-Mn-base alloy systems. The origins of this counter-intuitive feature are traced back to two effects, i.e., the magnetic ordering effect and the low cohesive energy of Mn solute element. Different aspects of interfacial phase stability and GB co-segregation in ternary alloys are investigated as well. We show that the concentration gradient energy contributions reduce segregation level but increase grain boundary solubility limit, stabilizing the GB against a chemical decomposition. Density-based grain boundary phase diagrams offer guidelines for systematic investigation of interfacial phase changes with applications to microstructure defects engineering.
For more than half a century, spinodal decomposition has been a key phenomenon in considering the formation of secondary phases in alloys. The most prominent aspect of the spinodal phenomenon is the lack of an energy barrier on its transformation pathway, offering an alternative to the nucleation and growth mechanism. The classical description of spinodal decomposition often neglects the influence of defects, such as grain boundaries, on the transformation because the innate ability for like-atoms to cluster is assumed to lead the process. Nevertheless, in nanocrystalline alloys, with a high population of grain boundaries with diverse characters, the structurally heterogeneous landscape can greatly influence the chemical decomposition behavior. Combining atom-probe tomography, precession electron diffraction and density-based phase-field simulations, we address how grain boundaries contribute to the temporal evolution of chemical decomposition within the miscibility gap of a Pt-Au nanocrystalline system. We found that grain boundaries can actually have their own miscibility gaps profoundly altering the spinodal decomposition in nanocrystalline alloys. A complex realm of multiple interfacial states, ranging from competitive grain boundary segregation to barrier-free low-dimensional interfacial decomposition, occurs with a dependency upon the grain boundary character.
Grain boundaries (GBs) are planar lattice defects that govern the properties of many types of polycrystalline materials. Hence, their structures have been investigated in great detail. However, much less is known about their chemical features, owing to the experimental difficulties to probe these features at the atomic length scale inside bulk material specimens. Atom probe tomography (APT) is a tool capable of accomplishing this task, with an ability to quantify chemical characteristics at near-atomic scale. Using APT data sets, we present here a machine-learning-based approach for the automated quantification of chemical features of GBs. We trained a convolutional neural network (CNN) using twenty thousand synthesized images of grain interiors, GBs, or triple junctions. Such a trained CNN automatically detects the locations of GBs from APT data. Those GBs are then subjected to compositional mapping and analysis, including revealing their in-plane chemical decoration patterns. We applied this approach to experimentally obtained APT data sets pertaining to three case studies, namely, Ni-P, Pt-Au, and Al-Zn-Mg-Cu alloys. In the first case, we extracted GB specific segregation features as a function of misorientation and coincidence site lattice character. Secondly, we revealed interfacial excesses and in-plane chemical features that could not have been found by standard compositional analyses. Lastly, we tracked the temporal evolution of chemical decoration from early-stage solute GB segregation in the dilute limit to interfacial phase separation, characterized by the evolution of complex composition patterns. This machine-learning-based approach provides quantitative, unbiased, and automated access to GB chemical analyses, serving as an enabling tool for new discoveries related to interface thermodynamics, kinetics, and the associated chemistry-structure-property relations.