[Auxiliary functions functions More...

#include <tutorials/adv-0/src/PlasticOpsGeneric.hpp>

Inheritance diagram for PlasticOps::OpCalculatePlasticSurfaceImpl< DIM, GAUSS, DomainEleOp >:

Collaboration diagram for PlasticOps::OpCalculatePlasticSurfaceImpl< DIM, GAUSS, DomainEleOp >:

Public Member Functions
	OpCalculatePlasticSurfaceImpl (const std::string field_name, boost::shared_ptr< CommonData > common_data_ptr)

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

Protected Attributes
boost::shared_ptr< CommonData >	commonDataPtr

Detailed Description

template<int DIM, typename DomainEleOp>
struct PlasticOps::OpCalculatePlasticSurfaceImpl< DIM, GAUSS, DomainEleOp >

[Auxiliary functions functions

Definition at line 387 of file PlasticOpsGeneric.hpp.

Constructor & Destructor Documentation

◆ OpCalculatePlasticSurfaceImpl()

template<int DIM, typename DomainEleOp >

PlasticOps::OpCalculatePlasticSurfaceImpl< DIM, GAUSS, DomainEleOp >::OpCalculatePlasticSurfaceImpl	(	const std::string	field_name,
		boost::shared_ptr< CommonData >	common_data_ptr
	)

Definition at line 399 of file PlasticOpsGeneric.hpp.

     : DomainEleOp(field_name, DomainEleOp::OPROW),
       commonDataPtr(common_data_ptr) {
   // Operator is only executed for vertices
   std::fill(&DomainEleOp::doEntities[MBEDGE],
             &DomainEleOp::doEntities[MBMAXTYPE], false);
 }

Member Function Documentation

◆ doWork()

template<int DIM, typename DomainEleOp >

MoFEMErrorCode PlasticOps::OpCalculatePlasticSurfaceImpl< DIM, GAUSS, DomainEleOp >::doWork	(	int	side,
		EntityType	type,
		EntData &	data
	)

Examples: PlasticOpsGeneric.hpp.

Definition at line 409 of file PlasticOpsGeneric.hpp.

                                               {
   MoFEMFunctionBegin;
  
   FTensor::Index<'i', DIM> i;
   FTensor::Index<'j', DIM> j;
  
   const size_t nb_gauss_pts = commonDataPtr->mStressPtr->size2();
   auto t_stress =
       getFTensor2SymmetricFromMat<DIM>(*(commonDataPtr->mStressPtr));
   auto t_plastic_strain =
       getFTensor2SymmetricFromMat<DIM>(commonDataPtr->plasticStrain);
  
   commonDataPtr->plasticSurface.resize(nb_gauss_pts, false);
   commonDataPtr->plasticFlow.resize(size_symm, nb_gauss_pts, false);
   auto t_flow = getFTensor2SymmetricFromMat<DIM>(commonDataPtr->plasticFlow);
  
   auto &params = commonDataPtr->blockParams;
  
   for (auto &f : commonDataPtr->plasticSurface) {
  
     f = platsic_surface(
  
         deviator(
             t_stress, trace(t_stress),
             kinematic_hardening(t_plastic_strain, params[CommonData::C1_k]))
  
     );
  
     auto t_flow_tmp =
         plastic_flow(f,
  
                      deviator(t_stress, trace(t_stress),
                               kinematic_hardening(t_plastic_strain,
                                                   params[CommonData::C1_k])),
  
                      diff_deviator(diff_tensor(FTensor::Number<DIM>())));
     t_flow(i, j) = t_flow_tmp(i, j);
  
     ++t_plastic_strain;
     ++t_flow;
     ++t_stress;
   }
  
   MoFEMFunctionReturn(0);
 }