OgdenElasticity.hpp

Go to the documentation of this file.

#pragma once
 
#include <polyfem/assembler/GenericElastic.hpp>
#include <polyfem/assembler/MatParams.hpp>
 
namespace polyfem::assembler
{
    class UnconstrainedOgdenElasticity : public GenericElastic<UnconstrainedOgdenElasticity>
    {
    public:
        UnconstrainedOgdenElasticity();
 
        // sets material params
        void add_multimaterial(const int index, const json &params, const Units &units) override;
 
        const GenericMatParams &alphas() const { return alphas_; }
        const GenericMatParams &mus() const { return mus_; }
        const GenericMatParams &Ds() const { return Ds_; }
 
        std::string name() const override { return "UnconstrainedOgden"; }
        std::map<std::string, ParamFunc> parameters() const override;
 
        template <typename T>
        T elastic_energy(
            const RowVectorNd &p,
            const double t,
            const int el_id,
            const DefGradMatrix<T> &def_grad) const;
 
    private:
        GenericMatParams alphas_;
        GenericMatParams mus_;
        GenericMatParams Ds_;
    };
 
    class IncompressibleOgdenElasticity : public GenericElastic<IncompressibleOgdenElasticity>
    {
    public:
        IncompressibleOgdenElasticity();
 
        // sets material params
        void add_multimaterial(const int index, const json &params, const Units &units) override;
 
        const GenericMatParams &coefficients() const { return coefficients_; }
        const GenericMatParams &expoenents() const { return expoenents_; }
        const GenericMatParam &bulk_modulus() const { return bulk_modulus_; }
 
        int num_terms() const { return coefficients_.size(); }
 
        std::string name() const override { return "IncompressibleOgden"; }
        std::map<std::string, ParamFunc> parameters() const override;
 
        template <typename T>
        T elastic_energy(
            const RowVectorNd &p,
            const double t,
            const int el_id,
            const DefGradMatrix<T> &def_grad) const;
 
    private:
        GenericMatParams coefficients_;
        GenericMatParams expoenents_;
        GenericMatParam bulk_modulus_;
    };
} // namespace polyfem::assembler
 
#include "OgdenElasticity.tpp"