muscleTissue_Martins.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_MUSCLE_TISSUE_MARTINS_HH
#define RFFGEN_MUSCLE_TISSUE_MARTINS_HH

#include "../../LinearAlgebra/modifiedPrincipalInvariants.hh"
#include "../../LinearAlgebra/modifiedMixedInvariants.hh"
#include "../../LinearAlgebra/strainTensor.hh"
#include "../../LinearAlgebra/tensorProduct.hh"
#include "../../CMath/exp.hh"
#include "../volumetricPenaltyFunctions.hh"
#include "../../generate.hh"

namespace RFFGen
{
  namespace MuscleTissueDetail
  {
    template < class Matrix , int offset = LinearAlgebra::dimension<Matrix>() >
    auto generateIncompressibleMuscleTissue_Martins(double c, double b, double A, double a, const Matrix& M, const Matrix& F )
    {
      using CMath::exp;
      using namespace LinearAlgebra;
      auto i1 = ShiftedFirstModifiedPrincipalInvariant<Matrix,offset>();
      auto i6 = ShiftedThirdModifiedMixedInvariant<Matrix>(F,M);

      return ( c * ( exp( b * i1 ) - 1 ) +
               A * ( exp( a * ( i6^2 ) ) - 1 )
             ) << LeftCauchyGreenStrainTensor<Matrix>(F);
    }
  }
  template < class Matrix , int offset = LinearAlgebra::dimension<Matrix>()>
  auto incompressibleMuscleTissue_Martins(double c, double b, double A, double a, const Matrix& M, const Matrix& F)
  {
    return MuscleTissueDetail::generateIncompressibleMuscleTissue_Martins<Matrix,offset>(c,b,A,a,M,F);
  }

  template < class Matrix , int offset = LinearAlgebra::dimension<Matrix>()>
  auto incompressibleMuscleTissue_Martins(const Matrix& M, const Matrix& F)
  {
    return incompressibleMuscleTissue_Martins<Matrix,offset>(0.387, 23.46, 0.584, 12.43, M, F);
  }

  template < class Inflation , class Compression , class Matrix , int offset = LinearAlgebra::dimension<Matrix>()>
  auto compressibleMuscleTissue_Martins(double c, double b, double A, double a, double d0, double d1, const Matrix& M, const Matrix& F)
  {
    return MuscleTissueDetail::generateIncompressibleMuscleTissue_Martins<Matrix,offset>(c,b,A,a,M,F) + volumetricPenalty<Inflation,Compression>(d0,d1,F);
  }

  template < class Inflation , class Compression , class Matrix , int offset = LinearAlgebra::dimension<Matrix>()>
  auto compressibleMuscleTissue_Martins(double d0, double d1, const Matrix& M, const Matrix& F)
  {
    return compressibleMuscleTissue_Martins<Inflation,Compression,Matrix,offset>(0.387, 23.46, 0.584, 12.43, d0, d1, M, F);
  }
}

#endif // RFFGEN_MUSCLE_TISSUE_MARTINS_HH