SeepageOps::OpDomainRhsHydrostaticStress< DIM > Struct Template Reference

#include "tutorials/adv-5_poroelasticity/src/SeepageOps.hpp"

Inheritance diagram for SeepageOps::OpDomainRhsHydrostaticStress< DIM >:

Collaboration diagram for SeepageOps::OpDomainRhsHydrostaticStress< DIM >:

Public Member Functions
	OpDomainRhsHydrostaticStress (std::string field_name1, boost::shared_ptr< VectorDouble > p_ptr, std::function< double(const double, const double, const double)> multiplier_func)
MoFEMErrorCode	iNtegrate (DataForcesAndSourcesCore::EntData &data)

Private Attributes
boost::shared_ptr< VectorDouble >	pPtr
std::function< double(const double, const double, const double)>	multiplier_func

Detailed Description

template<int DIM>
struct SeepageOps::OpDomainRhsHydrostaticStress< DIM >

Examples

mofem/tutorials/adv-5_poroelasticity/seepage.cpp.

Definition at line 22 of file SeepageOps.hpp.

Constructor & Destructor Documentation

OpDomainRhsHydrostaticStress()

template<int DIM>

SeepageOps::OpDomainRhsHydrostaticStress< DIM >::OpDomainRhsHydrostaticStress ( std::string  field_name1, boost::shared_ptr< VectorDouble >  p_ptr, std::function< double(const double, const double, const double)>  multiplier_func  )

inline

Definition at line 25 of file SeepageOps.hpp.

      : AssemblyDomainEleOp(field_name1, field_name1, DomainEleOp::OPROW),
        pPtr(p_ptr), multiplier_func(multiplier_func) {}

Member Function Documentation

iNtegrate()

template<int DIM>

MoFEMErrorCode SeepageOps::OpDomainRhsHydrostaticStress< DIM >::iNtegrate ( DataForcesAndSourcesCore::EntData &  data )

inline

Examples

mofem/tutorials/adv-5_poroelasticity/src/SeepageOps.hpp.

Definition at line 31 of file SeepageOps.hpp.

                                                                  {
    MoFEMFunctionBegin;
 
    const int nb_dofs = data.getIndices().size();
 
    if (nb_dofs) {
      
      auto &nf = AssemblyDomainEleOp::locF; //access to local vector
      // get element measure
      const double measure = getMeasure();
 
      // get number of integration points
      const int nb_integration_points = getGaussPts().size2();
      // get integration weights
      auto t_w = getFTensor0IntegrationWeight();
 
      // get base function
      auto t_base_diff = data.getFTensor1DiffN<DIM>();
 
      //constexpr double g_acceleration = 9.81;
 
      FTensor::Index<'i', DIM> i;
 
      auto t_p = getFTensor0FromVec(*pPtr);
      for (int gg = 0; gg != nb_integration_points; gg++) {
        
        auto t_nf = getFTensor1FromPtr<DIM>(&nf[0]);
 
        const double biot = multiplier_func(0.0, 0.0, 0.0);
        const double a = t_w * measure * t_p * biot;
 
        for (int rr = 0; rr != nb_dofs / DIM; rr++) {  //loop over dofs (rr = alpha)
          t_nf(i) -= t_base_diff(i) * a;  //loop over dimensions
 
          // move to the next base function
          ++t_base_diff; // moves the pointer to the next shape function
          ++t_nf;
        }
 
        // move to the weight of the next integration point
        ++t_w;
        ++t_p;
      }
   
    }
 
    MoFEMFunctionReturn(0);
  }

Member Data Documentation

multiplier_func

template<int DIM>

std::function<double(const double, const double, const double)> SeepageOps::OpDomainRhsHydrostaticStress< DIM >::multiplier_func

private

Examples

mofem/tutorials/adv-5_poroelasticity/src/SeepageOps.hpp.

Definition at line 82 of file SeepageOps.hpp.

pPtr

template<int DIM>

boost::shared_ptr<VectorDouble> SeepageOps::OpDomainRhsHydrostaticStress< DIM >::pPtr

private

Examples

mofem/tutorials/adv-5_poroelasticity/src/SeepageOps.hpp.

Definition at line 81 of file SeepageOps.hpp.

---

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

• tutorials/adv-5_poroelasticity/src/SeepageOps.hpp