InertialForceBeam Class Reference

#include <InertialForceBeam.h>

Inheritance diagram for InertialForceBeam:

Public Member Functions
	InertialForceBeam (const InputParameters &parameters)
virtual void	computeResidual () override
virtual void	computeJacobian () override
virtual void	computeOffDiagJacobian (MooseVariableFEBase &jvar) override

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual Real	computeQpResidual () override

Private Attributes
const bool	_has_beta
	Booleans for validity of params. More...
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
	Density of the beam. More...
unsigned int	_nrot
	Number of coupled rotational variables. More...
unsigned int	_ndisp
	Number of coupled displacement variables. More...
std::vector< unsigned int >	_rot_num
	Variable numbers corresponding to rotational variables. More...
std::vector< unsigned int >	_disp_num
	Variable numbers corresponding to displacement variables. More...
std::vector< unsigned int >	_vel_num
	Variable numbers corresponding to velocity aux variables. More...
std::vector< unsigned int >	_accel_num
	Variable numbers corresponding to acceleraion aux variables. More...
std::vector< unsigned int >	_rot_vel_num
	Variable numbers corresponding to rotational velocity aux variables. More...
std::vector< unsigned int >	_rot_accel_num
	Variable numbers corresponding to rotational acceleration aux variables. More...
const VariableValue &	_area
	Coupled variable for beam cross-sectional area. More...
const VariableValue &	_Ay
	Coupled variable for first moment of area of beam in y direction, i.e., integral of y*dA over the cross-section. More...
const VariableValue &	_Az
	Coupled variable for first moment of area of beam in z direction, i.e., integral of z*dA over the cross-section. More...
const VariableValue &	_Ix
	Coupled variable for second moment of area of beam in x direction, i.e., integral of (y^2+z^2)*dA over the cross-section. More...
const VariableValue &	_Iy
	Coupled variable for second moment of area of beam in y direction, i.e., integral of y^2*dA over the cross-section. More...
const VariableValue &	_Iz
	Coupled variable for second momemnt of area of beam in z direction, i.e., integral of z^2*dA over the cross-section. More...
const Real	_beta
	Newmark time integration parameter. More...
const Real	_gamma
	Newmark time integraion parameter. More...
const MaterialProperty< Real > &	_eta
	Mass proportional Rayleigh damping parameter. More...
const Real	_alpha
	HHT time integration parameter. More...
const MaterialProperty< RankTwoTensor > &	_original_local_config
	Rotational transformation from global to initial beam local coordinate system. More...
const MaterialProperty< Real > &	_original_length
	Initial length of beam. More...
const unsigned int	_component
	Direction along which residual is calculated. More...
RealVectorValue	_vel_old_0
	Old translational and rotational velocities at the two nodes of the beam in the global coordinate system. More...
RealVectorValue	_vel_old_1
RealVectorValue	_rot_vel_old_0
RealVectorValue	_rot_vel_old_1
RealVectorValue	_vel_0
	Current translational and rotational velocities at the two nodes of the beam in the global coordinate system. More...
RealVectorValue	_vel_1
RealVectorValue	_rot_vel_0
RealVectorValue	_rot_vel_1
RealVectorValue	_accel_0
	Current translational and rotational accelerations at the two nodes of the beam in the global coordinate system. More...
RealVectorValue	_accel_1
RealVectorValue	_rot_accel_0
RealVectorValue	_rot_accel_1
RealVectorValue	_local_vel_old_0
	Old translational and rotational velocities at the two nodes of the beam in the initial beam local coordinate system. More...
RealVectorValue	_local_vel_old_1
RealVectorValue	_local_rot_vel_old_0
RealVectorValue	_local_rot_vel_old_1
RealVectorValue	_local_vel_0
	Current translational and rotational velocities at the two nodes of the beam in the initial beam local coordinate system. More...
RealVectorValue	_local_vel_1
RealVectorValue	_local_rot_vel_0
RealVectorValue	_local_rot_vel_1
RealVectorValue	_local_accel_0
	Current translational and rotational accelerations at the two nodes of the beam in the initial beam local coordinate system. More...
RealVectorValue	_local_accel_1
RealVectorValue	_local_rot_accel_0
RealVectorValue	_local_rot_accel_1
std::vector< RealVectorValue >	_local_force
	Forces and moments at the two end nodes of the beam in the initial beam local configuration. More...
std::vector< RealVectorValue >	_local_moment
RealVectorValue	_global_force_0
	Forces and moments at the two end nodes of the beam in the global coordinate system. More...
RealVectorValue	_global_force_1
RealVectorValue	_global_moment_0
RealVectorValue	_global_moment_1
const VariableValue *	_du_dot_du
	Coupled variable for du_dot_du calculated by time integrator. More...
const VariableValue *	_du_dotdot_du
	Coupled variable for du_dotdot_du calculated by time integrator. More...

Detailed Description

Definition at line 22 of file InertialForceBeam.h.

Constructor & Destructor Documentation

InertialForceBeam()

InertialForceBeam::InertialForceBeam

(

const InputParameters &

parameters

)

