Dislocation based modelling of deformation and fracture in real engineering alloys

Lead Research Organisation: University of Oxford
Department Name: Materials

Abstract

If loaded a small amount, a metal will deform elastically, returning to its original shape when the load is removed. However if the load exceeds some value, then permanent deformation occurs, known as plasticity. Plasticity is far more complex to understand than elasticity as it involves breaking lines of atomic bonds in the metal. These lines of broken atomic bonds are called dislocations. This is analogous to the motion of a caterpillar: which does not attempt to move its whole body forward simultaneously; instead it incrementally moves its body forward in a wave of motion sweeping through the caterpillar's body. Metals contains a huge number of dislocations: these lines sweep through the metal allowing atomic planes to slip over each over, causing the metal to be permanently deformed. When metal is loaded, new dislocations are nucleated and some become trapped at obstacles. However, if the load is applied too quickly or the metal is too cold, then the dislocation lines do not have time to nucleate and move: instead whole planes of atoms are ripped apart, fracturing the metal.

In a nuclear reactor, the fuel rods are cladded in a zirconium alloy: over time, hydrogen from water used to cool the fuel rods, diffuses into the zirconium and is attracted to dislocation lines and to any small cracks or notches in the metal. If the hydrogen concentration becomes too high, hydrogen atoms will clump together to form precipitates which block dislocation motion and can easily fracture.

It is this complex interaction between, dislocations, diffusion, precipitate formation and fracture which I aim to simulate on a computer. This is possible by utilising the power of modern graphics cards (developed to play video games) which allow massively parallel simulations to be performed easily and at little cost. Even then it is only possible to simulate a very small volume of material. Traditional mechanical tests (bending or compressing pieces of metal) were always performed on large specimens, several millimetres in size, meaning it was simply not possible to simulate all the dislocations in the sample explicitly.

In the last decade it has become possible to perform mechanical tests on samples that are only a few microns in size. The samples are so small, that by utilizing the power of modern graphics cards, it will be possible to simulate the experiment including every dislocation in the material explicitly, and watch how they interact with each other and with multiple precipitates. Being able to simulate an entire experiment at this level of detail is unprecedented and it will provide new insights into the details of what exactly goes on when metal deforms plastically and fractures.

The fundamental new insights gained during the project will be used to develop more accurate engineering design rules for industry and involves close collaboration with scientists and engineers at Lawrence Livermore National Laboratory in California, Imperial College London, Culham Centre for Fusion Energy in Oxfordshire, The National Physical Laboratory in Teddington and Rolls-Royce in Derby.

Planned Impact

Who will benefit from this research?
The project will provide new insight and enable engineering design rules to be developed which will benefit the immediate project partners (Rolls-Royce, the National Physics Laboratory, Imperial College London, Culham Centre for Fusion Energy and Lawrence Livermore National Laboratory) as well as researchers in nuclear materials in Oxford, in both the host Department of Materials and in the Department of Engineering Science. The modelling capability developed during the project will extend the range of problems that can be simulated with dislocation dynamics to include complex alloys and hydrogen diffusion. By publishing the code and training material online, researchers throughout the world will be able to benefit, particularly universities and research groups which currently have little or no dislocation modelling capability. Experimental researchers in micromechanics will be immediate users of these models. Researchers in materials for nuclear reactors will benefit; for example the model will be able to accurately simulate a range of alloys used in current or future reactors including: Zr alloys, Ni alloys, steels and ODS alloys. In the longer term the model could be used to simulate stress corrosion cracking or in virtual prototyping and alloy development. The model will also provide a link between the nano-scale, electronic structure and molecular dynamics simulations (for example research at the Thomas Young Centre in London) and larger scale crystal plasticity finite element modelling (e.g. Fionn Dunne's group at Imperial College, Anish Roy's at Loughborough University) and industrial simulation of whole components (e.g. at Rolls-Royce).

