arc_length_interface.cpp File Reference

#include <MoFEM.hpp>
#include <algorithm>
#include <limits>
#include <ArcLengthTools.hpp>
#include <CohesiveInterfaceElement.hpp>
#include <HookeOps.hpp>
#include <InterfaceGapArcLengthControl.hpp>
#include <CohesiveNaturalBoundaryBC.hpp>

Go to the source code of this file.

Classes
struct	BoundaryBCs
	Boundary conditions marker. More...
struct	CohesiveElement::AssembleRhsVectors
struct	ArcLengthInterfaceExample

Namespaces
namespace	CohesiveElement

Typedefs
using	DomainEle = PipelineManager::ElementsAndOpsByDim< SPACE_DIM >::DomainEle
using	DomainEleOp = DomainEle::UserDataOperator
using	BoundaryEle = PipelineManager::ElementsAndOpsByDim< SPACE_DIM >::BoundaryEle
using	PostProcEle = PostProcBrokenMeshInMoab< VolumeElementForcesAndSourcesCore >

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

Variables
constexpr int	SPACE_DIM = 3
static char	help []

Typedef Documentation

BoundaryEle

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

Definition at line 19 of file arc_length_interface.cpp.

DomainEle

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

Definition at line 17 of file arc_length_interface.cpp.

DomainEleOp

using DomainEleOp = DomainEle::UserDataOperator

Definition at line 18 of file arc_length_interface.cpp.

PostProcEle

using PostProcEle = PostProcBrokenMeshInMoab<VolumeElementForcesAndSourcesCore>

Definition at line 21 of file arc_length_interface.cpp.

Function Documentation

main()

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

[Check]

Definition at line 1087 of file arc_length_interface.cpp.

                                 {
 
  const string default_options = "-ksp_type fgmres \n"
                                 "-pc_type lu \n"
                                 "-pc_factor_mat_solver_type mumps\n"
                                 "-mat_mumps_icntl_20 0\n"
                                 "-ksp_monitor \n"
                                 "-ksp_atol 1e-10 \n"
                                 "-ksp_rtol 1e-10 \n"
                                 "-snes_monitor \n"
                                 "-snes_type newtonls \n"
                                 "-snes_linesearch_type l2 \n"
                                 "-snes_linesearch_monitor \n"
                                 "-snes_max_it 16 \n"
                                 "-snes_atol 1e-8 \n"
                                 "-snes_rtol 1e-8 \n"
                                 "-snes_converged_reason \n";
 
  string param_file = "param_file.petsc";
  if (!static_cast<bool>(ifstream(param_file))) {
    std::ofstream file(param_file.c_str(), std::ios::ate);
    if (file.is_open()) {
      file << default_options;
      file.close();
    }
  }
  MoFEM::Core::Initialize(&argc, &argv, param_file.c_str(), help);
 
  auto core_log = logging::core::get();
  core_log->add_sink(
      LogManager::createSink(LogManager::getStrmWorld(), "ARC_LENGTH"));
  LogManager::setLog("ARC_LENGTH");
  MOFEM_LOG_TAG("ARC_LENGTH", "ArcLength");
 
  try {
    moab::Core mb_instance;
    moab::Interface &moab = mb_instance;
 
    ArcLengthInterfaceExample example(moab);
    CHKERR example.runProblem();
  }
  CATCH_ERRORS;
 
  MoFEM::Core::Finalize();
 
  return 0;
}

Variable Documentation

help

char help[]

static

Initial value:

= "\
 -my_file mesh file name\n\
 -my_sr reduction of step size\n\
 -my_its_d desired number of steps\n\
 -my_ms maximal number of steps\n\
 -gamma arc-length step adaptation exponent\n\
 -min_arc_length_step minimum arc-length step\n\
 -max_arc_length_step maximum arc-length step\n\
 -field_eval_coords x,y,z coordinates where displacement and stress are evaluated\n\n"

Definition at line 33 of file arc_length_interface.cpp.

SPACE_DIM

constexpr int SPACE_DIM = 3

constexpr

Definition at line 16 of file arc_length_interface.cpp.