EshelbianPlasticity::HMHNeohookean Struct Reference

Inheritance diagram for EshelbianPlasticity::HMHNeohookean:

Collaboration diagram for EshelbianPlasticity::HMHNeohookean:

Classes
struct	BlockData
struct	OpGetScale
struct	OpJacobian
struct	OpSpatialPhysical
struct	OpSpatialPhysical_du_du

Public Member Functions
	HMHNeohookean (MoFEM::Interface &m_field, const double c10, const double K)
auto	getMaterialParameters (EntityHandle ent)
VolUserDataOperator *	returnOpSetScale (boost::shared_ptr< double > scale_ptr, boost::shared_ptr< PhysicalEquations > physics_ptr)
virtual EshelbianPlasticity::OpJacobian *	returnOpJacobian (const int tag, const bool eval_rhs, const bool eval_lhs, boost::shared_ptr< DataAtIntegrationPts > data_ptr, boost::shared_ptr< PhysicalEquations > physics_ptr)
MoFEMErrorCode	getOptions ()
MoFEMErrorCode	extractBlockData (Sev sev)
MoFEMErrorCode	extractBlockData (std::vector< const CubitMeshSets * > meshset_vec_ptr, Sev sev)
MoFEMErrorCode	recordTape (const int tape, DTensor2Ptr *t_h_ptr)
virtual VolUserDataOperator *	returnOpSpatialPhysical (const std::string &field_name, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const double alpha_u)
VolUserDataOperator *	returnOpSpatialPhysical_du_du (std::string row_field, std::string col_field, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const double alpha)
Public Member Functions inherited from EshelbianPlasticity::PhysicalEquations
	PhysicalEquations ()=delete
	PhysicalEquations (const int size_active, const int size_dependent)
virtual	~PhysicalEquations ()=default
virtual VolUserDataOperator *	returnOpCalculateEnergy (boost::shared_ptr< DataAtIntegrationPts > data_ptr, boost::shared_ptr< double > total_energy_ptr)
virtual VolUserDataOperator *	returnOpCalculateStretchFromStress (boost::shared_ptr< DataAtIntegrationPts > data_ptr, boost::shared_ptr< PhysicalEquations > physics_ptr)
virtual VolUserDataOperator *	returnOpCalculateVarStretchFromStress (boost::shared_ptr< DataAtIntegrationPts > data_ptr, boost::shared_ptr< PhysicalEquations > physics_ptr)
DTensor2Ptr	get_P ()
DTensor3Ptr	get_P_dh0 ()
DTensor3Ptr	get_P_dh1 ()
DTensor3Ptr	get_P_dh2 ()
DTensor2Ptr	get_h ()

Static Public Attributes
static constexpr int	numberOfActiveVariables = 0
static constexpr int	numberOfDependentVariables = 0

Private Attributes
MoFEM::Interface &	mField
double	c10_default
double	K_default
double	alphaGradU
std::vector< BlockData >	blockData
double	eqScaling = 1.

Additional Inherited Members
Public Types inherited from EshelbianPlasticity::PhysicalEquations
typedef FTensor::Tensor1< adouble, 3 >	ATensor1
typedef FTensor::Tensor2< adouble, 3, 3 >	ATensor2
typedef FTensor::Tensor3< adouble, 3, 3, 3 >	ATensor3
typedef FTensor::Tensor1< double, 3 >	DTensor1
typedef FTensor::Tensor2< double, 3, 3 >	DTensor2
typedef FTensor::Tensor3< double, 3, 3, 3 >	DTensor3
typedef FTensor::Tensor0< FTensor::PackPtr< double *, 1 > >	DTensor0Ptr
typedef FTensor::Tensor2< FTensor::PackPtr< double *, 1 >, 3, 3 >	DTensor2Ptr
typedef FTensor::Tensor3< FTensor::PackPtr< double *, 1 >, 3, 3, 3 >	DTensor3Ptr
Static Public Member Functions inherited from EshelbianPlasticity::PhysicalEquations
template<int S>
static DTensor2Ptr	get_VecTensor2 (std::vector< double > &v)
template<int S>
static DTensor0Ptr	get_VecTensor0 (std::vector< double > &v)
template<int S0>
static DTensor3Ptr	get_vecTensor3 (std::vector< double > &v, const int nba)
Public Attributes inherited from EshelbianPlasticity::PhysicalEquations
std::vector< double >	activeVariables
std::vector< double >	dependentVariablesPiola
std::vector< double >	dependentVariablesPiolaDirevatives

Detailed Description

Definition at line 14 of file HMHNeohookean.cpp.

