ComputeGBMisorientationType.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://mooseframework.inl.gov
//*
//* All rights reserved, see COPYRIGHT for full restrictions
//* https://github.com/idaholab/moose/blob/master/COPYRIGHT
//*
//* Licensed under LGPL 2.1, please see LICENSE for details
//* https://www.gnu.org/licenses/lgpl-2.1.html
 
#pragma once
 
#include "Material.h"
#include "GrainTracker.h"
#include "EBSDReader.h"
 
class ComputeGBMisorientationType : public Material
{
public:
  static InputParameters validParams();
 
  ComputeGBMisorientationType(const InputParameters & parameters);
 
protected:
  virtual void computeQpProperties() override;
 
  virtual unsigned int getTotalLineNum() const;
  virtual unsigned int getLineNum(unsigned int grain_i, unsigned int grain_j);
  virtual Real getTripleJunctionType();
  void rotationSymmetryToQuaternion(const double O[3][3], Eigen::Quaternion<Real> & q);
  void defineSymmetryOperator();
  double getMisorientationFromQuaternion(const Eigen::Quaternion<Real> qi,
                                         const Eigen::Quaternion<Real> qj);
  void getMisorientationAngles();
 
  const GrainTracker & _grain_tracker;
 
  const EBSDReader & _ebsd_reader;
 
  std::vector<double> _misorientation_angles;
 
  std::vector<unsigned int> _gb_pairs;
  std::vector<Real> _gb_op_pairs;
 
  const unsigned int _op_num;
  const std::vector<const VariableValue *> _vals;
 
  const Real _angle_threshold;
 
  std::vector<Eigen::Quaternion<Real>> _sym_quat;
  int _o_sym = 24;
  std::vector<EulerAngles> _euler_angle;
  std::vector<Eigen::Quaternion<Real>> _quat_angle;
 
  ADMaterialProperty<Real> & _gb_type;
};