Example 144: Auckland's rectangular cartesian anatomical mesh of canine ventricles

Reads geometric information from ipfiles for Auckland's rectangular Cartesian mesh of canine ventricles and exports node and element files to cmgui.

 Heart H27 nodal parameters. 
 (Note: Fibre and sheet angles are taken from the prolate spheroidal
        mesh.  As node positions are slightly different, refitting the fibre
        and sheet data may improve the model.
        Geometry d2/ds1/ds2 cross derivatives at node 41 have been hand
        editted to make the mesh smooth in the endo-apical region) 

 There are 99 nodes which are used in 60 elements using the 
 rectangular Cartesian coordinate system (x,y,z). 
 See Nielsen, LeGrice, Smaill and Hunter "Mathematical model of 
 geometry and fibrous structure of the heart" AJP 260 H1365-H1378, 1991.

 Coordinates:
   x-axis is main axis of heart from base to apex,
   y-axis goes through middle of RV,
   z-axis completes right-handed orthogonal system.

 Bases used are:
   Basis nb=1 is tricubic (used for x,y,z),
   Basis nb=2 is trilinear (unused for constant imbrication field),
   Basis nb=3 is bilinear-cubic (used for fibre field),
   Basis nb=8 is linear-bicubic (used for sheet field),
   Bases nb=4,5,6,7 are for faces.

 At nodes:
   Xj(1) is the x coord,
   Xj(2) is the y coord,
   Xj(3) is the z coord.
   Derivatives are wrt arc length,
   direction 1 is the direction of the 1st local element coord,
   direction 2 is the 2nd local element coord,
   direction 3 is the 3rd local element coord.

 Nodes on the epicardium are, in layers from base to apex and in
 clockwise order (viewed from base towards apex) starting at z=0:
     Layer 1           Layer 2           Layer 3           Apex
  node   version    node   version    node   version    node   version
       xyz fib she       xyz fib she       xyz fib she       xyz fib she 
    1   1   1   1     5   1   1   1     9   1   1   1    13   2   1   1
   44   1   1   1    43   1   1   1    42   1   1   1    13   3   1   1
    2   1   1   1     6   1   1   1    10   1   1   1    13   4   1   1
   56   1   1   1    55   1   1   1    54   1   1   1    13   5   1   1
   47   1   1   1    46   1   1   1    45   1   1   1    13   6   1   1
    3   1   1   1     7   1   1   1    11   1   1   1    13   7   1   1
   50   1   1   1    49   1   1   1    48   1   1   1    13   8   1   1
   59   1   1   1    58   1   1   1    57   1   1   1    13   9   1   1
    4   1   1   1     8   1   1   1    12   1   1   1    13  10   1   1
   53   1   1   1    52   1   1   1    51   1   1   1    13   1   1   1

 Nodes in the LV midwall are, in layers from base to apex and in
 clockwise order (viewed from base towards apex):
     Layer 1           Layer 2           Layer 3           Apex
  node   version    node   version    node   version    node   version
       xyz fib she       xyz fib she       xyz fib she       xyz fib she 
  (RV)              (RV)               24  2,1  1   1    28   2   1   1
  (RV)              (RV)               78  2,1  1   1    28   3   1   1
   15  3,1  1   1    20   1   1   1    25  2,1  1   1    28   4   1   1
   90   1   1   1    89   1   1   1    88   1   1   1    28   5   1   1
   93   1   1   1    92   1   1   1    91   1   1   1    28   6   1   1
   16   1   1   1    21   1   1   1    26   1   1   1    28   7   1   1
   96   1   1   1    95   1   1   1    94   1   1   1    28   8   1   1
   99   1   1   1    96   1   1   1    97   1   1   1    28   9   1   1
   17  4,2  1   1    22   2   1   1    27  2,1  1   1    28  10   1   1
  (RV)              (RV)               81  2,1  1   1    28   1   1   1

 Nodes on the RV free wall endo are, in layers from base to apex and in
 clockwise order (viewed from base towards apex):
     Layer 1           Layer 2           Layer 3
  node   version    node   version    node   version
       xyz fib she       xyz fib she       xyz fib she
   17   3   2   1    22   1   2   1    27   2   2   1
   83   1   1   1    82   1   1   1    81   2   2   1
   14   1   1   1    19   1   1   1    24   2   2   1
   80   1   1   1    79   1   1   1    78   2   2   1
   15   4   2   1    20   1   2   1    25   2   2   1

 Nodes on the RV septal wall endo are, in layers from base to apex and in
 clockwise order (viewed from base towards apex):
     Layer 1           Layer 2           Layer 3           
  node   version    node   version    node   version
       xyz fib she       xyz fib she       xyz fib she
   17   1   3   1    22   1   3   1    27   2   3   1
   87   1   1   1    22   1   1   1    81   2   3   1
   18   1   1   1    23   1   1   1    24   2   3   1
   85   1   1   1    84   1   1   1    78   2   3   1
   15   2   3   1    20   1   3   1    25   2   3   1

 Nodes on the LV endocardium are, in layers from base to apex and in
 clockwise order (viewed from base towards apex) starting at z=0:
     Layer 1           Layer 2           Layer 3           Apex
  node   version    node   version    node   version    node   version
       xyz fib she       xyz fib she       xyz fib she       xyz fib she 
   29   1   1   1    33   1   1   1    37   1   1   1    41   2   1   1
   62   1   1   1    61   1   1   1    60   1   1   1    41   3   1   1
   30   1   1   1    34   1   1   1    38   1   1   1    41   4   1   1
   74   1   1   1    73   1   1   1    72   1   1   1    41   5   1   1
   65   1   1   1    64   1   1   1    63   1   1   1    41   6   1   1
   31   1   1   1    35   1   1   1    39   1   1   1    41   7   1   1
   68   1   1   1    67   1   1   1    66   1   1   1    41   8   1   1
   77   1   1   1    76   1   1   1    75   1   1   1    41   9   1   1
   32   1   1   1    36   1   1   1    40   1   1   1    41  10   1   1
   71   1   1   1    70   1   1   1    69   1   1   1    41   1   1   1

 The outer ring of elements are:
   LV free wall (subepi)    RV free wall
        3.. 8                9,10, 1, 2
       13..18               19,20,11,12
       23..28,29,30,21,22  

 The inner ring of elements are:
   LV free wall (subendo)   septal wall
       33..38               39,40,31,32
       43..48               49,50,41,42
       53..58,59,60,51,52

 Nodal derivatives are multiplied by scale factors to provide derivatives wrt
 local element coordinates.  The scale factors are node based for first
 derivative continuity where possible, and calculated from the harmonic mean of
 the lengths of adjacent arcs in the corresponding direction.

