ADDensity< compute_stage > Class Template Reference

#include <ADDensity.h>

Inheritance diagram for ADDensity< compute_stage >:

Public Member Functions
	ADDensity (const InputParameters &parameters)

Protected Member Functions
virtual void	initQpStatefulProperties ()
virtual void	computeQpProperties ()

Private Member Functions
	ADMaterialProperty (Real) &_density

Private Attributes
Moose::CoordinateSystemType	_coord_system
std::vector< const ADVariableGradient * >	_grad_disp
const ADVariableValue &	_disp_r
const Real	_initial_density
	usingMaterialMembers

Detailed Description

template<ComputeStage compute_stage>
class ADDensity< compute_stage >

Definition at line 16 of file ADDensity.h.

Constructor & Destructor Documentation

ADDensity()

template<ComputeStage compute_stage>

ADDensity< compute_stage >::ADDensity

(

const InputParameters &

parameters

)

Definition at line 22 of file ADDensity.C.

  : ADMaterial<compute_stage>(parameters),
    _coord_system(getBlockCoordSystem()),
    _disp_r(adCoupledValue("displacements", 0)),
    _initial_density(getParam<Real>("density")),
    _density(declareADProperty<Real>("density"))
{
  // get coupled gradients
  const unsigned int ndisp = coupledComponents("displacements");
  _grad_disp.resize(ndisp);
  for (unsigned int i = 0; i < ndisp; ++i)
    _grad_disp[i] = &adCoupledGradient("displacements", i);
 
  // fill remaining components with zero
  _grad_disp.resize(3, &adZeroGradient());
}

Member Function Documentation

ADMaterialProperty()

template<ComputeStage compute_stage>

ADDensity< compute_stage >::ADMaterialProperty ( Real  ) &

private

computeQpProperties()

template<ComputeStage compute_stage>

void ADDensity< compute_stage >::computeQpProperties ( )

protectedvirtual

Definition at line 48 of file ADDensity.C.

{
  // rho * V = rho0 * V0
  // rho = rho0 * V0 / V
  // rho = rho0 / det(F)
  // rho = rho0 / det(grad(u) + 1)
 
  const auto Axx = (*_grad_disp[0])[_qp](0) + 1.0;
  const auto & Axy = (*_grad_disp[0])[_qp](1);
  const auto & Axz = (*_grad_disp[0])[_qp](2);
  const auto & Ayx = (*_grad_disp[1])[_qp](0);
  auto Ayy = (*_grad_disp[1])[_qp](1) + 1.0;
  const auto & Ayz = (*_grad_disp[1])[_qp](2);
  const auto & Azx = (*_grad_disp[2])[_qp](0);
  const auto & Azy = (*_grad_disp[2])[_qp](1);
  auto Azz = (*_grad_disp[2])[_qp](2) + 1.0;
 
  switch (_coord_system)
  {
    case Moose::COORD_XYZ:
      Azz = (*_grad_disp[2])[_qp](2) + 1.0;
      break;
 
    case Moose::COORD_RZ:
      if (_q_point[_qp](0) != 0.0)
        Azz = _disp_r[_qp] / _q_point[_qp](0) + 1.0;
      break;
 
    case Moose::COORD_RSPHERICAL:
      if (_q_point[_qp](0) != 0.0)
        Ayy = Azz = _disp_r[_qp] / _q_point[_qp](0) + 1.0;
      break;
  }
 
  const auto detF = Axx * Ayy * Azz + Axy * Ayz * Azx + Axz * Ayx * Azy - Azx * Ayy * Axz -
                    Azy * Ayz * Axx - Azz * Ayx * Axy;
  _density[_qp] = _initial_density / detF;
}

initQpStatefulProperties()

template<ComputeStage compute_stage>

void ADDensity< compute_stage >::initQpStatefulProperties ( )

protectedvirtual

Definition at line 41 of file ADDensity.C.

{
  _density[_qp] = _initial_density;
}

Member Data Documentation

_coord_system

template<ComputeStage compute_stage>

Moose::CoordinateSystemType ADDensity< compute_stage >::_coord_system

private

Definition at line 31 of file ADDensity.h.

_disp_r

template<ComputeStage compute_stage>

const ADVariableValue& ADDensity< compute_stage >::_disp_r

private

Definition at line 33 of file ADDensity.h.

_grad_disp

template<ComputeStage compute_stage>

std::vector<const ADVariableGradient *> ADDensity< compute_stage >::_grad_disp

private

Definition at line 32 of file ADDensity.h.

Referenced by ADDensity< compute_stage >::ADDensity().

_initial_density

template<ComputeStage compute_stage>

const Real ADDensity< compute_stage >::_initial_density

private

Definition at line 35 of file ADDensity.h.

usingMaterialMembers

template<ComputeStage compute_stage>

ADDensity< compute_stage >::usingMaterialMembers

private

Definition at line 38 of file ADDensity.h.

---

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

• misc/include/materials/ADDensity.h
• misc/src/materials/ADDensity.C