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

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

This model uses the same bases and connectivity as example 144, however the nodal positions are the same those in the prolate model example 141 (except for node 15 and 28 which have been moved slightly). The all the surface derivatives are fitted to surface data from the prolate model, except the lv endocardial surface which has been fitted from the original h27 data set.

 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.

 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 comfile run by this example is as follows:

#Example_148: 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;new_rc_heart;example     #99 nodes
fem define base;r;heart;example    
fem define element;r;new_rc_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:

#Testing example 148
fem list elements

Files used by this example are:

Name                 Modified     Size

example_148.com 29-Dec-2000 810 heart.exelem 03-Jan-2001 510k heart.exnode 03-Jan-2001 116k heart.ipbase 03-Jan-2001 11k heart.ipelem 03-Jan-2001 64k heart.ipelfb 03-Jan-2001 27k heart.ipfibr 30-Jan-2001 49k new_rc_heart.ipelem 03-Jan-2001 64k new_rc_heart.ipnode 03-Jan-2001 268k rc3d.ipcoor 03-Jan-2001 728 test_output.com 29-Dec-2000 39 viewheart.com 09-Nov-2001 624 work_files.tar.gz 03-Jan-2001 5.7M

Download the entire example:

Name                      Modified     Size

examples_1_14_148.tar.gz 18-Aug-2006 5.7M

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:01:36 20161
cm-debugSuccessSat Mar 5 00:05:02 20162
mips-irix
cmSuccessSun Aug 19 01:43:47 200712
cm-debugSuccessWed Aug 15 01:31:14 200718
cm-debug-clear-mallocSuccessSat Aug 18 01:38:02 200721
cm-debug-clear-malloc7SuccessMon Aug 20 01:34:37 200721
cm64SuccessSun Aug 19 01:43:53 200712
cm64-debugSuccessTue Aug 21 01:30:19 200718
cm64-debug-clear-mallocSuccessThu Apr 1 10:10:53 20049
rs6000-aix
cmSuccessWed Mar 4 01:08:47 20092
cm-debugSuccessMon Mar 2 01:08:06 20093
cm64SuccessWed Mar 4 01:07:59 20091
cm64-debugSuccessTue Mar 3 01:13:07 20093
x86_64-linux
cmSuccessSun Mar 6 00:01:01 20161
cm-debugSuccessSat Mar 5 00:01:39 20161

Testing status by file:


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

Input last modified: Wed Oct 2 10:37:27 2002


CMISS Help / Examples / 1 / 14 / 148