Search

The phase field fracture method is attracting significant interest. Phase field approaches have enabled predicting - on arbitrary geometries and dimensions - complex fracture phenomena such as crack branching, coalescence, deflection and nucleation. In this work, we present a simple and robust implementation of the phase field fracture method in the commercial finite element package Abaqus. The implementation exploits the analogy between the phase field evolution law and the heat transfer equation, enabling the use of Abaqus’ in-built features and circumventing the need for defining user elements. The framework is general, and is shown to accommodate different solution schemes (staggered and monolithic), as well as various constitutive choices for preventing damage under compression. The robustness and applicability of the numerical framework presented is demonstrated by addressing several 2D and 3D boundary value problems of particular interest. Focus is on the solution of paradigmatic case studies that are known to be particularly demanding from a convergence perspective. The results reveal that our phase field fracture implementation can be readily combined with other advanced computational features, such as contact, and deliver robust and precise solutions. The code developed can be downloaded from www.empaneda.com/codes .

We present a simple and robust implementation of the phase field fracture method in Abaqus. Unlike previous works, only a user material (UMAT) subroutine is used. This is achieved by exploiting the analogy between the phase field balance equation and heat transfer, which avoids the need for a user element mesh and enables taking advantage of Abaqus' in-built features. A unified theoretical framework and its implementation are presented, suitable for any arbitrary choice of crack density function and fracture driving force. Specifically, the framework is exemplified with the so-called AT1, AT2 and phase field-cohesive zone models (PF-CZM). Both staggered and monolithic solution schemes are handled. We demonstrate the potential and robustness of this new implementation by addressing several paradigmatic 2D and 3D boundary value problems. The numerical examples show how the current implementation can be used to reproduce numerical and experimental results from the literature, and efficiently capture advanced features such as complex crack trajectories, crack nucleation from arbitrary sites and contact problems. The code developed can be downloaded from www.empaneda.com/codes.

Experiments and finite element (FE) calculations were performed to study the raster angle–dependent fracture behaviour of acrylonitrile butadiene styrene (ABS) thermoplastic processed via fused filament fabrication (FFF) additive manufacturing (AM). The fracture properties of 3D-printed ABS were characterized based on the concept of essential work of fracture (EWF), utilizing double-edge-notched tension (DENT) specimens considering rectilinear infill patterns with different raster angles (0°, 90° and + 45/− 45°). The measurements showed that the resistance to fracture initiation of 3D-printed ABS specimens is substantially higher for the printing direction perpendicular to the crack plane (0° raster angle) as compared to that of the samples wherein the printing direction is parallel to the crack (90° raster angle), reporting EWF values of 7.24 kJ m−2 and 3.61 kJ m−2, respectively. A relatively high EWF value was also reported for the specimens with + 45/− 45° raster angle (7.40 kJ m−2). Strain field analysis performed via digital image correlation showed that connected plastic zones existed in the ligaments of the DENT specimens prior to the onset of fracture, and this was corroborated by SEM fractography which showed that fracture proceeded by a ductile mechanism involving void growth and coalescence followed by drawing and ductile tearing of fibrils. It was further shown that the raster angle–dependent strength and fracture properties of 3D-printed ABS can be predicted with an acceptable accuracy by a relatively simple FE model considering the anisotropic elasticity and failure properties of FFF specimens. The findings of this study offer guidelines for fracture-resistant design of AM-enabled thermoplastics. Graphical abstract

The phase field paradigm, in combination with a suitable variational structure, has opened a path for using Griffith's energy balance to predict the fracture of solids. These so-called phase field fracture methods have gained significant popularity over the past decade, and are now part of commercial finite element packages and engineering fitness- for-service assessments. Crack paths can be predicted, in arbitrary geometries and dimensions, based on a global energy minimization-without the need for ad hoc criteria. In this work, we review the fundamentals of phase field fracture methods and examine their capabilities in delivering predictions in agreement with the classical fracture mechanics theory pioneered by Griffith. The two most widely used phase field fracture models are implemented in the context of the finite element method, and several paradigmatic boundary value problems are addressed to gain insight into their predictive abilities across all cracking stages; both the initiation of growth and stable crack propagation are investigated. In addition, we examine the effectiveness of phase field models with an internal material length scale in capturing size effects and the transition flaw size concept. Our results show that phase field fracture methods satisfactorily approximate classical fracture mechanics predictions and can also reconcile stress and toughness criteria for fracture. The accuracy of the approximation is however dependent on modelling and constitutive choices; we provide a rationale for these differences and identify suitable approaches for delivering phase field fracture predictions that are in good agreement with well-established fracture mechanics paradigms. This article is part of a discussion meeting issue 'A cracking approach to inventing new tough materials: fracture stranger than friction'., Royal Commission for the 1851 Exhibition

