skinTissue_Hendriks.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_EXAMPLES_SKIN_HENDRIKS_HH
#define RFFGEN_EXAMPLES_SKIN_HENDRIKS_HH

#include "../../LinearAlgebra/strainTensor.hh"
#include "../../LinearAlgebra/unitMatrix.hh"
#include "../../LinearAlgebra/principalInvariants.hh"
#include "../../generate.hh"

namespace RFFGen
{
  template < class Matrix , int offset = LinearAlgebra::dimension<Matrix>()>
  auto incompressibleSkin_Hendriks(double c0, double c1, const Matrix& F)
  {
    using namespace LinearAlgebra;
    auto i1 = ShiftedFirstPrincipalInvariant<Matrix,offset>();
    auto i2 = ShiftedSecondPrincipalInvariant<Matrix,offset>();

    auto S = LeftCauchyGreenStrainTensor<Matrix>(F);
    return ( c0*i1 + c1*i1*i2 ) << S;

  }

  template <class Matrix , int offset = LinearAlgebra::dimension<Matrix>()>
  auto incompressibleSkin_Hendriks(const Matrix& F)
  {
    return incompressibleSkin_Hendriks<Matrix,offset>(9.4,82.,F);
  }


  template <class InflationPenalty, class CompressionPenalty, class Matrix, int offset = LinearAlgebra::dimension<Matrix>()>
  auto compressibleSkin_Hendriks(double c0, double c1, double d0, double d1, const Matrix& F)
  {
    return incompressibleSkin_Hendriks<Matrix,offset>(c0,c1,F) + volumetricPenalty<InflationPenalty,CompressionPenalty>(d0,d1,F);
  }

  template < class InflationPenalty , class CompressionPenalty , class Matrix , int offset = LinearAlgebra::dimension<Matrix>() >
  auto compressibleSkin_Hendriks(double d0, double d1, const Matrix& F)
  {
    return compressibleSkin_Hendriks<InflationPenalty,CompressionPenalty,Matrix,offset>(9.4,82.,d0,d1,F);
  }
}

#endif // RFFGEN_EXAMPLES_SKIN_HENDRIKS_HH