How will they benefit from this research?
The specific benefits to the UK are more accurate modelling capability of engineering alloys. This will enable more accurate engineering design rules to be extracted as well as constitutive laws for use in larger scale engineering models. The project aims to produce a model able to predict delayed hydride cracking and to develop design rules which Rolls-Royce will be able to use. Today dislocation models are not used in industry as they are not able to simulate alloys accurately. It is anticipated that this project will generate not only good science but also have a real and significant industrial impact. A key area of impact that this project will accomplish is the training of the next generation of skilled researchers in 3D discrete dislocation plasticity. This will be achieved through the training of a postdoctoral research assistant and doctoral students who will become knowledgeable in state of the art 3D coupled discrete dislocation and finite element modelling techniques. Interactions with undergraduates will also be developed through the undertaking of final year (Part II) projects to simulate reactor materials. The UK is lacking in these core skills and as such, this project will provide a means by which the UK can once again become a leader in the generation of new materials technology.

Publications

10 25 50
 
Description A 3D discrete dislocation plasticity code has been developed to enable more accurate simulation of finite domains. The code can be used to simulate a range of experiments: nano-indentation, micro cantilever bending, and micro pillar compressions. This allows direct comparison with and analysis of experimental data. The code is relatively easy to use and modify and is publicly available.
Exploitation Route The code could be used be other groups to simulate their tests or to perform virtual tests to understand how dislocation dynamics influences mechanical properties.
Sectors Aerospace

Defence and Marine

Energy

URL https://github.com/TarletonGroup/EasyDD
 
Description I was invited to become a member of the Rolls-Royce core materials technical focus group due to my expertise in materials modelling.
First Year Of Impact 2020
Sector Aerospace, Defence and Marine
Impact Types Economic

 
Description Joined Rolls-Royce core materials technical focus group
Geographic Reach National 
Policy Influence Type Membership of a guideline committee
 
Description Senior Research Fellowship
Amount £760,000 (GBP)
Organisation Royal Academy of Engineering 
Sector Charity/Non Profit
Country United Kingdom
Start 03/2022 
End 02/2027
 
Description Summer internship funding
Amount £2,300 (GBP)
Organisation Embassy of France, London 
Sector Public
Country United Kingdom
Start 04/2016 
End 08/2016
 
Title 3D Discrete Dislocation Plasticity Code 
Description Finite domain discrete dislocation plasticity code to allow easier simulations of micro mechanical tests such as nano-indentation, microcatilever beam bending and pillar compression. Written mainly in Matlab and based on original DDLab from Lawrence Livermore National lab. 
Type Of Material Computer model/algorithm 
Year Produced 2021 
Provided To Others? Yes  
Impact Code provided new insights into the relationship between dislocation structure evolution and mechanical properties. Several papers have resulted from the development and use of this code. 
URL https://github.com/TarletonGroup/EasyDD
 
Title Crystal Plasticity UMAT 
Description Crystal plasticity user material subroutine for Abaqus, based on the UEL by F. Dunne et al. 
Type Of Material Computer model/algorithm 
Year Produced 2021 
Provided To Others? Yes  
Impact The code has resulted in several publications and provided new insights into the plastic deformation of various materials. 
URL https://github.com/TarletonGroup/CrystalPlasticity
 