Constructor & Destructor Documentation

HMHNeohookean()

EshelbianPlasticity::HMHNeohookean::HMHNeohookean ( MoFEM::Interface &  m_field, const double  c10, const double  K  )

inline

Definition at line 19 of file HMHNeohookean.cpp.

      : PhysicalEquations(numberOfActiveVariables, numberOfDependentVariables),
        mField(m_field), c10_default(c10), K_default(K) {
 
    CHK_THROW_MESSAGE(getOptions(), "get options failed");
    CHK_THROW_MESSAGE(extractBlockData(Sev::inform),
                      "extract block data failed");
 
#ifdef NEOHOOKEAN_SCALING
    if (blockData.size()) {
      double ave_K = 0;
      for (auto b : blockData) {
        ave_K += b.K;
      }
      eqScaling = ave_K / blockData.size();
    }
#endif
 
    MOFEM_LOG("EP", Sev::inform) << "Neo-Hookean scaling = " << eqScaling;
 
    if (EshelbianCore::stretchSelector != StretchSelector::LOG) {
      CHK_THROW_MESSAGE(MOFEM_DATA_INCONSISTENCY,
                        "Stretch selector is not equal to LOG");
    } else {
      if (EshelbianCore::exponentBase != exp(1)) {
        CHK_THROW_MESSAGE(MOFEM_DATA_INCONSISTENCY,
                          "Exponent base is not equal to exp(1)");
      }
    }
  }

Member Function Documentation

extractBlockData() [1/2]

MoFEMErrorCode EshelbianPlasticity::HMHNeohookean::extractBlockData ( Sev  sev )

inline

Definition at line 128 of file HMHNeohookean.cpp.

                                           {
    return extractBlockData(
 
        mField.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(std::regex(
 
            (boost::format("%s(.*)") % "MAT_NEOHOOKEAN").str()
 
                )),
 
        sev);
  }

extractBlockData() [2/2]

MoFEMErrorCode EshelbianPlasticity::HMHNeohookean::extractBlockData ( std::vector< const CubitMeshSets * >  meshset_vec_ptr, Sev  sev  )

inline

Definition at line 141 of file HMHNeohookean.cpp.

                            {
    MoFEMFunctionBegin;
 
    for (auto m : meshset_vec_ptr) {
      MOFEM_LOG("EP", sev) << *m;
      std::vector<double> block_data;
      CHKERR m->getAttributes(block_data);
      if (block_data.size() != 2) {
        SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
                "Expected that block has two attribute");
      }
      auto get_block_ents = [&]() {
        Range ents;
        CHKERR mField.get_moab().get_entities_by_handle(m->meshset, ents, true);
        return ents;
      };
 
      double c10 = block_data[0];
      double K = block_data[1];
 
      blockData.push_back({c10, K, get_block_ents()});
 
      MOFEM_LOG("EP", sev) << "MatBlock Neo-Hookean c10 = "
                           << blockData.back().c10
                           << " K = " << blockData.back().K << " nb ents. = "
                           << blockData.back().blockEnts.size();
    }
    MoFEMFunctionReturn(0);
  }

getMaterialParameters()

auto EshelbianPlasticity::HMHNeohookean::getMaterialParameters ( EntityHandle  ent )

inline

Definition at line 50 of file HMHNeohookean.cpp.

                                               {
    for (auto &b : blockData) {
      if (b.blockEnts.find(ent) != b.blockEnts.end()) {
        return std::make_pair(b.c10, b.K);
      }
    }
    if (blockData.size() != 0)
      CHK_THROW_MESSAGE(MOFEM_DATA_INCONSISTENCY,
                        "Block not found for entity handle. If you mat set "
                        "block, set it to all elements");
    return std::make_pair(c10_default, K_default);
  }

getOptions()

MoFEMErrorCode EshelbianPlasticity::HMHNeohookean::getOptions ( )

inline

Definition at line 103 of file HMHNeohookean.cpp.

                              {
    MoFEMFunctionBegin;
    CHKERR PetscOptionsBegin(PETSC_COMM_WORLD, "neo_hookean_", "", "none");
 
    c10_default = 1;
    CHKERR PetscOptionsScalar("-c10", "C10", "", c10_default, &c10_default,
                              PETSC_NULL);
    K_default = 1;
    CHKERR PetscOptionsScalar("-K", "Bulk modulus K", "", K_default, &K_default,
                              PETSC_NULL);
 
    alphaGradU = 0;
    CHKERR PetscOptionsScalar("-viscosity_alpha_grad_u", "viscosity", "",
                              alphaGradU, &alphaGradU, PETSC_NULL);
 
    ierr = PetscOptionsEnd();
    CHKERRG(ierr);
 
    MOFEM_LOG_CHANNEL("WORLD");
    MOFEM_TAG_AND_LOG("WORLD", Sev::inform, "MatBlock Neo-Hookean (default)")
        << "c10 = " << c10_default << " K = " << K_default
        << " grad alpha u = " << alphaGradU;
    MoFEMFunctionReturn(0);
  }

