Inheritance diagram for EshelbianPlasticity::HMHNeohookean:

Collaboration diagram for EshelbianPlasticity::HMHNeohookean:

Classes
struct	BlockData

struct	OpGetScale

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)

MoFEMErrorCode	setParams (const int tag, EntityHandle ent, int gg)

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

UserDataOperator *	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)

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 ()=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 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.

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)

Static Public Attributes
static constexpr int	numberOfActiveVariables = 9

static constexpr int	numberOfDependentVariables = 9

Private Attributes
MoFEM::Interface &	mField

double	c10_default

double	K_default

double	alphaGradU

std::vector< BlockData >	blockData

double	eqScaling = 1.

PetscBool	adolCOn = PETSC_FALSE

ATensor2	th

ATensor2	tF

adouble	detF

ATensor2	tInvF

ATensor2	tCof

ATensor2	tP

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

Public Attributes inherited from EshelbianPlasticity::PhysicalEquations
std::vector< double >	activeVariables

std::vector< double >	dependentVariablesPiola

std::vector< double >	dependentVariablesPiolaDirevatives

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 71 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 (adolCOn == PETSC_FALSE) {
      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 201 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 214 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 104 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 170 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);
 
    CHKERR PetscOptionsBool("-adolc_on", "adolc ON", "", adolCOn, &adolCOn,
                            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;
    MOFEM_TAG_AND_LOG("WORLD", Sev::inform, "MatBlock Neo-Hookean (default)")
            << "ADOL-C On " << (adolCOn)
        ? "Yes"
        : "No";
    MoFEMFunctionReturn(0);
  }

◆ recordTape()

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

inlinevirtual

Implements EshelbianPlasticity::PhysicalEquations.

Definition at line 245 of file HMHNeohookean.cpp.

                                                                  {
    MoFEMFunctionBegin;
 
    FTENSOR_INDEX(3, i);
    FTENSOR_INDEX(3, j);
    FTENSOR_INDEX(3, I);
    FTENSOR_INDEX(3, J);
 
    auto large = [&]() {
      MoFEMFunctionBegin;
 
      constexpr auto t_kd = FTensor::Kronecker_Delta<int>();
 
      auto ih = get_h();
      if (t_h_ptr)
        ih(i, j) = (*t_h_ptr)(i, j);
      else {
        ih(i, j) = t_kd(i, j);
      }
 
      auto r_P = get_P();
 
      enableMinMaxUsingAbs();
      trace_on(tape);
 
      // Set active variables to ADOL-C
      th(i, j) <<= get_h()(i, j);
 
      tF(i, I) = th(i, I);
      CHKERR determinantTensor3by3(tF, detF);
      CHKERR invertTensor3by3(tF, detF, tInvF);
 
      // tCof(i, I) = detF * tInvF(I, i);
 
      // Psi = c10* (trace(tCof) + 2ln(det(F)) - 3) + (K/2.) * (log(detF))^2;
 
      auto p_c10 = mkparam(c10_default);
      auto p_K = mkparam(K_default);
      tP(i, I) =
          2. * p_c10 * (tF(i, I) - tInvF(i, I)) + p_K * log(detF) * tInvF(i, I);
 
      // Set dependent variables to ADOL-C
      tP(i, j) >>= r_P(i, j);
 
      trace_off();
 
      MoFEMFunctionReturn(0);
    };
 
    switch (EshelbianCore::gradApproximator) {
    case NO_H1_CONFIGURATION:
    case LARGE_ROT:
    case MODERATE_ROT:
      CHKERR large();
      break;
    default:
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "gradApproximator not handled");
    };
 
    MoFEMFunctionReturn(0);
  }

◆ returnOpJacobian()

UserDataOperator * 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

Reimplemented from EshelbianPlasticity::PhysicalEquations.

