#include <MoFEM.hpp>

Classes
struct	ElementsAndOps< DIM >

struct	ElementsAndOps< 2 >

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

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

struct	Plate

Typedefs
using	DomainEle = ElementsAndOps< SPACE_DIM >::DomainEle

using	DomainEleOp = DomainEle::UserDataOperator
	Finire element operator type.

using	BoundaryEle = ElementsAndOps< SPACE_DIM >::BoundaryEle

using	BoundaryEleOp = BoundaryEle::UserDataOperator

using	PostProcEle = PostProcBrokenMeshInMoab< DomainEle >

using	FaceSideEle = ElementsAndOps< SPACE_DIM >::FaceSideEle

using	FaceSideOp = ElementsAndOps< SPACE_DIM >::FaceSideOp

using	EntData = EntitiesFieldData::EntData
	Data on entities.

using	OpDomainGradGrad = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >::OpGradGrad< BASE_DIM, FIELD_DIM, SPACE_DIM >

using	OpDomainPlateStiffness = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::BiLinearForm< GAUSS >::OpGradGradSymTensorGradGrad< 1, 1, SPACE_DIM, 0 >

using	OpDomainPlateLoad = FormsIntegrators< DomainEleOp >::Assembly< PETSC >::LinearForm< GAUSS >::OpSource< BASE_DIM, FIELD_DIM >

Enumerations
enum	ElementSide { LEFT_SIDE = 0 , RIGHT_SIDE }

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

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)

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

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

Variables
static char	help [] = "...\n\n"

constexpr int	BASE_DIM = 1
	dimension of base

constexpr int	SPACE_DIM = 2
	dimension of space

constexpr int	FIELD_DIM = 1
	dimension of approx. field

constexpr auto	t_kd = FTensor::Kronecker_Delta_symmetric<int>()

constexpr double	lambda = 1

constexpr double	mu = 1
	lame parameter

constexpr double	t = 1
	plate stiffness

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

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

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

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

static double	penalty = 1e6

static double	phi

static double	nitsche = 1

static int	order = 4

auto	source

auto	plate_stiffness
	get fourth-order constitutive tensor

std::array< std::vector< VectorInt >, 2 >	indicesSideMap
	indices on rows for left hand-side

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

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

std::array< double, 2 >	areaMap

std::array< int, 2 >	senseMap

Typedef Documentation

◆ BoundaryEle

using BoundaryEle = ElementsAndOps<SPACE_DIM>::BoundaryEle

Definition at line 31 of file plate.cpp.

◆ BoundaryEleOp

using BoundaryEleOp = BoundaryEle::UserDataOperator

Definition at line 32 of file plate.cpp.

◆ DomainEle

using DomainEle = ElementsAndOps<SPACE_DIM>::DomainEle

Definition at line 29 of file plate.cpp.

◆ DomainEleOp

using DomainEleOp = DomainEle::UserDataOperator

Finire element operator type.

Definition at line 30 of file plate.cpp.

◆ EntData

using EntData = EntitiesFieldData::EntData

Data on entities.

Definition at line 41 of file plate.cpp.

◆ FaceSideEle

using FaceSideEle = ElementsAndOps<SPACE_DIM>::FaceSideEle

Definition at line 36 of file plate.cpp.

◆ FaceSideOp

using FaceSideOp = ElementsAndOps<SPACE_DIM>::FaceSideOp

Definition at line 37 of file plate.cpp.

◆ OpDomainGradGrad

using OpDomainGradGrad = FormsIntegrators<DomainEleOp>::Assembly< PETSC>::BiLinearForm<GAUSS>::OpGradGrad<BASE_DIM, FIELD_DIM, SPACE_DIM>

Definition at line 43 of file plate.cpp.

◆ OpDomainPlateLoad

using OpDomainPlateLoad = FormsIntegrators<DomainEleOp>::Assembly< PETSC>::LinearForm<GAUSS>::OpSource<BASE_DIM, FIELD_DIM>

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 49 of file plate.cpp.

◆ OpDomainPlateStiffness

using OpDomainPlateStiffness = FormsIntegrators<DomainEleOp>::Assembly<PETSC>::BiLinearForm< GAUSS>::OpGradGradSymTensorGradGrad<1, 1, SPACE_DIM, 0>

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 46 of file plate.cpp.

◆ PostProcEle

using PostProcEle = PostProcBrokenMeshInMoab<DomainEle>

