gsElasticityAssembler.h

Go to the documentation of this file.

 
#pragma once
 
#include <gsElasticity/gsBaseAssembler.h>
#include <gsElasticity/gsElasticityFunctions.h>
#include <gsElasticity/gsBaseUtils.h>
 
namespace gismo
{
 
template <class T>
class gsElasticityAssembler : public gsBaseAssembler<T>
{
public:
    typedef gsBaseAssembler<T> Base;
 
    gsElasticityAssembler(const gsMultiPatch<T> & patches,
                          const gsMultiBasis<T> & basis,
                          const gsBoundaryConditions<T> & bconditions,
                          const gsFunction<T> & body_force);
 
    gsElasticityAssembler(const gsMultiPatch<T> & patches,
                          const gsMultiBasis<T> & basisDisp,
                          const gsMultiBasis<T> & basisPres,
                          const gsBoundaryConditions<T> & bconditions,
                          const gsFunction<T> & body_force);
 
    static gsOptionList defaultOptions();
 
    virtual void refresh();
 
    //--------------------- SYSTEM ASSEMBLY ----------------------------------//
 
    virtual void assemble(bool saveEliminationMatrix);
 
    virtual void assemble() { assemble(false); };
 
    virtual bool assemble(const gsMatrix<T> & solutionVector,
                          const std::vector<gsMatrix<T> > & fixedDoFs);
protected:
    virtual void assemble(const gsMultiPatch<T> & displacement);
 
    virtual void assemble(const gsMultiPatch<T> & displacement, const gsMultiPatch<T> & pressure);
 
    using Base::assemble;
 
    //--------------------- SOLUTION CONSTRUCTION ----------------------------------//
 
public:
 
    using Base::constructSolution;
 
    virtual void constructSolution(const gsMatrix<T> & solVector,
                                   const std::vector<gsMatrix<T> > & fixedDoFs,
                                   gsMultiPatch<T> & displacement) const;
 
    virtual void constructSolution(const gsMatrix<T> & solVector,
                                   const std::vector<gsMatrix<T> > & fixedDoFs,
                                   gsMultiPatch<T> & displacement,
                                   gsMultiPatch<T> & pressure) const;
 
    virtual void constructPressure(const gsMatrix<T> & solVector,
                                   const std::vector<gsMatrix<T> > & fixedDoFs,
                                   gsMultiPatch<T> & pressure) const;
 
    //--------------------- SPECIALS ----------------------------------//
 
    virtual void constructCauchyStresses(const gsMultiPatch<T> & displacement,
                                         gsPiecewiseFunction<T> & result,
                                         stress_components::components component = stress_components::von_mises) const;
 
    virtual void constructCauchyStresses(const gsMultiPatch<T> & displacement,
                                         const gsMultiPatch<T> & pressure,
                                         gsPiecewiseFunction<T> & result,
                                         stress_components::components component = stress_components::von_mises) const;
 
protected:
    virtual void reserve();
 
protected:
    short_t m_dim;
 
    using Base::m_pde_ptr;
    using Base::m_bases;
    using Base::m_ddof;
    using Base::m_options;
    using Base::m_system;
    using Base::eliminationMatrix;
};
 
 
} // namespace gismo ends
 
#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsElasticityAssembler.hpp)
#endif