v0.13.0
Public Member Functions | Public Attributes | List of all members
NitscheMethod::OpLhsNormal Struct Reference

Calculate Nitsche method terms on left hand side. More...

#include <users_modules/homogenisation/src/NitscheMethod.hpp>

Inheritance diagram for NitscheMethod::OpLhsNormal:
[legend]
Collaboration diagram for NitscheMethod::OpLhsNormal:
[legend]

Public Member Functions

 OpLhsNormal (const std::string field_name, BlockData &nitsche_block_data, CommonData &nitsche_common_data, NonlinearElasticElement::BlockData &data, NonlinearElasticElement::CommonData &common_data, bool field_disp)
 
MoFEMErrorCode doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, DataForcesAndSourcesCore::EntData &row_data, DataForcesAndSourcesCore::EntData &col_data)
 
- Public Member Functions inherited from NitscheMethod::OpCommon
 OpCommon (const std::string field_name, BlockData &nitsche_block_data, CommonData &nitsche_common_data, NonlinearElasticElement::BlockData &data, NonlinearElasticElement::CommonData &common_data, bool field_disp, const char type)
 
virtual MoFEMErrorCode calculateP (int gg, int fgg, int ff)
 
MoFEMErrorCode getJac (DataForcesAndSourcesCore::EntData &data, int gg, MatrixDouble &jac)
 
MoFEMErrorCode getTractionVariance (int gg, int fgg, int ff, MatrixDouble &jac, MatrixDouble &trac)
 
- Public Member Functions inherited from NitscheMethod::OpBasicCommon
 OpBasicCommon (const std::string field_name, BlockData &nitsche_block_data, CommonData &nitsche_common_data, bool field_disp, const char type)
 
virtual MoFEMErrorCode getFaceRadius (int ff)
 
virtual MoFEMErrorCode getGammaH (double gamma, int gg)
 
- Public Member Functions inherited from MoFEM::VolumeElementForcesAndSourcesCoreBase::UserDataOperator
int getNumNodes ()
 get element number of nodes More...
 
const EntityHandlegetConn ()
 get element connectivity More...
 
double getVolume () const
 element volume (linear geometry) More...
 
double & getVolume ()
 element volume (linear geometry) More...
 
FTensor::Tensor2< double *, 3, 3 > & getJac ()
 get element Jacobian More...
 
FTensor::Tensor2< double *, 3, 3 > & getInvJac ()
 get element inverse Jacobian More...
 
VectorDoublegetCoords ()
 nodal coordinates More...
 
VolumeElementForcesAndSourcesCoreBasegetVolumeFE () const
 return pointer to Generic Volume Finite Element object More...
 
- Public Member Functions inherited from MoFEM::ForcesAndSourcesCore::UserDataOperator
 UserDataOperator (const FieldSpace space, const char type=OPSPACE, const bool symm=true)
 
 UserDataOperator (const std::string &field_name, const char type, const bool symm=true)
 
 UserDataOperator (const std::string &row_field_name, const std::string &col_field_name, const char type, const bool symm=true)
 
boost::shared_ptr< const NumeredEntFiniteElementgetNumeredEntFiniteElementPtr () const
 Return raw pointer to NumeredEntFiniteElement. More...
 
EntityHandle getFEEntityHandle () const
 Return finite element entity handle. More...
 
int getFEDim () const
 Get dimension of finite element. More...
 
EntityType getFEType () const
 Get dimension of finite element. More...
 
boost::weak_ptr< SideNumbergetSideNumberPtr (const int side_number, const EntityType type)
 Get the side number pointer. More...
 
EntityHandle getSideEntity (const int side_number, const EntityType type)
 Get the side entity. More...
 
int getNumberOfNodesOnElement () const
 Get the number of nodes on finite element. More...
 
MoFEMErrorCode getProblemRowIndices (const std::string filed_name, const EntityType type, const int side, VectorInt &indices) const
 Get row indices. More...
 
