Variational phase field model for dynamic brittle fracture

Simulating crack nucleation and propagation remains a challenging problematic because of the complexity of crack patterns observed in fracture mechanics experiments. Whereas some numerical methods aim at explicitly tracking the crack front evolution, an interesting alternative is offered by continuous approaches of brittle fracture which consist in representing the crack topology using a continuous field varying from 0 (sound material) to 1 (fully cracked material) across an internal length scale. This « phase field » approach benefits from a variational framework, strongly related to gradient damage models, and can be seen as a regularization of the variational approach to fracture developed by Francfort and Marigo in 1998. Moreover, it does not require any a priori knowledge of the crack path or topology, its evolution being driven only by energy minimization. Using such an approach combined to a finite-element discretization, the present work aims at providing some insights on crack propagation in a dynamic context. More specifically, crack branching (splitting of a single crack in two or more cracks) is a characteristic phenomenon of dynamic brittle fracture which still lacks a sound theoretical explanation. Numerical simulations will help us better understand some aspects of the branching phenomenon, especially the role played by material heterogeneities in the onset or delay of crack branching.

Constitutive behavior and fracture properties of tempered martensitic steels for nuclear applications

Raul Alejandro Bonadé

In this study, the tensile and fracture properties and the microstructure of the reduced-activation tempered martensitic steel Eurofer97 have been investigated. This technical alloy is a 9%Cr steel developed within the European fusion material research program. To a lesser extend, the plastic flow properties of a equiaxed ferritic Fe9%Cr model alloy were also studied for comparison with those of the tempered martensitic structure. The main objectives of this work are described hereafter: to correlate the microstructural features with the plastic flow properties measured by the tensile tests for both the Eurofer97 steel and the model alloy. The correlation established should be reflected in a physically-based model of plastic flow. to study the fracture properties of the Eurofer97 steel in details in the lower ductile-to-brittle transition region. to calculate with finite element modeling the stress fields at the crack tip. This information is further used in conjunction with a local approach model for quasi-cleavage to reconstruct the fracture toughness-temperature curve. Tensile tests were carried out at different imposed nominal strain rate at several temperatures from 77K up to 473K, on both the Eurofer97 steel and the Fe9%Cr model alloy. The temperature dependence of the yield stress was precisely determined. As expected for body-centered cubic (BCC) materials, a strong increase of the yield stress by decreasing temperature, below 200 K, was observed. At higher temperatures, the temperature dependence of the yield stress was found much weaker, being associated with the temperature dependence of the shear modulus. Efforts were made to analyze in details the post-yield behavior (strain-hardening) as a function of temperature. The post-yield behavior was modeled using the Kocks phenomenological model based on the competition of storage and annihilation of dislocations. While this model was originally developed for face-centered cubic (FCC) metals where the rate-controlling mechanism of dislocation motion is the dislocation-dislocation interaction, we used is model in the high temperature domain (T>200K) of BCC materials, to model the strain-hardening evolution of the Eurofer97 steel and the model alloy. The values obtained for the key parameters of the model, namely the dislocation mean free path and annihilation coefficient, were found consistent with the microstructural features. The parameters temperature-dependence observed was also consistent with the physics of two basic mechanisms of dislocation storage and annihilation (dynamical recovery). For the low temperature domain, the strain-hardening model was modified to account for the strong Peierls lattice friction, based on an original idea of Rauch. Transmission electron microscopy observations were done to characterize the undeformed microstructures and their evolution with strain. A clear correlation was established between the stress-dependence of the strain-hardening and the microstructures. Fracture toughness tests on the Eurofer97 steel were performed with 0.2T C(T) and 0.4T C(T) specimens in the lower transition. Five temperatures were selected at which many tests were repeated to determine the amplitude of the inherent scatter of quasi-cleavage. The temperature were chosen in such a way that the measure fracture toughness remain below 150 MPa m1/2, which correspond for the 0.4T C(T) at a M value larger than 70 at the highest temperature. Such a M value for 0.4T C(T) specimens is known to ensure enough constraint. An attempt to analyze the experimental data in the framework of the master-curve approach by following the ASTM E-1921 standard was done. It was clearly demonstrated that the assumed shape of the toughness-temperature curve as described in the ASTM E-1921 standard for the reactor pressure vessel steels is not adequate for the tempered martensitic Eurofer97 steel, which present a particularly steep transition. Our Eurofer97 fracture database was then compared to the existing one on another similar steel, the F82H. Differences in the amplitude of the scatter of both steels was found while the lower bound of the toughness-temperature curve, describing the 1% failure probability was shown to be the same. 2D finite element simulations of the compact tension specimens were performed at various temperatures using the constitutive laws determined previously. The stress field around the crack tip were calculated and used to determine a local criterion of quasi-cleavage. The criterion is defined by the attainment of a critical stress encompassing a critical area. The lower bound of the toughness-temperature curve was then successfully reconstructed by using this local criterion. Finally, the relationship between the critical area and the applied stress intensity factor for the C(T) specimen was shown to follow a power law whose coefficients are dependent on the real dimensions of the specimen. Such a relationship allows scaling the C(T) toughness data from one size to another in case of in-plane constraint loss.