Definition at line 34 of file plate.cpp.

Enumeration Type Documentation

◆ ElementSide

enum ElementSide

Enumerator
LEFT_SIDE
RIGHT_SIDE

Definition at line 92 of file plate.cpp.

92{ LEFT_SIDE = 0, RIGHT_SIDE };

LEFT_SIDE

@ LEFT_SIDE

Definition plate.cpp:92

RIGHT_SIDE

@ RIGHT_SIDE

Definition plate.cpp:92

Function Documentation

◆ get_diff2_ntensor() [1/2]

template<typename T >

auto get_diff2_ntensor ( T & base_mat )

inline

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 489 of file plate.cpp.

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

◆ get_diff2_ntensor() [2/2]

template<typename T >

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

inline

Definition at line 495 of file plate.cpp.

                                                           {
  double *ptr = &base_mat(gg, 4 * bb);
  return getFTensor2SymmetricLowerFromPtr<2>(ptr);
};

◆ get_diff_ntensor() [1/2]

template<typename T >

auto get_diff_ntensor ( T & base_mat )

inline

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 478 of file plate.cpp.

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

◆ get_diff_ntensor() [2/2]

template<typename T >

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

inline

Definition at line 484 of file plate.cpp.

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

◆ get_ntensor() [1/2]

template<typename T >

auto get_ntensor ( T & base_mat )

inline

Examples: mofem/tutorials/adv-3/level_set.cpp, and mofem/tutorials/vec-6/plate.cpp.

Definition at line 468 of file plate.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 473 of file plate.cpp.

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

◆ main()

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

[Run program]

[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]

[Plate]

Definition at line 390 of file plate.cpp.

                                 {
 
  // Initialisation of MoFEM/PETSc and MOAB data structures
  const char param_file[] = "param_file.petsc";
  MoFEM::Core::Initialize(&argc, &argv, param_file, help);
 
  // Add logging channel for example
  auto core_log = logging::core::get();
  core_log->add_sink(LogManager::createSink(LogManager::getStrmWorld(), "PL"));
  LogManager::setLog("PL");
  MOFEM_LOG_TAG("PL", "plate");
 
  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]
 
    //! [Plate]
    Plate ex(m_field);
    CHKERR ex.runProblem();
    //! [Plate]
  }
  CATCH_ERRORS;
 
  CHKERR MoFEM::Core::Finalize();
}

Variable Documentation

◆ areaMap

std::array<double, 2> areaMap

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 101 of file plate.cpp.

◆ BASE_DIM

constexpr int BASE_DIM = 1

constexpr

dimension of base

Definition at line 16 of file plate.cpp.

◆ diff2BaseSideMap

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

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 100 of file plate.cpp.

◆ diffBaseSideMap

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

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 98 of file plate.cpp.

◆ FIELD_DIM

constexpr int FIELD_DIM = 1

constexpr

dimension of approx. field

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 18 of file plate.cpp.

◆ help

char help[] = "...\n\n"

static

Definition at line 13 of file plate.cpp.

◆ i

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

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 60 of file plate.cpp.

◆ indicesSideMap

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

indices on rows for left hand-side

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 96 of file plate.cpp.

◆ j

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

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 61 of file plate.cpp.

◆ k

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

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 62 of file plate.cpp.

◆ l

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

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 63 of file plate.cpp.

◆ lambda

constexpr double lambda = 1

constexpr

Definition at line 56 of file plate.cpp.

◆ mu

constexpr double mu = 1

constexpr

lame parameter

Definition at line 57 of file plate.cpp.

◆ nitsche

double nitsche = 1

static

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 68 of file plate.cpp.

◆ order

int order = 4

static

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 69 of file plate.cpp.

◆ penalty

double penalty = 1e6

static

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 65 of file plate.cpp.

◆ phi

double phi

static

Initial value:

=

-1

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 66 of file plate.cpp.

◆ plate_stiffness

auto plate_stiffness

Initial value:

= []() {
  constexpr auto a = (SPACE_DIM * (SPACE_DIM + 1)) / 2;
  auto mat_D_ptr = boost::make_shared<MatrixDouble>(a * a, 1);
  auto t_D = getFTensor4DdgFromMat<2, 2, 0>(*(mat_D_ptr));
  constexpr double t3 = t * t * t;
  constexpr double A = mu * t3 / 12;
  constexpr double B = lambda * t3 / 12;
  t_D(i, j, k, l) =
      2 * B * ((t_kd(i, k) ^ t_kd(j, l)) / 4.) + A * t_kd(i, j) * t_kd(k, l);
  
  return mat_D_ptr;
}