We present a computational framework to explore the effect of microstructure and constituent properties upon the fracture toughness of fibre-reinforced polymer composites. To capture microscopic matrix cracking and fibre-matrix debonding, the framework couples a phase field fracture method and a cohesive zone model in the context of the finite element method. Virtual single-notched three point bending tests on fibre reinforced composites are conducted. The actual microstructure of the composite is simulated by an embedded cell in the fracture process zone, while the remaining area is homogenised to be an anisotropic elastic solid. A detailed comparison of the predicted results with experimental observations reveals that it is possible to accurately capture the crack path, interface debonding and load versus displacement response. The sensitivity of the crack growth resistance curve (R-curve) to the matrix fracture toughness and the fibre-matrix interface properties is determined. The influence of porosity upon the R-curve of fibre-reinforced composites is also explored, revealing a higher crack growth resistance with increasing void volume fraction. These results shed light into microscopic fracture mechanisms and set the basis for efficient design of high fracture toughness composites.

The phase field fracture method has emerged as a promising computational tool for modelling a variety of problems including, since recently, hydrogen embrittlement and stress corrosion cracking. In this work, we demonstrate the potential of phase field-based multi-physics models in transforming the engineering assessment and design of structural components in hydrogen-containing environments. First, we present a theoretical and numerical framework coupling deformation, diffusion and fracture, which accounts for inertia effects. Several constitutive choices are considered for the crack density function, including choices with and without an elastic phase in the damage response. The material toughness is defined as a function of the hydrogen content using an atomistically-informed hydrogen degradation law. The model is numerically implemented in 2D and 3D using the finite element method. The resulting computational framework is used to address a number of case studies of particular engineering interest. These are intended to showcase the model capabilities in: (i) capturing complex fracture phenomena, such as dynamic crack branching or void-crack interactions, (ii) simulating standardised tests for critical components, such as bolts, and (iii) enabling simulation-based paradigms such as Virtual Testing or Digital Twins by coupling model predictions with inspection data of large-scale engineering components. The evolution of defects under in-service conditions can be predicted, up to the ultimate failure. By reproducing the precise geometry of the defects, as opposed to re-characterising them as sharp cracks, phase field modelling enables more realistic and effective structural integrity assessments.

We present a generalised framework for resolving the electrochemistry-diffusion interface and modelling hydrogen transport near a crack tip. The adsorption and absorption kinetics are captured by means of Neumann-type generalised boundary conditions. The diffusion model includes the role of trapping, with a constant or evolving trap density, and the influence of the hydrostatic stress. Both conventional plasticity and strain gradient plasticity are used to model the mechanical behaviour of the solid. Notable differences are found in the estimated crack tip hydrogen concentrations when comparing with the common procedure of prescribing a constant hydrogen concentration at the crack surfaces., E. Martínez-Pañeda acknowledges financial support from EPSRC funding under grant no.EP/R010161/1, from the UKCRIC Coordination Node EPSRC grant number EP/R017727/1, which funds UKCRIC's ongoing coordination, and from Wolfson College Cambridge (Junior Research Fellowship). A. Díaz gratefully acknowledges financial support from the Ministry of Science, Innovation and Universities of Spain through grant RTI2018-096070-B-C33.

Failures attributed to hydrogen embrittlement are a major concern for metals so a better understanding of damage micro-mechanisms and hydrogen diffusion within the metal is needed. Local concentrations depend on transport phenomena including trapping effects, which are usually characterised by a temperature-programmed desorption method often referred to as Thermal Desorption Analysis (TDA). When the hydrogen is released from the specimen during the programmed heating, some desorption peaks are observed that are commonly related to detrapping energies by means of an analytical procedure. The limitations of this approach are revisited here and gaseous hydrogen charging at high temperatures is simulated. This popular procedure enables attaining high concentrations due to the higher solubility of hydrogen at high temperatures. However, the segregation behaviour of hydrogen into traps depends on charging time and temperature. This process and the subsequent cooling alter hydrogen distribution are numerically modelled; it is found that TDA spectra are strongly affected by the charging temperature and the charging time, both for weak and strong traps. However, the influence of ageing time at room temperature after cooling and before desorption is only appreciable for weak traps., Ministry of Economy and Competitiveness of Spain through grants MAT2014-58738 and RTI2018-096070-B-C33.

