BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

E. Martínez-Pañeda

E. Martínez-Pañeda contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
cond-mat.mtrl-sciComputational Engineering, Finance, and Sciencephysics.app-phphysics.chem-phmath.NANumerical Analysisphysics.comp-phphysics.flu-dynphysics.geo-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
16works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2026arXiv

Coupled thermo-chemo-mechanical phase field-based modelling of hydrogen-assisted cracking in girth welds

A new computational framework is presented to predict the structural integrity of welds in hydrogen transmission pipelines. The framework combines: (i) a thermo-mechanical weld process model, and (ii) a coupled deformation-diffusion-fracture phase field-based model that accounts for plasticity and hydrogen trapping, considering multiple trap types, with stationary and evolving trap densities. This enables capturing, for the first time, the interplay between residual stresses, trap creation, hydrogen transport, and fracture. The computational framework is particularised and applied to the study of weld integrity in X80 pipeline steel. The focus is on girth welds, as they are more complex due to their multi-pass nature. The weld process model enables identifying the dimensions and characteristics of the three weld regions: base metal, heat-affected zone, and weld metal, and these are treated distinctively. This is followed by virtual fracture experiments, which reveal a very good agreement with laboratory studies. Then, weld pipeline integrity is assessed, estimating critical failure pressures for a wide range of scenarios. Of particular interest is to assess the structural integrity implications of welding defects present in existing natural gas pipelines under consideration for hydrogen transport: pores, lack of penetration, imperfections, lack of fusion, root contraction, and undercutting. The results obtained in hydrogen-containing environments reveal an important role of the weld microstructure and the detrimental effect of weld defects that are likely to be present in existing natural gas pipelines, as they are considered safe in gas pipeline standards.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A coupled phase field formulation for modelling fatigue cracking in lithium-ion battery electrode particles

Electrode particle cracking is one of the main phenomena driving battery capacity degradation. Recent phase field fracture studies have investigated particle cracking behaviour. However, only the beginning of life has been considered and effects such as damage accumulation have been neglected. Here, a multi-physics phase field fatigue model has been developed to study crack propagation in battery electrode particles undergoing hundreds of cycles. In addition, we couple our electrochemo-mechanical formulation with X-ray CT imaging to simulate fatigue cracking of realistic particle microstructures. Using this modelling framework, non-linear crack propagation behaviour is predicted, leading to the observation of an exponential increase in cracked area with cycle number. Three stages of crack growth (slow, accelerating and unstable) are observed, with phenomena such as crack initialisation at concave regions and crack coalescence having a significant contribution to the resulting fatigue crack growth rates. The critical values of C-rate, particle size and initial crack length are determined, and found to be lower than those reported in the literature using static fracture models. Therefore, this work demonstrates the importance of considering fatigue damage in battery degradation models and provides insights on the control of fatigue crack propagation to alleviate battery capacity degradation.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A general framework for decomposing the phase field fracture driving force, particularised to a Drucker-Prager failure surface

Due to its computational robustness and versatility, the phase field fracture model has become the preferred tool for predicting a wide range of cracking phenomena. However, in its conventional form, its intrinsic tension-compression symmetry in damage evolution prevents its application to the modelling of compressive failures in brittle and quasi-brittle solids, such as concrete or rock materials. In this work, we present a general methodology for decomposing the phase field fracture driving force, the strain energy density, so as to reproduce asymmetrical tension-compression fracture behaviour. The generalised approach presented is particularised to the case of linear elastic solids and the Drucker-Prager failure criterion. The ability of the presented model to capture the compressive failure of brittle materials is showcased by numerically implementing the resulting strain energy split formulation and addressing four case studies of particular interest. Firstly, insight is gained into the capabilities of the model in predicting friction and dilatancy effects under shear loading. Secondly, virtual direct shear tests are conducted to assess fracture predictions under different pressure levels. Thirdly, a concrete cylinder is subjected to uniaxial and triaxial compression to investigate the influence of confinement. Finally, the localised failure of a soil slope is predicted and the results are compared with other formulations for the strain energy decomposition proposed in the literature. The results provide a good qualitative agreement with experimental observations and demonstrate the capabilities of phase field fracture methods to predict crack nucleation and growth under multi-axial loading in materials exhibiting asymmetric tension-compression fracture behaviour.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A generalised, multi-phase-field theory for dissolution-driven stress corrosion cracking and hydrogen embrittlement

