v0.16.3
Loading...
Searching...
No Matches
Public Types | Public Member Functions | Protected Attributes | List of all members
OpSpringLhsBc_dP Struct Reference

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

Inheritance diagram for OpSpringLhsBc_dP:
[legend]
Collaboration diagram for OpSpringLhsBc_dP:
[legend]

Public Types

using OP = typename FormsIntegrators< FaceUserDataOperator >::Assembly< A >::OpBrokenBase
 

Public Member Functions

 OpSpringLhsBc_dP (std::string row_field, boost::shared_ptr< std::vector< BrokenBaseSideData > > broken_base_side_data, boost::shared_ptr< SpringBcVec > bc_spring_ptr)
 
MoFEMErrorCode iNtegrate (EntData &row_data, EntData &col_data)
 

Protected Attributes

boost::shared_ptr< SpringBcVec > bcSpringPtr
 

Detailed Description

Examples
/home/lk58p/mofem_install/vanilla_dev_release/mofem-cephas/mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianPlasticity.cpp.

Definition at line 605 of file EshelbianOperators.hpp.

Member Typedef Documentation

◆ OP

using OpSpringLhsBc_dP::OP = typename FormsIntegrators<FaceUserDataOperator>::Assembly<A>::OpBrokenBase

Definition at line 608 of file EshelbianOperators.hpp.

Constructor & Destructor Documentation

◆ OpSpringLhsBc_dP()

OpSpringLhsBc_dP::OpSpringLhsBc_dP ( std::string  row_field,
boost::shared_ptr< std::vector< BrokenBaseSideData > >  broken_base_side_data,
boost::shared_ptr< SpringBcVec >  bc_spring_ptr 
)
inline

Definition at line 611 of file EshelbianOperators.hpp.

615 : OP(row_field, broken_base_side_data, false, false, nullptr),
616 bcSpringPtr(bc_spring_ptr) {}
boost::shared_ptr< SpringBcVec > bcSpringPtr
typename FormsIntegrators< FaceUserDataOperator >::Assembly< A >::OpBrokenBase OP

Member Function Documentation

◆ iNtegrate()

MoFEMErrorCode OpSpringLhsBc_dP::iNtegrate ( EntData &  row_data,
EntData &  col_data 
)
Examples
/home/lk58p/mofem_install/vanilla_dev_release/mofem-cephas/mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp.

Definition at line 2540 of file EshelbianOperators.cpp.

2541 {
2543
2544 int row_nb_dofs = row_data.getIndices().size();
2545 int col_nb_dofs = col_data.getIndices().size();
2546 auto &locMat = OP::locMat;
2547 locMat.resize(row_nb_dofs, col_nb_dofs, false);
2548 locMat.clear();
2549
2550 // get entity of face
2551 EntityHandle fe_ent = OP::getFEEntityHandle();
2552 // iterate over all boundary data
2553 for (auto &bc : (*bcSpringPtr)) {
2554 // check if finite element entity is part of boundary condition
2555 if (bc.faces.find(fe_ent) != bc.faces.end()) {
2556
2557 auto t_normal = OP::getFTensor1NormalsAtGaussPts();
2558 auto t_w = OP::getFTensor0IntegrationWeight();
2559
2563
2564 int nb_integration_pts = OP::getGaussPts().size2();
2565 int nb_base_functions = row_data.getN().size2();
2566 auto t_row_base = row_data.getFTensor0N();
2567
2569
2570 for (int gg = 0; gg != nb_integration_pts; ++gg) {
2571
2573 t_N(i) = t_normal(i);
2574 t_N.normalize();
2575
2576 int rr = 0;
2577 for (; rr != row_nb_dofs / SPACE_DIM; ++rr) {
2578 auto t_mat = getFTensor2FromArray<SPACE_DIM, SPACE_DIM, SPACE_DIM>(
2579 OP::locMat, SPACE_DIM * rr);
2580 auto t_col_base = col_data.getFTensor1N<3>(gg, 0);
2581 for (auto cc = 0; cc != col_nb_dofs / SPACE_DIM; ++cc) {
2582 t_mat(i, j) +=
2583 ((t_w * t_row_base) * (t_N(k) * t_col_base(k))) * 0.5 * t_kd(i, j);
2584 ++t_mat;
2585 ++t_col_base;
2586 }
2587 ++t_row_base;
2588 }
2589
2590 for (; rr != nb_base_functions; ++rr)
2591 ++t_row_base;
2592
2593 ++t_w;
2594 ++t_normal;
2595 }
2596
2597 locMat *= OP::getMeasure();
2598 }
2599 }
2601}
#define FTENSOR_INDEX(DIM, I)
constexpr int SPACE_DIM
Kronecker Delta class.
Tensor1< T, Tensor_Dim > normalize()
#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()
constexpr auto t_kd
FTensor::Index< 'i', SPACE_DIM > i
FTensor::Index< 'j', 3 > j
FTensor::Index< 'k', 3 > k
FTensor::Tensor0< FTensor::PackPtr< double *, 1 > > getFTensor0N(const FieldApproximationBase base)
Get base function as Tensor0.
MatrixDouble & getN(const FieldApproximationBase base)
get base functions this return matrix (nb. of rows is equal to nb. of Gauss pts, nb....
auto getFTensor1N(FieldApproximationBase base)
Get base functions for Hdiv/Hcurl spaces.
const VectorInt & getIndices() const
Get global indices of degrees of freedom on entity.

Member Data Documentation

◆ bcSpringPtr

boost::shared_ptr<SpringBcVec> OpSpringLhsBc_dP::bcSpringPtr
protected

Definition at line 621 of file EshelbianOperators.hpp.


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