DG cohesive zone model kernel for the Total Lagrangian formulation. More...

#include <CZMInterfaceKernelTotalLagrangian.h>

Inheritance diagram for CZMInterfaceKernelTotalLagrangian:

Public Types
typedef std::vector< int >	JvarMap

Public Member Functions
	CZMInterfaceKernelTotalLagrangian (const InputParameters &parameters)

virtual void	computeOffDiagJacobian (unsigned int jvar) override

unsigned int	mapJvarToCvar (unsigned int jvar)

int	mapJvarToCvar (unsigned int jvar, const JvarMap &jvar_map)

bool	mapJvarToCvar (unsigned int jvar, unsigned int &cvar)

const JvarMap &	getJvarMap ()

const JvarMap &	getParameterJvarMap (std::string parameter_name)

Static Public Member Functions
static InputParameters	validParams ()

Protected Member Functions
Real	computeDResidualDDisplacement (const unsigned int &component_j, const Moose::DGJacobianType &type) const override
	method computing the derivative of residual[_component] w.r.t displacement[component_j]

Real	JacLD (const unsigned int cc, const bool neighbor) const
	method computing the jacobian contribution due to rotations and area changes

Real	computeQpResidual (Moose::DGResidualType type) override

Real	computeQpJacobian (Moose::DGJacobianType type) override

Real	computeQpOffDiagJacobian (Moose::DGJacobianType type, unsigned int jvar) override

Protected Attributes
const MaterialProperty< RankThreeTensor > &	_dPK1traction_dF

const std::string	_base_name
	Base name of the material system that this kernel applies to.

const unsigned int	_component
	the displacement component this kernel is operating on (0=x, 1=y, 2 =z)

const unsigned int	_ndisp
	number of displacement components

std::vector< MooseVariable * >	_vars

const unsigned int	_n_args


std::vector< unsigned int >	_disp_var
	Coupled displacement component variable IDs.

std::vector< unsigned int >	_disp_neighbor_var


const MaterialProperty< RealVectorValue > &	_traction_global

const MaterialProperty< RankTwoTensor > &	_dtraction_djump_global

Private Attributes
const std::size_t	_jvar_max_size

JvarMap	_jvar_map

std::map< std::string, JvarMap >	_jvar_local_map

Detailed Description

DG cohesive zone model kernel for the Total Lagrangian formulation.

This kernel assummes the traction sepration law only depends from the displacement jump. One kernel is required for each displacement component

Definition at line 17 of file CZMInterfaceKernelTotalLagrangian.h.

Constructor & Destructor Documentation

◆ CZMInterfaceKernelTotalLagrangian()

CZMInterfaceKernelTotalLagrangian::CZMInterfaceKernelTotalLagrangian ( const InputParameters & parameters )

Definition at line 24 of file CZMInterfaceKernelTotalLagrangian.C.

  : CZMInterfaceKernelBase(parameters),
    _dPK1traction_dF(getMaterialPropertyByName<RankThreeTensor>(_base_name + "dPK1traction_dF"))
{
}

Member Function Documentation

◆ computeDResidualDDisplacement()

Real CZMInterfaceKernelTotalLagrangian::computeDResidualDDisplacement	(	const unsigned int &	component_j,
		const Moose::DGJacobianType &	type
	)		const

overrideprotectedvirtual

method computing the derivative of residual[_component] w.r.t displacement[component_j]

Implements CZMInterfaceKernelBase.

Definition at line 32 of file CZMInterfaceKernelTotalLagrangian.C.

