v0.16.0
Loading...
Searching...
No Matches
Public Member Functions | Private Attributes | List of all members
PlasticOps::OpCalculateConstraintsLhs_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp > Struct Template Reference

#include "tutorials/adv-0_plasticity/src/PlasticOpsGeneric.hpp"

Inheritance diagram for PlasticOps::OpCalculateConstraintsLhs_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >:
[legend]
Collaboration diagram for PlasticOps::OpCalculateConstraintsLhs_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >:
[legend]

Public Member Functions

 OpCalculateConstraintsLhs_dEPImpl (const std::string row_field_name, const std::string col_field_name, boost::shared_ptr< CommonData > common_data_ptr, boost::shared_ptr< MatrixDouble > mat_D_ptr)
 
MoFEMErrorCode iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)
 
 OpCalculateConstraintsLhs_dEPImpl (const std::string row_field_name, const std::string col_field_name, boost::shared_ptr< CommonData > common_data_ptr, boost::shared_ptr< MatrixDouble > mat_D_ptr)
 
MoFEMErrorCode iNtegrate (EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)
 

Private Attributes

boost::shared_ptr< CommonDatacommonDataPtr
 
boost::shared_ptr< MatrixDouble > mDPtr
 

Detailed Description

template<int DIM, typename AssemblyDomainEleOp>
struct PlasticOps::OpCalculateConstraintsLhs_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >

Definition at line 1165 of file PlasticOpsGeneric.hpp.

Constructor & Destructor Documentation

◆ OpCalculateConstraintsLhs_dEPImpl() [1/2]

template<int DIM, typename AssemblyDomainEleOp >
PlasticOps::OpCalculateConstraintsLhs_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::OpCalculateConstraintsLhs_dEPImpl ( const std::string  row_field_name,
const std::string  col_field_name,
boost::shared_ptr< CommonData common_data_ptr,
boost::shared_ptr< MatrixDouble >  mat_D_ptr 
)

Definition at line 1180 of file PlasticOpsGeneric.hpp.

1185 : AssemblyDomainEleOp(row_field_name, col_field_name,
1186 DomainEleOp::OPROWCOL),
1187 commonDataPtr(common_data_ptr), mDPtr(m_D_ptr) {
1188 AssemblyDomainEleOp::sYmm = false;
1189}
FormsIntegrators< DomainEleOp >::Assembly< A >::OpBase AssemblyDomainEleOp

◆ OpCalculateConstraintsLhs_dEPImpl() [2/2]

template<int DIM, typename AssemblyDomainEleOp >
PlasticOps::OpCalculateConstraintsLhs_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::OpCalculateConstraintsLhs_dEPImpl ( const std::string  row_field_name,
const std::string  col_field_name,
boost::shared_ptr< CommonData common_data_ptr,
boost::shared_ptr< MatrixDouble >  mat_D_ptr 
)

Member Function Documentation

◆ iNtegrate() [1/2]

template<int DIM, typename AssemblyDomainEleOp >
MoFEMErrorCode PlasticOps::OpCalculateConstraintsLhs_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::iNtegrate ( EntitiesFieldData::EntData &  row_data,
EntitiesFieldData::EntData &  col_data 
)

Definition at line 1203 of file PlasticOpsGeneric.hpp.

1205 {
1207
1208 FTensor::Index<'i', DIM> i;
1209 FTensor::Index<'j', DIM> j;
1210 constexpr auto size_symm = (DIM * (DIM + 1)) / 2;
1212
1213 const auto nb_integration_pts = AssemblyDomainEleOp::getGaussPts().size2();
1214 const auto nb_row_base_functions = row_data.getN().size2();
1215
1216 auto t_c_dstrain =
1217 getFTensor2SymmetricFromMat<SPACE_DIM>(commonDataPtr->resCdStrain);
1218 auto t_c_dplastic_strain =
1219 getFTensor2SymmetricFromMat<SPACE_DIM>(commonDataPtr->resCdPlasticStrain);
1220
1221 auto next = [&]() {
1222 ++t_c_dstrain;
1223 ++t_c_dplastic_strain;
1224 };
1225
1227
1228 auto t_w = AssemblyDomainEleOp::getFTensor0IntegrationWeight();
1229 auto t_row_base = row_data.getFTensor0N();
1230 for (auto gg = 0; gg != nb_integration_pts; ++gg) {
1231 const double alpha = AssemblyDomainEleOp::getMeasure() * t_w;
1232 ++t_w;
1233
1235 t_res_vec(L) =
1236 t_L(i, j, L) * (t_c_dplastic_strain(i, j) - t_c_dstrain(i, j));
1237 next();
1238
1239 auto t_mat = get_mat_scalar_dtensor_sym(AssemblyDomainEleOp::locMat,
1241 size_t rr = 0;
1242 for (; rr != AssemblyDomainEleOp::nbRows; ++rr) {
1243 const auto row_base = alpha * t_row_base;
1244 auto t_col_base = col_data.getFTensor0N(gg, 0);
1245 for (size_t cc = 0; cc != AssemblyDomainEleOp::nbCols / size_symm; cc++) {
1246 t_mat(L) += (row_base * t_col_base) * t_res_vec(L);
1247 ++t_mat;
1248 ++t_col_base;
1249 }
1250 ++t_row_base;
1251 }
1252 for (; rr != nb_row_base_functions; ++rr)
1253 ++t_row_base;
1254 }
1255
1257}
#define MoFEMFunctionBegin
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
#define MoFEMFunctionReturn(a)
Last executable line of each PETSc function used for error handling. Replaces return()
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'j', 3 > j
auto get_mat_scalar_dtensor_sym(MatrixDouble &mat, FTensor::Number< 2 >)
auto symm_L_tensor(FTensor::Number< DIM >)
constexpr auto size_symm
Definition plastic.cpp:42

◆ iNtegrate() [2/2]

template<int DIM, typename AssemblyDomainEleOp >
MoFEMErrorCode PlasticOps::OpCalculateConstraintsLhs_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::iNtegrate ( EntitiesFieldData::EntData &  row_data,
EntitiesFieldData::EntData &  col_data 
)

Member Data Documentation

◆ commonDataPtr

template<int DIM, typename AssemblyDomainEleOp >
boost::shared_ptr< CommonData > PlasticOps::OpCalculateConstraintsLhs_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::commonDataPtr
private

Definition at line 1175 of file PlasticOpsGeneric.hpp.

◆ mDPtr

template<int DIM, typename AssemblyDomainEleOp >
boost::shared_ptr< MatrixDouble > PlasticOps::OpCalculateConstraintsLhs_dEPImpl< DIM, GAUSS, AssemblyDomainEleOp >::mDPtr
private

Definition at line 1176 of file PlasticOpsGeneric.hpp.


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