Homec4science
Diffusion akantu (features/solver-rewrite) Browse

rAKA/test/test_model/patch_tests/patch_test_linear_solid_mechanics_fixture.hh
165c792569e1features/solver-rewrite

Last Change
Raw File

patch_test_linear_solid_mechanics_fixture.hh
View Options

/**
* @file patch_test_linear_solid_mechanics_fixture.hh
*
* @author Nicolas Richart <nicolas.richart@epfl.ch>
*
* @date creation: Tue Jan 30 2018
*
* @brief SolidMechanics patch tests fixture
*
* @section LICENSE
*
* Copyright (©) 2016-2018 EPFL (Ecole Polytechnique Fédérale de Lausanne)
* Laboratory (LSMS - Laboratoire de Simulation en Mécanique des Solides)
*
* Akantu is free software: you can redistribute it and/or modify it under the
* terms of the GNU Lesser General Public License as published by the Free
* Software Foundation, either version 3 of the License, or (at your option) any
* later version.
*
* Akantu is distributed in the hope that it will be useful, but WITHOUT ANY
* WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR
* A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more
* details.
*
* You should have received a copy of the GNU Lesser General Public License
* along with Akantu. If not, see <http://www.gnu.org/licenses/>.
*
*/
/* -------------------------------------------------------------------------- */
#include "patch_test_linear_fixture.hh"
#include "solid_mechanics_model.hh"
/* -------------------------------------------------------------------------- */
#ifndef __AKANTU_PATCH_TEST_LINEAR_SOLID_MECHANICS_FIXTURE_HH__
#define __AKANTU_PATCH_TEST_LINEAR_SOLID_MECHANICS_FIXTURE_HH__
/* -------------------------------------------------------------------------- */
template <typename tuple_>
class TestPatchTestSMMLinear
: public TestPatchTestLinear<std::tuple_element_t<0, tuple_>,
SolidMechanicsModel> {
using parent =
TestPatchTestLinear<std::tuple_element_t<0, tuple_>, SolidMechanicsModel>;
public:
static constexpr bool plane_strain = std::tuple_element_t<1, tuple_>::value;
void applyBC() override {
parent::applyBC();
auto & displacement = this->model->getDisplacement();
this->applyBConDOFs(displacement);
}
void checkAll() {
auto & displacement = this->model->getDisplacement();
auto & mat = this->model->getMaterial(0);
this->checkDOFs(displacement);
this->checkGradient(mat.getGradU(this->type), displacement);
this->checkResults(
[&](const Matrix<Real> & pstrain) {
Real nu = this->model->getMaterial(0).get("nu");
Real E = this->model->getMaterial(0).get("E");
auto strain = (pstrain + pstrain.transpose()) / 2.;
auto trace = strain.trace();
auto lambda = nu * E / ((1 + nu) * (1 - 2 * nu));
auto mu = E / (2 * (1 + nu));
if (not this->plane_strain) {
lambda = nu * E / (1 - nu * nu);
}
decltype(strain) stress(this->dim, this->dim);
if (this->dim == 1) {
stress(0, 0) = E * strain(0, 0);
} else {
for (UInt i = 0; i < this->dim; ++i)
for (UInt j = 0; j < this->dim; ++j)
stress(i, j) =
(i == j) * lambda * trace + 2 * mu * strain(i, j);
}
return stress;
},
mat.getStress(this->type), displacement);
}
};
template <typename tuple_>
constexpr bool TestPatchTestSMMLinear<tuple_>::plane_strain;
template <typename T> struct invalid_plan_stress : std::true_type {};
template <typename type, typename bool_c>
struct invalid_plan_stress<std::tuple<type, bool_c>>
: aka::bool_constant<ElementClass<type::value>::getSpatialDimension() !=
2 and
not bool_c::value> {};
using true_false =
std::tuple<aka::bool_constant<true>, aka::bool_constant<false>>;
template <typename T> using valid_types = aka::negation<invalid_plan_stress<T>>;
using types = gtest_list_t<
tuple_filter_t<valid_types, cross_product_t<TestElementTypes, true_false>>>;
TYPED_TEST_CASE(TestPatchTestSMMLinear, types);
#endif /* __AKANTU_PATCH_TEST_LINEAR_SOLID_MECHANICS_FIXTURE_HH__ */
c4science · Help