We present a phase field-based electro-chemo-mechanical formulation for modelling mechanics-enhanced corrosion and hydrogen-assisted cracking in elastic-plastic solids. A multi-phase-field approach is used to present, for the first time, a general framework for stress corrosion cracking, incorporating both anodic dissolution and hydrogen embrittlement mechanisms. We numerically implement our theory using the finite element method and defining as primary kinematic variables the displacement components, the phase field corrosion order parameter, the metal ion concentration, the phase field fracture order parameter and the hydrogen concentration. Representative case studies are addressed to showcase the predictive capabilities of the model in various materials and environments, attaining a promising agreement with benchmark tests and experimental observations. We show that the generalised formulation presented can capture, as a function of the environment, the interplay between anodic dissolution- and hydrogen-driven failure mechanisms; including the transition from one to the other, their synergistic action and their individual occurrence. Such a generalised framework can bring new insight into environment-material interactions and the understanding of stress corrosion cracking, as demonstrated here by providing the first simulation results for Gruhl's seminal experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A phase field electro-chemo-mechanical formulation for predicting void evolution at the Li-electrolyte interface in all-solid-state batteries

We present a mechanistic theory for predicting void evolution in the Li metal electrode during the charge and discharge of all-solid-state battery cells. A phase field formulation is developed to model vacancy annihilation and nucleation, and to enable the tracking of the void-Li metal interface. This is coupled with a viscoplastic description of Li deformation, to capture creep effects, and a mass transfer formulation accounting for substitutional (bulk and surface) Li diffusion and current-driven flux. Moreover, we incorporate the interaction between the electrode and the solid electrolyte, resolving the coupled electro-chemical-mechanical problem in both domains. This enables predicting the electrolyte current distribution and thus the emergence of local current 'hot spots', which act as precursors for dendrite formation and cell death. The theoretical framework is numerically implemented, and single and multiple void case studies are carried out to predict the evolution of voids and current hot spots as a function of the applied pressure, material properties and charge (magnitude and cycle history). For both plating and stripping, insight is gained into the interplay between bulk diffusion, Li dissolution and deposition, creep, and the nucleation and annihilation of vacancies. The model is shown to capture the main experimental observations, including not only key features of electrolyte current and void morphology but also the sensitivity to the applied current, the role of pressure in increasing the electrode-electrolyte contact area, and the dominance of creep over vacancy diffusion.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