We present a gradient-based theoretical framework for predicting hydrogen assisted fracture in elastic-plastic solids. The novelty of the model lies in the combination of: (i) stress-assisted diffusion of solute species, (ii) strain gradient plasticity, and (iii) a hydrogen-sensitive phase field fracture formulation, inspired by first principles calculations. The theoretical model is numerically implemented using a mixed finite element formulation and several boundary value problems are addressed to gain physical insight and showcase model predictions. The results reveal the critical role of plastic strain gradients in rationalising decohesion-based arguments and capturing the transition to brittle fracture observed in hydrogen-rich environments. Large crack tip stresses are predicted, which in turn raise the hydrogen concentration and reduce the fracture energy. The computation of the steady state fracture toughness as a function of the cohesive strength shows that cleavage fracture can be predicted in otherwise ductile metals using sensible values for the material parameters and the hydrogen concentration. In addition, we compute crack growth resistance curves in a wide variety of scenarios and demonstrate that the model can appropriately capture the sensitivity to: the plastic length scales, the fracture length scale, the loading rate and the hydrogen concentration. Model predictions are also compared with fracture experiments on a modern ultra-high strength steel, AerMet100. A promising agreement is observed with experimental measurements of threshold stress intensity factor Kth over a wide range of applied potentials.

Interface cracking is one of the most prominent failure modes in fibre reinforced polymer (FRP) composites. Recent trends in high-tech applications of FRP composites exploit the limits of the load bearing capacity, generally encompassing the development of notable nonlinear effects from geometrical and material signatures. In this investigation, we present a comprehensive assessment of the new Linear Elastic Brittle Interface Model (LEBIM) in geometrically nonlinear applications undergoing mixed-mode fracture conditions. This interface model for triggering fracture events is formulated through the advocation of continuum-like assumptions (for initial non-zero interface thickness) and allows the incorporation of the potential role of the in-plane deformation effects. The performance of the present interface model is demonstrated through the simulation of specimens with mixed-mode delamination, with special attention for its application in samples equipped with structured interfaces. Current predictions exhibit an excellent agreement with respect to experimental data, validating the proposed methodology., Composite Structures (2020)

In this work, we extend the recently proposed adaptive phase field method to model fracture in orthotropic functionally graded materials (FGMs). A recovery type error indicator combined with quadtree decomposition is employed for adaptive mesh refinement. The proposed approach is capable of capturing the fracture process with a localized mesh refinement that provides notable gains in computational efficiency. The implementation is validated against experimental data and other numerical experiments on orthotropic materials with different material orientations. The results reveal an increase in the stiffness and the maximum force with increasing material orientation angle. The study is then extended to the analysis of orthotropic FGMs. It is observed that, if the gradation in fracture properties is neglected, the material gradient plays a secondary role, with the fracture behaviour being dominated by the orthotropy of the material. However, when the toughness increases along the crack propagation path, a substantial gain in fracture resistance is observed.

We investigate the influence of microstructural traps on hydrogen diffusion and embrittlement in the presence of cyclic loads. A mechanistic, multi-trap model for hydrogen transport is developed, implemented into a finite element framework, and used to capture the variation of crack tip lattice and trapped hydrogen concentrations as a function of the loading frequency, the trap binding energies and the trap densities. We show that the maximum value attained by the lattice hydrogen concentration during the cyclic analysis exhibits a notable sensitivity to the ratio between the loading frequency and the effective diffusion coefficient. This is observed for both hydrogen pre-charged samples (closed-systems) and samples exposed to a permanent source of hydrogen (open-systems). Experiments are used to determine the critical concentration for embrittlement, by mapping the range of frequencies where the output is the same as testing in inert environments. We then quantitatively investigate and discuss the implications of developing materials with higher trap densities in mitigating embrittlement in the presence of cyclic loads. It is shown that, unlike the static case, increasing the density of “beneficial traps” is a viable strategy in designing alloys resistant to hydrogen assisted fatigue for both closed- and open-systems.

The ductile fracture process in porous metals due to growth and coalescence of micron scale voids is affected not only by the imposed stress state but also by the distribution of the voids and the material size effect. The objective of this study is to understand the interaction of the inter-void spacing (or ligaments) and the resultant gradient-induced material size effect on void coalescence for a range of imposed stress states. To this end, three-dimensional finite element calculations of unit cell models with a discrete void embedded in a strain gradient-enhanced material matrix are performed. The calculations are carried out for a range of initial inter-void ligament sizes and imposed stress states characterized by fixed values of the stress triaxiality and the Lode parameter. Our results show that in the absence of strain gradient effects on the material response, decreasing the inter-void ligament size results in an increase in the propensity for void coalescence. However, in a strain gradient-enhanced material matrix, the strain gradients harden the material in the inter-void ligament and decrease the effect of inter-void ligament size on the propensity for void coalescence.

