#include "tutorials/cor-7/src/ElasticityMixedFormulation.hpp"

Collaboration diagram for DataAtIntegrationPts:

Public Member Functions
	DataAtIntegrationPts (MoFEM::Interface &m_field)

MoFEMErrorCode	getParameters ()

MoFEMErrorCode	getBlockData (BlockData &data)

MoFEMErrorCode	setBlocks ()

	DataAtIntegrationPts (MoFEM::Interface &m_field)

	DataAtIntegrationPts ()

Public Attributes
boost::shared_ptr< MatrixDouble >	gradDispPtr

boost::shared_ptr< VectorDouble >	pPtr

FTensor::Ddg< double, 3, 3 >	tD

double	pOisson

double	yOung

double	lAmbda

double	mU

std::map< int, BlockData >	setOfBlocksData

SmartPetscObj< Vec >	petscVec

double	blockPermittivity

double	blockChrgDens

boost::shared_ptr< MatrixDouble >	smallStrainMat

boost::shared_ptr< MatrixDouble >	hMat

boost::shared_ptr< MatrixDouble >	FMat

boost::shared_ptr< MatrixDouble >	HMat

boost::shared_ptr< VectorDouble >	detHVec

boost::shared_ptr< MatrixDouble >	invHMat

boost::shared_ptr< MatrixDouble >	cauchyStressMat

boost::shared_ptr< MatrixDouble >	stiffnessMat

boost::shared_ptr< VectorDouble >	energyVec

boost::shared_ptr< MatrixDouble >	eshelbyStressMat

boost::shared_ptr< MatrixDouble >	eshelbyStress_dx

Range	forcesOnlyOnEntitiesRow

Range	forcesOnlyOnEntitiesCol

Private Attributes
MoFEM::Interface &	mField

Detailed Description

Examples: ElasticityMixedFormulation.hpp, EshelbianPlasticity.cpp, HookeElement.cpp, HookeElement.hpp, and elasticity_mixed_formulation.cpp.

Definition at line 20 of file ElasticityMixedFormulation.hpp.

Constructor & Destructor Documentation

◆ DataAtIntegrationPts() [1/3]

DataAtIntegrationPts::DataAtIntegrationPts ( MoFEM::Interface & m_field )

inline

Definition at line 33 of file ElasticityMixedFormulation.hpp.

                                                : mField(m_field) {
 
    // Setting default values for coeffcients
    gradDispPtr = boost::shared_ptr<MatrixDouble>(new MatrixDouble());
    pPtr = boost::shared_ptr<VectorDouble>(new VectorDouble());
 
    ierr = setBlocks();
    CHKERRABORT(PETSC_COMM_WORLD, ierr);
  }

◆ DataAtIntegrationPts() [2/3]

DataAtIntegrationPts::DataAtIntegrationPts ( MoFEM::Interface & m_field )

inline

Definition at line 62 of file electrostatics.hpp.

                                                {
    blockPermittivity = 0.0;
    blockChrgDens = 0.0;
    PetscInt ghosts[2] = {0, 1};
    if (!m_field.get_comm_rank())
      petscVec = createGhostVector(m_field.get_comm(), 2, 2, 0, ghosts);
    else
      petscVec = createGhostVector(m_field.get_comm(), 0, 2, 2, ghosts);
  }

◆ DataAtIntegrationPts() [3/3]

DataAtIntegrationPts::DataAtIntegrationPts ( )

inline

Examples: ElasticityMixedFormulation.hpp.

