#include <MoFEM.hpp>
#include <slepceps.h>
#include <slepcpep.h>
#include <OpPostProcElastic.hpp>

Classes
struct	ElementsAndOps< DIM >

struct	ElementsAndOps< 2 >

struct	ElementsAndOps< 3 >

struct	OpCalculateSideData
	Operator tp collect data from elements on the side of Edge/Face. More...

struct	OpLhsSkeleton
	Operator the left hand side matrix. More...

struct	FSI

Typedefs
using	EntData = EntitiesFieldData::EntData

using	DomainEle = ElementsAndOps< SPACE_DIM >::DomainEle

using	DomainEleOp = ElementsAndOps< SPACE_DIM >::DomainEleOp

using	PostProcEle = ElementsAndOps< SPACE_DIM >::PostProcEle

using	FaceSideEle = ElementsAndOps< SPACE_DIM >::FaceSideEle

using	FaceSideOp = ElementsAndOps< SPACE_DIM >::FaceSideOp

using	SkeletonEle = ElementsAndOps< SPACE_DIM >::SkeletonEle

using	SkeletonEleOp = ElementsAndOps< SPACE_DIM >::SkeletonEleOp

using	OpK = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >::OpGradSymTensorGrad< 1, SPACE_DIM, SPACE_DIM, 0 >

using	OpMass = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >::OpMass< 1, SPACE_DIM >

Enumerations
enum	ElementSide { LEFT_SIDE = 0 , RIGHT_SIDE , LEFT_SIDE = 0 , RIGHT_SIDE }

Functions
int	main (int argc, char *argv[])

template<typename T >
auto	get_ntensor (T &base_mat)

template<typename T >
auto	get_ntensor (T &base_mat, int gg, int bb)

template<typename T >
auto	get_diff_ntensor (T &base_mat)

template<typename T >
auto	get_diff_ntensor (T &base_mat, int gg, int bb)

Variables
constexpr int	SPACE_DIM

FTensor::Index< 'i', SPACE_DIM >	i

FTensor::Index< 'j', SPACE_DIM >	j

FTensor::Index< 'k', SPACE_DIM >	k

FTensor::Index< 'l', SPACE_DIM >	l

FTensor::Index< 'm', SPACE_DIM >	m

FTensor::Index< 'n', SPACE_DIM >	n

static double	penalty = 1e2

static double	phi = -1

double	rho_solid = 2.28

double	K_solid = 95 * 1e3

double	G_solid = 63 * 1e3

double	rho_fluid = 0.998

double	K_fluid = 2.2 * 1e3

double	G_fluid = 1e-6

double	MU_fluid = 1e-3

int	order = 1

double	eig_threshold = 1e-1

auto	poisson_ratio = 0.265

static char	help [] = "...\n\n"
	[Postprocess results] More...

std::array< std::vector< VectorInt >, 2 >	indicesUsSideMap

std::array< std::vector< VectorInt >, 2 >	indicesUfSideMap

std::array< std::vector< MatrixDouble >, 2 >	baseUfSideMap

std::array< std::vector< MatrixDouble >, 2 >	diffUfBaseSideMap

std::array< std::vector< MatrixDouble >, 2 >	baseUsSideMap

std::array< std::vector< MatrixDouble >, 2 >	diffUsBaseSideMap

std::array< EntityHandle, 2 >	sideHandle

std::array< double, 2 >	areaMap

std::array< int, 2 >	senseMap

Typedef Documentation

◆ DomainEle

using DomainEle = ElementsAndOps<SPACE_DIM>::DomainEle

Definition at line 46 of file fluid_structure_eigenproblem.cpp.

◆ DomainEleOp

using DomainEleOp = ElementsAndOps<SPACE_DIM>::DomainEleOp

Definition at line 47 of file fluid_structure_eigenproblem.cpp.

◆ EntData

using EntData = EntitiesFieldData::EntData

Definition at line 45 of file fluid_structure_eigenproblem.cpp.

◆ FaceSideEle

using FaceSideEle = ElementsAndOps<SPACE_DIM>::FaceSideEle

Definition at line 50 of file fluid_structure_eigenproblem.cpp.

◆ FaceSideOp

using FaceSideOp = ElementsAndOps<SPACE_DIM>::FaceSideOp

Examples: PoissonDiscontinousGalerkin.hpp.

Definition at line 51 of file fluid_structure_eigenproblem.cpp.

◆ OpK

using OpK = FormsIntegrators<DomainEleOp>::Assembly<PETSC>::BiLinearForm< GAUSS>::OpGradSymTensorGrad<1, SPACE_DIM, SPACE_DIM, 0>

Definition at line 56 of file fluid_structure_eigenproblem.cpp.

◆ OpMass

using OpMass = FormsIntegrators<DomainEleOp>::Assembly<PETSC>::BiLinearForm< GAUSS>::OpMass<1, SPACE_DIM>

