ADComputeIncrementalShellStrain.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 "DenseMatrix.h"
#include "Material.h"
#include "ADRankTwoTensorForward.h"
 
namespace libMesh
{
class QGauss;
}
 
class ADComputeIncrementalShellStrain : public Material
{
public:
  static InputParameters validParams();
 
  ADComputeIncrementalShellStrain(const InputParameters & parameters);
 
protected:
  virtual void initQpStatefulProperties() override;
  virtual void computeProperties() override;
 
  virtual void computeSolnVector();
 
  virtual void computeBMatrix();
 
  virtual void computeNodeNormal();
 
  virtual void updateGVectors() {}
 
  virtual void updatedxyz() {}
 
  virtual void computeGMatrix();
 
  virtual void computeLocalCoordinates();
 
  unsigned int _nrot;
 
  unsigned int _ndisp;
 
  std::vector<unsigned int> _rot_num;
 
  std::vector<unsigned int> _disp_num;
 
  const VariableValue & _thickness;
 
  const bool _large_strain;
 
  std::vector<ADMaterialProperty<RankTwoTensor> *> _strain_increment;
 
  std::vector<ADMaterialProperty<RankTwoTensor> *> _total_strain;
 
  std::vector<const MaterialProperty<RankTwoTensor> *> _total_strain_old;
 
  NonlinearSystemBase & _nonlinear_sys;
 
  std::vector<std::vector<unsigned int>> _soln_disp_index;
 
  std::vector<std::vector<unsigned int>> _soln_rot_index;
 
  ADDenseVector _soln_vector;
 
  ADDenseVector _strain_vector;
 
  std::vector<const Node *> _nodes;
 
  ADMaterialProperty<RealVectorValue> & _node_normal;
 
  const MaterialProperty<RealVectorValue> & _node_normal_old;
 
  std::unique_ptr<libMesh::QGauss> _t_qrule;
 
  std::vector<Point> _t_points;
 
  std::vector<Point> _2d_points;
 
  std::vector<std::vector<Real>> _dphidxi_map;
 
  std::vector<std::vector<Real>> _dphideta_map;
 
  std::vector<std::vector<Real>> _phi_map;
 
  std::vector<ADMaterialProperty<RealVectorValue> *> _dxyz_dxi;
 
  std::vector<ADMaterialProperty<RealVectorValue> *> _dxyz_deta;
 
  std::vector<ADMaterialProperty<RealVectorValue> *> _dxyz_dzeta;
 
  std::vector<const MaterialProperty<RealVectorValue> *> _dxyz_dxi_old;
 
  std::vector<const MaterialProperty<RealVectorValue> *> _dxyz_deta_old;
 
  std::vector<const MaterialProperty<RealVectorValue> *> _dxyz_dzeta_old;
 
  std::vector<ADRealVectorValue> _v1;
 
  std::vector<ADRealVectorValue> _v2;
 
  std::vector<ADMaterialProperty<DenseMatrix<Real>> *> _B;
 
  std::vector<const MaterialProperty<DenseMatrix<Real>> *> _B_old;
 
  std::vector<ADMaterialProperty<RankTwoTensor> *> _ge;
 
  std::vector<const MaterialProperty<RankTwoTensor> *> _ge_old;
 
  std::vector<ADMaterialProperty<Real> *> _J_map;
 
  std::vector<const MaterialProperty<Real> *> _J_map_old;
 
  std::vector<MaterialProperty<RankTwoTensor> *> _local_transformation_matrix;
  std::vector<const MaterialProperty<RankTwoTensor> *> _local_transformation_matrix_old;
 
  std::vector<MaterialProperty<RankTwoTensor> *> _covariant_transformation_matrix;
  std::vector<const MaterialProperty<RankTwoTensor> *> _covariant_transformation_matrix_old;
 
  std::vector<MaterialProperty<RankTwoTensor> *> _contravariant_transformation_matrix;
  std::vector<const MaterialProperty<RankTwoTensor> *> _contravariant_transformation_matrix_old;
 
  std::vector<MaterialProperty<RankTwoTensor> *> _total_global_strain;
 
  ADMaterialProperty<RankTwoTensor> * _transformation_matrix;
 
  ADRealVectorValue _x2;
  ADRealVectorValue _x3;
  const NumericVector<Number> * const & _sol;
  const NumericVector<Number> & _sol_old;
  const RealVectorValue _ref_first_local_dir;
  ADRealVectorValue _g3_a;
  ADRealVectorValue _g3_c;
  ADRealVectorValue _g3_b;
  ADRealVectorValue _g3_d;
  ADRealVectorValue _g1_a;
  ADRealVectorValue _g1_c;
  ADRealVectorValue _g2_b;
  ADRealVectorValue _g2_d;
  RankTwoTensor _unrotated_total_strain;
  ADRealVectorValue _e1;
  ADRealVectorValue _e2;
  ADRealVectorValue _e3;
};