{
  Real jacsd = _dtraction_djump_global[_qp](_component, component_j);
  Real jac = 0;
 
  switch (type)
  {
    case Moose::ElementElement: // Residual_sign -1  ddeltaU_ddisp sign -1;
      jac += _test[_i][_qp] * jacsd * _vars[component_j]->phiFace()[_j][_qp];
      jac -= _test[_i][_qp] * JacLD(component_j, /*neighbor=*/false);
      break;
    case Moose::ElementNeighbor: // Residual_sign -1  ddeltaU_ddisp sign 1;
      jac -= _test[_i][_qp] * jacsd * _vars[component_j]->phiFaceNeighbor()[_j][_qp];
      jac -= _test[_i][_qp] * JacLD(component_j, /*neighbor=*/true);
      break;
    case Moose::NeighborElement: // Residual_sign 1  ddeltaU_ddisp sign -1;
      jac -= _test_neighbor[_i][_qp] * jacsd * _vars[component_j]->phiFace()[_j][_qp];
      jac += _test_neighbor[_i][_qp] * JacLD(component_j, /*neighbor=*/false);
      break;
    case Moose::NeighborNeighbor: // Residual_sign 1  ddeltaU_ddisp sign 1;
      jac += _test_neighbor[_i][_qp] * jacsd * _vars[component_j]->phiFaceNeighbor()[_j][_qp];
      jac += _test_neighbor[_i][_qp] * JacLD(component_j, /*neighbor=*/true);
      break;
  }
 
  return jac;
}

◆ computeQpJacobian()

Real CZMInterfaceKernelBase::computeQpJacobian ( Moose::DGJacobianType type )

overrideprotectedinherited

Definition at line 79 of file CZMInterfaceKernelBase.C.

{
  // retrieve the diagonal Jacobian coefficient dependning on the displacement
  // component (_component) this kernel is working on
  return computeDResidualDDisplacement(_component, type);
}

◆ computeQpOffDiagJacobian()

Real CZMInterfaceKernelBase::computeQpOffDiagJacobian	(	Moose::DGJacobianType	type,
		unsigned int	jvar
	)

overrideprotectedinherited

Definition at line 87 of file CZMInterfaceKernelBase.C.

{
  // bail out if jvar is not coupled
  if (getJvarMap()[jvar] < 0)
    return 0.0;
 
  // Jacobian of the residul[_component] w.r.t to the coupled displacement
  // component[off_diag_component]
  for (unsigned int off_diag_component = 0; off_diag_component < _ndisp; ++off_diag_component)
  {
    if (jvar == _disp_var[off_diag_component])
      return computeDResidualDDisplacement(off_diag_component, type);
  }
  // this is the place where one should implement derivatives of the residual w.r.t. other variables
  return 0.0;
}

◆ computeQpResidual()

Real CZMInterfaceKernelBase::computeQpResidual ( Moose::DGResidualType type )

overrideprotectedinherited

Definition at line 61 of file CZMInterfaceKernelBase.C.

{
  Real r = _traction_global[_qp](_component);
 
  switch (type)
  {
    // [test_secondary-test_primary]*T where T represents the traction.
    case Moose::Element:
      r *= -_test[_i][_qp];
      break;
    case Moose::Neighbor:
      r *= _test_neighbor[_i][_qp];
      break;
  }
  return r;
}

◆ JacLD()

Real CZMInterfaceKernelTotalLagrangian::JacLD	(	const unsigned int	cc,
		const bool	neighbor
	)		const

protected

method computing the jacobian contribution due to rotations and area changes

Definition at line 62 of file CZMInterfaceKernelTotalLagrangian.C.

{
  Real jacld = 0;
  RealVectorValue phi;
  if (neighbor)
    phi = 0.5 * _vars[cc]->gradPhiFaceNeighbor()[_j][_qp];
  else
    phi = 0.5 * _vars[cc]->gradPhiFace()[_j][_qp];
 
  for (unsigned int j = 0; j < 3; j++)
    jacld += _dPK1traction_dF[_qp](_component, cc, j) * phi(j);
 
  mooseAssert(std::isfinite(jacld), "CZMInterfaceKernelTotalLagrangian JacLD is not finite");
  return jacld;
}

Referenced by computeDResidualDDisplacement().

◆ validParams()

InputParameters CZMInterfaceKernelTotalLagrangian::validParams ( )

static

Definition at line 15 of file CZMInterfaceKernelTotalLagrangian.C.

