Considers cleavage plane anisotropy in the crack propagation. More...

#include <ACInterfaceCleavageFracture.h>

Inheritance diagram for ACInterfaceCleavageFracture:

Public Types
typedef DerivativeMaterialPropertyNameInterface::SymbolName	SymbolName

Public Member Functions
	ACInterfaceCleavageFracture (const InputParameters &parameters)

virtual void	initialSetup ()

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialProperty (const std::string &name)

const GenericMaterialProperty< U, is_ad > &	getDefaultMaterialPropertyByName (const std::string &name)

void	validateDerivativeMaterialPropertyBase (const std::string &base)

const MaterialPropertyName	derivativePropertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	derivativePropertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	derivativePropertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	derivativePropertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

GenericMaterialProperty< U, is_ad > &	declarePropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, const SymbolName &c1, unsigned int v2, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivative (const std::string &base, unsigned int v1, unsigned int v2=libMesh::invalid_uint, unsigned int v3=libMesh::invalid_uint)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< VariableName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const std::vector< SymbolName > &c)

const GenericMaterialProperty< U, is_ad > &	getMaterialPropertyDerivativeByName (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2="", const SymbolName &c3="")

void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)

void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateCoupling (const MaterialPropertyName &base, const std::vector< VariableName > &c, bool validate_aux=true)

void	validateCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

void	validateNonlinearCoupling (const MaterialPropertyName &base, const VariableName &c1="", const VariableName &c2="", const VariableName &c3="")

const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	propertyName (const MaterialPropertyName &base, const std::vector< SymbolName > &c) const

const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	propertyNameFirst (const MaterialPropertyName &base, const SymbolName &c1) const

const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	propertyNameSecond (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2) const

const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

const MaterialPropertyName	propertyNameThird (const MaterialPropertyName &base, const SymbolName &c1, const SymbolName &c2, const SymbolName &c3) const

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
virtual Real	computeQpResidual ()

virtual Real	computeQpJacobian ()

Real	betaNablaPsi ()
	term with beta penalty

virtual Real	computeQpOffDiagJacobian (unsigned int jvar)

RealGradient	gradL ()

RealGradient	gradKappa ()

RealGradient	nablaLPsi ()
	the \( \nabla(L\psi) \) term

RealGradient	kappaNablaLPsi ()
	the \( \kappa\nabla(L\psi) \) term

Protected Attributes
const Real	_beta_penalty
	penalty for damage on planes not normal to the weak (favoured) cleavage plane (Clayton & Knap, 2015)

const RealVectorValue	_cleavage_plane_normal
	Plane normal to the weak cleavage plane: M in (Clayton & Knap, 2015)

const MaterialProperty< Real > &	_L
	Mobility.

const MaterialProperty< Real > &	_kappa
	Interfacial parameter.

const bool	_variable_L
	flag set if L is a function of non-linear variables in args

const MaterialProperty< Real > &	_dkappadop
	kappa derivative w.r.t. order parameter

std::vector< const MaterialProperty< Real > * >	_dkappadarg
	kappa derivative w.r.t. other coupled variables

std::vector< const VariableGradient * >	_gradarg
	Gradients for all coupled variables.


const MaterialProperty< Real > &	_dLdop
	Mobility derivatives w.r.t. order parameter.

const MaterialProperty< Real > &	_d2Ldop2


std::vector< const MaterialProperty< Real > * >	_dLdarg
	Mobility derivative w.r.t. other coupled variables.

std::vector< const MaterialProperty< Real > * >	_d2Ldargdop

std::vector< std::vector< const MaterialProperty< Real > * > >	_d2Ldarg2

Private Member Functions
bool	haveMaterialProperty (const std::string &prop_name)

std::vector< VariableName >	buildVariableVector (const VariableName &c1, const VariableName &c2, const VariableName &c3)

void	validateCouplingHelper (const MaterialPropertyName &base, const std::vector< VariableName > &c, const System &system, std::vector< VariableName > &missing)

bool	isNotObjectVariable (const VariableName &name)

Private Attributes
FEProblemBase &	_dmi_fe_problem

Detailed Description

Considers cleavage plane anisotropy in the crack propagation.

Definition at line 16 of file ACInterfaceCleavageFracture.h.

Constructor & Destructor Documentation

◆ ACInterfaceCleavageFracture()

ACInterfaceCleavageFracture::ACInterfaceCleavageFracture ( const InputParameters & parameters )

Definition at line 31 of file ACInterfaceCleavageFracture.C.

  : ACInterface(parameters),
    _beta_penalty(getParam<Real>("beta_penalty")),
    _cleavage_plane_normal(getParam<RealVectorValue>("cleavage_plane_normal"))
{
}

Member Function Documentation

◆ betaNablaPsi()

Real ACInterfaceCleavageFracture::betaNablaPsi ( )

protected

term with beta penalty

Definition at line 39 of file ACInterfaceCleavageFracture.C.

{
  return _beta_penalty * _L[_qp] * _kappa[_qp] * (_grad_u[_qp] * _cleavage_plane_normal) *
         (_grad_test[_i][_qp] * _cleavage_plane_normal);
}

Referenced by computeQpResidual().

◆ computeQpJacobian()

Real ACInterfaceCleavageFracture::computeQpJacobian ( )

protectedvirtual

dsum is the derivative \( \frac\partial{\partial \eta} \left( \nabla (L\psi) \right) \)

compute the derivative of the gradient of the mobility

Reimplemented from ACInterface.

Definition at line 52 of file ACInterfaceCleavageFracture.C.

