EshelbianPlasticity::HMHNeohookean Struct Reference

Inheritance diagram for EshelbianPlasticity::HMHNeohookean:

Collaboration diagram for EshelbianPlasticity::HMHNeohookean:

Classes
struct	BlockData
struct	OpJacobian
struct	OpSpatialPhysical
struct	OpSpatialPhysical_du_du
struct	OpSpatialPhysicalExternalStrain

Public Member Functions
	HMHNeohookean (MoFEM::Interface &m_field, const double c10, const double K)
auto	getMaterialParameters (EntityHandle ent)
UserDataOperator *	returnOpJacobian (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)
virtual VolUserDataOperator *	returnOpSpatialPhysical (const std::string &field_name, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const double alpha_u)
virtual VolUserDataOperator *	returnOpSpatialPhysicalExternalStrain (const std::string &field_name, boost::shared_ptr< DataAtIntegrationPts > data_ptr, boost::shared_ptr< ExternalStrainVec > external_strain_vec_ptr, std::map< std::string, boost::shared_ptr< ScalingMethod > > smv)
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 ()=default
virtual	~PhysicalEquations ()=default
virtual VolUserDataOperator *	returnOpTopoSpatialPhysical (const std::string &field_name, boost::shared_ptr< DataAtIntegrationPts > data_ptr, SmartPetscObj< Vec > assemble_vec, boost::shared_ptr< TopologicalData > topo_ptr, const double alpha_u, boost::shared_ptr< double > J_ptr)
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)

Static Public Member Functions
static double	fun_neohookean (double c10, double v)
	Definition of Neo-hookean function.
static double	fun_d_neohookean (double c10, double v)
	Definition of derivative of Neo-hookean function.
static double	fun_neohookean_bulk (double K, double tr)
	Definition of axiator of Neo-hookean function.
static double	fun_diff_neohookean_bulk (double K, double tr)
	Definition of derivative of axiator of Neo-hookean function.

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

Detailed Description

Definition at line 13 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 68 of file HMHNeohookean.cpp.

      : PhysicalEquations(), 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 149 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 162 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 atleast two attributes");
      }
      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);
  }

fun_d_neohookean()

static double EshelbianPlasticity::HMHNeohookean::fun_d_neohookean ( double  c10, double  v  )

inlinestatic

Definition of derivative of Neo-hookean function.

Parameters

c10
v

Returns

* double

Definition at line 38 of file HMHNeohookean.cpp.

                                                              {
    return
 
        (4.0 * c10) * EshelbianCore::d_f(2 * v);
  }

fun_diff_neohookean_bulk()

static double EshelbianPlasticity::HMHNeohookean::fun_diff_neohookean_bulk ( double  K, double  tr  )

inlinestatic

Definition of derivative of axiator of Neo-hookean function.

Parameters

K
tr

Returns

double

Definition at line 64 of file HMHNeohookean.cpp.

                                                                     {
    return K;//K * exp(tr) * (1. + tr); // K*2*ln(J)/J
  }

fun_neohookean()

static double EshelbianPlasticity::HMHNeohookean::fun_neohookean ( double  c10, double  v  )

inlinestatic

Definition of Neo-hookean function.

P(i, I) = 2. * c10_default * (tF(i, I) - tInvF(i, I)) + K_default * log(detF) * tInvF(i, I) tCof(i, I) = detF * tInvF(I, i) Psi = c10* (trace(tCof) + 2ln(det(F)) - 3) + (K/2.) * (log(detF))^2

Parameters

c10
v

Returns

double

Definition at line 27 of file HMHNeohookean.cpp.

                                                            {
    return (2.0 * c10) * (-1. + EshelbianCore::f(2 * v));
  }

fun_neohookean_bulk()

static double EshelbianPlasticity::HMHNeohookean::fun_neohookean_bulk ( double  K, double  tr  )

inlinestatic

Definition of axiator of Neo-hookean function.

Psi = (K/2.) * (log(J))^2 = (K/2.) * (log(exp tr H))^2 = (K/2.) * (tr H)^2

Parameters

c10
v

Returns

double

Definition at line 53 of file HMHNeohookean.cpp.

                                                                {
    return K * tr; 
  }

getMaterialParameters()

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

inline

Definition at line 99 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 125 of file HMHNeohookean.cpp.

                              {
    MoFEMFunctionBegin;
    PetscOptionsBegin(PETSC_COMM_WORLD, "neo_hookean_", "", "none");
 
    c10_default = 1;
    CHKERR PetscOptionsScalar("-c10", "C10", "", c10_default, &c10_default,
                              PETSC_NULLPTR);
    K_default = 1;
    CHKERR PetscOptionsScalar("-K", "Bulk modulus K", "", K_default, &K_default,
                              PETSC_NULLPTR);
 
    alphaGradU = 0;
    CHKERR PetscOptionsScalar("-viscosity_alpha_grad_u", "viscosity", "",
                              alphaGradU, &alphaGradU, PETSC_NULLPTR);
 
    PetscOptionsEnd();
 
    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);
  }

returnOpJacobian()

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

inlinevirtual

Reimplemented from EshelbianPlasticity::PhysicalEquations.

Definition at line 119 of file HMHNeohookean.cpp.

                                                                   {
    return (new OpJacobian(eval_rhs, eval_lhs, data_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 206 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 248 of file HMHNeohookean.cpp.

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

returnOpSpatialPhysicalExternalStrain()

virtual VolUserDataOperator * EshelbianPlasticity::HMHNeohookean::returnOpSpatialPhysicalExternalStrain ( const std::string &  field_name, boost::shared_ptr< DataAtIntegrationPts >  data_ptr, boost::shared_ptr< ExternalStrainVec >  external_strain_vec_ptr, std::map< std::string, boost::shared_ptr< ScalingMethod > >  smv  )

inlinevirtual

Reimplemented from EshelbianPlasticity::PhysicalEquations.

Definition at line 227 of file HMHNeohookean.cpp.

                                                             {
    return new OpSpatialPhysicalExternalStrain(field_name, data_ptr,
                                               external_strain_vec_ptr, smv);
  }

Member Data Documentation

alphaGradU

double EshelbianPlasticity::HMHNeohookean::alphaGradU

private

Definition at line 259 of file HMHNeohookean.cpp.

blockData

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

private

Definition at line 265 of file HMHNeohookean.cpp.

c10_default

double EshelbianPlasticity::HMHNeohookean::c10_default

private

Definition at line 257 of file HMHNeohookean.cpp.

eqScaling

double EshelbianPlasticity::HMHNeohookean::eqScaling = 1.

private

Definition at line 267 of file HMHNeohookean.cpp.

K_default

double EshelbianPlasticity::HMHNeohookean::K_default

private

Definition at line 258 of file HMHNeohookean.cpp.

mField

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

private

Definition at line 255 of file HMHNeohookean.cpp.

---

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

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