ImplicitCardiacMechanicsAssembler< DIM > Class Template Reference

#include <ImplicitCardiacMechanicsAssembler.hpp>

Inherits AbstractCardiacMechanicsAssembler< DIM >.

Collaboration diagram for ImplicitCardiacMechanicsAssembler< DIM >:
Collaboration graph
[legend]

List of all members.

Public Member Functions

 ImplicitCardiacMechanicsAssembler (ContractionModel contractionModel, QuadraticMesh< DIM > *pQuadMesh, std::string outputDirectory, std::vector< unsigned > &rFixedNodes, AbstractIncompressibleMaterialLaw< DIM > *pMaterialLaw=NULL)
virtual ~ImplicitCardiacMechanicsAssembler ()
std::vector< double > & rGetFibreStretches ()
void Solve (double time, double nextTime, double odeTimestep)

Private Member Functions

bool IsImplicitSolver ()
void GetActiveTensionAndTensionDerivs (double currentFibreStretch, unsigned currentQuadPointGlobalIndex, bool assembleJacobian, double &rActiveTension, double &rDerivActiveTensionWrtLambda, double &rDerivActiveTensionWrtDLambdaDt)

Private Attributes

std::vector< double > mStretchesLastTimeStep

Friends

class TestImplicitCardiacMechanicsAssembler

Detailed Description

template<unsigned DIM>
class ImplicitCardiacMechanicsAssembler< DIM >

Implicit Cardiac Mechanics Assembler

Solves cardiac mechanics implicitly (together with the contraction models for determining the active tension), taking in the intracellular calcium concentration. See CardiacElectroMechanicsProblem documentation for more detail.

Definition at line 51 of file ImplicitCardiacMechanicsAssembler.hpp.

Constructor & Destructor Documentation

template<unsigned DIM>
ImplicitCardiacMechanicsAssembler< DIM >::ImplicitCardiacMechanicsAssembler ( ContractionModel  contractionModel,
QuadraticMesh< DIM > *  pQuadMesh,
std::string  outputDirectory,
std::vector< unsigned > &  rFixedNodes,
AbstractIncompressibleMaterialLaw< DIM > *  pMaterialLaw = NULL 
) [inline]

Constructor

Parameters:
contractionModel The contraction model.
pQuadMesh A pointer to the mesh.
outputDirectory The output directory, relative to TEST_OUTPUT
rFixedNodes The fixed nodes
pMaterialLaw The material law for the tissue. Defaults to NULL, in which case a default material law is used.

Definition at line 35 of file ImplicitCardiacMechanicsAssembler.cpp.

References AbstractCardiacMechanicsAssembler< DIM >::mContractionModelSystems, ImplicitCardiacMechanicsAssembler< DIM >::mStretchesLastTimeStep, and AbstractCardiacMechanicsAssembler< DIM >::mTotalQuadPoints.

template<unsigned DIM>
ImplicitCardiacMechanicsAssembler< DIM >::~ImplicitCardiacMechanicsAssembler (  )  [inline, virtual]

Destructor

Definition at line 88 of file ImplicitCardiacMechanicsAssembler.cpp.

References AbstractCardiacMechanicsAssembler< DIM >::mContractionModelSystems.

Member Function Documentation

template<unsigned DIM>
void ImplicitCardiacMechanicsAssembler< DIM >::GetActiveTensionAndTensionDerivs ( double  currentFibreStretch,
unsigned  currentQuadPointGlobalIndex,
bool  assembleJacobian,
double &  rActiveTension,
double &  rDerivActiveTensionWrtLambda,
double &  rDerivActiveTensionWrtDLambdaDt 
) [inline, private, virtual]

A method called by AbstractCardiacMechanicsAssembler::AssembleOnElement(), providing the active tension (and other info) at a particular quadrature point. This version uses C to determine the current stretch and stretch rate, and integrates the contraction model ODEs to determine the active tension, and derivatives of the active tension with respect to stretch and stretch rate.

Parameters:
currentFibreStretch The stretch in the fibre direction
currentQuadPointGlobalIndex Quadrature point integrand currently being evaluated at in AssembleOnElement.
assembleJacobian A bool stating whether to assemble the Jacobian matrix.
rActiveTension The returned active tension
rDerivActiveTensionWrtLambda The returned dT_dLam, derivative of active tension wrt stretch
rDerivActiveTensionWrtDLambdaDt The returned dT_dLamDot, derivative of active tension wrt stretch rate

