CrackTipEnrichmentStressDivergenceTensors implements the residual and jacobian for enrichement displacement variables. More...

#include <CrackTipEnrichmentStressDivergenceTensors.h>

Inheritance diagram for CrackTipEnrichmentStressDivergenceTensors:

[legend]

Public Member Functions
	CrackTipEnrichmentStressDivergenceTensors (const InputParameters &parameters)

virtual void	computeJacobian () override

virtual void	computeOffDiagJacobian (MooseVariableFEBase &jvar) override

virtual unsigned int	crackTipEnrichementFunctionAtPoint (const Point &point, std::vector< Real > &B)
	calculate the enrichment function values at point More...

virtual unsigned int	crackTipEnrichementFunctionDerivativeAtPoint (const Point &point, std::vector< RealVectorValue > &dB)
	calculate the enrichment function derivatives at point More...

void	rotateFromCrackFrontCoordsToGlobal (const RealVectorValue &vector, RealVectorValue &rotated_vector, const unsigned int point_index)
	rotate a vector from crack front coordinate to global cooridate More...

Protected Member Functions
virtual Real	computeQpResidual () override

virtual Real	computeQpJacobian () override

virtual Real	computeQpOffDiagJacobian (unsigned int jvar) override

Protected Attributes
std::string	_base_name

const MaterialProperty< RankTwoTensor > &	_stress

const MaterialProperty< RankFourTensor > &	_Jacobian_mult

const MaterialProperty< RankTwoTensor > *	_deformation_gradient

const MaterialProperty< RankTwoTensor > *	_deformation_gradient_old

const MaterialProperty< RankTwoTensor > *	_rotation_increment

const unsigned int	_component
	displacement components More...

const unsigned int	_enrichment_component
	enrichment function components More...

unsigned int	_nenrich_disp
	Coupled enrichment displacement variables. More...

std::vector< unsigned int >	_enrich_disp_var

unsigned int	_ndisp
	Coupled displacement variables. More...

std::vector< unsigned int >	_disp_var

Assembly &	_assembly_undisplaced
	undisplaced problem More...

MooseVariable &	_var_undisplaced
	Reference to this Kernel's undisplaced MooseVariable object. More...


const VariablePhiGradient &	_grad_phi_undisplaced
	Shape and test functions on the undisplaced mesh. More...

const VariableTestGradient &	_grad_test_undisplaced

Private Attributes
std::vector< Real >	_B
	enrichment function value More...

std::vector< RealVectorValue >	_dBX
	derivatives of enrichment function respect to global cooridnate More...

std::vector< RealVectorValue >	_dBx
	derivatives of enrichment function respect to crack front cooridnate More...

std::vector< Real >	_BI
	enrichment function at node I More...

std::vector< Real >	_BJ
	enrichment function at node J More...

Detailed Description

CrackTipEnrichmentStressDivergenceTensors implements the residual and jacobian for enrichement displacement variables.

Definition at line 32 of file CrackTipEnrichmentStressDivergenceTensors.h.

Constructor & Destructor Documentation

◆ CrackTipEnrichmentStressDivergenceTensors()

CrackTipEnrichmentStressDivergenceTensors::CrackTipEnrichmentStressDivergenceTensors ( const InputParameters & parameters )

Definition at line 38 of file CrackTipEnrichmentStressDivergenceTensors.C.

   : ALEKernel(parameters),
     EnrichmentFunctionCalculation(&getUserObject<CrackFrontDefinition>("crack_front_definition")),
     _base_name(isParamValid("base_name") ? getParam<std::string>("base_name") + "_" : ""),
     _stress(getMaterialPropertyByName<RankTwoTensor>(_base_name + "stress")),
     _Jacobian_mult(getMaterialPropertyByName<RankFourTensor>(_base_name + "Jacobian_mult")),
     _component(getParam<unsigned int>("component")),
     _enrichment_component(getParam<unsigned int>("enrichment_component")),
     _nenrich_disp(coupledComponents("enrichment_displacements")),
     _ndisp(coupledComponents("displacements")),
     _B(4),
     _dBX(4),
     _dBx(4),
     _BI(4),
     _BJ(4)
 {
   _enrich_disp_var.resize(_nenrich_disp);
   for (unsigned int i = 0; i < _nenrich_disp; ++i)
     _enrich_disp_var[i] = coupled("enrichment_displacements", i);
 
   _disp_var.resize(_ndisp);
   for (unsigned int i = 0; i < _ndisp; ++i)
     _disp_var[i] = coupled("displacements", i);
 }

