This website, while active, is no longer being frequently updated. For the most recent information visit our Oxford website: https://mechmat.web.ox.ac.uk/
Please find below the source files and documentation for most of the codes tat we have developed during our research. Please remember to cite the corresponding paper(s) if you use any of these codes for research or industrial purposes. Do not hesitate to contact us if questions arise.
Codes available for download
- Julia code for phase field fracture and fatigue
- MATLAB code for electrochemical phase field hydrogen fractures
- COMSOL phase field code to predict material biodegradation
- COMSOL phase field code to predict concrete corrosion-induced cracking
- MATLAB code to predict corrosion under charge conservation conditions
- COMSOL implementation of electro-chemo-mechanical phase field corrosion
- Abaqus UEL subroutine for piezoresistive fracture
- MATLAB code to predict H uptake using lumped integration
- COMSOL-MATLAB code to estimate H diffusivity from isothermal TDS
- Abaqus UMAT for phase field fracture with a Drucker-Prager criterion
- COMSOL module to predict hydrogen uptake in metals
- BrazVal: A MATLAB App to assess the validity of the Brazilian test
- COMSOL model for phase field fatigue in Li-Ion battery materials
- Abaqus UMAT for phase field fracture
- Abaqus UEL subroutine for phase field corrosion
- Abaqus UMATHT subroutine for coupled hydrogen diffusion
- Abaqus UEL subroutine for distortion gradient plasticity
- Abaqus UEL for phase field fracture/fatigue and quasi-Newton
- Abaqus UMAT for power law hardening plasticity
- FEniCS script for phase field fracture
- Abaqus UEL subroutine for strain gradient plasticity
- Abaqus UEL subroutine for phase field fracture and hydrogen diffusion
- Abaqus UEL subroutine for phase field fracture
- Control algorithm for convergence problems in crack propagation studies
- Matlab script for hydrogen assisted cracking predictions
- Abaqus2Matlab: a suitable tool for finite element post-processing
- Abaqus USDFLD subroutine for Functionally Graded Materials
- Matlab non-linear X-FEM code for MSG plasticity
- Abaqus UEL subroutine for cohesive elements
- Abaqus UMAT subroutine for CMSG plasticity
- Abaqus UELMAT for phase field fracture
Julia code for phase field fracture and fatigue
Julia code for efficient fracture and fatigue predictions using the phase field method.
MATLAB code for electrochemical phase field hydrogen fractures
MATLAB code to simulate hydrogen embrittlement using phase field fracture and electrolyte modelling.
COMSOL phase field code to predict material biodegradation
COMSOL implementation of a phase field model for prediting biocorrosion.
COMSOL phase field code to predict concrete corrosion-induced cracking
MATLAB code to predict corrosion under charge conservation conditions
Finite element MATLAB code to simulate natural corrosion rates in the absence of external current sources.
COMSOL implementation of electro-chemo-mechanical phase field corrosion
COMSOL implementation of phase field corrosion accounting for the role of mechanics and electrolyte electrochemistry.
Abaqus UEL subroutine for piezoresistive fracture
Abaqus UEL subroutine to simulate coupled deformation-electrical-phase field fracture behaviour in CNT-based composites.
MATLAB code to predict H uptake using lumped integration
Finite element MATLAB code for electrochemical reactions, using lumped integration for efficiency and robustness, and particularised to the case of hydrogen uptake.
COMSOL-MATLAB code to estimate H diffusivity from isothermal TDS
COMSOL-MATLAB files to fit the output of an isothermal TDS experiment so as to obtain the diffusivity of hydrogen in metals.
Abaqus UMAT for phase field fracture with a Drucker-Prager criterion
ABAQUS UMAT user subroutine for implementing a generalised version of phase field fracture allowing for any fracture driving force split (including Drucker-Prager) and any choice of crack density function.
COMSOL module to predict hydrogen uptake in metals
COMSOL Physics builder to predict hydrogen uptake through an electro-chemo-mechanics framework that resolves the HER.
BrazVal: A MATLAB App to assess the validity of the Brazilian test
MATLAB App that enables assessing the validity of the Brazilian test and provides the material strength from test data.
COMSOL model for phase field fatigue in Li-Ion battery materials
Abaqus UMAT for phase field fracture
ABAQUS user material (UMAT) user subroutine for implementing phase field fracture without the need for user elements by taking advantage of the heat transfer equation.
Abaqus UEL subroutine for phase field corrosion
Abaqus UMATHT subroutine for coupled hydrogen diffusion
Abaqus UEL subroutine for distortion gradient plasticity
ABAQUS user-element subroutine (UEL) for Gurtin (2004) distortion gradient plasticity formulation, incorporating the role of the plastic spin as well as dissipative and energetic higher order contributions.
Abaqus UEL for phase field fracture/fatigue and quasi-Newton
ABAQUS user-element subroutine (UEL) with phase field fatigue (and fracture). Includes an input file for using Quasi-Newton to run in a robust monolithic manner and a document with detailed instructions.
Abaqus UMAT for power law hardening plasticity
ABAQUS user-material subroutine (UMAT) with a plane strain/3D implementation of von Mises plasticity with isotropic power law hardening. Includes documentation.
FEniCS script for phase field fracture
FEniCS Python script with a staggered implementation of the phase field fracture method, suitable for 2D and 3D case studies. Includes detailed documenation, with a line-by-line of code explanation, facilitating the learning of both FEniCS and phase field fracture methods.
Abaqus UEL subroutine for strain gradient plasticity
ABAQUS user-element subroutine (UEL) for Gudmundson (2004) higher order strain gradient plasticity formulation, including both energetic and dissipative length scales.
Abaqus UEL subroutine for phase field fracture and hydrogen diffusion
ABAQUS user-element subroutine (UEL) for the coupled deformation – hydrogen transport – phase field fracture scheme presented in the associated paper.
Abaqus UEL subroutine for phase field fracture
ABAQUS user-element subroutine (UEL) with a robust phase field formulation for fracture. Several staggered integration schemes are available, as detailed in the accompanying documentation. (See also the code above for a monolithic quasi-Newton implementation)
Control algorithm for convergence problems in crack propagation studies
ABAQUS input files with a control algorithm to overcome convergence problems in cohesive zone modelling of crack propagation. The code is largely inspired by the work by Segurado and Llorca (2004) on particle fracture in composites.
Matlab script for hydrogen assisted cracking predictions
Matlab script to predict the variation with the applied potential of the cracking threshold and the stage II crack growth rate. The model builds on Gerberich’s decohesion dislocation-based model, strain gradient plasticity and advanced electrochemistry. See details in the associated paper.
Abaqus2Matlab: a suitable tool for finite element post-processing
A toolbox to connect Abaqus and Matlab, enabling the usage of the statistical analysis, image processing, integrated graph-plotting and mathematical optimization capabilities of Matlab to post-process the outcome of advanced finite element calculations. See
Abaqus USDFLD subroutine for Functionally Graded Materials
ABAQUS USDFLD Subroutine for the implementation of a continuous variation of the material elastic properties between integration points.
Matlab non-linear X-FEM code for MSG plasticity
Abaqus UEL subroutine for cohesive elements
ABAQUS user-element subroutine (UEL) with a cohesive zone formulation. Includes the dependence of the cohesive strength on the total hydrogen concentration and the effect of cyclic loading.
Abaqus UMAT subroutine for CMSG plasticity
ABAQUS User Material (UMAT) subroutine with the constitutive formulation of the conventional mechanism-based strain gradient (CMSG) plasticity theory.
Abaqus UELMAT for phase field fracture
i, E. Martínez-Pañeda. A generalised phase field model for fatigue crack growth in elastic-plastic solids with an efficient monolithic solver. Computer Methods in Applied Mechanics and Engineering 388, 114286 (2022)Contact us
Write to e.martinez-paneda@imperial.ac.uk for establishing a collaboration or to inquire for vacancies.