Examples: eigen_elastic.cpp, heat_method.cpp, plastic.cpp, and scalar_check_approximation.cpp.

Definition at line 58 of file fluid_structure_eigenproblem.cpp.

◆ PostProcEle

using PostProcEle = ElementsAndOps<SPACE_DIM>::PostProcEle

Definition at line 48 of file fluid_structure_eigenproblem.cpp.

◆ SkeletonEle

using SkeletonEle = ElementsAndOps<SPACE_DIM>::SkeletonEle

Definition at line 53 of file fluid_structure_eigenproblem.cpp.

◆ SkeletonEleOp

using SkeletonEleOp = ElementsAndOps<SPACE_DIM>::SkeletonEleOp

Definition at line 54 of file fluid_structure_eigenproblem.cpp.

Enumeration Type Documentation

◆ ElementSide

enum ElementSide

Enumerator
LEFT_SIDE
RIGHT_SIDE
LEFT_SIDE
RIGHT_SIDE

Definition at line 1014 of file fluid_structure_eigenproblem.cpp.

1014{ LEFT_SIDE = 0, RIGHT_SIDE };

LEFT_SIDE

@ LEFT_SIDE

Definition: fluid_structure_eigenproblem.cpp:1014

RIGHT_SIDE

@ RIGHT_SIDE

Definition: fluid_structure_eigenproblem.cpp:1014

Function Documentation

◆ get_diff_ntensor() [1/2]

template<typename T >

auto get_diff_ntensor ( T & base_mat )

inline

Examples: plate.cpp.

Definition at line 1026 of file fluid_structure_eigenproblem.cpp.

                                                                {
  double *ptr = &*base_mat.data().begin();
  return getFTensor1FromPtr<SPACE_DIM>(ptr);
};

◆ get_diff_ntensor() [2/2]

template<typename T >

auto get_diff_ntensor	(	T &	base_mat,
		int	gg,
		int	bb
	)

inline

Definition at line 1032 of file fluid_structure_eigenproblem.cpp.

                                                          {
  double *ptr = &base_mat(gg, SPACE_DIM * bb);
  return getFTensor1FromPtr<SPACE_DIM>(ptr);
};

◆ get_ntensor() [1/2]

template<typename T >

auto get_ntensor ( T & base_mat )

inline

Examples: level_set.cpp, and plate.cpp.

Definition at line 1016 of file fluid_structure_eigenproblem.cpp.

                                                           {
  return FTensor::Tensor0<FTensor::PackPtr<double *, 1>>(
      &*base_mat.data().begin());
};

◆ get_ntensor() [2/2]

template<typename T >

auto get_ntensor	(	T &	base_mat,
		int	gg,
		int	bb
	)

inline

Definition at line 1021 of file fluid_structure_eigenproblem.cpp.

                                                                           {
  double *ptr = &base_mat(gg, bb);
  return FTensor::Tensor0<FTensor::PackPtr<double *, 1>>(ptr);
};

◆ main()

int main	(	int	argc,
		char *	argv[]
	)

[Register MoFEM discrete manager in PETSc]

[Register MoFEM discrete manager in PETSc

[Create MoAB]

< mesh database

< mesh database interface

[Create MoAB]

[Create MoFEM]

< finite element database

< finite element database insterface

[Create MoFEM]

[FSI]

Definition at line 882 of file fluid_structure_eigenproblem.cpp.

                                 {
 
  // Initialisation of MoFEM/PETSc and MOAB data structures
  const char param_file[] = "param_file.petsc";
  SlepcInitialize(&argc, &argv, param_file, help);
  MoFEM::Core::Initialize(&argc, &argv, param_file, help);
 
  // Add logging channel for FSI
  auto core_log = logging::core::get();
  core_log->add_sink(LogManager::createSink(LogManager::getStrmWorld(), "FSI"));
  LogManager::setLog("FSI");
  MOFEM_LOG_TAG("FSI", "FSI");
 
  try {
 
    //! [Register MoFEM discrete manager in PETSc]
    DMType dm_name = "DMMOFEM";
    CHKERR DMRegister_MoFEM(dm_name);
    //! [Register MoFEM discrete manager in PETSc
 
    //! [Create MoAB]
    moab::Core mb_instance;              ///< mesh database
    moab::Interface &moab = mb_instance; ///< mesh database interface
    //! [Create MoAB]
 
    //! [Create MoFEM]
    MoFEM::Core core(moab);           ///< finite element database
    MoFEM::Interface &m_field = core; ///< finite element database insterface
    //! [Create MoFEM]
 
    //! [FSI]
    FSI ex(m_field);
    CHKERR ex.runProblem();
    //! [FSI]
  }
  CATCH_ERRORS;
 
  SlepcFinalize();
  CHKERR MoFEM::Core::Finalize();
}