Definition at line 75 of file InertialForceBeam.C.

  : TimeKernel(parameters),
    _has_beta(isParamValid("beta")),
    _has_gamma(isParamValid("gamma")),
    _has_velocities(isParamValid("velocities")),
    _has_rot_velocities(isParamValid("rotational_velocities")),
    _has_accelerations(isParamValid("accelerations")),
    _has_rot_accelerations(isParamValid("rotational_accelerations")),
    _has_Ix(isParamValid("Ix")),
    _density(getMaterialProperty<Real>("density")),
    _nrot(coupledComponents("rotations")),
    _ndisp(coupledComponents("displacements")),
    _rot_num(_nrot),
    _disp_num(_ndisp),
    _vel_num(_ndisp),
    _accel_num(_ndisp),
    _rot_vel_num(_nrot),
    _rot_accel_num(_nrot),
    _area(coupledValue("area")),
    _Ay(coupledValue("Ay")),
    _Az(coupledValue("Az")),
    _Ix(_has_Ix ? coupledValue("Ix") : _zero),
    _Iy(coupledValue("Iy")),
    _Iz(coupledValue("Iz")),
    _beta(_has_beta ? getParam<Real>("beta") : 0.1),
    _gamma(_has_gamma ? getParam<Real>("gamma") : 0.1),
    _eta(getMaterialProperty<Real>("eta")),
    _alpha(getParam<Real>("alpha")),
    _original_local_config(getMaterialPropertyByName<RankTwoTensor>("initial_rotation")),
    _original_length(getMaterialPropertyByName<Real>("original_length")),
    _component(getParam<unsigned int>("component")),
    _local_force(2),
    _local_moment(2)
{
  // Checking for consistency between the length of the provided rotations and displacements vector
  if (_ndisp != _nrot)
    mooseError("InertialForceBeam: The number of variables supplied in 'displacements' and "
               "'rotations' must match.");
 
  if (_has_beta && _has_gamma && _has_velocities && _has_accelerations && _has_rot_velocities &&
      _has_rot_accelerations)
  {
    if ((coupledComponents("velocities") != _ndisp) ||
        (coupledComponents("accelerations") != _ndisp) ||
        (coupledComponents("rotational_velocities") != _ndisp) ||
        (coupledComponents("rotational_accelerations") != _ndisp))
      mooseError(
          "InertialForceBeam: The number of variables supplied in 'velocities', "
          "'accelerations', 'rotational_velocities' and 'rotational_accelerations' must match "
          "the number of displacement variables.");
 
    // fetch coupled velocities and accelerations
    for (unsigned int i = 0; i < _ndisp; ++i)
    {
      MooseVariable * vel_variable = getVar("velocities", i);
      _vel_num[i] = vel_variable->number();
 
      MooseVariable * accel_variable = getVar("accelerations", i);
      _accel_num[i] = accel_variable->number();
 
      MooseVariable * rot_vel_variable = getVar("rotational_velocities", i);
      _rot_vel_num[i] = rot_vel_variable->number();
 
      MooseVariable * rot_accel_variable = getVar("rotational_accelerations", i);
      _rot_accel_num[i] = rot_accel_variable->number();
    }
  }
  else if (!_has_beta && !_has_gamma && !_has_velocities && !_has_accelerations &&
           !_has_rot_velocities && !_has_rot_accelerations)
  {
    _du_dot_du = &coupledDotDu("displacements", 0);
    _du_dotdot_du = &coupledDotDotDu("displacements", 0);
  }
  else
    mooseError("InertialForceBeam: Either all or none of `beta`, `gamma`, `velocities`, "
               "`accelerations`, `rotational_velocities` and `rotational_accelerations` should be "
               "provided as input.");
 
  // fetch coupled displacements and rotations
  for (unsigned int i = 0; i < _ndisp; ++i)
  {
    MooseVariable * disp_variable = getVar("displacements", i);
    _disp_num[i] = disp_variable->number();
 
    MooseVariable * rot_variable = getVar("rotations", i);
    _rot_num[i] = rot_variable->number();
  }
}

Member Function Documentation

computeJacobian()

void InertialForceBeam::computeJacobian ( )

overridevirtual

Definition at line 425 of file InertialForceBeam.C.

{
  prepareMatrixTag(_assembly, _var.number(), _var.number());
 
  mooseAssert(_beta > 0.0, "InertialForceBeam: Beta parameter should be positive.");
 
  Real factor = 0.0;
  if (_has_beta)
    factor = 1.0 / (_beta * _dt * _dt) + _eta[0] * (1.0 + _alpha) * _gamma / _beta / _dt;
  else
    factor = (*_du_dotdot_du)[0] + _eta[0] * (1.0 + _alpha) * (*_du_dot_du)[0];
 
  for (unsigned int i = 0; i < _test.size(); ++i)
  {
    for (unsigned int j = 0; j < _phi.size(); ++j)
    {
      if (_component < 3)
        _local_ke(i, j) = (i == j ? 1.0 / 3.0 : 1.0 / 6.0) * _density[0] * _area[0] *
                          _original_length[0] * factor;
      else if (_component > 2)
      {
        RankTwoTensor I;
        if (_has_Ix)
          I(0, 0) = _Ix[0];
        else
          I(0, 0) = _Iy[0] + _Iz[0];
        I(1, 1) = _Iz[0];
        I(2, 2) = _Iy[0];
 
        // conversion from local config to global coordinate system
        RankTwoTensor Ig = _original_local_config[0].transpose() * I * _original_local_config[0];
 
        _local_ke(i, j) = (i == j ? 1.0 / 3.0 : 1.0 / 6.0) * _density[0] *
                          Ig(_component - 3, _component - 3) * _original_length[0] * factor;
      }
    }
  }
 
  accumulateTaggedLocalMatrix();
 
  if (_has_diag_save_in)
  {
    unsigned int rows = _local_ke.m();
    DenseVector<Number> diag(rows);
    for (unsigned int i = 0; i < rows; ++i)
      diag(i) = _local_ke(i, i);
 
    Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
    for (unsigned int i = 0; i < _diag_save_in.size(); ++i)
      _diag_save_in[i]->sys().solution().add_vector(diag, _diag_save_in[i]->dofIndices());
  }
}

Referenced by computeOffDiagJacobian().

computeOffDiagJacobian()

void InertialForceBeam::computeOffDiagJacobian ( MooseVariableFEBase &  jvar )

overridevirtual

Definition at line 479 of file InertialForceBeam.C.

