Plate Struct Reference

Collaboration diagram for Plate:

Public Member Functions
	Plate (MoFEM::Interface &m_field)
MoFEMErrorCode	runProblem ()
	[Output results] More...

Private Member Functions
MoFEMErrorCode	readMesh ()
	[Read mesh] More...
MoFEMErrorCode	setupProblem ()
	[Read mesh] More...
MoFEMErrorCode	boundaryCondition ()
	[Set up problem] More...
MoFEMErrorCode	assembleSystem ()
	[Boundary condition] More...
MoFEMErrorCode	solveSystem ()
	[Push operators to pipeline] More...
MoFEMErrorCode	outputResults ()
	[Solve system] More...

Private Attributes
MoFEM::Interface &	mField

Detailed Description

Examples

plate.cpp.

Definition at line 140 of file plate.cpp.

Constructor & Destructor Documentation

Plate()

Plate::Plate ( MoFEM::Interface &  m_field )

inline

Definition at line 142 of file plate.cpp.

: mField(m_field) {}

Member Function Documentation

assembleSystem()

MoFEMErrorCode Plate::assembleSystem ( )

private

[Boundary condition]

[Push operators to pipeline]

calculate jacobian

Examples

plate.cpp.

Definition at line 234 of file plate.cpp.

                                     {
  MoFEMFunctionBegin;
 
  auto pipeline_mng = mField.getInterface<PipelineManager>();
 
  auto rule_lhs = [](int, int, int p) -> int { return 2 * p; };
  auto rule_rhs = [](int, int, int p) -> int { return 2 * p; };
  auto rule_2 = [this](int, int, int) { return 2 * order; };
 
  CHKERR pipeline_mng->setDomainLhsIntegrationRule(rule_lhs);
  CHKERR pipeline_mng->setDomainRhsIntegrationRule(rule_rhs);
 
  CHKERR pipeline_mng->setSkeletonLhsIntegrationRule(rule_2);
  CHKERR pipeline_mng->setSkeletonRhsIntegrationRule(rule_2);
  CHKERR pipeline_mng->setBoundaryLhsIntegrationRule(rule_2);
  CHKERR pipeline_mng->setBoundaryRhsIntegrationRule(rule_2);
 
  auto det_ptr = boost::make_shared<VectorDouble>();
  auto jac_ptr = boost::make_shared<MatrixDouble>();
  auto inv_jac_ptr = boost::make_shared<MatrixDouble>();
  auto base_mass_ptr = boost::make_shared<MatrixDouble>();
  auto data_l2_ptr = boost::make_shared<EntitiesFieldData>(MBENTITYSET);
 
  auto mat_D_ptr = plate_stiffness();
 
  auto push_ho_direcatives = [&](auto &pipeline) {
    pipeline.push_back(new OpBaseDerivativesMass<BASE_DIM>(
        base_mass_ptr, data_l2_ptr, AINSWORTH_LEGENDRE_BASE, L2));
    pipeline.push_back(new OpBaseDerivativesNext<BASE_DIM>(
        BaseDerivatives::SecondDerivative, base_mass_ptr, data_l2_ptr,
        AINSWORTH_LEGENDRE_BASE, H1));
  };
 
  /**
   * @brief calculate jacobian
   *
   */
  auto push_jacobian = [&](auto &pipeline) {
    pipeline.push_back(new OpSetHOWeightsOnFace());
    pipeline.push_back(new OpCalculateHOJac<SPACE_DIM>(jac_ptr));
    pipeline.push_back(
        new OpInvertMatrix<SPACE_DIM>(jac_ptr, det_ptr, inv_jac_ptr));
    // push first base derivatives tp physical element shape
    pipeline.push_back(new OpSetInvJacH1ForFace<1>(inv_jac_ptr));
    // push second base derivatives tp physical element shape
    pipeline.push_back(new OpSetInvJacH1ForFace<2>(inv_jac_ptr));
  };
 
  push_ho_direcatives(pipeline_mng->getOpDomainLhsPipeline());
  push_jacobian(pipeline_mng->getOpDomainLhsPipeline());
 
  pipeline_mng->getOpDomainLhsPipeline().push_back(
      new OpDomainPlateStiffness("U", "U", mat_D_ptr));
  // pipeline_mng->getOpDomainLhsPipeline().push_back(new OpDomainGradGrad(
  // "U", "U", [](const double, const double, const double) { return 1; }));
 
  pipeline_mng->getOpDomainRhsPipeline().push_back(
      new OpDomainPlateLoad("U", source));
 
  // Push operators to the Pipeline for Skeleton
  auto side_fe_ptr = boost::make_shared<FaceSideEle>(mField);
  push_ho_direcatives(side_fe_ptr->getOpPtrVector());
  push_jacobian(side_fe_ptr->getOpPtrVector());
  side_fe_ptr->getOpPtrVector().push_back(new OpCalculateSideData("U", "U"));
 
  // Push operators to the Pipeline for Skeleton
  pipeline_mng->getOpSkeletonLhsPipeline().push_back(
      new OpH1LhsSkeleton(side_fe_ptr, mat_D_ptr));
 
  MoFEMFunctionReturn(0);
}

