adiposeTissue_SommerHolzapfel.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_ADIPOSE_TISSUE_SOMMER_HOLZAPFEL_HH
#define RFFGEN_ADIPOSE_TISSUE_SOMMER_HOLZAPFEL_HH

#include "../../LinearAlgebra/principalInvariants.hh"
#include "../../LinearAlgebra/mixedInvariants.hh"
#include "../../LinearAlgebra/strainTensor.hh"
#include "../../LinearAlgebra/tensorProduct.hh"
#include "../../CMath/exp.hh"
#include "../../generate.hh"

namespace RFFGen
{
  namespace Detail
  {
    template < class Matrix , int offset = LinearAlgebra::dimension<Matrix>() >
    auto generateIncompressibleAdiposeTissue_SommerHolzapfel(double cCells, double k1, double k2, double kappa, const Matrix& M, const Matrix& F)
    {
      using CMath::exp;
      using namespace LinearAlgebra;
      auto i1 = FirstPrincipalInvariant<Matrix>(F);
      auto si1 = ShiftedFirstPrincipalInvariant<Matrix,offset>(F);
      auto i4 = FirstMixedInvariant<Matrix>(F,M);

      auto aniso = kappa*i1 + (1-3*kappa)*i4 - 1;
      return ( cCells*si1 +
               (k1/k2)*( exp( k2*(aniso^2) ) - 1)
             ) << LeftCauchyGreenStrainTensor<Matrix>(F);
    }
  }
  template < class Matrix , int offset = LinearAlgebra::dimension<Matrix>()>
  auto incompressibleAdiposeTissue_SommerHolzapfel(double cCells, double k1, double k2, double kappa, const Matrix& M, const Matrix& F)
  {
    return Detail::generateIncompressibleAdiposeTissue_SommerHolzapfel<Matrix,offset>(cCells,k1,k2,kappa,M,F);
  }

  template < class Matrix , int offset = LinearAlgebra::dimension<Matrix>()>
  auto incompressibleAdiposeTissue_SommerHolzapfel(const Matrix& M, const Matrix& F)
  {
    return incompressibleAdiposeTissue_SommerHolzapfel<Matrix,offset>(0.15,0.8,47.3,0.09,M,F);
  }


  template <class Inflation, class Compression, class Matrix , int offset = LinearAlgebra::dimension<Matrix>()>
  auto compressibleAdiposeTissue_SommerHolzapfel(double cCells, double k1, double k2, double kappa, double d0, double d1, const Matrix& M, const Matrix& F)
  {
    return Detail::generateIncompressibleAdiposeTissue_SommerHolzapfel<Matrix,offset>(cCells,k1,k2,kappa,M,F) + volumetricPenalty<Inflation,Compression>(d0,d1,F);
  }

  template <class Inflation, class Compression, class Matrix , int offset = LinearAlgebra::dimension<Matrix>()>
  auto compressibleAdiposeTissue_SommerHolzapfel(double d0, double d1, const Matrix& M, const Matrix& F)
  {
    return compressibleAdiposeTissue_SommerHolzapfel<Inflation,Compression,Matrix,offset>(0.15,0.8,47.3,0.09,d0,d1,M,F);
  }
}

#endif // RFFGEN_ADIPOSE_TISSUE_SOMMER_HOLZAPFEL_HH