EnrichmentFunctionCalculation Class Reference

Perform calculation of enrichment function values and derivatives. More...

#include <EnrichmentFunctionCalculation.h>

Inheritance diagram for EnrichmentFunctionCalculation:

Public Member Functions
	EnrichmentFunctionCalculation (const CrackFrontDefinition *crack_front_definition)
virtual unsigned int	crackTipEnrichementFunctionAtPoint (const Point &point, std::vector< Real > &B)
	calculate the enrichment function values at point
virtual unsigned int	crackTipEnrichementFunctionDerivativeAtPoint (const Point &point, std::vector< RealVectorValue > &dB)
	calculate the enrichment function derivatives at point
void	rotateFromCrackFrontCoordsToGlobal (const RealVectorValue &vector, RealVectorValue &rotated_vector, const unsigned int point_index)
	rotate a vector from crack front coordinate to global cooridate

Private Attributes
const CrackFrontDefinition &	_crack_front_definition
Real	_r
Real	_theta

Detailed Description

Perform calculation of enrichment function values and derivatives.

Definition at line 17 of file EnrichmentFunctionCalculation.h.

Constructor & Destructor Documentation

EnrichmentFunctionCalculation()

EnrichmentFunctionCalculation::EnrichmentFunctionCalculation

(

const CrackFrontDefinition *

crack_front_definition

)

Definition at line 12 of file EnrichmentFunctionCalculation.C.

  : _crack_front_definition(*crack_front_definition)
{
}

Member Function Documentation

crackTipEnrichementFunctionAtPoint()

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

virtual

calculate the enrichment function values at point

Returns

the closest crack front index

Definition at line 19 of file EnrichmentFunctionCalculation.C.

{
  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;
}

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

crackTipEnrichementFunctionDerivativeAtPoint()

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

virtual

calculate the enrichment function derivatives at point

Returns

the closest crack front index

Definition at line 43 of file EnrichmentFunctionCalculation.C.

{
  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;
}

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

rotateFromCrackFrontCoordsToGlobal()

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

rotate a vector from crack front coordinate to global cooridate

Parameters

rotated_vector

rotated vector

Definition at line 78 of file EnrichmentFunctionCalculation.C.

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

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

Member Data Documentation

_crack_front_definition

const CrackFrontDefinition& EnrichmentFunctionCalculation::_crack_front_definition

private

Definition at line 43 of file EnrichmentFunctionCalculation.h.

Referenced by crackTipEnrichementFunctionAtPoint(), crackTipEnrichementFunctionDerivativeAtPoint(), and rotateFromCrackFrontCoordsToGlobal().

_r

Real EnrichmentFunctionCalculation::_r

private

Definition at line 44 of file EnrichmentFunctionCalculation.h.

Referenced by crackTipEnrichementFunctionAtPoint(), and crackTipEnrichementFunctionDerivativeAtPoint().

_theta

Real EnrichmentFunctionCalculation::_theta

private

Definition at line 45 of file EnrichmentFunctionCalculation.h.

Referenced by crackTipEnrichementFunctionAtPoint(), and crackTipEnrichementFunctionDerivativeAtPoint().

---

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

• xfem/include/base/EnrichmentFunctionCalculation.h
• xfem/src/base/EnrichmentFunctionCalculation.C