#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, HookeElement.hpp, and HookeInternalStressElement.hpp.

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 52 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.

Definition at line 43 of file ElasticityMixedFormulation.hpp.

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

◆ setBlocks()

MoFEMErrorCode DataAtIntegrationPts::setBlocks ( )

inline

Examples: ElasticityMixedFormulation.hpp.

Definition at line 86 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);
   }