dynamic_first_order_con_law.cpp File Reference

#include <MoFEM.hpp>
#include <MatrixFunction.hpp>

Go to the source code of this file.

Classes
struct	PostProcEleByDim< 2 >
struct	PostProcEleByDim< 3 >
struct	OpCalculateFStab< DIM_0, DIM_1 >
struct	OpCalculatePiola< DIM_0, DIM_1 >
struct	OpCalculateDisplacement< DIM >
struct	OpCalculatePiolaIncompressibleNH< DIM_0, DIM_1 >
struct	OpCalculateDeformationGradient< DIM_0, DIM_1 >
struct	TSPrePostProc
	Set of functions called by PETSc solver used to refine and update mesh. More...
struct	LinMomTimeScale
struct	CommonData
struct	Example
	[Example] More...
struct	Example::DynamicFirstOrderConsSinusTimeScale
struct	Example::DynamicFirstOrderConsConstantTimeScale
struct	Monitor
	[Push operators to pipeline] More...

Typedefs
using	DomainEle = PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::DomainEle
using	PostProcEle = PostProcBrokenMeshInMoab<DomainEle>
using	PostProcFaceEle
using	BoundaryEle
using	PostProcEleDomain = PostProcEleByDim<SPACE_DIM>::PostProcEleDomain
using	SideEle = PostProcEleByDim<SPACE_DIM>::SideEle
using	PostProcEleBdy = PostProcEleByDim<SPACE_DIM>::PostProcEleBdy
using	OpMassV
using	OpMassF
using	OpInertiaForce
using	OpBodyForce
using	DomainNaturalBC
using	OpBodyForceVector
using	OpGradTimesTensor2
using	OpRhsTestPiola
using	OpGradTimesPiola
using	BoundaryNaturalBC
using	OpForce = BoundaryNaturalBC::OpFlux<NaturalForceMeshsets, 1, SPACE_DIM>

Functions
template<typename T >
double	trace (FTensor::Tensor2< T, 2, 2 > &t_stress)
template<typename T >
double	trace (FTensor::Tensor2< T, 3, 3 > &t_stress)
int	main (int argc, char *argv[])

Variables
constexpr int	SPACE_DIM
constexpr double	omega = 1.
	Save field DOFS on vertices/tags.
constexpr double	young_modulus = 1.
constexpr double	poisson_ratio = 0.
double	bulk_modulus_K = young_modulus / (3. * (1. - 2. * poisson_ratio))
double	shear_modulus_G = young_modulus / (2. * (1. + poisson_ratio))
double	mu = young_modulus / (2. * (1. + poisson_ratio))
double	lamme_lambda
static boost::weak_ptr< TSPrePostProc >	tsPrePostProc
static char	help [] = "...\n\n"
	[Check]

Typedef Documentation

BoundaryEle

using BoundaryEle

Initial value:

 
    PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::BoundaryEle

Definition at line 33 of file dynamic_first_order_con_law.cpp.

BoundaryNaturalBC

using BoundaryNaturalBC

Initial value:

 
    NaturalBC<BoundaryEleOp>::Assembly<PETSC>::LinearForm<GAUSS>

Definition at line 85 of file dynamic_first_order_con_law.cpp.

DomainEle

using DomainEle = PipelineManager::ElementsAndOpsByDim<SPACE_DIM>::DomainEle

Definition at line 27 of file dynamic_first_order_con_law.cpp.

DomainNaturalBC

using DomainNaturalBC

Initial value:

 
    NaturalBC<DomainEleOp>::Assembly<PETSC>::LinearForm<GAUSS>

Definition at line 67 of file dynamic_first_order_con_law.cpp.

OpBodyForce

using OpBodyForce

Initial value:

 FormsIntegrators<DomainEleOp>::Assembly<PETSC>::LinearForm<
    GAUSS>::OpBaseTimesVector<1, SPACE_DIM, 0>

Definition at line 64 of file dynamic_first_order_con_law.cpp.

OpBodyForceVector

using OpBodyForceVector

Initial value:

 
    DomainNaturalBC::OpFlux<NaturalMeshsetTypeVectorScaling<BLOCKSET>, 1,
                            SPACE_DIM>

Definition at line 69 of file dynamic_first_order_con_law.cpp.

OpForce

using OpForce = BoundaryNaturalBC::OpFlux<NaturalForceMeshsets, 1, SPACE_DIM>

Definition at line 87 of file dynamic_first_order_con_law.cpp.

OpGradTimesPiola

using OpGradTimesPiola

Initial value:

 
    FormsIntegrators<DomainEleOp>::Assembly<AssemblyType::PETSC>::LinearForm<
        IntegrationType::GAUSS>::OpGradTimesTensor<1, SPACE_DIM, SPACE_DIM>

Definition at line 81 of file dynamic_first_order_con_law.cpp.

OpGradTimesTensor2

using OpGradTimesTensor2

Initial value:

 
    FormsIntegrators<DomainEleOp>::Assembly<AssemblyType::PETSC>::LinearForm<
        IntegrationType::GAUSS>::OpGradTimesTensor<1, SPACE_DIM, SPACE_DIM>

Examples

dynamic_first_order_con_law.cpp.

Definition at line 73 of file dynamic_first_order_con_law.cpp.

OpInertiaForce

using OpInertiaForce

Initial value:

 FormsIntegrators<DomainEleOp>::Assembly<
    PETSC>::LinearForm<GAUSS>::OpBaseTimesVector<1, SPACE_DIM, 1>

Examples

plastic.cpp.

Definition at line 61 of file dynamic_first_order_con_law.cpp.

OpMassF

using OpMassF

