Implementation of the non-stateful exponential traction separation law proposed by Salehani, Mohsen Khajeh and Irani, Nilgoon 2018. More...

#include <SalehaniIrani3DCTraction.h>

Inheritance diagram for SalehaniIrani3DCTraction:

[legend]

Public Member Functions
	SalehaniIrani3DCTraction (const InputParameters &parameters)

Protected Member Functions
virtual RealVectorValue	computeTraction () override
	method returning the traction in the interface coordinate system. More...

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

virtual void	computeQpProperties () override

Protected Attributes
const std::vector< Real >	_delta_u0
	the displacement jump associated to the maximum traction More...

const std::vector< Real >	_max_allowable_traction
	the vector representing the maximum allowed traction in each direction More...

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

Implementation of the non-stateful exponential traction separation law proposed by Salehani, Mohsen Khajeh and Irani, Nilgoon 2018.

Definition at line 20 of file SalehaniIrani3DCTraction.h.

Constructor & Destructor Documentation

◆ SalehaniIrani3DCTraction()

SalehaniIrani3DCTraction::SalehaniIrani3DCTraction ( const InputParameters & parameters )

Definition at line 33 of file SalehaniIrani3DCTraction.C.

   : CZMMaterialBase(parameters),
     _delta_u0({getParam<Real>("normal_gap_at_maximum_normal_traction"),
                std::sqrt(2) * getParam<Real>("tangential_gap_at_maximum_shear_traction"),
                std::sqrt(2) * getParam<Real>("tangential_gap_at_maximum_shear_traction")}),
     _max_allowable_traction({getParam<Real>("maximum_normal_traction"),
                              getParam<Real>("maximum_shear_traction"),
                              getParam<Real>("maximum_shear_traction")})
 {
 }

Member Function Documentation

◆ computeQpProperties()

void CZMMaterialBase::computeQpProperties ( )

overrideprotectedvirtualinherited

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());
 }

◆ computeTraction()

RealVectorValue SalehaniIrani3DCTraction::computeTraction ( )

overrideprotectedvirtual

method returning the traction in the interface coordinate system.

Implements CZMMaterialBase.

Definition at line 45 of file SalehaniIrani3DCTraction.C.

 {
  
   // The convention for ordering the traction is N, T, S, where N is the normal direction, and T and
   // S are two arbitrary tangential directions.
   RealVectorValue traction_local;
  
   // temporary containers for auxiliary calculations
   Real aa, x, exp_x, a_i, b_i;
   // index variable to avoid multiple redefinitions
   unsigned int i;
  
   x = 0;
   for (i = 0; i < 3; i++)
   {
     aa = _displacement_jump[_qp](i) / _delta_u0[i];
     if (i > 0)
       aa *= aa; // square for shear component
  
     x += aa;
   }
  
   exp_x = std::exp(-x);
  
   for (i = 0; i < 3; i++)
   {
     if (i == 0)
       aa = std::exp(1);
     else
       aa = std::sqrt(2 * std::exp(1));
  
     a_i = _max_allowable_traction[i] * aa;
     b_i = _displacement_jump[_qp](i) / _delta_u0[i];
     traction_local(i) = a_i * b_i * exp_x;
   }
  
   return traction_local;
 }

◆ computeTractionDerivatives()

RankTwoTensor SalehaniIrani3DCTraction::computeTractionDerivatives ( )

overrideprotectedvirtual

method returning the traction derivitaves wrt local displacement jump.

Implements CZMMaterialBase.

Definition at line 85 of file SalehaniIrani3DCTraction.C.

 {
   RankTwoTensor traction_jump_derivatives_local;
  
   // this function computes the partial derivates of Tn[0][:], Tt[1][:], Ts[2][:]
   // wrt dun, dut, dus
   // T_i = a_i*b_i*exp(-x) with:
   // a_i = \sigma_i,max * (\alpha_i*e)^{1/\alpha_i} with \alpha_i = 1 for i==n
   // \alpha_i = 2 for i!=n
   // b_i = \delta_u,i / \delta_0,i
   // x = sum_i=1^3{(\delta_u,i / \delta_0,i)^\alpha_i}  with \alpha_i = 1 for i==n
   // \alpha_i = 2 for i!=n
  
   // dTi_duj = a_i * ( dBi_duj * exp(-x) + b_i * exp(-x) * dx_duj  )
   //         = a_i * ( exp(-x) * (dBi_duj + b_i * dx_duj ) )
  
   // temporary containers for auxiliary calculations
   Real aa, exp_x, x;
   // index variables to avoid multiple redefinitions
   unsigned int i, j;
  
   // compute x and the exponential term
   aa = 0;
   x = 0;
   for (i = 0; i < 3; i++)
   {
     aa = _displacement_jump[_qp](i) / _delta_u0[i];
     if (i > 0)
       aa *= aa;
     x += aa;
   }
   exp_x = std::exp(-x);
  
   // compute partial derivatives in local coordiante wrt the displacement jump
   //            | dTn/dun dTn/dut dTn/dus |
   // dTi_duj  = | dTt/dun dTt/dut dTt/dus | = _traction_derivative[i][j]
   //            | dTs/dun dTs/dut dTs/dus |
   Real a_i, b_i;
   Real dbi_dui, dx_duj;
  
   for (i = 0; i < 3; i++)
   {
     if (i == 0) // alpha = 1
       a_i = std::exp(1);
     else // alpha = 2
       a_i = std::sqrt(2 * std::exp(1));
  
     a_i *= _max_allowable_traction[i];
     b_i = _displacement_jump[_qp](i) / _delta_u0[i];
  
     for (j = 0; j < 3; j++)
     {
  
       dbi_dui = 0;
       if (i == j)
         dbi_dui = 1 / _delta_u0[j];
  
       if (j == 0) // alpha = 1
         dx_duj = 1. / _delta_u0[j];
       else // alpha = 2
         dx_duj = 2. * _displacement_jump[_qp](j) / (_delta_u0[j] * _delta_u0[j]);
  
       traction_jump_derivatives_local(i, j) =
           a_i * exp_x * (dbi_dui - b_i * dx_duj); // the minus sign is due to exp(-x)
     }
   }
  
   return traction_jump_derivatives_local;
 }

