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.
#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
fem list element total
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
Name Modified Size
examples_1_14_144.tar.gz 18-Aug-2006 115k
| Status | Tested | Real time (s) | |
| i686-linux | |||
| cm | Success | Sun Mar 6 00:01:35 2016 | 1 |
| cm-debug | Success | Sat Mar 5 00:01:54 2016 | 1 |
| mips-irix | |||
| cm | Success | Sun Aug 19 01:32:54 2007 | 6 |
| cm-debug | Success | Wed Aug 15 01:27:29 2007 | 10 |
| cm-debug-clear-malloc | Success | Sat Aug 18 01:31:43 2007 | 15 |
| cm-debug-clear-malloc7 | Success | Mon Aug 20 01:28:58 2007 | 14 |
| cm64 | Success | Sun Aug 19 01:32:56 2007 | 6 |
| cm64-debug | Success | Tue Aug 21 01:26:49 2007 | 11 |
| cm64-debug-clear-malloc | Success | Thu Apr 1 10:10:27 2004 | 5 |
| rs6000-aix | |||
| cm | Success | Wed Mar 4 01:08:41 2009 | 1 |
| cm-debug | Success | Mon Mar 2 01:07:18 2009 | 2 |
| cm64 | Success | Wed Mar 4 01:06:49 2009 | 1 |
| cm64-debug | Success | Tue Mar 3 01:12:04 2009 | 3 |
| x86_64-linux | |||
| cm | Success | Sun Mar 6 00:01:01 2016 | 0 |
| cm-debug | Success | Sat Mar 5 00:01:12 2016 | 0 |
| i686-linux | |||
| Success | cm: | cmiss_test.log.retain. | |
| Success | cm-debug: | cmiss_test.log.retain. | |
| mips-irix | |||
| Success | cm: | cmiss_test.log.retain. | |
| Success | cm-debug: | cmiss_test.log.retain. | |
| Success | cm-debug-clear-malloc: | cmiss_test.log.retain. | |
| Success | cm-debug-clear-malloc7: | cmiss_test.log.retain. | |
| Success | cm64: | cmiss_test.log.retain. | |
| Success | cm64-debug: | cmiss_test.log.retain. | |
| Success | cm64-debug-clear-malloc: | cmiss_test.log.retain. | |
| rs6000-aix | |||
| Success | cm: | cmiss_test.log.retain. | |
| Success | cm-debug: | cmiss_test.log.retain. | |
| Success | cm64: | cmiss_test.log.retain. | |
| Success | cm64-debug: | cmiss_test.log.retain. | |
| x86_64-linux | |||
| Success | cm: | cmiss_test.log.retain. | |
| Success | cm-debug: | cmiss_test.log.retain. |
| i686-linux | |||
| Success | cm: | ndiff test: no significant differences with generic answer. | |
| Success | cm-debug: | ndiff test: no significant differences with generic answer. | |
| mips-irix | |||
| Success | cm: | ndiff test: no significant differences with generic answer. | |
| Success | cm-debug: | ndiff test: no significant differences with generic answer. | |
| Success | cm-debug-clear-malloc: | ndiff test: no significant differences with generic answer. | |
| Success | cm-debug-clear-malloc7: | ndiff test: no significant differences with generic answer. | |
| Success | cm64: | ndiff test: no significant differences with generic answer. | |
| Success | cm64-debug: | ndiff test: no significant differences with generic answer. | |
| Success | cm64-debug-clear-malloc: | ndiff test: no significant differences with generic answer. | |
| rs6000-aix | |||
| Success | cm: | ndiff test: no significant differences with generic answer. | |
| Success | cm-debug: | ndiff test: no significant differences with generic answer. | |
| Success | cm64: | ndiff test: no significant differences with generic answer. | |
| Success | cm64-debug: | ndiff test: no significant differences with generic answer. | |
| x86_64-linux | |||
| Success | cm: | ndiff test: no significant differences with generic answer. | |
| Success | cm-debug: | ndiff test: no significant differences with generic answer. |
Html last generated: Sun Mar 6 05:50:07 2016
Input last modified: Wed Oct 2 10:36:46 2002