Initial value:

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

Examples

dynamic_first_order_con_law.cpp.

Definition at line 58 of file dynamic_first_order_con_law.cpp.

OpMassV

using OpMassV

Initial value:

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

Examples

dynamic_first_order_con_law.cpp.

Definition at line 56 of file dynamic_first_order_con_law.cpp.

OpRhsTestPiola

using OpRhsTestPiola

Initial value:

 
    FormsIntegrators<DomainEleOp>::Assembly<AssemblyType::PETSC>::LinearForm<
        IntegrationType::GAUSS>::OpBaseTimesVector<1, SPACE_DIM * SPACE_DIM, 1>

Examples

dynamic_first_order_con_law.cpp.

Definition at line 77 of file dynamic_first_order_con_law.cpp.

PostProcEle

using PostProcEle = PostProcBrokenMeshInMoab<DomainEle>

Definition at line 29 of file dynamic_first_order_con_law.cpp.

PostProcEleBdy

using PostProcEleBdy = PostProcEleByDim<SPACE_DIM>::PostProcEleBdy

Definition at line 54 of file dynamic_first_order_con_law.cpp.

PostProcEleDomain

using PostProcEleDomain = PostProcEleByDim<SPACE_DIM>::PostProcEleDomain

Definition at line 52 of file dynamic_first_order_con_law.cpp.

PostProcFaceEle

using PostProcFaceEle

Initial value:

 
    PostProcBrokenMeshInMoab<FaceElementForcesAndSourcesCore>

Definition at line 30 of file dynamic_first_order_con_law.cpp.

SideEle

using SideEle = PostProcEleByDim<SPACE_DIM>::SideEle

Definition at line 53 of file dynamic_first_order_con_law.cpp.

Function Documentation

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 interface

[Create MoFEM]

[Example]

[Example]

Definition at line 1190 of file dynamic_first_order_con_law.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(), "EXAMPLE"));
  LogManager::setLog("EXAMPLE");
  MOFEM_LOG_TAG("EXAMPLE", "example");
 
  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 interface
    //! [Create MoFEM]
 
    //! [Example]
    Example ex(m_field);
    CHKERR ex.runProblem();
    //! [Example]
  }
  CATCH_ERRORS;
 
  CHKERR MoFEM::Core::Finalize();
}

trace() [1/2]

template<typename T >

double trace ( FTensor::Tensor2< T, 2, 2 > & t_stress )

inline

Examples

dynamic_first_order_con_law.cpp.

Definition at line 14 of file dynamic_first_order_con_law.cpp.

                                                                             {
  constexpr double third = boost::math::constants::third<double>();
  return (t_stress(0, 0) + t_stress(1, 1));
};

trace() [2/2]

template<typename T >

double trace ( FTensor::Tensor2< T, 3, 3 > & t_stress )

inline

Definition at line 19 of file dynamic_first_order_con_law.cpp.

                                                                             {
  return (t_stress(0, 0) + t_stress(1, 1) + t_stress(2, 2));
};

Variable Documentation

bulk_modulus_K

double bulk_modulus_K = young_modulus / (3. * (1. - 2. * poisson_ratio))

Examples

HenckyOps.hpp, PlasticOps.hpp, dynamic_first_order_con_law.cpp, eigen_elastic.cpp, mofem/tutorials/vec-7/adjoint.cpp, and seepage.cpp.

Definition at line 95 of file dynamic_first_order_con_law.cpp.

help

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

static

[Check]

Definition at line 1188 of file dynamic_first_order_con_law.cpp.

lamme_lambda

double lamme_lambda

Initial value:

= young_modulus * poisson_ratio /
                      ((1. + poisson_ratio) * (1. - 2. * poisson_ratio))

Examples

dynamic_first_order_con_law.cpp.

Definition at line 98 of file dynamic_first_order_con_law.cpp.

mu

double mu = young_modulus / (2. * (1. + poisson_ratio))

Examples

ADOLCPlasticityMaterialModels.hpp, ElasticityMixedFormulation.hpp, NavierStokesElement.cpp, Remodeling.cpp, Remodeling.hpp, dynamic_first_order_con_law.cpp, plate.cpp, and shallow_wave.cpp.

Definition at line 97 of file dynamic_first_order_con_law.cpp.

omega

double omega = 1.

constexpr

Save field DOFS on vertices/tags.

Examples

dynamic_first_order_con_law.cpp, and shallow_wave.cpp.

Definition at line 92 of file dynamic_first_order_con_law.cpp.

poisson_ratio

double poisson_ratio = 0.

constexpr

Definition at line 94 of file dynamic_first_order_con_law.cpp.

shear_modulus_G

double shear_modulus_G = young_modulus / (2. * (1. + poisson_ratio))

Examples

HenckyOps.hpp, PlasticOps.hpp, dynamic_first_order_con_law.cpp, eigen_elastic.cpp, mofem/tutorials/vec-7/adjoint.cpp, and seepage.cpp.

Definition at line 96 of file dynamic_first_order_con_law.cpp.

SPACE_DIM

int SPACE_DIM

constexpr

Initial value:

=
    EXECUTABLE_DIMENSION

Examples

dynamic_first_order_con_law.cpp.

Definition at line 23 of file dynamic_first_order_con_law.cpp.

tsPrePostProc

boost::weak_ptr<TSPrePostProc> tsPrePostProc

static

Examples

dynamic_first_order_con_law.cpp, and free_surface.cpp.

Definition at line 404 of file dynamic_first_order_con_law.cpp.

young_modulus

double young_modulus = 1.

constexpr

Definition at line 93 of file dynamic_first_order_con_law.cpp.