ConstantGrainForceAndTorque.C

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
 
#include "ConstantGrainForceAndTorque.h"
 
registerMooseObject("PhaseFieldApp", ConstantGrainForceAndTorque);
 
template <>
InputParameters
validParams<ConstantGrainForceAndTorque>()
{
  InputParameters params = validParams<GeneralUserObject>();
  params.addClassDescription("Userobject for calculating force and torque acting on a grain");
  params.addParam<std::vector<Real>>("force", "force acting on grains");
  params.addParam<std::vector<Real>>("torque", "torque acting on grains");
  return params;
}
 
ConstantGrainForceAndTorque::ConstantGrainForceAndTorque(const InputParameters & parameters)
  : GrainForceAndTorqueInterface(),
    GeneralUserObject(parameters),
    _F(getParam<std::vector<Real>>("force")),
    _M(getParam<std::vector<Real>>("torque")),
    _grain_num(_F.size() / 3),
    _ncomp(6 * _grain_num),
    _force_values(_grain_num),
    _torque_values(_grain_num)
{
}
 
void
ConstantGrainForceAndTorque::initialize()
{
  unsigned int total_dofs = _subproblem.es().n_dofs();
  for (unsigned int i = 0; i < _grain_num; ++i)
  {
    _force_values[i](0) = _F[3 * i + 0];
    _force_values[i](1) = _F[3 * i + 1];
    _force_values[i](2) = _F[3 * i + 2];
    _torque_values[i](0) = _M[3 * i + 0];
    _torque_values[i](1) = _M[3 * i + 1];
    _torque_values[i](2) = _M[3 * i + 2];
  }
 
  if (_fe_problem.currentlyComputingJacobian())
  {
    _c_jacobians.assign(6 * _grain_num * total_dofs, 0.0);
    _eta_jacobians.resize(_grain_num);
    for (unsigned int i = 0; i < _grain_num; ++i)
      _eta_jacobians[i].assign(6 * _grain_num * total_dofs, 0.0);
  }
}
 
const std::vector<RealGradient> &
ConstantGrainForceAndTorque::getForceValues() const
{
  return _force_values;
}
 
const std::vector<RealGradient> &
ConstantGrainForceAndTorque::getTorqueValues() const
{
  return _torque_values;
}
 
const std::vector<Real> &
ConstantGrainForceAndTorque::getForceCJacobians() const
{
  return _c_jacobians;
}
 
const std::vector<std::vector<Real>> &
ConstantGrainForceAndTorque::getForceEtaJacobians() const
{
  return _eta_jacobians;
}