CohesiveNaturalBoundaryBC.hpp

Go to the documentation of this file.

/**
 * @file CohesiveNaturalBoundaryBC.hpp
 * @brief Natural force and pressure boundary conditions
 */
 
#ifndef COHESIVE_INTERFACE_NATURAL_BOUNDARY_BC_HPP
#define COHESIVE_INTERFACE_NATURAL_BOUNDARY_BC_HPP
 
template <int BASE_DIM, int FIELD_DIM, AssemblyType A, IntegrationType I,
          typename OpBase>
struct AddFluxToRhsPipelineImpl<
 
    OpFluxRhsImpl<BoundaryBCs, BASE_DIM, FIELD_DIM, A, I, OpBase>, A, I, OpBase
 
    > {
 
  AddFluxToRhsPipelineImpl() = delete;
 
  using T =
      typename NaturalBC<OpBase>::template Assembly<A>::template LinearForm<I>;
  using OpForce =
      typename T::template OpFlux<NaturalForceMeshsets, 1, SPACE_DIM>;
 
  static MoFEMErrorCode add(
 
      boost::ptr_deque<ForcesAndSourcesCore::UserDataOperator> &pipeline,
      MoFEM::Interface &m_field, const std::string field_name,
      std::vector<boost::shared_ptr<ScalingMethod>> scaling_methods, Sev sev
 
  ) {
    MoFEMFunctionBegin;
    CHKERR T::template AddFluxToPipeline<OpForce>::add(
        pipeline, m_field, field_name, scaling_methods, "FORCE", "PRESSURE",
        sev);
    MoFEMFunctionReturn(0);
  }
 
  static MoFEMErrorCode add(
      boost::ptr_deque<ForcesAndSourcesCore::UserDataOperator> &pipeline,
      MoFEM::Interface &m_field, const std::string field_name, Sev sev) {
    return add(pipeline, m_field, field_name, {}, sev);
  }
};
 
#endif // COHESIVE_INTERFACE_NATURAL_BOUNDARY_BC_HPP