Title Synchrotron Diffraction During Stress Relaxation in CP Ti (grade 4) 
Description These data support the associated paper: Cold Creep of Titanium: Analysis of stress relaxation using synchrotron diffraction and crystal plasticity simulations
Yi Xiong, Phani Karamched, Chi-Toan Nguyen, David M Collins, Christopher M Magazzeni, Edmund Tarleton, Angus J Wilkinson
Acta Materialia (2020) vol. 199, 561-577
https://doi.org/10.1016/j.actamat.2020.08.010 The Authors' Accepted Manuscript version of the paper is available open access on arXiv:
https://arxiv.org/ftp/arxiv/papers/2003/2003.01682.pdf This dataset arises from an in situ stress relaxation experiment on commercially pure (grade 4) Ti undertaken at the Diamond Light Source, beamline ID12, as part of experiment EE17222. The sample was loaded to just beyond the yield point, and then held at constant strain for 5 minutes over which time the stress relaxed. The sample was then reloaded elastically and a further period of stress relaxation at fixed total strain undertaken. in total five periods of stress relaxation were imposed. Throughout the mechanical testing cycle powder diffraction patterns were recorded in the transmission geometry, at 1 second intervals using a 2d Pixium detector held 1097 mm from the sample. The beam energy was determined to be 79.79 keV. Diffraction Patterns are contained as 16 bit TIF files bundled into the Patterns_72995.zip file. Macroscopic mechanical test data are in the excel file MechTest_2995.xlsx. Small EBSD map in #.ctf format converted from Bruker #bcf file in CPg4.ctf file. Matlab + MTEX script to load and make simple plots from EBSD data file in CPg4Ti_EBSDmap.m file. Details of the MTEX orientation analysis package from Matlab can be found and freely downloaded at:
https://mtex-toolbox.github.io/ 
Type Of Material Database/Collection of data 
Year Produced 2020 
Provided To Others? Yes  
URL https://zenodo.org/record/3960528
 
Title Synchrotron Diffraction During Stress Relaxation in CP Ti (grade 4) 
Description These data support the associated paper: Cold Creep of Titanium: Analysis of stress relaxation using synchrotron diffraction and crystal plasticity simulations
Yi Xiong, Phani Karamched, Chi-Toan Nguyen, David M Collins, Christopher M Magazzeni, Edmund Tarleton, Angus J Wilkinson
Acta Materialia (2020) vol. 199, 561-577
https://doi.org/10.1016/j.actamat.2020.08.010 The Authors' Accepted Manuscript version of the paper is available open access on arXiv:
https://arxiv.org/ftp/arxiv/papers/2003/2003.01682.pdf This dataset arises from an in situ stress relaxation experiment on commercially pure (grade 4) Ti undertaken at the Diamond Light Source, beamline ID12, as part of experiment EE17222. The sample was loaded to just beyond the yield point, and then held at constant strain for 5 minutes over which time the stress relaxed. The sample was then reloaded elastically and a further period of stress relaxation at fixed total strain undertaken. in total five periods of stress relaxation were imposed. Throughout the mechanical testing cycle powder diffraction patterns were recorded in the transmission geometry, at 1 second intervals using a 2d Pixium detector held 1097 mm from the sample. The beam energy was determined to be 79.79 keV. Diffraction Patterns are contained as 16 bit TIF files bundled into the Patterns_72995.zip file. Macroscopic mechanical test data are in the excel file MechTest_2995.xlsx. Small EBSD map in #.ctf format converted from Bruker #bcf file in CPg4.ctf file. Matlab + MTEX script to load and make simple plots from EBSD data file in CPg4Ti_EBSDmap.m file. Details of the MTEX orientation analysis package from Matlab can be found and freely downloaded at:
https://mtex-toolbox.github.io/ 
Type Of Material Database/Collection of data 
Year Produced 2020 
Provided To Others? Yes  
URL https://zenodo.org/record/3960529
 
Description Summer School (Hong Kong) 
Form Of Engagement Activity Participation in an activity, workshop or similar
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Schools
Results and Impact 100 (15-16 year olds) attended classes on Materials Engineering science over a week in Hong Kong. Age appropriate materials was developed including, hands on activities (model bridge building, building a speaker and a Chladni plate demonstration). This inspired and introduced a subject outside of the standard school curriculum and gave students the opportunity for discussion. The course was the third most popular out of eight subjects.
Year(s) Of Engagement Activity 2016,2018
URL http://www.st-annes.ox.ac.uk/current/st-annes-international/summer-schools/2016-hong-kong
 
Description TMS conference (California) 
Form Of Engagement Activity A talk or presentation
Part Of Official Scheme? No
Geographic Reach International
Primary Audience Postgraduate students
Results and Impact I gave a talk at the international conference TMS showing some preliminary results, to raise awareness of the new modelling techniques i am developing. Approximately 50 people attended and this has produced a new collaboration with an experimental researcher.
Year(s) Of Engagement Activity 2017
URL http://www.tms.org/tms2017