Computes the inertial torque, which is density * displacement x acceleration (a cross-product is used). More...

#include <InertialTorque.h>

Inheritance diagram for InertialTorque:

Public Member Functions
	InertialTorque (const InputParameters &parameters)

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual Real	computeQpResidual () override

virtual Real	computeQpJacobian () override

virtual Real	computeQpOffDiagJacobian (unsigned int jvar) override

Private Attributes
const MaterialProperty< Real > &	_density
	density More...

const Real	_beta
	Newmark beta parameter. More...

const Real	_gamma
	Newmark gamma parameter. More...

const MaterialProperty< Real > &	_eta
	Rayleigh-damping eta parameter. More...

const Real	_alpha
	HHT alpha parameter. More...

const unsigned	_component
	Component of the cross-product desired. More...

const unsigned	_ndisp
	Number of displacement variables. This must be 3. More...

std::vector< unsigned >	_disp_num
	MOOSE internal variable numbers corresponding to the displacments. More...

std::vector< const VariableValue * >	_disp
	Displacements. More...

std::vector< const VariableValue * >	_disp_old
	Old value of displacements. More...

std::vector< const VariableValue * >	_vel_old
	Old value of velocities. More...

std::vector< const VariableValue * >	_accel_old
	Old value of accelerations. More...

std::vector< Real >	_accel
	Acceleration (instantiating this vector avoids re-creating a new vector every residual calculation) More...

std::vector< Real >	_vel
	Velocity (instantiating this vector avoids re-creating a new vector every residual calculation) More...

std::vector< Real >	_daccel
	Derivative of acceleration with respect to displacement. More...

std::vector< Real >	_dvel
	Derivative of velocity with respect to displacement. More...

Detailed Description

Computes the inertial torque, which is density * displacement x acceleration (a cross-product is used).

Newmark time integration is used, and parameters are included that allow Rayleigh damping and HHT time integration

Definition at line 28 of file InertialTorque.h.

Constructor & Destructor Documentation

◆ InertialTorque()

InertialTorque::InertialTorque ( const InputParameters & parameters )

Definition at line 49 of file InertialTorque.C.

   : TimeKernel(parameters),
     _density(getMaterialProperty<Real>("density")),
     _beta(getParam<Real>("beta")),
     _gamma(getParam<Real>("gamma")),
     _eta(getMaterialProperty<Real>("eta")),
     _alpha(getParam<Real>("alpha")),
     _component(getParam<unsigned int>("component")),
     _ndisp(coupledComponents("displacements")),
     _disp_num(_ndisp),
     _disp(_ndisp),
     _disp_old(_ndisp),
     _vel_old(_ndisp),
     _accel_old(_ndisp),
     _accel(_ndisp),
     _vel(_ndisp),
     _daccel(_ndisp),
     _dvel(_ndisp)
 {
   if (_ndisp != 3 || coupledComponents("velocities") != 3 ||
       coupledComponents("accelerations") != 3)
     mooseError(
         "InertialTorque: This Kernel is only defined for 3-dimensional simulations so 3 "
         "displacement variables, 3 velocity variables and 3 acceleration variables are needed");
   for (unsigned i = 0; i < _ndisp; ++i)
   {
     _disp_num[i] = coupled("displacements", i);
     _disp[i] = &coupledValue("displacements", i);
     _disp_old[i] = &coupledValueOld("displacements", i);
     _vel_old[i] = &coupledValueOld("velocities", i);
     _accel_old[i] = &coupledValueOld("accelerations", i);
   }
 }

Member Function Documentation

◆ computeQpJacobian()

Real InertialTorque::computeQpJacobian ( )

overrideprotectedvirtual

Definition at line 106 of file InertialTorque.C.

 {
   return InertialTorque::computeQpOffDiagJacobian(_var.number());
 }

◆ computeQpOffDiagJacobian()

Real InertialTorque::computeQpOffDiagJacobian ( unsigned int jvar )

overrideprotectedvirtual

