Homec4science
Diffusion µSpectre prototype implementation (add_eigenstrain) Browse

muSpectre_proto/src/materials/material.hh
447e9272747badd_eigenstrain

Last Change
Raw File

material.hh
View Options

/**
* file material.hh
*
* @author Till Junge <till.junge@epfl.ch>
*
* @date 01 May 2017
*
* @brief Base class for materials (constitutive models)
*
* @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 <vector>
#include <unsupported/Eigen/CXX11/Tensor>
#include "common/common.hh"
namespace muSpectre {
template<Dim_t spatial_dimension, Dim_t material_dimension>
class MaterialBase {
public:
using Index = Index_t<spatial_dimension>;
//! Gradients are second-order tensors per pixel
using SecondArray = Eigen::Tensor<Real, spatial_dimension + 2>;
//! Tangent moduli are fourth-order tensors per pixel
using FourthArray = Eigen::Tensor<Real, spatial_dimension + 4>;
//! Default constructor
MaterialBase();
//! take responsibility for a pixel
void add_pixel(const Index && local_index);
virtual void compute_First_Piola_Kirchhoff_stress(const SecondArray & F,
SecondArray & P) = 0;
virtual void compute_tangent_moduli(const SecondArray & F,
FourthArray & K) = 0;
protected:
//! container of pixel indices using this material
std::vector<Index> indices;
private:
};
} // muSpectre
c4science · Help