{
  mooseAssert(_beta > 0.0, "InertialForceBeam: Beta parameter should be positive.");
 
  Real factor = 0.0;
  if (_has_beta)
    factor = 1.0 / (_beta * _dt * _dt) + _eta[0] * (1.0 + _alpha) * _gamma / _beta / _dt;
  else
    factor = (*_du_dotdot_du)[0] + _eta[0] * (1.0 + _alpha) * (*_du_dot_du)[0];
 
  size_t jvar_num = jvar.number();
  if (jvar_num == _var.number())
    computeJacobian();
  else
  {
    unsigned int coupled_component = 0;
    bool disp_coupled = false;
    bool rot_coupled = false;
 
    for (unsigned int i = 0; i < _ndisp; ++i)
    {
      if (jvar_num == _disp_num[i] && _component > 2)
      {
        coupled_component = i;
        disp_coupled = true;
        break;
      }
    }
 
    for (unsigned int i = 0; i < _nrot; ++i)
    {
      if (jvar_num == _rot_num[i])
      {
        coupled_component = i + 3;
        rot_coupled = true;
        break;
      }
    }
 
    prepareMatrixTag(_assembly, _var.number(), jvar_num);
 
    if (disp_coupled || rot_coupled)
    {
      for (unsigned int i = 0; i < _test.size(); ++i)
      {
        for (unsigned int j = 0; j < _phi.size(); ++j)
        {
          if (_component < 3 && coupled_component > 2)
          {
            RankTwoTensor A;
            A(0, 1) = _Az[0];
            A(0, 2) = -_Ay[0];
            A(1, 0) = -_Az[0];
            A(2, 0) = _Ay[0];
 
            // conversion from local config to global coordinate system
            const RankTwoTensor Ag =
                _original_local_config[0].transpose() * A * _original_local_config[0];
 
            _local_ke(i, j) += (i == j ? 1.0 / 3.0 : 1.0 / 6.0) * _density[0] *
                               Ag(_component, coupled_component - 3) * _original_length[0] * factor;
          }
          else if (_component > 2 && coupled_component < 3)
          {
            RankTwoTensor A;
            A(0, 1) = -_Az[0];
            A(0, 2) = _Ay[0];
            A(1, 0) = _Az[0];
            A(2, 0) = -_Ay[0];
 
            // conversion from local config to global coordinate system
            const RankTwoTensor Ag =
                _original_local_config[0].transpose() * A * _original_local_config[0];
 
            _local_ke(i, j) += (i == j ? 1.0 / 3.0 : 1.0 / 6.0) * _density[0] *
                               Ag(_component - 3, coupled_component) * _original_length[0] * factor;
          }
          else if (_component > 2 && coupled_component > 2)
          {
            RankTwoTensor I;
            if (_has_Ix)
              I(0, 0) = _Ix[0];
            else
              I(0, 0) = _Iy[0] + _Iz[0];
            I(1, 1) = _Iz[0];
            I(2, 2) = _Iy[0];
 
            // conversion from local config to global coordinate system
            const RankTwoTensor Ig =
                _original_local_config[0].transpose() * I * _original_local_config[0];
 
            _local_ke(i, j) += (i == j ? 1.0 / 3.0 : 1.0 / 6.0) * _density[0] *
                               Ig(_component - 3, coupled_component - 3) * _original_length[0] *
                               factor;
          }
        }
      }
    }
 
    accumulateTaggedLocalMatrix();
  }
}

computeQpResidual()

virtual Real InertialForceBeam::computeQpResidual ( )

inlineoverrideprotectedvirtual

Definition at line 37 of file InertialForceBeam.h.

{ return 0.0; };

computeResidual()

void InertialForceBeam::computeResidual ( )

overridevirtual

Definition at line 165 of file InertialForceBeam.C.