get fourth-order constitutive tensor

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 79 of file plate.cpp.

                            {
  constexpr auto a = (SPACE_DIM * (SPACE_DIM + 1)) / 2;
  auto mat_D_ptr = boost::make_shared<MatrixDouble>(a * a, 1);
  auto t_D = getFTensor4DdgFromMat<2, 2, 0>(*(mat_D_ptr));
  constexpr double t3 = t * t * t;
  constexpr double A = mu * t3 / 12;
  constexpr double B = lambda * t3 / 12;
  t_D(i, j, k, l) =
      2 * B * ((t_kd(i, k) ^ t_kd(j, l)) / 4.) + A * t_kd(i, j) * t_kd(k, l);
  // t_D(i, j, k, l) = (t_kd(i, k) ^ t_kd(j, l)) / 4.;
  return mat_D_ptr;
};

◆ senseMap

std::array<int, 2> senseMap

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 102 of file plate.cpp.

◆ source

auto source

Initial value:

= [](const double x, const double y, const double) {
  return cos(2 * x * M_PI) * sin(2 * y * M_PI);
}

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 71 of file plate.cpp.

                                                               {
  return cos(2 * x * M_PI) * sin(2 * y * M_PI);
};

◆ SPACE_DIM

constexpr int SPACE_DIM = 2

constexpr

dimension of space

Examples: mofem/tutorials/vec-6/plate.cpp.

Definition at line 17 of file plate.cpp.

◆ t

constexpr double t = 1

constexpr

plate stiffness

Examples: mofem/atom_tests/gauss_points_on_outer_product.cpp, mofem/atom_tests/hanging_node_approx.cpp, mofem/atom_tests/hcurl_divergence_operator_2d.cpp, mofem/atom_tests/matrix_function.cpp, mofem/atom_tests/prism_elements_from_surface.cpp, mofem/atom_tests/remove_entities_from_problem.cpp, mofem/atom_tests/remove_entities_from_problem_not_partitioned.cpp, mofem/tutorials/adv-0/plastic.cpp, mofem/tutorials/adv-1/src/ContactOps.hpp, mofem/tutorials/adv-3/level_set.cpp, mofem/tutorials/max-1/lorentz_force.cpp, mofem/tutorials/scl-11/src/PoissonDiscontinousGalerkin.hpp, mofem/tutorials/scl-6/heat_equation.cpp, mofem/tutorials/scl-7/wave_equation.cpp, mofem/tutorials/vec-5/free_surface.cpp, mofem/tutorials/vec-6/plate.cpp, mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianOperators.cpp, and mofem/users_modules/eshelbian_plasticity/src/impl/EshelbianPlasticity.cpp.

Definition at line 58 of file plate.cpp.

◆ t_kd

constexpr auto t_kd = FTensor::Kronecker_Delta_symmetric<int>()

constexpr

Definition at line 53 of file plate.cpp.

Classes

Typedefs

Enumerations

Functions

Variables

Typedef Documentation

◆ BoundaryEle

◆ BoundaryEleOp

◆ DomainEle

◆ DomainEleOp

◆ EntData

◆ FaceSideEle

◆ FaceSideOp

◆ OpDomainGradGrad

◆ OpDomainPlateLoad

◆ OpDomainPlateStiffness

◆ PostProcEle

Enumeration Type Documentation

◆ ElementSide

Function Documentation

◆ get_diff2_ntensor() [1/2]

◆ get_diff2_ntensor() [2/2]

◆ get_diff_ntensor() [1/2]

◆ get_diff_ntensor() [2/2]

◆ get_ntensor() [1/2]

◆ get_ntensor() [2/2]

◆ main()

Variable Documentation

◆ areaMap

◆ BASE_DIM

◆ diff2BaseSideMap

◆ diffBaseSideMap

◆ FIELD_DIM

◆ help

◆ i

◆ indicesSideMap

◆ j

◆ k

◆ l

◆ lambda

◆ mu

◆ nitsche

◆ order

◆ penalty

◆ phi

◆ plate_stiffness

◆ senseMap

◆ source

◆ SPACE_DIM

◆ t

◆ t_kd