This is the base Material class for implementing a traction separation material model. More...

#include <CZMMaterialBase.h>

Inheritance diagram for CZMMaterialBase:

Public Member Functions
	CZMMaterialBase (const InputParameters &parameters)

Protected Member Functions
virtual void	computeQpProperties () override

virtual RealVectorValue	computeTraction ()=0
	method returning the traction in the interface coordinate system. More...

virtual RankTwoTensor	computeTractionDerivatives ()=0
	method returning the traction derivitaves wrt local displacement jump. More...

Protected Attributes
const MooseArray< Point > &	_normals
	normal to the interface More...

const unsigned int	_ndisp
	number of displacement components More...


std::vector< const VariableValue * >	_disp
	the coupled displacement and neighbor displacement values More...

std::vector< const VariableValue * >	_disp_neighbor


MaterialProperty< RealVectorValue > &	_displacement_jump_global
	the displacement jump in global and local coordiante More...

MaterialProperty< RealVectorValue > &	_displacement_jump


MaterialProperty< RealVectorValue > &	_traction_global
	the value of the traction in global and local coordinates More...

MaterialProperty< RealVectorValue > &	_traction


MaterialProperty< RankTwoTensor > &	_traction_derivatives_global
	the traction's derivatives wrt the displacement jump in global and local coordinates More...

MaterialProperty< RankTwoTensor > &	_traction_derivatives

Detailed Description

This is the base Material class for implementing a traction separation material model.

The responsibility of this class is to rotate the displacement jump from global to local coordinate and rotate back traction and traction derivatives. The local coordinate system assumes the following order: opening, tangential1, tangential2. Note that tangential1, tangential2 are arbitrary and therefore the interface assumes an in-plane isotropic behavior. By overriding computeTraction and computeTractionDerivatives in aderived class, different traction separation laws can be implemented. The computeTraction and computeTractionDerivatives methods assumes calculations are performed in the local frame. CZM laws should always be implemented in 3D even if they are going to be used in 2D or 1D simulations. This class assumes small deformations and that the traction separation law is only dependent upon the the displacement jump.

Definition at line 29 of file CZMMaterialBase.h.

Constructor & Destructor Documentation

◆ CZMMaterialBase()

CZMMaterialBase::CZMMaterialBase ( const InputParameters & parameters )

Definition at line 26 of file CZMMaterialBase.C.

   : InterfaceMaterial(parameters),
     _normals(_assembly.normals()),
     _ndisp(coupledComponents("displacements")),
     _disp(_ndisp),
     _disp_neighbor(_ndisp),
     _displacement_jump_global(declareProperty<RealVectorValue>("displacement_jump_global")),
     _displacement_jump(declareProperty<RealVectorValue>("displacement_jump")),
     _traction_global(declareProperty<RealVectorValue>("traction_global")),
     _traction(declareProperty<RealVectorValue>("traction")),
     _traction_derivatives_global(declareProperty<RankTwoTensor>("traction_derivatives_global")),
     _traction_derivatives(declareProperty<RankTwoTensor>("traction_derivatives"))
 {
   if (_ndisp > 3 || _ndisp < 1)
     mooseError("the CZM material requires 1, 2 or 3 displacement variables");
  
   if (getParam<bool>("use_displaced_mesh") == true)
     mooseError("This material cannot be used with use_displaced_mesh = true");
  
   // initializing the displacement vectors
   for (unsigned int i = 0; i < _ndisp; ++i)
   {
     _disp[i] = &coupledValue("displacements", i);
     _disp_neighbor[i] = &coupledNeighborValue("displacements", i);
   }
 }

Member Function Documentation

◆ computeQpProperties()

void CZMMaterialBase::computeQpProperties ( )

overrideprotectedvirtual

Definition at line 54 of file CZMMaterialBase.C.

 {
  
   RealTensorValue RotationGlobalToLocal =
       RotationMatrix::rotVec1ToVec2(_normals[_qp], RealVectorValue(1, 0, 0));
  
   // computing the displacement jump
   for (unsigned int i = 0; i < _ndisp; i++)
     _displacement_jump_global[_qp](i) = (*_disp_neighbor[i])[_qp] - (*_disp[i])[_qp];
   for (unsigned int i = _ndisp; i < 3; i++)
     _displacement_jump_global[_qp](i) = 0;
  
   // rotate the displacement jump to the local coordiante system
   _displacement_jump[_qp] = RotationGlobalToLocal * _displacement_jump_global[_qp];
  
   // compute local traction
   _traction[_qp] = computeTraction();
  
   // compute local traction derivatives wrt the displacement jump
   _traction_derivatives[_qp] = computeTractionDerivatives();
  
   // rotate local traction and derivatives to the global coordinate system
   _traction_global[_qp] = RotationGlobalToLocal.transpose() * _traction[_qp];
   _traction_derivatives_global[_qp] = _traction_derivatives[_qp];
   _traction_derivatives_global[_qp].rotate(RotationGlobalToLocal.transpose());
 }