|
| v0.14.0
|
|
| OpSpatialPhysical (const std::string &field_name, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const double alpha_u) |
|
MoFEMErrorCode | integrate (EntData &data) |
|
MoFEMErrorCode | integrateHencky (EntData &data) |
|
MoFEMErrorCode | integratePolyconvexHencky (EntData &data) |
|
| OpAssembleVolume (const std::string &field, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const char type) |
|
| OpAssembleVolume (std::string row_field, std::string col_field, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const char type, const bool assemble_symmetry) |
|
| OpAssembleBasic (const std::string &field_name, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const char type) |
|
| OpAssembleBasic (std::string row_field, std::string col_field, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const char type, const bool assemble_symmetry, ScaleOff scale_off=[]() { return 1;}) |
|
virtual MoFEMErrorCode | integrate (int row_side, EntityType row_type, EntData &data) |
|
virtual MoFEMErrorCode | integrate (EntData &row_data, EntData &col_data) |
|
virtual MoFEMErrorCode | assemble (EntData &data) |
|
virtual MoFEMErrorCode | assemble (int row_side, EntityType row_type, EntData &data) |
|
virtual MoFEMErrorCode | assemble (int row_side, int col_side, EntityType row_type, EntityType col_type, EntData &row_data, EntData &col_data) |
|
MoFEMErrorCode | doWork (int side, EntityType type, EntData &data) |
|
MoFEMErrorCode | doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntData &row_data, EntData &col_data) |
|
Definition at line 79 of file Hencky.cpp.
◆ OpSpatialPhysical()
EshelbianPlasticity::HMHHencky::OpSpatialPhysical::OpSpatialPhysical |
( |
const std::string & |
field_name, |
|
|
boost::shared_ptr< DataAtIntegrationPts > |
data_ptr, |
|
|
const double |
alpha_u |
|
) |
| |
Definition at line 287 of file Hencky.cpp.
294 "get polyconvex option failed");
◆ integrate()
MoFEMErrorCode EshelbianPlasticity::HMHHencky::OpSpatialPhysical::integrate |
( |
EntData & |
data | ) |
|
|
virtual |
◆ integrateHencky()
MoFEMErrorCode EshelbianPlasticity::HMHHencky::OpSpatialPhysical::integrateHencky |
( |
EntData & |
data | ) |
|
Definition at line 307 of file Hencky.cpp.
314 int nb_integration_pts = data.
getN().size1();
315 auto v = getVolume();
316 auto t_w = getFTensor0IntegrationWeight();
317 auto t_approx_P_adjont_log_du =
318 getFTensor1FromMat<size_symm>(
dataAtPts->adjointPdUAtPts);
319 auto t_log_stretch_h1 =
320 getFTensor2SymmetricFromMat<3>(
dataAtPts->logStretchTotalTensorAtPts);
322 getFTensor2SymmetricFromMat<3>(
dataAtPts->logStretchDotTensorAtPts);
324 auto t_D = getFTensor4DdgFromPtr<3, 3, 0>(&*
dataAtPts->matD.data().begin());
330 auto get_ftensor2 = [](
auto &
v) {
332 &
v[0], &
v[1], &
v[2], &
v[3], &
v[4], &
v[5]);
335 int nb_base_functions = data.
getN().size2();
337 for (
int gg = 0; gg != nb_integration_pts; ++gg) {
339 auto t_nf = get_ftensor2(
nF);
343 t_D(
i,
j,
k,
l) * (t_log_stretch_h1(
k,
l) +
alphaU * t_dot_log_u(
k,
l));
346 a * (t_approx_P_adjont_log_du(
L) - t_L(
i,
j,
L) * t_T(
i,
j));
349 for (; bb != nb_dofs / 6; ++bb) {
350 t_nf(
L) += t_row_base_fun * t_residual(
L);
354 for (; bb != nb_base_functions; ++bb)
358 ++t_approx_P_adjont_log_du;
◆ integratePolyconvexHencky()
MoFEMErrorCode EshelbianPlasticity::HMHHencky::OpSpatialPhysical::integratePolyconvexHencky |
( |
EntData & |
data | ) |
|
Definition at line 366 of file Hencky.cpp.
373 int nb_integration_pts = data.
getN().size1();
374 auto v = getVolume();
375 auto t_w = getFTensor0IntegrationWeight();
376 auto t_approx_P_adjont_log_du =
377 getFTensor1FromMat<size_symm>(
dataAtPts->adjointPdUAtPts);
378 auto t_log_stretch_h1 =
379 getFTensor2SymmetricFromMat<3>(
dataAtPts->logStretchTotalTensorAtPts);
381 getFTensor2SymmetricFromMat<3>(
dataAtPts->logStretchDotTensorAtPts);
383 auto t_D = getFTensor4DdgFromPtr<3, 3, 0>(&*
dataAtPts->matD.data().begin());
389 auto get_ftensor2 = [](
auto &
v) {
391 &
v[0], &
v[1], &
v[2], &
v[3], &
v[4], &
v[5]);
394 constexpr
double nohat_k = 1. / 4;
395 constexpr
double hat_k = 1. / 8;
399 constexpr
double third = boost::math::constants::third<double>();
403 int nb_base_functions = data.
getN().size2();
405 for (
int gg = 0; gg != nb_integration_pts; ++gg) {
407 auto t_nf = get_ftensor2(
nF);
409 double log_det = t_log_stretch_h1(
i,
i);
410 double log_det2 = log_det * log_det;
413 double dev_norm2 = t_dev(
i,
j) * t_dev(
i,
j);
416 auto A = 2 *
mu * std::exp(nohat_k * dev_norm2);
417 auto B =
lambda * std::exp(hat_k * log_det2) * log_det;
420 A * (t_dev(
k,
l) * t_diff_deviator(
k,
l,
i,
j))
432 a * (t_approx_P_adjont_log_du(
L) - t_L(
i,
j,
L) * t_T(
i,
j));
435 for (; bb != nb_dofs /
size_symm; ++bb) {
436 t_nf(
L) += t_row_base_fun * t_residual(
L);
440 for (; bb != nb_base_functions; ++bb)
444 ++t_approx_P_adjont_log_du;
◆ alphaU
const double EshelbianPlasticity::HMHHencky::OpSpatialPhysical::alphaU |
|
private |
◆ polyConvex
PetscBool EshelbianPlasticity::HMHHencky::OpSpatialPhysical::polyConvex = PETSC_FALSE |
|
private |
The documentation for this struct was generated from the following file:
- users_modules/eshelbian_plasticit/src/impl/Hencky.cpp
#define CHK_MOAB_THROW(err, msg)
Check error code of MoAB function and throw MoFEM exception.
constexpr static auto size_symm
VectorDouble nF
local right hand side vector
FTensor::Tensor0< FTensor::PackPtr< double *, 1 > > getFTensor0N(const FieldApproximationBase base)
Get base function as Tensor0.
MoFEMErrorCode integrateHencky(EntData &data)
#define CHKERR
Inline error check.
const VectorInt & getIndices() const
Get global indices of dofs on entity.
boost::shared_ptr< DataAtIntegrationPts > dataAtPts
data at integration pts
auto diff_deviator(FTensor::Ddg< double, 3, 3 > &&t_diff_stress)
FTensor::Index< 'i', SPACE_DIM > i
constexpr auto field_name
const double v
phase velocity of light in medium (cm/ns)
MatrixDouble & getN(const FieldApproximationBase base)
get base functions this return matrix (nb. of rows is equal to nb. of Gauss pts, nb....
MoFEMErrorCode integratePolyconvexHencky(EntData &data)
FTensor::Index< 'j', 3 > j
Kronecker Delta class symmetric.
OpAssembleVolume(const std::string &field, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const char type)
FTensor::Index< 'k', 3 > k
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
FTensor::Index< 'l', 3 > l
PetscErrorCode PetscOptionsGetBool(PetscOptions *, const char pre[], const char name[], PetscBool *bval, PetscBool *set)