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 ----------------------------------//

    virtual bool assemble(const gsMatrix<T> & solutionVector,
                          const std::vector<gsMatrix<T> > & fixedDoFs);

    //--------------------- SPECIALS ----------------------------------//

    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