In this work, distortion gradient plasticity is used to gain insight into material deformation ahead of a crack tip. This also constitutes the first fracture mechanics analysis of gradient plasticity theories adopting Nye's tensor as primal kinematic variable. First, the asymptotic nature of crack tip fields is analytically investigated. We show that an inner elastic region exists, adjacent to the crack tip, where elastic strains dominate plastic strains and Cauchy stresses follow the linear elastic stress singularity. This finding is verified by detailed finite element analyses using a new numerical framework, which builds upon a viscoplastic constitutive law that enables capturing both rate-dependent and rate-independent behaviour in a computationally efficient manner. Numerical analysis is used to gain further insight into the stress elevation predicted by distortion gradient plasticity, relative to conventional J2 plasticity, and the influence of the plastic spin under both mode I and mixed-mode fracture conditions. It is found that Nye's tensor contributions have a weaker effect in elevating the stresses in the plastic region, while predicting the same asymptotic behaviour as constitutive choices based on the plastic strain gradient tensor. A minor sensitivity to X, the parameter governing the dissipation due to the plastic spin, is observed. Finally, distortion gradient plasticity and suitable higher order boundary conditions are used to appropriately model the phenomenon of brittle failure along elastic-plastic material interfaces. We reproduce paradigmatic experiments on niobium-sapphire interfaces and show that the combination of strain gradient hardening and dislocation blockage leads to interface crack tip stresses that are larger than the theoretical lattice strength, rationalising cleavage in the presence of plasticity at bi-material interfaces., Comment: International Journal of Engineering Science (2020)

The onset of sub-critical crack growth during slow strain rate tensile testing (SSRT) is assessed through a combined experimental and modeling approach. A systematic comparison of the extent of intergranular fracture and expected hydrogen ingress suggests that hydrogen diffusion alone is insufficient to explain the intergranular fracture depths observed after SSRT experiments in a Ni–Cu superalloy. Simulations of these experiments using a new phase field formulation indicate that crack initiation occurs as low as 40% of the time to failure. The implications of such sub-critical crack growth on the validity and interpretation of SSRT metrics are then explored.

The electrochemical permeation test is one of the most used methods for characterising hydrogen diffusion in metals. The flux of hydrogen atoms registered in the oxidation cell might be fitted to obtain apparent diffusivities. The magnitude of this coefficient has a decisive influence on the kinetics of fracture or fatigue phenomena assisted by hydrogen and depends largely on hydrogen retention in microstructural traps. In order to improve the numerical fitting of diffusion coefficients, a permeation test has been reproduced using FEM simulations considering two approaches: a continuum 1D model in which the trap density, binding energy and the input lattice concentrations are critical variables and a polycrystalline model where trapping at grain boundaries is simulated explicitly including a segregation factor and a diffusion coefficient different from that of the interior of the grain. Results show that the continuum model captures trapping delay, but it should be modified to model the trapping influence on the steady state flux. Permeation behaviour might be classified according to different regimes depending on deviation from Fickian diffusion. Polycrystalline synthetic permeation shows a strong influence of segregation on output flux magnitude. This approach is able to simulate also the short-circuit diffusion phenomenon. The comparison between different grain sizes and grain boundary thicknesses by means of the fitted apparent diffusivity shows the relationships between the registered flux and the characteristic parameters of traps., The authors gratefully acknowledge financial support from the project MINECO Refs: MAT2014-58738-C3-2-R and RTI2018-096070-B-C33. E. Martínez-Pañeda acknowledges financial support from the People Programme (Marie Curie Actions) of the European Union’s Seventh Framework Programme (FP7/2007-2013) under REA Grant Agreement No. 609405 (COFUNDPostdocDTU).

The shear strength of a pre-cracked sandwich layer is predicted, assuming that the layer is linear elastic or elastic-plastic, with yielding characterized either by the J2 plasticity theory or by a strip-yield model. The substrates are elastic and of dissimilar modulus to that of the layer. Two geometries are analyzed: (i) a semi-infinite crack in a sandwich layer, subjected to a remote mode II K-field and (ii) a center-cracked sandwich plate of finite width under remote shear stress. For the semi-infinite crack, the near-tip stress field is determined as a function of elastic mismatch, and crack tip plasticity is either prevented (the elastic case) or duly accounted for (the elastic-plastic case). Analytical and numerical solutions are then obtained for the center-cracked sandwich plate of the finite width. First, a mode II K-calibration is obtained for a finite crack in the elastic sandwich layer. Second, the analysis is extended to account for crack tip plasticity via a mode II strip-yield model of finite strength and finite toughness. The analytical predictions are verified by finite element simulations, and a failure map is constructed in terms of specimen geometry and crack length.