EPFL2006

Finite element modeling and experimental study of brittle fracture in tempered martensitic steels for thermonuclear fusion applications

Pablo Federico Mueller

In this work we have studied brittle fracture in high-chromium reduced activation tempered martensitic steels foreseen as structural materials for thermonuclear fusion reactors. Developing the adequate materials that can withstand the severe irradiation conditions of the burning plasma in a future fusion reactor is one of the major challenges to be solved in order to make profit from the great advantages of thermonuclear fusion as an energy source. High-chromium tempered martensitic steels such as F82H and the most advanced version Eurofer97 are among the main candidate materials for structural applications in future fusion power plants due to low irradiation-induced swelling, good mechanical and thermal properties, and reasonably fast radioactive decay. The most concerning drawback of these kind of steels is irradiation embrittlement, which is manifested by a ductile-to-brittle transition temperature shift to higher temperatures after irradiation whose amplitude depends on the irradiation conditions (temperature, neutron flux, neutron fluence, etc). The aim of this work was to study and model brittle fracture in the ductile-brittle transition region of this kind of steels in the as-received unirradiated conditions. It is necessary to be able to transfer laboratory specimen fracture data to real components and structures in order to assess the performance of these steels in the different operating and transient conditions they could find during the operation life of a fusion reactor. In order to do so, the specimen geometry effects and specimen size effects on measured fracture toughness need to be properly understood, taken into account and predicted with an appropriate model. In particular, specimen size effect on measured toughness is a major concern for the nuclear materials research community owing to the limited irradiation volume in current and planed materials irradiation facilities. The main results of this PhD work are summarized below. The microstructure of Eurofer97 and F82H has been characterized and compared by means of optical microscopy, scanning electron microscopy, transmission electron microscopy and energy-dispersive X-ray spectroscopy in order to identify microstructural features that could play a role in the measured fracture toughness. Both steels have similar but slightly different chemical composition and final heat-treatments but the prior austenitic grain size measured in F82H is approximately 8 times larger than in Eurofer97. It was shown that the alloying element Tantalum, added to stabilize the austenite grain size, played a different role in both steels. After a careful analysis of the particles present in both steels, it was found that Tantalum in Eurofer97 formed carbides of an average size around 100 nanometers. In contrast in F82H it did not form small carbides but formed big oxide inclusions with a size up to 30 µm. These large particles do not effectively pin the grain boundaries. The different behavior of Tantalum in these steels is believed to be mainly a consequence of the larger content of Oxygen present and the smaller amount of Aluminum in F82H compared to Eurofer97. The Master-Curve ASTM-E1921 standard is a method initially developed to determine the ductile-to-brittle transition reference temperature, T0, in fission reactor ferritic steels from a small number of experiments. In this work the applicability of the Master-Curve method to reduced activation tempered martensitic steels such as Eurofer97 and F82H was studied in detail. Fracture tests with pre-cracked sub-sized compact tension specimens (three different sizes, 0.18T, 0.35T and 0.87T of Eurofer97) were carried out in the temperature range [-196 °C, -40 °C]. The toughness-temperature behavior and scatter were shown to deviate from the ASTM-E1921 standard predictions near the lower shelf. Using the method of maximum likelihood, the athermal component of the Master-Curve was calculated to better fit our fracture toughness data from the lower to the middle transition region. We showed that these Master-Curve adjustments are necessary to make the T0 values obtained near the lower shelf with 0.35T size C(T) specimens consistent with those obtained in the middle transition region with 0.87T C(T) specimens. The ASTM-E1921 specimen size limitations, setting the maximum toughness measurable with a given specimen size, were found to be too lenient for this kind of steels. This problem was especially evident in fracture toughness data of Eurofer97 and F82H obtained in the upper transition temperature range with two different specimen sizes. Thus a more stringent specimen size requirement was proposed to avoid inconsistent transition temperature determinations. A promising local fracture model with the potential of predicting cleavage fracture toughness was studied. Finite element simulations were undertaken for compact specimens, notched specimens and tensile specimens of Eurofer97 steel tested from 20 °C down to -197 °C. Three and two dimensional as well as axisymmetric simulations were run in order to calculate the stress and strain fields at the onset of brittle fracture. For each tested temperature, the calculated load-displacement curves were found to reproduce very well those of the experiments. A local approach fracture criterion was studied. This criterion states that when the maximum principal stress is larger than a critical stress within a critical volume ahead of the crack tip, notch or neck, cleavage is triggered leading to macroscopic fracture of the specimen. It was shown that this model is able to predict the minimum fracture load of the notched specimens with the same values of critical stress and critical volume that were calibrated to predict the lower bound of fracture in compact specimens. This local approach model was also successfully used to predict the strong size effect observed experimentally in pre-cracked compact tension specimens in the upper transition region. The critical fracture stress determined in the standard tensile specimens was found higher than that of the fracture specimens. It was suggested that this difference stems from a significant difference in the stress state between the different specimens (triaxiality). Finally, the comparison between the overall fracture behavior in the transition between F82H and Eurofer97 steels indicated that these two materials are quite similar. A difference in the reference temperature T0 of about 20 °C was found, with a nominal value of -100 °C and -80 °C for F82H and Eurofer97 respectively. However, some excessive scatter in the toughness data was found in F82H with more data points than expected lying below the 1% tolerance bound; this was not observed in Eurofer97. From an engineering point of view, the entire fracture databases of these two steels are encompassed practically by the same lower tolerance bound.

