#include <ExplicitCardiacMechanicsAssembler.hpp>
Inherits AbstractCardiacMechanicsAssembler< DIM >.
Public Member Functions | |
ExplicitCardiacMechanicsAssembler (ContractionModel contractionModel, QuadraticMesh< DIM > *pQuadMesh, std::string outputDirectory, std::vector< unsigned > &rFixedNodes, AbstractIncompressibleMaterialLaw< DIM > *pMaterialLaw=NULL) | |
virtual | ~ExplicitCardiacMechanicsAssembler () |
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) |
Friends | |
class | TestExplicitCardiacMechanicsAssembler |
Explicit cardiac mechanics assembler for solving electromechanic problems where the contraction model is not stretch-rate-dependent (for those the implicit assembler is needed).
The general explicit solution procedure is to do, each timestep: (0) [solve the electrics and interpolate Ca and voltage onto quad points (i) pass Ca and voltage to the contraction models (ii) pass the fibre stretch to the contraction models in case this is needed. (iii) integrate the contraction models in order to get the active tension (iv) solve for the deformation using this active tension.
Definition at line 49 of file ExplicitCardiacMechanicsAssembler.hpp.
ExplicitCardiacMechanicsAssembler< DIM >::ExplicitCardiacMechanicsAssembler | ( | ContractionModel | contractionModel, | |
QuadraticMesh< DIM > * | pQuadMesh, | |||
std::string | outputDirectory, | |||
std::vector< unsigned > & | rFixedNodes, | |||
AbstractIncompressibleMaterialLaw< DIM > * | pMaterialLaw = NULL | |||
) | [inline] |
Constructor
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 38 of file ExplicitCardiacMechanicsAssembler.cpp.
References AbstractCardiacMechanicsAssembler< DIM >::mContractionModelSystems, and AbstractCardiacMechanicsAssembler< DIM >::mTotalQuadPoints.
ExplicitCardiacMechanicsAssembler< DIM >::~ExplicitCardiacMechanicsAssembler | ( | ) | [inline, virtual] |
Destructor
Definition at line 99 of file ExplicitCardiacMechanicsAssembler.cpp.
References AbstractCardiacMechanicsAssembler< DIM >::mContractionModelSystems.
void ExplicitCardiacMechanicsAssembler< DIM >::GetActiveTensionAndTensionDerivs | ( | double | currentFibreStretch, | |
unsigned | currentQuadPointGlobalIndex, | |||
bool | assembleJacobian, | |||
double & | rActiveTension, | |||
double & | rDerivActiveTensionWrtLambda, | |||
double & | rDerivActiveTensionWrtDLambdaDt | |||
) | [inline, private, virtual] |
Get the active tension and other info at the given quadrature point. This is an explicit assembler so just sets the active tension, it doesn't set the derivatives. It stores the stretch for the next timestep.
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. Unset in this explicit solver. | |
rDerivActiveTensionWrtDLambdaDt | The returned dT_dLamDot, derivative of active tension wrt stretch rate. Unset in this explicit solver. |
Implements AbstractCardiacMechanicsAssembler< DIM >.
Definition at line 108 of file ExplicitCardiacMechanicsAssembler.cpp.
References AbstractCardiacMechanicsAssembler< DIM >::mContractionModelSystems, and AbstractCardiacMechanicsAssembler< DIM >::mStretches.
bool ExplicitCardiacMechanicsAssembler< DIM >::IsImplicitSolver | ( | ) | [inline, private, virtual] |
This solver is an explicit solver (overloaded pure method)
Implements AbstractCardiacMechanicsAssembler< DIM >.
Definition at line 55 of file ExplicitCardiacMechanicsAssembler.hpp.
void ExplicitCardiacMechanicsAssembler< 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 explicitly over the time range using the ODE timestep provided then the active tension used to solve for the deformation
time | the current time | |
nextTime | the next time | |
odeTimestep | the ODE timestep |
Implements AbstractCardiacMechanicsAssembler< DIM >.
Definition at line 129 of file ExplicitCardiacMechanicsAssembler.cpp.
References NonlinearElasticityAssembler< DIM >::AssembleSystem(), AbstractCardiacMechanicsAssembler< DIM >::mContractionModelSystems, AbstractCardiacMechanicsAssembler< DIM >::mCurrentTime, AbstractCardiacMechanicsAssembler< DIM >::mNextTime, AbstractCardiacMechanicsAssembler< DIM >::mOdeTimestep, and AbstractCardiacMechanicsAssembler< DIM >::mStretches.