Additive manufacturing is becoming increasingly popular in academia and industry. Accordingly, there has been a growing interest in characterizing 3D printed samples to determine their structural integrity behaviour. We employ the Essential Work of Fracture (EWF) to investigate the mechanical response of polymer sheets obtained through additive manufacturing. Our goal is twofold; first, we aim at gaining insight into the role of fibre reinforcement on the fracture resistance of additively manufactured polymer sheets. Deeply double-edge notched tensile (DDEN-T) tests are conducted on four different polymers: Onyx, a crystalline, nylon-reinforced polymer, and three standard polymers used in additive manufacturing – PLA, PP and ABS. Results show that fibre-reinforcement translates into a notable increase in fracture resistance, with the fracture energy of Onyx being an order of magnitude higher than that reported for non-reinforced polymers. On the other hand, we propose the use of a miniature test specimen, the deeply double-edge notched small punch specimens (DDEN-SP), to characterize the mechanical response using a limited amount of material. The results obtained exhibit good alignment with the DDEN-T data, suggesting the suitability of the DDEN-SP test for measuring fracture properties of additively manufactured polymers in a cost-effective manner., Ministry of Education, Government of Castile-Leon, which started in 2018 under the support of the Recognized Research Groups of the Public Universities of Castile and Leon (project: BU033G18)

We present a phase field formulation for fracture in functionally graded materials (FGMs). The model builds upon homogenization theory and accounts for the spatial variation of elastic and fracture properties. Several paradigmatic case studies are addressed to demonstrate the potential of the proposed modelling framework. Specifically, we (i) gain insight into the crack growth resistance of FGMs by conducting numerical experiments over a wide range of material gradation profiles and orientations, (ii) accurately reproduce the crack trajectories observed in graded photodegradable copolymers and glass-filled epoxy FGMs, (iii) benchmark our predictions with results from alternative numerical methodologies, and (iv) model complex crack paths and failure in three dimensional functionally graded solids. The suitability of phase field fracture methods in capturing the crack deflections intrinsic to crack tip mode-mixity due to material gradients is demonstrated. Material gradient profiles that prevent unstable fracture and enhance crack growth resistance are identified: this provides the foundation for the design of fracture resistant FGMs. The finite element code developed can be downloaded from www.empaneda.com/codes., E. Martínez-Pañeda acknowledges financial support from the Royal Commission for the 1851 Exhibition through their Research Fellowship programme (RF496/2018).

Thick adhesive layers have potential structural application in ship construction for the joining of a composite superstructure to a steel hull. The purpose of this study is to develop a mechanics model for the adhesive fracture of such lap joints under shear loading. Modified Thick-Adherend-Shear-Test (TAST) specimens made from a MMA-based adhesive and steel adherents are designed and fabricated. Crack initiation and growth of these joints is measured and monitored by Digital Image Correlation (DIC). An attempt is made to use a cohesive zone model to predict the magnitude of shear strain across the adhesive layer both at crack initiation and at peak load, and to predict the extent of crack growth as a function of shear strain across the adhesive layer. The ability of a cohesive zone model to predict several features of specimen failure is assessed for the case of an adhesive layer of high shear ductility.

© 2017, Springer-Verlag Berlin Heidelberg. Strain gradient plasticity theories are being widely used for fracture assessment, as they provide a richer description of crack tip fields by incorporating the influence of geometrically necessary dislocations. Characterizing the behavior at the small scales involved in crack tip deformation requires, however, the use of a very refined mesh within microns to the crack. In this work a novel and efficient gradient-enhanced numerical framework is developed by means of the extended finite element method (X-FEM). A mechanism-based gradient plasticity model is employed and the approximation of the displacement field is enriched with the stress singularity of the gradient-dominated solution. Results reveal that the proposed numerical methodology largely outperforms the standard finite element approach. The present work could have important implications on the use of microstructurally-motivated models in large scale applications. The non-linear X-FEM code developed in MATLAB can be downloaded from www.empaneda.com/codes.

The mode I crack tip asymptotic response of a solid characterised by strain gradient plasticity is investigated. It is found that elastic strains dominate plastic strains near the crack tip, and thus the Cauchy stress and the strain state are given asymptotically by the elastic K -field. This crack tip elastic zone is embedded within an annular elasto-plastic zone. This feature is predicted by both a crack tip asymptotic analysis and a finite element computation . When small scale yielding applies, three distinct regimes exist: an outer elastic K field, an intermediate elasto-plastic field, and an inner elastic K field. The inner elastic core significantly influences the crack opening profile. Crack tip plasticity is suppressed when the material length scale l of the gradient theory is on the order of the plastic zone size estimation, as dictated by the remote stress intensity factor . A generalized J -integral for strain gradient plasticity is stated and used to characterise the asymptotic response ahead of a short crack. Finite element analysis of a cracked three point bend specimen reveals that the crack tip elastic zone persists in the presence of bulk plasticity and an outer J-field.

