neohookean_auxiliary_log_stress_atom.cpp File Reference

Verify the five-coordinate Neo-Hookean material map. More...

#include <MoFEM.hpp>
#include <NeoHookeanLogarithmicMaterial.hpp>

Go to the source code of this file.

Typedefs
using	Material = EshelbianPlasticity::NeoHookeanLogarithmicMaterial

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

Variables
static char	help []

Detailed Description

Verify the five-coordinate Neo-Hookean material map.

Definition in file neohookean_auxiliary_log_stress_atom.cpp.

Typedef Documentation

Material

using Material = EshelbianPlasticity::NeoHookeanLogarithmicMaterial

Definition at line 9 of file neohookean_auxiliary_log_stress_atom.cpp.

Function Documentation

main()

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

Definition at line 499 of file neohookean_auxiliary_log_stress_atom.cpp.

                                 {
  MoFEM::Core::Initialize(&argc, &argv, nullptr, help);
  try {
    moab::Core moab;
    MoFEM::Core core(moab);
    MoFEM::Interface &m_field = core;
    auto *json_config = m_field.getInterface<JsonConfigManager>();
    auto *meshsets_manager = m_field.getInterface<MeshsetsManager>();
    // Only material parameters are needed; this atom has no element geometry.
    CHKERR meshsets_manager->addMeshset(BLOCKSET, 1001, "MAT_NEOHOOKEAN");
    CHKERR meshsets_manager->setMeshsetFromFile();
    const auto block =
        json_config->getParamsFromBlockset("MAT_NEOHOOKEAN", 1001);
    if (block.size() != 2 || !block.count("c10") || !block.count("k"))
      SETERRQ(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
              "Atom requires MAT_NEOHOOKEAN JSON block 1001 with c10,k");
    const Material::Parameters parameters{block.at("c10"), block.at("k")};
    CHKERR Material::validateParameters(parameters);
    CHKERR checkNear(parameters.c10, 1.7, 1.e-14, "reference C10");
    CHKERR checkNear(parameters.bulkModulus, 8.5, 1.e-14, "reference K");
    PetscReal fd_step = 0., tolerance = 0.;
    CHKERR PetscOptionsGetReal(nullptr, nullptr, "-auxiliary_atom_fd_step",
                               &fd_step, nullptr);
    CHKERR PetscOptionsGetReal(nullptr, nullptr, "-auxiliary_atom_tolerance",
                               &tolerance, nullptr);
    if (!std::isfinite(fd_step) || fd_step <= 0. || !std::isfinite(tolerance) ||
        tolerance <= 0.)
      SETERRQ(PETSC_COMM_SELF, MOFEM_ATOM_TEST_INVALID,
              "Atom requires positive finite JSON finite-difference settings");
    CHKERR checkReference(parameters, tolerance);
    CHKERR checkReference({parameters.c10, 0.1}, tolerance);
    CHKERR checkSpectralStates(parameters, fd_step, tolerance);
    CHKERR checkForwardReferences(parameters, fd_step, tolerance);
    CHKERR checkLargeStress(parameters, tolerance);
    CHKERR checkSmallDeviator(parameters, tolerance);
    CHKERR checkEnergyScaling(tolerance);
    CHKERR checkConditioning(parameters, tolerance);
    CHKERR checkInvalidInputs(parameters);
    CHKERR PetscPrintf(
        PETSC_COMM_WORLD,
        "Neo-Hookean auxiliary logarithmic stress atom passed\n");
  }
  CATCH_ERRORS;
  MoFEM::Core::Finalize();
  return 0;
}

Variable Documentation

help

char help[]

static

Initial value:

=
    "Neo-Hookean auxiliary logarithmic stress material atom.\n"
    "Use -json_config neohookean_auxiliary_log_stress_atom.json.\n"

Definition at line 495 of file neohookean_auxiliary_log_stress_atom.cpp.