gsMaterialMatrixLinear.h

Go to the documentation of this file.

#pragma once

#include <gsKLShell/src/gsMaterialMatrixBase.h>
#include <gsKLShell/src/gsMaterialMatrixBaseDim.h>
#include <gsKLShell/src/gsMaterialMatrixUtils.h>
#include <gsIO/gsOptionList.h>
#include <gsCore/gsFuncData.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:

    gsMaterialMatrixLinear( const gsMaterialMatrixLinear<dim,T> & other);

    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;

    gsMatrix<T> eval3D_dmatrix(const index_t patch, const gsMatrix<T> & u, const gsMatrix<T>& z, enum MaterialOutput out = MaterialOutput::Generic) const override;

    gsMatrix<T> eval3D_vector (const index_t patch, const gsMatrix<T> & u, const gsMatrix<T>& z, enum MaterialOutput out = MaterialOutput::Generic) const override;

    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;

    gsMatrix<T> eval3D_CauchyStress(const index_t patch, const gsMatrix<T> & u, const gsMatrix<T> & z, enum MaterialOutput out) const override;

    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