The mesh may be viewed after running this example by running the CMGUI comfile viewheart.


The comfile run by this example is as follows:

#Example_144: Reads all ipfiles for Auckland's rectangular Cartesian
#full anatomical model of the ventricles (geometry fibres and sheets). 
#Exports of node and element files to cmgui.

fem
fem define coord;r;rc3d;example           #3D rectangular Cartesian coordinates
fem define node;r;heart;example           #99 nodes
fem define base;r;heart;example    
fem define element;r;heart;example        #60 tricubic-Hermite elements
fem define fibre;r;heart;example          #fibre angle: bilinear/cubic-Hermite
fem define element;r;heart;example fibre  #imbrication angle: zero everywhere 
                                          #sheet angle: linear/bicubic-Hermite
fem export nodes;heart as heart           #export nodes to cmgui
fem export elements;heart as heart        #export elements to cmgui

Additional testing commands:

fem list element total

Files used by this example are:

Name              Modified     Size

example_144.com 10-Apr-2000 801 heart-fit.ipnode 10-Apr-2000 268k heart.ipbase 10-Apr-2000 11k heart.ipelem 10-Apr-2000 64k heart.ipelfb 10-Apr-2000 27k heart.ipfibr 30-Jan-2001 49k heart.ipnode 10-Apr-2000 267k rc3d.ipcoor 10-Apr-2000 728 test_output.com 10-Apr-2000 23 viewheart.com 09-Nov-2001 624

Download the entire example:

Name                      Modified     Size

examples_1_14_144.tar.gz 18-Aug-2006 115k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:01:35 20161
cm-debugSuccessSat Mar 5 00:01:54 20161
mips-irix
cmSuccessSun Aug 19 01:32:54 20076
cm-debugSuccessWed Aug 15 01:27:29 200710
cm-debug-clear-mallocSuccessSat Aug 18 01:31:43 200715
cm-debug-clear-malloc7SuccessMon Aug 20 01:28:58 200714
cm64SuccessSun Aug 19 01:32:56 20076
cm64-debugSuccessTue Aug 21 01:26:49 200711
cm64-debug-clear-mallocSuccessThu Apr 1 10:10:27 20045
rs6000-aix
cmSuccessWed Mar 4 01:08:41 20091
cm-debugSuccessMon Mar 2 01:07:18 20092
cm64SuccessWed Mar 4 01:06:49 20091
cm64-debugSuccessTue Mar 3 01:12:04 20093
x86_64-linux
cmSuccessSun Mar 6 00:01:01 20160
cm-debugSuccessSat Mar 5 00:01:12 20160

Testing status by file:


Html last generated: Sun Mar 6 05:50:07 2016

Input last modified: Wed Oct 2 10:36:46 2002


CMISS Help / Examples / 1 / 14 / 144