recordTape()

MoFEMErrorCode EshelbianPlasticity::HMHNeohookean::recordTape ( const int  tape, DTensor2Ptr *  t_h_ptr  )

inlinevirtual

Implements EshelbianPlasticity::PhysicalEquations.

Definition at line 172 of file HMHNeohookean.cpp.

                                                                  {
    MoFEMFunctionBegin;
    MoFEMFunctionReturn(0);
  }

returnOpJacobian()

virtual EshelbianPlasticity::OpJacobian* EshelbianPlasticity::HMHNeohookean::returnOpJacobian ( const int  tag, const bool  eval_rhs, const bool  eval_lhs, boost::shared_ptr< DataAtIntegrationPts >  data_ptr, boost::shared_ptr< PhysicalEquations >  physics_ptr  )

inlinevirtual

Implements EshelbianPlasticity::PhysicalEquations.

Definition at line 97 of file HMHNeohookean.cpp.

                                                                   {
    return (new OpJacobian(tag, eval_rhs, eval_lhs, data_ptr, physics_ptr));
  }

returnOpSetScale()

VolUserDataOperator* EshelbianPlasticity::HMHNeohookean::returnOpSetScale ( boost::shared_ptr< double >  scale_ptr, boost::shared_ptr< PhysicalEquations >  physics_ptr  )

inlinevirtual

Reimplemented from EshelbianPlasticity::PhysicalEquations.

Definition at line 84 of file HMHNeohookean.cpp.

                                                                   {
 
    return new OpGetScale(scale_ptr, physics_ptr);
  };

returnOpSpatialPhysical()

virtual VolUserDataOperator* EshelbianPlasticity::HMHNeohookean::returnOpSpatialPhysical ( const std::string &  field_name, boost::shared_ptr< DataAtIntegrationPts >  data_ptr, const double  alpha_u  )

inlinevirtual

Reimplemented from EshelbianPlasticity::PhysicalEquations.

Definition at line 190 of file HMHNeohookean.cpp.

                                                {
    return new OpSpatialPhysical(field_name, data_ptr, alpha_u);
  }

returnOpSpatialPhysical_du_du()

VolUserDataOperator* EshelbianPlasticity::HMHNeohookean::returnOpSpatialPhysical_du_du ( std::string  row_field, std::string  col_field, boost::shared_ptr< DataAtIntegrationPts >  data_ptr, const double  alpha  )

inlinevirtual

Reimplemented from EshelbianPlasticity::PhysicalEquations.

Definition at line 206 of file HMHNeohookean.cpp.

                                                                          {
    return new OpSpatialPhysical_du_du(row_field, col_field, data_ptr, alpha);
  }

Member Data Documentation

alphaGradU

double EshelbianPlasticity::HMHNeohookean::alphaGradU

private

Definition at line 217 of file HMHNeohookean.cpp.

blockData

std::vector<BlockData> EshelbianPlasticity::HMHNeohookean::blockData

private

Definition at line 223 of file HMHNeohookean.cpp.

c10_default

double EshelbianPlasticity::HMHNeohookean::c10_default

private

Definition at line 215 of file HMHNeohookean.cpp.

eqScaling

double EshelbianPlasticity::HMHNeohookean::eqScaling = 1.

private

Definition at line 225 of file HMHNeohookean.cpp.

K_default

double EshelbianPlasticity::HMHNeohookean::K_default

private

Definition at line 216 of file HMHNeohookean.cpp.

mField

MoFEM::Interface& EshelbianPlasticity::HMHNeohookean::mField

private

Definition at line 213 of file HMHNeohookean.cpp.

numberOfActiveVariables

constexpr int EshelbianPlasticity::HMHNeohookean::numberOfActiveVariables = 0

staticconstexpr

Definition at line 16 of file HMHNeohookean.cpp.

numberOfDependentVariables

constexpr int EshelbianPlasticity::HMHNeohookean::numberOfDependentVariables = 0

staticconstexpr

Definition at line 17 of file HMHNeohookean.cpp.

---

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

• users_modules/eshelbian_plasticit/src/impl/HMHNeohookean.cpp