Classes
struct	AnalyticalDisp

struct	AnalyticalForces

struct	AnalyticalOptions

struct	ConstantArea
	Constant area constrains. More...

struct	CPMeshCut

struct	CPSolvers

struct	CrackFrontSingularBase

struct	CrackPropagation

struct	FirendVolumeOnSide

struct	GetSmoothingElementsSkin

struct	GriffithForceElement
	Implementation of Griffith element. More...

struct	MWLSApprox

struct	OpAleLhsWithDensitySingularElement_dX_dX

struct	OpAleLhsWithDensitySingularElement_dx_dX

struct	OpAnalyticalMaterialTraction

struct	OpAnalyticalSpatialTraction

struct	OpGetCrackFrontCommonDataAtGaussPts

struct	OpGetCrackFrontDataAtGaussPts

struct	OpGetCrackFrontDataGradientAtGaussPts

struct	OpGetDensityFieldForTesting
	Op to generate artificial density field. More...

struct	OpLhsBoneExplicitDerivariveWithHooke_dX
	Calculate explicit derivative of energy. More...

struct	OpLhsBoneExplicitDerivariveWithHooke_dx

struct	OpPostProcDisplacements

struct	OpPrint

struct	OpRhsBoneExplicitDerivariveWithHooke
	Calculate explicit derivative of energy. More...

struct	OpSetTagRangeOnSkin
	Mark crack surfaces on skin. More...

struct	OpTransfromSingularBaseFunctions

class	TwoType

Typedefs
typedef CrackFrontSingularBase< NonlinearElasticElement::MyVolumeFE, VolumeElementForcesAndSourcesCore >	CrackFrontElement

Enumerations
enum	Materials { HOOKE , KIRCHHOFF , NEOHOOKEAN , BONEHOOKE , LASTOP }

enum	AdolcTags { ELASTIC_TAG = 1 , MATERIAL_TAG , GRIFFITH_FORCE_TAG , GRIFFITH_CONSTRAINS_TAG , CONSTANT_AREA_TAG , ARC_LENGTH_TAG , FRONT_TANGENT , EXTERIOR_DERIVATIVE_TAG , SMOOTHING_TAG , SURFACE_SLIDING_TAG , EDGE_SLIDING_TAG }
	Tapes numbers used by ADOL-C. More...

enum	tangent_tests { MWLS_STRESS_TAN , MWLS_DENSITY_TAN , MWLS_GRIFFITH_TAN , LASTTAN }
	Names of tangent matrices tests. More...

Functions
MoFEMErrorCode	clean_pcomms (moab::Interface &moab, boost::shared_ptr< WrapMPIComm > moab_comm_wrap)

MoFEMErrorCode	broadcast_entities (moab::Interface &moab, moab::Interface &moab_tmp, const int from_proc, Range &entities, const bool adjacencies, const bool tags)

static double	calMax (double a, double b, double r)

static double	diffCalMax_a (double a, double b, double r)

MoFEMErrorCode	clean_pcomms (moab::Interface &moab, boost::shared_ptr< WrapMPIComm > moab_comm_wrap)

static MoFEMErrorCode	elastic_snes_rhs (SNES snes, Vec x, Vec f, void *ctx)

static MoFEMErrorCode	elastic_snes_mat (SNES snes, Vec x, Mat A, Mat B, void *ctx)

static MoFEMErrorCode	propagation_snes_rhs (SNES snes, Vec x, Vec f, void *ctx)

static MoFEMErrorCode	propagation_snes_mat (SNES snes, Vec x, Mat A, Mat B, void *ctx)

Variables
static std::map< long int, MatrixDouble >	mapRefCoords

const char *	materials_list []

Typedef Documentation

◆ CrackFrontElement

typedef CrackFrontSingularBase<NonlinearElasticElement::MyVolumeFE, VolumeElementForcesAndSourcesCore> FractureMechanics::CrackFrontElement

Definition at line 402 of file CrackFrontElement.hpp.

Enumeration Type Documentation

◆ AdolcTags

enum FractureMechanics::AdolcTags

Tapes numbers used by ADOL-C.