Member Function Documentation

◆ computeJacobian()

void ALEKernel::computeJacobian ( )

overridevirtualinherited

Reimplemented in StressDivergenceTensors.

Definition at line 35 of file ALEKernel.C.

Referenced by StressDivergenceTensors::computeJacobian().

 {
   _fe_problem.prepareShapes(_var.number(), _tid);
   Kernel::computeJacobian();
 }

◆ computeOffDiagJacobian()

void ALEKernel::computeOffDiagJacobian ( MooseVariableFEBase & jvar )

overridevirtualinherited

Reimplemented in StressDivergenceTensors.

Definition at line 42 of file ALEKernel.C.

Referenced by StressDivergenceTensors::computeOffDiagJacobian().

 {
   _fe_problem.prepareShapes(jvar.number(), _tid);
   Kernel::computeOffDiagJacobian(jvar);
 }

◆ computeQpJacobian()

Real CrackTipEnrichmentStressDivergenceTensors::computeQpJacobian ( )

overrideprotectedvirtual

Definition at line 84 of file CrackTipEnrichmentStressDivergenceTensors.C.

 {
   crackTipEnrichementFunctionAtPoint(*_current_elem->get_node(_i), _BI);
   crackTipEnrichementFunctionAtPoint(*_current_elem->get_node(_j), _BJ);
 
   crackTipEnrichementFunctionAtPoint(_q_point[_qp], _B);
 
   unsigned int crack_front_point_index =
       crackTipEnrichementFunctionDerivativeAtPoint(_q_point[_qp], _dBx);
 
   for (unsigned int i = 0; i < 4; ++i)
     rotateFromCrackFrontCoordsToGlobal(_dBx[i], _dBX[i], crack_front_point_index);
 
   RealVectorValue grad_B(_dBX[_enrichment_component]);
 
   RealVectorValue grad_test =
       _grad_test[_i][_qp] * (_B[_enrichment_component] - _BI[_enrichment_component]) +
       _test[_i][_qp] * grad_B;
   RealVectorValue grad_phi =
       _grad_phi[_j][_qp] * (_B[_enrichment_component] - _BJ[_enrichment_component]) +
       _phi[_j][_qp] * grad_B;
 
   return ElasticityTensorTools::elasticJacobian(
       _Jacobian_mult[_qp], _component, _component, grad_test, grad_phi);
 }

◆ computeQpOffDiagJacobian()

Real CrackTipEnrichmentStressDivergenceTensors::computeQpOffDiagJacobian ( unsigned int jvar )

overrideprotectedvirtual

Definition at line 111 of file CrackTipEnrichmentStressDivergenceTensors.C.

 {
   unsigned int coupled_component = 0;
   unsigned int coupled_enrichment_component = 0;
   bool active(false);
   bool active_enrich(false);
 
   for (unsigned int i = 0; i < _enrich_disp_var.size(); ++i)
   {
     if (jvar == _enrich_disp_var[i])
     {
       coupled_component = i / 4;
       coupled_enrichment_component = i % 4;
       active_enrich = true;
     }
   }
 
   for (unsigned int i = 0; i < _disp_var.size(); ++i)
   {
     if (jvar == _disp_var[i])
     {
       coupled_component = i;
       active = true;
     }
   }
 
   if (active_enrich)
   {
     crackTipEnrichementFunctionAtPoint(*_current_elem->get_node(_i), _BI);
     crackTipEnrichementFunctionAtPoint(*_current_elem->get_node(_j), _BJ);
 
     crackTipEnrichementFunctionAtPoint(_q_point[_qp], _B);
     unsigned int crack_front_point_index =
         crackTipEnrichementFunctionDerivativeAtPoint(_q_point[_qp], _dBx);
 
     for (unsigned int i = 0; i < 4; ++i)
       rotateFromCrackFrontCoordsToGlobal(_dBx[i], _dBX[i], crack_front_point_index);
 
     RealVectorValue grad_B_test(_dBX[_enrichment_component]);
     RealVectorValue grad_B_phi(_dBX[coupled_enrichment_component]);
 
     RealVectorValue grad_test =
         _grad_test[_i][_qp] * (_B[_enrichment_component] - _BI[_enrichment_component]) +
         _test[_i][_qp] * grad_B_test;
     RealVectorValue grad_phi = _grad_phi[_j][_qp] * (_B[coupled_enrichment_component] -
                                                      _BJ[coupled_enrichment_component]) +
                                _phi[_j][_qp] * grad_B_phi;
 
     return ElasticityTensorTools::elasticJacobian(
         _Jacobian_mult[_qp], _component, coupled_component, grad_test, grad_phi);
   }
   else if (active)
   {
     crackTipEnrichementFunctionAtPoint(*_current_elem->get_node(_i), _BI);
 
     crackTipEnrichementFunctionAtPoint(_q_point[_qp], _B);
     unsigned int crack_front_point_index =
         crackTipEnrichementFunctionDerivativeAtPoint(_q_point[_qp], _dBx);
 
     for (unsigned int i = 0; i < 4; ++i)
       rotateFromCrackFrontCoordsToGlobal(_dBx[i], _dBX[i], crack_front_point_index);
 
     RealVectorValue grad_B_test(_dBX[_enrichment_component]);
 
     RealVectorValue grad_test =
         _grad_test[_i][_qp] * (_B[_enrichment_component] - _BI[_enrichment_component]) +
         _test[_i][_qp] * grad_B_test;
 
     return ElasticityTensorTools::elasticJacobian(
         _Jacobian_mult[_qp], _component, coupled_component, grad_test, _grad_phi[_j][_qp]);
   }
 
   return 0;
 }