{
  prepareVectorTag(_assembly, _var.number());
 
  if (_dt != 0.0)
  {
    // fetch the two end nodes for _current_elem
    std::vector<const Node *> node;
    for (unsigned int i = 0; i < 2; ++i)
      node.push_back(_current_elem->node_ptr(i));
 
    // Fetch the solution for the two end nodes at time t
    NonlinearSystemBase & nonlinear_sys = _fe_problem.getNonlinearSystemBase();
 
    if (_has_beta)
    {
      const NumericVector<Number> & sol = *nonlinear_sys.currentSolution();
      const NumericVector<Number> & sol_old = nonlinear_sys.solutionOld();
 
      AuxiliarySystem & aux = _fe_problem.getAuxiliarySystem();
      const NumericVector<Number> & aux_sol_old = aux.solutionOld();
 
      mooseAssert(_beta > 0.0, "InertialForceBeam: Beta parameter should be positive.");
      mooseAssert(_eta[0] >= 0.0, "InertialForceBeam: Eta parameter should be non-negative.");
 
      for (unsigned int i = 0; i < _ndisp; ++i)
      {
        // obtain delta displacement
        unsigned int dof_index_0 = node[0]->dof_number(nonlinear_sys.number(), _disp_num[i], 0);
        unsigned int dof_index_1 = node[1]->dof_number(nonlinear_sys.number(), _disp_num[i], 0);
        const Real disp_0 = sol(dof_index_0) - sol_old(dof_index_0);
        const Real disp_1 = sol(dof_index_1) - sol_old(dof_index_1);
 
        dof_index_0 = node[0]->dof_number(nonlinear_sys.number(), _rot_num[i], 0);
        dof_index_1 = node[1]->dof_number(nonlinear_sys.number(), _rot_num[i], 0);
        const Real rot_0 = sol(dof_index_0) - sol_old(dof_index_0);
        const Real rot_1 = sol(dof_index_1) - sol_old(dof_index_1);
 
        // obtain new translational and rotational velocities and accelerations using newmark-beta
        // time integration
        _vel_old_0(i) = aux_sol_old(node[0]->dof_number(aux.number(), _vel_num[i], 0));
        _vel_old_1(i) = aux_sol_old(node[1]->dof_number(aux.number(), _vel_num[i], 0));
        const Real accel_old_0 = aux_sol_old(node[0]->dof_number(aux.number(), _accel_num[i], 0));
        const Real accel_old_1 = aux_sol_old(node[1]->dof_number(aux.number(), _accel_num[i], 0));
 
        _rot_vel_old_0(i) = aux_sol_old(node[0]->dof_number(aux.number(), _rot_vel_num[i], 0));
        _rot_vel_old_1(i) = aux_sol_old(node[1]->dof_number(aux.number(), _rot_vel_num[i], 0));
        const Real rot_accel_old_0 =
            aux_sol_old(node[0]->dof_number(aux.number(), _rot_accel_num[i], 0));
        const Real rot_accel_old_1 =
            aux_sol_old(node[1]->dof_number(aux.number(), _rot_accel_num[i], 0));
 
        _accel_0(i) =
            1. / _beta * (disp_0 / (_dt * _dt) - _vel_old_0(i) / _dt - accel_old_0 * (0.5 - _beta));
        _accel_1(i) =
            1. / _beta * (disp_1 / (_dt * _dt) - _vel_old_1(i) / _dt - accel_old_1 * (0.5 - _beta));
        _rot_accel_0(i) =
            1. / _beta *
            (rot_0 / (_dt * _dt) - _rot_vel_old_0(i) / _dt - rot_accel_old_0 * (0.5 - _beta));
        _rot_accel_1(i) =
            1. / _beta *
            (rot_1 / (_dt * _dt) - _rot_vel_old_1(i) / _dt - rot_accel_old_1 * (0.5 - _beta));
 
        _vel_0(i) = _vel_old_0(i) + (_dt * (1 - _gamma)) * accel_old_0 + _gamma * _dt * _accel_0(i);
        _vel_1(i) = _vel_old_1(i) + (_dt * (1 - _gamma)) * accel_old_1 + _gamma * _dt * _accel_1(i);
        _rot_vel_0(i) = _rot_vel_old_0(i) + (_dt * (1 - _gamma)) * rot_accel_old_0 +
                        _gamma * _dt * _rot_accel_0(i);
        _rot_vel_1(i) = _rot_vel_old_1(i) + (_dt * (1 - _gamma)) * rot_accel_old_1 +
                        _gamma * _dt * _rot_accel_1(i);
      }
    }
    else
    {
      if (!nonlinear_sys.solutionUDot())
        mooseError("InertialForceBeam: Time derivative of solution (`u_dot`) is not stored. Please "
                   "set uDotRequested() to true in FEProblemBase before requesting `u_dot`.");
 
      if (!nonlinear_sys.solutionUDotOld())
        mooseError("InertialForceBeam: Old time derivative of solution (`u_dot_old`) is not "
                   "stored. Please set uDotOldRequested() to true in FEProblemBase before "
                   "requesting `u_dot_old`.");
 
      if (!nonlinear_sys.solutionUDotDot())
        mooseError("InertialForceBeam: Second time derivative of solution (`u_dotdot`) is not "
                   "stored. Please set uDotDotRequested() to true in FEProblemBase before "
                   "requesting `u_dotdot`.");
 
      const NumericVector<Number> & vel = *nonlinear_sys.solutionUDot();
      const NumericVector<Number> & vel_old = *nonlinear_sys.solutionUDotOld();
      const NumericVector<Number> & accel = *nonlinear_sys.solutionUDotDot();
 
      for (unsigned int i = 0; i < _ndisp; ++i)
      {
        // translational velocities and accelerations
        unsigned int dof_index_0 = node[0]->dof_number(nonlinear_sys.number(), _disp_num[i], 0);
        unsigned int dof_index_1 = node[1]->dof_number(nonlinear_sys.number(), _disp_num[i], 0);
        _vel_0(i) = vel(dof_index_0);
        _vel_1(i) = vel(dof_index_1);
        _vel_old_0(i) = vel_old(dof_index_0);
        _vel_old_1(i) = vel_old(dof_index_1);
        _accel_0(i) = accel(dof_index_0);
        _accel_1(i) = accel(dof_index_1);
 
        // rotational velocities and accelerations
        dof_index_0 = node[0]->dof_number(nonlinear_sys.number(), _rot_num[i], 0);
        dof_index_1 = node[1]->dof_number(nonlinear_sys.number(), _rot_num[i], 0);
        _rot_vel_0(i) = vel(dof_index_0);
        _rot_vel_1(i) = vel(dof_index_1);
        _rot_vel_old_0(i) = vel_old(dof_index_0);
        _rot_vel_old_1(i) = vel_old(dof_index_1);
        _rot_accel_0(i) = accel(dof_index_0);
        _rot_accel_1(i) = accel(dof_index_1);
      }
    }
 
    // transform translational and rotational velocities and accelerations to the initial local
    // configuration of the beam
    _local_vel_old_0 = _original_local_config[0] * _vel_old_0;
    _local_vel_old_1 = _original_local_config[0] * _vel_old_1;
    _local_vel_0 = _original_local_config[0] * _vel_0;
    _local_vel_1 = _original_local_config[0] * _vel_1;
    _local_accel_0 = _original_local_config[0] * _accel_0;
    _local_accel_1 = _original_local_config[0] * _accel_1;
 
    _local_rot_vel_old_0 = _original_local_config[0] * _rot_vel_old_0;
    _local_rot_vel_old_1 = _original_local_config[0] * _rot_vel_old_1;
    _local_rot_vel_0 = _original_local_config[0] * _rot_vel_0;
    _local_rot_vel_1 = _original_local_config[0] * _rot_vel_1;
    _local_rot_accel_0 = _original_local_config[0] * _rot_accel_0;
    _local_rot_accel_1 = _original_local_config[0] * _rot_accel_1;
 
    // local residual
    for (unsigned int i = 0; i < _ndisp; ++i)
    {
      if (_component < 3)
      {
        _local_force[0](i) = _density[0] * _area[0] * _original_length[0] / 3.0 *
                             (_local_accel_0(i) + _local_accel_1(i) / 2.0 +
                              _eta[0] * (1.0 + _alpha) * (_local_vel_0(i) + _local_vel_1(i) / 2.0) -
                              _alpha * _eta[0] * (_local_vel_old_0(i) + _local_vel_old_1(i) / 2.0));
        _local_force[1](i) = _density[0] * _area[0] * _original_length[0] / 3.0 *
                             (_local_accel_1(i) + _local_accel_0(i) / 2.0 +
                              _eta[0] * (1.0 + _alpha) * (_local_vel_1(i) + _local_vel_0(i) / 2.0) -
                              _alpha * _eta[0] * (_local_vel_old_1(i) + _local_vel_old_0(i) / 2.0));
      }
 
      if (_component > 2)
      {
        Real I = _Iy[0] + _Iz[0];
        if (_has_Ix && (i == 0))
          I = _Ix[0];
        if (i == 1)
          I = _Iz[0];
        else if (i == 2)
          I = _Iy[0];
 
        _local_moment[0](i) =
            _density[0] * I * _original_length[0] / 3.0 *
            (_local_rot_accel_0(i) + _local_rot_accel_1(i) / 2.0 +
             _eta[0] * (1.0 + _alpha) * (_local_rot_vel_0(i) + _local_rot_vel_1(i) / 2.0) -
             _alpha * _eta[0] * (_local_rot_vel_old_0(i) + _local_rot_vel_old_1(i) / 2.0));
        _local_moment[1](i) =
            _density[0] * I * _original_length[0] / 3.0 *
            (_local_rot_accel_1(i) + _local_rot_accel_0(i) / 2.0 +
             _eta[0] * (1.0 + _alpha) * (_local_rot_vel_1(i) + _local_rot_vel_0(i) / 2.0) -
             _alpha * _eta[0] * (_local_rot_vel_old_1(i) + _local_rot_vel_old_0(i) / 2.0));
      }
    }
 
    // If Ay or Az are non-zero, contribution of rotational accelerations to translational forces
    // and vice versa have to be added
    if (_component < 3)
    {
      _local_force[0](0) +=
          _density[0] * _original_length[0] / 3.0 *
          (_Az[0] * (_local_rot_accel_0(1) + _local_rot_accel_1(1) / 2.0 +
                     _eta[0] * (1.0 + _alpha) * (_local_rot_vel_0(1) + _local_rot_vel_1(1) / 2.0) -
                     _alpha * _eta[0] * (_local_rot_vel_old_0(1) + _local_rot_vel_old_1(1) / 2.0)) -
           _Ay[0] * (_local_rot_accel_0(2) + _local_rot_accel_1(2) / 2.0 +
                     _eta[0] * (1.0 + _alpha) * (_local_rot_vel_0(2) + _local_rot_vel_1(2) / 2.0) -
                     _alpha * _eta[0] * (_local_rot_vel_old_0(2) + _local_rot_vel_old_1(2) / 2.0)));
      _local_force[1](0) +=
          _density[0] * _original_length[0] / 3.0 *
          (_Az[0] * (_local_rot_accel_1(1) + _local_rot_accel_0(1) / 2.0 +
                     _eta[0] * (1.0 + _alpha) * (_local_rot_vel_1(1) + _local_rot_vel_0(1) / 2.0) -
                     _alpha * _eta[0] * (_local_rot_vel_old_1(1) + _local_rot_vel_old_0(1) / 2.0)) -
           _Ay[0] * (_local_rot_accel_1(2) + _local_rot_accel_0(2) / 2.0 +
                     _eta[0] * (1.0 + _alpha) * (_local_rot_vel_1(2) + _local_rot_vel_0(2) / 2.0) -
                     _alpha * _eta[0] * (_local_rot_vel_old_1(2) + _local_rot_vel_old_0(2) / 2.0)));
 
      _local_force[0](1) +=
          -_density[0] * _original_length[0] / 3.0 * _Az[0] *
          (_local_rot_accel_0(0) + _local_rot_accel_1(0) / 2.0 +
           _eta[0] * (1.0 + _alpha) * (_local_rot_vel_0(0) + _local_rot_vel_1(0) / 2.0) -
           _alpha * _eta[0] * (_local_rot_vel_old_0(0) + _local_rot_vel_old_1(0) / 2.0));
      _local_force[1](1) +=
          -_density[0] * _original_length[0] / 3.0 * _Az[0] *
          (_local_rot_accel_1(0) + _local_rot_accel_0(0) / 2.0 +
           _eta[0] * (1.0 + _alpha) * (_local_rot_vel_1(0) + _local_rot_vel_0(0) / 2.0) -
           _alpha * _eta[0] * (_local_rot_vel_old_1(0) + _local_rot_vel_old_0(0) / 2.0));
 
      _local_force[0](2) +=
          _density[0] * _original_length[0] / 3.0 * _Ay[0] *
          (_local_rot_accel_0(0) + _local_rot_accel_1(0) / 2.0 +
           _eta[0] * (1.0 + _alpha) * (_local_rot_vel_0(0) + _local_rot_vel_1(0) / 2.0) -
           _alpha * _eta[0] * (_local_rot_vel_old_0(0) + _local_rot_vel_old_1(0) / 2.0));
      _local_force[1](2) +=
          _density[0] * _original_length[0] / 3.0 * _Ay[0] *
          (_local_rot_accel_1(0) + _local_rot_accel_0(0) / 2.0 +
           _eta[0] * (1.0 + _alpha) * (_local_rot_vel_1(0) + _local_rot_vel_0(0) / 2.0) -
           _alpha * _eta[0] * (_local_rot_vel_old_1(0) + _local_rot_vel_old_0(0) / 2.0));
    }
    else
    {
      _local_moment[0](0) += _density[0] * _original_length[0] / 3.0 *
                             (-_Az[0] * (_local_accel_0(1) + _local_accel_1(1) / 2.0) +
                              _Ay[0] * (_local_accel_0(1) + _local_accel_1(1) / 2.0));
      _local_moment[1](0) += _density[0] * _original_length[0] / 3.0 *
                             (-_Az[0] * (_local_accel_1(1) + _local_accel_0(1) / 2.0) +
                              _Ay[0] * (_local_accel_1(1) + _local_accel_0(1) / 2.0));
 
      _local_moment[0](1) += _density[0] * _original_length[0] / 3.0 * _Az[0] *
                             (_local_accel_0(0) + _local_accel_1(0) / 2.0);
      _local_moment[1](1) += _density[0] * _original_length[0] / 3.0 * _Az[0] *
                             (_local_accel_1(0) + _local_accel_0(0) / 2.0);
 
      _local_moment[0](2) += -_density[0] * _original_length[0] / 3.0 * _Ay[0] *
                             (_local_accel_0(0) + _local_accel_1(0) / 2.0);
      _local_moment[1](2) += -_density[0] * _original_length[0] / 3.0 * _Ay[0] *
                             (_local_accel_1(0) + _local_accel_0(0) / 2.0);
    }
 
    // Global force and moments
    if (_component < 3)
    {
      _global_force_0 = _original_local_config[0] * _local_force[0];
      _global_force_1 = _original_local_config[0] * _local_force[1];
      _local_re(0) = _global_force_0(_component);
      _local_re(1) = _global_force_1(_component);
    }
    else
    {
      _global_moment_0 = _original_local_config[0] * _local_moment[0];
      _global_moment_1 = _original_local_config[0] * _local_moment[1];
      _local_re(0) = _global_moment_0(_component - 3);
      _local_re(1) = _global_moment_1(_component - 3);
    }
  }
 
  accumulateTaggedLocalResidual();
 
  if (_has_save_in)
  {
    Threads::spin_mutex::scoped_lock lock(Threads::spin_mtx);
    for (unsigned int i = 0; i < _save_in.size(); ++i)
      _save_in[i]->sys().solution().add_vector(_local_re, _save_in[i]->dofIndices());
  }
}

