5.5 Materialmodellierung
Filtern
Dokumenttyp
- Vortrag (45)
- Zeitschriftenartikel (44)
- Beitrag zu einem Tagungsband (12)
- Posterpräsentation (4)
- Forschungsdatensatz (3)
- Buchkapitel (2)
- Dissertation (1)
- Forschungsbericht (1)
Sprache
- Englisch (98)
- Deutsch (13)
- Französisch (1)
Schlagworte
- CALPHAD (9)
- Creep (9)
- Nickel-base superalloys (8)
- Fatigue (7)
- LCF (7)
- Thermomechanical fatigue (7)
- Additive manufacturing (6)
- Crack growth (6)
- Crystal plasticity (6)
- Anisotropy (5)
Organisationseinheit der BAM
- 5 Werkstofftechnik (112)
- 5.5 Materialmodellierung (112)
- 5.2 Metallische Hochtemperaturwerkstoffe (78)
- 5.1 Mikrostruktur Design und Degradation (15)
- 7 Bauwerkssicherheit (11)
- 3 Gefahrgutumschließungen; Energiespeicher (7)
- 9 Komponentensicherheit (7)
- 3.0 Abteilungsleitung und andere (5)
- 8 Zerstörungsfreie Prüfung (5)
- 8.5 Röntgenbildgebung (5)
Paper des Monats
- ja (3)
Eingeladener Vortrag
- nein (45)
Lifetime aspects including fatigue failure of concrete structures were traditionally only of minor importance. Because of the growing interest in maxing out the capacities of concrete, its fatigue failure under compression has become an issue. A variety of interacting phenomena such as e.g. loss of prestress, degradation due to chemical reactions or creep and shrinkage influence the fatigue resistance. Failure due to cyclic loads is generally not instantaneous, but characterized by a steady damage accumulation. Therefore, a reliable numerical model to predict the performance of concrete over its lifetime is required, which accurately captures order effects and full three-dimensional stress states.
Many constitutive models for concrete are currently available, which are applicable for specific loading regimes, different time scales and different resolution scales.
However, a key limitation of those models is that they generally do not address issues related to fatigue on a structural level. Very few models can be found in the literature that reproduce deterioration of concrete under repeated loading-unloading cycles. This is due to the computational effort necessary to explicitly resolve every cycle which exceeds the currently available computational resources. The limitation can only be overcome by the application of multiscale methods in time.
The objective of the paper is the development of numerical methods for the simulation of concrete under fatigue loading using temporal multiscale methods.
First, a continuum damage model for concrete is developed with a focus on fatigue under compressive stresses [1]. This includes the possibility to model stress redistributions and capture size effects. In contrast to cycle based approaches, where damage is accumulated based on the number of full stress cycles, a strain based approach is developed that can capture cyclic degradation under variable loading cycles including different amplitudes and loading frequencies. The model is designed to represent failure under static loading as a particular case of fatigue failure after a single loading cycle. As a consequence, most of the material parameters can be deduced from static tests. Only a limit set of additional constitutive parameters is required to accurately describe the evolution under fatigue loading. Another advantage of the proposed model is the possibility to directly incorporate other multi-physics effects such as creep and shrinkage or thermal loading on the constitutive level.
Second, a multiscale approach in time is presented to enable structural computations of fatigue failure with a reduced computational effort. The damage rate within the short time scale corresponding to a single cycle is computed based on a Fourier based approach [2]. This evolution equation is then solved on the long time scale using different implicit and explicit time integration schemes. Their performance and some limitations for specific loading regimes is discussed.
Finally, the developed methods will be validated and compared to experimental data.
[1] Vitaliy Kindrachuk, Marc Thiele, Jörg F. Unger. Constitutive modeling of creep-fatigue interaction for normal strength concrete under compression, International Journal of Fatigue, 78:81-94, 2015
[2] Vitaliy Kindrachuk, Jörg F. Unger. A Fourier transformation-based temporal integration scheme for viscoplastic solids subjected to fatigue deterioration, International Journal of Fatigue, 100:215-228, 2017
Lifetime aspects including fatigue failure of concrete structures were traditionally only of minor importance. Because of the growing interest in maxing out the capacities of concrete, its fatigue failure under compression has become an issue. A variety of interacting phenomena such as e.g. loss of prestress, degradation due to chemical reactions or creep and shrinkage influence the fatigue resistance. Failure due to cyclic loads is generally not instantaneous, but characterized by a steady damage accumulation. Therefore, a reliable numerical model to predict the performance of concrete over its lifetime is required, which accurately captures order effects and full three-dimensional stress states.
Many constitutive models for concrete are currently available, which are applicable for specific loading regimes, different time scales and different resolution scales.
However, a key limitation of those models is that they generally do not address issues related to fatigue on a structural level. Very few models can be found in the literature that reproduce deterioration of concrete under repeated loading-unloading cycles. This is due to the computational effort necessary to explicitly resolve every cycle which exceeds the currently available computational resources. The limitation can only be overcome by the application of multiscale methods in time.
The objective of the paper is the development of numerical methods for the simulation of concrete under fatigue loading using temporal multiscale methods.
First, a continuum damage model for concrete is developed with a focus on fatigue under compressive stresses. This includes the possibility to model stress redistributions and capture size effects. In contrast to cycle based approaches, where damage is accumulated based on the number of full stress cycles, a strain based approach is developed that can capture cyclic degradation under variable loading cycles including different amplitudes and loading frequencies. The model is designed to represent failure under static loading as a particular case of fatigue failure after a single loading cycle. As a consequence, most of the material parameters can be deduced from static tests. Only a limit set of additional constitutive parameters is required to accurately describe the evolution under fatigue loading. Another advantage of the proposed model is the possibility to directly incorporate other multi-physics effects such as creep and shrinkage or thermal loading on the constitutive level.
Second, a multiscale approach in time is presented to enable structural computations of fatigue failure with a reduced computational effort. The damage rate within the short time scale corresponding to a single cycle is computed based on a Fourier based approach. This evolution equation is then solved on the long time scale using different implicit and explicit time integration schemes. Their performance and some limitations for specific loading regimes is discussed.
Finally, the developed methods will be validated and compared to experimental data.
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.
A key limitation of the most constitutive models that reproduce a Degradation of quasi-brittle materials is that they generally do not address issues related to fatigue. One reason is the huge computational costs to resolve each load cycle on the structural level. The goal of this paper is the development of a temporal Integration scheme, which significantly increases the computational efficiency of the finite element method in comparison to conventional temporal integrations.
The essential constituent of the fatigue model is an implicit gradient-enhanced formulation of the damage rate. The evolution of the field variables is computed as amultiscale Fourier series in time.On a microchronological scale attributed to single cycles, the initial boundary value problem is approximated by linear BVPs with respect to the Fourier coefficients. Using the adaptive cycle jump concept, the obtained damage rates are transferred to a coarsermacrochronological scale associated with the duration of material deterioration. The performance of the developedmethod is hence improved due to an efficient numerical treatment of the microchronological problem in combination with the cycle jump technique on the macrochronological scale. Validation examples demonstrate the convergence of the obtained solutions to the reference simulations while significantly reducing the computational costs.
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.
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 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 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.
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.
A phenomenological criterion for an optical assessment of PE-HD fracture surfaces obtained from FNCT
(2021)
The full-notch creep test (FNCT) is a common test method to evaluate the environmental stress cracking (ESC) behavior of high-density polyethylene (PE-HD), e.g. for container materials. The test procedure as specified in ISO 16770 provides a comparative measure of the resistance against ESC using the time to failure of PE-HD specimens under constant mechanical load in a well-defined liquid test environment. Since the craze-crack damage mechanism underlying the ESC phenomenon is associated with brittle failure, the occurrence of a predominantly brittle fracture surface is a prerequisite to consider an FNCT measurement as representative for ESC, i.e. a time to failure dominated by craze-crack propagation.
The craze-crack propagation continuously reduces the effective residual cross-sectional area of the specimen during the test, which results in a corresponding increase of the effective mechanical stress. Thus, a transition to ductile shear deformation is inevitable at later stages of the test, leading usually to a pronounced central ligament.
Therefore, an optical evaluation of FNCT fracture surfaces concerning their brittleness is essential. An enhanced imaging analysis of FNCT fracture surfaces enables a detailed assessment of craze-crack Propagation during ESC. In this study, laser scanning microscopy (LSM) was employed to evaluate whether FNCT fracture surfaces are representative with respect to craze-crack propagation and ESC. Based on LSM height data, a phenomenological criterion is proposed to assess the validity of distinct FNCT measurements. This criterion is
supposed to facilitate a quick evaluation of FNCT results in practical routine testing. Its applicability is verified on a sample basis for seven different commercial PE-HD container materials.
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.
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.
The full-notch creep test (FNCT) is a common method to evaluate the environmental stress cracking (ESC) behavior of high-density polyethylene (PE-HD) container materials . The test procedure as specified in ISO 16770 provides a comparative measure of the resistance against ESC using the time to failure of specimens mechanically loaded in a well-defined liquid environment. Since the craze-crack damage mechanism underlying the ESC process is associated with brittle failure, the occurrence of globally brittle fracture surfaces is a prerequisite to consider an FNCT measurement as representative for ESC . Therefore, an optical evaluation of FNCT fracture surfaces concerning their brittleness is essential. Due to the experimental setup, an inevitable increase of the true mechanical stress and the associated appearance of small ductile parts on fracture surfaces is induced in any case. Hence, an FNCT experiment is considered as 'valid', if the corresponding fracture surface is predominantly brittle . Based on laser scanning microscopy (LSM) height data of FNCT fracture surfaces , a universal and easy-to-use phenomenological criterion was developed to assess the validity of distinct FNCT experiments. This criterion is supposed to facilitate a quick evaluation of FNCT results in practical routine testing.
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.
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.
Druckbehälter werden vor ihrer Inbetriebnahme obligatorisch einer hydraulischen Prüfung bis zum sog. Prüfdruck belastet. Damit sollen fehlerhaft gefertigte Individuen erkannt und ausgesondert werden. Mit dieser Prüfung werden aber auch diejenigen Einzelbehälter ausgesondert, die ohne einem nachvollziehbaren Fertigungsfehler, nur auf Basis der statistischen Streuung ggf. eine Festigkeit unterhalb des Prüfdrucks aufweisen.
Da der Prüfdruck im Fall der Wasserstoffspeicher mindestens 20% über dem maximal zulässigem Betriebsdruck liegt, wäre damit ein Versagen im Betrieb ausgenommen, gäbe es die verschiedenen Alterungs- bzw. Ermüdungseffekt nicht.
Um den Einfluss der erstmaligen Prüfung auf eine angenommene ausfallfreie Zeit bewerten zu können, wird hier die Alterung und die erstmalige Prüfung in Kombination analysiert. Hierzu wird auf Basis der Restfestigkeitsdaten aus einer umfangreichen Prüfkampagne ein Ansatz für die Beschreibung der Alterung ermittelt und auf andere Lastzustände übertragen.
So kann für das Baumuster, das der o.g. Prüfkampagne zugrunde lag, mithilfe einer Monte-Carlo-Simulation gezeigt werden, dass mindestens 10.000 LW vor einem ersten Alterungsversagen mit 1 aus 1 Mio. In Abhängigkeit der Anzahl der Individuen, die bei der erstmaligen Prüfung versagen, und der Streuung der gesamten Population, kann diese Mindestfestigkeit auf realistisch 50.000 Lastwechsel steigen.
Damit ist gezeigt, dass aufgrund der erstmaligen Prüfung eine Mindestfestigkeit erzeugt wird, die stark vom Alterungsverhalten abhängt. Im Ergebnis heißt dies, dass in den Vorschriften aktuell die Bedeutung des Mindestberstdrucks in der Baumusterprüfung überschätzt und die möglichst betriebsbegleitende Erfassung der Alterung unterschätzt sind.
Druckbehälter werden vor ihrer Inbetriebnahme obligatorisch einer hydraulischen Prüfung bis zum sog. Prüfdruck belastet. Damit sollen fehlerhaft gefertigte Individuen erkannt und ausgesondert werden. Mit dieser Prüfung werden aber auch diejenigen Einzelbehälter ausgesondert, die ohne einem nachvollziehbaren Fertigungsfehler, nur auf Basis der statistischen Streuung ggf. eine Festigkeit unterhalb des Prüfdrucks aufweisen.
Da der Prüfdruck im Fall der Wasserstoffspeicher mindestens 20% über dem maximal zulässigem Betriebsdruck liegt, wäre damit ein Versagen im Betrieb ausgenommen, gäbe es die verschiedenen Alterungs- bzw. Ermüdungseffekt nicht.
Um den Einfluss der erstmaligen Prüfung auf eine angenommene ausfallfreie Zeit bewerten zu können, wird hier die Alterung und die erstmalige Prüfung in Kombination analysiert. Hierzu wird auf Basis der Restfestigkeitsdaten aus einer umfangreichen Prüfkampagne ein Ansatz für die Beschreibung der Alterung ermittelt und auf andere Lastzustände übertragen.
So kann für das Baumuster, das der o.g. Prüfkampagne zugrunde lag, mithilfe einer Monte-Carlo-Simulation gezeigt werden, dass mindestens 10.000 LW vor einem ersten Alterungsversagen mit 1 aus 1 Mio. In Abhängigkeit der Anzahl der Individuen, die bei der erstmaligen Prüfung versagen, und der Streuung der gesamten Population, kann diese Mindestfestigkeit auf realistisch 50.000 Lastwechsel steigen.
Damit ist gezeigt, dass aufgrund der erstmaligen Prüfung eine Mindestfestigkeit erzeugt wird, die stark vom Alterungsverhalten abhängt. Im Ergebnis heißt dies, dass in den Vorschriften aktuell die Bedeutung des Mindestberstdrucks in der Baumusterprüfung überschätzt und die möglichst betriebsbegleitende Erfassung der Alterung unterschätzt sind.
Aero-engine turbine disks are safety-relevant components which are operated under high thermal and mechanical stress conditions. The actual part qualification and certification procedures make use of spin-tests conducted on production-similar disks. While these tests provide, on the one hand, a reliable definition of the critical conditions for real components, on the other hand they represent a relevant cost item for engine manufacturers. The aim of this work is to present two alternative burst speed assessment methods under development based on the Failure Assessment Diagram (FAD) and a global stability criterion, respectively. In the scope of the fracture mechanics assessment, the failure modes hoop-burst and rim-peeling are investigated with semicircular surface cracks modelled at the critical regions on the turbine disk. The comparison of the predicted critical rotational speed shows good agreement between the assessment methods.
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.
Along with the desire for developing novel multi-principal element alloys, also known as high-entropy alloys, the concern about their safe application is also increasingly growing. This relates to the alloys’ phase stability, in particular, the control required over unexpected phase decompositions resulting from solute segregation at grain boundaries. Yet, the mechanisms of co-segregation and grain boundary phase decomposition in multi-component alloys are rather challenging to explore. In fact, quantitative investigation of grain boundary behaviors is mostly conducted for binary and a few ternary alloys. In this work, we apply the recently introduced CALPHAD-integrated density-based formalism [RSC Advances 10 (2020) 26728-26741] for considering co-segregation phenomena in alloys with an arbitrary number of components —the term ‘co-segregation’ here refers to co-evolution and any mutual interplay among the solute atoms during their interaction with a grain boundary. Quaternary Fe-Co-Mn-Cr alloy system is studied. We present two major advances beyond previous results: First, a co-segregation-induced multi-component grain boundary spinodal decomposition is quantitatively simulated for the first time. We found that in addition to its low cohesive energy and asymmetrical mixing enthalpy due to magnetic ordering, Mn plays a leading role in triggering interfacial phase decomposition by having a relatively large, concentration-dependent atomic mobility. Second, as an alternative to grain boundary phase diagrams proposed for binary and ternary alloys, we introduce the concept of co-segregation maps for grain boundary segregation screening and design in multi-component alloys. Applying the co-segregation maps, the nonlinear Mn and Cr co-segregation are discussed. Depicted on the alloying composition and phase space, the co-segregation maps enable the required insights to guide a safer, more controlled design of high-entropy alloys.
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.
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.
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.
The durability of concrete structures and its performance over the lifetime is strongly influenced by many interacting phenomena such as e.g. mechanical degradation due to fatigue loading, loss of prestress, degradation due to chemical reactions or creep and shrinkage. Failure due to cyclic loading is generally not instantaneous, but characterized by a steady damage accumulation.
Many constitutive models for concrete are currently available, which are applicable for specific loading regimes, different time scales and different resolution scales. A key limitation is that the models often do not address issues related to fatigue on a structural level. Very few models can be found in the literature that reproduce deterioration of concrete under repeated loading-unloading cycles.
The objective of this paper is the presentation of numerical methods for the simulation of concrete under fatigue loading using a temporal multiscale method.
First, a continuum damage model for concrete is developed with a focus on fatigue under compressive stresses. This includes the possibility to model stress redistributions and capture size effects. In contrast to cycle based approaches, where damage is accumulated based on the number of full stress cycles, a strain based approach is developed that can capture cyclic degradation under variable loading cycles including different amplitudes and loading frequencies. Second, a multiscale approach in time is presented to enable structural computations of fatigue failure with a reduced computational effort. The damage rate within the short time scale corresponding to a single cycle is computed based on a Fourier based approach. This evolution equation is then solved on the long time scale using different time integration schemes.
A density-based phase field model is developed where the free energy functional is explicitly linked with molecular dynamics and is referred to as the Molecular Phase Field Method (MoPF). MoPF simulations involve expressing interatomic potentials in terms of density to form a density based free energy functional. Inputs to this functional are taken from atomistics such that the phase field density profile matches the corresponding density profile from atomistic simulations. We analyze our results by comparing the MoPF calculated excess interfacial energies with excess interfacial energies calculated using molecular dynamics associated with several nickel grain boundaries. Additionally, a comparison is made between our results and the interfacial energies of a \Sigma7 boundary across a variety of FCC systems simulated using density functional theory. The MoPF method is able to successfully predict grain boundary free energy trends between grain boundary and material types offering an atomistically informed mesoscale formulation for studying grain boundary physics.
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.
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.
Porosity in single-crystal nickel-base superalloys is removed by hot isostatic pressing (HIP) at temperatures above gamma’-solvus where the material is very soft and ductile. For example, single-crystal nickel-base superalloy CMSX-4 is HIPed at temperature 1288 °C, which is slightly higher than the gamma’-solvus temperature of this alloy equal to about 1280 °C. It is assumed that pore shrinking during HIP is mostly due to dislocation creep. Such a modelling of HIP of CMSX-4 was started in our group on the base of results of creep tests of [001] single-crystals at 1288 °C [1]. However, it was found later [2] that the alloy CMSX-4 shows very strong creep anisotropy at 1288 °C. Therefore, for calibration of the creep law, creep tests of different orientations under different stress levels are required at the HIP temperature. This was the main task of present work. Single-crystals of CMSX-4 of axial orientations [001], [011], [123] and [111] were cast by VIAM Moscow and tested by BAM Berlin under creep conditions at 1288 °C and stress levels between 4 MPa and 16 MPa. At all stress levels, the creep rate increases by an order of magnitude when changing the orientation from [001] to [111] with [011] and [123] orientations in between. Such a character of creep anisotropy corresponds to the orientation dependence of the Schmid factor for octahedral glide.
The crystal viscoplasticity model developed in [1] was improved to better represent the time induced softening observed during creep. The creep tests for different stresses and orientations as well as pore closure were simulated. The results of pore closure simulation are compared with measurements of porosity decrease during Hiping.
The 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.
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.
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.
9-12% Cr ferritic-martensitic stainless steels are widely used as high temperature construction materials in power plants due to their excellent creep and oxidation resistance. The growing share of renewable energy sources in power generation forces many of these plants into more flexible operation with frequent load shifts or shutdowns. These cyclic operation profiles constitute a major lifetime issue. The present contribution reports on current findings obtained in a multidisciplinary project which combines cyclic mechanical and cyclic oxidation testing of different 9-12% Cr grades with detailed microstructural analyses.
Mechanical analyses are carried out on P92 and P91 steel grades to give an overview of softening phenomena and lifetimes obtained in isothermal cyclic loading (low cycle fatigue, LCF), non-isothermal cyclic loading (thermo-mechanical fatigue, TMF), and service-like combinations of creep and fatigue periods.
Complementary microstructural investigations by scanning and transmission electron microscopy plus EBSD are used for phase identification, substrate/oxide interface characterization and quantification of the microstructure evolution under cyclic conditions.
The current competitive situation on electricity markets forces power plants into cyclic operation regimes with frequent load shifts and starts/shutdowns. In the present work, the cyclic mechanical behavior of ferritic-martensitic 9-12 % Cr steels under isothermal and thermomechanical loading was investigated for the example of grade P92 material. A continuous softening was observed under all loading conditions. The introduction of hold periods to the applied cycles reduced material lifetime, with most prominent effects at technologically relevant small strain levels. The microstructural characterization reveals a coarsening of the original “martensitic” lath-type microstructure to a structure with polygonal subgrains and reduced dislocation density. The microstructural data forms the input for a physically-based modelling approach.
Cyclic mechanical performance and microstructure evolution of P92 under LCF and TMF conditions
(2019)
9-12% Cr ferritic-martensitic stainless steels are widely used as high temperature construction materials in fossil fueled power plants due to their excellent creep and oxidation resistance, but changes in electricity markets during the last two decades have considerably changed the typical working conditions of these facilities. The growing contribution of renewable energy sources in power generation forces most of these plants into flexible operation with frequent load shifts or shutdowns. These cyclic operation profiles constitute a major lifetime issue, raising the question which fundamental processes govern the reaction of ferritic-martensitic steels to cyclic load and temperature variations.
The present contribution reports on current findings obtained in a multidisciplinary project funded by German Ministry of Education and Research (BMBF) which combines cyclic mechanical and cyclic oxidation testing of different 9-12% Cr grades with detailed microstructural analyses and related micromechanical modeling.
In this contribution, an overview will be given on the results obtained in the mechanical testing programme of the project. Mechanical analyses were carried out on P91 and (mainly) P92 steel grades, particularly looking at softening phenomena and lifetimes obtained in isothermal cyclic loading (low cycle fatigue, LCF), non-isothermal cyclic loading (thermo-mechanical fatigue, TMF), and service-like combinations of fatigue and creep/relaxation periods. For this purpose, cylindrical specimens were extracted from thick-walled steam pipes, orthogonal to the pipe axis, and subjected to strain controlled cyclic loading (± 0.2 to ±0.5 % mechanical strain). Temperature intervals of TMF tests were chosen as either 300-620°C or 500-620°C, resembling so-called warm or hot start conditions of a power plant.
The test results will be presented and discussed with a focus on the impact of hold periods during testing (combined creep/relaxation-fatigue conditions) on mechanical softening, lifetime and formation of cracks. The findings will be complemented by results on the modification of the hierarchical ferritic-martensitic microstructure under different loading scenarios.
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.
This contribution contains the raw data used to compare experimental results with thermodynamic calculations using the CALPHAD method, which is related to the research article “The AlMo0.5NbTa0.5TiZr refractory high entropy superalloy: experimental findings and comparison with calculations using the CALPHAD method” [1] , and therefore this article can be used as a basis for interpreting the data contained therein. The AlMo0.5NbTa0.5TiZr refractory superalloy was characterized in the cast and annealed condition (1400 °C for 24 h) in order to measure grain size and to identify and measure the size and area fraction of the phases present. The raw data of this article include X-ray diffraction (XRD) measurements, microstructural characterization by scanning and transmission electron microscopy (SEM and TEM), and elemental analysis by energy dispersive X-ray spectroscopy (EDX). XRD includes the determination of phases and the lattice parameters (A2, B2, and hexagonal structure). Microstructural analysis by scanning and transmission electron microscopy includes (1) identification of composition, size, and volume fraction of the present phases and (2) determination of grain size. Based on these experimental data, it is possible to identify similarities and discrepancies with the data calculated using the CALPHAD method for the alloy under study in Ref. [1] , which provides the basis for better and more efficient development of reliable databases.
This contribution contains the raw data used to compare experimental results with thermodynamic calculations using the CALPHAD method, which is related to the research article “The AlMo0.5NbTa0.5TiZr refractory high entropy superalloy: experimental findings and comparison with calculations using the CALPHAD method” [1] , and therefore this article can be used as a basis for interpreting the data contained therein. The AlMo0.5NbTa0.5TiZr refractory superalloy was characterized in the cast and annealed condition (1400 °C for 24 h) in order to measure grain size and to identify and measure the size and area fraction of the phases present. The raw data of this article include X-ray diffraction (XRD) measurements, microstructural characterization by scanning and transmission electron microscopy (SEM and TEM), and elemental analysis by energy dispersive X-ray spectroscopy (EDX). XRD includes the determination of phases and the lattice parameters (A2, B2, and hexagonal structure). Microstructural analysis by scanning and transmission electron microscopy includes (1) identification of composition, size, and volume fraction of the present phases and (2) determination of grain size. Based on these experimental data, it is possible to identify similarities and discrepancies with the data calculated using the CALPHAD method for the alloy under study in Ref. [1] , which provides the basis for better and more efficient development of reliable databases.
The dataset provided in this repository comprises data obtained from a series of full-notch creep tests (FNCT) performed on selected high-density polyethylene (PE-HD) materials (for further details, see section 1 Materials in this document) in accordance with the corresponding standard ISO 16770.
The FNCT is one of the mechanical testing procedures used to characterize polymer materials with respect to their environmental stress cracking (ESC) behavior. It is widely applied for PE-HD materials, that are predominantly used for pipe and container applications. It is based on the determination of the time to failure for a test specimen under constant mechanical load in a well-defined and temperature controlled liquid environment. The test device used here also allows for continuous monitoring of applied force, specimen elongation and temperature.
Phase diagrams are the roadmaps for designing bulk phases. Similar to bulk, grain boundaries can possess various phases, but their phase diagrams remain largely unknown. Using a recently introduced density-based model, here we devise a strategy for computing multi-component grain boundary phase diagrams based on available bulk (CALPHAD) thermodynamic data. Fe-Mn-Cr, Fe-Mn-Ni, Fe-Mn-Co, Fe-Cr-Ni and Fe-Cr-Co alloy systems, as important ternary bases for several trending steels and high-entropy alloys, are studied. We found that despite its solute segregation enrichment, a grain boundary can have lower solubility limit than its corresponding bulk, promoting an interfacial chemical decomposition upon solute segregation. This is revealed here for the Fe-Mn-base alloy systems. The origins of this counter-intuitive feature are traced back to two effects, i.e., the magnetic ordering effect and the low cohesive energy of Mn solute element. Different aspects of interfacial phase stability and GB co-segregation in ternary alloys are investigated as well. We show that the concentration gradient energy contributions reduce segregation level but increase grain boundary solubility limit, stabilizing the GB against a chemical decomposition. Density-based grain boundary phase diagrams offer guidelines for systematic investigation of interfacial phase changes with applications to microstructure defects engineering.
Grain boundary (GB) chemical and structural variations can significantly influence materials performance. The former is generally ascribed to the structural gradient between the grain and GB. While GB segregation may be accompanied by chemical and structural variations, clear insights about the GB’s thermodynamic phase behaviour upon coupling between its chemistry and structure is lacking. Using the CALPHAD integrated density-based phase field model, we study the co-evolution of GB’s structure and segregation in Fe-Mn alloys. We found that the GB segregation transition is amplified if its structure can respond to chemical variation. Additionally, the coupling between GB structural and segregation evolution was found to enable co-existence of the spinodally formed low- and high-Mn phases within the GB. In the light of atomistic simulations, we expand on investigating the correlation between the parameters that characterise the GB density map with GB properties.
Systematic microstructure design requires reliable thermodynamic descriptions and phase diagrams of each and all microstructure elements. While such descriptions are well established for most bulk phases, thermodynamic assessment of crystal defects is greatly challenged by their individualistic aspects. In this talk, we present a density-based thermodynamic concept to describe defects based on available bulk thermodynamic data. Here dealing with grain boundaries (GBs), we apply this concept to compute GB (phase) diagram. Applications to segregation engineering of GBs in bulk and nanocrystalline alloys will be presented. We further develop this model to include the effect of elastic interactions due to atom size mismatch and obtain the corresponding GB (phase) diagram for the ternary Al-Cu-Li system.
Determination of Distribution Function used in MCS on Safety Analysis of Hydrogen Pressure Vessel
(2019)
The test data of static burst strength and load cycle strength of composite pressure vessels are often described by GAUSSian normal or WEIBULL distribution function to perform safety analyses. The goodness of assumed distribution function plays a significant role in the inferential statistics to predict the population properties by using limited test data. Often, GAUSSian and WEIBULL probability nets are empirical methods used to validate the distribution function; Anderson-Darling and KolmogorovSmirnov tests are the mostly favorable approaches for Goodness of Fit. However, the different approaches used to determine the parameters of distribution function lead mostly to different conclusions for safety assessments. In this study, six different methods are investigated to show the variations on the rates for accepting the composite pressure vessels according to GTR No. 13 life test procedure. The six methods are: a) NormLog based method, b) Least squares regression, c) Weighted least squares regression, d) A linear approach based on good linear unbiased estimators, e) Maximum likelihood estimation and f) The method of moments estimation. In addition, various approaches of ranking function are considered. In the study, Monte Carlo simulations are conducted to generate basic populations based on the distribution functions which are determined using different methods. Then the samples are extracted randomly from a population and evaluated to obtain acceptance rate. Here, the “populations” and “samples” are corresponding to the burst strength or load cycle strength of the pressure vessels made from composite material and a plastic liner (type 4) for the storage of hydrogen. To the end, the results are discussed, and the best reliable methods are proposed.
Determination of Distribution Function used in MCS on Safety Analysis of Hydrogen Pressure Vessel
(2019)
The test data of static burst strength and load cycle strength of composite pressure vessels are often described by GAUSSian normal or WEIBULL distribution function to perform safety analyses. The goodness of assumed distribution function plays a significant role in the inferential statistics to predict the population properties by using limited test data. Often, GAUSSian and WEIBULL probability nets are empirical methods used to validate the distribution function; Anderson-Darling and Kolmogorov-Smirnov tests are the mostly favorable approaches for Goodness of Fit. However, the different approaches used to determine the parameters of distribution function lead mostly to different conclusions for safety assessments. In this study, six different methods are investigated to show the variations on the rates for accepting the composite pressure vessels according to GTR No. 13 life test procedure. The six methods are: a) Norm-Log based method, b) Least squares regression, c) Weighted least squares regression, d) A linear approach based on good linear unbiased estimators, e) Maximum likelihood estimation and f) The method of moments estimation. In addition, various approaches of ranking function are considered. In the study, Monte Carlo simulations are conducted to generate basic populations based on the distribution functions which are determined using different methods. Then the samples are extracted randomly from a population and evaluated to obtain acceptance rate. Here, the “populations” and “samples” are corresponding to the burst strength or load cycle strength of the pressure vessels made from composite material and a plastic liner (type 4) for the storage of hydrogen. To the end, the results are discussed, and the best reliable methods are proposed.
The employment of renewable energy resources, which are naturally intermittent, for electricity generation has altered the working conditions of conventional power plants from continuous (baseload) to cyclic or flexible operation. For a long time, 9-12% Cr ferritic-martensitic stainless steels have been widely used in power plants due to their favorable characteristics such as high creep strength at high temperatures and oxidation and corrosion resistance. The components of power plants are subjected to long term cyclic loadings including fatigue and creep-fatigue at high temperatures. As ferritic-martensitic steels are known to exhibit cyclic softening when subjected to such loading scenarios, it is crucial to study the material response in such conditions. Since it is impossible to test the material behavior exactly as the operation conditions of power plants, due to technical difficulties and cost issues, it is necessary to develop physically-based material models that can predict the material behavior in more realistic situations.
In recent years, many material models have been proposed to describe the behavior of 9-12% Cr ferritic-martensitic stainless steels, which follow phenomenological or physically-based approaches. Phenomenological models provide a stress-strain relation based on empirical observations although they usually lack physical background. To alleviate this drawback and to allow for more flexibility and wider ranges of strain-rate and temperature, physically-based models are suggested. In this approach, microstructural evolution, dislocation movement and/or kinetics of plastic deformation processes are included. The physically-based models allow for a better extrapolation from the experiment results to other operation conditions and their material constants can be interpreted physically. Compared to the former approach, the material behavior can be described more accurately and flexibly and the number of material constants is less in general. In the presented work, the existing micromechanical models developed for P92 steel are compared and extended to allow for new dislocation-based strengthening/cyclic mechanisms. Their performance is assessed in the light of mechanical test data from creep-fatigue and thermo-mechanical fatigue experiments and detailed characterizations of the microstructure evolution in the fatigued material.
The formation of pores in CMSX-4 nickel based superalloys is detrimental to the service life of the material. A way to avoid the problem is to treat the superalloys under Hot Isostatic Pressing (HIP), which enables a large volume fraction of pores to be annihilated. This paper aims to understand the contribution of plastic activity related to the gliding of dislocations on the pore annihilation. Simulations based on a phase-field model of dislocation are performed and make it possible to consider the strong anisotropy of the CMSX-4 under HIP conditions in conjunction to the strong elastic heterogeneity introduced by the pore. For pores with a radius of few micrometers, it is shown that edge parts of dislocation lines that present an extra half atomic plane oriented towards the pore are stacked above and under it in the direction which is perpendicular to their slip-planes, causing an increase of the number of dislocation along the four octahedral directions of the FCC single crystal which intersect the pore center. Results are streamlined within the isotropic elastic theory of dislocations. Effects of elastic anisotropy and dislocation reactions are also investigated in order to specify what would be the dislocation configuration around a pore in CMSX-4 under HIP conditions. Notably, the elastic anisotropy is shown to significantly modify the arrangement of dislocations close to the pore equator. Simulations also allow for the characterization of pore/dislocation interactions when dislocations are involved in Low Angle Boundaries as experimentally observed.