{
  InputParameters params = CZMInterfaceKernelBase::validParams();
  params.addClassDescription("Calculate residual contribution for balancing the traction across an "
                             "interface (used in the cohesive zone method).");
  params.set<std::string>("traction_global_name") = "PK1traction";
  return params;
}

Member Data Documentation

◆ _base_name

const std::string CZMInterfaceKernelBase::_base_name

protectedinherited

Base name of the material system that this kernel applies to.

Definition at line 35 of file CZMInterfaceKernelBase.h.

◆ _component

const unsigned int CZMInterfaceKernelBase::_component

protectedinherited

the displacement component this kernel is operating on (0=x, 1=y, 2 =z)

Definition at line 38 of file CZMInterfaceKernelBase.h.

Referenced by CZMInterfaceKernelSmallStrain::computeDResidualDDisplacement(), computeDResidualDDisplacement(), CZMInterfaceKernelBase::computeQpJacobian(), CZMInterfaceKernelBase::computeQpResidual(), and JacLD().

◆ _disp_neighbor_var

std::vector<unsigned int> CZMInterfaceKernelBase::_disp_neighbor_var

protectedinherited

Definition at line 46 of file CZMInterfaceKernelBase.h.

Referenced by CZMInterfaceKernelBase::CZMInterfaceKernelBase().

◆ _disp_var

std::vector<unsigned int> CZMInterfaceKernelBase::_disp_var

protectedinherited

Coupled displacement component variable IDs.

Definition at line 45 of file CZMInterfaceKernelBase.h.

Referenced by CZMInterfaceKernelBase::computeQpOffDiagJacobian(), and CZMInterfaceKernelBase::CZMInterfaceKernelBase().

◆ _dPK1traction_dF

const MaterialProperty<RankThreeTensor>& CZMInterfaceKernelTotalLagrangian::_dPK1traction_dF

protected

Definition at line 32 of file CZMInterfaceKernelTotalLagrangian.h.

Referenced by JacLD().

◆ _dtraction_djump_global

const MaterialProperty<RankTwoTensor>& CZMInterfaceKernelBase::_dtraction_djump_global

protectedinherited

Definition at line 55 of file CZMInterfaceKernelBase.h.

Referenced by CZMInterfaceKernelSmallStrain::computeDResidualDDisplacement(), and computeDResidualDDisplacement().

◆ _ndisp

const unsigned int CZMInterfaceKernelBase::_ndisp

protectedinherited

number of displacement components

Definition at line 41 of file CZMInterfaceKernelBase.h.

Referenced by CZMInterfaceKernelBase::computeQpOffDiagJacobian(), and CZMInterfaceKernelBase::CZMInterfaceKernelBase().

◆ _traction_global

const MaterialProperty<RealVectorValue>& CZMInterfaceKernelBase::_traction_global

protectedinherited

Definition at line 54 of file CZMInterfaceKernelBase.h.

Referenced by CZMInterfaceKernelBase::computeQpResidual().

◆ _vars

std::vector<MooseVariable *> CZMInterfaceKernelBase::_vars

protectedinherited

Definition at line 50 of file CZMInterfaceKernelBase.h.

Referenced by CZMInterfaceKernelSmallStrain::computeDResidualDDisplacement(), computeDResidualDDisplacement(), CZMInterfaceKernelBase::CZMInterfaceKernelBase(), and JacLD().

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

solid_mechanics/include/interfacekernels/CZMInterfaceKernelTotalLagrangian.h
solid_mechanics/src/interfacekernels/CZMInterfaceKernelTotalLagrangian.C

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Protected Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ CZMInterfaceKernelTotalLagrangian()

Member Function Documentation

◆ computeDResidualDDisplacement()

◆ computeQpJacobian()

◆ computeQpOffDiagJacobian()

◆ computeQpResidual()

◆ JacLD()

◆ validParams()

Member Data Documentation

◆ _base_name

◆ _component

◆ _disp_neighbor_var

◆ _disp_var

◆ _dPK1traction_dF

◆ _dtraction_djump_global

◆ _ndisp

◆ _traction_global

◆ _vars