Inflation of a Passive Prolate Spheroid

Description

  • This example guides you through setting up and running a physical and material nonlinear problem of the passive (non-contracting) inflation of a truncated ellipsoid using pressure (natural) boundary conditions. This mesh is a rough approximation of a left ventricle.
  • The mesh consists of three elements that are wrapped onto themselves in a prolate spheroidal coordinate system. The basis functions are linear in the circumferential direction and cubic in the longitudinal and radial directions (linear-cubic-cubic). The material chosen is a transversely isotropic exponential strain energy function. The material is stiffer in the fiber direction than perpendicular to them. The fibers whose orientation is defined by fiber angles in the node form with respect to the circumferential direction vary linearly in radial direction across the wall from -37 degrees on the epicardium to +84 degrees on the endocardium. The boundary conditions assure that all rigid body motions are suppressed. Pressure is prescribed on the endocardium to (partly) simulate the passive filling of a ventricle (diastole).
  • The final (refined) model nodes and elements are defined according to the figure below:
  • Stress and strain solution tables for various pressure loads and complete model files are provided at the end of this tutorial.

Building and solving the model

Start Continuity

  • Launch the Continuity Client
  • On the About Continuity startup screen

    • Leave the mesh checkbox checked under Use Modules:

    • In addition, check the biomechanics checkbox

  • Click OK to bring up the main window

Create mesh

Refine the Mesh

To get a sufficiently converged result using linear elements, it is necessary to use multiple elements. Therefore, we will refine our single element into many elements.

Add biomechanics data

  • Load the required biomechanics model from the database
    • File→Library→Search...
    • In the window near the top, enter 'lagrangian' and hit return.

    • From the listed models select BM_TI_of_Lagrangian_strains_comp_sympy by right-clicking on it and selecting 'Load'

    • When the warning window display, select the third choice: 'Retain current problem but overwrite the following objects: [dims, renderer, matEquations]'
  • Biomechanics→Update→Initial Conditions with undeformed nodes

  • Biomechanics→Edit→Boundary Conditions...

    • Click on the Deformed Coordinates 1 tab

      • Click the Insert Nodes button three times

      • Edit the `Node(s)' and 'Derivative' parameters of the three nodes based on the table below
        • Node in Nodes List

          Node(s)

          Derivative

          Value

          1

          APEX

          wrt s(2)

          0

          2

          APEX

          wrt s(2)s(3)

          0

          3

          APEX

          Value

          0

    • Click on the Deformed Coordinates 2 tab

      • Click the Insert Nodes button six times

      • Edit the `Node(s)' and 'Derivative' parameters of the three nodes based on the table below
        • Node in Nodes List

          Node(s)

          Derivative

          Value

          1

          APEX

          Value

          0

          2

          APEX

          wrt s(3)

          0

          3

          APEX

          wrt s(2)s(3)

          0

          4

          BASE

          Value

          0

          5

          BASE

          wrt s(3)

          0

          6

          BASE

          wrt s(2)s(3)

          0

    • Click on the Deformed Coordinates 3 tab

      • Click the Insert Nodes button once

      • Enter BASE_EPI$ in the `Node(s)' text field

      • For Derivative, choose Value. Under Value enter 0.

    • Click on the External Pressure tab

      • For Element Number, enter ENDO, and hit the Enter key on your keyboard

      • For the Select Pressure Type drop-down list, choose Incremental

      • In the Inner surface box on the bottom, enter 1.0 for the Specify Pressure Increment field

    • It is a good idea to now go back through the Boundary Conditions Form to double check this parameters you just set up
    • Click the OK button

Solve the biomechanics

Calculate and render stress and strain

  • Mesh→Render→Elements...

    • Click the lines radio button

    • This time select the deformed radio button

    • Click Render to display deformed mesh lines

  • View→Show→OpenMesh...

    • Click on 3. element lines3 in the list on the left, and enter 1,0,0 in the R,G,B entry field to change the mesh lines from blue to red.

    • Press [return] and close the window
    • The mesh should look like the screenshot below
  • Biomechanics→Render→Surface...

    • Change the pop-up menu choice after At Xi to 1 and change Location to 0.5

    • Select the deformed radio button

    • under Variables select T Cauchy Stress Tensor

    • Click OK to render the fiber (circumferential) stress, OR select [2,2] to render the radial stress.

Solution files for various pressure loads

Pre-built model

This cont6 file contains all data and parameters for this problem: bm2.cont6 (original), bm2_refined.cont6 (refined)