ElasticTieMeshTutorial Struct Reference

Inheritance diagram for ElasticTieMeshTutorial:

Collaboration diagram for ElasticTieMeshTutorial:

Public Member Functions
MoFEMErrorCode	checkResults () override
	[Postprocess results]
Public Member Functions inherited from ElasticTieMeshExample
	ElasticExample (MoFEM::Interface &m_field)
Public Member Functions inherited from ElasticExample
	ElasticExample (MoFEM::Interface &m_field)
virtual	~ElasticExample ()=default
MoFEMErrorCode	runProblem ()
	[Run problem]

Private Member Functions
MoFEMErrorCode	checkStressError (const int test)

Additional Inherited Members
Protected Types inherited from ElasticExample
using	Pipeline = boost::ptr_deque< ForcesAndSourcesCore::UserDataOperator >
Protected Member Functions inherited from ElasticTieMeshExample
MoFEMErrorCode	readMesh () override
	[Run problem]
MoFEMErrorCode	setupProblem () override
	[Read mesh]
MoFEMErrorCode	setUpSolver (SmartPetscObj< KSP > solver) override
	[Push operators to pipeline]
MoFEMErrorCode	kspSetUpAndSolve (SmartPetscObj< KSP > solver) override
Protected Member Functions inherited from ElasticExample
virtual MoFEMErrorCode	boundaryCondition ()
	[Set up problem]
virtual MoFEMErrorCode	assembleSystem ()
virtual MoFEMErrorCode	solveSystem ()
	[Solve]
virtual MoFEMErrorCode	outputResults ()
	[Solve]
virtual MoFEMErrorCode	addDomainMaterialOps (Pipeline &lhs_pipeline, Pipeline &rhs_pipeline)
	Add constitutive operators to the domain assembly pipelines.
Protected Attributes inherited from ElasticTieMeshExample
boost::shared_ptr< ElasticTie::CommonData >	tieData = nullptr
Range	physicalBoundaryFaces
Protected Attributes inherited from ElasticExample
MoFEM::Interface &	mField
boost::shared_ptr< MatrixDouble >	vectorFieldPtr = nullptr

Detailed Description

Examples

mofem/tutorials/vec-11_elastic_tie_mesh/elastic_tie_mesh.cpp.

Definition at line 39 of file elastic_tie_mesh.cpp.

Member Function Documentation

checkResults()

MoFEMErrorCode ElasticTieMeshTutorial::checkResults ( )

overridevirtual

[Postprocess results]

[Check]

Reimplemented from ElasticTieMeshExample.

Examples

mofem/tutorials/vec-11_elastic_tie_mesh/elastic_tie_mesh.cpp.

Definition at line 135 of file elastic_tie_mesh.cpp.

                                                    {
  MoFEMFunctionBegin;
 
  int test = 0;
  CHKERR PetscOptionsGetInt(PETSC_NULLPTR, "", "-test", &test, PETSC_NULLPTR);
  if (test == 10 || test == 11)
    CHKERR checkStressError(test);
  MoFEMFunctionReturn(0);
}

checkStressError()

MoFEMErrorCode ElasticTieMeshTutorial::checkStressError ( const int  test )

private

Examples

mofem/tutorials/vec-11_elastic_tie_mesh/elastic_tie_mesh.cpp.

