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.