Example 811b: Porcine ventricular model

The complete model of a pig heart. Geometry and fibres from specimen MRP01, sheets from MRP02.


The comfile run by this example is as follows:

set echo;
#==========================================================================================================
print "\033[0;30;42m ==================================================================== \033[0m\n"; 
print "\033[0;30;42m         Nodes, Elements, Fibres & Sheets for Pig heart Model         \033[0m\n"; 
print "\033[0;30;42m ==================================================================== \033[0m\n"; 
#==========================================================================================================
set echo;
#
fem define parameters;r;MRP01;example;
fem define coordinates;r;MRP01;example;
fem define region;r;MRP01;example;
#
fem define node;r;MRP01;example;
fem define nodes;r;MRP01_base;example region 2;
#
fem define base;r;MRP01;example;                                 # 1..8 Geometry and fibres
fem define;add bases;r;3Linear-2PressAuxXi3_4x4x4Gauss;example;  # 9 Pressure
#
fem define elem;r;MRP01;example;
fem define elements;r;MRP01_base;example region 2;
#
fem define base;r;MRP01_readse;example;                          #redefine bases to read in scale factors 
fem define;add bases;r;3Linear-2PressAuxXi3_4x4x4Gauss;example; 
fem define line;r;MRP01;example;
fem define line;r;MRP01_base;example region 2;
#
fem update nodes hanging nodes 144,146 element 47 direction 1;
#
fem define fibre;r;MRP01;example;
fem define elem;r;MRP01;example fibre;
#
#
#
fem group nodes 5..8,11..12,15..16,19..20,22,24,27..28,33..34,36,38,40,66..67,69..70 as EPI_NODES1;
fem group nodes 72..73,75..76,78,80,83..84,86,88,97..98,101..102 as EPI_NODES2;
fem group nodes 104,106,108,110,112,114,116,118..119,128..136,146..147,30 as EPI_NODES3;
fem group nodes EPI_NODES1,EPI_NODES2,EPI_NODES3 as EPI_NODES;
#
fem group nodes 1,3,9,13,17,20..21,23,25..26,29,31..32,35,37,39,41..42,45..46,49..50,53..54,57 as LV_NODES1;
fem group nodes 59,61,63,81..82,85,87,89,91,93,95..96,103,105,107,120,122,124..127 as LV_NODES2;
fem group nodes LV_NODES1,LV_NODES2 as LV_NODES;
#
fem group nodes 2,4,10,14,43,47,51,55,65,68,71,74,77,79,92,94,99..100,109,111,113,115,117,123 as RV_FREE_NODES;
fem group nodes 2,4,10,14,43..44,47..48,51..52,55..56,58,60,62,64,90,92,94,99..100,109,121,123 as RV_SEPT_NODES;
#
fem group faces 6,21,43,183,210,214,218,239,244,248,252,323 as EPI_FACES1;
fem group faces 11,25,47,120,124,133,187,195,199,203,206 as EPI_FACES2;
fem group faces 16,29,51,101,106,111,116,129,158,163,172 as EPI_FACES3;
fem group faces 34,55,137,179,277,296,299,302,305,308,310 as EPI_FACES4;
fem group faces 38,58,140,176,261,274,312,314,316,318..319 as EPI_FACES5;
fem group faces EPI_FACES1,EPI_FACES2,EPI_FACES3,EPI_FACES4,EPI_FACES5 as EPI_FACES;
#
fem group faces 20,24,28,33,37,42,46,50,54,57,63,68,73,78,83,87,91,95,128,132,136,139,144,148 as LV_FACES1;
fem group faces 152,157,167,175,178,182,186,190,221,225,230,234,247,255,265,268,280,283,286,289,292 as LV_FACES2;
fem group faces LV_FACES1,LV_FACES2 as LV_FACES;
#
fem group faces 100,105,110,115,119,123,162,198,202,205,209,213,217,238,243,322 as RV_FREE_FACES;
fem group faces 64,69,74,79,84,88,92,96,145,149,153,168,191,222,226 as RV_SEPT_FACES;
fem group faces RV_FREE_FACES,RV_SEPT_FACES as RV_FACES;
#
#
#
fem export nodes;MRP01 as MRP01 region all;
fem export elements;MRP01 as MRP01 region 1;
fem export elements;MRP01_base as base_skel region 2;