Definition at line 112 of file InertialTorque.C.

 {
   for (unsigned int j = 0; j < _ndisp; ++j)
   {
     _accel[j] = 1.0 / _beta *
                 ((((*_disp[j])[_qp] - (*_disp_old[j])[_qp]) / (_dt * _dt)) -
                  (*_vel_old[j])[_qp] / _dt - (*_accel_old[j])[_qp] * (0.5 - _beta));
     _daccel[j] = 1.0 / _beta / _dt / _dt;
     _vel[j] = (*_vel_old[j])[_qp] + (_dt * (1 - _gamma)) * (*_accel_old[j])[_qp] +
               _gamma * _dt * _accel[j];
     _dvel[j] = _gamma * _dt * _daccel[j];
     _accel[j] =
         _accel[j] + _vel[j] * _eta[_qp] * (1 + _alpha) - _alpha * _eta[_qp] * (*_vel_old[j])[_qp];
     _daccel[j] = _daccel[j] + _dvel[j] * _eta[_qp] * (1 + _alpha);
   }
  
   for (unsigned v = 0; v < _ndisp; ++v)
   {
     if (v == _component)
       continue; // PermutationTensor will kill this
     if (jvar == _disp_num[v])
     {
       Real the_sum = 0.0;
       for (unsigned int k = 0; k < _ndisp; ++k)
         the_sum += PermutationTensor::eps(_component, v, k) * _accel[k];
       for (unsigned int j = 0; j < _ndisp; ++j)
         the_sum += PermutationTensor::eps(_component, j, v) * (*_disp[j])[_qp] * _daccel[v];
       return _test[_i][_qp] * _density[_qp] * the_sum * _phi[_j][_qp];
     }
   }
   return 0.0;
 }

Referenced by computeQpJacobian().

◆ computeQpResidual()

Real InertialTorque::computeQpResidual ( )

overrideprotectedvirtual

Definition at line 84 of file InertialTorque.C.

 {
   for (unsigned int j = 0; j < _ndisp; ++j)
   {
     // Newmark and damping
     _accel[j] = 1.0 / _beta *
                 ((((*_disp[j])[_qp] - (*_disp_old[j])[_qp]) / (_dt * _dt)) -
                  (*_vel_old[j])[_qp] / _dt - (*_accel_old[j])[_qp] * (0.5 - _beta));
     _vel[j] = (*_vel_old[j])[_qp] + (_dt * (1 - _gamma)) * (*_accel_old[j])[_qp] +
               _gamma * _dt * _accel[j];
     _accel[j] =
         _accel[j] + _vel[j] * _eta[_qp] * (1 + _alpha) - _alpha * _eta[_qp] * (*_vel_old[j])[_qp];
   }
  
   Real the_sum = 0.0;
   for (unsigned int j = 0; j < _ndisp; ++j)
     for (unsigned int k = 0; k < _ndisp; ++k)
       the_sum += PermutationTensor::eps(_component, j, k) * (*_disp[j])[_qp] * _accel[k];
   return _test[_i][_qp] * _density[_qp] * the_sum;
 }

◆ validParams()

InputParameters InertialTorque::validParams ( )

static

Definition at line 20 of file InertialTorque.C.

 {
   InputParameters params = TimeKernel::validParams();
   params.addClassDescription("Kernel for inertial torque: density * displacement x acceleration");
   params.set<bool>("use_displaced_mesh") = false;
   params.addRequiredRangeCheckedParam<unsigned int>(
       "component",
       "component<3",
       "The component of the Cosserat rotation Variable prescribed to "
       "this Kernel (0 for x, 1 for y, 2 for z)");
   params.addRequiredCoupledVar("displacements", "The 3 displacement variables");
   params.addRequiredCoupledVar("velocities", "The 3 velocity variables");
   params.addRequiredCoupledVar("accelerations", "The 3 acceleration variables");
   params.addRequiredParam<Real>("beta", "beta parameter for Newmark Time integration");
   params.addRequiredParam<Real>("gamma", "gamma parameter for Newmark Time integration");
   params.addParam<MaterialPropertyName>(
       "density", "density", "Name of Material Property that provides the density");
   params.addParam<MaterialPropertyName>("eta",
                                         0.0,
                                         "Name of material property or a constant real "
                                         "number defining the eta parameter for the "
                                         "Rayleigh damping.");
   params.addParam<Real>("alpha",
                         0,
                         "alpha parameter for mass dependent numerical damping induced "
                         "by HHT time integration scheme");
   return params;
 }

Member Data Documentation

◆ _accel

std::vector<Real> InertialTorque::_accel

private

Acceleration (instantiating this vector avoids re-creating a new vector every residual calculation)

Definition at line 82 of file InertialTorque.h.