validParams()

InputParameters InertialForceBeam::validParams ( )

static

Definition at line 24 of file InertialForceBeam.C.

{
  InputParameters params = TimeKernel::validParams();
  params.addClassDescription("Calculates the residual for the inertial force/moment and the "
                             "contribution of mass dependent Rayleigh damping and HHT time "
                             "integration scheme.");
  params.set<bool>("use_displaced_mesh") = true;
  params.addRequiredCoupledVar(
      "rotations",
      "The rotational variables appropriate for the simulation geometry and coordinate system");
  params.addRequiredCoupledVar(
      "displacements",
      "The displacement variables appropriate for the simulation geometry and coordinate system");
  params.addCoupledVar("velocities", "Translational velocity variables");
  params.addCoupledVar("accelerations", "Translational acceleration variables");
  params.addCoupledVar("rotational_velocities", "Rotational velocity variables");
  params.addCoupledVar("rotational_accelerations", "Rotational acceleration variables");
  params.addRangeCheckedParam<Real>(
      "beta", "beta>0.0", "beta parameter for Newmark Time integration");
  params.addRangeCheckedParam<Real>(
      "gamma", "gamma>0.0", "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.addRangeCheckedParam<Real>("alpha",
                                    0.0,
                                    "alpha >= -0.3333 & alpha <= 0.0",
                                    "Alpha parameter for mass dependent numerical damping induced "
                                    "by HHT time integration scheme");
  params.addParam<MaterialPropertyName>(
      "density",
      "density",
      "Name of Material Property  or a constant real number defining the density of the beam.");
  params.addRequiredCoupledVar("area", "Variable containing cross-section area");
  params.addCoupledVar("Ay", 0.0, "Variable containing first moment of area about y axis");
  params.addCoupledVar("Az", 0.0, "Variable containing first moment of area about z axis");
  params.addCoupledVar("Ix",
                       "Variable containing second moment of area about x axis. Defaults to Iy+Iz");
  params.addRequiredCoupledVar("Iy", "Variable containing second moment of area about y axis");
  params.addRequiredCoupledVar("Iz", "Variable containing second moment of area about z axis");
  params.addRequiredRangeCheckedParam<unsigned int>(
      "component",
      "component<6",
      "An integer corresponding to the direction "
      "the variable this kernel acts in. (0 for disp_x, "
      "1 for disp_y, 2 for disp_z, 3 for rot_x, 4 for rot_y and 5 for rot_z)");
  return params;
}

