This Material calculates the force density acting on a particle/grain due to interaction between particles. More...

#include <ExternalForceDensityMaterial.h>

Inheritance diagram for ExternalForceDensityMaterial:

Public Member Functions
	ExternalForceDensityMaterial (const InputParameters &parameters)

Protected Member Functions
virtual void	computeQpProperties ()

Private Attributes
const Function &	_force_x

const Function &	_force_y

const Function &	_force_z

const VariableValue &	_c
	concentration field considered to be the density of particles More...

VariableName	_c_name

const Real	_k
	stiffness constant More...

unsigned int	_op_num

std::vector< const VariableValue * >	_vals

std::vector< VariableName >	_vals_name

MaterialProperty< std::vector< RealGradient > > &	_dF
	force density material More...

MaterialProperty< std::vector< RealGradient > > &	_dFdc
	first order derivative of force density material w.r.t c More...

std::vector< MaterialProperty< std::vector< RealGradient > > * >	_dFdeta

Detailed Description

This Material calculates the force density acting on a particle/grain due to interaction between particles.

Definition at line 26 of file ExternalForceDensityMaterial.h.

Constructor & Destructor Documentation

◆ ExternalForceDensityMaterial()

ExternalForceDensityMaterial::ExternalForceDensityMaterial ( const InputParameters & parameters )

Definition at line 31 of file ExternalForceDensityMaterial.C.

   : DerivativeMaterialInterface<Material>(parameters),
     _force_x(getFunction("force_x")),
     _force_y(getFunction("force_y")),
     _force_z(getFunction("force_z")),
     _c(coupledValue("c")),
     _c_name(getVar("c", 0)->name()),
     _k(getParam<Real>("k")),
     _op_num(coupledComponents(
         "etas")),   // determine number of grains from the number of names passed in.
     _vals(_op_num), // Size variable arrays
     _vals_name(_op_num),
     _dF(declareProperty<std::vector<RealGradient>>("force_density_ext")),
     _dFdc(declarePropertyDerivative<std::vector<RealGradient>>("force_density_ext", _c_name)),
     _dFdeta(_op_num)
 {
   // Loop through grains and load coupled variables into the arrays
   for (unsigned int i = 0; i < _op_num; ++i)
   {
     _vals[i] = &coupledValue("etas", i);
     _vals_name[i] = getVar("etas", i)->name();
     _dFdeta[i] =
         &declarePropertyDerivative<std::vector<RealGradient>>("force_density_ext", _vals_name[i]);
   }
 }

Member Function Documentation

◆ computeQpProperties()

void ExternalForceDensityMaterial::computeQpProperties ( )

protectedvirtual

Definition at line 58 of file ExternalForceDensityMaterial.C.

 {
   _dF[_qp].resize(_op_num);
   _dFdc[_qp].resize(_op_num);
  
   for (unsigned int i = 0; i < _op_num; ++i)
   {
     _dF[_qp][i](0) = _k * _c[_qp] * _force_x.value(_t, _q_point[_qp]) * (*_vals[i])[_qp];
     _dF[_qp][i](1) = _k * _c[_qp] * _force_y.value(_t, _q_point[_qp]) * (*_vals[i])[_qp];
     _dF[_qp][i](2) = _k * _c[_qp] * _force_z.value(_t, _q_point[_qp]) * (*_vals[i])[_qp];
  
     _dFdc[_qp][i](0) = _k * _force_x.value(_t, _q_point[_qp]) * (*_vals[i])[_qp];
     _dFdc[_qp][i](1) = _k * _force_y.value(_t, _q_point[_qp]) * (*_vals[i])[_qp];
     _dFdc[_qp][i](2) = _k * _force_z.value(_t, _q_point[_qp]) * (*_vals[i])[_qp];
   }
  
   for (unsigned int i = 0; i < _op_num; ++i)
   {
     (*_dFdeta[i])[_qp].resize(_op_num);
     for (unsigned int j = 0; j < _op_num; ++j)
     {
       (*_dFdeta[i])[_qp][j](0) = _k * _c[_qp] * _force_x.value(_t, _q_point[_qp]);
       (*_dFdeta[i])[_qp][j](1) = _k * _c[_qp] * _force_y.value(_t, _q_point[_qp]);
       (*_dFdeta[i])[_qp][j](2) = _k * _c[_qp] * _force_z.value(_t, _q_point[_qp]);
     }
   }
 }