◆ computeQpResidual()

Real CrackTipEnrichmentStressDivergenceTensors::computeQpResidual ( )

overrideprotectedvirtual

Definition at line 65 of file CrackTipEnrichmentStressDivergenceTensors.C.

 {
   crackTipEnrichementFunctionAtPoint(*_current_elem->get_node(_i), _BI);
 
   crackTipEnrichementFunctionAtPoint(_q_point[_qp], _B);
   unsigned int crack_front_point_index =
       crackTipEnrichementFunctionDerivativeAtPoint(_q_point[_qp], _dBx);
 
   for (unsigned int i = 0; i < 4; ++i)
     rotateFromCrackFrontCoordsToGlobal(_dBx[i], _dBX[i], crack_front_point_index);
 
   RealVectorValue grad_B(_dBX[_enrichment_component]);
 
   return _stress[_qp].row(_component) *
          (_grad_test[_i][_qp] * (_B[_enrichment_component] - _BI[_enrichment_component]) +
           _test[_i][_qp] * grad_B);
 }

◆ crackTipEnrichementFunctionAtPoint()

unsigned int EnrichmentFunctionCalculation::crackTipEnrichementFunctionAtPoint	(	const Point &	point,
		std::vector< Real > &	B
	)

virtualinherited

calculate the enrichment function values at point

Returns: the closest crack front index

Definition at line 19 of file EnrichmentFunctionCalculation.C.

Referenced by ComputeCrackTipEnrichmentSmallStrain::computeProperties(), computeQpJacobian(), computeQpOffDiagJacobian(), ComputeCrackTipEnrichmentSmallStrain::computeQpProperties(), and computeQpResidual().

 {
   unsigned int crack_front_point_index =
       _crack_front_definition->calculateRThetaToCrackFront(point, _r, _theta);
 
   if (MooseUtils::absoluteFuzzyEqual(_r, 0.0))
     mooseError("EnrichmentFunctionCalculation: the distance between a point and the crack "
                "tip/front is zero.");
 
   Real st = std::sin(_theta);
   Real st2 = std::sin(_theta / 2.0);
   Real ct2 = std::cos(_theta / 2.0);
   Real sr = std::sqrt(_r);
 
   B[0] = sr * st2;
   B[1] = sr * ct2;
   B[2] = sr * st2 * st;
   B[3] = sr * ct2 * st;
 
   return crack_front_point_index;
 }

◆ crackTipEnrichementFunctionDerivativeAtPoint()

unsigned int EnrichmentFunctionCalculation::crackTipEnrichementFunctionDerivativeAtPoint	(	const Point &	point,
		std::vector< RealVectorValue > &	dB
	)

virtualinherited

calculate the enrichment function derivatives at point

Returns: the closest crack front index