Member Data Documentation

◆ _delta_u0

const std::vector<Real> SalehaniIrani3DCTraction::_delta_u0

protected

the displacement jump associated to the maximum traction

Definition at line 31 of file SalehaniIrani3DCTraction.h.

Referenced by computeTraction(), and computeTractionDerivatives().

◆ _disp

std::vector<const VariableValue *> CZMMaterialBase::_disp

protectedinherited

the coupled displacement and neighbor displacement values

Definition at line 45 of file CZMMaterialBase.h.

Referenced by CZMMaterialBase::computeQpProperties(), and CZMMaterialBase::CZMMaterialBase().

◆ _disp_neighbor

std::vector<const VariableValue *> CZMMaterialBase::_disp_neighbor

protectedinherited

Definition at line 46 of file CZMMaterialBase.h.

Referenced by CZMMaterialBase::computeQpProperties(), and CZMMaterialBase::CZMMaterialBase().

◆ _displacement_jump

MaterialProperty<RealVectorValue>& CZMMaterialBase::_displacement_jump

protectedinherited

Definition at line 58 of file CZMMaterialBase.h.

Referenced by CZMMaterialBase::computeQpProperties(), computeTraction(), and computeTractionDerivatives().

◆ _displacement_jump_global

MaterialProperty<RealVectorValue>& CZMMaterialBase::_displacement_jump_global

protectedinherited

the displacement jump in global and local coordiante

Definition at line 57 of file CZMMaterialBase.h.

Referenced by CZMMaterialBase::computeQpProperties().

◆ _max_allowable_traction

const std::vector<Real> SalehaniIrani3DCTraction::_max_allowable_traction

protected

the vector representing the maximum allowed traction in each direction

Definition at line 34 of file SalehaniIrani3DCTraction.h.

Referenced by computeTraction(), and computeTractionDerivatives().

◆ _ndisp

const unsigned int CZMMaterialBase::_ndisp

protectedinherited

number of displacement components

Definition at line 41 of file CZMMaterialBase.h.

Referenced by CZMMaterialBase::computeQpProperties(), and CZMMaterialBase::CZMMaterialBase().

◆ _normals

const MooseArray<Point>& CZMMaterialBase::_normals

protectedinherited

normal to the interface

Definition at line 38 of file CZMMaterialBase.h.

Referenced by CZMMaterialBase::computeQpProperties().

◆ _traction

MaterialProperty<RealVectorValue>& CZMMaterialBase::_traction

protectedinherited

Definition at line 64 of file CZMMaterialBase.h.

Referenced by CZMMaterialBase::computeQpProperties().

◆ _traction_derivatives

MaterialProperty<RankTwoTensor>& CZMMaterialBase::_traction_derivatives

protectedinherited

Definition at line 70 of file CZMMaterialBase.h.

Referenced by CZMMaterialBase::computeQpProperties().

◆ _traction_derivatives_global

MaterialProperty<RankTwoTensor>& CZMMaterialBase::_traction_derivatives_global

protectedinherited

the traction's derivatives wrt the displacement jump in global and local coordinates

Definition at line 69 of file CZMMaterialBase.h.

Referenced by CZMMaterialBase::computeQpProperties().

◆ _traction_global

MaterialProperty<RealVectorValue>& CZMMaterialBase::_traction_global

protectedinherited

the value of the traction in global and local coordinates

Definition at line 63 of file CZMMaterialBase.h.

Referenced by CZMMaterialBase::computeQpProperties().

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

tensor_mechanics/include/materials/czm/SalehaniIrani3DCTraction.h
tensor_mechanics/src/materials/czm/SalehaniIrani3DCTraction.C

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Constructor & Destructor Documentation

◆ SalehaniIrani3DCTraction()

Member Function Documentation

◆ computeQpProperties()

◆ computeTraction()

◆ computeTractionDerivatives()

Member Data Documentation

◆ _delta_u0

◆ _disp

◆ _disp_neighbor

◆ _displacement_jump

◆ _displacement_jump_global

◆ _max_allowable_traction

◆ _ndisp

◆ _normals

◆ _traction

◆ _traction_derivatives

◆ _traction_derivatives_global

◆ _traction_global