Revision e75c84b4d4c8a26c254c972e0ad3b35075980133 authored by Thomas Fischer on 18 June 2021, 06:22:59 UTC, committed by Thomas Fischer on 18 June 2021, 15:56:34 UTC
1 parent 3908c72
TestFeHEX20.h
/**
* \copyright
* Copyright (c) 2012-2021, OpenGeoSys Community (http://www.opengeosys.org)
* Distributed under a Modified BSD License.
* See accompanying file LICENSE.txt or
* http://www.opengeosys.org/project/license
*
* \file
* Created on July 13, 2017, 2:40 PM
*/
#pragma once
#include "MeshLib/Elements/Hex.h"
#include "NumLib/Fem/FiniteElement/C0IsoparametricElements.h"
namespace FeTestData
{
class TestFeHEX20
{
public:
using ShapeFunction = NumLib::ShapeHex20;
// Fe type information
template <template <typename> class ShapeMatrixPolicy_>
using FeType = NumLib::FeHEX20<ShapeMatrixPolicy_>;
using MeshElementType = MeshLib::Hex20;
static const unsigned dim = 3; // MeshElementType::dimension;
static const unsigned e_nnodes = MeshElementType::n_all_nodes;
static const unsigned n_sample_pt_order2 = 2 * 2 * 2;
static const unsigned n_sample_pt_order3 = 3 * 3 * 3;
static const unsigned global_dim = MeshElementType::dimension;
const double perturbation = 0.2;
/// create a 20 node hexahedral element with curved edges
MeshElementType* createMeshElement()
{
auto** nodes = new MeshLib::Node*[e_nnodes];
// bottom
nodes[0] = new MeshLib::Node(0.0, 0.0, 0.0);
nodes[1] = new MeshLib::Node(1.0, 0.0, 0.0);
nodes[2] = new MeshLib::Node(1.0, 1.0, 0.0);
nodes[3] = new MeshLib::Node(0.0, 1.0, 0.0);
// top
nodes[4] = new MeshLib::Node(0.0, 0.0, 1.0);
nodes[5] = new MeshLib::Node(1.0, 0.0, 1.0);
nodes[6] = new MeshLib::Node(1.0, 1.0, 1.0);
nodes[7] = new MeshLib::Node(0.0, 1.0, 1.0);
// mid-edge nodes on bottom
nodes[8] = new MeshLib::Node(0.5, perturbation, 0.0);
nodes[9] = new MeshLib::Node(1.0 + perturbation, 0.5, 0.0);
nodes[10] = new MeshLib::Node(0.5, 1.0 + perturbation, 0.0);
nodes[11] = new MeshLib::Node(perturbation, 0.5, 0.0);
// mid-edge nodes on top
nodes[12] = new MeshLib::Node(0.5, perturbation, 1.0);
nodes[13] = new MeshLib::Node(1.0 + perturbation, 0.5, 1.0);
nodes[14] = new MeshLib::Node(0.5, 1.0 + perturbation, 1.0);
nodes[15] = new MeshLib::Node(perturbation, 0.5, 1.0);
// mid-edge nodes between the bottom and the top
nodes[16] = new MeshLib::Node(perturbation, 0.0, 0.5);
nodes[17] = new MeshLib::Node(1.0 + perturbation, 0.0, 0.5);
nodes[18] = new MeshLib::Node(1.0 + perturbation, 1.0, 0.5);
nodes[19] = new MeshLib::Node(perturbation, 1.0, 0.5);
return new MeshElementType(nodes);
}
double getVolume() const { return 1.0; }
};
} // namespace FeTestData
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...