#include <InertialForce.h>

Inheritance diagram for InertialForce:

Public Member Functions
	InertialForce (const InputParameters &parameters)

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual Real	computeQpResidual ()

virtual Real	computeQpJacobian ()

Private Attributes
const MaterialProperty< Real > &	_density

const VariableValue *	_u_old

const VariableValue *	_vel_old

const VariableValue *	_accel_old

const bool	_has_beta

const bool	_has_gamma

const Real	_beta

const Real	_gamma

const bool	_has_velocity

const bool	_has_acceleration

const MaterialProperty< Real > &	_eta

const Real	_alpha

const VariableValue *	_u_dot

const VariableValue *	_u_dotdot

const VariableValue *	_u_dot_old

const VariableValue *	_du_dot_du

const VariableValue *	_du_dotdot_du

Detailed Description

Definition at line 21 of file InertialForce.h.

Constructor & Destructor Documentation

◆ InertialForce()

InertialForce::InertialForce ( const InputParameters & parameters )

Definition at line 45 of file InertialForce.C.

   : TimeKernel(parameters),
     _density(getMaterialProperty<Real>("density")),
     _has_beta(isParamValid("beta")),
     _has_gamma(isParamValid("gamma")),
     _beta(_has_beta ? getParam<Real>("beta") : 0.1),
     _gamma(_has_gamma ? getParam<Real>("gamma") : 0.1),
     _has_velocity(isParamValid("velocity")),
     _has_acceleration(isParamValid("acceleration")),
     _eta(getMaterialProperty<Real>("eta")),
     _alpha(getParam<Real>("alpha"))
 {
   if (_has_beta && _has_gamma && _has_velocity && _has_acceleration)
   {
     _vel_old = &coupledValueOld("velocity");
     _accel_old = &coupledValueOld("acceleration");
     _u_old = &valueOld();
   }
   else if (!_has_beta && !_has_gamma && !_has_velocity && !_has_acceleration)
   {
     _u_dot = &(_var.uDot());
     _u_dotdot = &(_var.uDotDot());
     _u_dot_old = &(_var.uDotOld());
     _du_dot_du = &(_var.duDotDu());
     _du_dotdot_du = &(_var.duDotDotDu());
   }
   else
     mooseError(
         "InertialForce: Either all or none of `beta`, `gamma`, `velocity`, and `acceleration` "
         "should be provided as input.");
 }

Member Function Documentation

◆ computeQpJacobian()

Real InertialForce::computeQpJacobian ( )

protectedvirtual

Definition at line 98 of file InertialForce.C.

 {
   if (_dt == 0)
     return 0;
   else if (_has_beta)
     return _test[_i][_qp] * _density[_qp] / (_beta * _dt * _dt) * _phi[_j][_qp] +
            _eta[_qp] * (1 + _alpha) * _test[_i][_qp] * _density[_qp] * _gamma / _beta / _dt *
                _phi[_j][_qp];
   else
     return _test[_i][_qp] * _density[_qp] * (*_du_dotdot_du)[_qp] * _phi[_j][_qp] +
            _eta[_qp] * (1 + _alpha) * _test[_i][_qp] * _density[_qp] * (*_du_dot_du)[_qp] *
                _phi[_j][_qp];
 }

◆ computeQpResidual()

Real InertialForce::computeQpResidual ( )

protectedvirtual

Definition at line 78 of file InertialForce.C.

 {
   if (_dt == 0)
     return 0;
   else if (_has_beta)
   {
     Real accel = 1. / _beta *
                  (((_u[_qp] - (*_u_old)[_qp]) / (_dt * _dt)) - (*_vel_old)[_qp] / _dt -
                   (*_accel_old)[_qp] * (0.5 - _beta));
     Real vel = (*_vel_old)[_qp] + (_dt * (1 - _gamma)) * (*_accel_old)[_qp] + _gamma * _dt * accel;
     return _test[_i][_qp] * _density[_qp] *
            (accel + vel * _eta[_qp] * (1 + _alpha) - _alpha * _eta[_qp] * (*_vel_old)[_qp]);
   }
   else
     return _test[_i][_qp] * _density[_qp] *
            ((*_u_dotdot)[_qp] + (*_u_dot)[_qp] * _eta[_qp] * (1.0 + _alpha) -
             _alpha * _eta[_qp] * (*_u_dot_old)[_qp]);
 }

◆ validParams()

InputParameters InertialForce::validParams ( )

static

Definition at line 18 of file InertialForce.C.

 {
   InputParameters params = TimeKernel::validParams();
   params.addClassDescription("Calculates the residual for the inertial force "
                              "($M \\cdot acceleration$) and the contribution of mass"
                              " dependent Rayleigh damping and HHT time "
                              " integration scheme ($\\eta \\cdot M \\cdot"
                              " ((1+\\alpha)velq2-\\alpha \\cdot vel-old) $)");
   params.set<bool>("use_displaced_mesh") = true;
   params.addCoupledVar("velocity", "velocity variable");
   params.addCoupledVar("acceleration", "acceleration variable");
   params.addParam<Real>("beta", "beta parameter for Newmark Time integration");
   params.addParam<Real>("gamma", "gamma parameter for Newmark Time integration");
   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");
   params.addParam<MaterialPropertyName>(
       "density", "density", "Name of Material Property that provides the density");
   return params;
 }