MoFEMErrorCode getProblemColIndices (const std::string filed_name, const EntityType type, const int side, VectorInt &indices) const
 Get col indices. More...
 
const FEMethodgetFEMethod () const
 Return raw pointer to Finite Element Method object. More...
 
int getOpType () const
 Get operator types. More...
 
void setOpType (const OpType type)
 Set operator type. More...
 
void addOpType (const OpType type)
 Add operator type. More...
 
int getNinTheLoop () const
 get number of finite element in the loop More...
 
int getLoopSize () const
 get size of elements in the loop More...
 
const std::string & getFEName () const
 Get name of the element. More...
 
const PetscData::SwitchesgetDataCtx () const
 
const KspMethod::KSPContext getKSPCtx () const
 
const SnesMethod::SNESContext getSNESCtx () const
 
const TSMethod::TSContext getTSCtx () const
 
Vec getKSPf () const
 
Mat getKSPA () const
 
Mat getKSPB () const
 
Vec getSNESf () const
 
Vec getSNESx () const
 
Mat getSNESA () const
 
Mat getSNESB () const
 
Vec getTSu () const
 
Vec getTSu_t () const
 
Vec getTSu_tt () const
 
Vec getTSf () const
 
Mat getTSA () const
 
Mat getTSB () const
 
int getTSstep () const
 
double getTStime () const
 
double getTSa () const
 
double getTSaa () const
 
MatrixDoublegetGaussPts ()
 matrix of integration (Gauss) points for Volume Element More...
 
auto getFTensor0IntegrationWeight ()
 Get integration weights. More...
 
MatrixDoublegetCoordsAtGaussPts ()
 Gauss points and weight, matrix (nb. of points x 3) More...
 
auto getFTensor1CoordsAtGaussPts ()
 Get coordinates at integration points assuming linear geometry. More...
 
double getMeasure () const
 get measure of element More...
 
double & getMeasure ()
 get measure of element More...
 
MoFEMErrorCode loopSide (const string &fe_name, ForcesAndSourcesCore *side_fe, const size_t dim, const EntityHandle ent_for_side=0, const int verb=QUIET, const LogManager::SeverityLevel sev=Sev::noisy)
 User call this function to loop over elements on the side of face. This function calls finite element with is operator to do calculations. More...
 
MoFEMErrorCode loopThis (const string &fe_name, ForcesAndSourcesCore *parent_fe, const int verb=QUIET, const LogManager::SeverityLevel sev=Sev::noisy)
 User call this function to loop over parent elements. This function calls finite element with is operator to do calculations. More...
 
MoFEMErrorCode loopParent (const string &fe_name, ForcesAndSourcesCore *parent_fe, const int verb=QUIET, const LogManager::SeverityLevel sev=Sev::noisy)
 User call this function to loop over parent elements. This function calls finite element with is operator to do calculations. More...
 
MoFEMErrorCode loopChildren (const string &fe_name, ForcesAndSourcesCore *child_fe, const int verb=QUIET, const LogManager::SeverityLevel sev=Sev::noisy)
 User call this function to loop over parent elements. This function calls finite element with is operator to do calculations. More...
 
- Public Member Functions inherited from MoFEM::DataOperator
 DataOperator (const bool symm=true)
 
virtual ~DataOperator ()=default
 
virtual MoFEMErrorCode doWork (int row_side, int col_side, EntityType row_type, EntityType col_type, EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)
 Operator for bi-linear form, usually to calculate values on left hand side. More...
 
virtual MoFEMErrorCode opLhs (EntitiesFieldData &row_data, EntitiesFieldData &col_data)
 
virtual MoFEMErrorCode doWork (int side, EntityType type, EntitiesFieldData::EntData &data)
 Operator for linear form, usually to calculate values on right hand side. More...
 
virtual MoFEMErrorCode opRhs (EntitiesFieldData &data, const bool error_if_no_base=false)
 