Enumerator
ELASTIC_TAG
MATERIAL_TAG
GRIFFITH_FORCE_TAG
GRIFFITH_CONSTRAINS_TAG
CONSTANT_AREA_TAG
ARC_LENGTH_TAG
FRONT_TANGENT
EXTERIOR_DERIVATIVE_TAG
SMOOTHING_TAG
SURFACE_SLIDING_TAG
EDGE_SLIDING_TAG

Definition at line 36 of file CrackPropagation.hpp.

               {
  ELASTIC_TAG = 1,
  MATERIAL_TAG,
  GRIFFITH_FORCE_TAG,
  GRIFFITH_CONSTRAINS_TAG,
  CONSTANT_AREA_TAG,
  ARC_LENGTH_TAG,
  FRONT_TANGENT,
  EXTERIOR_DERIVATIVE_TAG,
  SMOOTHING_TAG,
  SURFACE_SLIDING_TAG,
  EDGE_SLIDING_TAG
};

◆ Materials

enum FractureMechanics::Materials

Enumerator
HOOKE
KIRCHHOFF
NEOHOOKEAN
BONEHOOKE
LASTOP

Definition at line 23 of file CrackFrontElement.hpp.

23{ HOOKE, KIRCHHOFF, NEOHOOKEAN, BONEHOOKE, LASTOP };

FractureMechanics::BONEHOOKE

@ BONEHOOKE

Definition: CrackFrontElement.hpp:23

FractureMechanics::LASTOP

@ LASTOP

Definition: CrackFrontElement.hpp:23

FractureMechanics::HOOKE

@ HOOKE

Definition: CrackFrontElement.hpp:23

FractureMechanics::KIRCHHOFF

@ KIRCHHOFF

Definition: CrackFrontElement.hpp:23

FractureMechanics::NEOHOOKEAN

@ NEOHOOKEAN

Definition: CrackFrontElement.hpp:23

◆ tangent_tests

enum FractureMechanics::tangent_tests

Names of tangent matrices tests.

Enumerator
MWLS_STRESS_TAN
MWLS_DENSITY_TAN
MWLS_GRIFFITH_TAN
LASTTAN

Definition at line 54 of file CrackPropagation.hpp.

                   {
  MWLS_STRESS_TAN,
  MWLS_DENSITY_TAN,
  MWLS_GRIFFITH_TAN,
  LASTTAN
};

Function Documentation

◆ broadcast_entities()

MoFEMErrorCode FractureMechanics::broadcast_entities	(	moab::Interface &	moab,
		moab::Interface &	moab_tmp,
		const int	from_proc,
		Range &	entities,
		const bool	adjacencies,
		const bool	tags
	)