The fracture toughness for a mode I/II crack propagating in a ductile material has been subject to numerous investigations. However, the influence of the material hardening law has received very limited attention, with isotropic hardening being the default choice if cyclic loads are absent. The present work extends the existing studies of monotonic mode I/II steady-state crack propagation with the goal to compare the predictions from an isotropic hardening model with that of a kinematic hardening model. The work is conducted through a purpose-built steady-state framework that directly delivers the steady-state solution. In order to provide a fracture criterion, a cohesive zone model is adopted and embedded at the crack tip in the steady-state framework, while a control algorithm for the far-field, that significantly reduces the number of equilibrium iterations is employed to couple the far-field loading to the correct crack tip opening. Results show that the steady-state fracture toughness (shielding ratio) obtained for a kinematic hardening material is larger than for the corresponding isotropic hardening case. The difference between the isotropic and kinematic model is tied to the non-proportional loading conditions and reverse plasticity. This also explains the vanishing difference in the shielding ratio when considering mode II crack propagation as the non-proportional loading is less pronounced and the reverse plasticity is absent.

In this work a novel methodology for fracture toughness characterization by means of the small punch test (SPT) is presented. Notched specimens are employed and fracture resistance is assessed through a critical value of the notch mouth displacement δ SPT . Finite element simulations and interrupted experiments are used to track the evolution of δ SPT as a function of the punch displacement. The onset of crack propagation is identified by means of a ductile damage model and the outcome is compared to the crack tip opening displacement estimated from conventional tests at crack initiation. The proposed numerical–experimental scheme is examined with two different grades of CrMoV steel and the differences in material toughness captured. Limitations and uncertainties arising from the different damage phenomena observed in the lowest toughness material examined are thoroughly discussed.

We present a probabilistic framework for brittle fracture that builds upon Weibull statistics and strain gradient plasticity . The constitutive response is given by the mechanism-based strain gradient plasticity theory, aiming to accurately characterize crack tip stresses by accounting for the role of plastic strain gradients in elevating local strengthening ahead of cracks. It is shown that gradients of plastic strain elevate the Weibull stress and the probability of failure for a given choice of the threshold stress and the Weibull parameters. The statistical framework presented is used to estimate failure probabilities across temperatures in ferritic steels . The framework has the capability to estimate the three statistical parameters present in the Weibull-type model without any prior assumptions. The calibration against experimental data shows important differences in the values obtained for strain gradient plasticity and conventional J2 plasticity. Moreover, local probability maps show that potential damage initiation sites are much closer to the crack tip in the case of gradient-enhanced plasticity. Finally, the fracture response across the ductile-to-brittle regime is investigated by computing the cleavage resistance curves with increasing temperature. Gradient plasticity predictions appear to show a better agreement with the experiments.

We investigate the potential of quasi-Newton methods in facilitating convergence of monolithic solution schemes for phase field fracture modelling. Several paradigmatic boundary value problems are addressed, spanning the fields of quasi-static fracture, fatigue damage and dynamic cracking. The finite element results obtained reveal the robustness of quasi-Newton monolithic schemes, with convergence readily attained under both stable and unstable cracking conditions. Moreover, since the solution method is unconditionally stable, very significant computational gains are observed relative to the widely used staggered solution schemes. In addition, a new adaptive time increment scheme is presented to further reduce the computational cost while allowing to accurately resolve sudden changes in material behavior, such as unstable crack growth. Computation times can be reduced by several orders of magnitude, with the number of load increments required by the corresponding staggered solution being up to 3000 times higher. Quasi-Newton monolithic solution schemes can be a key enabler for large scale phase field fracture simulations. Implications are particularly relevant for the emerging field of phase field fatigue, as results show that staggered cycle-by-cycle calculations are prohibitive in mid or high cycle fatigue. The finite element codes are available to download from www.empaneda.com/codes .

In recent years, the pre-notched or pre-cracked small punch test (P-SPT) has been successfully used to estimate the fracture properties of metallic materials for cases in which there is not sufficient material to identify these properties from standard tests, such as CT or SENB specimens. The P-SPT basically consists of deforming a pre-notched miniature specimen, whose edges are firmly gripped by a die, using a high strength punch. The novelty of this paper lies in the estimation of fracture properties using dog-bone-shaped specimens with different confinement levels. With these specimens, three confinement variations have been studied. The results obtained enable the establishment of a variation of fracture properties depending on the level of confinement of each miniature specimen and selection of the most appropriate confinement for this goal., Ministry of Economy and Competitiveness of Spain through grant MAT2014-58738-C3