Definition at line 49 of file elastic_tie_mesh.cpp.

                                                                      {
  MoFEMFunctionBegin;
 
  auto simple = mField.getInterface<Simple>();
 
  auto norm_fe = boost::make_shared<DomainEle>(mField);
  norm_fe->getRuleHook = [](int, int, int approx_order) {
    return 2 * approx_order + 1;
  };
 
  CHKERR AddHOOps<SPACE_DIM, SPACE_DIM, SPACE_DIM>::add(
      norm_fe->getOpPtrVector(), {H1}, "GEOMETRY");
 
  auto common_ptr = HookeOps::commonDataFactory<SPACE_DIM, GAUSS, DomainEle>(
      mField, norm_fe->getOpPtrVector(), "U", "MAT_ELASTIC", Sev::verbose);
 
  auto exact_stress_ptr = boost::make_shared<MatrixDouble>();
 
  constexpr double uniaxial_stress_xx = 1e3 * 0.3 / 3.0;
  auto uniaxial_stress = [](const double, const double, const double) {
    MatrixDouble stress((SPACE_DIM * (SPACE_DIM + 1)) / 2, 1);
    stress.clear();
    auto t_stress = getFTensor2SymmetricFromMat<
        SPACE_DIM, -1, DataLayoutTraits<DataLayout::CoeffsByGauss>>(stress);
    t_stress(0, 0) = uniaxial_stress_xx;
    return stress;
  };
 
  const AnalyticalSolutions::HollowCylinderUnderRadialPressure<SPACE_DIM>
      lame_solution(0.5, 1.0, 2.0, 1.0, 1e3, 0.3);
  auto lame_stress = [&lame_solution](const double x, const double y,
                                      const double z) {
    return lame_solution.stress(x, y, z);
  };
 
  MatrixFunc stress_func = (test == 11) ? MatrixFunc(lame_stress)
                                        : MatrixFunc(uniaxial_stress);
  norm_fe->getOpPtrVector().push_back(
      new OpGetTensor2SymmetricfromFunc<SPACE_DIM, SPACE_DIM>(
          exact_stress_ptr, stress_func));
 
  enum Norms { STRESS_ERROR_L2 = 0, STRESS_EXACT_L2, LAST_NORM };
  auto norms_vec =
      createVectorMPI(mField.get_comm(),
                      (mField.get_comm_rank() == 0) ? LAST_NORM : 0,
                      LAST_NORM);
  CHKERR VecZeroEntries(norms_vec);
 
  norm_fe->getOpPtrVector().push_back(
      new OpCalcNormL2Tensor2Symmetric<SPACE_DIM>(
          exact_stress_ptr, norms_vec, STRESS_EXACT_L2));
  norm_fe->getOpPtrVector().push_back(
      new OpCalcNormL2Tensor2Symmetric<SPACE_DIM>(
          common_ptr->getMatCauchyStress(), norms_vec, STRESS_ERROR_L2,
          exact_stress_ptr));
 
  CHKERR DMoFEMLoopFiniteElements(simple->getDM(), simple->getDomainFEName(),
                                  norm_fe);
  CHKERR VecAssemblyBegin(norms_vec);
  CHKERR VecAssemblyEnd(norms_vec);
 
  if (mField.get_comm_rank() == 0) {
    const double *norms;
    CHKERR VecGetArrayRead(norms_vec, &norms);
 
    const double error_l2 = std::sqrt(norms[STRESS_ERROR_L2]);
    const double exact_l2 = std::sqrt(norms[STRESS_EXACT_L2]);
 
    MOFEM_LOG_C("WORLD", Sev::inform, "STRESS_ERROR_L2 = %.16e\n", error_l2);
    MOFEM_LOG_C("WORLD", Sev::inform, "STRESS_EXACT_L2 = %.16e\n", exact_l2);
 
    CHKERR VecRestoreArrayRead(norms_vec, &norms);
 
    double max_error_l2 = 1e-2;
    CHKERR PetscOptionsGetReal(PETSC_NULLPTR, "", "-stress_error_l2_tol",
                               &max_error_l2, PETSC_NULLPTR);
    if (error_l2 >= max_error_l2) {
      SETERRQ(PETSC_COMM_SELF, MOFEM_DATA_INCONSISTENCY,
              "Stress L2 error %.16e exceeds tolerance %.16e", error_l2,
              max_error_l2);
    }
  }
 
  MoFEMFunctionReturn(0);
}

---

The documentation for this struct was generated from the following file:

• tutorials/vec-11_elastic_tie_mesh/elastic_tie_mesh.cpp