Implements AbstractCardiacMechanicsAssembler< DIM >.

Definition at line 136 of file ImplicitCardiacMechanicsAssembler.cpp.

References AbstractContractionModel::GetNextActiveTension(), AbstractCardiacMechanicsAssembler< DIM >::mContractionModelSystems, AbstractCardiacMechanicsAssembler< DIM >::mCurrentTime, AbstractCardiacMechanicsAssembler< DIM >::mNextTime, AbstractCardiacMechanicsAssembler< DIM >::mOdeTimestep, AbstractCardiacMechanicsAssembler< DIM >::mStretches, ImplicitCardiacMechanicsAssembler< DIM >::mStretchesLastTimeStep, AbstractContractionModel::RunDoNotUpdate(), and AbstractContractionModel::SetStretchAndStretchRate().

template<unsigned DIM>
bool ImplicitCardiacMechanicsAssembler< DIM >::IsImplicitSolver (  )  [inline, private, virtual]

The current stretch ratio (in the fibre direction). Note the indexing: the i-th entry corresponds to the i-th global quad point, when looping over elements and then quad points This solver is an implicit solver (overloaded pure method)

Implements AbstractCardiacMechanicsAssembler< DIM >.

Definition at line 69 of file ImplicitCardiacMechanicsAssembler.hpp.

template<unsigned DIM>
std::vector< double > & ImplicitCardiacMechanicsAssembler< DIM >::rGetFibreStretches (  )  [inline]

Get lambda (the stretch ratio). NOTE: the i-th entry of this vector is assumed to be the i-th quad point obtained by looping over cells in the obvious way and then looping over quad points. These quad points, in the same order, can be obtained by using the QuadraturePointsGroup class.

Definition at line 98 of file ImplicitCardiacMechanicsAssembler.cpp.

References AbstractCardiacMechanicsAssembler< DIM >::mStretches.

template<unsigned DIM>
void ImplicitCardiacMechanicsAssembler< DIM >::Solve ( double  time,
double  nextTime,
double  odeTimestep 
) [inline, virtual]

Solve for the deformation using quasi-static nonlinear elasticity. (not dynamic nonlinear elasticity, despite the times taken in - just ONE deformation is solved for. The cell models are integrated implicitly over the time range using the ODE timestep provided, as part of the solve, and updated at the end once the solution has been found, as is lambda.

Parameters:
time the current time
nextTime the next time
odeTimestep the ODE timestep

Implements AbstractCardiacMechanicsAssembler< DIM >.

Definition at line 105 of file ImplicitCardiacMechanicsAssembler.cpp.

References NonlinearElasticityAssembler< DIM >::AssembleSystem(), AbstractNonlinearElasticityAssembler< DIM >::GetNumNewtonIterations(), AbstractCardiacMechanicsAssembler< DIM >::mContractionModelSystems, AbstractCardiacMechanicsAssembler< DIM >::mCurrentTime, AbstractCardiacMechanicsAssembler< DIM >::mNextTime, AbstractCardiacMechanicsAssembler< DIM >::mOdeTimestep, AbstractCardiacMechanicsAssembler< DIM >::mStretches, and ImplicitCardiacMechanicsAssembler< DIM >::mStretchesLastTimeStep.

Member Data Documentation

template<unsigned DIM>
std::vector<double> ImplicitCardiacMechanicsAssembler< DIM >::mStretchesLastTimeStep [private]

The stretch in the fibre direction at the last timestep, in order to compute the stretch rate. Note the indexing: the i-th entry corresponds to the i-th global quad point, when looping over elements and then quad points

Definition at line 61 of file ImplicitCardiacMechanicsAssembler.hpp.

Referenced by ImplicitCardiacMechanicsAssembler< DIM >::GetActiveTensionAndTensionDerivs(), ImplicitCardiacMechanicsAssembler< DIM >::ImplicitCardiacMechanicsAssembler(), and ImplicitCardiacMechanicsAssembler< DIM >::Solve().

The documentation for this class was generated from the following files:
Generated by  doxygen 1.6.2