We develop a strain gradient plasticity formulation for composite materials with spatially varying volume fractions to characterize size effects in functionally graded materials (FGMs). The model is grounded on the mechanism-based strain gradient plasticity theory and effective properties are determined by means of a linear homogenization scheme. Several paradigmatic boundary value problems are numerically investigated to gain insight into the strengthening effects associated with plastic strain gradients and geometrically necessary dislocations (GNDs). The analysis of bending in micro-size functionally graded foils shows a notably stiffer response with diminishing thickness. Micro-hardness measurements from indentation reveal a significant increase with decreasing indenter size. And large dislocation densities in the vicinity of the crack substantially elevate stresses in cracked FGM components. We comprehensively assess the influence of the length scale parameter and material gradation profile to accurately characterize the micro-scale response and identify regimes of GNDs relevance in FGMs., Comment: Composite Structures (2018)

© 2016 Acta Materialia Inc. Finite element analysis of stress about a blunt crack tip, emphasizing finite strain and phenomenological and mechanism-based strain gradient plasticity (SGP) formulations, is integrated with electrochemical assessment of occluded-crack tip hydrogen (H) solubility and two H-decohesion models to predict hydrogen environment assisted crack growth properties. SGP elevates crack tip geometrically necessary dislocation density and flow stress, with enhancement declining with increasing alloy strength. Elevated hydrostatic stress promotes high-trapped H concentration for crack tip damage; it is imperative to account for SGP in H cracking models. Predictions of the threshold stress intensity factor and H-diffusion limited Stage II crack growth rate agree with experimental data for a high strength austenitic Ni-Cu superalloy (Monel®K-500) and two modern ultra-high strength martensitic steels (AerMet™100 and Ferrium™M54) stressed in 0.6 M NaCl solution over a range of applied potential. For Monel®K-500, KTHis accurately predicted versus cathodic potential using either classical or gradient-modified formulations; however, Stage II growth rate is best predicted by a SGP description of crack tip stress that justifies a critical distance of 1 μm. For steel, threshold and growth rate are best predicted using high-hydrostatic stress that exceeds 6 to 8 times alloy yield strength and extends 1 μm ahead of the crack tip. This stress is nearly achieved with a three-length phenomenological SGP formulation, but additional stress enhancement is needed, perhaps due to tip geometry or slip-microstructure.

The sensitivity of crack growth resistance to the choice of isotropic or kinematic hardening is investigated. Monotonic mode I crack advance under small scale yielding conditions is modelled via a cohesive zone formulation endowed with a traction-separation law. R-curves are computed for materials that exhibit linear or power law hardening. Kinematic hardening leads to an enhanced crack growth resistance relative to isotropic hardening. Moreover, kinematic hardening requires greater crack extension to achieve the steady state. These differences are traced to the non-proportional loading of material elements near the crack tip as the crack advances. The sensitivity of the R-curve to the cohesive zone properties and to the level of material strain hardening is explored for both isotropic and kinematic hardening., Journal of Applied Mechanics (2018)

We investigate the influence of gradient-enhanced dislocation hardening on the mechanics of notch-induced failure. The role of geometrically necessary dislocations (GNDs) in enhancing cracking is assessed by means of a mechanism-based strain gradient plasticity theory. Both stationary and propagating cracks from notch-like defects are investigated through the finite element method. A cohesive zone formulation incorporating monotonic and cyclic damage contributions is employed to address both loading conditions. Computations are performed for a very wide range of length scale parameters and numerous geometries are addressed, covering the main types of notches. Results reveal a strong influence of the plastic strain gradients in all the scenarios considered. Transitional combinations of notch angle, radius and length scale parameter are identified that establish the regimes of GNDs-relevance, laying the foundations for the rational application of gradient plasticity models in damage assessment of notched components.

Ductile damage modeling within the Small Punch Test (SPT) is extensively investigated. The capabilities of the SPT to reliably estimate fracture and damage properties are thoroughly discussed and emphasis is placed on the use of notched specimens. First, different notch profiles are analyzed and constraint conditions quantified. The role of the notch shape is comprehensively examined from both triaxiality and notch fabrication perspectives. Afterwards, a methodology is presented to extract the micromechanical-based ductile damage parameters from the load-displacement curve of notched SPT samples. Furthermore, Gurson-Tvergaard-Needleman model predictions from a top-down approach are employed to gain insight into the mechanisms governing crack initiation and subsequent propagation in small punch experiments. An accurate assessment of micromechanical toughness parameters from the SPT is of tremendous relevance when little material is available., The authors gratefully acknowledge financial support from the Ministry of Economy and Competitiveness of Spain through grant MAT2014-58738-C3.

© 2017 Elsevier Ltd A suitable piece of software is presented to connect Abaqus, a sophisticated finite element package, with Matlab, the most comprehensive program for mathematical analysis. This interface between these well-known codes not only benefits from the image processing and the integrated graph-plotting features of Matlab but also opens up new opportunities in results post-processing, statistical analysis and mathematical optimization, among many other possibilities. The software architecture and usage are appropriately described and two problems of particular engineering significance are addressed to demonstrate its capabilities. Firstly, the software is employed to assess cleavage fracture through a novel 3-parameter Weibull probabilistic framework. Then, its potential to create and train neural networks is used to identify damage parameters through a hybrid experimental–numerical scheme, and model crack propagation in structural materials by means of a cohesive zone approach. The source code, detailed documentation and a large number of tutorials can be freely downloaded from www.abaqus2matlab.com.

