InertialForceBeam.h

Go to the documentation of this file.

//* This file is part of the MOOSE framework
//* https://www.mooseframework.org
//*
//* 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;
 
template <>
InputParameters validParams<InertialForceBeam>();
 
class InertialForceBeam : public TimeKernel
{
public:
  static InputParameters validParams();
 
  InertialForceBeam(const InputParameters & parameters);
 
  virtual void computeResidual() override;
 
  virtual void computeJacobian() override;
 
  virtual void computeOffDiagJacobian(MooseVariableFEBase & jvar) override;
  using Kernel::computeOffDiagJacobian;
 
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;
};