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;
};