EshelbianPlasticity::OpCalculateRotationAndSpatialGradient Struct Reference

#include <users_modules/eshelbian_plasticty/src/EshelbianPlasticity.hpp>

Inheritance diagram for EshelbianPlasticity::OpCalculateRotationAndSpatialGradient:

Collaboration diagram for EshelbianPlasticity::OpCalculateRotationAndSpatialGradient:

Public Member Functions
	OpCalculateRotationAndSpatialGradient (const std::string &field_name, boost::shared_ptr< DataAtIntegrationPts > &data_ptr)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Public Attributes
boost::shared_ptr< DataAtIntegrationPts >	dataAtPts
	data at integration pts More...

Detailed Description

Definition at line 513 of file EshelbianPlasticity.hpp.

Constructor & Destructor Documentation

OpCalculateRotationAndSpatialGradient()

EshelbianPlasticity::OpCalculateRotationAndSpatialGradient::OpCalculateRotationAndSpatialGradient ( const std::string &  field_name, boost::shared_ptr< DataAtIntegrationPts > &  data_ptr  )

inline

Definition at line 516 of file EshelbianPlasticity.hpp.

      : VolUserDataOperator(field_name, OPROW), dataAtPts(data_ptr) {}

Member Function Documentation

doWork()

MoFEMErrorCode EshelbianPlasticity::OpCalculateRotationAndSpatialGradient::doWork	(	int	side,
		EntityType	type,
		EntData &	data
	)

Examples

EshelbianOperators.cpp.

Definition at line 166 of file EshelbianOperators.cpp.

                                                                            {
  MoFEMFunctionBegin;
  if (type != MBTET)
    MoFEMFunctionReturnHot(0);
  int nb_integration_pts = data.getN().size1();
  FTensor::Index<'i', 3> i;
  FTensor::Index<'j', 3> j;
  FTensor::Index<'k', 3> k;
  FTensor::Index<'l', 3> l;
 
  dataAtPts->hAtPts.resize(9, nb_integration_pts, false);
  dataAtPts->rotMatAtPts.resize(9, nb_integration_pts, false);
  dataAtPts->diffRotMatAtPts.resize(27, nb_integration_pts, false);
  dataAtPts->streachTensorAtPts.resize(6, nb_integration_pts, false);
  dataAtPts->diffStreachTensorAtPts.resize(36, nb_integration_pts, false);
  dataAtPts->eigenVals.resize(3, nb_integration_pts, false);
  dataAtPts->eigenVecs.resize(9, nb_integration_pts, false);
  dataAtPts->nbUniq.resize(nb_integration_pts, false);
 
  auto t_h = getFTensor2FromMat<3, 3>(dataAtPts->hAtPts);
  auto t_omega = getFTensor1FromMat<3>(dataAtPts->rotAxisAtPts);
  auto t_R = getFTensor2FromMat<3, 3>(dataAtPts->rotMatAtPts);
  auto t_diff_R = getFTensor3FromMat(dataAtPts->diffRotMatAtPts);
  auto t_log_u =
      getFTensor2SymmetricFromMat<3>(dataAtPts->logStreachTensorAtPts);
  auto t_u = getFTensor2SymmetricFromMat<3>(dataAtPts->streachTensorAtPts);
 
  auto t_diff_u =
      getFTensor4DdgFromMat<3, 3, 1>(dataAtPts->diffStreachTensorAtPts);
  auto t_eigen_vals = getFTensor1FromMat<3>(dataAtPts->eigenVals);
  auto t_eigen_vecs = getFTensor2FromMat<3, 3>(dataAtPts->eigenVecs);
  auto &nbUniq = dataAtPts->nbUniq;
 
  for (int gg = 0; gg != nb_integration_pts; ++gg) {
 
    auto t0_diff = get_diff_rotation_form_vector(t_omega);
    auto t0 = get_rotation_form_vector(t_omega);
 
    t_diff_R(i, j, k) = t0_diff(i, j, k);
    t_R(i, j) = t0(i, j);
 
    FTensor::Tensor1<double, 3> eig;
    FTensor::Tensor2<double, 3, 3> eigen_vec;
    CHKERR get_eigen_val_and_proj_lapack(t_log_u, eig, eigen_vec);
 
    t_eigen_vals(i) = eig(i);
    t_eigen_vecs(i, j) = eigen_vec(i, j);
 
    // rare case when two eigen values are equal
    nbUniq[gg] = get_uniq_nb(eig);
    if (nbUniq[gg] == 2)
      sort_eigen_vals(eig, eigen_vec);
 
    auto t_u_tmp = EigenMatrix::getMat(t_eigen_vals, t_eigen_vecs, f);
    auto t_diff_u_tmp =
        EigenMatrix::getDiffMat(t_eigen_vals, t_eigen_vecs, f, d_f, nbUniq[gg]);
    for (int ii = 0; ii != 3; ++ii) {
      for (int jj = ii; jj != 3; ++jj) {
        for (int kk = 0; kk != 3; ++kk) {
          for (int ll = 0; ll < kk; ++ll) {
            t_diff_u_tmp(ii, jj, kk, ll) *= 2;
          }
        }
      }
    }
 
    t_u(i, j) = t_u_tmp(i, j);
    t_diff_u(i, j, k, l) = t_diff_u_tmp(i, j, k, l);
    t_h(i, j) = t_R(i, k) * t_u(k, j);
 
    ++t_h;
    ++t_R;
    ++t_diff_R;
    ++t_log_u;
    ++t_u;
    ++t_diff_u;
    ++t_eigen_vals;
    ++t_eigen_vecs;
    ++t_omega;
  }
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

dataAtPts

boost::shared_ptr<DataAtIntegrationPts> EshelbianPlasticity::OpCalculateRotationAndSpatialGradient::dataAtPts

data at integration pts

Examples

EshelbianOperators.cpp.

Definition at line 515 of file EshelbianPlasticity.hpp.

---

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

• users_modules/eshelbian_plasticty/src/EshelbianPlasticity.hpp
• users_modules/eshelbian_plasticty/src/impl/EshelbianOperators.cpp