Definition at line 96 of file HookeElement.hpp.

                           {
 
      smallStrainMat = boost::shared_ptr<MatrixDouble>(new MatrixDouble());
      hMat = boost::shared_ptr<MatrixDouble>(new MatrixDouble());
      FMat = boost::shared_ptr<MatrixDouble>(new MatrixDouble());
 
      HMat = boost::shared_ptr<MatrixDouble>(new MatrixDouble());
      detHVec = boost::shared_ptr<VectorDouble>(new VectorDouble());
      invHMat = boost::shared_ptr<MatrixDouble>(new MatrixDouble());
 
      cauchyStressMat = boost::shared_ptr<MatrixDouble>(new MatrixDouble());
      stiffnessMat = boost::shared_ptr<MatrixDouble>(new MatrixDouble());
      energyVec = boost::shared_ptr<VectorDouble>(new VectorDouble());
      eshelbyStressMat = boost::shared_ptr<MatrixDouble>(new MatrixDouble());
 
      eshelbyStress_dx = boost::shared_ptr<MatrixDouble>(new MatrixDouble());
    }

Member Function Documentation

◆ getBlockData()

MoFEMErrorCode DataAtIntegrationPts::getBlockData ( BlockData & data )

inline

Examples: ElasticityMixedFormulation.hpp.

Definition at line 51 of file ElasticityMixedFormulation.hpp.

                                               {
    MoFEMFunctionBegin;
 
    yOung = data.yOung;
    pOisson = data.pOisson;
    lAmbda = (yOung * pOisson) / ((1. + pOisson) * (1. - 2. * pOisson));
    mU = yOung / (2. * (1. + pOisson));
 
    FTensor::Index<'i', 3> i;
    FTensor::Index<'j', 3> j;
    FTensor::Index<'k', 3> k;
    FTensor::Index<'l', 3> l;
 
    tD(i, j, k, l) = 0.;
 
    tD(0, 0, 0, 0) = 2 * mU;
    tD(0, 1, 0, 1) = mU;
    tD(0, 1, 1, 0) = mU;
    tD(0, 2, 0, 2) = mU;
    tD(0, 2, 2, 0) = mU;
    tD(1, 0, 0, 1) = mU;
    tD(1, 0, 1, 0) = mU;
    tD(1, 1, 1, 1) = 2 * mU;
    tD(1, 2, 1, 2) = mU;
    tD(1, 2, 2, 1) = mU;
    tD(2, 0, 0, 2) = mU;
    tD(2, 0, 2, 0) = mU;
    tD(2, 1, 1, 2) = mU;
    tD(2, 1, 2, 1) = mU;
    tD(2, 2, 2, 2) = 2 * mU;
 
    MoFEMFunctionReturn(0);
  }

◆ getParameters()

MoFEMErrorCode DataAtIntegrationPts::getParameters ( )

inline

Examples: ElasticityMixedFormulation.hpp, and elasticity_mixed_formulation.cpp.

Definition at line 43 of file ElasticityMixedFormulation.hpp.

                                 {
    MoFEMFunctionBegin; // They will be overwriten by BlockData
    PetscOptionsBegin(PETSC_COMM_WORLD, "", "Problem", "none");
 
    PetscOptionsEnd();
    MoFEMFunctionReturn(0);
  }

◆ setBlocks()

MoFEMErrorCode DataAtIntegrationPts::setBlocks ( )

inline

Examples: ElasticityMixedFormulation.hpp.

Definition at line 85 of file ElasticityMixedFormulation.hpp.

                             {
    MoFEMFunctionBegin;
    for (_IT_CUBITMESHSETS_BY_BCDATA_TYPE_FOR_LOOP_(
             mField, BLOCKSET | MAT_ELASTICSET, it)) {
      Mat_Elastic mydata;
      CHKERR it->getAttributeDataStructure(mydata);
      int id = it->getMeshsetId();
      EntityHandle meshset = it->getMeshset();
      CHKERR mField.get_moab().get_entities_by_type(
          meshset, MBTET, setOfBlocksData[id].tEts, true);
      setOfBlocksData[id].iD = id;
      setOfBlocksData[id].yOung = mydata.data.Young;
      setOfBlocksData[id].pOisson = mydata.data.Poisson;
    }
    MoFEMFunctionReturn(0);
  }