Example 5g4: LV inflation using incompressible pole-zero material law

An elliptically shaped mesh in rectangular Cartesian coordinates is inflated.

This example is the same problem as example 5g3, but instead of using cellML, it uses the back-end version of the pole-zero law. The point of this example is thus to produce the same outputs as in example 5g3 so it is possible to verify that they produce the same results.

Since this example is not using cellML it may not belong under the 5g examples, but I leave that for others to decide and change. If this example is moved, make sure to update this description and the reference to this example in example 5g3.

Created by Espen Remme May 2004.
 


The comfile run by this example is as follows:

# Example 5g4

# If the example path is not set, default to current directory
if(!defined $example)
{
  $example = "./";
}

fem def param;r;full_big;example
fem define coor;r;rcCoord;example
fem define base;r;basisfuncs;example
#
fem define nodes;r;rc3;example
fem define elements;r;rc;example
fem define fibre;r;simple;example
fem define element;r;simple;example fibre
  
# Material parameter groups
fem group elem 1..4 as apexElems
fem group elem 5..12 as betweenElems
fem group elem 13..16 as baseElems
fem group elem 5..16 as notApexElems

# Boundary Condition groups
fem group nodes 1,18 as fixedApexNodes
fem group nodes 14..17,31..34 as fixedBaseNodes
fem group nodes 14,16 as y_axisNodes
fem group nodes 15,17 as z_axisNodes

# Residual strains groups
fem group nodes 6..13 as nodes095
fem group nodes 1..5,14..22,31..34 as nodes100
fem group nodes 23..30 as nodes105

$name = "heartOriginal"
$file = $name
fem export nodes;$file as $name offset 1000
fem export elements;$file as $name offset_elem 1000

fem define equation;r;equation;example lock
fem define material;r;polezero_gauss;example
fem define acti;r;calcium000;example

fem define initial;r;x_Base_node_restricted;example
fem define solve;r;solve;example

fem list elements total

fem solve increment 1 iter 30 error 10e-8

fem list elements deformed total

$name = "heartInflated"
$file = $name
fem export nodes;$file field as $name offset 0000
fem export elements;$file field as $name offset_elem 0000

fem list strain;strain fibre
fem list stress;stress fibre total 

$file = "heartStress"
fem update gauss stress fibre total Cauchy
fem export gauss;$file yg as gauss_stress offset 0

$file = "heartStrain"
fem update gauss strain fibre extension_ratios
fem export gauss;$file yg as gauss_strain

fem update geometry from solution
$file = "heartInflated"
fem define nodes;w;$file

Files used by this example are:

Name                           Modified     Size

example_5g4.com 03-May-2004 1.8k basisfuncs.ipbase 03-May-2004 16k calcium000.ipacti 03-May-2004 758 equation.ipequa 03-May-2004 2.1k full_big.ippara 03-May-2004 5.8k gaussDraw.com 03-May-2004 1.3k polezero_gauss.ipmate 03-May-2004 17k rc.ipelem 03-May-2004 15k rc3.ipnode 03-May-2004 75k rcCoord.ipcoor 03-May-2004 678 simple.ipelfb 03-May-2004 4.2k simple.ipfibr 03-May-2004 6.6k solve.ipsolv 16-Aug-2010 2.2k solve.ipsolv.old 13-Apr-2007 2.1k x_Base_node_restricted.ipinit 03-May-2004 9.8k

Download the entire example:

Name                      Modified     Size

examples_5_5g_5g4.tar.gz 14-Aug-2014 64k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:18:30 201634
cm-debugSuccessSat Mar 5 01:55:45 2016191
mips-irix
cmSuccessSun Aug 19 06:14:50 2007292
cm-debugSuccessWed Aug 8 05:18:01 20071613
cm-debug-clear-mallocFailureSat Aug 4 06:33:51 20071615
last breakSat Aug 4 06:06:00 20071615
last successFri Jun 29 04:23:00 20071783
cm-debug-clear-malloc7SuccessMon Aug 20 05:32:44 20071844
cm64SuccessThu Aug 16 03:55:02 2007309
cm64-debugSuccessFri Aug 17 04:08:47 20071753
cm64-debug-clear-mallocSuccessTue Feb 1 09:54:56 2005719
rs6000-aix
cmSuccessWed Mar 4 01:22:33 200937
cm-debugSuccessMon Mar 2 02:10:03 2009580
cm64SuccessWed Mar 4 01:22:37 200939
cm64-debugSuccessTue Mar 3 02:15:47 2009562
x86_64-linux
cmSuccessSun Mar 6 00:02:09 201619
cm-debugSuccessSat Mar 5 00:07:46 201686

Testing status by file:


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

Input last modified: Thu Aug 14 10:03:36 2014


CMISS Help / Examples / 5 / 5g / 5g4