bool getSymm () const
 Get if operator uses symmetry of DOFs or not. More...
 
void setSymm ()
 set if operator is executed taking in account symmetry More...
 
void unSetSymm ()
 unset if operator is executed for non symmetric problem More...
 

Public Attributes

MatrixDouble kMatrix
 
MatrixDouble kMatrix0
 
MatrixDouble kMatrix1
 
std::vector< MatrixDouble > kMatrixFace
 
std::vector< MatrixDouble > kMatrixFace0
 
std::vector< MatrixDouble > kMatrixFace1
 
- Public Attributes inherited from NitscheMethod::OpCommon
NonlinearElasticElement::BlockDatadAta
 
NonlinearElasticElement::CommonDatacommonData
 
VectorDouble dIsp
 
VectorDouble tRaction
 
MatrixDouble jAc_row
 
MatrixDouble jAc_col
 
MatrixDouble tRac_v
 
MatrixDouble tRac_u
 
- Public Attributes inherited from NitscheMethod::OpBasicCommon
BlockDatanitscheBlockData
 
CommonDatanitscheCommonData
 
bool fieldDisp
 
double faceRadius
 
double gammaH
 
- Public Attributes inherited from MoFEM::ForcesAndSourcesCore::UserDataOperator
char opType
 
std::string rowFieldName
 
std::string colFieldName
 
FieldSpace sPace
 
- Public Attributes inherited from MoFEM::DataOperator
DoWorkLhsHookFunType doWorkLhsHook
 
DoWorkRhsHookFunType doWorkRhsHook
 
bool sYmm
 If true assume that matrix is symmetric structure. More...
 
std::array< bool, MBMAXTYPE > doEntities
 If true operator is executed for entity. More...
 
booldoVertices
 \deprectaed If false skip vertices More...
 
booldoEdges
 \deprectaed If false skip edges More...
 
booldoQuads
 \deprectaed More...
 
booldoTris
 \deprectaed More...
 
booldoTets
 \deprectaed More...
 
booldoPrisms
 \deprectaed More...
 

Additional Inherited Members

- Public Types inherited from MoFEM::ForcesAndSourcesCore::UserDataOperator
enum  OpType { OPROW = 1 << 0 , OPCOL = 1 << 1 , OPROWCOL = 1 << 2 , OPSPACE = 1 << 3 }
 Controls loop over entities on element. More...
 
- Public Types inherited from MoFEM::DataOperator
using DoWorkLhsHookFunType = boost::function< MoFEMErrorCode(DataOperator *op_ptr, int row_side, int col_side, EntityType row_type, EntityType col_type, EntitiesFieldData::EntData &row_data, EntitiesFieldData::EntData &col_data)>
 
using DoWorkRhsHookFunType = boost::function< MoFEMErrorCode(DataOperator *op_ptr, int side, EntityType type, EntitiesFieldData::EntData &data)>
 
- Static Public Attributes inherited from MoFEM::ForcesAndSourcesCore::UserDataOperator
static const char *const OpTypeNames []
 
- Protected Member Functions inherited from MoFEM::VolumeElementForcesAndSourcesCoreBase::UserDataOperator
MoFEMErrorCode setPtrFE (ForcesAndSourcesCore *ptr)
 
- Protected Member Functions inherited from MoFEM::ForcesAndSourcesCore::UserDataOperator
ForcesAndSourcesCoregetPtrFE () const
 
ForcesAndSourcesCoregetSidePtrFE () const
 
ForcesAndSourcesCoregetRefinePtrFE () const
 
- Protected Attributes inherited from MoFEM::ForcesAndSourcesCore::UserDataOperator
ForcesAndSourcesCoreptrFE
 

Detailed Description

Calculate Nitsche method terms on left hand side.

Definition at line 513 of file NitscheMethod.hpp.

Constructor & Destructor Documentation

◆ OpLhsNormal()

