LV inflation using incompressible pole-zero material law

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:

	Status	Tested	Real time (s)
i686-linux
cm	Success	Sun Mar 6 00:18:30 2016	34
cm-debug	Success	Sat Mar 5 01:55:45 2016	191
mips-irix
cm	Success	Sun Aug 19 06:14:50 2007	292
cm-debug	Success	Wed Aug 8 05:18:01 2007	1613
cm-debug-clear-malloc	Failure	Sat Aug 4 06:33:51 2007	1615
	last break	Sat Aug 4 06:06:00 2007	1615
	last success	Fri Jun 29 04:23:00 2007	1783
cm-debug-clear-malloc7	Success	Mon Aug 20 05:32:44 2007	1844
cm64	Success	Thu Aug 16 03:55:02 2007	309
cm64-debug	Success	Fri Aug 17 04:08:47 2007	1753
cm64-debug-clear-malloc	Success	Tue Feb 1 09:54:56 2005	719
rs6000-aix
cm	Success	Wed Mar 4 01:22:33 2009	37
cm-debug	Success	Mon Mar 2 02:10:03 2009	580
cm64	Success	Wed Mar 4 01:22:37 2009	39
cm64-debug	Success	Tue Mar 3 02:15:47 2009	562
x86_64-linux
cm	Success	Sun Mar 6 00:02:09 2016	19
cm-debug	Success	Sat Mar 5 00:07:46 2016	86

Testing status by file:

signal

mips-irix

SIGKILL cm-debug-clear-malloc: exit due to SIGKILL signal.

cmiss_test.log

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.

heartInflated.ipnode

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.
Missing	cm-debug-clear-malloc:	output file not generated for ndiff; 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.

strain.opstra

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.
Missing	cm-debug-clear-malloc:	output file not generated for ndiff; 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.

stress.opstre

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.
Missing	cm-debug-clear-malloc:	output file not generated for ndiff; 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:29 2016

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

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

mips-irix
SIGKILL	cm-debug-clear-malloc:	exit due to SIGKILL signal.