Member Data Documentation

_accel_0

RealVectorValue InertialForceBeam::_accel_0

private

Current translational and rotational accelerations at the two nodes of the beam in the global coordinate system.

Definition at line 149 of file InertialForceBeam.h.

Referenced by computeResidual().

_accel_1

RealVectorValue InertialForceBeam::_accel_1

private

Definition at line 149 of file InertialForceBeam.h.

Referenced by computeResidual().

_accel_num

std::vector<unsigned int> InertialForceBeam::_accel_num

private

Variable numbers corresponding to acceleraion aux variables.

Definition at line 68 of file InertialForceBeam.h.

Referenced by computeResidual(), and InertialForceBeam().

_alpha

const Real InertialForceBeam::_alpha

private

HHT time integration parameter.

Definition at line 119 of file InertialForceBeam.h.

Referenced by computeJacobian(), computeOffDiagJacobian(), and computeResidual().

_area

const VariableValue& InertialForceBeam::_area

private

Coupled variable for beam cross-sectional area.

Definition at line 77 of file InertialForceBeam.h.

Referenced by computeJacobian(), and computeResidual().

_Ay

const VariableValue& InertialForceBeam::_Ay

private

Coupled variable for first moment of area of beam in y direction, i.e., integral of y*dA over the cross-section.

Definition at line 83 of file InertialForceBeam.h.

Referenced by computeOffDiagJacobian(), and computeResidual().

_Az

const VariableValue& InertialForceBeam::_Az

private

Coupled variable for first moment of area of beam in z direction, i.e., integral of z*dA over the cross-section.

Definition at line 89 of file InertialForceBeam.h.

Referenced by computeOffDiagJacobian(), and computeResidual().

_beta

const Real InertialForceBeam::_beta

private

Newmark time integration parameter.

