EshelbianPlasticity::HMHHencky Struct Reference

Inheritance diagram for EshelbianPlasticity::HMHHencky:

Collaboration diagram for EshelbianPlasticity::HMHHencky:

Classes
struct	BlockData
struct	OpHenckyJacobian

Public Member Functions
	HMHHencky (MoFEM::Interface &m_field, const double E, const double nu)
MoFEMErrorCode	recordTape (const int tag, DTensor2Ptr *t_h)
virtual 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 (boost::shared_ptr< DataAtIntegrationPts > data_ptr)
MoFEMErrorCode	extractBlockData (Sev sev)
MoFEMErrorCode	extractBlockData (std::vector< const CubitMeshSets * > meshset_vec_ptr, Sev sev)
MoFEMErrorCode	getMatDPtr (boost::shared_ptr< MatrixDouble > mat_D_ptr, boost::shared_ptr< MatrixDouble > mat_axiator_D_ptr, boost::shared_ptr< MatrixDouble > mat_deviator_D_ptr, double bulk_modulus_K, double shear_modulus_G)
Public Member Functions inherited from EshelbianPlasticity::PhysicalEquations
	PhysicalEquations ()=delete
	PhysicalEquations (const int size_active, const int size_dependent)
virtual	~PhysicalEquations ()=default
DTensor2Ptr	get_P ()
DTensor3Ptr	get_P_dh0 ()
DTensor3Ptr	get_P_dh1 ()
DTensor3Ptr	get_P_dh2 ()
DTensor2Ptr	get_h ()

Private Attributes
MoFEM::Interface &	mField
std::vector< BlockData >	blockData
double	E
double	nu

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 465 of file EshelbianADOL-C.cpp.

Constructor & Destructor Documentation

HMHHencky()

EshelbianPlasticity::HMHHencky::HMHHencky ( MoFEM::Interface &  m_field, const double  E, const double  nu  )

inline

Definition at line 467 of file EshelbianADOL-C.cpp.

      : PhysicalEquations(0, 0), mField(m_field), E(E), nu(nu) {}

Member Function Documentation

extractBlockData() [1/2]

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

inline

Definition at line 547 of file EshelbianADOL-C.cpp.

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

extractBlockData() [2/2]

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

inline

Definition at line 560 of file EshelbianADOL-C.cpp.

                            {
    MoFEMFunctionBegin;
 
    for (auto m : meshset_vec_ptr) {
      MOFEM_TAG_AND_LOG("WORLD", sev, "MatBlock") << *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 young_modulus = block_data[0];
      double poisson_ratio = block_data[1];
      double bulk_modulus_K = young_modulus / (3 * (1 - 2 * poisson_ratio));
      double shear_modulus_G = young_modulus / (2 * (1 + poisson_ratio));
 
      MOFEM_TAG_AND_LOG("WORLD", sev, "MatBlock")
          << "E = " << young_modulus << " nu = " << poisson_ratio;
 
      blockData.push_back({bulk_modulus_K, shear_modulus_G, get_block_ents()});
    }
    MOFEM_LOG_CHANNEL("WORLD");
    MoFEMFunctionReturn(0);
  }

getMatDPtr()

MoFEMErrorCode EshelbianPlasticity::HMHHencky::getMatDPtr ( boost::shared_ptr< MatrixDouble >  mat_D_ptr, boost::shared_ptr< MatrixDouble >  mat_axiator_D_ptr, boost::shared_ptr< MatrixDouble >  mat_deviator_D_ptr, double  bulk_modulus_K, double  shear_modulus_G  )

inline

[Calculate elasticity tensor]

[Calculate elasticity tensor]

Definition at line 592 of file EshelbianADOL-C.cpp.

                                                                           {
    MoFEMFunctionBegin;
    //! [Calculate elasticity tensor]
    auto set_material_stiffness = [&]() {
      FTensor::Index<'i', SPACE_DIM> i;
      FTensor::Index<'j', SPACE_DIM> j;
      FTensor::Index<'k', SPACE_DIM> k;
      FTensor::Index<'l', SPACE_DIM> l;
      constexpr auto t_kd = FTensor::Kronecker_Delta_symmetric<int>();
      double A = (SPACE_DIM == 2)
                     ? 2 * shear_modulus_G /
                           (bulk_modulus_K + (4. / 3.) * shear_modulus_G)
                     : 1;
      auto t_D = getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM, 0>(*mat_D_ptr);
      auto t_axiator_D =
          getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM, 0>(*mat_axiator_D_ptr);
      auto t_deviator_D =
          getFTensor4DdgFromMat<SPACE_DIM, SPACE_DIM, 0>(*mat_deviator_D_ptr);
      t_axiator_D(i, j, k, l) = A *
                                (bulk_modulus_K - (2. / 3.) * shear_modulus_G) *
                                t_kd(i, j) * t_kd(k, l);
      t_deviator_D(i, j, k, l) =
          2 * shear_modulus_G * ((t_kd(i, k) ^ t_kd(j, l)) / 4.);
      t_D(i, j, k, l) = t_axiator_D(i, j, k, l) + t_deviator_D(i, j, k, l);
    };
    //! [Calculate elasticity tensor]
    constexpr auto size_symm = (SPACE_DIM * (SPACE_DIM + 1)) / 2;
    mat_D_ptr->resize(size_symm * size_symm, 1);
    mat_axiator_D_ptr->resize(size_symm * size_symm, 1);
    mat_deviator_D_ptr->resize(size_symm * size_symm, 1);
    set_material_stiffness();
    MoFEMFunctionReturn(0);
  }

getOptions()

MoFEMErrorCode EshelbianPlasticity::HMHHencky::getOptions ( boost::shared_ptr< DataAtIntegrationPts >  data_ptr )

inline

Definition at line 532 of file EshelbianADOL-C.cpp.

                                                                            {
    MoFEMFunctionBegin;
    CHKERR PetscOptionsBegin(PETSC_COMM_WORLD, "hencky_", "", "none");
 
    CHKERR PetscOptionsScalar("-young_modulus", "Young modulus", "", E, &E,
                              PETSC_NULL);
    CHKERR PetscOptionsScalar("-poisson_ratio", "poisson ratio", "", nu, &nu,
                              PETSC_NULL);
 
    ierr = PetscOptionsEnd();
    CHKERRG(ierr);
 
    MoFEMFunctionReturn(0);
  }

recordTape()

MoFEMErrorCode EshelbianPlasticity::HMHHencky::recordTape ( const int  tag, DTensor2Ptr *  t_h  )

inlinevirtual

Implements EshelbianPlasticity::PhysicalEquations.

Definition at line 470 of file EshelbianADOL-C.cpp.

{ return 0; }

returnOpJacobian()

virtual OpJacobian* EshelbianPlasticity::HMHHencky::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 525 of file EshelbianADOL-C.cpp.

                                                                    {
    return (new OpHenckyJacobian(
        data_ptr, boost::dynamic_pointer_cast<HMHHencky>(physics_ptr)));
  }

Member Data Documentation

blockData

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

private

Definition at line 637 of file EshelbianADOL-C.cpp.

E

double EshelbianPlasticity::HMHHencky::E

private

Definition at line 639 of file EshelbianADOL-C.cpp.

mField

MoFEM::Interface& EshelbianPlasticity::HMHHencky::mField

private

Definition at line 630 of file EshelbianADOL-C.cpp.

nu

double EshelbianPlasticity::HMHHencky::nu

private

Definition at line 640 of file EshelbianADOL-C.cpp.

---

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

• users_modules/eshelbian_plasticit/src/impl/EshelbianADOL-C.cpp