Referenced by computeQpOffDiagJacobian(), and computeQpResidual().

◆ _accel_old

std::vector<const VariableValue *> InertialTorque::_accel_old

private

Old value of accelerations.

Definition at line 79 of file InertialTorque.h.

Referenced by computeQpOffDiagJacobian(), computeQpResidual(), and InertialTorque().

◆ _alpha

const Real InertialTorque::_alpha

private

HHT alpha parameter.

Definition at line 54 of file InertialTorque.h.

Referenced by computeQpOffDiagJacobian(), and computeQpResidual().

◆ _beta

const Real InertialTorque::_beta

private

Newmark beta parameter.

Definition at line 45 of file InertialTorque.h.

Referenced by computeQpOffDiagJacobian(), and computeQpResidual().

◆ _component

const unsigned InertialTorque::_component

private

Component of the cross-product desired.

This kernel will calculate the _component component of density * displacement x acceleration

Definition at line 61 of file InertialTorque.h.

Referenced by computeQpOffDiagJacobian(), and computeQpResidual().

◆ _daccel

std::vector<Real> InertialTorque::_daccel

private

Derivative of acceleration with respect to displacement.

Definition at line 88 of file InertialTorque.h.

Referenced by computeQpOffDiagJacobian().

◆ _density

const MaterialProperty<Real>& InertialTorque::_density

private

density

Definition at line 42 of file InertialTorque.h.

Referenced by computeQpOffDiagJacobian(), and computeQpResidual().

◆ _disp

std::vector<const VariableValue *> InertialTorque::_disp

private

Displacements.

Definition at line 70 of file InertialTorque.h.

Referenced by computeQpOffDiagJacobian(), computeQpResidual(), and InertialTorque().

◆ _disp_num

std::vector<unsigned> InertialTorque::_disp_num

private

MOOSE internal variable numbers corresponding to the displacments.

Definition at line 67 of file InertialTorque.h.

Referenced by computeQpOffDiagJacobian(), and InertialTorque().

◆ _disp_old

std::vector<const VariableValue *> InertialTorque::_disp_old

private

Old value of displacements.

Definition at line 73 of file InertialTorque.h.

Referenced by computeQpOffDiagJacobian(), computeQpResidual(), and InertialTorque().

◆ _dvel

std::vector<Real> InertialTorque::_dvel

private

Derivative of velocity with respect to displacement.

Definition at line 91 of file InertialTorque.h.

Referenced by computeQpOffDiagJacobian().

◆ _eta

const MaterialProperty<Real>& InertialTorque::_eta

private

Rayleigh-damping eta parameter.

Definition at line 51 of file InertialTorque.h.

Referenced by computeQpOffDiagJacobian(), and computeQpResidual().

◆ _gamma

const Real InertialTorque::_gamma

private

Newmark gamma parameter.

Definition at line 48 of file InertialTorque.h.

Referenced by computeQpOffDiagJacobian(), and computeQpResidual().

◆ _ndisp

const unsigned InertialTorque::_ndisp

private

Number of displacement variables. This must be 3.

Definition at line 64 of file InertialTorque.h.

Referenced by computeQpOffDiagJacobian(), computeQpResidual(), and InertialTorque().

◆ _vel

std::vector<Real> InertialTorque::_vel

private

Velocity (instantiating this vector avoids re-creating a new vector every residual calculation)

Definition at line 85 of file InertialTorque.h.

Referenced by computeQpOffDiagJacobian(), and computeQpResidual().

◆ _vel_old

std::vector<const VariableValue *> InertialTorque::_vel_old

private

Old value of velocities.

Definition at line 76 of file InertialTorque.h.

Referenced by computeQpOffDiagJacobian(), computeQpResidual(), and InertialTorque().

The documentation for this class was generated from the following files:

tensor_mechanics/include/kernels/InertialTorque.h
tensor_mechanics/src/kernels/InertialTorque.C

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ InertialTorque()

Member Function Documentation

◆ computeQpJacobian()

◆ computeQpOffDiagJacobian()

◆ computeQpResidual()

◆ validParams()

Member Data Documentation

◆ _accel

◆ _accel_old

◆ _alpha

◆ _beta

◆ _component

◆ _daccel

◆ _density

◆ _disp

◆ _disp_num

◆ _disp_old

◆ _dvel

◆ _eta

◆ _gamma

◆ _ndisp

◆ _vel

◆ _vel_old