InertialForceBeam.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 "TimeKernel.h"
#include "Material.h"
#include "RankTwoTensor.h"

// Forward Declarations

class InertialForceBeam : public TimeKernel
{
public:
  static InputParameters validParams();

  InertialForceBeam(const InputParameters & parameters);

  virtual void computeResidual() override;

  virtual void computeJacobian() override;

  virtual void computeOffDiagJacobian(unsigned int jvar) override;

protected:
  virtual Real computeQpResidual() override { return 0.0; };

private:
  const bool _has_beta;
  const bool _has_gamma;
  const bool _has_velocities;
  const bool _has_rot_velocities;
  const bool _has_accelerations;
  const bool _has_rot_accelerations;
  const bool _has_Ix;

  const MaterialProperty<Real> & _density;

  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;

  const VariableValue & _area;

  const VariableValue & _Ay;

  const VariableValue & _Az;

  const VariableValue & _Ix;

  const VariableValue & _Iy;

  const VariableValue & _Iz;

  const Real _beta;

  const Real _gamma;

  const MaterialProperty<Real> & _eta;

  const Real _alpha;

  const MaterialProperty<RankTwoTensor> & _original_local_config;

  const MaterialProperty<Real> & _original_length;

  const unsigned int _component;

  RealVectorValue _vel_old_0, _vel_old_1, _rot_vel_old_0, _rot_vel_old_1;

  RealVectorValue _vel_0, _vel_1, _rot_vel_0, _rot_vel_1;

  RealVectorValue _accel_0, _accel_1, _rot_accel_0, _rot_accel_1;

  RealVectorValue _local_vel_old_0, _local_vel_old_1, _local_rot_vel_old_0, _local_rot_vel_old_1;

  RealVectorValue _local_vel_0, _local_vel_1, _local_rot_vel_0, _local_rot_vel_1;

  RealVectorValue _local_accel_0, _local_accel_1, _local_rot_accel_0, _local_rot_accel_1;

  std::vector<RealVectorValue> _local_force, _local_moment;

  RealVectorValue _global_force_0, _global_force_1, _global_moment_0, _global_moment_1;

  const VariableValue * _du_dot_du;

  const VariableValue * _du_dotdot_du;
};