v0.15.0
Loading...
Searching...
No Matches
EshelbianPlasticity::OpDispBcImpl< A, GAUSS > Struct Template Reference

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

Inheritance diagram for EshelbianPlasticity::OpDispBcImpl< A, GAUSS >:
[legend]
Collaboration diagram for EshelbianPlasticity::OpDispBcImpl< A, GAUSS >:
[legend]

Public Types

using OP
 

Public Member Functions

 OpDispBcImpl (boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_base_side_data, boost::shared_ptr< BcDispVec > &bc_disp_ptr, std::map< std::string, boost::shared_ptr< ScalingMethod > > smv, boost::shared_ptr< Range > ents_ptr=nullptr)
 

Protected Member Functions

MoFEMErrorCode iNtegrate (EntData &data)
 

Protected Attributes

boost::shared_ptr< BcDispVecbcDispPtr
 
std::map< std::string, boost::shared_ptr< ScalingMethod > > scalingMethodsMap
 

Detailed Description

template<AssemblyType A>
struct EshelbianPlasticity::OpDispBcImpl< A, GAUSS >

Definition at line 420 of file EshelbianPlasticity.hpp.

Member Typedef Documentation

◆ OP

template<AssemblyType A>
using EshelbianPlasticity::OpDispBcImpl< A, GAUSS >::OP
Initial value:
typename FormsIntegrators<FaceUserDataOperator>::Assembly<

Definition at line 423 of file EshelbianPlasticity.hpp.

Constructor & Destructor Documentation

◆ OpDispBcImpl()

template<AssemblyType A>
EshelbianPlasticity::OpDispBcImpl< A, GAUSS >::OpDispBcImpl ( boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_base_side_data,
boost::shared_ptr< BcDispVec > & bc_disp_ptr,
std::map< std::string, boost::shared_ptr< ScalingMethod > > smv,
boost::shared_ptr< Range > ents_ptr = nullptr )
inline

Definition at line 426 of file EshelbianPlasticity.hpp.

Member Function Documentation

◆ iNtegrate()

template<AssemblyType A>
MoFEMErrorCode OpDispBcImpl< A, GAUSS >::iNtegrate ( EntData & data)
protected

Definition at line 1423 of file EshelbianOperators.cpp.

1423 {
1425 // get entity of face
1426 EntityHandle fe_ent = OP::getFEEntityHandle();
1427 // iterate over all boundary data
1428 for (auto &bc : (*bcDispPtr)) {
1429 // check if finite element entity is part of boundary condition
1430 if (bc.faces.find(fe_ent) != bc.faces.end()) {
1431 int nb_dofs = data.getIndices().size();
1432
1433 int nb_integration_pts = OP::getGaussPts().size2();
1434 auto t_normal = OP::getFTensor1NormalsAtGaussPts();
1435 auto t_w = OP::getFTensor0IntegrationWeight();
1436 int nb_base_functions = data.getN().size2() / 3;
1437 auto t_row_base_fun = data.getFTensor1N<3>();
1438
1441
1442 double scale = 1;
1443 if (scalingMethodsMap.find(bc.blockName) != scalingMethodsMap.end()) {
1445 scale *= scalingMethodsMap.at(bc.blockName)
1446 ->getScale(EshelbianCore::dynamicTime);
1447 } else {
1448 scale *= scalingMethodsMap.at(bc.blockName)
1449 ->getScale(OP::getFEMethod()->ts_t);
1450 }
1451 } else {
1452 MOFEM_LOG("SELF", Sev::warning)
1453 << "No scaling method found for " << bc.blockName;
1454 }
1455
1456 // get bc data
1457 FTensor::Tensor1<double, 3> t_bc_disp(bc.vals[0], bc.vals[1], bc.vals[2]);
1458 t_bc_disp(i) *= scale;
1459
1460 for (int gg = 0; gg != nb_integration_pts; ++gg) {
1461 auto t_nf = getFTensor1FromPtr<3>(&*OP::locF.begin());
1462 int bb = 0;
1463 for (; bb != nb_dofs / SPACE_DIM; ++bb) {
1464 t_nf(i) +=
1465 t_w * (t_row_base_fun(j) * t_normal(j)) * t_bc_disp(i) * 0.5;
1466 ++t_nf;
1467 ++t_row_base_fun;
1468 }
1469 for (; bb != nb_base_functions; ++bb)
1470 ++t_row_base_fun;
1471
1472 ++t_w;
1473 ++t_normal;
1474 }
1475 }
1476 }
1478}
#define FTENSOR_INDEX(DIM, I)
constexpr int SPACE_DIM
#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()
#define MOFEM_LOG(channel, severity)
Log.
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'j', 3 > j
FTensor::Tensor1< FTensor::PackPtr< double *, S >, DIM > getFTensor1FromPtr(double *ptr)
Make Tensor1 from pointer.
double scale
Definition plastic.cpp:123
std::map< std::string, boost::shared_ptr< ScalingMethod > > scalingMethodsMap
MatrixDouble & getN(const FieldApproximationBase base)
get base functions this return matrix (nb. of rows is equal to nb. of Gauss pts, nb....
FTensor::Tensor1< FTensor::PackPtr< double *, Tensor_Dim >, Tensor_Dim > getFTensor1N(FieldApproximationBase base)
Get base functions for Hdiv/Hcurl spaces.
const VectorInt & getIndices() const
Get global indices of dofs on entity.

Member Data Documentation

◆ bcDispPtr

template<AssemblyType A>
boost::shared_ptr<BcDispVec> EshelbianPlasticity::OpDispBcImpl< A, GAUSS >::bcDispPtr
protected

Definition at line 436 of file EshelbianPlasticity.hpp.

◆ scalingMethodsMap

template<AssemblyType A>
std::map<std::string, boost::shared_ptr<ScalingMethod> > EshelbianPlasticity::OpDispBcImpl< A, GAUSS >::scalingMethodsMap
protected

Definition at line 437 of file EshelbianPlasticity.hpp.


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