SymmAnisotropicElasticityTensor.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* 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 "SymmElasticityTensor.h"
 
class SymmAnisotropicElasticityTensor : public SymmElasticityTensor
{
public:
  SymmAnisotropicElasticityTensor();
  SymmAnisotropicElasticityTensor(std::vector<Real> & init_list, bool all_21);
  SymmAnisotropicElasticityTensor(const SymmAnisotropicElasticityTensor & a);
  SymmAnisotropicElasticityTensor & operator=(const SymmAnisotropicElasticityTensor & a) = default;
 
  virtual ~SymmAnisotropicElasticityTensor() = default;
 
  void setFirstEulerAngle(const Real a1);
 
  void setSecondEulerAngle(const Real a2);
 
  void setThirdEulerAngle(const Real a3);
 
  void setMaterialConstantc11(const Real c11);
 
  void setMaterialConstantc12(const Real c12);
 
  void setMaterialConstantc44(const Real c44);
 
  virtual void rotate(const Real a1, const Real a2, const Real a3);
 
  Real firstEulerAngle();
 
  Real secondEulerAngle();
 
  Real thirdEulerAngle();
 
  // Debugging functions
  void show_dt_matrix();
  void show_r_matrix();
 
protected:
  DenseMatrix<Real> _dmat;        // 9 x 9 Material Matrix
  DenseMatrix<Real> _qdmat;       // Rotated Material Matrix
  DenseMatrix<Real> _dt;          // 6 x 6 Material Matrix
  DenseMatrix<Real> _qdt;         // Rotated Material Matrix
  DenseMatrix<Real> _r;           // Rotational Matrix
  DenseMatrix<Real> _q;           // Q = R (dyadic) R
  DenseMatrix<Real> _qt;          // Q Transpose
  std::vector<Real> _euler_angle; // Stores Euler angeles
 
  DenseMatrix<Real> _trans_d6_to_d9;
  // Transformation matrix from a 6 x 6 to 9 to 9
 
  DenseMatrix<Real> _trans_d9_to_d6;
  // Transformation matrix from a 9 x 9 to 6 x 6
 
  Real _c11, _c12, _c44; // Material Constants
 
  void form_r_matrix();
  void initialize_material_dt_matrix();
  void initialize_anisotropic_material_dt_matrix();
  void form_rotational_q_matrix();
  void form_transformation_t_matrix();
  void form_transformed_material_dmat_matrix();
  void form_transformed_material_dt_matrix();
  void form_rotated_material_qdmat_matrix();
 
  virtual void calculateEntries(unsigned int qp);
};