5 Werkstofftechnik
Filtern
Dokumenttyp
- Zeitschriftenartikel (44)
- Vortrag (44)
- Beitrag zu einem Tagungsband (12)
- Posterpräsentation (4)
- Forschungsdatensatz (3)
- Buchkapitel (2)
- Dissertation (1)
- Forschungsbericht (1)
Sprache
- Englisch (97)
- Deutsch (13)
- Französisch (1)
Schlagworte
- CALPHAD (9)
- Creep (9)
- Nickel-base superalloys (8)
- Fatigue (7)
- LCF (7)
- Crack growth (6)
- Crystal plasticity (6)
- Thermomechanical fatigue (6)
- Additive manufacturing (5)
- Anisotropy (5)
Organisationseinheit der BAM
- 5.5 Materialmodellierung (111) (entfernen)
Paper des Monats
- ja (2)
Eingeladener Vortrag
- nein (44)
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.
The current trend towards cyclic, “flexible” operation of fossil-fueled power plants constitutes a major issue regarding lifetime and operational safety of the respective installations and their components, as was outlined in our complementary contribution (Part 1). The present contribution reports on the investigation of the microstructure evolution in cyclically loaded ferritic-martensitic steels and its representation in a physically-based micromechanical model.
For this purpose, specimens of P92 steel grade from the mechanical test programme outlined in our companion contribution (Part 1) were analyzed by scanning electron microscopy (SEM), including backscattered diffraction (EBSD) mapping, and transmission electron microscopy (TEM). A novel method was implemented to improve angular resolution of EBSD scans. Additionally, a correlative microscopy approach was developed and used to correlate EBSD and TEM measurements on the same locations of thick regions of electron transparent specimens. By applying these techniques, a detailed quantitative microstructure description of the as-received material condition, namely in terms of subgrain morphology and dislocation density/distributions, was established. Comparisons of as-received and cyclically loaded conditions from tests interrupted at different stages of lifetime indicate a rapid redistribution of in-grain dislocations with a strong interaction between mobile dislocations and low angle grain boundaries (LABs).
The proposed micromechanical model is formulated in a viscoplastic self-consistent (VPSC) scheme, which is a mean-field approach that allows us to include the crystal details at the level of slip systems while avoiding the considerable computational costs of full-field approaches (such as the classical crystal plasticity finite element analysis). Being physically-based, the model uses dislocation densities and includes the interaction between dislocations, e.g. annihilation of mobile dislocations, and evolution of microstructure, e.g. the grain coarsening. Particularly, the constitutive laws for dislocation evolution and interaction between dislocations and low angle boundaries are calibrated based on two-dimensional discrete dislocation dynamic (2D DDD) simulations, which are performed at a micro-/meso-scale. The results of the beforementioned EBSD experiments are considered as a direct input, involving e.g. the amount of geometrically necessary dislocations, average misorientations and grain characteristics.
The embrittlement of metallic alloys by liquid metals leads to catastrophic material failure and severely impacts their structural integrity. The weakening of grain boundaries by the ingress of liquid metal and preceding segregation in the solid are thought to promote early fracture. However, the potential of balancing between the segregation of cohesion-enhancing interstitial solutes and embrittling elements inducing grain boundary decohesion is not understood. Here, we unveil the mechanisms of how boron segregation mitigates the detrimental effects of the prime embrittler, zinc, in a Σ5 [0 0 1] tilt grain boundary in α −Fe (4 at.% Al). Zinc forms nanoscale segregation patterns inducing structurally and compositionally complex grain boundary states. Ab-initio simulations reveal that boron hinders zinc segregation and compensates for the zinc induced loss in grain boundary cohesion. Our work sheds new light on how interstitial solutes intimately modify grain boundaries, thereby opening pathways to use them as dopants for preventing disastrous material failure.
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.
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.
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.
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.
The large surface-to-volume ratio of nanoparticles is understood to be the source of many interesting phenomena. The melting temperature of nanoparticles is shown to dramatically reduce compared to bulk material. Yet, at temperatures below this reduced melting point, a liquid-like atomic arrangement on the surface of nanoparticles is still anticipated to influence its properties. To understand such surface effects, here, we study the coalescence of Au nanoparticles of various sizes using molecular dynamics simulations. Analysis of the potential energy and Lindemann index distribution across the nanoparticles reveals that high-energy, high-mobility surface atoms can enable the coalescence of nanoparticles at temperatures much lower than their corresponding melting point. The smaller the nanoparticles, the larger the difference between their melting and coalescence temperatures. For small enough particles and/or elevated enough temperatures, we found that the coalescence leads to a melting transition of the two nominally solid nanoparticles, here discussed in relation to the heat released due to the surface reduction upon the coalescence and the size dependence of latent heat. Such discontinuous melting transitions can lead to abrupt changes in the properties of nanoparticles, important for their applications at intermediate temperatures.
A main source of current challenges in materials science and engineering is the ever-increasing complexity in materials chemistry and processing resulting in complex microstructures, making the assessment of process-microstructure-property-performance relations difficult, even unmanageable. Here the computational materials science is facing the same situation. In this talk, I share a viewpoint that the complexities in chemistry, processing and microstructures can be circumscribed by integrating existing knowledges of bulk thermodynamics and kinetics to the unknown thermodynamics and kinetics of microstructure elements. To this end, I discuss several successful examples on grain boundary segregation engineering how this scientific advance can be conducted. A roadmap is proposed, beginning to form on generalizing the concept of phase diagrams.
Engineering grain boundaries demands a quantitative description of both their segregation and specific phase behavior. Recently I have proposed a density-based model for grain boundary thermodynamics that enables CALPHAD integrated derivation of grain boundary phase diagrams, broadly applied now in studying various alloys. Combining this model with experimental investigations, in this talk, new aspects of interfacial segregation and phase transformation revealed in polycrystalline alloys are discussed. The effect of elastic interaction on grain boundary phase behavior is incorporated. We consider Al alloys and novel high-entropy alloys and discuss a general strategy for grain boundary engineering.
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.
Several recent studies have revealed that the phase decomposition in medium- and high-entropy alloys is triggered by solute segregation at grain boundaries. The chemically complex nature of these alloys already clues that such segregation in materials should be something more than the interaction between solutes and grain boundary, but also affected by the solute-solute interactions during the segregation. In fact, experiments indicate the significance of such solute-solute interaction. Yet, the mechanisms of co-segregation in multi-component alloys are rather challenging to explore and rarely studied quantitatively. Recently, a CALPHAD-integrated density-based phase-field model has been proposed for studying grain boundary phenomena. Several applications of this model have shown its expansive capacity for investigating grain boundary segregation and phase behavior. In this talk, the applications of this model to studying co-segregation phenomena and grain boundary phase diagrams in multi-component alloys are presented ––the term ‘co-segregation’ here refers to co-evolution and any mutual interplay among the solutes and grain boundary during the segregation. In particular, the iron-based ternary and quaternary alloys are discussed. I show how a grain boundary may have its own miscibility gap and how this immiscibility can influence the co-segregation behavior. As an alternative to grain boundary phase diagrams, rather suited for binary and ternary alloys, a new concept of co-segregation maps for screening and segregation design in multi-component alloys is presented. Applying the co-segregation maps, the nonlinear Mn and Cr co-segregation in Fe-Co-Mn-Cr is discussed.
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.
A giant Zn segregation transition is revealed using CALPHAD-integrated density-based modeling of segregation into Fe grain boundaries (GBs). The results show that above a threshold of only a few atomic percent Zn in the alloy, a substantial amount of up to 60 at.% Zn can segregate to the GB. We found that the amount of segregation abruptly increases with decreasing temperature, while the Zn content in the alloy required for triggering the segregation transition decreases. Direct evidence of the Zn segregation transition is obtained using high-resolution scanning transmission electron microscopy. Base on the model, we trace the origin of the segregation transition back to the low cohesive energy of Zn and a miscibility gap in Fe-Zn GB, arising from the magnetic ordering effect, which is confirmed by ab-initio calculations. We also show that the massive Zn segregation resulting from the segregation transition greatly assists with liquid wetting and reduces the work of separation along the GB. The current predictions suggest that control over Zn segregation, by both alloy design and optimizing the galvanization and welding processes, may offer preventive strategies against liquid metal embrittlement.
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.
Recent advances in Deep Learning (DL) have significantly impacted the field of materials modelling. DL tools have been recently considered as promising tools to address the complex relationships among processing, microstructure and property of materials. The thermal stability of polycrystalline materials is a highly interesting and complex problem that could be addressed using DL techniques. The Multi-Phase-Field (MPF) method has emerged as a powerful tool for addressing grain growth phenomena from multiple perspectives. Unlike sharp-interface based methods, the MPF approach bypasses the need for detailed information on individual grains. In this work, we use DL to address issues related to MPF simulations of grain growth including numerical efficiency, computing speed, and resource consumption. This presents specific challenges for high-performance computing (HPC) due to the large datasets and complex computations required by both MPF and DL methods. We study various 3D microstructure settings with the goal of accelerating the simulation process while exploring different physical effects. In particular, the impact of grain boundary and triple junction energies on grain growth are to be investigated. The results will be presented in terms of the evolving size and shape distribution of the grains.
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.
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.
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.
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.
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.
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.
Extended X-ray absorption fine structure (EXAFS) conducted on an equiatomic MoNbTaW bcc medium-entropy alloy that was annealed at 2273 K reveals unexpectedly small 1st and 2nd shell element-specific lattice distortions. An experimental size-mismatch parameter, δexp, is determined to be ca. 50% lower than the corresponding calculated value. Around W, short-range order (SRO) preferring 4d elements in the 1st and 2nd shells persists. A Nb-W ordering is found, which is reminiscent of ordering emerging at lower temperatures in the B2(Mo,W;Ta,Nb)- and B32(Nb,W)-phases. With high-temperature ordering preferences in fcc also foreshadowing low-temperature phase, these findings suggest a general feature of high-temperature SRO.
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.
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.
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.
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.
Inconel 718 (IN718) is a traditional age-hardenable nickel-based alloy that has been increasingly processed by additive manufacturing (AM) in recent years. In the as-solidified condition, IN718 exhibits chemical segregation and the undesired Laves phase, requiring a solution annealing (SA) prior to aging. The material produced by AM does not respond to the established thermal routines in the same way as conventionally produced IN718, and there is still no consensus on which routine yields optimal results. This work aims to provide a fundamental understanding of the heat treatment (HT) response by continuously monitoring the microstructural evolution during SA via time-resolved synchrotron x-ray diffraction, complemented by ex-situ scanning electron microscopy (SEM). The samples were produced by laser powder bed fusion to a geometry of 10x20x90 mm³, from which Ø1x5 mm³ cylindric specimens were extracted. Two different scanning strategies – incremental 67° rotations, Rot, and alternating 0°/67° tracks, Alt – were used, leading to two different as-built conditions. 1-hour SAs were carried out in the beamline ID22 of the ESRF at 50 KeV. Two SA temperatures, SA1 = 1020 °C, and SA2 = 1080 °C were tested for each scanning strategy. Data were processed using the software PDIndexer. In the as-built state, all samples showed typical subgrain columnar cell structures with predominant Nb/Mo segregation and Laves phase at the cell walls, as seen by SEM. The Alt scan induced higher intensity on the Laves peaks than the Rot scan, suggesting a greater content of Laves. Chemical homogenization in the SA was largely achieved during the heating ramp (Fig. 1). SA2 eliminated the Laves peaks just before reaching 1080 °C, and mitigated differences between Rot and Alt samples. On the other hand, SA1 induced the precipitation of the generally detrimental δ phase, also observed by SEM. Furthermore, the Rot scan showed higher δ peak intensities than the Alt scan, indicating a higher content of δ in the latter. No signs of recrystallization were observed in any of the investigated SAs.
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.