gsThinShellFunctions.h

Go to the documentation of this file.

#pragma once

#include <gsKLShell/src/gsThinShellUtils.h>
#include <gsKLShell/src/gsMaterialMatrixEval.h>
#include <gsKLShell/src/gsMaterialMatrixIntegrate.h>
#include <gsKLShell/src/gsMaterialMatrixLinear.h>
#include <gsKLShell/src/gsMaterialMatrixNonlinear.h>

#include <gsAssembler/gsExprAssembler.h>


namespace gismo
{

struct stress_type
{
    enum type
    {
        displacement       = -1,
        von_mises          = 0,  
        von_mises_membrane = 1,  
        von_mises_flexural = 2,  
        membrane           = 3,  
        membrane_force     = 4,  
        membrane_PK2       = 5,  
        membrane_force_PK2 = 6,  
        flexural           = 7,  
        flexural_PK2       = 8,  
        flexural_moment    = 9,  
        flexural_moment_PK2= 10, 
        total              = 11, 
        membrane_strain    = 12, 
        flexural_strain    = 13, 
        principal_membrane_strain    = 14,  
        principal_flexural_strain    = 15,  
        principal_stretch  = 16, 
        principal_stress   = 17,  
        principal_stress_membrane  = 18,  
        principal_stress_flexural  = 19,  
        principal_stretch_dir1  = 81,  
        principal_stretch_dir2  = 82,  
        principal_stretch_dir3  = 83,  
        principal_stress_dir1  = 91,  
        principal_stress_dir2  = 92,  
        principal_stress_dir3  = 93,  
        tension_field  = 99,
    };
};

template <class T>
class gsShellStressFunction : public gsFunction<T>
{
public:

    gsShellStressFunction(const gsFunctionSet<T> & geometry,
                           const gsFunctionSet<T> & deformed,
                           const gsMaterialMatrixContainer<T> & mm,
                           index_t patch,
                           stress_type::type type)
        : m_patches(&geometry),
          m_defpatches(&deformed),
          m_materialMatrices(mm),
          m_patchID(patch),
          m_stress_type(type)
    {

    }

    ~gsShellStressFunction()
    {
        // delete m_patches;
        // delete m_defpatches;
    }

    virtual short_t domainDim() const
    {
        return 2;
    }

    virtual short_t targetDim() const
    {
        switch (m_stress_type)
        {
            default:
                return 0;
                break;
            case stress_type::displacement :
                return 3;
                break;

            case stress_type::membrane :
                return 3;
                break;

            case stress_type::flexural :
                return 3;
                break;

            case stress_type::membrane_PK2 :
                return 3;
                break;

            case stress_type::flexural_PK2 :
                return 3;
                break;

            case stress_type::membrane_force :
                return 3;
                break;

            case stress_type::flexural_moment :
                return 3;
                break;

            case stress_type::membrane_force_PK2 :
                return 3;
                break;

            case stress_type::flexural_moment_PK2 :
                return 3;
                break;

            // TO BE IMPLEMENTED
            // -------------------------------------
            case stress_type::von_mises :
                return 2;
                break;

            case stress_type::von_mises_membrane :
                return 1;
                break;

            case stress_type::von_mises_flexural :
                return 1;
                break;
            // -------------------------------------

            case stress_type::membrane_strain :
                return 3;
                break;
            case stress_type::principal_membrane_strain :
                return 3;
                break;
            case stress_type::principal_flexural_strain :
                return 3;
                break;

            case stress_type::flexural_strain :
                return 3;
                break;
            case stress_type::principal_stretch :
                return 3;
                break;
            case stress_type::principal_stress:
                return 2;
                break;
            case stress_type::principal_stress_membrane :
                return 2;
                break;
            case stress_type::principal_stress_flexural :
                return 2;
                break;
            case stress_type::principal_stretch_dir1 :
                return 3;
                break;
            case stress_type::principal_stretch_dir2 :
                return 3;
                break;
            case stress_type::principal_stretch_dir3 :
                return 3;
                break;
            case stress_type::principal_stress_dir1 :
                return 3;
                break;
            case stress_type::principal_stress_dir2 :
                return 3;
                break;
            case stress_type::principal_stress_dir3 :
                return 3;
            case stress_type::tension_field :
                return 1;
                break;
            /*
                DEFAULT includes:
                stress_type::membrane;
                stress_type::flexural;
                stress_type::membrane_strain;
                stress_type::flexural_strain;
                stress_type::principal_stretch;
                stress_type::principal_stretch_dir1;
                stress_type::principal_stretch_dir2;
                stress_type::principal_stretch_dir3
            */
        }
    }

    virtual void eval_into(const gsMatrix<T> & u, gsMatrix<T> & result) const;

protected:

    typedef gsExprAssembler<>::geometryMap geometryMap;
    typedef gsExprAssembler<>::variable    variable;
    typedef gsExprAssembler<>::space       space;
    typedef gsExprAssembler<>::solution    solution;

    const gsFunctionSet<T> * m_patches;
    const gsFunctionSet<T> * m_defpatches;
    const gsMaterialMatrixContainer<T> & m_materialMatrices;
    index_t m_patchID;
    stress_type::type m_stress_type;

}; // class definition ends


} // namespace ends


#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsThinShellFunctions.hpp)
#endif