Definition at line 43 of file EnrichmentFunctionCalculation.C.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), ComputeCrackTipEnrichmentSmallStrain::computeQpProperties(), and computeQpResidual().

 {
   unsigned int crack_front_point_index =
       _crack_front_definition->calculateRThetaToCrackFront(point, _r, _theta);
 
   if (MooseUtils::absoluteFuzzyEqual(_r, 0.0))
     mooseError("EnrichmentFunctionCalculation: the distance between a point and the crack "
                "tip/front is zero.");
 
   Real st = std::sin(_theta);
   Real ct = std::cos(_theta);
   Real st2 = std::sin(_theta / 2.0);
   Real ct2 = std::cos(_theta / 2.0);
   Real st15 = std::sin(1.5 * _theta);
   Real ct15 = std::cos(1.5 * _theta);
   Real sr = std::sqrt(_r);
 
   dB[0](0) = -0.5 / sr * st2;
   dB[0](1) = 0.5 / sr * ct2;
   dB[0](2) = 0.0;
   dB[1](0) = 0.5 / sr * ct2;
   dB[1](1) = 0.5 / sr * st2;
   dB[1](2) = 0.0;
   dB[2](0) = -0.5 / sr * st15 * st;
   dB[2](1) = 0.5 / sr * (st2 + st15 * ct);
   dB[2](2) = 0.0;
   dB[3](0) = -0.5 / sr * ct15 * st;
   dB[3](1) = 0.5 / sr * (ct2 + ct15 * ct);
   dB[3](2) = 0.0;
 
   return crack_front_point_index;
 }

◆ rotateFromCrackFrontCoordsToGlobal()

void EnrichmentFunctionCalculation::rotateFromCrackFrontCoordsToGlobal	(	const RealVectorValue &	vector,
		RealVectorValue &	rotated_vector,
		const unsigned int	point_index
	)

inherited

rotate a vector from crack front coordinate to global cooridate

Parameters

rotated_vector rotated vector

Definition at line 78 of file EnrichmentFunctionCalculation.C.

Referenced by computeQpJacobian(), computeQpOffDiagJacobian(), ComputeCrackTipEnrichmentSmallStrain::computeQpProperties(), and computeQpResidual().

 {
   rotated_vector = _crack_front_definition->rotateFromCrackFrontCoordsToGlobal(vector, point_index);
 }

Member Data Documentation

◆ _assembly_undisplaced

Assembly& ALEKernel::_assembly_undisplaced

protectedinherited

undisplaced problem

Definition at line 33 of file ALEKernel.h.

◆ _B

std::vector<Real> CrackTipEnrichmentStressDivergenceTensors::_B

private

enrichment function value

Definition at line 67 of file CrackTipEnrichmentStressDivergenceTensors.h.

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

◆ _base_name

std::string CrackTipEnrichmentStressDivergenceTensors::_base_name

protected

Definition at line 43 of file CrackTipEnrichmentStressDivergenceTensors.h.

◆ _BI

std::vector<Real> CrackTipEnrichmentStressDivergenceTensors::_BI

private

enrichment function at node I

Definition at line 73 of file CrackTipEnrichmentStressDivergenceTensors.h.

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

◆ _BJ

std::vector<Real> CrackTipEnrichmentStressDivergenceTensors::_BJ

private

enrichment function at node J

Definition at line 75 of file CrackTipEnrichmentStressDivergenceTensors.h.

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

◆ _component

const unsigned int CrackTipEnrichmentStressDivergenceTensors::_component

protected

displacement components

Definition at line 53 of file CrackTipEnrichmentStressDivergenceTensors.h.

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

◆ _dBX

std::vector<RealVectorValue> CrackTipEnrichmentStressDivergenceTensors::_dBX

private

derivatives of enrichment function respect to global cooridnate

Definition at line 69 of file CrackTipEnrichmentStressDivergenceTensors.h.

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

◆ _dBx

std::vector<RealVectorValue> CrackTipEnrichmentStressDivergenceTensors::_dBx

private

derivatives of enrichment function respect to crack front cooridnate

Definition at line 71 of file CrackTipEnrichmentStressDivergenceTensors.h.

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

◆ _deformation_gradient

const MaterialProperty<RankTwoTensor>* CrackTipEnrichmentStressDivergenceTensors::_deformation_gradient

protected

Definition at line 48 of file CrackTipEnrichmentStressDivergenceTensors.h.

◆ _deformation_gradient_old

const MaterialProperty<RankTwoTensor>* CrackTipEnrichmentStressDivergenceTensors::_deformation_gradient_old

protected

