Formation of grain boundaries in ductile single crystals under plane-strain simple shear: MATLAB implementation
Authors/Creators
-
1.
Ton Duc Thang University
- 2. Institute for Advanced Study in Technology
- 3. Faculty of Mathematics and Statistics
Description
MATLAB implementation and the scripts that reproduce every figure and table of the paper "Formation of grain boundaries in ductile single crystals under plane-strain simple shear: a block-coordinate finite element method" by K. C. Le and T. D. Nguyen.
The code minimises the continuum-dislocation energy over the deformation and the plastic slip by a block-coordinate scheme: an ADMM solve for the non-smooth slip block and a Levenberg-regularised Newton solve for the deformation block, with a polyconvex Ciarlet-Geymonat elastic energy.
Running reproduce_all.m regenerates every figure and prints every table in a few seconds from the saved results. Twelve further scripts verify the propositions and the algebra of the paper, each printing PASS or FAIL.
Requires MATLAB R2025b; no toolboxes.