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.
#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
#Testing example 148 fem list elements
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
Name Modified Size
examples_1_14_148.tar.gz 18-Aug-2006 5.7M
| Status | Tested | Real time (s) | |
| i686-linux | |||
| cm | Success | Sun Mar 6 00:01:36 2016 | 1 |
| cm-debug | Success | Sat Mar 5 00:05:02 2016 | 2 |
| mips-irix | |||
| cm | Success | Sun Aug 19 01:43:47 2007 | 12 |
| cm-debug | Success | Wed Aug 15 01:31:14 2007 | 18 |
| cm-debug-clear-malloc | Success | Sat Aug 18 01:38:02 2007 | 21 |
| cm-debug-clear-malloc7 | Success | Mon Aug 20 01:34:37 2007 | 21 |
| cm64 | Success | Sun Aug 19 01:43:53 2007 | 12 |
| cm64-debug | Success | Tue Aug 21 01:30:19 2007 | 18 |
| cm64-debug-clear-malloc | Success | Thu Apr 1 10:10:53 2004 | 9 |
| rs6000-aix | |||
| cm | Success | Wed Mar 4 01:08:47 2009 | 2 |
| cm-debug | Success | Mon Mar 2 01:08:06 2009 | 3 |
| cm64 | Success | Wed Mar 4 01:07:59 2009 | 1 |
| cm64-debug | Success | Tue Mar 3 01:13:07 2009 | 3 |
| x86_64-linux | |||
| cm | Success | Sun Mar 6 00:01:01 2016 | 1 |
| cm-debug | Success | Sat Mar 5 00:01:39 2016 | 1 |
| 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:08 2016
Input last modified: Wed Oct 2 10:37:27 2002