#include <SchmidCostaExponentialLaw2d.hpp>

Inheritance diagram for SchmidCostaExponentialLaw2d:

Collaboration diagram for SchmidCostaExponentialLaw2d:

Public Member Functions
	SchmidCostaExponentialLaw2d ()

void	ComputeStressAndStressDerivative (c_matrix< double, 2, 2 > &rC, c_matrix< double, 2, 2 > &rInvC, double pressure, c_matrix< double, 2, 2 > &rT, FourthOrderTensor< 2, 2, 2, 2 > &rDTdE, bool computeDTdE)

double	GetA ()

std::vector< std::vector< double > >	GetB ()

double	GetZeroStrainPressure ()

Public Member Functions inherited from AbstractIncompressibleMaterialLaw< 2 >
	AbstractIncompressibleMaterialLaw ()

virtual	~AbstractIncompressibleMaterialLaw ()

Public Member Functions inherited from AbstractMaterialLaw< DIM >
	AbstractMaterialLaw ()

virtual	~AbstractMaterialLaw ()

virtual void	ComputeStressAndStressDerivative (c_matrix< double, DIM, DIM > &rC, c_matrix< double, DIM, DIM > &rInvC, double pressure, c_matrix< double, DIM, DIM > &rT, FourthOrderTensor< DIM, DIM, DIM, DIM > &rDTdE, bool computeDTdE)=0

void	ComputeCauchyStress (c_matrix< double, DIM, DIM > &rF, double pressure, c_matrix< double, DIM, DIM > &rSigma)

void	Compute1stPiolaKirchoffStress (c_matrix< double, DIM, DIM > &rF, double pressure, c_matrix< double, DIM, DIM > &rS)

void	Compute2ndPiolaKirchoffStress (c_matrix< double, DIM, DIM > &rC, double pressure, c_matrix< double, DIM, DIM > &rT)

virtual void	ScaleMaterialParameters (double scaleFactor)

void	SetChangeOfBasisMatrix (c_matrix< double, DIM, DIM > &rChangeOfBasisMatrix)

void	ResetToNoChangeOfBasisMatrix ()

Private Attributes
double	mA

std::vector< std::vector< double > >	mB

c_matrix< double, 2, 2 >	mIdentity

Additional Inherited Members
Protected Member Functions inherited from AbstractMaterialLaw< DIM >
void	ComputeTransformedDeformationTensor (c_matrix< double, DIM, DIM > &rC, c_matrix< double, DIM, DIM > &rInvC, c_matrix< double, DIM, DIM > &rCTransformed, c_matrix< double, DIM, DIM > &rInvCTransformed)

void	TransformStressAndStressDerivative (c_matrix< double, DIM, DIM > &rT, FourthOrderTensor< DIM, DIM, DIM, DIM > &rDTdE, bool transformDTdE)

Protected Attributes inherited from AbstractMaterialLaw< DIM >
c_matrix< double, DIM, DIM > *	mpChangeOfBasisMatrix

Detailed Description

A 2d version of the material law in Costa, Holmes, McCulloch "Modelling Cardiac Mechanical Properties in Three Dimensions" Philo. Trans. R. Soc.

W = a[exp(Q)-1]/2 where Q = bff*Eff^2 + bfs*Efs^2 + bsf*Esf^2 + bss*Ess^2

where the parameters are taken from the fitting in Schmid,Nash,Young,Hunter "Myocardial Material Parameter Estimation - A Comparative Study for Simple Shear" Transactions of the ASME.

Note, by default, the fibre direction is assumed to be THE X-DIRECTION, and the sheet direction the Y-DIRECTION (ie sheets in the XY plane). Call SetChangeOfBasisMatrix() before ComputeStressAndStressDerivative(), with the matrix P = [fibre_vec, sheet_vec, normal_vec] if this is not the case.

Definition at line 58 of file SchmidCostaExponentialLaw2d.hpp.

Constructor & Destructor Documentation

◆ SchmidCostaExponentialLaw2d()

SchmidCostaExponentialLaw2d::SchmidCostaExponentialLaw2d ( )

Constructor.

Definition at line 38 of file SchmidCostaExponentialLaw2d.cpp.

References mA, mB, and mIdentity.

Member Function Documentation

◆ ComputeStressAndStressDerivative()

void SchmidCostaExponentialLaw2d::ComputeStressAndStressDerivative	(	c_matrix< double, 2, 2 > &	rC,
		c_matrix< double, 2, 2 > &	rInvC,
		double	pressure,
		c_matrix< double, 2, 2 > &	rT,
		FourthOrderTensor< 2, 2, 2, 2 > &	rDTdE,
		bool	computeDTdE
	)

Compute the (2nd Piola Kirchoff) stress T and the stress derivative dT/dE for a given strain.

NOTE: the strain E is not expected to be passed in, instead the Lagrangian deformation tensor C is required (recall, E = 0.5(C-I))

dT/dE is a fourth-order tensor, where dT/dE[M][N][P][Q] = dT^{MN}/dE_{PQ}

Parameters

rC	The Lagrangian deformation tensor (F^T F)
rInvC	The inverse of C. Should be computed by the user. (Change this?)
pressure	the current pressure
rT	the stress will be returned in this parameter
rDTdE	the stress derivative will be returned in this parameter, assuming the final parameter is true
computeDTdE	a boolean flag saying whether the stress derivative is required or not.