CohesiveZoneModelTools Namespace Reference

Functions
template<bool is_ad = false>
GenericRealVectorValue< is_ad >	computedVnormdV (const GenericRealVectorValue< is_ad > &V)
template<bool is_ad = false>
GenericRankFourTensor< is_ad >	computedRdF (const GenericRankTwoTensor< is_ad > &R, const GenericRankTwoTensor< is_ad > &U)
	compute the derivative of the rotation matrix, R=FU^-1, using Chen and Wheeler 1993
RankFourTensor	computedFinversedF (const RankTwoTensor &F_inv)
	compute the derivative of F^-1 w.r.t. to F
template<bool is_ad = false>
GenericRealVectorValue< is_ad >	computeNormalComponents (const GenericRealVectorValue< is_ad > &normal, const GenericRealVectorValue< is_ad > &vector)
	compute the normal componets of a vector
template<bool is_ad = false>
GenericRealVectorValue< is_ad >	computeTangentComponents (const GenericRealVectorValue< is_ad > &normal, const GenericRealVectorValue< is_ad > &vector)
	compute the tangent componets of a vector
template<bool is_ad = false>
GenericRankTwoTensor< is_ad >	computeReferenceRotation (const GenericRealVectorValue< is_ad > &normal, const unsigned int mesh_dimension)
	compute the czm reference rotation based on the normal in the undeformed configuration and the mesh dimension
template<bool is_ad = false>
GenericReal< is_ad >	computeAreaRatio (const RankTwoTensor &FinvT, const Real &J, const RealVectorValue &N)
	compute the area ratio betweeen the deformed and undeformed configuration, and its derivatives w.r.t.
template<bool is_ad = false>
GenericRankTwoTensor< is_ad >	computeDAreaRatioDF (const GenericRankTwoTensor< is_ad > &FinvT, const GenericRealVectorValue< is_ad > &N, const GenericReal< is_ad > &J, const GenericRankFourTensor< is_ad > &DFinv_DF)

Function Documentation

computeAreaRatio()

template<bool is_ad = false>

GenericReal< is_ad > CohesiveZoneModelTools::computeAreaRatio	(	const RankTwoTensor &	FinvT,
		const Real &	J,
		const RealVectorValue &	N
	)

compute the area ratio betweeen the deformed and undeformed configuration, and its derivatives w.r.t.

the deformation gradient, F

Definition at line 61 of file CohesiveZoneModelTools.h.

{
  return J * (FinvT * N).norm();
}

Referenced by CZMComputeGlobalTractionTotalLagrangian::computeEquilibriumTracionAndDerivatives().

computeDAreaRatioDF()

template<bool is_ad = false>

GenericRankTwoTensor< is_ad > CohesiveZoneModelTools::computeDAreaRatioDF	(	const GenericRankTwoTensor< is_ad > &	FinvT,
		const GenericRealVectorValue< is_ad > &	N,
		const GenericReal< is_ad > &	J,
		const GenericRankFourTensor< is_ad > &	DFinv_DF
	)

Definition at line 68 of file CohesiveZoneModelTools.h.

{
  const auto Fitr_N = FinvT * N;
  const auto Fitr_N_norm = Fitr_N.norm();
  GenericRankTwoTensor<is_ad> R2temp;
  for (unsigned int l = 0; l < 3; l++)
    for (unsigned int m = 0; m < 3; m++)
    {
      R2temp(l, m) = 0;
      for (unsigned int i = 0; i < 3; i++)
        for (unsigned int j = 0; j < 3; j++)
          R2temp(l, m) += Fitr_N(i) * DFinv_DF(j, i, l, m) * N(j);
 
      R2temp(l, m) *= J / Fitr_N_norm;
    }
 
  return J * FinvT * Fitr_N_norm + R2temp;
}

Referenced by CZMComputeGlobalTractionTotalLagrangian::computeAreaRatioAndDerivatives().

computedFinversedF()

RankFourTensor CohesiveZoneModelTools::computedFinversedF

(

const RankTwoTensor &

F_inv

)

compute the derivative of F^-1 w.r.t. to F

