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

Inheritance diagram for EshelbianPlasticity::OpSpatialEquilibrium:

Collaboration diagram for EshelbianPlasticity::OpSpatialEquilibrium:

Public Member Functions
	OpSpatialEquilibrium (const std::string &field_name, boost::shared_ptr< DataAtIntegrationPts > &data_ptr, const double alpha, const double rho)

MoFEMErrorCode	integrate (EntData &data)

Public Member Functions inherited from EshelbianPlasticity::OpAssembleVolume
	OpAssembleVolume (const std::string &field, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const char type)

	OpAssembleVolume (const std::string &row_field, const std::string &col_field, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const char type, const bool assemble_symmetry)

Public Member Functions inherited from EshelbianPlasticity::OpAssembleBasic< VolUserDataOperator >
	OpAssembleBasic (const std::string &field_name, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const char type)

	OpAssembleBasic (const std::string &row_field, const std::string &col_field, boost::shared_ptr< DataAtIntegrationPts > data_ptr, const char type, const bool assemble_symmetry)

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)

Public Attributes
const double	alphaW

const double	alphaRho

Public Attributes inherited from EshelbianPlasticity::OpAssembleBasic< VolUserDataOperator >
const bool	assembleSymmetry

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

VectorDouble	nF
	local right hand side vector More...

MatrixDouble	K
	local tangent matrix More...

MatrixDouble	transposeK

Detailed Description

Definition at line 768 of file EshelbianPlasticity.hpp.

Constructor & Destructor Documentation

◆ OpSpatialEquilibrium()

EshelbianPlasticity::OpSpatialEquilibrium::OpSpatialEquilibrium	(	const std::string &	field_name,
		boost::shared_ptr< DataAtIntegrationPts > &	data_ptr,
		const double	alpha,
		const double	rho
	)

Definition at line 771 of file EshelbianPlasticity.hpp.

774 : OpAssembleVolume(field_name, data_ptr, OPROW), alphaW(alpha),

775 alphaRho(rho) {}

Member Function Documentation

◆ integrate()

MoFEMErrorCode EshelbianPlasticity::OpSpatialEquilibrium::integrate ( EntData & data )

virtual

Reimplemented from EshelbianPlasticity::OpAssembleBasic< VolUserDataOperator >.

Definition at line 252 of file EshelbianOperators.cpp.

                                                             {
   MoFEMFunctionBegin;
   int nb_dofs = data.getIndices().size();
   int nb_integration_pts = data.getN().size1();
   auto v = getVolume();
   auto t_w = getFTensor0IntegrationWeight();
   auto t_div_P = getFTensor1FromMat<3>(dataAtPts->divPAtPts);
   auto t_s_dot_w = getFTensor1FromMat<3>(dataAtPts->wDotAtPts);
   if (dataAtPts->wDotDotAtPts.size1() == 0 &&
       dataAtPts->wDotDotAtPts.size2() != nb_integration_pts) {
     dataAtPts->wDotDotAtPts.resize(3, nb_integration_pts);
     dataAtPts->wDotDotAtPts.clear();
   }
   auto t_s_dot_dot_w = getFTensor1FromMat<3>(dataAtPts->wDotDotAtPts);
  
   int nb_base_functions = data.getN().size2();
   auto t_row_base_fun = data.getFTensor0N();
  
   FTensor::Index<'i', 3> i;
   auto get_ftensor1 = [](auto &v) {
     return FTensor::Tensor1<FTensor::PackPtr<double *, 3>, 3>(&v[0], &v[1],
                                                               &v[2]);
   };
  
   for (int gg = 0; gg != nb_integration_pts; ++gg) {
     double a = v * t_w;
     auto t_nf = get_ftensor1(nF);
     int bb = 0;
     for (; bb != nb_dofs / 3; ++bb) {
       t_nf(i) += a * t_row_base_fun * t_div_P(i);
       t_nf(i) -= a * t_row_base_fun * alphaW * t_s_dot_w(i);
       t_nf(i) -= a * t_row_base_fun * alphaRho * t_s_dot_dot_w(i);
       ++t_nf;
       ++t_row_base_fun;
     }
     for (; bb != nb_base_functions; ++bb)
       ++t_row_base_fun;
     ++t_w;
     ++t_div_P;
     ++t_s_dot_w;
     ++t_s_dot_dot_w;
   }
  
   MoFEMFunctionReturn(0);
 }

Member Data Documentation

◆ alphaRho

const double EshelbianPlasticity::OpSpatialEquilibrium::alphaRho

Definition at line 770 of file EshelbianPlasticity.hpp.

◆ alphaW

const double EshelbianPlasticity::OpSpatialEquilibrium::alphaW

Definition at line 769 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

Public Member Functions

Public Attributes