OpCalculateReactionForces Struct Reference

#include "users_modules/eshelbian_plasticity/src/EshelbianOperators.hpp"

Inheritance diagram for OpCalculateReactionForces:

Collaboration diagram for OpCalculateReactionForces:

Public Member Functions
	OpCalculateReactionForces (boost::shared_ptr< DataAtIntegrationPts > data_ptr, std::string block_name, Range block_entities, std::array< double, 6 > &reaction_vec)
MoFEMErrorCode	doWork (int side, EntityType type, EntData &data)

Private Attributes
boost::shared_ptr< DataAtIntegrationPts >	dataAtPts
	data at integration pts
Range	blockEntities
std::string	blockName
std::array< double, 6 > &	reactionVec

Detailed Description

Definition at line 316 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

OpCalculateReactionForces()

OpCalculateReactionForces::OpCalculateReactionForces ( boost::shared_ptr< DataAtIntegrationPts >  data_ptr, std::string  block_name, Range  block_entities, std::array< double, 6 > &  reaction_vec  )

inline

Definition at line 317 of file EshelbianOperators.hpp.

      : FaceUserDataOperator(NOSPACE, OPSPACE), dataAtPts(data_ptr),
        blockEntities(block_entities), blockName(block_name),
        reactionVec(reaction_vec) {}

Member Function Documentation

doWork()

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

Examples

/home/lk58p/mofem_install/vanilla_dev_release/mofem-cephas/mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 975 of file EshelbianOperators.cpp.

                                                                {
  MoFEMFunctionBegin;
  if (blockEntities.find(getFEEntityHandle()) == blockEntities.end()) {
    MoFEMFunctionReturnHot(0);
  };
  FTENSOR_INDEX(SPACE_DIM, i);
  FTENSOR_INDEX(SPACE_DIM, j);
  FTENSOR_INDEX(SPACE_DIM, k);
  int nb_integration_pts = getGaussPts().size2();
  auto t_w = getFTensor0IntegrationWeight();
  auto t_traction = dataAtPts->getFTensorTraction(nb_integration_pts);
  auto t_coords = getFTensor1CoordsAtGaussPts();
  auto t_spatial_disp = dataAtPts->getFTensorSmallWL2(nb_integration_pts);
 
  FTensor::Tensor1<double, 3> t_coords_spatial{0., 0., 0.};
  // Offset for center of mass. Can be added in the future.
  FTensor::Tensor1<double, 3> t_off{0.0, 0.0, 0.0};
  FTensor::Tensor1<double, 3> loc_reaction_forces{0., 0., 0.};
  FTensor::Tensor1<double, 3> loc_moment_forces{0., 0., 0.};
 
  for (auto gg = 0; gg != nb_integration_pts; ++gg) {
    double a = t_w * getMeasure();
    loc_reaction_forces(i) += a*t_traction(i);
    t_coords_spatial(i) = t_coords(i) + t_spatial_disp(i);
    // t_coords_spatial(i) -= t_off(i);
    loc_moment_forces(i) +=
        (a * (FTensor::levi_civita<double>(i, j, k) * t_coords_spatial(j))) *
        t_traction(k);
    ++t_coords;
    ++t_spatial_disp;
    ++t_w;
    ++t_traction;
  }
 
  reactionVec[0] += loc_reaction_forces(0);
  reactionVec[1] += loc_reaction_forces(1);
  reactionVec[2] += loc_reaction_forces(2);
  reactionVec[3] += loc_moment_forces(0);
  reactionVec[4] += loc_moment_forces(1);
  reactionVec[5] += loc_moment_forces(2);
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

blockEntities

Range OpCalculateReactionForces::blockEntities

private

Definition at line 328 of file EshelbianOperators.hpp.

blockName

std::string OpCalculateReactionForces::blockName

private

Definition at line 329 of file EshelbianOperators.hpp.

dataAtPts

boost::shared_ptr<DataAtIntegrationPts> OpCalculateReactionForces::dataAtPts

private

data at integration pts

Definition at line 327 of file EshelbianOperators.hpp.

reactionVec

std::array<double, 6>& OpCalculateReactionForces::reactionVec

private

Definition at line 330 of file EshelbianOperators.hpp.

---

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

• users_modules/eshelbian_plasticity/src/EshelbianOperators.hpp
• users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp