FINAL ORIENTATION OF EACH GRAINS

[Anmol-Bakshi]

I am using this command to run my simulation :
docker run --rm -it --volume "G:\Anmol\DAMASK\S11_Trial-20:/wd" --env OMP_NUM_THREADS=8 damaskmultiphysics/damask-grid:latest --load tension1Z.yaml --geom scaled_simulation6_mat-id_0-9.vti --material material.yaml

How to get/extract the final orientation of each grains from the .hdf5 file?

material.yaml.txt
scaled_simulation6_mat-id_0-9.vti.txt
tension1Z.yaml.txt

[Anmol-Bakshi]

Same like this :
matplotlib inline
import damask
import numpy as np
import matplotlib.pyplot as plt

adjust to your situation, file needs to exist

result_file = 'r-value_Marc/r-value_texture.hdf5'
result = damask.Result(result_file)
result.add_strain('F_p','U')
result.add_equivalent_Mises('epsilon_U^0.0(F_p)')

I am trying to add orientation "O" in my result, but it is not getting added. How to do that?

[MarDiehl]

the add_XXX functions calculate additional data based on existing data, e.g. the strain is calculated from F or F_p which are requested in your material configuration. The orientation O is also a base output, just request it in your configuration, e.g. with output: [F, P, F_e, F_p, L_p, O] under mechanical.

[Anmol-Bakshi]

I did that :
homogenization:
fullfield:
N_constituents: 1
mechanical:
type: pass
output: ['F', 'P', 'O']
And I tried with all three mechanical solvers : spectral_basic | spectral_polarization | FEM
But when I am doing the following :
docker run --rm -it -v "G:\Anmol\DAMASK\S11_Trial-20:/wd" damaskmultiphysics/python-damask /bin/bash
root@3c43c9d1c4b3:/wd# python3
Python 3.12.4 (main, Jun 20 2024, 05:25:54) [GCC 12.2.0] on linux
Type "help", "copyright", "credits" or "license" for more information.

import damask
r = damask.Result('scaled_simulation6_mat-id_0-9_tension1Z_material.hdf5')
r

Here, 'O' is not being showed as one of the output.

Kindly guide me where am I going wrong and how to resolve this?

[eisenlohr]

The homogenization level does not have 'O' as a valid output. This can only be provided at the phase level, i.e., per component at each grid point. Hence, you would have to include 'O' in the list of outputs under

phase:
  theName:
    mechanical:
      output: ['O',...]

[Anmol-Bakshi]

Thank you for your suggestion and I did that :
homogenization:
Taylor:
N_constituents: 1
mechanical:
type: isostrain
output: ['F', 'P', 'material']
phase:
Iron:
references: ['https://en.wikipedia.org/wiki/Iron']
lattice: cI
rho: 7874.0
mechanical:
elastic:
type: Hooke
references: ['W. F. Hosford. The mechanics of crystals and textured polycrystals. Oxford University Press, New York-Oxford 1993. https://doi.org/10.1002/crat.2170290414']
C_11: 233500000000.0
C_12: 135500000000.0
C_44: 118000000000.0
plastic:
type: phenopowerlaw
output: ['xi_sl', 'gamma_sl', 'O']
N_sl: [12, 12]

But still I got only these in .hdf5 file :
increment_900 (1.0 s)
phase
Iron
mechanical
gamma_sl / 1: plastic shear
xi_sl / Pa: resistance against plastic slip
homogenization
Taylor
mechanical
F / 1: deformation gradient
P / Pa: first Piola-Kirchhoff stress
epsilon_U^0.0(F) / 1: Seth-Hill strain tensor of order 0.0 based on right stretch tensor of F (deformation gradient)
epsilon_V^0.0(F) / 1: Seth-Hill strain tensor of order 0.0 based on left stretch tensor of F (deformation gradient)
sigma / Pa: Cauchy stress calculated from P (first Piola-Kirchhoff stress) and F (deformation gradient)

Now how to resolve it?

[eisenlohr]

You put the 'O' into the plastic output list. That will not work. It needs to be (directly!) under mechanical as I stated above...

By the way, could you change your formatting to reflect proper YAML indentation when posting?