boundaryCondition()

MoFEMErrorCode Plate::boundaryCondition ( )

private

[Set up problem]

[Boundary condition]

Examples

plate.cpp.

Definition at line 199 of file plate.cpp.

                                        {
  MoFEMFunctionBegin;
 
  // get edges and vertices on body skin
  auto get_skin = [&]() {
    Range body_ents;
    MOAB_THROW(
        mField.get_moab().get_entities_by_dimension(0, SPACE_DIM, body_ents));
    Skinner skin(&mField.get_moab());
    Range skin_ents;
    MOAB_THROW(skin.find_skin(0, body_ents, false, skin_ents));
    Range verts;
    MOAB_THROW(mField.get_moab().get_connectivity(skin_ents, verts, true));
    skin_ents.merge(verts);
    return skin_ents;
  };
 
  // remove dofs on skin edges from priblem
  auto remove_dofs_from_problem = [&](Range &&skin) {
    MoFEMFunctionBegin;
    auto problem_mng = mField.getInterface<ProblemsManager>();
    auto simple = mField.getInterface<Simple>();
    CHKERR problem_mng->removeDofsOnEntities(simple->getProblemName(), "U",
                                             skin, 0, 1);
    MoFEMFunctionReturn(0);
  };
 
  //  it make plate simply supported
  CHKERR remove_dofs_from_problem(get_skin());
 
  MoFEMFunctionReturn(0);
}

outputResults()

MoFEMErrorCode Plate::outputResults ( )

private

[Solve system]

[Output results]

Examples

plate.cpp.

Definition at line 335 of file plate.cpp.

                                    {
  MoFEMFunctionBegin;
 
  auto pipeline_mng = mField.getInterface<PipelineManager>();
  pipeline_mng->getDomainLhsFE().reset();
  pipeline_mng->getSkeletonRhsFE().reset();
  pipeline_mng->getSkeletonLhsFE().reset();
  pipeline_mng->getBoundaryRhsFE().reset();
  pipeline_mng->getBoundaryLhsFE().reset();
 
  auto post_proc_fe = boost::make_shared<PostProcEle>(mField);
 
  auto u_ptr = boost::make_shared<VectorDouble>();
  post_proc_fe->getOpPtrVector().push_back(
      new OpCalculateScalarFieldValues("U", u_ptr));
 
  using OpPPMap = OpPostProcMapInMoab<2, 2>;
 
  post_proc_fe->getOpPtrVector().push_back(
 
      new OpPPMap(post_proc_fe->getPostProcMesh(),
                  post_proc_fe->getMapGaussPts(),
 
                  {{"U", u_ptr}},
 
                  {}, {}, {}
 
                  )
 
  );
 
  pipeline_mng->getDomainRhsFE() = post_proc_fe;
  CHKERR pipeline_mng->loopFiniteElements();
  CHKERR post_proc_fe->writeFile("out_result.h5m");
 
  MoFEMFunctionReturn(0);
}

