gsMaterialMatrixLinear.h

Go to the documentation of this file.

 
#pragma once
 
#include <gsKLShell/src/gsMaterialMatrixBaseDim.h>
#include <gsKLShell/src/gsMaterialMatrixUtils.h>
#include <gsIO/gsOptionList.h>
 
namespace gismo
{
 
 
template <  short_t dim,
            class T     >
class gsMaterialMatrixLinear : public gsMaterialMatrixBaseDim<dim,T>
{
public:
 
    typedef T Scalar_t;
 
    GISMO_OVERRIDE_CLONE_FUNCTION(gsMaterialMatrixLinear);
 
    using Base = gsMaterialMatrixBaseDim<dim,T>;
 
    typedef typename Base::function_ptr function_ptr;
 
    enum {Linear=1};
 
    gsMaterialMatrixLinear();
 
    gsMaterialMatrixLinear(   const gsFunctionSet<T> & mp,
                        const gsFunctionSet<T> & thickness);
 
    gsMaterialMatrixLinear(   const gsFunctionSet<T> & mp,
                        const gsFunctionSet<T> & thickness,
                        const gsFunctionSet<T> & YoungsModulus,
                        const gsFunctionSet<T> & PoissonRatio);
 
    gsMaterialMatrixLinear(   const gsFunctionSet<T> & mp,
                        const gsFunctionSet<T> & thickness,
                        const gsFunctionSet<T> & YoungsModulus,
                        const gsFunctionSet<T> & PoissonRatio,
                        const gsFunctionSet<T> & Density);
 
protected:
    gsMaterialMatrixLinear(   const gsFunctionSet<T> * mp,
                        const gsFunctionSet<T> * thickness,
                        const gsFunctionSet<T> & YoungsModulus,
                        const gsFunctionSet<T> & PoissonRatio,
                        const gsFunctionSet<T> * Density);
public:
    gsMaterialMatrixLinear(   const gsFunctionSet<T> & mp,
                        const gsFunctionSet<T> & thickness,
                        const std::vector<gsFunctionSet<T> *> &pars);
 
    gsMaterialMatrixLinear(   const gsFunctionSet<T> & mp,
                        const gsFunctionSet<T> & thickness,
                        const std::vector<gsFunctionSet<T> *> &pars,
                        const gsFunctionSet<T> & Density);
 
    gsMaterialMatrixLinear(    const gsFunctionSet<T> & thickness,
                        const std::vector<gsFunctionSet<T> *> &pars);
 
    gsMaterialMatrixLinear(   const gsFunctionSet<T> & thickness,
                        const std::vector<gsFunctionSet<T> *> &pars,
                        const gsFunctionSet<T> & Density);
 
protected:
    gsMaterialMatrixLinear(   const gsFunctionSet<T> * mp,
                        const gsFunctionSet<T> * thickness,
                        const std::vector<gsFunctionSet<T> *> &pars,
                        const gsFunctionSet<T> * Density);
public:
 
    inline enum MatIntegration isMatIntegrated() const override {return MatIntegration::Constant; }
 
    inline enum MatIntegration isVecIntegrated() const override {return MatIntegration::Constant; }
 
    void defaultOptions() override;
 
    gsMatrix<T> eval3D_matrix (const index_t patch, const gsMatrix<T> & u, const gsMatrix<T>& z, enum MaterialOutput out = MaterialOutput::Generic) const override;
    using Base::eval3D_matrix;
 
    gsMatrix<T> eval3D_dmatrix(const index_t patch, const gsMatrix<T> & u, const gsMatrix<T>& z, enum MaterialOutput out = MaterialOutput::Generic) const override;
    using Base::eval3D_dmatrix;
 
    gsMatrix<T> eval3D_vector (const index_t patch, const gsMatrix<T> & u, const gsMatrix<T>& z, enum MaterialOutput out = MaterialOutput::Generic) const override;
    using Base::eval3D_vector;
 
    gsMatrix<T> eval3D_CauchyVector (const index_t patch, const gsMatrix<T> & u, const gsMatrix<T>& z, enum MaterialOutput out = MaterialOutput::Generic) const override;
 
