ADInertialForceShell.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 "ADTimeKernel.h"
#include "Material.h"
#include "RankTwoTensor.h"
#include "libmesh/dense_vector.h"
 
namespace libMesh
{
class QGauss;
}
 
struct PosRotVectors
{
  std::array<ADRealVectorValue, 4> pos;
  std::array<ADRealVectorValue, 4> rot;
};
 
class ADInertialForceShell : public ADTimeKernel
{
public:
  static InputParameters validParams();
 
  ADInertialForceShell(const InputParameters & parameters);
 
protected:
  virtual ADReal computeQpResidual() override { return 0.0; };
 
  virtual void computeResidual() override;
  virtual void computeResidualsForJacobian() override;
  virtual void computeShellInertialForces(const MooseArray<ADReal> & _ad_coord,
                                          const MooseArray<ADReal> & _ad_JxW);
 
private:
  unsigned int _nrot;
 
  unsigned int _ndisp;
 
  std::vector<unsigned int> _rot_num;
 
  std::vector<unsigned int> _disp_num;
 
  std::vector<unsigned int> _vel_num;
 
  std::vector<unsigned int> _accel_num;
 
  std::vector<unsigned int> _rot_vel_num;
 
  std::vector<unsigned int> _rot_accel_num;
 
  ADRankTwoTensor _original_local_config;
 
  const unsigned int _component;
 
  PosRotVectors _vel;
 
  PosRotVectors _old_vel;
 
  PosRotVectors _accel;
 
  PosRotVectors _local_vel;
 
  PosRotVectors _local_old_vel;
 
  PosRotVectors _local_accel;
 
  std::array<ADRealVectorValue, 4> _local_force, _local_moment;
 
  std::array<ADRealVectorValue, 4> _global_force, _global_moment;
 
  std::vector<std::vector<Real>> _dphidxi_map;
 
  std::vector<std::vector<Real>> _dphideta_map;
 
  std::vector<std::vector<Real>> _phi_map;
 
  std::vector<const Node *> _nodes;
 
  std::vector<Point> _2d_points;
 
  std::vector<ADReal> _2d_weights;
 
  std::vector<ADRealVectorValue> _v1;
 
  std::vector<ADRealVectorValue> _v2;
 
  ADRealVectorValue _x2;
 
  ADRealVectorValue _x3;
 
  std::vector<ADRealVectorValue> _node_normal;
 
  std::array<ADRealVectorValue, 2> _0g1_vectors;
 
  std::array<ADRealVectorValue, 2> _0g2_vectors;
 
  std::array<ADRealVectorValue, 2> _0g3_vectors;
 
  std::array<ADRealVectorValue, 2> _0g4_vectors;
 
  const MaterialProperty<Real> & _eta;
 
  const ADMaterialProperty<RankTwoTensor> & _transformation_matrix;
 
  const ADMaterialProperty<Real> & _J_map;
 
  const ADReal _thickness;
 
  const MaterialProperty<Real> & _density;
 
  const Real _alpha;
};