NitscheMethod::OpLhsNormal::OpLhsNormal ( const std::string  field_name,
BlockData nitsche_block_data,
CommonData nitsche_common_data,
NonlinearElasticElement::BlockData data,
NonlinearElasticElement::CommonData common_data,
bool  field_disp 
)

Definition at line 515 of file NitscheMethod.hpp.

520  : OpCommon(field_name, nitsche_block_data, nitsche_common_data, data,
521  common_data, field_disp, UserDataOperator::OPROWCOL) {}
OpCommon(const std::string field_name, BlockData &nitsche_block_data, CommonData &nitsche_common_data, NonlinearElasticElement::BlockData &data, NonlinearElasticElement::CommonData &common_data, bool field_disp, const char type)

Member Function Documentation

◆ doWork()

MoFEMErrorCode NitscheMethod::OpLhsNormal::doWork ( int  row_side,
int  col_side,
EntityType  row_type,
EntityType  col_type,
DataForcesAndSourcesCore::EntData row_data,
DataForcesAndSourcesCore::EntData col_data 
)

Definition at line 526 of file NitscheMethod.hpp.

529  {
531 
532  if (dAta.tEts.find(getNumeredEntFiniteElementPtr()->getEnt()) ==
533  dAta.tEts.end()) {
535  }
536  if (row_data.getIndices().size() == 0)
538  if (col_data.getIndices().size() == 0)
540  int nb_dofs_row = row_data.getIndices().size();
541  int nb_dofs_col = col_data.getIndices().size();
542  double gamma = nitscheBlockData.gamma;
543  double phi = nitscheBlockData.phi;
544 
545  kMatrix0.resize(nb_dofs_row, nb_dofs_col, false);
546  kMatrix1.resize(nb_dofs_row, nb_dofs_col, false);
547  kMatrix.resize(nb_dofs_row, nb_dofs_col, false);
548  kMatrix.clear();
549 
550  try {
551 
552  int gg = 0;
553  for (int ff = 0; ff < 4; ff++) {
554  if (!nitscheCommonData.facesFePtr[ff])
555  continue;
556  int nb_face_gauss_pts = nitscheCommonData.faceGaussPts[ff].size2();
557  ierr = getFaceRadius(ff);
558  CHKERRG(ierr);
559  kMatrix0.clear();
560  kMatrix1.clear();
561  for (int fgg = 0; fgg < nb_face_gauss_pts; fgg++, gg++) {
562  ierr = getGammaH(gamma, gg);
563  CHKERRG(ierr);
564  double val = getGaussPts()(3, gg);
565  ierr = getJac(row_data, gg, jAc_row);
566  CHKERRG(ierr);
567  ierr = getTractionVariance(gg, fgg, ff, jAc_row, tRac_v);
568  CHKERRG(ierr);
569  ierr = getJac(col_data, gg, jAc_col);
570  CHKERRG(ierr);
571  ierr = getTractionVariance(gg, fgg, ff, jAc_col, tRac_u);
572  CHKERRG(ierr);
574  3, ublas::shallow_array_adaptor<double>(
575  3, &nitscheCommonData.faceNormals[ff](fgg, 0)));
576  double area = cblas_dnrm2(3, &normal[0], 1);
577 
578  ierr = calculateP(gg, fgg, ff);
579  CHKERRG(ierr);
580  MatrixDouble &P = nitscheCommonData.P[ff][fgg];
581 
582  // P
583  for (int dd1 = 0; dd1 < nb_dofs_row / 3; dd1++) {
584  double n_row = row_data.getN()(gg, dd1);
585  for (int dd2 = 0; dd2 < nb_dofs_col / 3; dd2++) {
586  double n_col = col_data.getN()(gg, dd2);
587  for (int dd3 = 0; dd3 < 3; dd3++) {
588  for (int dd4 = 0; dd4 < 3; dd4++) {
589  if (!P(dd3, dd4))
590  continue;
591  kMatrix0(3 * dd1 + dd3, 3 * dd2 + dd4) +=
592  val * n_row * P(dd3, dd4) * n_col * area;
593  }
594  }
595  }
596  }
597  for (int dd1 = 0; dd1 < nb_dofs_row / 3; dd1++) {
598  double n_row = row_data.getN()(gg, dd1);
599  for (int dd2 = 0; dd2 < nb_dofs_col; dd2++) {
600  for (int dd3 = 0; dd3 < 3; dd3++) {
601  double t = cblas_ddot(3, &P(dd3, 0), 1, &tRac_u(0, dd2),
602  tRac_u.size2());
603  kMatrix1(3 * dd1 + dd3, dd2) -= val * n_row * t;
604  }
605  }
606  }
607  for (int dd1 = 0; dd1 < nb_dofs_row; dd1++) {
608  for (int dd2 = 0; dd2 < nb_dofs_col / 3; dd2++) {
609  double n_col = col_data.getN()(gg, dd2);
610  for (int dd3 = 0; dd3 < 3; dd3++) {
611  double t = cblas_ddot(3, &P(0, dd3), 3, &tRac_v(0, dd1),
612  tRac_v.size2());
613  kMatrix1(dd1, 3 * dd2 + dd3) -= phi * val * t * n_col;
614  }
615  }
616  }
617  // dP
618  if (nitscheCommonData.dP[col_type].empty() != 0) {
619  if (nitscheCommonData.dP[col_type].size() == 4) {
620  if (nitscheCommonData.dP[col_type][ff].size() ==
621  (unsigned int)nb_face_gauss_pts) {
622  MatrixDouble &dP = nitscheCommonData.dP[col_type][ff][fgg];
623  if (dP.size1() == 3 &&
624  dP.size2() == (unsigned int)nb_dofs_col) {
625  VectorDouble &u =
626  (commonData
628  for (int dd1 = 0; dd1 < nb_dofs_row / 3; dd1++) {
629  double n_row = row_data.getN()(gg, dd1);
630  for (int dd2 = 0; dd2 < nb_dofs_col / 3; dd2++) {
631  for (int dd3 = 0; dd3 < 3; dd3++) {
632  for (int dd4 = 0; dd4 < 3; dd4++) {
633  kMatrix0(3 * dd1 + dd3, 3 * dd2 + dd4) +=
634  val * n_row * u[dd4] * dP(dd4, 3 * dd2 + dd4);
635  }
636  }
637  }
638  }
639  VectorDouble t =
640  prod(trans(commonData.sTress[gg]),
641  normal); // this is Piola-Kirchhoff I stress
642  for (int dd1 = 0; dd1 < nb_dofs_row / 3; dd1++) {
643  double n_row = row_data.getN()(gg, dd1);
644  for (int dd2 = 0; dd2 < nb_dofs_col / 3; dd2++) {
645  for (int dd3 = 0; dd3 < 3; dd3++) {
646  for (int dd4 = 0; dd4 < 3; dd4++) {
647  kMatrix1(3 * dd1 + dd3, 3 * dd2 + dd4) +=
648  val * n_row * t[dd4] * dP(dd4, 3 * dd2 + dd4);
649  }
650  }
651  }
652  }
653  }
654  }
655  }
656  }
657  }
658 
659  kMatrix0 /= gammaH;
660  noalias(kMatrix) += kMatrix0 + kMatrix1;
661  }
662 
663  if (gg != (int)getGaussPts().size2()) {
664  SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
665  "wrong number of gauss pts");
666  }
667 
668  ierr = MatSetValues(
669  getFEMethod()->snes_B, nb_dofs_row, &row_data.getIndices()[0],
670  nb_dofs_col, &col_data.getIndices()[0], &kMatrix(0, 0), ADD_VALUES);
671  CHKERRG(ierr);
672 
673  } catch (const std::exception &ex) {
674  std::ostringstream ss;
675  ss << "throw in method: " << ex.what() << std::endl;
676  SETERRQ(PETSC_COMM_SELF, 1, ss.str().c_str());
677  }
678 
680  }
#define MoFEMFunctionReturnHot(a)
Last executable line of each PETSc function used for error handling. Replaces return()
Definition: definitions.h:460
#define CHKERRG(n)
Check error code of MoFEM/MOAB/PETSc function.
Definition: definitions.h:496
@ MOFEM_DATA_INCONSISTENCY
Definition: definitions.h:44
#define MoFEMFunctionBeginHot
First executable line of each MoFEM function, used for error handling. Final line of MoFEM functions ...
Definition: definitions.h:453
FTensor::Tensor1< double, SPACE_DIM > normal(FTensor::Tensor1< T1, 3 > &t_coords, FTensor::Tensor1< T2, SPACE_DIM > &t_disp)
Definition: ContactOps.hpp:77
static MoFEMErrorCodeGeneric< PetscErrorCode > ierr
Definition: Exceptions.hpp:87
VectorShallowArrayAdaptor< double > VectorAdaptor
Definition: Types.hpp:126
UBlasMatrix< double > MatrixDouble
Definition: Types.hpp:88
UBlasVector< double > VectorDouble
Definition: Types.hpp:79
constexpr double t
plate stiffness
Definition: plate.cpp:76
static double phi
boost::shared_ptr< const NumeredEntFiniteElement > getNumeredEntFiniteElementPtr() const
Return raw pointer to NumeredEntFiniteElement.
const FEMethod * getFEMethod() const
Return raw pointer to Finite Element Method object.
MatrixDouble & getGaussPts()
matrix of integration (Gauss) points for Volume Element
FTensor::Tensor2< double *, 3, 3 > & getJac()
get element Jacobian
double gamma
Penalty term, see .
double phi
Nitsche method parameter, see .
std::vector< boost::shared_ptr< const NumeredEntFiniteElement > > facesFePtr
std::vector< MatrixDouble > faceNormals
std::map< EntityType, std::vector< std::vector< MatrixDouble > > > dP
derivative of projection matrix in respect DoFs This is EntityType, face, gauss point at face....
std::vector< std::vector< MatrixDouble > > P
projection matrix
std::vector< MatrixDouble > faceGaussPts
virtual MoFEMErrorCode getFaceRadius(int ff)
virtual MoFEMErrorCode getGammaH(double gamma, int gg)
NonlinearElasticElement::CommonData & commonData
virtual MoFEMErrorCode calculateP(int gg, int fgg, int ff)
NonlinearElasticElement::BlockData & dAta
MoFEMErrorCode getTractionVariance(int gg, int fgg, int ff, MatrixDouble &jac, MatrixDouble &trac)
Range tEts
constrains elements in block set
std::map< std::string, std::vector< VectorDouble > > dataAtGaussPts
std::vector< MatrixDouble3by3 > sTress

Member Data Documentation

◆ kMatrix

MatrixDouble NitscheMethod::OpLhsNormal::kMatrix

Definition at line 523 of file NitscheMethod.hpp.

◆ kMatrix0

MatrixDouble NitscheMethod::OpLhsNormal::kMatrix0

Definition at line 523 of file NitscheMethod.hpp.

◆ kMatrix1

MatrixDouble NitscheMethod::OpLhsNormal::kMatrix1

Definition at line 523 of file NitscheMethod.hpp.

◆ kMatrixFace

std::vector<MatrixDouble> NitscheMethod::OpLhsNormal::kMatrixFace

Definition at line 524 of file NitscheMethod.hpp.

◆ kMatrixFace0

std::vector<MatrixDouble> NitscheMethod::OpLhsNormal::kMatrixFace0

Definition at line 524 of file NitscheMethod.hpp.

◆ kMatrixFace1

std::vector<MatrixDouble> NitscheMethod::OpLhsNormal::kMatrixFace1

Definition at line 524 of file NitscheMethod.hpp.


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