    gsMatrix<T> eval3D_matrix_C(const gsMatrix<T> & Cmat, const index_t patch, const gsVector<T>& u, const T z, enum MaterialOutput out = MaterialOutput::Generic) const override;
    gsMatrix<T> eval3D_vector_C(const gsMatrix<T> & Cmat, const index_t patch, const gsVector<T> & u, const T z, enum MaterialOutput out = MaterialOutput::Generic) const override;
 
    gsMatrix<T> eval3D_pstress(const index_t patch, const gsMatrix<T> & u, const gsMatrix<T>& z, enum MaterialOutput out = MaterialOutput::Generic) const override;
 
    gsMatrix<T> eval3D_pstressDir(const index_t patch, const gsMatrix<T> & u, const gsMatrix<T>& z, enum MaterialOutput out = MaterialOutput::Generic) const override;
 
    gsMatrix<T> eval3D_CauchyPStress(const index_t patch, const gsMatrix<T> & u, const gsMatrix<T>& z, enum MaterialOutput out = MaterialOutput::Generic) const override;
 
    gsMatrix<T> eval3D_stress(const index_t patch, const gsMatrix<T> & u, const gsMatrix<T> & z, enum MaterialOutput out) const override;
    using Base::eval3D_stress;
 
    virtual gsMatrix<T> eval3D_CauchyStress(const index_t patch, const gsMatrix<T> & u, const gsMatrix<T> & z, enum MaterialOutput out) const override;
    using Base::eval3D_CauchyStress;
 
    gsMatrix<T> eval3D_detF (const index_t patch, const gsMatrix<T> & u, const gsMatrix<T>& z, enum MaterialOutput out = MaterialOutput::Generic) const override;
 
    void setYoungsModulus(const gsFunctionSet<T> & YoungsModulus) override { Base::setParameter(0,YoungsModulus); }
 
    const function_ptr getYoungsModulus() const override { return Base::getParameter(0); }
 
    void setPoissonsRatio(const gsFunctionSet<T> & PoissonsRatio) override { Base::setParameter(1,PoissonsRatio); }
 
    const function_ptr getPoissonsRatio() const override { return Base::getParameter(1); }
 
    std::ostream &print(std::ostream &os) const override;
 
    gsMatrix<T> S(const gsMatrix<T> & strain) const override;
 
    gsMatrix<T> C(const gsMatrix<T> & strain) const override;
 
    gsMatrix<T> S(const gsMatrix<T> & C, const index_t patch, const gsMatrix<T> & u, const gsMatrix<T> & z, enum MaterialOutput out) const;
 
    gsMatrix<T> C(const gsMatrix<T> & C, const index_t patch, const gsMatrix<T> & u, const gsMatrix<T> & z, enum MaterialOutput out) const;
 
    gsMatrix<T> dC(const gsMatrix<T> & C, const index_t patch, const gsMatrix<T> & u, const gsMatrix<T> & z, enum MaterialOutput out) const;
 
public:
    typedef memory::shared_ptr< gsMaterialMatrixLinear > Ptr;
 
    typedef memory::unique_ptr< gsMaterialMatrixLinear > uPtr;
 
protected:
    void _initialize();
 
protected:
    T _Cijkl  (const index_t i, const index_t j, const index_t k, const index_t l) const;
 
    T _Sij    (const index_t i, const index_t j, const gsMatrix<T> & z) const;
 
 
    using Base::_getMetric;
 
    gsMatrix<T> _E      (const T z, enum MaterialOutput out) const;
 
protected:
    // constructor
    using Base::m_thickness;
    using Base::m_pars;
    using Base::m_density;
 
    // Geometric data
    using Base::m_data;
 
    using Base::m_options;
 
};
 
#ifdef GISMO_WITH_PYBIND11
 
  void pybind11_init_gsMaterialMatrixLinear2(pybind11::module &m);
  void pybind11_init_gsMaterialMatrixLinear3(pybind11::module &m);
 
#endif // GISMO_WITH_PYBIND11
 
} // namespace
 
 
#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsMaterialMatrixLinear.hpp)
#endif