Definition at line 159 of file HMHNeohookean.cpp.

                                                                   {
    if (adolCOn) {
      return PhysicalEquations::returnOpJacobian(tag, eval_rhs, eval_lhs,
                                                 data_ptr, physics_ptr);
    } else {
      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 146 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 321 of file HMHNeohookean.cpp.

                                                {
    if (adolCOn) {
      return PhysicalEquations::returnOpSpatialPhysical(field_name, data_ptr,
                                                        alpha_u);
    } else {
      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 366 of file HMHNeohookean.cpp.

                                                                          {
    if (adolCOn) {
      return PhysicalEquations::returnOpSpatialPhysical_du_du(
          row_field, col_field, data_ptr, alpha);
    } else {
      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 347 of file HMHNeohookean.cpp.

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

◆ setParams()

MoFEMErrorCode EshelbianPlasticity::HMHNeohookean::setParams	(	const int	tag,
		EntityHandle	ent,
		int	gg
	)

inlinevirtual

Reimplemented from EshelbianPlasticity::PhysicalEquations.

Definition at line 117 of file HMHNeohookean.cpp.

                                                                    {
    MoFEMFunctionBegin;
    auto [def_c10, def_K] = getMaterialParameters(ent);
    std::vector<double> params = {def_c10, def_K};
    set_param_vec(tag, params.size(), params.data());
    MoFEMFunctionReturn(0);
  };

Member Data Documentation

◆ adolCOn

PetscBool EshelbianPlasticity::HMHNeohookean::adolCOn = PETSC_FALSE

private

Definition at line 393 of file HMHNeohookean.cpp.

◆ alphaGradU

double EshelbianPlasticity::HMHNeohookean::alphaGradU

private

Definition at line 382 of file HMHNeohookean.cpp.

◆ blockData

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

private

Definition at line 388 of file HMHNeohookean.cpp.

◆ c10_default

double EshelbianPlasticity::HMHNeohookean::c10_default

private

Definition at line 380 of file HMHNeohookean.cpp.

◆ detF

adouble EshelbianPlasticity::HMHNeohookean::detF

private

Definition at line 397 of file HMHNeohookean.cpp.

◆ eqScaling

double EshelbianPlasticity::HMHNeohookean::eqScaling = 1.

private

Definition at line 390 of file HMHNeohookean.cpp.

◆ K_default

double EshelbianPlasticity::HMHNeohookean::K_default

private

Definition at line 381 of file HMHNeohookean.cpp.

◆ mField

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

private

Definition at line 378 of file HMHNeohookean.cpp.

◆ numberOfActiveVariables

constexpr int EshelbianPlasticity::HMHNeohookean::numberOfActiveVariables = 9

staticconstexpr

Definition at line 68 of file HMHNeohookean.cpp.

◆ numberOfDependentVariables

constexpr int EshelbianPlasticity::HMHNeohookean::numberOfDependentVariables = 9

staticconstexpr

Definition at line 69 of file HMHNeohookean.cpp.

◆ tCof

ATensor2 EshelbianPlasticity::HMHNeohookean::tCof

private

Definition at line 399 of file HMHNeohookean.cpp.

◆ tF

ATensor2 EshelbianPlasticity::HMHNeohookean::tF

private

Definition at line 396 of file HMHNeohookean.cpp.

◆ th

ATensor2 EshelbianPlasticity::HMHNeohookean::th

private

Definition at line 395 of file HMHNeohookean.cpp.

◆ tInvF

ATensor2 EshelbianPlasticity::HMHNeohookean::tInvF

private

Definition at line 398 of file HMHNeohookean.cpp.

◆ tP

ATensor2 EshelbianPlasticity::HMHNeohookean::tP

private

Definition at line 400 of file HMHNeohookean.cpp.

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

users_modules/eshelbian_plasticity/src/impl/HMHNeohookean.cpp

Classes

Public Member Functions

Static Public Member Functions

Static Public Attributes

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ HMHNeohookean()

Member Function Documentation

◆ extractBlockData() [1/2]

◆ extractBlockData() [2/2]

◆ fun_d_neohookean()

◆ fun_diff_neohookean_bulk()

◆ fun_neohookean()

◆ fun_neohookean_bulk()

◆ getMaterialParameters()

◆ getOptions()

◆ recordTape()

◆ returnOpJacobian()

◆ returnOpSetScale()

◆ returnOpSpatialPhysical()

◆ returnOpSpatialPhysical_du_du()

◆ returnOpSpatialPhysicalExternalStrain()

◆ setParams()

Member Data Documentation

◆ adolCOn

◆ alphaGradU

◆ blockData

◆ c10_default

◆ detF

◆ eqScaling

◆ K_default

◆ mField

◆ numberOfActiveVariables

◆ numberOfDependentVariables

◆ tCof

◆ tF

◆ th

◆ tInvF

◆ tP