EPFL2009

A meso-scale computational approach to dynamic failure in concrete

Leonardo Snozzi

The aim of the present thesis is to provide a numerical tool for the mechanical modeling of dynamic failure in concrete. This widely used construction material is characterized by a non-linear failure behavior, which is consequently difficult to describe with macroscopic quantities. In order to grasp the mechanical damage response of this heterogeneous brittle material, our approach is based on a meso-mechanical point of view. The whole research work is, therefore, mainly focused on the understanding of the mechanisms that are linked with the material’s heterogeneous composition. The advantage of this level of observation is that it allows to represent the most important concrete constituents (e.g. aggregates and cement paste), thus facilitating the physical identification of the material parameters of the model and of the mechanisms (interaction between matrix and inclusions) that characterize its constitutive behavior. For this purpose we exploit the capabilities of a two-dimensional finite-element model. The onset of fracture is explicitly modeled using the well-known cohesive approach. We first investigate the dynamic tensile response of concrete specimens. Different simulations are carried out to assess the influence of aggregates properties on peak strength and dissipated fracture energy at different strain rates. The model aims to explain the strain-rate strengthening through structural effects. However, to capture the full extent observed experimentally, our results suggest that is not possible to discard the combination of inertial with material rate hardening mechanisms. Next, in order to account for crack-interactions as well as path dependent behavior, the model is enriched with the introduction of an explicit contact algorithm and a mode-dependent fracture energy. To demonstrate the capability of the proposed approach to provide accurate results, the model is first applied to two benchmark tests in masonry engineering. Afterwards, the developed framework is applied to reproduce dynamic compressive failure of meso-scale concrete samples. Simulations involving different strain rates and levels of lateral confinement are conducted. An energetic analysis shows that dissipation of energy through fracture and friction is an increasing function of the applied confinement and strain rate. Our results underline thus the importance of capturing frictional mechanisms, which appear to dissipate a raising amount of frictional energy with increasing strain and applied pressure. Finally, a multi-scale computational framework is developed to up scale the obtained fine scale response to the coarse scale. The selected approach considers concrete macroscopically as being homogeneous and to behave linear elastically aside from the propagating cracks. The cohesive macroscopic tractions are however extracted from the response of meso-scale representative volume elements, which are activated on the fly when a macroscopic discontinuity propagates. Parallel simulation shows the capability of the model to predict the structural response of a tested unit including the physical mechanisms occurring at the fine scale.

EPFL2013