#include <users_modules/fracture_mechanics/src/GriffithForceElement.hpp>

Inheritance diagram for FractureMechanics::GriffithForceElement::AuxFunctions< TYPE >:

Collaboration diagram for FractureMechanics::GriffithForceElement::AuxFunctions< TYPE >:

Public Member Functions
	AuxFunctions ()

MoFEMErrorCode	calculateAreaAndNormal (const MatrixDouble &diff_n)

MoFEMErrorCode	calculateMatA (const MatrixDouble &diff_n)

MoFEMErrorCode	calculateGriffithForce (const double gc, const double beta, const MatrixDouble &diff_n)

Public Attributes
ublas::vector< TYPE >	referenceCoords

ublas::vector< TYPE >	currentCoords

ublas::vector< TYPE >	griffithForce

FTensor::Tensor1< TYPE, 3 >	t_Normal

FTensor::Tensor1< TYPE, 3 >	t_Tangent1_ksi

FTensor::Tensor1< TYPE, 3 >	t_Tangent2_eta

TYPE	currentArea

FTensor::Index< 'i', 3 >	i

FTensor::Index< 'j', 3 >	j

FTensor::Index< 'k', 3 >	k

FTensor::Number< 0 >	N0

FTensor::Number< 1 >	N1

boost::function< FTensor::Tensor1< TYPE *, 3 >(ublas::vector< TYPE > &v)>	get_ftensor_from_vec

boost::function< FTensor::Tensor1< FTensor::PackPtr< const double *, 2 >, 2 >(const MatrixDouble &m)>	get_ftensor_from_mat_2d

Detailed Description

template<typename TYPE>
struct FractureMechanics::GriffithForceElement::AuxFunctions< TYPE >

\[ A = \| \mathbf{N} \| = \left\| \textrm{Spin}\left[ \frac{\partial\mathbf{N}\mathbf{X}}{\partial\xi} \right] \frac{\partial\mathbf{N}\mathbf{X}}{\partial\eta} \right\| \]

Definition at line 123 of file GriffithForceElement.hpp.

Constructor & Destructor Documentation

◆ AuxFunctions()

template<typename TYPE >

FractureMechanics::GriffithForceElement::AuxFunctions< TYPE >::AuxFunctions ( )

inline

Definition at line 147 of file GriffithForceElement.hpp.

                   {
 
      get_ftensor_from_vec = [](ublas::vector<TYPE> &v) {
        return FTensor::Tensor1<TYPE *, 3>(&v(0), &v(1), &v(2), 3);
      };
 
      get_ftensor_from_mat_2d = [](const MatrixDouble &m) {
        return FTensor::Tensor1<FTensor::PackPtr<const double *, 2>, 2>(
            &m(0, 0), &m(0, 1));
      };
    }

Member Function Documentation

◆ calculateAreaAndNormal()

template<typename TYPE >

MoFEMErrorCode FractureMechanics::GriffithForceElement::AuxFunctions< TYPE >::calculateAreaAndNormal ( const MatrixDouble & diff_n )

inline

Definition at line 159 of file GriffithForceElement.hpp.

                                                                      {
      MoFEMFunctionBeginHot;
 
      auto t_coords = get_ftensor_from_vec(currentCoords);
      auto t_diff = get_ftensor_from_mat_2d(diff_n);
      t_Tangent1_ksi(i) = (TYPE)0;
      t_Tangent2_eta(i) = (TYPE)0;
      for (int nn = 0; nn != 3; ++nn) {
        t_Tangent1_ksi(i) += t_coords(i) * t_diff(N0);
        t_Tangent2_eta(i) += t_coords(i) * t_diff(N1);
        ++t_diff;
        ++t_coords;
      }
 
      t_Normal(j) =
          FTensor::levi_civita(i, j, k) * t_Tangent1_ksi(k) * t_Tangent2_eta(i);
      currentArea = 0.5 * sqrt(t_Normal(i) * t_Normal(i));
 
      MoFEMFunctionReturnHot(0);
    }

◆ calculateGriffithForce()

template<typename TYPE >

MoFEMErrorCode FractureMechanics::GriffithForceElement::AuxFunctions< TYPE >::calculateGriffithForce	(	const double	gc,
		const double	beta,
		const MatrixDouble &	diff_n
	)

inline

Definition at line 201 of file GriffithForceElement.hpp.

                                                                      {
      MoFEMFunctionBeginHot;
 
      CHKERR calculateAreaAndNormal(diff_n);
      CHKERR calculateMatA(diff_n);
 
      for (int dd = 0; dd != 9; dd++)
        griffithForce[dd] *= gc * beta;
 
      MoFEMFunctionReturnHot(0);
    }

◆ calculateMatA()

template<typename TYPE >

MoFEMErrorCode FractureMechanics::GriffithForceElement::AuxFunctions< TYPE >::calculateMatA ( const MatrixDouble & diff_n )

inline

Definition at line 180 of file GriffithForceElement.hpp.

                                                             {
      MoFEMFunctionBeginHot;
 
      griffithForce.resize(9, false);
      std::fill(griffithForce.begin(), griffithForce.end(), 0);
      auto t_griffith = get_ftensor_from_vec(griffithForce);
      auto t_diff = get_ftensor_from_mat_2d(diff_n);
 
      for (int dd = 0; dd != 3; dd++) {
 
        t_griffith(i) += FTensor::levi_civita(i, j, k) * t_Tangent1_ksi(k) *
                         t_diff(N1) * t_Normal(j) * 0.25 / currentArea;
        t_griffith(i) -= FTensor::levi_civita(i, j, k) * t_Tangent2_eta(k) *
                         t_diff(N0) * t_Normal(j) * 0.25 / currentArea;
        ++t_diff;
        ++t_griffith;
      }
 
      MoFEMFunctionReturnHot(0);
    }