Definition at line 16 of file CohesiveZoneModelTools.C.

{
  usingTensorIndices(i_, j_, k_, l_);
  return -F_inv.times<i_, k_, j_, l_>(F_inv.transpose());
}

Referenced by CZMComputeGlobalTractionTotalLagrangian::computeAreaRatioAndDerivatives().

computedRdF()

template<bool is_ad = false>

GenericRankFourTensor< is_ad > CohesiveZoneModelTools::computedRdF	(	const GenericRankTwoTensor< is_ad > &	R,
		const GenericRankTwoTensor< is_ad > &	U
	)

compute the derivative of the rotation matrix, R=FU^-1, using Chen and Wheeler 1993

Definition at line 30 of file CohesiveZoneModelTools.h.

{
  const auto Uhat = U.trace() * RankTwoTensor::Identity() - U;
  unsigned int k, l, m, n, p, q;
  const auto Uhat_det = Uhat.det();
 
  GenericRankFourTensor<is_ad> dR_dF;
  for (k = 0; k < 3; k++)
    for (l = 0; l < 3; l++)
      for (m = 0; m < 3; m++)
        for (n = 0; n < 3; n++)
        {
          dR_dF(k, l, m, n) = 0.;
          for (p = 0; p < 3; p++)
            for (q = 0; q < 3; q++)
              dR_dF(k, l, m, n) +=
                  R(k, p) * (Uhat(p, q) * R(m, q) * Uhat(n, l) - Uhat(p, n) * R(m, q) * Uhat(q, l));
 
          dR_dF(k, l, m, n) /= Uhat_det;
        }
 
  return dR_dF;
}

Referenced by CZMComputeGlobalTractionTotalLagrangian::computeAreaRatioAndDerivatives().

computedVnormdV()

template<bool is_ad = false>

GenericRealVectorValue< is_ad > CohesiveZoneModelTools::computedVnormdV

(

const GenericRealVectorValue< is_ad > &

V

)

Definition at line 22 of file CohesiveZoneModelTools.h.

{
  return V / V.norm();
}

computeNormalComponents()

template<bool is_ad = false>

GenericRealVectorValue< is_ad > CohesiveZoneModelTools::computeNormalComponents	(	const GenericRealVectorValue< is_ad > &	normal,
		const GenericRealVectorValue< is_ad > &	vector
	)

compute the normal componets of a vector

Definition at line 94 of file CohesiveZoneModelTools.h.

{
  return (normal * vector) * normal;
}

Referenced by CZMRealVectorScalar::computeQpProperties().

computeReferenceRotation()

template<bool is_ad = false>

GenericRankTwoTensor< is_ad > CohesiveZoneModelTools::computeReferenceRotation	(	const GenericRealVectorValue< is_ad > &	normal,
		const unsigned int	mesh_dimension
	)

compute the czm reference rotation based on the normal in the undeformed configuration and the mesh dimension

Definition at line 112 of file CohesiveZoneModelTools.h.

{
  GenericRankTwoTensor<is_ad> rot;
  switch (mesh_dimension)
  {
    case 3:
      rot = RotationMatrix::rotVec1ToVec2<is_ad>(GenericRealVectorValue<false>(1, 0, 0), normal);
      break;
    case 2:
      rot = RotationMatrix::rotVec2DToX<is_ad>(normal).transpose();
      break;
    case 1:
      rot = RankTwoTensor::Identity();
      break;
    default:
      mooseError("computeReferenceRotation: mesh_dimension value should be 1, 2 or, 3. You "
                 "provided " +
                 std::to_string(mesh_dimension));
  }
  return rot;
}

computeTangentComponents()

template<bool is_ad = false>

GenericRealVectorValue< is_ad > CohesiveZoneModelTools::computeTangentComponents	(	const GenericRealVectorValue< is_ad > &	normal,
		const GenericRealVectorValue< is_ad > &	vector
	)

compute the tangent componets of a vector

Definition at line 103 of file CohesiveZoneModelTools.h.

{
  return vector - computeNormalComponents<is_ad>(normal, vector);
}

Referenced by CZMRealVectorScalar::computeQpProperties().