Volume Length Quality. More...

#include "tools/src/VolumeLengthQuality.hpp"

Inheritance diagram for VolumeLengthQuality< TYPE >:

Collaboration diagram for VolumeLengthQuality< TYPE >:

Public Member Functions
	VolumeLengthQuality (VolumeLengthQualityType type, double alpha, double gamma, const std::string &field_name="MESH_NODE_POSITIONS")

MoFEMErrorCode	getMaterialOptions ()

MoFEMErrorCode	getEdgesFromElemCoords ()

MoFEMErrorCode	calculateLrms ()
	Calculate mean element edge length.

MoFEMErrorCode	calculateQ ()
	Calculate Q.

virtual MoFEMErrorCode	calculateP_PiolaKirchhoffI ()
	Volume Length Quality.

	VolumeLengthQuality (VolumeLengthQualityType type, double alpha, double gamma)

MoFEMErrorCode	getMaterialOptions ()

MoFEMErrorCode	getEdgesFromElemCoords ()

MoFEMErrorCode	calculateLrms ()
	Calculate mean element edge length.

MoFEMErrorCode	calculateQ ()
	Calculate Q.

virtual MoFEMErrorCode	calculateP_PiolaKirchhoffI (const NonlinearElasticElement::BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	Volume Length Quality.

	VolumeLengthQuality (VolumeLengthQualityType type, double alpha, double gamma)

MoFEMErrorCode	getMaterialOptions ()

MoFEMErrorCode	getEdgesFromElemCoords ()

MoFEMErrorCode	calculateLrms ()
	Calculate mean element edge length.

MoFEMErrorCode	calculateQ ()
	Calculate Q.

virtual MoFEMErrorCode	calculateP_PiolaKirchhoffI (const NonlinearElasticElement::BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	Volume Length Quality.

Public Member Functions inherited from Smoother::CalculatePK1< TYPE >
	CalculatePK1 (const std::string field_name)

Public Member Functions inherited from MoFEM::VolumeElementForcesAndSourcesCore::UserDataOperator
int	getNumNodes ()
	get element number of nodes

const EntityHandle *	getConn ()
	get element connectivity

double	getVolume () const
	element volume (linear geometry)

double &	getVolume ()
	element volume (linear geometry)

FTensor::Tensor2< double *, 3, 3 > &	getJac ()
	get element Jacobian

FTensor::Tensor2< double *, 3, 3 > &	getInvJac ()
	get element inverse Jacobian

VectorDouble &	getCoords ()
	nodal coordinates

VolumeElementForcesAndSourcesCore *	getVolumeFE () const
	return pointer to Generic Volume Finite Element object

Public Member Functions inherited from MoFEM::ForcesAndSourcesCore::UserDataOperator
	UserDataOperator (const FieldSpace space, const char type=OPSPACE, const bool symm=true)
	Constructor for operators working on finite element spaces.

	UserDataOperator (const std::string field_name, const char type, const bool symm=true)
	Constructor for operators working on a single field.

	UserDataOperator (const std::string row_field_name, const std::string col_field_name, const char type, const bool symm=true)
	Constructor for operators working on two fields (bilinear forms)

boost::shared_ptr< const NumeredEntFiniteElement >	getNumeredEntFiniteElementPtr () const
	Return raw pointer to NumeredEntFiniteElement.

EntityHandle	getFEEntityHandle () const
	Return finite element entity handle.

int	getFEDim () const
	Get dimension of finite element.

EntityType	getFEType () const
	Get dimension of finite element.

boost::weak_ptr< SideNumber >	getSideNumberPtr (const int side_number, const EntityType type)
	Get the side number pointer.

EntityHandle	getSideEntity (const int side_number, const EntityType type)
	Get the side entity.

int	getNumberOfNodesOnElement () const
	Get the number of nodes on finite element.

MoFEMErrorCode	getProblemRowIndices (const std::string filed_name, const EntityType type, const int side, VectorInt &indices) const
	Get row indices.

MoFEMErrorCode	getProblemColIndices (const std::string filed_name, const EntityType type, const int side, VectorInt &indices) const
	Get col indices.

const FEMethod *	getFEMethod () const
	Return raw pointer to Finite Element Method object.

int	getOpType () const
	Get operator types.

void	setOpType (const OpType type)
	Set operator type.

void	addOpType (const OpType type)
	Add operator type.

int	getNinTheLoop () const
	get number of finite element in the loop

int	getLoopSize () const
	get size of elements in the loop

std::string	getFEName () const
	Get name of the element.

ForcesAndSourcesCore *	getPtrFE () const

ForcesAndSourcesCore *	getSidePtrFE () const

ForcesAndSourcesCore *	getRefinePtrFE () const

const PetscData::Switches &	getDataCtx () const

const KspMethod::KSPContext	getKSPCtx () const

const SnesMethod::SNESContext	getSNESCtx () const

const TSMethod::TSContext	getTSCtx () const

Vec	getKSPf () const

Mat	getKSPA () const

Mat	getKSPB () const

Vec	getSNESf () const

Vec	getSNESx () const

Mat	getSNESA () const

Mat	getSNESB () const

Vec	getTSu () const

Vec	getTSu_t () const

Vec	getTSu_tt () const

Vec	getTSf () const

Mat	getTSA () const

Mat	getTSB () const

int	getTSstep () const

double	getTStime () const

double	getTStimeStep () const

double	getTSa () const

double	getTSaa () const

MatrixDouble &	getGaussPts ()
	matrix of integration (Gauss) points for Volume Element

auto	getFTensor0IntegrationWeight ()
	Get integration weights.

MatrixDouble &	getCoordsAtGaussPts ()
	Gauss points and weight, matrix (nb. of points x 3)

auto	getFTensor1CoordsAtGaussPts ()
	Get coordinates at integration points assuming linear geometry.

double	getMeasure () const
	get measure of element

double &	getMeasure ()
	get measure of element

MoFEM::Interface &	getMField ()

moab::Interface &	getMoab ()

virtual boost::weak_ptr< ForcesAndSourcesCore >	getSubPipelinePtr () const

MoFEMErrorCode	loopSide (const string &fe_name, ForcesAndSourcesCore side_fe, const size_t dim, const EntityHandle ent_for_side=0, boost::shared_ptr< Range > fe_range=nullptr, const int verb=QUIET, const LogManager::SeverityLevel sev=Sev::noisy, AdjCache adj_cache=nullptr)
	User calls this function to loop over elements on the side of face. This function calls finite element with its operator to do calculations.

MoFEMErrorCode	loopThis (const string &fe_name, ForcesAndSourcesCore *this_fe, const int verb=QUIET, const LogManager::SeverityLevel sev=Sev::noisy)
	User calls this function to loop over the same element using a different set of integration points. This function calls finite element with its operator to do calculations.

MoFEMErrorCode	loopParent (const string &fe_name, ForcesAndSourcesCore *parent_fe, const int verb=QUIET, const LogManager::SeverityLevel sev=Sev::noisy)
	User calls this function to loop over parent elements. This function calls finite element with its operator to do calculations.

MoFEMErrorCode	loopChildren (const string &fe_name, ForcesAndSourcesCore *child_fe, const int verb=QUIET, const LogManager::SeverityLevel sev=Sev::noisy)
	User calls this function to loop over parent elements. This function calls finite element with its operator to do calculations.

MoFEMErrorCode	loopRange (const string &fe_name, ForcesAndSourcesCore *range_fe, boost::shared_ptr< Range > fe_range, const int verb=QUIET, const LogManager::SeverityLevel sev=Sev::noisy)
	Iterate over range of elements.

Public Member Functions inherited from MoFEM::DataOperator
	DataOperator (const bool symm=true)

virtual	~DataOperator ()=default

virtual MoFEMErrorCode	doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)
	Operator for bi-linear form, usually to calculate values on left hand side.

virtual MoFEMErrorCode	opLhs (EntitiesFieldData &row_data, EntitiesFieldData &col_data)

virtual MoFEMErrorCode	doWork (int side, EntityType type, EntitiesFieldData::EntData &data)
	Operator for linear form, usually to calculate values on right hand side.

virtual MoFEMErrorCode	opRhs (EntitiesFieldData &data, const bool error_if_no_base=false)

bool	getSymm () const
	Get if operator uses symmetry of DOFs or not.

void	setSymm ()
	set if operator is executed taking in account symmetry

void	unSetSymm ()
	unset if operator is executed for non symmetric problem

Public Member Functions inherited from NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI< TYPE >
	FunctionsToCalculatePiolaKirchhoffI ()

virtual	~FunctionsToCalculatePiolaKirchhoffI ()=default

MoFEMErrorCode	calculateC_CauchyDeformationTensor ()

MoFEMErrorCode	calculateE_GreenStrain ()

MoFEMErrorCode	calculateS_PiolaKirchhoffII ()

virtual MoFEMErrorCode	calculateCauchyStress (const BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	Function overload to implement user material.

virtual MoFEMErrorCode	setUserActiveVariables (int &nb_active_variables)
	add additional active variables

virtual MoFEMErrorCode	setUserActiveVariables (VectorDouble &activeVariables)
	Add additional independent variables More complex physical models depend on gradient of defamation and some additional variables. For example can depend on temperature. This function adds additional independent variables to the model.

virtual MoFEMErrorCode	calculateElasticEnergy (const BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	Calculate elastic energy density.

virtual MoFEMErrorCode	calculatesIGma_EshelbyStress (const BlockData block_data, boost::shared_ptr< const NumeredEntFiniteElement > fe_ptr)
	Calculate Eshelby stress.

virtual MoFEMErrorCode	getDataOnPostProcessor (std::map< std::string, std::vector< VectorDouble > > &field_map, std::map< std::string, std::vector< MatrixDouble > > &grad_map)
	Do operations when pre-process.

Public Attributes
int	tYpe

double	aLpha

double	gAmma

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

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

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

VectorDouble	coordsEdges

double	lrmsSquared0

VectorDouble	deltaChi

ublas::vector< TYPE >	deltaX

ublas::matrix< TYPE >	Q

ublas::matrix< TYPE >	dXdChiT

TYPE	lrmsSquared

TYPE	q

TYPE	b

TYPE	detF

TYPE	currentVolume

TYPE	tMp

Public Attributes inherited from Smoother::CalculatePK1< TYPE >
FTensor::Index< 'i', 3 >	i

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

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

MatrixBoundedArray< TYPE, 9 >	F

MatrixBoundedArray< TYPE, 9 >	invF

MatrixBoundedArray< TYPE, 9 >	P

FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_F

FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_invF

FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_P

int	gG = 0
	Gauss point number.

CommonData *	commonDataPtr = nullptr
	field values at Gauss pts.)

MoFEM::VolumeElementForcesAndSourcesCore::UserDataOperator *	opPtr = nullptr
	pointer to finite element tetrahedral operator

Public Attributes inherited from MoFEM::ForcesAndSourcesCore::UserDataOperator
char	opType

std::string	rowFieldName

std::string	colFieldName

FieldSpace	sPace

Public Attributes inherited from MoFEM::DataOperator
DoWorkLhsHookFunType	doWorkLhsHook

DoWorkRhsHookFunType	doWorkRhsHook

bool	sYmm
	If true assume that matrix is symmetric structure.

std::array< bool, MBMAXTYPE >	doEntities
	If true operator is executed for entity.

bool &	doVertices
	\deprectaed If false skip vertices

bool &	doEdges
	\deprectaed If false skip edges

bool &	doQuads
	\deprectaed

bool &	doTris
	\deprectaed

bool &	doTets
	\deprectaed

bool &	doPrisms
	\deprectaed

Public Attributes inherited from NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI< TYPE >
FTensor::Index< 'i', 3 >	i

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

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

double	lambda

double	mu

MatrixBoundedArray< TYPE, 9 >	F

MatrixBoundedArray< TYPE, 9 >	C

MatrixBoundedArray< TYPE, 9 >	E

MatrixBoundedArray< TYPE, 9 >	S

MatrixBoundedArray< TYPE, 9 >	invF

MatrixBoundedArray< TYPE, 9 >	P

MatrixBoundedArray< TYPE, 9 >	sIGma

MatrixBoundedArray< TYPE, 9 >	h

MatrixBoundedArray< TYPE, 9 >	H

MatrixBoundedArray< TYPE, 9 >	invH

MatrixBoundedArray< TYPE, 9 >	sigmaCauchy

FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_F

FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_C

FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_E

FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_S

FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_invF

FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_P

FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_sIGma

FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_h

FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_H

FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_invH

FTensor::Tensor2< FTensor::PackPtr< TYPE *, 0 >, 3, 3 >	t_sigmaCauchy

TYPE	J

TYPE	eNergy

TYPE	detH

TYPE	detF

int	gG
	Gauss point number.

CommonData *	commonDataPtr

MoFEM::VolumeElementForcesAndSourcesCore::UserDataOperator *	opPtr
	pointer to finite element tetrahedral operator

Additional Inherited Members
Public Types inherited from MoFEM::ForcesAndSourcesCore::UserDataOperator
enum	OpType { OPROW = 1 << 0 , OPCOL = 1 << 1 , OPROWCOL = 1 << 2 , OPSPACE = 1 << 3 , OPLAST = 1 << 3 }
	Controls loop over entities on element. More...

using	AdjCache = std::map< EntityHandle, std::vector< boost::weak_ptr< NumeredEntFiniteElement > > >

Public Types inherited from MoFEM::DataOperator
using	DoWorkLhsHookFunType = boost::function< MoFEMErrorCode(DataOperator *op_ptr, int row_side, int col_side, EntityType row_type, EntityType col_type, EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)>

using	DoWorkRhsHookFunType = boost::function< MoFEMErrorCode(DataOperator *op_ptr, int side, EntityType type, EntitiesFieldData::EntData &data)>

Static Public Attributes inherited from MoFEM::ForcesAndSourcesCore::UserDataOperator
static const char *const	OpTypeNames []

Protected Member Functions inherited from MoFEM::VolumeElementForcesAndSourcesCore::UserDataOperator
MoFEMErrorCode	setPtrFE (ForcesAndSourcesCore *ptr)

Static Protected Member Functions inherited from Smoother::CalculatePK1< TYPE >
static auto	resizeAndSet (MatrixBoundedArray< TYPE, 9 > &m)

Static Protected Member Functions inherited from NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI< TYPE >
static auto	resizeAndSet (MatrixBoundedArray< TYPE, 9 > &m)

Protected Attributes inherited from MoFEM::ForcesAndSourcesCore::UserDataOperator
ForcesAndSourcesCore *	ptrFE

Detailed Description

template<typename TYPE>
struct VolumeLengthQuality< TYPE >

Volume Length Quality.

Definition at line 34 of file VolumeLengthQuality.hpp.

Constructor & Destructor Documentation

◆ VolumeLengthQuality() [1/3]

template<typename TYPE >

VolumeLengthQuality< TYPE >::VolumeLengthQuality	(	VolumeLengthQualityType	type,
		double	alpha,
		double	gamma,
		const std::string &	field_name = `"MESH_NODE_POSITIONS"`
	)

inline

Definition at line 46 of file VolumeLengthQuality.hpp.

      : Smoother::CalculatePK1<TYPE>(field_name), tYpe(type), aLpha(alpha),
        gAmma(gamma) {
    ierr = getMaterialOptions();
    CHKERRABORT(PETSC_COMM_WORLD, ierr);
  }

◆ VolumeLengthQuality() [2/3]

template<typename TYPE >

VolumeLengthQuality< TYPE >::VolumeLengthQuality	(	VolumeLengthQualityType	type,
		double	alpha,
		double	gamma
	)

inline

Definition at line 45 of file VolumeLengthQuality.hpp.

      : NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI<TYPE>(),
        tYpe(type), aLpha(alpha), gAmma(gamma) {
    ierr = getMaterialOptions();
    CHKERRABORT(PETSC_COMM_WORLD, ierr);
  }

◆ VolumeLengthQuality() [3/3]

template<typename TYPE >

VolumeLengthQuality< TYPE >::VolumeLengthQuality	(	VolumeLengthQualityType	type,
		double	alpha,
		double	gamma
	)

inline

Definition at line 45 of file VolumeLengthQuality.hpp.

      : NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI<TYPE>(),
        tYpe(type), aLpha(alpha), gAmma(gamma) {
    ierr = getMaterialOptions();
    CHKERRABORT(PETSC_COMM_WORLD, ierr);
  }

Member Function Documentation

◆ calculateLrms() [1/3]

template<typename TYPE >

MoFEMErrorCode VolumeLengthQuality< TYPE >::calculateLrms ( )

inline

Calculate mean element edge length.

\[ \Delta \boldsymbol\chi = \boldsymbol\chi^1 - \boldsymbol\chi^2 \]

\[ \Delta X = \mathbf{F} \Delta \boldsymbol\chi \]

\[ l_\textrm{rms} = \sqrt{\frac{1}{6} \sum_{i=0}^6 l_i^2 } = L_\textrm{rms}dl_\textrm{rms} \]

Definition at line 130 of file VolumeLengthQuality.hpp.

                                 {
    MoFEMFunctionBegin;
    if (deltaChi.size() != 3) {
      deltaChi.resize(3);
      deltaX.resize(3);
      dXdChiT.resize(3, 3);
    }
    lrmsSquared = 0;
    lrmsSquared0 = 0;
    std::fill(dXdChiT.data().begin(), dXdChiT.data().end(), 0);
    for (int ee = 0; ee < 6; ee++) {
      for (int dd = 0; dd < 3; dd++) {
        deltaChi[dd] = coordsEdges[6 * ee + dd] - coordsEdges[6 * ee + 3 + dd];
      }
      std::fill(deltaX.begin(), deltaX.end(), 0);
      for (int ii = 0; ii != 3; ++ii)
        for (int jj = 0; jj != 3; ++jj)
          deltaX(ii) += this->F(ii, jj) * deltaChi(jj);
 
      for (int dd = 0; dd < 3; dd++) {
        lrmsSquared += (1. / 6.) * deltaX[dd] * deltaX[dd];
        lrmsSquared0 += (1. / 6.) * deltaChi[dd] * deltaChi[dd];
      }
      for (int ii = 0; ii != 3; ++ii)
        for (int jj = 0; jj != 3; ++jj)
          dXdChiT(ii, jj) += deltaX(ii) * deltaChi(jj);
    }
    MoFEMFunctionReturn(0);
  }

◆ calculateLrms() [2/3]

template<typename TYPE >

MoFEMErrorCode VolumeLengthQuality< TYPE >::calculateLrms ( )

inline

Calculate mean element edge length.

\[ \Delta \boldsymbol\chi = \boldsymbol\chi^1 - \boldsymbol\chi^2 \]

\[ \Delta X = \mathbf{F} \Delta \boldsymbol\chi \]

\[ l_\textrm{rms} = \sqrt{\frac{1}{6} \sum_{i=0}^6 l_i^2 } = L_\textrm{rms}dl_\textrm{rms} \]

Definition at line 127 of file VolumeLengthQuality.hpp.

                                 {
    MoFEMFunctionBegin;
    if (deltaChi.size() != 3) {
      deltaChi.resize(3);
      deltaX.resize(3);
      dXdChiT.resize(3, 3);
    }
    lrmsSquared = 0;
    lrmsSquared0 = 0;
    std::fill(dXdChiT.data().begin(), dXdChiT.data().end(), 0);
    for (int ee = 0; ee < 6; ee++) {
      for (int dd = 0; dd < 3; dd++) {
        deltaChi[dd] = coordsEdges[6 * ee + dd] - coordsEdges[6 * ee + 3 + dd];
      }
      std::fill(deltaX.begin(), deltaX.end(), 0);
      for (int ii = 0; ii != 3; ++ii)
        for (int jj = 0; jj != 3; ++jj)
          deltaX(ii) += this->F(ii, jj) * deltaChi(jj);
 
      for (int dd = 0; dd < 3; dd++) {
        lrmsSquared += (1. / 6.) * deltaX[dd] * deltaX[dd];
        lrmsSquared0 += (1. / 6.) * deltaChi[dd] * deltaChi[dd];
      }
      for (int ii = 0; ii != 3; ++ii)
        for (int jj = 0; jj != 3; ++jj)
          dXdChiT(ii, jj) += deltaX(ii) * deltaChi(jj);
    }
    MoFEMFunctionReturn(0);
  }

◆ calculateLrms() [3/3]

template<typename TYPE >

MoFEMErrorCode VolumeLengthQuality< TYPE >::calculateLrms ( )

inline

Calculate mean element edge length.

\[ \Delta \boldsymbol\chi = \boldsymbol\chi^1 - \boldsymbol\chi^2 \]

\[ \Delta X = \mathbf{F} \Delta \boldsymbol\chi \]

\[ l_\textrm{rms} = \sqrt{\frac{1}{6} \sum_{i=0}^6 l_i^2 } = L_\textrm{rms}dl_\textrm{rms} \]

Definition at line 128 of file VolumeLengthQuality.hpp.

                                 {
    MoFEMFunctionBegin;
    if (deltaChi.size() != 3) {
      deltaChi.resize(3);
      deltaX.resize(3);
      dXdChiT.resize(3, 3);
    }
    lrmsSquared = 0;
    lrmsSquared0 = 0;
    std::fill(dXdChiT.data().begin(), dXdChiT.data().end(), 0);
    for (int ee = 0; ee < 6; ee++) {
      for (int dd = 0; dd < 3; dd++) {
        deltaChi[dd] = coordsEdges[6 * ee + dd] - coordsEdges[6 * ee + 3 + dd];
      }
      std::fill(deltaX.begin(), deltaX.end(), 0);
      for (int ii = 0; ii != 3; ++ii)
        for (int jj = 0; jj != 3; ++jj)
          deltaX(ii) += this->F(ii, jj) * deltaChi(jj);
 
      for (int dd = 0; dd < 3; dd++) {
        lrmsSquared += (1. / 6.) * deltaX[dd] * deltaX[dd];
        lrmsSquared0 += (1. / 6.) * deltaChi[dd] * deltaChi[dd];
      }
      for (int ii = 0; ii != 3; ++ii)
        for (int jj = 0; jj != 3; ++jj)
          dXdChiT(ii, jj) += deltaX(ii) * deltaChi(jj);
    }
    MoFEMFunctionReturn(0);
  }

◆ calculateP_PiolaKirchhoffI() [1/3]

template<typename TYPE >

virtual MoFEMErrorCode VolumeLengthQuality< TYPE >::calculateP_PiolaKirchhoffI ( )

inlinevirtual

Volume Length Quality.

Based on: Three‐dimensional brittle fracture: configurational‐force‐driven crack propagation International Journal for Numerical Methods in Engineering 97 (7), 531-550

\[ \mathcal{B}(a)=\frac{a}{2(1-\gamma)}-\ln{(a-\gamma)} \]

\[ q = q_0 b, \quad q_0 = 6\sqrt{2}\frac{V_0}{L^3_\textrm{rms}}, \quad b = \frac{\textrm{det}(\mathbf{F})}{\textrm{d}l^3_\textrm{rms}} \]

\[ \mathbf{P} = \mathcal{B}(a)\mathbf{Q}, \]

where \(a\) depending on problem could be \(q\) or \(b\).

Reimplemented from Smoother::CalculatePK1< TYPE >.

Definition at line 207 of file VolumeLengthQuality.hpp.

                                                      {
    MoFEMFunctionBegin;
 
    CHKERR getEdgesFromElemCoords();
    detF = determinantTensor3by3(this->F);
    if (this->invF.size1() != 3)
      this->invF.resize(3, 3);
 
    CHKERR invertTensor3by3(this->F, detF, this->invF);
    CHKERR calculateLrms();
    CHKERR calculateQ();
 
    double lrms03 = lrmsSquared0 * sqrt(lrmsSquared0);
    b = detF / (lrmsSquared * sqrt(lrmsSquared) / lrms03);
 
    currentVolume = detF * this->opPtr->getVolume();
    q = 6. * sqrt(2.) * currentVolume / (lrmsSquared * sqrt(lrmsSquared));
 
    if (this->P.size1() != 3)
      this->P.resize(3, 3);
 
    switch (tYpe) {
    case QUALITY:
      // Only use for testing, simple quality gradient
      noalias(this->P) = q * Q;
      break;
    case BARRIER_AND_QUALITY:
      // This is used form mesh smoothing
      tMp = q / (1.0 - gAmma) - 1.0 / (q - gAmma);
      noalias(this->P) = tMp * Q;
      break;
    case BARRIER_AND_CHANGE_QUALITY:
      // Works well with Arbitrary Lagrangian Formulation
      tMp = b / (1.0 - gAmma) - 1.0 / (b - gAmma);
      noalias(this->P) = tMp * Q;
      break;
    case BARRIER_AND_CHANGE_QUALITY_SCALED_BY_VOLUME:
      // When scaled by volume works well with ALE and face flipping.
      // Works well with smooth crack propagation
      tMp = currentVolume;
      tMp *= b / (1.0 - gAmma) - 1.0 / (b - gAmma);
      noalias(this->P) = tMp * Q;
      break;
    }
 
    // Divide by volume, to make units as they should be
    this->P *= aLpha / this->opPtr->getVolume();
 
    MoFEMFunctionReturn(0);
  }

◆ calculateP_PiolaKirchhoffI() [2/3]

template<typename TYPE >

virtual MoFEMErrorCode VolumeLengthQuality< TYPE >::calculateP_PiolaKirchhoffI	(	const NonlinearElasticElement::BlockData	block_data,
		boost::shared_ptr< const NumeredEntFiniteElement >	fe_ptr
	)

inlinevirtual

Volume Length Quality.

Based on: Three‐dimensional brittle fracture: configurational‐force‐driven crack propagation International Journal for Numerical Methods in Engineering 97 (7), 531-550

\[ \mathcal{B}(a)=\frac{a}{2(1-\gamma)}-\ln{(a-\gamma)} \]

\[ q = q_0 b, \quad q_0 = 6\sqrt{2}\frac{V_0}{L^3_\textrm{rms}}, \quad b = \frac{\textrm{det}(\mathbf{F})}{\textrm{d}l^3_\textrm{rms}} \]

\[ \mathbf{P} = \mathcal{B}(a)\mathbf{Q}, \]

where \(a\) depending on problem could be \(q\) or \(b\).

Reimplemented from NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI< TYPE >.

Definition at line 204 of file VolumeLengthQuality.hpp.

                                                             {
    MoFEMFunctionBegin;
 
    CHKERR getEdgesFromElemCoords();
    detF = determinantTensor3by3(this->F);
    if (this->invF.size1() != 3)
      this->invF.resize(3, 3);
 
    CHKERR invertTensor3by3(this->F, detF, this->invF);
    CHKERR calculateLrms();
    CHKERR calculateQ();
 
    double lrms03 = lrmsSquared0 * sqrt(lrmsSquared0);
    b = detF / (lrmsSquared * sqrt(lrmsSquared) / lrms03);
 
    currentVolume = detF * this->opPtr->getVolume();
    q = 6. * sqrt(2.) * currentVolume / (lrmsSquared * sqrt(lrmsSquared));
 
    if (this->P.size1() != 3)
      this->P.resize(3, 3);
 
    switch (tYpe) {
    case QUALITY:
      // Only use for testing, simple quality gradient
      noalias(this->P) = q * Q;
      break;
    case BARRIER_AND_QUALITY:
      // This is used form mesh smoothing
      tMp = q / (1.0 - gAmma) - 1.0 / (q - gAmma);
      noalias(this->P) = tMp * Q;
      break;
    case BARRIER_AND_CHANGE_QUALITY:
      // Works well with Arbitrary Lagrangian Formulation
      tMp = b / (1.0 - gAmma) - 1.0 / (b - gAmma);
      noalias(this->P) = tMp * Q;
      break;
    case BARRIER_AND_CHANGE_QUALITY_SCALED_BY_VOLUME:
      // When scaled by volume works well with ALE and face flipping.
      // Works well with smooth crack propagation
      tMp = currentVolume;
      tMp *= b / (1.0 - gAmma) - 1.0 / (b - gAmma);
      noalias(this->P) = tMp * Q;
      break;
    }
 
    // Divide by volume, to make units as they should be
    this->P *= aLpha / this->opPtr->getVolume();
 
    MoFEMFunctionReturn(0);
  }

◆ calculateP_PiolaKirchhoffI() [3/3]

template<typename TYPE >

virtual MoFEMErrorCode VolumeLengthQuality< TYPE >::calculateP_PiolaKirchhoffI	(	const NonlinearElasticElement::BlockData	block_data,
		boost::shared_ptr< const NumeredEntFiniteElement >	fe_ptr
	)

inlinevirtual

Volume Length Quality.

Based on: Three‐dimensional brittle fracture: configurational‐force‐driven crack propagation International Journal for Numerical Methods in Engineering 97 (7), 531-550

\[ \mathcal{B}(a)=\frac{a}{2(1-\gamma)}-\ln{(a-\gamma)} \]

\[ q = q_0 b, \quad q_0 = 6\sqrt{2}\frac{V_0}{L^3_\textrm{rms}}, \quad b = \frac{\textrm{det}(\mathbf{F})}{\textrm{d}l^3_\textrm{rms}} \]

\[ \mathbf{P} = \mathcal{B}(a)\mathbf{Q}, \]

where \(a\) depending on problem could be \(q\) or \(b\).

Reimplemented from NonlinearElasticElement::FunctionsToCalculatePiolaKirchhoffI< TYPE >.

Definition at line 205 of file VolumeLengthQuality.hpp.

                                                             {
    MoFEMFunctionBegin;
 
    CHKERR getEdgesFromElemCoords();
    detF = determinantTensor3by3(this->F);
    if (this->invF.size1() != 3)
      this->invF.resize(3, 3);
 
    CHKERR invertTensor3by3(this->F, detF, this->invF);
    CHKERR calculateLrms();
    CHKERR calculateQ();
 
    double lrms03 = lrmsSquared0 * sqrt(lrmsSquared0);
    b = detF / (lrmsSquared * sqrt(lrmsSquared) / lrms03);
 
    currentVolume = detF * this->opPtr->getVolume();
    q = 6. * sqrt(2.) * currentVolume / (lrmsSquared * sqrt(lrmsSquared));
 
    if (this->P.size1() != 3)
      this->P.resize(3, 3);
 
    switch (tYpe) {
    case QUALITY:
      // Only use for testing, simple quality gradient
      noalias(this->P) = q * Q;
      break;
    case BARRIER_AND_QUALITY:
      // This is used form mesh smoothing
      tMp = q / (1.0 - gAmma) - 1.0 / (q - gAmma);
      noalias(this->P) = tMp * Q;
      break;
    case BARRIER_AND_CHANGE_QUALITY:
      // Works well with Arbitrary Lagrangian Formulation
      tMp = b / (1.0 - gAmma) - 1.0 / (b - gAmma);
      noalias(this->P) = tMp * Q;
      break;
    case BARRIER_AND_CHANGE_QUALITY_SCALED_BY_VOLUME:
      // When scaled by volume works well with ALE and face flipping.
      // Works well with smooth crack propagation
      tMp = currentVolume;
      tMp *= b / (1.0 - gAmma) - 1.0 / (b - gAmma);
      noalias(this->P) = tMp * Q;
      break;
    }
 
    // Divide by volume, to make units as they should be
    this->P *= aLpha / this->opPtr->getVolume();
 
    MoFEMFunctionReturn(0);
  }

◆ calculateQ() [1/3]

template<typename TYPE >

MoFEMErrorCode VolumeLengthQuality< TYPE >::calculateQ ( )

inline

Calculate Q.

\[ \mathbf{Q} = \mathbf{F}^{-\mathsf{T}} - \frac{1}{2} \frac{1}{l^2_\textrm{rms}} \sum_i^6 \Delta\mathbf{X}_i \Delta\boldsymbol\chi_i^\mathsf{T} \]

Definition at line 174 of file VolumeLengthQuality.hpp.

                              {
    MoFEMFunctionBegin;
    Q.resize(3, 3, false);
    auto t_Q = getFTensor2FromArray3by3(Q, FTensor::Number<0>(), 0);
    auto t_invF = getFTensor2FromArray3by3(this->invF, FTensor::Number<0>(), 0);
    auto t_dXdChiT = getFTensor2FromArray3by3(dXdChiT, FTensor::Number<0>(), 0);
    t_Q(i, j) = t_invF(j, i) - 0.5 * t_dXdChiT(i, j) / lrmsSquared;
    MoFEMFunctionReturn(0);
  }

◆ calculateQ() [2/3]

template<typename TYPE >

MoFEMErrorCode VolumeLengthQuality< TYPE >::calculateQ ( )

inline

Calculate Q.

\[ \mathbf{Q} = \mathbf{F}^{-\mathsf{T}} - \frac{1}{2} \frac{1}{l^2_\textrm{rms}} \sum_i^6 \Delta\mathbf{X}_i \Delta\boldsymbol\chi_i^\mathsf{T} \]

Definition at line 171 of file VolumeLengthQuality.hpp.

                              {
    MoFEMFunctionBegin;
    Q.resize(3, 3, false);
    auto t_Q = getFTensor2FromArray3by3(Q, FTensor::Number<0>(), 0);
    auto t_invF = getFTensor2FromArray3by3(this->invF, FTensor::Number<0>(), 0);
    auto t_dXdChiT = getFTensor2FromArray3by3(dXdChiT, FTensor::Number<0>(), 0);
    t_Q(i, j) = t_invF(j, i) - 0.5 * t_dXdChiT(i, j) / lrmsSquared;
    MoFEMFunctionReturn(0);
  }

◆ calculateQ() [3/3]

template<typename TYPE >

MoFEMErrorCode VolumeLengthQuality< TYPE >::calculateQ ( )

inline

Calculate Q.

\[ \mathbf{Q} = \mathbf{F}^{-\mathsf{T}} - \frac{1}{2} \frac{1}{l^2_\textrm{rms}} \sum_i^6 \Delta\mathbf{X}_i \Delta\boldsymbol\chi_i^\mathsf{T} \]

Definition at line 172 of file VolumeLengthQuality.hpp.

                              {
    MoFEMFunctionBegin;
    Q.resize(3, 3, false);
    auto t_Q = getFTensor2FromArray3by3(Q, FTensor::Number<0>(), 0);
    auto t_invF = getFTensor2FromArray3by3(this->invF, FTensor::Number<0>(), 0);
    auto t_dXdChiT = getFTensor2FromArray3by3(dXdChiT, FTensor::Number<0>(), 0);
    t_Q(i, j) = t_invF(j, i) - 0.5 * t_dXdChiT(i, j) / lrmsSquared;
    MoFEMFunctionReturn(0);
  }

◆ getEdgesFromElemCoords() [1/3]

template<typename TYPE >

MoFEMErrorCode VolumeLengthQuality< TYPE >::getEdgesFromElemCoords ( )

inline

Get coordinates of edges using cannonical element numeration

Definition at line 83 of file VolumeLengthQuality.hpp.

                                          {
    MoFEMFunctionBegin;
    if (coordsEdges.empty()) {
      coordsEdges.resize(6 * 2 * 3, false);
    }
    cblas_dcopy(3, &this->opPtr->getCoords()[0 * 3], 1,
                &coordsEdges[0 * 3 * 2 + 0], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[1 * 3], 1,
                &coordsEdges[0 * 3 * 2 + 3], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[0 * 3], 1,
                &coordsEdges[1 * 3 * 2 + 0], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[2 * 3], 1,
                &coordsEdges[1 * 3 * 2 + 3], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[0 * 3], 1,
                &coordsEdges[2 * 3 * 2 + 0], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[3 * 3], 1,
                &coordsEdges[2 * 3 * 2 + 3], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[1 * 3], 1,
                &coordsEdges[3 * 3 * 2 + 0], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[2 * 3], 1,
                &coordsEdges[3 * 3 * 2 + 3], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[1 * 3], 1,
                &coordsEdges[4 * 3 * 2 + 0], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[3 * 3], 1,
                &coordsEdges[4 * 3 * 2 + 3], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[2 * 3], 1,
                &coordsEdges[5 * 3 * 2 + 0], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[3 * 3], 1,
                &coordsEdges[5 * 3 * 2 + 3], 1);
    MoFEMFunctionReturn(0);
  }

◆ getEdgesFromElemCoords() [2/3]

template<typename TYPE >

MoFEMErrorCode VolumeLengthQuality< TYPE >::getEdgesFromElemCoords ( )

inline

Get coordinates of edges using cannonical element numeration

Definition at line 80 of file VolumeLengthQuality.hpp.

                                          {
    MoFEMFunctionBegin;
    if (coordsEdges.empty()) {
      coordsEdges.resize(6 * 2 * 3, false);
    }
    cblas_dcopy(3, &this->opPtr->getCoords()[0 * 3], 1,
                &coordsEdges[0 * 3 * 2 + 0], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[1 * 3], 1,
                &coordsEdges[0 * 3 * 2 + 3], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[0 * 3], 1,
                &coordsEdges[1 * 3 * 2 + 0], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[2 * 3], 1,
                &coordsEdges[1 * 3 * 2 + 3], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[0 * 3], 1,
                &coordsEdges[2 * 3 * 2 + 0], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[3 * 3], 1,
                &coordsEdges[2 * 3 * 2 + 3], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[1 * 3], 1,
                &coordsEdges[3 * 3 * 2 + 0], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[2 * 3], 1,
                &coordsEdges[3 * 3 * 2 + 3], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[1 * 3], 1,
                &coordsEdges[4 * 3 * 2 + 0], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[3 * 3], 1,
                &coordsEdges[4 * 3 * 2 + 3], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[2 * 3], 1,
                &coordsEdges[5 * 3 * 2 + 0], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[3 * 3], 1,
                &coordsEdges[5 * 3 * 2 + 3], 1);
    MoFEMFunctionReturn(0);
  }

◆ getEdgesFromElemCoords() [3/3]

template<typename TYPE >

MoFEMErrorCode VolumeLengthQuality< TYPE >::getEdgesFromElemCoords ( )

inline

Get coordinates of edges using cannonical element numeration

Definition at line 81 of file VolumeLengthQuality.hpp.

                                          {
    MoFEMFunctionBegin;
    if (coordsEdges.empty()) {
      coordsEdges.resize(6 * 2 * 3, false);
    }
    cblas_dcopy(3, &this->opPtr->getCoords()[0 * 3], 1,
                &coordsEdges[0 * 3 * 2 + 0], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[1 * 3], 1,
                &coordsEdges[0 * 3 * 2 + 3], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[0 * 3], 1,
                &coordsEdges[1 * 3 * 2 + 0], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[2 * 3], 1,
                &coordsEdges[1 * 3 * 2 + 3], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[0 * 3], 1,
                &coordsEdges[2 * 3 * 2 + 0], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[3 * 3], 1,
                &coordsEdges[2 * 3 * 2 + 3], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[1 * 3], 1,
                &coordsEdges[3 * 3 * 2 + 0], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[2 * 3], 1,
                &coordsEdges[3 * 3 * 2 + 3], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[1 * 3], 1,
                &coordsEdges[4 * 3 * 2 + 0], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[3 * 3], 1,
                &coordsEdges[4 * 3 * 2 + 3], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[2 * 3], 1,
                &coordsEdges[5 * 3 * 2 + 0], 1);
    cblas_dcopy(3, &this->opPtr->getCoords()[3 * 3], 1,
                &coordsEdges[5 * 3 * 2 + 3], 1);
    MoFEMFunctionReturn(0);
  }

◆ getMaterialOptions() [1/3]

template<typename TYPE >

MoFEMErrorCode VolumeLengthQuality< TYPE >::getMaterialOptions ( )

inline

Definition at line 54 of file VolumeLengthQuality.hpp.

                                      {
    MoFEMFunctionBeginHot;
    PetscOptionsBegin(PETSC_COMM_WORLD, "",
                             "Get VolumeLengthQuality material options",
                             "none");
    CHKERR PetscOptionsEList(
        "-volume_length_type", "Volume length quality type", "",
        VolumeLengthQualityTypeNames, LASTOP_VOLUMELENGTHQUALITYTYPE,
        VolumeLengthQualityTypeNames[tYpe], &tYpe, PETSC_NULLPTR);
    CHKERR PetscOptionsScalar("-volume_length_alpha",
                              "volume length alpha parameter", "", aLpha,
                              &aLpha, PETSC_NULLPTR);
    CHKERR PetscOptionsScalar("-volume_length_gamma",
                              "volume length parameter (barrier)", "", gAmma,
                              &gAmma, PETSC_NULLPTR);
    PetscOptionsEnd();
    MoFEMFunctionReturnHot(0);
  }

◆ getMaterialOptions() [2/3]

template<typename TYPE >

MoFEMErrorCode VolumeLengthQuality< TYPE >::getMaterialOptions ( )

inline

Definition at line 52 of file VolumeLengthQuality.hpp.

                                      {
    MoFEMFunctionBeginHot;
    PetscOptionsBegin(PETSC_COMM_WORLD, "",
                      "Get VolumeLengthQuality material options", "none");
    CHKERR PetscOptionsEList(
        "-volume_length_type", "Volume length quality type", "",
        VolumeLengthQualityTypeNames, LASTOP_VOLUMELENGTHQUALITYTYPE,
        VolumeLengthQualityTypeNames[tYpe], &tYpe, PETSC_NULLPTR);
    CHKERR PetscOptionsScalar("-volume_length_alpha",
                              "volume length alpha parameter", "", aLpha,
                              &aLpha, PETSC_NULLPTR);
    CHKERR PetscOptionsScalar("-volume_length_gamma",
                              "volume length parameter (barrier)", "", gAmma,
                              &gAmma, PETSC_NULLPTR);
    PetscOptionsEnd();
    MoFEMFunctionReturnHot(0);
  }

◆ getMaterialOptions() [3/3]

template<typename TYPE >

MoFEMErrorCode VolumeLengthQuality< TYPE >::getMaterialOptions ( )

inline

Definition at line 52 of file VolumeLengthQuality.hpp.

                                      {
    MoFEMFunctionBeginHot;
    PetscOptionsBegin(PETSC_COMM_WORLD, "",
                             "Get VolumeLengthQuality material options",
                             "none");
    CHKERR PetscOptionsEList(
        "-volume_length_type", "Volume length quality type", "",
        VolumeLengthQualityTypeNames, LASTOP_VOLUMELENGTHQUALITYTYPE,
        VolumeLengthQualityTypeNames[tYpe], &tYpe, PETSC_NULLPTR);
    CHKERR PetscOptionsScalar("-volume_length_alpha",
                              "volume length alpha parameter", "", aLpha,
                              &aLpha, PETSC_NULLPTR);
    CHKERR PetscOptionsScalar("-volume_length_gamma",
                              "volume length parameter (barrier)", "", gAmma,
                              &gAmma, PETSC_NULLPTR);
    PetscOptionsEnd();
    MoFEMFunctionReturnHot(0);
  }

Member Data Documentation

◆ aLpha

template<typename TYPE >

double VolumeLengthQuality< TYPE >::aLpha

Definition at line 39 of file VolumeLengthQuality.hpp.

◆ b

template<typename TYPE >

TYPE VolumeLengthQuality< TYPE >::b

Definition at line 79 of file VolumeLengthQuality.hpp.

◆ coordsEdges

template<typename TYPE >

VectorDouble VolumeLengthQuality< TYPE >::coordsEdges

Definition at line 73 of file VolumeLengthQuality.hpp.

◆ currentVolume

template<typename TYPE >

TYPE VolumeLengthQuality< TYPE >::currentVolume

Definition at line 79 of file VolumeLengthQuality.hpp.

◆ deltaChi

template<typename TYPE >

VectorDouble VolumeLengthQuality< TYPE >::deltaChi

Definition at line 75 of file VolumeLengthQuality.hpp.

◆ deltaX

template<typename TYPE >

ublas::vector< TYPE > VolumeLengthQuality< TYPE >::deltaX

Definition at line 77 of file VolumeLengthQuality.hpp.

◆ detF

template<typename TYPE >

TYPE VolumeLengthQuality< TYPE >::detF

Definition at line 79 of file VolumeLengthQuality.hpp.

◆ dXdChiT

template<typename TYPE >

ublas::matrix< TYPE > VolumeLengthQuality< TYPE >::dXdChiT

Definition at line 78 of file VolumeLengthQuality.hpp.

◆ gAmma

template<typename TYPE >

double VolumeLengthQuality< TYPE >::gAmma

Definition at line 40 of file VolumeLengthQuality.hpp.

◆ i

template<typename TYPE >

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

Definition at line 42 of file VolumeLengthQuality.hpp.

◆ j

template<typename TYPE >

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

Definition at line 43 of file VolumeLengthQuality.hpp.

◆ k

template<typename TYPE >

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

Definition at line 44 of file VolumeLengthQuality.hpp.

◆ lrmsSquared

template<typename TYPE >

TYPE VolumeLengthQuality< TYPE >::lrmsSquared

Definition at line 79 of file VolumeLengthQuality.hpp.

◆ lrmsSquared0

template<typename TYPE >

double VolumeLengthQuality< TYPE >::lrmsSquared0

Definition at line 74 of file VolumeLengthQuality.hpp.

◆ Q

template<typename TYPE >

ublas::matrix< TYPE > VolumeLengthQuality< TYPE >::Q

Definition at line 78 of file VolumeLengthQuality.hpp.

◆ q

template<typename TYPE >

TYPE VolumeLengthQuality< TYPE >::q

Definition at line 79 of file VolumeLengthQuality.hpp.

◆ tMp

template<typename TYPE >

TYPE VolumeLengthQuality< TYPE >::tMp

Definition at line 79 of file VolumeLengthQuality.hpp.

◆ tYpe

template<typename TYPE >

int VolumeLengthQuality< TYPE >::tYpe

Definition at line 38 of file VolumeLengthQuality.hpp.

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

tools/src/VolumeLengthQuality.hpp
users_modules/basic_finite_elements/nonlinear_elastic_materials/src/VolumeLengthQuality.hpp
users_modules/basic_finite_elements/src/VolumeLengthQuality.hpp

Public Member Functions

Public Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ VolumeLengthQuality() [1/3]

◆ VolumeLengthQuality() [2/3]

◆ VolumeLengthQuality() [3/3]

Member Function Documentation

◆ calculateLrms() [1/3]

◆ calculateLrms() [2/3]

◆ calculateLrms() [3/3]

◆ calculateP_PiolaKirchhoffI() [1/3]

◆ calculateP_PiolaKirchhoffI() [2/3]

◆ calculateP_PiolaKirchhoffI() [3/3]

◆ calculateQ() [1/3]

◆ calculateQ() [2/3]

◆ calculateQ() [3/3]

◆ getEdgesFromElemCoords() [1/3]

◆ getEdgesFromElemCoords() [2/3]

◆ getEdgesFromElemCoords() [3/3]

◆ getMaterialOptions() [1/3]

◆ getMaterialOptions() [2/3]

◆ getMaterialOptions() [3/3]

Member Data Documentation

◆ aLpha

◆ b

◆ coordsEdges

◆ currentVolume

◆ deltaChi

◆ deltaX

◆ detF

◆ dXdChiT

◆ gAmma

◆ i

◆ j

◆ k

◆ lrmsSquared

◆ lrmsSquared0

◆ Q

◆ q

◆ tMp

◆ tYpe