readMesh()

MoFEMErrorCode Plate::readMesh ( )

private

[Read mesh]

Examples

plate.cpp.

Definition at line 158 of file plate.cpp.

                               {
  MoFEMFunctionBegin;
 
  auto simple = mField.getInterface<Simple>();
  CHKERR simple->getOptions();
  CHKERR simple->loadFile();
 
  MoFEMFunctionReturn(0);
}

runProblem()

MoFEMErrorCode Plate::runProblem

(

)

[Output results]

[Run program]

Examples

plate.cpp.

Definition at line 375 of file plate.cpp.

                                 {
  MoFEMFunctionBegin;
 
  CHKERR readMesh();
  CHKERR setupProblem();
  CHKERR boundaryCondition();
  CHKERR assembleSystem();
  CHKERR solveSystem();
  // CHKERR checkResults();
  CHKERR outputResults();
 
  MoFEMFunctionReturn(0);
}

setupProblem()

MoFEMErrorCode Plate::setupProblem ( )

private

[Read mesh]

[Set up problem]

Examples

plate.cpp.

Definition at line 170 of file plate.cpp.

                                   {
  MoFEMFunctionBegin;
 
  auto simple = mField.getInterface<Simple>();
 
  CHKERR PetscOptionsGetInt(PETSC_NULL, "", "-order", &order, PETSC_NULL);
  CHKERR PetscOptionsGetScalar(PETSC_NULL, "", "-penalty", &penalty,
                               PETSC_NULL);
  CHKERR PetscOptionsGetScalar(PETSC_NULL, "", "-phi", &phi, PETSC_NULL);
  CHKERR PetscOptionsGetScalar(PETSC_NULL, "", "-nitsche", &nitsche,
                               PETSC_NULL);
 
  MOFEM_LOG("WORLD", Sev::inform) << "Set order: " << order;
  MOFEM_LOG("WORLD", Sev::inform) << "Set penalty: " << penalty;
  MOFEM_LOG("WORLD", Sev::inform) << "Set phi: " << phi;
  MOFEM_LOG("WORLD", Sev::inform) << "Set nitche: " << nitsche;
 
  CHKERR simple->addDomainField("U", H1, AINSWORTH_LEGENDRE_BASE, FIELD_DIM);
  CHKERR simple->addSkeletonField("U", H1, AINSWORTH_LEGENDRE_BASE, FIELD_DIM);
  CHKERR simple->addBoundaryField("U", H1, AINSWORTH_LEGENDRE_BASE, FIELD_DIM);
 
  CHKERR simple->setFieldOrder("U", order);
  CHKERR simple->setUp();
 
  MoFEMFunctionReturn(0);
}

solveSystem()

MoFEMErrorCode Plate::solveSystem ( )

private

[Push operators to pipeline]

[Solve system]

Examples

plate.cpp.

Definition at line 308 of file plate.cpp.

                                  {
  MoFEMFunctionBegin;
 
  auto pipeline_mng = mField.getInterface<PipelineManager>();
  auto simple = mField.getInterface<Simple>();
 
  auto ksp_solver = pipeline_mng->createKSP();
  CHKERR KSPSetFromOptions(ksp_solver);
  CHKERR KSPSetUp(ksp_solver);
 
  // Create RHS and solution vectors
  auto dm = simple->getDM();
  auto F = createDMVector(dm);
  auto D = vectorDuplicate(F);
 
  CHKERR KSPSolve(ksp_solver, F, D);
 
  // Scatter result data on the mesh
  CHKERR VecGhostUpdateBegin(D, INSERT_VALUES, SCATTER_FORWARD);
  CHKERR VecGhostUpdateEnd(D, INSERT_VALUES, SCATTER_FORWARD);
  CHKERR DMoFEMMeshToLocalVector(dm, D, INSERT_VALUES, SCATTER_REVERSE);
 
  MoFEMFunctionReturn(0);
}

Member Data Documentation

mField

MoFEM::Interface& Plate::mField

private

Definition at line 154 of file plate.cpp.

---

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

• tutorials/vec-6/plate.cpp