Homec4science
Diffusion µSpectre (feat/add_lam_pixels) Browse

muSpectre/src/cell/cell_laminate.hh
a2b173105256feat/add_lam_pixels

Last Change
Raw File

cell_laminate.hh
View Options

/**
* @file cell_split.hh
*
* @author Ali Falsafi <ali.falsafi@epfl.ch>
*
* @date 10 Sep 2018
*
* @brief Base class representing a unit cell able to handle
* split material assignments
*
* Copyright © 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.
*/
#ifndef SRC_CELL_CELL_LAMINATE_HH_
#define SRC_CELL_CELL_LAMINATE_HH_
#include "common/muSpectre_common.hh"
#include "libmugrid/ccoord_operations.hh"
#include "libmugrid/field.hh"
#include "libmugrid/cpp_compliance.hh"
#include "materials/material_base.hh"
#include "projection/projection_base.hh"
#include "cell/cell_traits.hh"
#include "cell/cell_base.hh"
#include "cell/cell_split.hh"
namespace muSpectre {
//! DimS spatial dimension (dimension of problem
//! DimM material_dimension (dimension of constitutive law)
/* This class handles the cells that has splitly assigned material to their
* pixels */
template <Dim_t DimS, Dim_t DimM = DimS>
class CellLaminate : public CellSplit<DimS, DimM> {
public:
using Parent = CellSplit<DimS, DimM>;
//! global field collection
using FieldCollection_t = muGrid::GlobalFieldCollection<DimS>;
using Projection_t = ProjectionBase<DimS, DimM>;
//! projections handled through `std::unique_ptr`s
using Projection_ptr = std::unique_ptr<Projection_t>;
using StrainField_t =
muGrid::TensorField<FieldCollection_t, Real, secondOrder, DimM>;
//! expected type for stress fields
using StressField_t =
muGrid::TensorField<FieldCollection_t, Real, secondOrder, DimM>;
//! expected type for tangent stiffness fields
using TangentField_t =
muGrid::TensorField<FieldCollection_t, Real, fourthOrder, DimM>;
//! combined stress and tangent field
using FullResponse_t =
std::tuple<const StressField_t &, const TangentField_t &>;
//! ref to constant vector
using ConstVector_ref = typename Parent::ConstVector_ref;
//! Default constructor
CellLaminate() = delete;
//! constructor using sizes and resolution
CellLaminate(Projection_ptr projection,
SplitCell is_cell_split = SplitCell::laminate);
//! Copy constructor
CellLaminate(const CellLaminate & other) = delete;
//! Move constructor
CellLaminate(CellLaminate && other) = default;
//! Destructor
virtual ~CellLaminate() = default;
//! Copy assignment operator
CellLaminate & operator=(const CellLaminate & other) = delete;
//! Move assignment operator
CellLaminate & operator=(CellLaminate && other) = default;
/* ---------------------------------------------------------------------- */
protected:
// full resppnse is composed of the stresses and tangent matrix
FullResponse_t evaluate_stress_tangent(StrainField_t & F) override;
// std::array<ConstVector_ref, 2> evaluate_stress_tangent() override;
// ConstVector_ref evaluate_stress() override;
};
} // namespace muSpectre
#endif // SRC_CELL_CELL_LAMINATE_HH_
c4science · Help