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.