Definition at line 49 of file CrackTipEnrichmentStressDivergenceTensors.h.

◆ _disp_var

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

protected

Definition at line 63 of file CrackTipEnrichmentStressDivergenceTensors.h.

Referenced by computeQpOffDiagJacobian(), and CrackTipEnrichmentStressDivergenceTensors().

◆ _enrich_disp_var

std::vector<unsigned int> CrackTipEnrichmentStressDivergenceTensors::_enrich_disp_var

protected

Definition at line 59 of file CrackTipEnrichmentStressDivergenceTensors.h.

Referenced by computeQpOffDiagJacobian(), and CrackTipEnrichmentStressDivergenceTensors().

◆ _enrichment_component

const unsigned int CrackTipEnrichmentStressDivergenceTensors::_enrichment_component

protected

enrichment function components

Definition at line 55 of file CrackTipEnrichmentStressDivergenceTensors.h.

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

◆ _grad_phi_undisplaced

const VariablePhiGradient& ALEKernel::_grad_phi_undisplaced

protectedinherited

Shape and test functions on the undisplaced mesh.

Definition at line 39 of file ALEKernel.h.

Referenced by StressDivergenceTensors::computeQpJacobian(), and StressDivergenceTensors::computeQpOffDiagJacobian().

◆ _grad_test_undisplaced

const VariableTestGradient& ALEKernel::_grad_test_undisplaced

protectedinherited

Definition at line 40 of file ALEKernel.h.

◆ _Jacobian_mult

const MaterialProperty<RankFourTensor>& CrackTipEnrichmentStressDivergenceTensors::_Jacobian_mult

protected

Definition at line 46 of file CrackTipEnrichmentStressDivergenceTensors.h.

Referenced by computeQpJacobian(), and computeQpOffDiagJacobian().

◆ _ndisp

unsigned int CrackTipEnrichmentStressDivergenceTensors::_ndisp

protected

Coupled displacement variables.

Definition at line 62 of file CrackTipEnrichmentStressDivergenceTensors.h.

Referenced by CrackTipEnrichmentStressDivergenceTensors().

◆ _nenrich_disp

unsigned int CrackTipEnrichmentStressDivergenceTensors::_nenrich_disp

protected

Coupled enrichment displacement variables.

Definition at line 58 of file CrackTipEnrichmentStressDivergenceTensors.h.

Referenced by CrackTipEnrichmentStressDivergenceTensors().

◆ _rotation_increment

const MaterialProperty<RankTwoTensor>* CrackTipEnrichmentStressDivergenceTensors::_rotation_increment

protected

Definition at line 50 of file CrackTipEnrichmentStressDivergenceTensors.h.

◆ _stress

const MaterialProperty<RankTwoTensor>& CrackTipEnrichmentStressDivergenceTensors::_stress

protected

Definition at line 45 of file CrackTipEnrichmentStressDivergenceTensors.h.

Referenced by computeQpResidual().

◆ _var_undisplaced

MooseVariable& ALEKernel::_var_undisplaced

protectedinherited

Reference to this Kernel's undisplaced MooseVariable object.

Definition at line 36 of file ALEKernel.h.

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

xfem/include/kernels/CrackTipEnrichmentStressDivergenceTensors.h
xfem/src/kernels/CrackTipEnrichmentStressDivergenceTensors.C

Public Member Functions

Protected Member Functions

Protected Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ CrackTipEnrichmentStressDivergenceTensors()

Member Function Documentation

◆ computeJacobian()

◆ computeOffDiagJacobian()

◆ computeQpJacobian()

◆ computeQpOffDiagJacobian()

◆ computeQpResidual()

◆ crackTipEnrichementFunctionAtPoint()

◆ crackTipEnrichementFunctionDerivativeAtPoint()

◆ rotateFromCrackFrontCoordsToGlobal()

Member Data Documentation

◆ _assembly_undisplaced

◆ _B

◆ _base_name

◆ _BI

◆ _BJ

◆ _component

◆ _dBX

◆ _dBx

◆ _deformation_gradient

◆ _deformation_gradient_old

◆ _disp_var

◆ _enrich_disp_var

◆ _enrichment_component

◆ _grad_phi_undisplaced

◆ _grad_test_undisplaced

◆ _Jacobian_mult

◆ _ndisp

◆ _nenrich_disp

◆ _rotation_increment

◆ _stress

◆ _var_undisplaced