Homec4science
Diffusion µSpectre (18_OSX_fixes) Browse

muSpectre/tests/test_material_linear_elastic_generic.cc
3046b397fb8618_OSX_fixes

Last Change
Raw File

test_material_linear_elastic_generic.cc
View Options

/**
* @file test_material_linear_elastic_generic.cc
*
* @author Till Junge <till.junge@altermail.ch>
*
* @date 21 Sep 2018
*
* @brief test for the generic linear elastic law
*
* Copyright © 2018 Till Junge
*
* µSpectre 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, 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 Lesser General Public License
* along with µSpectre; see the file COPYING. If not, write to the
* Free Software Foundation, Inc., 59 Temple Place - Suite 330,
* * Boston, MA 02111-1307, USA.
*
* Additional permission under GNU GPL version 3 section 7
*
* If you modify this Program, or any covered work, by linking or combining it
* with proprietary FFT implementations or numerical libraries, containing parts
* covered by the terms of those libraries' licenses, the licensors of this
* Program grant you additional permission to convey the resulting work.
*/
#include "tests.hh"
#include "materials/material_linear_elastic_generic.cc"
#include "materials/materials_toolbox.hh"
#include <boost/mpl/list.hpp>
namespace muSpectre {
BOOST_AUTO_TEST_SUITE(material_linear_elastic_generic);
template <Dim_t Dim>
struct MatFixture
{
using Mat_t = MaterialLinearElasticGeneric<Dim, Dim>;
using T2_t = Eigen::Matrix<Real, Dim, Dim>;
using T4_t = T4Mat<Real, Dim>;
using V_t = Eigen::Matrix<Real, vsize(Dim), vsize(Dim)>;
constexpr static Real lambda{2}, mu{1.5};
constexpr static Real get_lambda() {return lambda;}
constexpr static Real get_mu() {return mu;}
constexpr static Real young{mu*(3*lambda + 2*mu)/(lambda + mu)};
constexpr static Real poisson{lambda/(2*(lambda + mu))};
using Hooke = MatTB::Hooke<Dim, T2_t, T4_t>;
MatFixture():
C_voigt{get_C_voigt()},
mat("material", this->C_voigt)
{}
static V_t get_C_voigt() {
V_t C{};
C.setZero();
C.template topLeftCorner<Dim, Dim>().setConstant(get_lambda());
C.template topLeftCorner<Dim, Dim>() += 2*get_mu()*T2_t::Identity();
constexpr Dim_t Rest{vsize(Dim)-Dim};
using Rest_t = Eigen::Matrix<Real, Rest, Rest>;
C.template bottomRightCorner<Rest, Rest>() += get_mu()*Rest_t::Identity();
return C;
}
V_t C_voigt;
Mat_t mat;
};
using mats = boost::mpl::list <MatFixture< twoD>,
MatFixture<threeD>>;
BOOST_FIXTURE_TEST_CASE_TEMPLATE(C_test, Fix, mats, Fix) {
const auto ref_C{Fix::Hooke::compute_C_T4(Fix::get_lambda(), Fix::get_mu())};
Real error{(ref_C-Fix::mat.get_C()).norm()};
BOOST_CHECK_LT(error, tol);
if (not (error < tol)) {
std::cout << "ref:" << std::endl << ref_C << std::endl;
std::cout << "new:" << std::endl << Fix::mat.get_C() << std::endl;
}
}
BOOST_AUTO_TEST_SUITE_END();
} // muSpectre
c4science · Help