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