Homec4science
Diffusion µSpectre prototype implementation (master) Browse

muSpectre_proto/tests/test_material_hyperelastic.cc
20b7cd3d088cmaster

Last Change
Raw File

test_material_hyperelastic.cc
View Options

/**
* file test_material_hyperelastic.cc
*
* @author Till Junge <till.junge@epfl.ch>
*
* @date 01 May 2017
*
* @brief Test the hyper-elastic material used in Section 4 of
* de Geus et al. / Comput. Methods in Appl. Mech. Engrg.
* 318, (2017), 412–430 https://doi.org/10.1016/j.cma.2016.12.032
*
* @section LICENCE
*
* Copyright (C) 2017 Till Junge
*
* µSpectre is free software; you can redistribute it and/or
* modify it under the terms of the GNU General Public License as
* published by the Free Software Foundation, either version 3, or (at
* your option) any later version.
*
* µSpectre 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
* General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with GNU Emacs; see the file COPYING. If not, write to the
* Free Software Foundation, Inc., 59 Temple Place - Suite 330,
* Boston, MA 02111-1307, USA.
*/
#include "materials/material_hyper_elastic.hh"
#include <boost/test/unit_test.hpp>
#include <boost/mpl/list.hpp>
namespace muSpectre {
const Real tol = 1e-14*100; //it's in percent
const Real E = 10.;
const Real nu = .3;
template<Dim_t s_dim, Dim_t m_dim>
struct mat_fixture: public MaterialHyperElastic<s_dim, m_dim>{
mat_fixture():
MaterialHyperElastic<s_dim, m_dim>(E, nu){
};
const static Dim_t DimS = s_dim;
const static Dim_t DimM = m_dim;
};
typedef boost::mpl::list<mat_fixture<2, 2>,
mat_fixture<2, 3>,
mat_fixture<3, 3>> test_types;
BOOST_FIXTURE_TEST_CASE_TEMPLATE(material_hyper_elastic_test, F, test_types, F) {
BOOST_CHECK_CLOSE(F::lambda, E*nu/((1+nu)*(1-2*nu)), tol);
BOOST_CHECK_CLOSE(F::mu, E/2/(1+nu), tol);
}
BOOST_FIXTURE_TEST_CASE_TEMPLATE(test_elastic_tensor, F, test_types, F) {
Eigen::Matrix<Real, vsize(F::DimM), vsize(F::DimM)> Cv;
F::VoigtConv::fourth_to_voigt(F::C, Cv);
auto & E = F::Young;
auto & nu = F::Poisson;
auto prefactor = E/(1+nu)/(1-2*nu);
const auto diff{vsize(F::DimM) - F::DimM};
for (Dim_t i = 0; i < F::DimM; ++i) {
BOOST_CHECK_CLOSE(Cv(i, i ), prefactor*(1-nu), 1e-12);
for (Dim_t j = i + 1; j < F::DimM; ++j ) {
BOOST_CHECK_CLOSE(Cv(i, j), prefactor*nu, 1e-12);
BOOST_CHECK_CLOSE(Cv(j, i), prefactor*nu, 1e-12);
}
}
for (Dim_t i = 0; i < diff; ++i) {
BOOST_CHECK_CLOSE(Cv(i+F::DimM, i+F::DimM), prefactor*(1-2*nu)/2, 1e-12);
for (Dim_t j = i + 1; j < diff; ++j ) {
BOOST_CHECK_CLOSE(Cv(i+F::DimM, j+F::DimM), 0, 1e-12);
}
for (Dim_t j = 0; j < diff; ++j) {
BOOST_CHECK_CLOSE(Cv(i+F::DimM, j ), 0, 1e-12);
BOOST_CHECK_CLOSE(Cv(i, j+F::DimM), 0, 1e-12);
}
}
}
} // muSpectre
c4science · Help