CZMInterfaceKernel Class Reference

DG kernel implementing cohesive zone models (CZM) for a 1D/2D/3D traction separation laws based on the displacement jump. More...

#include <CZMInterfaceKernel.h>

Inheritance diagram for CZMInterfaceKernel:

Public Member Functions
	CZMInterfaceKernel (const InputParameters &parameters)

Protected Member Functions
virtual Real	computeQpResidual (Moose::DGResidualType type)
virtual Real	computeQpJacobian (Moose::DGJacobianType type)
virtual Real	computeQpOffDiagJacobian (Moose::DGJacobianType type, unsigned int jvar)

Protected Attributes
const unsigned int	_component
	the displacement component this kernel is operating on (0=x, 1=y, 2 =z) More...
const unsigned int	_ndisp
	number of displacement components More...
std::vector< unsigned int >	_disp_var
	Coupled displacement component variable IDs. More...
std::vector< unsigned int >	_disp_neighbor_var
const MaterialProperty< RealVectorValue > &	_traction_global
const MaterialProperty< RankTwoTensor > &	_traction_derivatives_global

Detailed Description

DG kernel implementing cohesive zone models (CZM) for a 1D/2D/3D traction separation laws based on the displacement jump.

This kernel operates only on a single displacement compenent. One kernel is required for each mesh dimension.

Definition at line 24 of file CZMInterfaceKernel.h.

Constructor & Destructor Documentation

CZMInterfaceKernel()

CZMInterfaceKernel::CZMInterfaceKernel

(

const InputParameters &

parameters

)

Definition at line 34 of file CZMInterfaceKernel.C.

  : InterfaceKernel(parameters),
    _component(getParam<unsigned int>("component")),
    _ndisp(coupledComponents("displacements")),
    _disp_var(_ndisp),
    _disp_neighbor_var(_ndisp),
    _traction_global(getMaterialPropertyByName<RealVectorValue>("traction_global")),
    _traction_derivatives_global(
        getMaterialPropertyByName<RankTwoTensor>("traction_derivatives_global"))
{
  if (getParam<bool>("use_displaced_mesh") == true)
    mooseError("CZMInterfaceKernel cannot be used with use_displaced_mesh = true");
 
  for (unsigned int i = 0; i < _ndisp; ++i)
  {
    _disp_var[i] = coupled("displacements", i);
    _disp_neighbor_var[i] = coupled("displacements", i);
  }
}

Member Function Documentation

computeQpJacobian()

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

protectedvirtual

Definition at line 74 of file CZMInterfaceKernel.C.

{
  // retrieve the diagonal Jacobian coefficient dependning on the displacement
  // component (_component) this kernel is working on
  Real jac = _traction_derivatives_global[_qp](_component, _component);
 
  switch (type)
  {
    case Moose::ElementElement: // Residual_sign -1  ddeltaU_ddisp sign -1;
      jac *= _test[_i][_qp] * _phi[_j][_qp];
      break;
    case Moose::ElementNeighbor: // Residual_sign -1  ddeltaU_ddisp sign 1;
      jac *= -_test[_i][_qp] * _phi_neighbor[_j][_qp];
      break;
    case Moose::NeighborElement: // Residual_sign 1  ddeltaU_ddisp sign -1;
      jac *= -_test_neighbor[_i][_qp] * _phi[_j][_qp];
      break;
    case Moose::NeighborNeighbor: // Residual_sign 1  ddeltaU_ddisp sign 1;
      jac *= _test_neighbor[_i][_qp] * _phi_neighbor[_j][_qp];
      break;
  }
  return jac;
}

computeQpOffDiagJacobian()

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

protectedvirtual

Definition at line 99 of file CZMInterfaceKernel.C.

{
 
  // find the displacement component associated to jvar
  unsigned int off_diag_component;
  for (off_diag_component = 0; off_diag_component < _ndisp; off_diag_component++)
    if (_disp_var[off_diag_component] == jvar)
      break;
 
  mooseAssert(off_diag_component < _ndisp,
              "CZMInterfaceKernel::computeQpOffDiagJacobian wrong offdiagonal variable");
 
  Real jac = _traction_derivatives_global[_qp](_component, off_diag_component);
 
  switch (type)
  {
    case Moose::ElementElement: // Residual_sign -1  ddeltaU_ddisp sign -1;
      jac *= _test[_i][_qp] * _phi[_j][_qp];
      break;
    case Moose::ElementNeighbor: // Residual_sign -1  ddeltaU_ddisp sign 1;
      jac *= -_test[_i][_qp] * _phi_neighbor[_j][_qp];
      break;
    case Moose::NeighborElement: // Residual_sign 1  ddeltaU_ddisp sign -1;
      jac *= -_test_neighbor[_i][_qp] * _phi[_j][_qp];
      break;
    case Moose::NeighborNeighbor: // Residual_sign 1  ddeltaU_ddisp sign 1;
      jac *= _test_neighbor[_i][_qp] * _phi_neighbor[_j][_qp];
      break;
  }
  return jac;
}

computeQpResidual()

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

protectedvirtual

Definition at line 55 of file CZMInterfaceKernel.C.

{
 
  Real r = _traction_global[_qp](_component);
 
  switch (type)
  {
    // [test_slave-test_master]*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;
}

Member Data Documentation

_component

const unsigned int CZMInterfaceKernel::_component

protected

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

Definition at line 35 of file CZMInterfaceKernel.h.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), and computeQpResidual().

_disp_neighbor_var

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

protected

Definition at line 43 of file CZMInterfaceKernel.h.

Referenced by CZMInterfaceKernel().

_disp_var

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

protected

Coupled displacement component variable IDs.

Definition at line 42 of file CZMInterfaceKernel.h.

Referenced by computeQpOffDiagJacobian(), and CZMInterfaceKernel().

_ndisp

const unsigned int CZMInterfaceKernel::_ndisp

protected

number of displacement components

Definition at line 38 of file CZMInterfaceKernel.h.

Referenced by computeQpOffDiagJacobian(), and CZMInterfaceKernel().

_traction_derivatives_global

const MaterialProperty<RankTwoTensor>& CZMInterfaceKernel::_traction_derivatives_global

protected

Definition at line 49 of file CZMInterfaceKernel.h.

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

_traction_global

const MaterialProperty<RealVectorValue>& CZMInterfaceKernel::_traction_global

protected

Definition at line 48 of file CZMInterfaceKernel.h.

Referenced by computeQpResidual().

---

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

• tensor_mechanics/include/interfacekernels/CZMInterfaceKernel.h
• tensor_mechanics/src/interfacekernels/CZMInterfaceKernel.C