An assessment of phase field fracture: crack initiation and growth

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 minimisation - without the need for \textit{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.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cohesive zone modelling of hydrogen assisted fatigue crack growth: the role of trapping

We investigate the influence of microstructural traps in hydrogen-assisted fatigue crack growth. To this end, a new formulation combining multi-trap stress-assisted diffusion, mechanism-based strain gradient plasticity and a hydrogen- and fatigue-dependent cohesive zone model is presented and numerically implemented. The results show that the ratio of loading frequency to effective diffusivity governs fatigue crack growth behaviour. Increasing the density of \emph{beneficial} traps, not involved in the fracture process, results in lower fatigue crack growth rates. The combinations of loading frequency and carbide trap densities that minimise embrittlement susceptibility are identified, providing the foundation for a rational design of hydrogen-resistant alloys.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Computational Assessment of Biomass Dust Explosions in the 20L Sphere

Determination of the explosion severity parameters of biomass is crucial for the safety management and dust explosion risk assessment of biomass-processing industries. These are commonly determined following experimental tests in the 20L sphere according to the international standards. Recently, CFD simulations have emerged as a reliable alternative to predict the explosion behavior with good accuracy and reduced labor and capital. In this work, numerical simulations of biomass dust explosions are conducted with the open-source CFD code OpenFOAM. The multi-phase (gas-solid) flow is treated in an Eulerian-Lagrangian framework, using a two-way coupling regime and considering the reactions of biomass conversion (moisture evaporation, devolatilization, and char oxidation), the combustion of volatile gases, and convective and radiative heat transfer. The model is validated with pressure-time and concentration-dependent experimental measurements of two biomass samples. Results suggest that the characteristics of the cold-flow (ı.e. turbulence levels, actual dust concentration, spatial distribution of the dust cloud, and turbophoresis effect) govern the course of the explosion process, and depend strongly on particle size, dust concentration, and ignition delay time effects. These findings may be relevant in the design of better dust explosion testing devices and to the reexamination of the guidelines for the operation of the experiment. Finally, a thorough discussion on the explosion pressures, degree of biomass conversion, flame temperature, flame propagation patterns, and the dust agglomeration effect is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Computational modelling of hydrogen assisted fracture in polycrystalline materials

We present a combined phase field and cohesive zone formulation for hydrogen embrittlement that resolves the polycrystalline microstructure of metals. Unlike previous studies, our deformation-diffusion-fracture modelling framework accounts for hydrogen-microstructure interactions and explicitly captures the interplay between bulk (transgranular) fracture and intergranular fracture, with the latter being facilitated by hydrogen through mechanisms such as grain boundary decohesion. We demonstrate the potential of the theoretical and computational formulation presented by simulating inter- and trans-granular cracking in relevant case studies. Firstly, verification calculations are conducted to show how the framework predicts the expected qualitative trends. Secondly, the model is used to simulate recent experiments on pure Ni and a Ni-Cu superalloy that have attracted particular interest. We show that the model is able to provide a good quantitative agreement with testing data and yields a mechanistic rationale for the experimental observations.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Griffith-based analysis of crack initiation location in a Brazilian test

The Brazilian test has been extremely popular while prompting significant debate. The main source of controversy is rooted in its indirect nature; the material tensile strength is inferred upon assuming that cracking initiates at the centre of the sample. Here, we use the Griffith criterion and finite element analysis to map the conditions (jaws geometry and material properties) that result in the nucleation of a centre crack. Unlike previous studies, we do not restrict ourselves to evaluating the stress state at the disk centre; the failure envelope of the generalised Griffith criterion is used to establish the crack nucleation location. We find that the range of conditions where the Brazilian test is valid is much narrower than previously assumed, with current practices and standards being inappropriate for a wide range of rock-like materials. The results obtained are used to develop a protocol that experimentalists can follow to obtain a valid estimate of the material tensile strength. This is showcased with specific case studies and examples of valid and invalid tests from the literature. Furthermore, the uptake of this protocol is facilitated by providing a MATLAB App that determines the validity of the experiment for arbitrary test conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Phase field fracture predictions of microscopic bridging behaviour of composite materials

We investigate the role of microstructural bridging on the fracture toughness of composite materials. To achieve this, a new computational framework is presented that integrates phase field fracture and cohesive zone models to simulate fibre breakage, matrix cracking and fibre-matrix debonding. The composite microstructure is represented by an embedded cell at the vicinity of the crack tip, whilst the rest of the sample is modelled as an anisotropic elastic solid. The model is first validated against experimental data of transverse matrix cracking from single-notched three-point bending tests. Then, the model is extended to predict the influence of grain bridging, brick-and-mortar microstructure and 3D fibre bridging on crack growth resistance. The results show that these microstructures are very efficient in enhancing the fracture toughness via fibre-matrix debonding, fibre breakage and crack deflection. In particular, the 3D fibre bridging effect can increase the energy dissipated at failure by more than three orders of magnitude, relative to that of the bulk matrix; well in excess of the predictions obtained from the rule of mixtures. These results shed light on microscopic bridging mechanisms and provide a virtual tool for developing high fracture toughness composites.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Creep Behaviour and Tensile Response of Adhesively Bonded Polyethylene Joints: Single-Lap and Double-Strap

The static and time-dependent behaviours of adhesively bonded polyethylene Double-Strap (DS) joints were investigated to assess the viability of this joint configuration relative to the Single-Lap (SL) joints. Both experiments and finite element simulations are conducted. First, we individually characterise the tensile and creep behaviour of the adhesive and adherent materials; an epoxy-based adhesive and polyethylene, respectively. This information is used to develop suitable constitutive models that are then implemented in the commercial finite element package ABAQUS by means of user material subroutines, UMATs. The numerical models are used to design the creep tests on the adhesive joints. Afterwards, an extensive experimental campaign is conducted where we characterise the static and creep behaviour of two joint configurations, SL and DS joints, and three selected values of the overlap length. In regard to the static case, results reveal an increase in the failure load with increasing overlap length, of up to 10% for an overlap length of 39 mm. Also, slightly better performance is observed for the SL joint configuration. For the creep experiments, we show that the DS adhesive joint configuration leads to much shorter elongations, relative to the SL joints. These differences diminish with increasing overlap length but remain substantial in all cases. In both joint configurations, the elongation increases with decreasing overlap length. For instance, increasing the overlap length to 39 mm led to a 50% and a 30% reduction in elongation for SL and DS joints, respectively. Moreover, the numerical predictions show a good agreement with the experiments. The stress redistribution is investigated and it is found that the shear stress is highly sensitive to the testing time, with differences being more noticeable for the DS joint system.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Fracture in distortion gradient plasticity

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Influence of charging conditions on simulated temperature-programmed desorption for hydrogen in metals

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.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Numerical study of interface cracking in composite structures using a novel geometrically nonlinear Linear Elastic Brittle Interface Model: mixed-mode fracture conditions and application to structured interfaces

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.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Mode II fracture of an elastic-plastic sandwich layer

The shear strength of a pre-cracked sandwich layer is predicted, assuming that the layer is linear elastic or elastic-plastic, with yielding characterized by either 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 analysed: (i) a semi-infinite crack in a sandwich layer, subjected to a remote mode II K-field and (ii) a centre-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 is duly accounted for (the elastic-plastic case). Analytical and numerical solutions are then obtained for the centre-cracked sandwich plate of 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 of 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.

0 citations0 reviewsNo code linkedOpen access