{
  RealGradient dsum =
      (_dkappadop[_qp] * _L[_qp] + _kappa[_qp] * _dLdop[_qp]) * _phi[_j][_qp] * _grad_test[_i][_qp];
 
  if (_variable_L)
  {
    RealGradient dgradL =
        _grad_phi[_j][_qp] * _dLdop[_qp] + _grad_u[_qp] * _phi[_j][_qp] * _d2Ldop2[_qp];
 
    for (unsigned int i = 0; i < _n_args; ++i)
      dgradL += (*_gradarg[i])[_qp] * _phi[_j][_qp] * (*_d2Ldargdop[i])[_qp];
 
    dsum += (_kappa[_qp] * dgradL + _dkappadop[_qp] * _phi[_j][_qp] * gradL()) * _test[_i][_qp];
  }
 
  return (1 + _beta_penalty) * _grad_phi[_j][_qp] * kappaNablaLPsi() + _grad_u[_qp] * dsum -
         _beta_penalty * _L[_qp] * _kappa[_qp] * (_grad_u[_qp] * _cleavage_plane_normal) *
             (_grad_phi[_j][_qp] * _cleavage_plane_normal);
}

◆ computeQpOffDiagJacobian()

Real ACInterface::computeQpOffDiagJacobian ( unsigned int jvar )

protectedvirtualinherited

Reimplemented in ACInterface2DMultiPhase1, and ACInterface2DMultiPhase2.

Definition at line 128 of file ACInterface.C.

{
  // get the coupled variable jvar is referring to
  const unsigned int cvar = mapJvarToCvar(jvar);
 
  // dsum is the derivative \f$ \frac\partial{\partial \eta} \left( \nabla (L\psi) \right) \f$
  RealGradient dsum = ((*_dkappadarg[cvar])[_qp] * _L[_qp] + _kappa[_qp] * (*_dLdarg[cvar])[_qp]) *
                      _phi[_j][_qp] * _grad_test[_i][_qp];
 
  // compute the derivative of the gradient of the mobility
  if (_variable_L)
  {
    RealGradient dgradL = _grad_phi[_j][_qp] * (*_dLdarg[cvar])[_qp] +
                          _grad_u[_qp] * _phi[_j][_qp] * (*_d2Ldargdop[cvar])[_qp];
 
    for (unsigned int i = 0; i < _n_args; ++i)
      dgradL += (*_gradarg[i])[_qp] * _phi[_j][_qp] * (*_d2Ldarg2[cvar][i])[_qp];
 
    dsum += (_kappa[_qp] * dgradL + _dkappadop[_qp] * _phi[_j][_qp] * gradL()) * _test[_i][_qp];
  }
 
  return _grad_u[_qp] * dsum;
}

◆ computeQpResidual()

Real ACInterfaceCleavageFracture::computeQpResidual ( )

protectedvirtual

Reimplemented from ACInterface.

Definition at line 46 of file ACInterfaceCleavageFracture.C.

{
  return (1 + _beta_penalty) * _grad_u[_qp] * kappaNablaLPsi() - betaNablaPsi();
}

◆ gradKappa()

RealGradient ACInterface::gradKappa ( )

protectedinherited

◆ gradL()

RealGradient ACInterface::gradL ( )

protectedinherited

Definition at line 73 of file ACInterface.C.

{
  RealGradient g = _grad_u[_qp] * _dLdop[_qp];
  for (unsigned int i = 0; i < _n_args; ++i)
    g += (*_gradarg[i])[_qp] * (*_dLdarg[i])[_qp];
  return g;
}

Referenced by ACInterface::computeQpJacobian(), computeQpJacobian(), ACInterface::computeQpOffDiagJacobian(), and ACInterface::nablaLPsi().

◆ initialSetup()

void ACInterface::initialSetup ( )

virtualinherited

Reimplemented in ACInterfaceChangedVariable.

Definition at line 66 of file ACInterface.C.

{
  validateCoupling<Real>("mob_name");
  validateCoupling<Real>("kappa_name");
}

◆ kappaNablaLPsi()

RealGradient ACInterface::kappaNablaLPsi ( )

protectedinherited

the \( \kappa\nabla(L\psi) \) term

Definition at line 94 of file ACInterface.C.

{
  return _kappa[_qp] * nablaLPsi();
}

Referenced by ACInterface::computeQpJacobian(), computeQpJacobian(), ACInterface::computeQpResidual(), and computeQpResidual().

◆ nablaLPsi()

RealGradient ACInterface::nablaLPsi ( )

protectedinherited

the \( \nabla(L\psi) \) term

Definition at line 82 of file ACInterface.C.

{
  // sum is the product rule gradient \f$ \nabla (L\psi) \f$
  RealGradient sum = _L[_qp] * _grad_test[_i][_qp];
 
  if (_variable_L)
    sum += gradL() * _test[_i][_qp];
 
  return sum;
}

Referenced by ACInterface2DMultiPhase1::computeQpJacobian(), ACInterface2DMultiPhase2::computeQpJacobian(), ACInterface2DMultiPhase1::computeQpOffDiagJacobian(), ACInterface2DMultiPhase2::computeQpOffDiagJacobian(), ACInterface2DMultiPhase1::computeQpResidual(), and ACInterface::kappaNablaLPsi().

◆ validParams()

InputParameters ACInterfaceCleavageFracture::validParams ( )

static

Definition at line 17 of file ACInterfaceCleavageFracture.C.

{
  InputParameters params = ACInterface::validParams();
  params.addClassDescription("Gradient energy Allen-Cahn Kernel where crack propagation along weak"
                             "cleavage plane is preferred");
  params.addRequiredParam<Real>(
      "beta_penalty",
      "penalty to penalize fracture on planes not normal to one cleavage plane normal which is "
      "normal to weak cleavage plane. Setting beta=0 results in isotropic damage.");
  params.addRequiredParam<RealVectorValue>("cleavage_plane_normal",
                                           "Normal to the weak cleavage plane");
  return params;
}