neoHooke.hh

Go to the documentation of this file.

/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
/*                                                                             */
/*   This file is part of the C++-library RFFGen.                              */
/*   Copyright 2015 Lars Lubkoll                                               */
/*                                                                             */
/*   RFFGen is free software: you can redistribute it and/or modify            */
/*   it under the terms of the GNU General Public License as published by      */
/*   the Free Software Foundation, either version 3 of the License, or         */
/*   (at your option) any later version.                                       */
/*                                                                             */
/*   RFFGen is distributed in the hope that it will be useful,                 */
/*   but WITHOUT ANY WARRANTY; without even the implied warranty of            */
/*   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the             */
/*   GNU General Public License for more details.                              */
/*                                                                             */
/*   You should have received a copy of the GNU General Public License         */
/*   along with RFFGen.  If not, see <http://www.gnu.org/licenses/>.           */
/*                                                                             */
/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */

#ifndef RFFGEN_NEO_HOOKE_HH
#define RFFGEN_NEO_HOOKE_HH

#include "../../LinearAlgebra/invariants.hh"
#include "../../LinearAlgebra/strainTensor.hh"
#include "../../LinearAlgebra/unitMatrix.hh"
#include "../../finalize.hh"
#include "../volumetricPenaltyFunctions.hh"
#include "../../generate.hh"

namespace RFFGen
{
  template < class Matrix , int offset = LinearAlgebra::dimension<Matrix>() >
  auto incompressibleNeoHooke(double c, const Matrix& F)
  {
    using namespace LinearAlgebra;
    using Inv = typename InvariantTraits<Invariant::PRINCIPAL>::template ShiftedFirstInvariant<Matrix,offset>;
    auto S = LeftCauchyGreenStrainTensor<Matrix>(F);
    auto f = c*Inv(S()) << S;
    return Finalize<decltype(f)>(f);
  }

  template < class Matrix , int offset = LinearAlgebra::dimension<Matrix>() >
  auto modifiedIncompressibleNeoHooke(double c, const Matrix& F)
  {
    using namespace LinearAlgebra;
    using Inv = typename InvariantTraits<Invariant::MODIFIED>::template ShiftedFirstInvariant<Matrix,offset>;
    auto S = LeftCauchyGreenStrainTensor<Matrix>(F);
    return c * Inv(S()) << S;
  }

  namespace NeoHookeDetail
  {
    template <class InflationPenalty, class CompressionPenalty, class Matrix, LinearAlgebra::Invariant I = LinearAlgebra::Invariant::PRINCIPAL , int offset = LinearAlgebra::dimension<Matrix>() >
    auto compressibleNeoHookeImpl(double c, double d0, double d1, const Matrix& F)
    {
      using namespace LinearAlgebra;
      using Inv = typename InvariantTraits<I>::template ShiftedFirstInvariant<Matrix,offset>;
      auto S = LeftCauchyGreenStrainTensor<Matrix>(F);
      return ( c*Inv(F) << S ) + volumetricPenalty<InflationPenalty,CompressionPenalty>(d0,d1,F);
      //return Finalize<decltype(f)>(f);
    }
  }
  template <class InflationPenalty, class CompressionPenalty, class Matrix , int offset = LinearAlgebra::dimension<Matrix>() >
  auto compressibleNeoHooke(double c, double d0, double d1, const Matrix& F)
  {
    return NeoHookeDetail::compressibleNeoHookeImpl<InflationPenalty,CompressionPenalty,Matrix,LinearAlgebra::Invariant::PRINCIPAL,offset>(c,d0,d1,F);
  }

  template <class InflationPenalty, class CompressionPenalty, class Matrix , int offset = LinearAlgebra::dimension<Matrix>() >
  auto modifiedCompressibleNeoHooke(double c, double d0, double d1, const Matrix& F)
  {
    return NeoHookeDetail::compressibleNeoHookeImpl<InflationPenalty,CompressionPenalty,Matrix,LinearAlgebra::Invariant::MODIFIED,offset>(c,d0,d1,F);
  }
}

#endif // RFFGEN_NEO_HOOKE_HH