EshelbianPlasticity::OpCalculateReactionForces Struct Reference

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

Inheritance diagram for EshelbianPlasticity::OpCalculateReactionForces:

Collaboration diagram for EshelbianPlasticity::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 324 of file EshelbianPlasticity.hpp.

Constructor & Destructor Documentation

OpCalculateReactionForces()

EshelbianPlasticity::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 325 of file EshelbianPlasticity.hpp.

Member Function Documentation

doWork()

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

Examples

EshelbianOperators.cpp.

Definition at line 913 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 = getFTensor1FromMat<SPACE_DIM>(dataAtPts->tractionAtPts);
  auto t_coords = getFTensor1CoordsAtGaussPts();
  auto t_spatial_disp = getFTensor1FromMat<SPACE_DIM>(dataAtPts->wL2AtPts);
 
  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) {
    loc_reaction_forces(i) += (t_traction(i)) * t_w * getMeasure();
    t_coords_spatial(i) = t_coords(i) + t_spatial_disp(i);
    // t_coords_spatial(i) -= t_off(i);
    loc_moment_forces(i) +=
        (FTensor::levi_civita<double>(i, j, k) * t_coords_spatial(j)) *
        t_traction(k) * t_w * getMeasure();
    ++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 EshelbianPlasticity::OpCalculateReactionForces::blockEntities

private

Examples

EshelbianOperators.cpp.

Definition at line 336 of file EshelbianPlasticity.hpp.

blockName

std::string EshelbianPlasticity::OpCalculateReactionForces::blockName

private

Definition at line 337 of file EshelbianPlasticity.hpp.

dataAtPts

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

private

data at integration pts

Examples

EshelbianOperators.cpp.

Definition at line 335 of file EshelbianPlasticity.hpp.

reactionVec

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

private

Examples

EshelbianOperators.cpp.

Definition at line 338 of file EshelbianPlasticity.hpp.

---

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

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