Deformation of the full prolate heart model
# CMISS (cm) command file for full ventricular mechanics.
#
# Prolate spheroidal coordinates are used to define the mesh consists
# using 120 elements with tricubicHermite interpolation for lambda and
# trilinear interpolation for mu, theta and the hydrostatic pressure.
# Fibre and sheet angles are interpolated with bilinear-cubicHermite
# and bicubicHermite-linear bases, resp. (imbrication angles are all
# zero). Equations of finite deformation elasticity govern the
# boundary value problem and the material response is described by the
# pole-zero constitutive law, with prescribed residual strains from
# Omens et al (??). Active fibre tension is modelled using the
# Hunter-McCulloch-ter Keurs equations for active myocyte contraction
# (n.b. zero for passive inflation). LV and RV cavities (endocardial
# surfaces) are inflated to 1.0 and 0.2 kPa, resp., while lambda and
# its in-wall-plane derivs are fixed for all epicardial nodes.
#
# Created: Martyn Nash October 2000
#
set echo;
#==========================================================================================================
print "\033[0;30;42m ============================================================ \033[0m\n";
print "\033[0;30;42m Execution Flow \033[0m\n";
print "\033[0;30;42m ============================================================ \033[0m\n";
#==========================================================================================================
set echo;
#
set num_threads 12;
#
fem define parameters;r;full_big;
#
$OUTPUT_DIR = "output/";
if( ! -d ${OUTPUT_DIR})
{
mkdir ${OUTPUT_DIR};
}
#
$TRUE = 1;
$FALSE = 0;
#
# Solve for or read in deformations
#
$MAX_PRESS = 10;
$SOLVE_INFLATION = $TRUE;
$SOLVE_ISOVOLUMIC_CONTRACTION = $TRUE;
$SOLVE_EJECTION = $TRUE;
$SOLVE_ISOVOLUMIC_RELAXATION = $FALSE;
#
if ( $SOLVE_ISOVOLUMIC_CONTRACTION || $SOLVE_EJECTION || $SOLVE_ISOVOLUMIC_RELAXATION )
{
$COUPLED_PROBLEM = $TRUE;
} else {
$COUPLED_PROBLEM = $FALSE;
}
set echo;
#==========================================================================================================
print "\033[0;30;42m ============================================================ \033[0m\n";
print "\033[0;30;42m Regions & Classes \033[0m\n";
print "\033[0;30;42m ============================================================ \033[0m\n";
#==========================================================================================================
set echo;
#
$WALL = 1;
$LV_CAVITY = 2;
$RV_CAVITY = 3;
set echo;
#==========================================================================================================
print "\033[0;30;42m ============================================================ \033[0m\n";
print "\033[0;30;42m Geometry \033[0m\n";
print "\033[0;30;42m ============================================================ \033[0m\n";
#==========================================================================================================
set echo;
#
fem define coord;r;fullheart_refmu region $WALL;
fem define node;r region $WALL;
fem define base;r;fullheart;
fem define element;r;fullheart_refmu region $WALL;
fem define fibre;r;fullheart_refmu region $WALL;
fem define element;r;fullheart_refmu fibre region $WALL;
#
#
# Define element groups
#
fem group elements 63..68 as base_lvfree_epi region $WALL;
fem group elements 93..98 as base_lvfree_endo region $WALL;
fem group elements 91,92,99,100 as base_septal region $WALL;
fem group elements 61,62,69,70 as base_rv region $WALL;
fem group elements base_lvfree_epi,base_lvfree_endo,base_septal,base_rv as base region $WALL;
#
fem group elements 3..8 as equatora_lvfree_epi region $WALL;
fem group elements 33..38 as equatora_lvfree_endo region $WALL;
fem group elements 31,32,39,40 as equatora_septal region $WALL;
fem group elements 1,2,9,10 as equatora_rv region $WALL;
fem group elements equatora_lvfree_epi,equatora_lvfree_endo,equatora_septal,equatora_rv as equatora region $WALL;
#
fem group elements 73..78 as equatorb_lvfree_epi region $WALL;
fem group elements 103..108 as equatorb_lvfree_endo region $WALL;
fem group elements 101,102,109,110 as equatorb_septal region $WALL;
fem group elements 71,72,79,80 as equatorb_rv region $WALL;
fem group elements equatorb_lvfree_epi,equatorb_lvfree_endo as equatorb_lvfree region $WALL;
fem group elements equatorb_lvfree,equatorb_septal,equatorb_rv as equatorb region $WALL;
#
fem group elements 13..18 as equatorc_lvfree_epi region $WALL;
fem group elements 43..48 as equatorc_lvfree_endo region $WALL;
fem group elements 41,42,49,50 as equatorc_septal region $WALL;
fem group elements 11,12,19,20 as equatorc_rv region $WALL;
fem group elements equatorc_lvfree_epi,equatorc_lvfree_endo as equatorc_lvfree region $WALL;
fem group elements equatorc_lvfree,equatorc_septal,equatorc_rv as equatorc region $WALL;
#
fem group elements equatora_lvfree_epi,equatorb_lvfree_epi,equatorc_lvfree_epi as equator_lvfree_epi region $WALL;
fem group elements equatora_lvfree_endo,equatorb_lvfree_endo,equatorc_lvfree_endo as equator_lvfree_endo region $WALL;
fem group elements equatora_septal,equatorb_septal,equatorc_septal as equator_septal region $WALL;
fem group elements equatora_rv,equatorb_rv,equatorc_rv as equator_rv region $WALL;
fem group elements equator_lvfree_epi,equator_lvfree_endo,equator_septal,equator_rv as equator region $WALL;
#
fem group elements 81..90 as apexa_epi region $WALL;
fem group elements 111..120 as apexa_endo region $WALL;
fem group elements apexa_epi,apexa_endo as apexa region $WALL;
#
fem group elements 21..30 as apexb_epi region $WALL;
fem group elements 51..60 as apexb_endo region $WALL;
fem group elements apexb_epi,apexb_endo as apexb region $WALL;
#
fem group elements apexa_epi,apexb_epi as apex_epi_elems region $WALL;
fem group elements apexa_endo,apexb_endo as apex_endo_elems region $WALL;
fem group elements apex_endo_elems,apex_epi_elems as apex_elems region $WALL;
#
fem group elements base_lvfree_endo,equator_lvfree_endo as lvfree_endo_elems region $WALL;
fem group elements base_lvfree_epi,equator_lvfree_epi as lvfree_epi_elems region $WALL;
fem group elements base_rv,equator_rv as rv_elems region $WALL;
fem group elements base_septal,equator_septal as septal_elems region $WALL;
#
#
# Define node groups
#
fem group nodes xi2=1 xi3=1 external elements base_lvfree_epi,base_rv as base_epi_nodes region $WALL;
fem group nodes xi2=1 xi3=0 external elements base_lvfree_endo,base_septal as base_lv_endo_nodes region $WALL;
fem group nodes xi2=1 xi3=0 external elements base_rv as base_rvfree_endo_nodes region $WALL;
fem group nodes xi2=1 xi3=1 external elements base_septal as base_rvsept_endo_nodes region $WALL;
fem group nodes base_lv_endo_nodes,base_rvfree_endo_nodes,base_rvsept_endo_nodes as base_endo_nodes region $WALL;
fem group nodes xi2=1 external elements base as base_nodes region $WALL;
#
# LV endo node groups
fem group nodes xi2=1 xi3=0 external elements equatora_lvfree_endo,equatora_septal as eqa_lv_endo_nodes region $WALL;
fem group nodes xi2=1 xi3=0 external elements equatorb_lvfree_endo,equatorb_septal as eqb_lv_endo_nodes region $WALL;
fem group nodes xi2=1 xi3=0 external elements equatorc_lvfree_endo,equatorc_septal as eqc_lv_endo_nodes region $WALL;
fem group nodes xi2=1 xi3=0 external elements apexa_endo as apexa_lv_endo_nodes region $WALL;
fem group nodes xi2=1 xi3=0 external elements apexb_endo as apexb_lv_endo_nodes region $WALL;
fem group nodes eqa_lv_endo_nodes,eqb_lv_endo_nodes,eqc_lv_endo_nodes,apexa_lv_endo_nodes as equator_lv_endo_nodes region $WALL;
#
# Epi node groups
fem group nodes xi2=1 xi3=1 external elements equatora_lvfree_epi,equatora_rv as eqa_epi_nodes region $WALL;
fem group nodes xi2=1 xi3=1 external elements equatorb_lvfree_epi,equatorb_rv as eqb_epi_nodes region $WALL;
fem group nodes xi2=1 xi3=1 external elements equatorc_lvfree_epi,equatorc_rv as eqc_epi_nodes region $WALL;
fem group nodes xi2=1 xi3=1 external elements apexa_epi as apexa_epi_nodes region $WALL;
fem group nodes xi2=1 xi3=1 external elements apexb_epi as apexb_epi_nodes region $WALL;
fem group nodes eqa_epi_nodes,eqb_epi_nodes,eqc_epi_nodes,apexa_epi_nodes as equator_epi_nodes region $WALL;
#
fem group nodes xi3=0 external elements equator_rv as equator_rvfree_endo_nodes region $WALL;
fem group nodes xi3=1 external elements equator_septal as equator_rvsept_endo_nodes region $WALL;
fem group nodes equator_rvfree_endo_nodes,equator_rvsept_endo_nodes as equator_rv_endo_nodes region $WALL;
fem group nodes xi3=1 elements equator_lvfree_endo,apex_endo_elems as equator_lv_mid_nodes region $WALL;
fem group nodes xi3=1 external elements lvfree_epi_elems,rv_elems,apex_epi_elems as epi_nodes region $WALL;
#
fem group nodes xi2=0 elements apexb as apex_nodes region $WALL;
fem group nodes xi2=0 xi3=0 external elements apexa_endo as near_apex_endo_nodes region $WALL;
fem group nodes xi2=0 xi3=1 external elements apexa_epi as near_apex_epi_nodes region $WALL;
fem group nodes equator_lv_mid_nodes,equator_rv_endo_nodes as equator_mid_nodes_temp region $WALL;
fem group nodes equator_mid_nodes_temp,near_apex_endo_nodes,near_apex_epi_nodes as equator_mid_nodes region $WALL;
#
set echo;
#==========================================================================================================
print "\033[0;30;42m ============================================================ \033[0m\n";
print "\033[0;30;42m Dependent variable/material information \033[0m\n";
print "\033[0;30;42m ============================================================ \033[0m\n";
#==========================================================================================================
set echo;
#
fem define equation;r;finelas_tch_incomp region $WALL lock;
fem define material;r;orth_incomp_active_residstrain region $WALL;
fem define active;r;active0_00 region $WALL;
fem define initial;r;fullheart_inflate_pericard region $WALL;
fem define solve;r;newton region $WALL;
set echo;
#==========================================================================================================
print "\033[0;30;42m ============================================================ \033[0m\n";
print "\033[0;30;42m Inflation \033[0m\n";
print "\033[0;30;42m ============================================================ \033[0m\n";
#==========================================================================================================
if ( $SOLVE_INFLATION )
{
print "\033[0;30;43m Increase cavity pressures incrementally to simulate diastole \033[0m\n";
#
set output;${OUTPUT_DIR}fullheart_inflate on;
$NAME1 = "fullheart_def_00";
$NAME2 = "_active_0000";
$ITERS = 50;
$ERROR_TOLERANCE = 0.001;
#
for $PRESS ( 0..$MAX_PRESS )
{
print "\033[0;30;43m Pressure = $PRESS \033[0m\n";
$NUMZEROS=2-length($PRESS);
$PRESSNUM="";
for $Z ( 1..$NUMZEROS )
{
$PRESSNUM=$PRESSNUM."0";
}
$PRESSNUM=$PRESSNUM.$PRESS;
if( $PRESS == 0 )
{
$INCREM = 0.0;
} else
{
$INCREM = 0.1;
}
fem solve increment $INCREM iterate $ITERS error $ERROR_TOLERANCE;
$FILENAME = ${NAME1}.${PRESSNUM}.${NAME2};
fem define initial;w;${OUTPUT_DIR}.${FILENAME} region $WALL;
fem export nodes;${OUTPUT_DIR}.${FILENAME} as heart;
fem export nodes;${OUTPUT_DIR}${FILENAME}_def field as heart;
fem export elements;${OUTPUT_DIR}.${FILENAME} as wall;
fem export elements;${OUTPUT_DIR}${FILENAME}_def field as wall;
}
set output off;
} else {
print "\033[0;30;43m Reading in end-diastolic state \033[0m\n";
#
$FILENAME = "fullheart_def_00".${MAX_PRESS}."_active_0000";
fem define initial;r;${OUTPUT_DIR}.${FILENAME} region $WALL;
fem define solve;r;newton region $WALL;
fem export nodes;${OUTPUT_DIR}.${FILENAME} as heart;
fem export nodes;${OUTPUT_DIR}${FILENAME}_def field as heart;
fem export elements;${OUTPUT_DIR}.${FILENAME} as wall;
fem export elements;${OUTPUT_DIR}${FILENAME}_def field as wall;
}
set echo;
if ( $COUPLED_PROBLEM ) { # Set up coupled regions
set echo;
#==========================================================================================================
print "\033[0;30;42m ============================================================ \033[0m\n";
print "\033[0;30;42m Set up Coupled Regions \033[0m\n";
print "\033[0;30;42m ============================================================ \033[0m\n";
#==========================================================================================================
#
# Save scale factors for the ventricular wall elements
#
fem define base;r;fullheart_readse; #redefine bases to read in scale factors
fem def line;w;fullheart_refmu region $WALL; #write out scale factors for wall elements
fem define base;r;fullheart; #redefine bases to calculate scale factors
#
# Define geometry and dependent variable information for LV and RV cavities
#
fem define region;r;coupled;
fem define coord;r;lvcavity_refmu region $LV_CAVITY;
fem define node;r region $LV_CAVITY;
fem define coord;r;rvcavity_refmu region $RV_CAVITY;
fem define node;r region $RV_CAVITY;
fem define element;r;lvcavity_refmu region $LV_CAVITY;
fem define element;r;rvcavity_refmu region $RV_CAVITY;
fem group elem 121..180 as lv_cavity_elems region $LV_CAVITY;
fem group elem 151..160 as lv_cavity_base_elems region $LV_CAVITY;
fem group nodes xi2=1 xi3=1 elem lv_cavity_base_elems as lv_cavity_base_nodes region $LV_CAVITY;
fem group elem 181..198 as rv_cavity_elems region $RV_CAVITY;
fem define equation;r;coupled region $WALL,$LV_CAVITY,$RV_CAVITY lock;
fem define material;r;cavity region $LV_CAVITY;
fem define material;r;cavity region $RV_CAVITY;
fem define initial;r;lvcavity_refmu region $LV_CAVITY;
fem define initial;r;rvcavity_refmu region $RV_CAVITY;
#
# Define cavity/wall coupling and solution information
#
fem define coupling;r;coupled;
fem define solve;r;coupled coupled region $WALL,$LV_CAVITY,$RV_CAVITY;
fem update solution coupled source_region $WALL;
$NHlambda = 1;
fem update solution cavity_reference average 196 in $NHlambda node lv_cavity_base_nodes region $LV_CAVITY;
fem update solution cavity_reference average 203 in $NHlambda node 14,18 region $RV_CAVITY;
#
# Read back in scale factors for the ventricular wall elements
#
fem define base;r;fullheart_readse; #redefine bases to read in scale factors
fem def line;r;fullheart_refmu region $WALL; #read in scale factors for wall elements
fem define base;r;fullheart; #redefine bases to calculate scale factors
fem export nodes;${OUTPUT_DIR}heart as heart region all;
fem export elements;${OUTPUT_DIR}lv_cavity as lv_cavity region $LV_CAVITY;
fem export elements;${OUTPUT_DIR}rv_cavity as rv_cavity region $RV_CAVITY;
fem export elements;${OUTPUT_DIR}wall as wall region $WALL;
fem export nodes;${OUTPUT_DIR}heart_def field as heart region all;
fem export elements;${OUTPUT_DIR}lv_cavity_def field as lv_cavity region $LV_CAVITY;
fem export elements;${OUTPUT_DIR}rv_cavity_def field as rv_cavity region $RV_CAVITY;
fem export elements;${OUTPUT_DIR}wall_def field as wall region $WALL;
set echo;
}
if ( $SOLVE_ISOVOLUMIC_CONTRACTION ) {
set echo;
#==========================================================================================================
print "\033[0;30;42m ============================================================ \033[0m\n";
print "\033[0;30;42m Solving Isovolumic Contraction \033[0m\n";
print "\033[0;30;42m ============================================================ \033[0m\n";
#==========================================================================================================
set output;${OUTPUT_DIR}steps_contract_infl0_${MAX_PRESS} on;
$NAME1 = "active0_";
$NAME2 = "fullheart_def_00".${MAX_PRESS}."_active_00";
$INCREM = 0.0;
$ITER = 20;
$ERROR_TOLERANCE = 0.001;
#
# Isovolumically contract the ventricles by incrementally
# increasing intracellular calcium concentration and solving.
if ( $MAX_PRESS == 10 ) {
$MAX_CA = 24;
} elsif ( $MAX_PRESS == 30 ) {
$MAX_CA = 30;
} else {
fem quit # required activation level unknown
}
for $CA ( 1..$MAX_CA )
{
$NUMZEROS=2-length($CA);
$CANUM="";
for $Z ( 1..$NUMZEROS )
{
$CANUM=$CANUM."0";
}
$CANUM=$CANUM.$CA
print "\033[0;30;43m Calcium Level = ${CA} \033[0m\n";
$FILENAME=${NAME2}.${CANUM}
fem define active;r;${NAME1}.${CANUM} region $WALL;
fem solve increment $INCREM iter $ITER coupled error $ERROR_TOLERANCE;
fem define initial;w;${OUTPUT_DIR}.${FILENAME} region $WALL,$LV_CAVITY,$RV_CAVITY;
fem export nodes;${OUTPUT_DIR}.${FILENAME} as heart region all;
fem export nodes;${OUTPUT_DIR}${FILENAME}_def field as heart region all;
}
set out off
set echo;
}
if ( $SOLVE_EJECTION ) {
set echo;
#==========================================================================================================
print "\033[0;30;42m ============================================================ \033[0m\n";
print "\033[0;30;42m Solving Ejection \033[0m\n";
print "\033[0;30;42m ============================================================ \033[0m\n";
#==========================================================================================================
if ( !$SOLVE_ISOVOLUMIC_CONTRACTION ) {
print "\033[0;30;43m Reading in state at the end of isovolumic contraction \033[0m\n";
$CA = 24;
fem define active;r;"active0_".${CA} region $WALL;
fem define initial;r;${OUTPUT_DIR}."fullheart_def_00".${MAX_PRESS}."_active_00".${CA} nofix region $WALL,$LV_CAVITY,$RV_CAVITY;
}
$CA = 24;
set out;${OUTPUT_DIR}."steps_eject_infl0_".${MAX_PRESS}."_act0_".${CA} on;
$NAME1 = "cavity_";
$NAME2 = "fullheart_def_00".${MAX_PRESS}."_active_00".${CA}."_cavk_";
$INCREM = 0.0;
$ITERS = 20;
$ERROR_TOLERANCE = 0.001;
#
# Simulate ventricular ejection by incrementally decreasing the
# cavity impedances.
$CAVITY_K =640;
for $I ( 1..26 )
{
$CAVITY_K2 = int($CAVITY_K*1000);
$NUMZEROS=8-length($CAVITY_K2)
$CAVITY_K_NUMBER="";
for $Z ( 1..$NUMZEROS )
{
$CAVITY_K_NUMBER=$CAVITY_K_NUMBER."0";
}
$CAVITY_K_NUMBER=$CAVITY_K_NUMBER.$CAVITY_K2;
$FILENAME=${NAME2}.${CAVITY_K_NUMBER};
print "\033[0;30;43m Cavity Stiffness Number ${CAVITY_K_NUMBER} \033[0m\n";
fem define material;r;${NAME1}.${CAVITY_K_NUMBER} region $LV_CAVITY;
fem define material;r;${NAME1}.${CAVITY_K_NUMBER} region $RV_CAVITY;
fem solve increment $INCREM iterate $ITERS coupled error $ERROR_TOLERANCE;
fem define initial;w;${OUTPUT_DIR}.${NAME2}.${CAVITY_K_NUMBER} region $WALL,$LV_CAVITY,$RV_CAVITY;
fem export nodes;${OUTPUT_DIR}.${FILENAME} as heart region all;
fem export nodes;${OUTPUT_DIR}${FILENAME}_def field as heart region all;
if( $CAVITY_K > 10 )
{
$CAVITY_K = $CAVITY_K / 2.0;
} elsif( $CAVITY_K > 4 )
{
$CAVITY_K = $CAVITY_K - 1.0;
} else
{
$CAVITY_K = $CAVITY_K - 0.3;
}
}
set out off;
set echo;
}
if ( $SOLVE_ISOVOLUMIC_RELAXATION ) {
set echo;
#==========================================================================================================
print "\033[0;30;42m ============================================================ \033[0m\n";
print "\033[0;30;42m Solving Isovolumic Relaxation \033[0m\n";
print "\033[0;30;42m ============================================================ \033[0m\n";
#==========================================================================================================
set echo;
}
Name Modified Size
example_8111.com 01-Feb-2001 20k active0_00.ipacti 03-Mar-2004 777 active0_01.ipacti 03-Mar-2004 777 active0_02.ipacti 03-Mar-2004 777 active0_03.ipacti 03-Mar-2004 777 active0_04.ipacti 03-Mar-2004 777 active0_05.ipacti 03-Mar-2004 777 active0_06.ipacti 03-Mar-2004 777 active0_07.ipacti 03-Mar-2004 777 active0_08.ipacti 03-Mar-2004 777 active0_09.ipacti 03-Mar-2004 777 active0_10.ipacti 03-Mar-2004 777 active0_11.ipacti 03-Mar-2004 777 active0_12.ipacti 03-Mar-2004 777 active0_13.ipacti 03-Mar-2004 777 active0_14.ipacti 03-Mar-2004 777 active0_15.ipacti 03-Mar-2004 777 active0_16.ipacti 03-Mar-2004 777 active0_17.ipacti 03-Mar-2004 777 active0_18.ipacti 03-Mar-2004 777 active0_19.ipacti 03-Mar-2004 777 active0_20.ipacti 03-Mar-2004 777 active0_21.ipacti 03-Mar-2004 777 active0_22.ipacti 03-Mar-2004 777 active0_23.ipacti 03-Mar-2004 777 active0_24.ipacti 03-Mar-2004 777 active0_25.ipacti 03-Mar-2004 777 active0_26.ipacti 03-Mar-2004 777 active0_27.ipacti 03-Mar-2004 777 active0_28.ipacti 03-Mar-2004 777 active0_29.ipacti 03-Mar-2004 777 active0_30.ipacti 03-Mar-2004 777 cavity.ipmate 01-Feb-2001 390 cavity_00000100.ipmate 01-Feb-2001 386 cavity_00000400.ipmate 01-Feb-2001 386 cavity_00000700.ipmate 01-Feb-2001 386 cavity_00001000.ipmate 01-Feb-2001 386 cavity_00001300.ipmate 01-Feb-2001 386 cavity_00001600.ipmate 01-Feb-2001 386 cavity_00001900.ipmate 01-Feb-2001 386 cavity_00002200.ipmate 01-Feb-2001 386 cavity_00002500.ipmate 01-Feb-2001 386 cavity_00002800.ipmate 01-Feb-2001 386 cavity_00003100.ipmate 01-Feb-2001 386 cavity_00003400.ipmate 01-Feb-2001 386 cavity_00003700.ipmate 01-Feb-2001 386 cavity_00004000.ipmate 01-Feb-2001 386 cavity_00005000.ipmate 01-Feb-2001 386 cavity_00006000.ipmate 01-Feb-2001 386 cavity_00007000.ipmate 01-Feb-2001 386 cavity_00008000.ipmate 01-Feb-2001 386 cavity_00009000.ipmate 01-Feb-2001 386 cavity_00010000.ipmate 01-Feb-2001 388 cavity_00020000.ipmate 01-Feb-2001 388 cavity_00040000.ipmate 01-Feb-2001 388 cavity_00080000.ipmate 01-Feb-2001 388 cavity_00160000.ipmate 01-Feb-2001 390 cavity_00320000.ipmate 01-Feb-2001 390 cavity_00640000.ipmate 01-Feb-2001 390 coupled.ipcoup 21-Aug-2002 465 coupled.ipcoup.old 01-Feb-2001 414 coupled.ipequa 02-May-2004 2.2k coupled.ipregi 01-Feb-2001 138 coupled.ipsolv 16-Aug-2010 3.2k coupled.ipsolv.old 13-Apr-2007 3.0k coupled.irequa 02-May-2004 5.6k coupled.irsolv 16-Aug-2010 3.3k coupled.irsolv.old 13-Apr-2007 3.1k finelas_tch_incomp.ipequa 02-May-2004 2.2k full_big.ippara 12-Nov-2002 5.8k fullheart.ipbase 01-Feb-2001 16k fullheart_inflate_pericard.ipinit 01-Feb-2001 2.8k fullheart_readse.ipbase 01-Feb-2001 16k fullheart_refmu.exelem 01-Feb-2001 741k fullheart_refmu.exnode 01-Feb-2001 105k fullheart_refmu.ipcoor 01-Feb-2001 758 fullheart_refmu.ipelem 01-Feb-2001 87k fullheart_refmu.ipelfb 01-Feb-2001 56k fullheart_refmu.ipfibr 01-Feb-2001 94k fullheart_refmu.ipline 01-Feb-2001 1.1M fullheart_refmu.ipnode 01-Feb-2001 157k lvcavity_refmu.ipcoor 01-Feb-2001 779 lvcavity_refmu.ipelem 01-Feb-2001 59k lvcavity_refmu.ipinit 01-Feb-2001 2.7k lvcavity_refmu.ipnode 01-Feb-2001 60k newton.ipsolv 16-Aug-2010 2.7k newton.ipsolv.old 13-Apr-2007 2.5k orth_incomp_active_residstrain.ipmate 24-May-2001 6.3k output/ 27-Feb-2004 - rvcavity_refmu.ipcoor 01-Feb-2001 779 rvcavity_refmu.ipelem 01-Feb-2001 12k rvcavity_refmu.ipinit 01-Feb-2001 1.5k rvcavity_refmu.ipnode 01-Feb-2001 30k
Name Modified Size
examples_8_81_811_8111.tar.gz 17-Aug-2010 10.0M
Html last generated: Sun Mar 6 05:50:32 2016
Input last modified: Mon Aug 16 11:22:37 2010