ElasticAdaptiveExample.hpp File Reference

elastic example with adaptive order refinement More...

#include <ElasticExample.hpp>
#include <AdaptiveOrderRef.hpp>
#include <HighOrderGeometryRef.hpp>

Go to the source code of this file.

Classes
struct	ElasticAdaptiveExample

Enumerations
enum	IndicInds { ERROR_TOT = 0 , COUNT_ENTS , ERROR_NORM , LAST_IND }

Variables
auto	sqr = [](double v) { return v * v; }
MoFEM::VectorFunc	analyticalDisplacement

Detailed Description

elastic example with adaptive order refinement

Copyright

Copyright (c) 2026

Definition in file ElasticAdaptiveExample.hpp.

Enumeration Type Documentation

IndicInds

enum IndicInds

Enumerator
ERROR_TOT
COUNT_ENTS
ERROR_NORM
LAST_IND

Definition at line 40 of file ElasticAdaptiveExample.hpp.

{ ERROR_TOT = 0, COUNT_ENTS, ERROR_NORM, LAST_IND };

Variable Documentation

analyticalDisplacement

MoFEM::VectorFunc analyticalDisplacement

inline

Definition at line 207 of file ElasticAdaptiveExample.hpp.

                                                               {
  double E = 1e3;
  double nu = 0.3;
 
  double r = std::sqrt(x * x + y * y);
 
  double a = 0.5;
  double b = 1.0;
 
  double pa = 2;
  double pb = 1;
 
  double A = (pa * sqr(a) - pb * sqr(b)) / (sqr(b) - sqr(a));
  double B = sqr(a) * sqr(b) * (pa - pb) / (sqr(b) - sqr(a));
 
  VectorDouble u;
  u.resize(SPACE_DIM);
  std::vector<double> v_u;
 
  double u_r = (1. + nu) / E * ((1. - 2. * nu) * A * r + B / r);
 
  if (SPACE_DIM == 2)
    v_u = {u_r / r * x, u_r / r * y};
  else
    v_u = {u_r / r * x, u_r / r * y, 0.0};
 
  for (int i = 0; i < SPACE_DIM; ++i)
    u(i) = v_u[i];
 
  return u;
};

sqr

auto sqr = [](double v) { return v * v; }

inline

Examples

mofem/tutorials/mix-0_mixed_poisson/mixed_poisson.cpp.

Definition at line 205 of file ElasticAdaptiveExample.hpp.

{ return v * v; };