gsMuscleAssembler.h

Go to the documentation of this file.

 
#pragma once
 
#include <gsElasticity/gsElasticityAssembler.h>
 
namespace gismo
{
 
template <class T>
class gsMuscleAssembler : public gsElasticityAssembler<T>
{
public:
    typedef gsElasticityAssembler<T> Base;
 
    gsMuscleAssembler(const gsMultiPatch<T> & patches,
                      const gsMultiBasis<T> & basisDisp,
                      const gsMultiBasis<T> & basisPres,
                      const gsBoundaryConditions<T> & bconditions,
                      const gsFunction<T> & body_force,
                      const gsPiecewiseFunction<T> & tendonMuscleDistribution,
                      const gsVector<T> & fiberDirection);
 
    //--------------------- SYSTEM ASSEMBLY ----------------------------------//
 
    using Base::assemble;
    virtual bool assemble(const gsMatrix<T> & solutionVector,
                          const std::vector<gsMatrix<T> > & fixedDoFs);
 
    //--------------------- SPECIALS ----------------------------------//
 
    using Base::constructCauchyStresses;
    virtual void constructCauchyStresses(const gsMultiPatch<T> & displacement,
                                         const gsMultiPatch<T> & pressure,
                                         gsPiecewiseFunction<T> & result,
                                         stress_components::components component = stress_components::von_mises) const;
 
protected:
    using Base::m_options;
    using Base::m_pde_ptr;
    using Base::m_system;
    gsPiecewiseFunction<T> const & muscleTendon;
    gsVector<T> const & fiberDir;
};
 
 
} // namespace gismo ends
 
#ifndef GISMO_BUILD_LIB
#include GISMO_HPP_HEADER(gsMuscleAssembler.hpp)
#endif