Density Plot with Pandas

• Contributor: Martin Diehl (https://martin-diehl.net)

• DAMASK version: 3.0.2

• Prerequisites (data): DADF5 file with 1. Piola-Kirchhoff stress (‘P’) and deformation gradient (‘F’)

%matplotlib inline

import damask
import pandas as pd

# adjust to your situation, file needs to exist
result_file = 'densityplot/20grains16x16x16_tensionX.hdf5'
idx = [0,0] # xx-component

result = damask.Result(result_file).view(
    increments=-1,            # only view last increment
    homogenizations=False,    # do not consider homogenization results
)

df = pd.DataFrame({'F_11':result.get('F')[:,idx[0],idx[1]],
                   'P_11':result.get('P')[:,idx[0],idx[1]],})
label = f'{idx[0]+1}{idx[1]+1}'
_ = df.plot.hexbin(x=f'F_{label}',y=f'P_{label}',gridsize=20,
                   cmap=damask.Colormap.from_predefined('strain'))