Definition at line 147 of file CrackPropagation.cpp.

                                                                           {
  ErrorCode result = MB_SUCCESS;
  int success;
  int buff_size;
 
  auto add_verts = [&](Range &sent_ents) {
    // Get the verts adj to these entities, since we'll have to send those
    // too
 
    // First check sets
    std::pair<Range::const_iterator, Range::const_iterator> set_range =
        sent_ents.equal_range(MBENTITYSET);
    ErrorCode result = MB_SUCCESS, tmp_result;
    for (Range::const_iterator rit = set_range.first; rit != set_range.second;
         ++rit) {
      tmp_result = moab_tmp.get_entities_by_type(*rit, MBVERTEX, sent_ents);
      CHK_MOAB_THROW(tmp_result, "Failed to get contained verts");
    }
 
    // Now non-sets
    Range tmp_ents;
    std::copy(sent_ents.begin(), set_range.first, range_inserter(tmp_ents));
    result = moab_tmp.get_adjacencies(tmp_ents, 0, false, sent_ents,
                                      moab::Interface::UNION);
    CHK_MOAB_THROW(result, "Failed to get vertices adj to ghosted ents");
 
    return result;
  };
 
  ParallelComm *pcomm = ParallelComm::get_pcomm(&moab_tmp, MYPCOMM_INDEX);
  auto &procConfig = pcomm->proc_config();
  const int MAX_BCAST_SIZE = (1 << 28);
 
  ParallelComm::Buffer buff(ParallelComm::INITIAL_BUFF_SIZE);
  buff.reset_ptr(sizeof(int));
  if ((int)procConfig.proc_rank() == from_proc) {
 
    result = add_verts(entities);
    CHK_MOAB_THROW(result, "Failed to add adj vertices");
 
    buff.reset_ptr(sizeof(int));
    result = pcomm->pack_buffer(entities, adjacencies, tags, false, -1, &buff);
    CHK_MOAB_THROW(result,
                   "Failed to compute buffer size in broadcast_entities");
    buff.set_stored_size();
    buff_size = buff.buff_ptr - buff.mem_ptr;
  }
 
  success =
      MPI_Bcast(&buff_size, 1, MPI_INT, from_proc, procConfig.proc_comm());
  if (MPI_SUCCESS != success) {
    THROW_MESSAGE("MPI_Bcast of buffer size failed");
  }
 
  if (!buff_size) // No data
    return MB_SUCCESS;
 
  if ((int)procConfig.proc_rank() != from_proc)
    buff.reserve(buff_size);
 
  size_t offset = 0;
  while (buff_size) {
    int sz = std::min(buff_size, MAX_BCAST_SIZE);
    success = MPI_Bcast(buff.mem_ptr + offset, sz, MPI_UNSIGNED_CHAR, from_proc,
                        procConfig.proc_comm());
    if (MPI_SUCCESS != success) {
      THROW_MESSAGE("MPI_Bcast of buffer failed");
    }
 
    offset += sz;
    buff_size -= sz;
  }
 
  std::vector<std::vector<EntityHandle>> dum1a, dum1b;
  std::vector<std::vector<int>> dum1p;
  std::vector<EntityHandle> dum2, dum4;
  std::vector<unsigned int> dum3;
  buff.reset_ptr(sizeof(int));
  if ((int)procConfig.proc_rank() == from_proc) {
    ParallelComm *pcomm = ParallelComm::get_pcomm(&moab, MYPCOMM_INDEX);
    result = pcomm->unpack_buffer(buff.buff_ptr, false, from_proc, -1, dum1a,
                                  dum1b, dum1p, dum2, dum2, dum3, dum4);
  } else {
    result = pcomm->unpack_buffer(buff.buff_ptr, false, from_proc, -1, dum1a,
                                  dum1b, dum1p, dum2, dum2, dum3, dum4);
  }
  CHK_MOAB_THROW(result, "Failed to unpack buffer in broadcast_entities");
  std::copy(dum4.begin(), dum4.end(), range_inserter(entities));
 
  return MB_SUCCESS;
};

◆ calMax()

static double FractureMechanics::calMax	(	double	a,
		double	b,
		double	r
	)

static

Definition at line 24 of file GriffithForceElement.hpp.

                                                   {
  return (a + b + (1 / r) * pow(fabs(a - b), r)) / 2;
}

◆ clean_pcomms() [1/2]

MoFEMErrorCode FractureMechanics::clean_pcomms	(	moab::Interface &	moab,
		boost::shared_ptr< WrapMPIComm >	moab_comm_wrap
	)

◆ clean_pcomms() [2/2]

MoFEMErrorCode FractureMechanics::clean_pcomms	(	moab::Interface &	moab,
		boost::shared_ptr< WrapMPIComm >	moab_comm_wrap
	)

Definition at line 126 of file CrackPropagation.cpp.

                                                                         {
  MoFEMFunctionBegin;
  std::vector<ParallelComm *> list_pcomms;
  ParallelComm::get_all_pcomm(&moab, list_pcomms);
  if (list_pcomms[MYPCOMM_INDEX]) {
    CHKERR moab.tag_delete(list_pcomms[MYPCOMM_INDEX]->pstatus_tag());
    CHKERR moab.tag_delete(list_pcomms[MYPCOMM_INDEX]->sharedps_tag());
    CHKERR moab.tag_delete(list_pcomms[MYPCOMM_INDEX]->sharedhs_tag());
    CHKERR moab.tag_delete(list_pcomms[MYPCOMM_INDEX]->sharedp_tag());
    CHKERR moab.tag_delete(list_pcomms[MYPCOMM_INDEX]->sharedh_tag());
  }
  for (auto p : list_pcomms) {
    delete p;
  }
  int pcomm_id = MYPCOMM_INDEX;
  list_pcomms[MYPCOMM_INDEX] =
      new ParallelComm(&moab, moab_comm_wrap->get_comm(), &pcomm_id);
  MoFEMFunctionReturn(0);
}

◆ diffCalMax_a()

static double FractureMechanics::diffCalMax_a	(	double	a,
		double	b,
		double	r
	)

static

