OGS
ComputeElementVolumeNumerically.cpp
Go to the documentation of this file.
1
12#include "ComputeElementVolumeNumerically.h"
13
14#include <typeinfo>
15
16#include "MeshLib/Elements/Element.h"
17#include "MeshLib/Elements/Hex.h"
18#include "MeshLib/Elements/Line.h"
19#include "MeshLib/Elements/Prism.h"
20#include "MeshLib/Elements/Pyramid.h"
21#include "MeshLib/Elements/Quad.h"
22#include "MeshLib/Elements/Tet.h"
23#include "MeshLib/Elements/Tri.h"
24#include "MeshLib/MeshEnums.h"
25#include "NumLib/Fem/InitShapeMatrices.h"
26#include "NumLib/Fem/Integration/IntegrationMethodRegistry.h"
27#include "NumLib/Fem/ShapeFunction/ShapeHex20.h"
28#include "NumLib/Fem/ShapeFunction/ShapeHex8.h"
29#include "NumLib/Fem/ShapeFunction/ShapeLine2.h"
30#include "NumLib/Fem/ShapeFunction/ShapeLine3.h"
31#include "NumLib/Fem/ShapeFunction/ShapePrism15.h"
32#include "NumLib/Fem/ShapeFunction/ShapePrism6.h"
33#include "NumLib/Fem/ShapeFunction/ShapePyra13.h"
34#include "NumLib/Fem/ShapeFunction/ShapePyra5.h"
35#include "NumLib/Fem/ShapeFunction/ShapeQuad4.h"
36#include "NumLib/Fem/ShapeFunction/ShapeQuad8.h"
37#include "NumLib/Fem/ShapeFunction/ShapeQuad9.h"
38#include "NumLib/Fem/ShapeFunction/ShapeTet10.h"
39#include "NumLib/Fem/ShapeFunction/ShapeTet4.h"
40#include "NumLib/Fem/ShapeFunction/ShapeTri3.h"
41#include "NumLib/Fem/ShapeFunction/ShapeTri6.h"
42#include "NumLib/Fem/ShapeMatrixPolicy.h"
43
44namespace MeshToolsLib
45{
46template <typename ShapeFunction>
47double computeElementVolumeNumerically(MeshLib::Element const& e)
48{
49 // Space dimension is set to 3 in case that 1D or 2D element is inclined.
50 constexpr int space_dim = 3;
51 using ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, space_dim>;
52
53 // Integration order is set to 3:
54 auto const& integration_method =
55 NumLib::IntegrationMethodRegistry::template getIntegrationMethod<
56 typename ShapeFunction::MeshElement>(NumLib::IntegrationOrder{3});
57
58 auto const shape_function_data =
59 NumLib::initShapeMatrices<ShapeFunction, ShapeMatricesType, space_dim>(
60 e, false /*is_axially_symmetric*/, integration_method);
61
62 auto const n_integration_points = integration_method.getNumberOfPoints();
63 double volume = 0.0;
64 for (unsigned ip = 0; ip < n_integration_points; ++ip)
65 {
66 auto const weight = integration_method.getWeightedPoint(ip).getWeight();
67 volume += shape_function_data[ip].detJ * weight;
68 }
69
70 return volume;
71}
72
73double computeElementVolumeNumerically(MeshLib::Element const& e)
74{
75 switch (e.getCellType())
76 {
77 case MeshLib::CellType::LINE2:
78 return computeElementVolumeNumerically<NumLib::ShapeLine2>(e);
79 case MeshLib::CellType::LINE3:
80 return computeElementVolumeNumerically<NumLib::ShapeLine3>(e);
81 case MeshLib::CellType::TRI3:
82 return computeElementVolumeNumerically<NumLib::ShapeTri3>(e);
83 case MeshLib::CellType::TRI6:
84 return computeElementVolumeNumerically<NumLib::ShapeTri6>(e);
85 case MeshLib::CellType::QUAD4:
86 return computeElementVolumeNumerically<NumLib::ShapeQuad4>(e);
87 case MeshLib::CellType::QUAD8:
88 return computeElementVolumeNumerically<NumLib::ShapeQuad8>(e);
89 case MeshLib::CellType::QUAD9:
90 return computeElementVolumeNumerically<NumLib::ShapeQuad9>(e);
91 case MeshLib::CellType::TET4:
92 return computeElementVolumeNumerically<NumLib::ShapeTet4>(e);
93 case MeshLib::CellType::HEX8:
94 return computeElementVolumeNumerically<NumLib::ShapeHex8>(e);
95 case MeshLib::CellType::HEX20:
96 return computeElementVolumeNumerically<NumLib::ShapeHex20>(e);
97 case MeshLib::CellType::TET10:
98 return computeElementVolumeNumerically<NumLib::ShapeTet10>(e);
99 case MeshLib::CellType::PRISM6:
100 return computeElementVolumeNumerically<NumLib::ShapePrism6>(e);
101 case MeshLib::CellType::PRISM15:
102 return computeElementVolumeNumerically<NumLib::ShapePrism15>(e);
103 case MeshLib::CellType::PYRAMID5:
104 return computeElementVolumeNumerically<NumLib::ShapePyra5>(e);
105 case MeshLib::CellType::PYRAMID13:
106 return computeElementVolumeNumerically<NumLib::ShapePyra13>(e);
107 default:
108 OGS_FATAL(
109 "Numerical volume calculation is not available for element "
110 "with type {}. ",
111 MeshLib::CellType2String(e.getCellType()));
112 }
113}
114
115} // namespace MeshToolsLib
Element.h
Definition of the Element class.
OGS_FATAL
#define OGS_FATAL(...)
Definition Error.h:26
Hex.h
Definition of the Hex class.
Line.h
Definition of the Line class.
MeshEnums.h
Definition of mesh-related Enumerations.
Prism.h
Definition of the Prism class.
Pyramid.h
Definition of the Pyramid class.
Quad.h
Definition of the Quad class.
Tet.h
Definition of the Tet class.
Tri.h
Definition of the Tri class.
MeshLib::Element::getCellType
virtual CellType getCellType() const =0
MeshLib::CellType2String
std::string CellType2String(const CellType t)
Given a MeshElemType this returns the appropriate string.
Definition MeshEnums.cpp:157
MeshToolsLib::computeElementVolumeNumerically
double computeElementVolumeNumerically(MeshLib::Element const &e)
Definition ComputeElementVolumeNumerically.cpp:47
NumLib::initShapeMatrices
std::vector< typename ShapeMatricesType::ShapeMatrices, Eigen::aligned_allocator< typename ShapeMatricesType::ShapeMatrices > > initShapeMatrices(MeshLib::Element const &e, bool const is_axially_symmetric, IntegrationMethod const &integration_method)
Definition InitShapeMatrices.h:53