We present a compelling finite element framework to model hydrogen assisted fatigue by means of a hydrogen- and cycle-dependent cohesive zone formulation. The model builds upon: (i) appropriate environmental boundary conditions, (ii) a coupled mechanical and hydrogen diffusion response, driven by chemical potential gradients, (iii) a mechanical behavior characterized by finite deformation J2 plasticity, (iv) a phenomenological trapping model, (v) an irreversible cohesive zone formulation for fatigue, grounded on continuum damage mechanics, and (vi) a traction-separation law dependent on hydrogen coverage calculated from first principles. The computations show that the present scheme appropriately captures the main experimental trends; namely, the sensitivity of fatigue crack growth rates to the loading frequency and the environment. The role of yield strength, work hardening, and constraint conditions in enhancing crack growth rates as a function of the frequency is thoroughly investigated. The results reveal the need to incorporate additional sources of stress elevation, such as gradient-enhanced dislocation hardening, to attain a quantitative agreement with the experiments.

© 2016 Elsevier Ltd A novel general purpose Finite Element framework is presented to study small-scale metal plasticity. A distinct feature of the adopted distortion gradient plasticity formulation, with respect to strain gradient plasticity theories, is the constitutive inclusion of the plastic spin, as proposed by Gurtin (2004) through the prescription of a free energy dependent on Nye's dislocation density tensor. The proposed numerical scheme is developed by following and extending the mathematical principles established by Fleck and Willis (2009). The modeling of thin metallic foils under bending reveals a significant influence of the plastic shear strain and spin due to a mechanism associated with the higher-order boundary conditions allowing dislocations to exit the body. This mechanism leads to an unexpected mechanical response in terms of bending moment versus curvature, dependent on the foil length, if either viscoplasticity or isotropic hardening are included in the model. In order to study the effect of dissipative higher-order stresses, the mechanical response under non-proportional loading is also investigated., Dr. Andrea Panteghini and Prof. Samuel Forest are acknowledged for helpful discussions. The authors gratefully acknowledge financial support from the Danish Council for Independent Research under the research career programme Sapere Aude in the project “Higher Order Theories in Solid Mechanics”. E. Martínez-Pañeda also acknowledges financial support from the Ministry of Science and Innovation of Spain through grant MAT2011-28796-CO3-03, and the University of Oviedo through grant UNOV-13-PF and an excellence mobility grant within the International Campus of Excellence programme. L. Bardella additionally acknowledges financial support from the Italian Ministry of Education, University, and Research (MIUR).

In this work hydrogen diffusion towards the fracture process zone is examined accounting for local hardening due to geometrically necessary dislocations (GNDs) by means of strain gradient plasticity (SGP). Finite element computations are performed within the finite deformation theory to characterize the gradient-enhanced stress elevation and subsequent diffusion of hydrogen towards the crack tip. Results reveal that GNDs, absent in conventional plasticity predictions, play a fundamental role on hydrogen transport ahead of a crack. SGP estimations provide a good agreement with experimental measurements of crack tip deformation and high levels of lattice hydrogen concentration are predicted within microns to the crack tip. The important implications of the results in the understanding of hydrogen embrittlement mechanisms are thoroughly discussed. (C) 2016 Hydrogen Energy Publications LLC. Published by Elsevier Ltd. All rights reserved.

All rights reserved. In this work a general framework for damage and fracture assessment including the effect of strain gradients is provided. Both mechanism-based and phenomenological strain gradient plasticity (SGP) theories are implemented numerically using finite deformation theory and crack tip fields are investigated. Differences and similarities between the two approaches within continuum SGP modeling are highlighted and discussed. Local strain hardening promoted by geometrically necessary dislocations (GNDs) in the vicinity of the crack leads to much higher stresses, relative to classical plasticity predictions. These differences increase significantly when large strains are taken into account, as a consequence of the contribution of strain gradients to the work hardening of the material. The magnitude of stress elevation at the crack tip and the distance ahead of the crack where GNDs significantly alter the stress distributions are quantified. The SGP dominated zone extends over meaningful physical lengths that could embrace the critical distance of several damage mechanisms, being particularly relevant for hydrogen assisted cracking models. A major role of a certain length parameter is observed in the multiple parameter version of the phenomenological SGP theory. Since this also dominates the mechanics of indentation testing, results suggest that length parameters characteristic of mode I fracture should be inferred from nanoindentation.