Definition at line 28 of file GriffithForceElement.hpp.

                                                         {
  double sgn = ((a - b) == 0) ? 0 : (((a - b) < 0) ? -1 : 1);
  return (1 + pow(fabs(a - b), r - 1) * sgn * (+1)) / 2;
}

◆ elastic_snes_mat()

static MoFEMErrorCode FractureMechanics::elastic_snes_mat	(	SNES	snes,
		Vec	x,
		Mat	A,
		Mat	B,
		void *	ctx
	)

static

Definition at line 8227 of file CrackPropagation.cpp.

                                                  {
  auto arc_snes_ctx =
      reinterpret_cast<CrackPropagation::ArcLengthSnesCtx *>(ctx);
  MoFEMFunctionBegin;
  CHKERR SnesMat(snes, x, A, B, ctx);
  CHKERR MatDiagonalScale(B, arc_snes_ctx->getVecDiagM(), PETSC_NULL);
  MoFEMFunctionReturn(0);
};

◆ elastic_snes_rhs()

static MoFEMErrorCode FractureMechanics::elastic_snes_rhs	(	SNES	snes,
		Vec	x,
		Vec	f,
		void *	ctx
	)

static

Definition at line 8218 of file CrackPropagation.cpp.

                                                                           {
  auto arc_snes_ctx =
      reinterpret_cast<CrackPropagation::ArcLengthSnesCtx *>(ctx);
  MoFEMFunctionBegin;
  CHKERR SnesRhs(snes, x, f, ctx);
  CHKERR VecPointwiseMult(f, f, arc_snes_ctx->getVecDiagM());
  MoFEMFunctionReturn(0);
}

◆ propagation_snes_mat()

static MoFEMErrorCode FractureMechanics::propagation_snes_mat	(	SNES	snes,
		Vec	x,
		Mat	A,
		Mat	B,
		void *	ctx
	)

static

Definition at line 8491 of file CrackPropagation.cpp.

                                                      {
  auto arc_snes_ctx =
      reinterpret_cast<CrackPropagation::ArcLengthSnesCtx *>(ctx);
  MoFEMFunctionBegin;
 
  CHKERR SnesMat(snes, x, A, B, ctx);
 
  CHKERR MatDiagonalScale(B, arc_snes_ctx->getVecDiagM(), PETSC_NULL);
  CHKERR VecPointwiseMult(arc_snes_ctx->getArcPtr()->F_lambda,
                          arc_snes_ctx->getArcPtr()->F_lambda,
                          arc_snes_ctx->getVecDiagM());
 
  MoFEMFunctionReturn(0);
};

◆ propagation_snes_rhs()

static MoFEMErrorCode FractureMechanics::propagation_snes_rhs	(	SNES	snes,
		Vec	x,
		Vec	f,
		void *	ctx
	)

static

Definition at line 8482 of file CrackPropagation.cpp.

                                                                               {
  auto arc_snes_ctx =
      reinterpret_cast<CrackPropagation::ArcLengthSnesCtx *>(ctx);
  MoFEMFunctionBegin;
  CHKERR SnesRhs(snes, x, f, ctx);
  CHKERR VecPointwiseMult(f, f, arc_snes_ctx->getVecDiagM());
  MoFEMFunctionReturn(0);
};

Variable Documentation

◆ mapRefCoords

std::map<long int, MatrixDouble> FractureMechanics::mapRefCoords

static

Definition at line 53 of file CrackFrontElement.cpp.

◆ materials_list

const char* FractureMechanics::materials_list[]

Initial value:

= {"HOOKE", "KIRCHHOFF", "NEOHOOKEAN",

"BONEHOOKE"}

Definition at line 492 of file CrackPropagation.cpp.

Classes

Typedefs

Enumerations

Functions

Variables

Typedef Documentation

◆ CrackFrontElement

Enumeration Type Documentation

◆ AdolcTags

◆ Materials

◆ tangent_tests

Function Documentation

◆ broadcast_entities()

◆ calMax()

◆ clean_pcomms() [1/2]

◆ clean_pcomms() [2/2]

◆ diffCalMax_a()

◆ elastic_snes_mat()

◆ elastic_snes_rhs()

◆ propagation_snes_mat()

◆ propagation_snes_rhs()

Variable Documentation

◆ mapRefCoords

◆ materials_list