This code calculates and plots schmid factors, plane traces and opening angles in a deformed sample Load your EBSD file in .txt format from TSL OIM or other Pick a grain and a lattice (e.g. BCC, FCC, HCP not implemented yet!), the programm will:
- compute all Schmid factors, plot and save them
- compute all plane traces, plot and save them
- compute all opening angles in the reference frame of the slip band, and save them Conventions: loading direction = horizontal (message me if you'd like me to make it a parameter of the code) Euler angles = Bunge's convention (angles are all in degrees, though there's a function in crystallography.py you can change to convert them if needed, message me if you'd like this to be a parameter in the code)
Download both files and change the path to your grain file. You want to edit the analyzeGrain.py file only:
- 1/ pick your lattice
- 2/ pick your grain # of interest
- 3/ hit 'run' and check the results!
Here's an example of the code's output, for a grain of Euler angles (194.3,83.2,68.1) in a BCC lattice, grain '48' in 'grainFile2.txt'
Schmid factors for {110} planes:
{112} planes:
{123} planes:
{134} planes:

Traces of {110} planes:
{112} planes:
{123} planes:
{134} planes:

Opening ('gamma') angles are to be considered when you are using a discontinuity-tolerant DIC code (e.g. Heaviside DIC, as in F. Bourdin, JC. Stinville et al., Measurements of Plastic Localization by Heaviside-Digital Image Correlation, Acta Materialia 157, pp.307:325 (2018), doi: 10.1016/j.actamat.2018.07.013, https://www.researchgate.net/publication/326406014_Measurements_of_Plastic_Localization_by_Heaviside-Digital_Image_Correlation).
This angle is calculated in the reference frame of the slip band and enables to determine the slip direction among discrete possibilities, once the slip plane is known. See above publication for more details, as well as MA. Charpagne, JC. Stinville et al. Materials Characterization (2020), details to come.
Convention for the gamma angle:
(see above publication for more details).
Gamma angles are saved in a .txt file when you run the code.