Created by: Martyn Nash and David Nickerson June 2000

WARNING: This example does not run correctly!

Dynamic HMT was not originally with compression in mind (model validation was based on tension experiments), but in any 3D incompressible model, some structures will experience compressive deformations. As a result, the history terms in the HMT fading memory description (Q), become unstable and hence the model solutions become unreliable.

Sets up an LV mechanics model coupled to a HMT cell model, from which the active fibre tension is derived. The LV, made up of two high order prolate elements, has been previously inflated (passively) to 1kPa. The solution phase for active contraction uses a [Ca]i pulse (defined using a time variable) to drive the HMT model, which computes the active tension at all computational points, causing the LV to contract.