Definition at line 110 of file InertialForceBeam.h.

Referenced by computeJacobian(), computeOffDiagJacobian(), and computeResidual().

_component

const unsigned int InertialForceBeam::_component

private

Direction along which residual is calculated.

Definition at line 131 of file InertialForceBeam.h.

Referenced by computeJacobian(), computeOffDiagJacobian(), and computeResidual().

_density

const MaterialProperty<Real>& InertialForceBeam::_density

private

Density of the beam.

Definition at line 50 of file InertialForceBeam.h.

Referenced by computeJacobian(), computeOffDiagJacobian(), and computeResidual().

_disp_num

std::vector<unsigned int> InertialForceBeam::_disp_num

private

Variable numbers corresponding to displacement variables.

Definition at line 62 of file InertialForceBeam.h.

Referenced by computeOffDiagJacobian(), computeResidual(), and InertialForceBeam().

_du_dot_du

const VariableValue* InertialForceBeam::_du_dot_du

private

Coupled variable for du_dot_du calculated by time integrator.

Definition at line 184 of file InertialForceBeam.h.

Referenced by computeJacobian(), computeOffDiagJacobian(), and InertialForceBeam().

_du_dotdot_du

const VariableValue* InertialForceBeam::_du_dotdot_du

private

Coupled variable for du_dotdot_du calculated by time integrator.

Definition at line 189 of file InertialForceBeam.h.

Referenced by InertialForceBeam().

_eta

const MaterialProperty<Real>& InertialForceBeam::_eta

private

Mass proportional Rayleigh damping parameter.

Definition at line 116 of file InertialForceBeam.h.

Referenced by computeJacobian(), computeOffDiagJacobian(), and computeResidual().

_gamma

const Real InertialForceBeam::_gamma

private

Newmark time integraion parameter.

Definition at line 113 of file InertialForceBeam.h.

Referenced by computeJacobian(), computeOffDiagJacobian(), and computeResidual().

_global_force_0

RealVectorValue InertialForceBeam::_global_force_0

private

Forces and moments at the two end nodes of the beam in the global coordinate system.

Definition at line 179 of file InertialForceBeam.h.

Referenced by computeResidual().

_global_force_1

RealVectorValue InertialForceBeam::_global_force_1

private

Definition at line 179 of file InertialForceBeam.h.

Referenced by computeResidual().

_global_moment_0

RealVectorValue InertialForceBeam::_global_moment_0

private

Definition at line 179 of file InertialForceBeam.h.

Referenced by computeResidual().

_global_moment_1

RealVectorValue InertialForceBeam::_global_moment_1

private

Definition at line 179 of file InertialForceBeam.h.

Referenced by computeResidual().

_has_accelerations

const bool InertialForceBeam::_has_accelerations

private

Definition at line 45 of file InertialForceBeam.h.

Referenced by InertialForceBeam().

_has_beta

const bool InertialForceBeam::_has_beta

private

Booleans for validity of params.

Definition at line 37 of file InertialForceBeam.h.

Referenced by computeJacobian(), computeOffDiagJacobian(), computeResidual(), and InertialForceBeam().

_has_gamma

const bool InertialForceBeam::_has_gamma

private

Definition at line 42 of file InertialForceBeam.h.

Referenced by InertialForceBeam().

_has_Ix

const bool InertialForceBeam::_has_Ix

private

Definition at line 47 of file InertialForceBeam.h.

Referenced by computeJacobian(), computeOffDiagJacobian(), and computeResidual().

_has_rot_accelerations

const bool InertialForceBeam::_has_rot_accelerations

private

Definition at line 46 of file InertialForceBeam.h.

Referenced by InertialForceBeam().

_has_rot_velocities

const bool InertialForceBeam::_has_rot_velocities

private

Definition at line 44 of file InertialForceBeam.h.

Referenced by InertialForceBeam().

_has_velocities

const bool InertialForceBeam::_has_velocities

private

Definition at line 43 of file InertialForceBeam.h.

Referenced by InertialForceBeam().

_Ix

const VariableValue& InertialForceBeam::_Ix

private

Coupled variable for second moment of area of beam in x direction, i.e., integral of (y^2+z^2)*dA over the cross-section.

Definition at line 95 of file InertialForceBeam.h.

Referenced by computeJacobian(), computeOffDiagJacobian(), and computeResidual().

_Iy

const VariableValue& InertialForceBeam::_Iy

private

Coupled variable for second moment of area of beam in y direction, i.e., integral of y^2*dA over the cross-section.

Definition at line 101 of file InertialForceBeam.h.

Referenced by computeJacobian(), computeOffDiagJacobian(), and computeResidual().

_Iz

const VariableValue& InertialForceBeam::_Iz

private

Coupled variable for second momemnt of area of beam in z direction, i.e., integral of z^2*dA over the cross-section.

Definition at line 107 of file InertialForceBeam.h.

Referenced by computeJacobian(), computeOffDiagJacobian(), and computeResidual().

_local_accel_0

RealVectorValue InertialForceBeam::_local_accel_0

private

Current translational and rotational accelerations at the two nodes of the beam in the initial beam local coordinate system.

Definition at line 167 of file InertialForceBeam.h.

Referenced by computeResidual().

_local_accel_1

RealVectorValue InertialForceBeam::_local_accel_1

private

Definition at line 167 of file InertialForceBeam.h.

Referenced by computeResidual().

_local_force

std::vector<RealVectorValue> InertialForceBeam::_local_force

private

Forces and moments at the two end nodes of the beam in the initial beam local configuration.

Definition at line 173 of file InertialForceBeam.h.

Referenced by computeResidual().

_local_moment

std::vector<RealVectorValue> InertialForceBeam::_local_moment

private

Definition at line 173 of file InertialForceBeam.h.

Referenced by computeResidual().

_local_rot_accel_0

RealVectorValue InertialForceBeam::_local_rot_accel_0

private

Definition at line 167 of file InertialForceBeam.h.

Referenced by computeResidual().

