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