Additional testing commands:

#Testing example 811b
fem list nodes region all;
fem list elements region all;
#
fem export nodes;MRP01 as MRP01 region all;
fem export elements;MRP01 as MRP01 region 1;
fem export elements;MRP01_base as base_skel region 2;

Files used by this example are:

Name                                    Modified     Size

example_811b.com 05-Feb-2002 3.3k 3Linear-2PressAuxXi3_4x4x4Gauss.ipbase 05-Feb-2002 1.9k MRP01.exelem 10-Sep-2002 528k MRP01.exnode 10-Sep-2002 340k MRP01.ipbase 10-Sep-2002 11k MRP01.ipcoor 03-Feb-2002 688 MRP01.ipelem 11-Sep-2002 102k MRP01.ipelfb 11-Sep-2002 68k MRP01.ipfibr 11-Sep-2002 415k MRP01.ipline 10-Sep-2002 455k MRP01.ipnode 11-Sep-2002 469k MRP01.ippara 12-Nov-2002 5.9k MRP01.ipregi 03-Feb-2002 142 MRP01_base.exelem 10-Sep-2002 42k MRP01_base.ipelem 05-Feb-2002 12k MRP01_base.ipline 03-Feb-2002 36k MRP01_base.ipnode 03-Feb-2002 121k MRP01_readse.ipbase 10-Sep-2002 11k test_output.com 03-Feb-2002 224 view.com 03-Feb-2002 6.0k

Download the entire example:

Name                           Modified     Size

examples_8_81_811_811b.tar.gz 20-Aug-2014 558k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmFailureSun Mar 6 00:40:39 201623
last breakMon Aug 18 01:29:00 201470
cm-debugFailureSun Mar 6 01:24:17 201678
last breakMon Aug 18 02:26:00 2014128
mips-irix
cmFailureTue Aug 21 01:27:30 20074
last breakTue Apr 11 01:21:00 20064
last successSun May 18 02:27:00 2003203
cm-debugFailureTue Aug 21 01:27:39 20078
last breakTue Apr 11 01:24:00 20067
last successWed May 21 04:11:00 2003712
cm-debug-clear-mallocFailureTue Aug 21 01:58:26 200717
last breakTue Apr 11 01:48:00 200617
cm-debug-clear-malloc7FailureTue Aug 21 02:46:05 200717
last breakTue Apr 11 01:47:00 200616
cm64SuccessSun Aug 19 04:21:14 2007203
cm64-debugSuccessTue Aug 21 03:59:11 2007755
cm64-debug-clear-mallocSuccessThu Apr 1 10:51:50 2004335
rs6000-aix
cmFailureWed Mar 4 01:29:41 200926
last breakThu Feb 7 01:20:00 200826
cm-debugFailureWed Mar 4 02:16:33 2009243
last breakThu Feb 7 01:22:00 2008244
cm64FailureWed Mar 4 01:30:09 200928
last breakThu Feb 7 01:20:00 200828
cm64-debugFailureWed Mar 4 02:20:16 2009278
last breakThu Feb 7 01:31:00 2008272
x86_64-linux
cmFailureSun Mar 6 00:02:19 20168
last breakWed Aug 20 00:12:00 201415
cm-debugFailureSun Mar 6 00:05:23 201637
last breakWed Aug 20 00:38:00 201448

Testing status by file:


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

Input last modified: Tue Aug 19 10:42:53 2014


CMISS Help / Examples / 8 / 81 / 811 / 811b