v0.14.0
Classes | Public Member Functions | Private Attributes | List of all members
EshelbianPlasticity::HMHHencky Struct Reference
Inheritance diagram for EshelbianPlasticity::HMHHencky:
[legend]
Collaboration diagram for EshelbianPlasticity::HMHHencky:
[legend]

Classes

struct  BlockData
 
struct  OpCalculateEnergy
 Calculate energy density for Hencky material model. More...
 
struct  OpCalculateStretchFromStress
 
struct  OpHenckyJacobian
 
struct  OpSpatialPhysical
 
struct  OpSpatialPhysical_du_du
 

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)
 
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)
 
VolUserDataOperator * returnOpCalculateEnergy (boost::shared_ptr< DataAtIntegrationPts > data_ptr, boost::shared_ptr< double > total_energy_ptr)
 
VolUserDataOperator * returnOpCalculateStretchFromStress (boost::shared_ptr< DataAtIntegrationPts > data_ptr, boost::shared_ptr< PhysicalEquations > physics_ptr)
 
VolUserDataOperator * returnOpCalculateVarStretchFromStress (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)
 
MoFEMErrorCode getInvMatDPtr (boost::shared_ptr< MatrixDouble > mat_inv_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
 
virtual VolUserDataOperator * returnOpSetScale (boost::shared_ptr< double > scale_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 ()
 

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 12 of file HMHHencky.cpp.

Constructor & Destructor Documentation

◆ HMHHencky()

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

Definition at line 14 of file HMHHencky.cpp.

15  : 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 204 of file HMHHencky.cpp.

204  {
205  return extractBlockData(
206 
207  mField.getInterface<MeshsetsManager>()->getCubitMeshsetPtr(std::regex(
208 
209  (boost::format("%s(.*)") % "MAT_ELASTIC").str()
210 
211  )),
212 
213  sev);
214  }

◆ extractBlockData() [2/2]

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

Definition at line 217 of file HMHHencky.cpp.

218  {
219  MoFEMFunctionBegin;
220 
221  for (auto m : meshset_vec_ptr) {
222  MOFEM_TAG_AND_LOG("WORLD", sev, "MatBlock") << *m;
223  std::vector<double> block_data;
224  CHKERR m->getAttributes(block_data);
225  if (block_data.size() != 2) {
226  SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
227  "Expected that block has two attribute");
228  }
229  auto get_block_ents = [&]() {
230  Range ents;
231  CHKERR mField.get_moab().get_entities_by_handle(m->meshset, ents, true);
232  return ents;
233  };
234 
235  double young_modulus = block_data[0];
236  double poisson_ratio = block_data[1];
237  double bulk_modulus_K = young_modulus / (3 * (1 - 2 * poisson_ratio));
238  double shear_modulus_G = young_modulus / (2 * (1 + poisson_ratio));
239 
240  MOFEM_TAG_AND_LOG("WORLD", sev, "MatBlock")
241  << "E = " << young_modulus << " nu = " << poisson_ratio;
242 
243  blockData.push_back({bulk_modulus_K, shear_modulus_G, get_block_ents()});
244  }
245  MOFEM_LOG_CHANNEL("WORLD");
246  MoFEMFunctionReturn(0);
247  }

◆ getInvMatDPtr()

MoFEMErrorCode EshelbianPlasticity::HMHHencky::getInvMatDPtr ( boost::shared_ptr< MatrixDouble >  mat_inv_D_ptr,
double  bulk_modulus_K,
double  shear_modulus_G 
)
inline

[Calculate elasticity tensor]

[Calculate elasticity tensor]

Definition at line 283 of file HMHHencky.cpp.

284  {
285  MoFEMFunctionBegin;
286  //! [Calculate elasticity tensor]
287  auto set_material_compilance = [&]() {
288  FTENSOR_INDEX(SPACE_DIM, i);
289  FTENSOR_INDEX(SPACE_DIM, j);
290  FTENSOR_INDEX(SPACE_DIM, k);
291  FTENSOR_INDEX(SPACE_DIM, l);
292  constexpr auto t_kd = FTensor::Kronecker_Delta_symmetric<int>();
293  const double A = 1. / (2. * shear_modulus_G);
294  const double B =
295  (1. / (9. * bulk_modulus_K)) - (1. / (6. * shear_modulus_G));
296  auto t_inv_D = getFTensor4DdgFromPtr<SPACE_DIM, SPACE_DIM, 0>(
297  &*(mat_inv_D_ptr->data().begin()));
298  t_inv_D(i, j, k, l) =
299  A * ((t_kd(i, k) ^ t_kd(j, l)) / 4.) + B * t_kd(i, j) * t_kd(k, l);
300  };
301  //! [Calculate elasticity tensor]
302  constexpr auto size_symm = (SPACE_DIM * (SPACE_DIM + 1)) / 2;
303  mat_inv_D_ptr->resize(size_symm, size_symm, false);
304  set_material_compilance();
305  MoFEMFunctionReturn(0);
306  }

◆ 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 249 of file HMHHencky.cpp.

252  {
253  MoFEMFunctionBegin;
254  //! [Calculate elasticity tensor]
255  auto set_material_stiffness = [&]() {
256  FTENSOR_INDEX(SPACE_DIM, i);
257  FTENSOR_INDEX(SPACE_DIM, j);
258  FTENSOR_INDEX(SPACE_DIM, k);
259  FTENSOR_INDEX(SPACE_DIM, l);
260  constexpr auto t_kd = FTensor::Kronecker_Delta_symmetric<int>();
261  auto t_D = getFTensor4DdgFromPtr<SPACE_DIM, SPACE_DIM, 0>(
262  &*(mat_D_ptr->data().begin()));
263  auto t_axiator_D = getFTensor4DdgFromPtr<SPACE_DIM, SPACE_DIM, 0>(
264  &*mat_axiator_D_ptr->data().begin());
265  auto t_deviator_D = getFTensor4DdgFromPtr<SPACE_DIM, SPACE_DIM, 0>(
266  &*mat_deviator_D_ptr->data().begin());
267  t_axiator_D(i, j, k, l) = (bulk_modulus_K - (2. / 3.) * shear_modulus_G) *
268  t_kd(i, j) * t_kd(k, l);
269  t_deviator_D(i, j, k, l) =
270  2 * shear_modulus_G * ((t_kd(i, k) ^ t_kd(j, l)) / 4.);
271  t_D(i, j, k, l) = t_axiator_D(i, j, k, l) + t_deviator_D(i, j, k, l);
272  };
273  //! [Calculate elasticity tensor]
274  constexpr auto size_symm = (SPACE_DIM * (SPACE_DIM + 1)) / 2;
275  mat_D_ptr->resize(size_symm, size_symm, false);
276  mat_axiator_D_ptr->resize(size_symm, size_symm, false);
277  mat_deviator_D_ptr->resize(size_symm, size_symm, false);
278  set_material_stiffness();
279  MoFEMFunctionReturn(0);
280  }

◆ getOptions()

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

Definition at line 186 of file HMHHencky.cpp.

186  {
187  MoFEMFunctionBegin;
188  CHKERR PetscOptionsBegin(PETSC_COMM_WORLD, "hencky_", "", "none");
189 
190  CHKERR PetscOptionsScalar("-young_modulus", "Young modulus", "", E, &E,
191  PETSC_NULL);
192  CHKERR PetscOptionsScalar("-poisson_ratio", "poisson ratio", "", nu, &nu,
193  PETSC_NULL);
194 
195  MOFEM_LOG_CHANNEL("WORLD");
196  MOFEM_LOG("WORLD", Sev::verbose) << "Hencky: E = " << E << " nu = " << nu;
197 
198  ierr = PetscOptionsEnd();
199  CHKERRG(ierr);
200 
201  MoFEMFunctionReturn(0);
202  }

◆ recordTape()

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

Implements EshelbianPlasticity::PhysicalEquations.

Definition at line 17 of file HMHHencky.cpp.

17 { return 0; }

◆ returnOpCalculateEnergy()

VolUserDataOperator* EshelbianPlasticity::HMHHencky::returnOpCalculateEnergy ( boost::shared_ptr< DataAtIntegrationPts >  data_ptr,
boost::shared_ptr< double >  total_energy_ptr 
)
inlinevirtual

Reimplemented from EshelbianPlasticity::PhysicalEquations.

Definition at line 148 of file HMHHencky.cpp.

149  {
150  return new OpCalculateEnergy(data_ptr, total_energy_ptr);
151  }

◆ returnOpCalculateStretchFromStress()

VolUserDataOperator* EshelbianPlasticity::HMHHencky::returnOpCalculateStretchFromStress ( boost::shared_ptr< DataAtIntegrationPts >  data_ptr,
boost::shared_ptr< PhysicalEquations >  physics_ptr 
)
inlinevirtual

Reimplemented from EshelbianPlasticity::PhysicalEquations.

Definition at line 169 of file HMHHencky.cpp.

171  {
172  return new OpCalculateStretchFromStress(
173  data_ptr, data_ptr->getLogStretchTensorAtPts(),
174  data_ptr->getApproxPAtPts(),
175  boost::dynamic_pointer_cast<HMHHencky>(physics_ptr));
176  }

◆ returnOpCalculateVarStretchFromStress()

VolUserDataOperator* EshelbianPlasticity::HMHHencky::returnOpCalculateVarStretchFromStress ( boost::shared_ptr< DataAtIntegrationPts >  data_ptr,
boost::shared_ptr< PhysicalEquations >  physics_ptr 
)
inlinevirtual

Reimplemented from EshelbianPlasticity::PhysicalEquations.

Definition at line 178 of file HMHHencky.cpp.

180  {
181  return new OpCalculateStretchFromStress(
182  data_ptr, data_ptr->getVarLogStreachPts(), data_ptr->getVarPiolaPts(),
183  boost::dynamic_pointer_cast<HMHHencky>(physics_ptr));
184  }

◆ 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 70 of file HMHHencky.cpp.

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

◆ returnOpSpatialPhysical()

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

Reimplemented from EshelbianPlasticity::PhysicalEquations.

Definition at line 95 of file HMHHencky.cpp.

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

◆ returnOpSpatialPhysical_du_du()

VolUserDataOperator* EshelbianPlasticity::HMHHencky::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 117 of file HMHHencky.cpp.

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

Member Data Documentation

◆ blockData

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

Definition at line 316 of file HMHHencky.cpp.

◆ E

double EshelbianPlasticity::HMHHencky::E
private

Definition at line 318 of file HMHHencky.cpp.

◆ mField

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

Definition at line 309 of file HMHHencky.cpp.

◆ nu

double EshelbianPlasticity::HMHHencky::nu
private

Definition at line 319 of file HMHHencky.cpp.

The documentation for this struct was generated from the following file:
MoFEM::UnknownInterface::getInterface
MoFEMErrorCode getInterface(IFACE *&iface) const
Get interface reference to pointer of interface.
Definition: UnknownInterface.hpp:93
EshelbianPlasticity::size_symm
constexpr static auto size_symm
Definition: EshelbianAux.hpp:41
MOFEM_LOG_CHANNEL
#define MOFEM_LOG_CHANNEL(channel)
Set and reset channel.
Definition: LogManager.hpp:284
young_modulus
double young_modulus
Young modulus.
Definition: plastic.cpp:125
E
A
constexpr AssemblyType A
Definition: operators_tests.cpp:30
bulk_modulus_K
double bulk_modulus_K
Definition: dynamic_first_order_con_law.cpp:95
EshelbianPlasticity::HMHHencky::blockData
std::vector< BlockData > blockData
Definition: HMHHencky.cpp:316
FTENSOR_INDEX
#define FTENSOR_INDEX(DIM, I)
Definition: Templates.hpp:2013
EshelbianPlasticity::HMHHencky::nu
double nu
Definition: HMHHencky.cpp:319
CHKERR
#define CHKERR
Inline error check.
Definition: definitions.h:548
MoFEM::CoreInterface::get_moab
virtual moab::Interface & get_moab()=0
SPACE_DIM
constexpr int SPACE_DIM
Definition: child_and_parent.cpp:16
OpCalculateEnergy
Definition: HookeElement.hpp:164
poisson_ratio
double poisson_ratio
Poisson ratio.
Definition: plastic.cpp:126
EshelbianPlasticity::HMHHencky::E
double E
Definition: HMHHencky.cpp:318
EshelbianPlasticity::HMHHencky::mField
MoFEM::Interface & mField
Definition: HMHHencky.cpp:309
i
FTensor::Index< 'i', SPACE_DIM > i
Definition: hcurl_divergence_operator_2d.cpp:27
t_kd
constexpr auto t_kd
Definition: free_surface.cpp:139
EshelbianPlasticity::HMHHencky::extractBlockData
MoFEMErrorCode extractBlockData(Sev sev)
Definition: HMHHencky.cpp:204
field_name
constexpr auto field_name
Definition: poisson_2d_homogeneous.cpp:13
Range
MOFEM_TAG_AND_LOG
#define MOFEM_TAG_AND_LOG(channel, severity, tag)
Tag and log in channel.
Definition: LogManager.hpp:362
MOFEM_LOG
#define MOFEM_LOG(channel, severity)
Log.
Definition: LogManager.hpp:308
shear_modulus_G
double shear_modulus_G
Definition: dynamic_first_order_con_law.cpp:96
j
FTensor::Index< 'j', 3 > j
Definition: matrix_function.cpp:19
MoFEM::Exceptions::ierr
static MoFEMErrorCodeGeneric< PetscErrorCode > ierr
Definition: Exceptions.hpp:76
MOFEM_DATA_INCONSISTENCY
@ MOFEM_DATA_INCONSISTENCY
Definition: definitions.h:31
m
FTensor::Index< 'm', 3 > m
Definition: shallow_wave.cpp:80
FTensor::Kronecker_Delta_symmetric
Kronecker Delta class symmetric.
Definition: Kronecker_Delta.hpp:49
k
FTensor::Index< 'k', 3 > k
Definition: matrix_function.cpp:20
MoFEMFunctionReturn
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:429
CHKERRG
#define CHKERRG(n)
Check error code of MoFEM/MOAB/PETSc function.
Definition: definitions.h:496
EshelbianPlasticity::PhysicalEquations::PhysicalEquations
PhysicalEquations()=delete
MoFEMFunctionBegin
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:359
l
FTensor::Index< 'l', 3 > l
Definition: matrix_function.cpp:21
Generated by Doxygen 1.8.17 and hosted at University of Glasgow