_local_rot_accel_1

RealVectorValue InertialForceBeam::_local_rot_accel_1

private

Definition at line 167 of file InertialForceBeam.h.

Referenced by computeResidual().

_local_rot_vel_0

RealVectorValue InertialForceBeam::_local_rot_vel_0

private

Definition at line 161 of file InertialForceBeam.h.

Referenced by computeResidual().

_local_rot_vel_1

RealVectorValue InertialForceBeam::_local_rot_vel_1

private

Definition at line 161 of file InertialForceBeam.h.

Referenced by computeResidual().

_local_rot_vel_old_0

RealVectorValue InertialForceBeam::_local_rot_vel_old_0

private

Definition at line 155 of file InertialForceBeam.h.

Referenced by computeResidual().

_local_rot_vel_old_1

RealVectorValue InertialForceBeam::_local_rot_vel_old_1

private

Definition at line 155 of file InertialForceBeam.h.

Referenced by computeResidual().

_local_vel_0

RealVectorValue InertialForceBeam::_local_vel_0

private

Current translational and rotational velocities at the two nodes of the beam in the initial beam local coordinate system.

Definition at line 161 of file InertialForceBeam.h.

Referenced by computeResidual().

_local_vel_1

RealVectorValue InertialForceBeam::_local_vel_1

private

Definition at line 161 of file InertialForceBeam.h.

Referenced by computeResidual().

_local_vel_old_0

RealVectorValue InertialForceBeam::_local_vel_old_0

private

Old translational and rotational velocities at the two nodes of the beam in the initial beam local coordinate system.

Definition at line 155 of file InertialForceBeam.h.

Referenced by computeResidual().

_local_vel_old_1

RealVectorValue InertialForceBeam::_local_vel_old_1

private

Definition at line 155 of file InertialForceBeam.h.

Referenced by computeResidual().

_ndisp

unsigned int InertialForceBeam::_ndisp

private

Number of coupled displacement variables.

Definition at line 56 of file InertialForceBeam.h.

Referenced by computeOffDiagJacobian(), computeResidual(), and InertialForceBeam().

_nrot

unsigned int InertialForceBeam::_nrot

private

Number of coupled rotational variables.

Definition at line 53 of file InertialForceBeam.h.

Referenced by computeOffDiagJacobian(), and InertialForceBeam().

_original_length

const MaterialProperty<Real>& InertialForceBeam::_original_length

private

Initial length of beam.

Definition at line 128 of file InertialForceBeam.h.

Referenced by computeJacobian(), computeOffDiagJacobian(), and computeResidual().

_original_local_config

const MaterialProperty<RankTwoTensor>& InertialForceBeam::_original_local_config

private

Rotational transformation from global to initial beam local coordinate system.

Definition at line 125 of file InertialForceBeam.h.

Referenced by computeJacobian(), computeOffDiagJacobian(), and computeResidual().

_rot_accel_0

RealVectorValue InertialForceBeam::_rot_accel_0

private

Definition at line 149 of file InertialForceBeam.h.

Referenced by computeResidual().

_rot_accel_1

RealVectorValue InertialForceBeam::_rot_accel_1

private

Definition at line 149 of file InertialForceBeam.h.

Referenced by computeResidual().

_rot_accel_num

std::vector<unsigned int> InertialForceBeam::_rot_accel_num

private

Variable numbers corresponding to rotational acceleration aux variables.

Definition at line 74 of file InertialForceBeam.h.

Referenced by computeResidual(), and InertialForceBeam().

_rot_num

std::vector<unsigned int> InertialForceBeam::_rot_num

private

Variable numbers corresponding to rotational variables.

Definition at line 59 of file InertialForceBeam.h.

Referenced by computeOffDiagJacobian(), computeResidual(), and InertialForceBeam().

_rot_vel_0

RealVectorValue InertialForceBeam::_rot_vel_0

private

Definition at line 143 of file InertialForceBeam.h.

Referenced by computeResidual().

_rot_vel_1

RealVectorValue InertialForceBeam::_rot_vel_1

private

Definition at line 143 of file InertialForceBeam.h.

Referenced by computeResidual().

_rot_vel_num

std::vector<unsigned int> InertialForceBeam::_rot_vel_num

private

Variable numbers corresponding to rotational velocity aux variables.

Definition at line 71 of file InertialForceBeam.h.

Referenced by computeResidual(), and InertialForceBeam().

_rot_vel_old_0

RealVectorValue InertialForceBeam::_rot_vel_old_0

private

Definition at line 137 of file InertialForceBeam.h.

Referenced by computeResidual().

_rot_vel_old_1

RealVectorValue InertialForceBeam::_rot_vel_old_1

private

Definition at line 137 of file InertialForceBeam.h.

Referenced by computeResidual().

_vel_0

RealVectorValue InertialForceBeam::_vel_0

private

Current translational and rotational velocities at the two nodes of the beam in the global coordinate system.

Definition at line 143 of file InertialForceBeam.h.

Referenced by computeResidual().

_vel_1

RealVectorValue InertialForceBeam::_vel_1

private

Definition at line 143 of file InertialForceBeam.h.

Referenced by computeResidual().

_vel_num

std::vector<unsigned int> InertialForceBeam::_vel_num

private

Variable numbers corresponding to velocity aux variables.

Definition at line 65 of file InertialForceBeam.h.

Referenced by computeResidual(), and InertialForceBeam().

_vel_old_0

RealVectorValue InertialForceBeam::_vel_old_0

private

Old translational and rotational velocities at the two nodes of the beam in the global coordinate system.

Definition at line 137 of file InertialForceBeam.h.

Referenced by computeResidual().

_vel_old_1

RealVectorValue InertialForceBeam::_vel_old_1

private

Definition at line 137 of file InertialForceBeam.h.

Referenced by computeResidual().

---

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

• tensor_mechanics/include/kernels/InertialForceBeam.h
• tensor_mechanics/src/kernels/InertialForceBeam.C