xz shearing of a Mooney-Rivlin unit cube with embeded coronary mesh

Example 847: xz shearing of a Mooney-Rivlin unit cube with embeded coronary mesh

This example demonstrates the solution to the coupled flow/solid mechanics problem. The flow network, comprising of elastic tubes (four, 1-d cubic Hermite elements) is embedded in a homogeneous Mooney-Rivlin cube. The cube is deformed in simple shear in the x-z plane in steps over a period of 0.005s. The fluid flow pressure at the inlet of the tube is increased from 10.6 - 14.6 kPa linearly over the same period while the exit flow pressure is maintained constant at 10.6 kPa. For the theory and further details of this problem, see section 5.3 of the thesis "Coronary flow mechanics" by N.P.Smith.

Simulation of coupled flow/solid mechanics

The comfile run by this example is as follows:

# Note : region 1 - solid cube; region 2 - fluid flow
#        class 1 - Finite elasticity problem ; class 2 - Flow thru elastic tube

fem define para;r;cube_coro;example
fem reallocate
fem define region;r;cube_coro;example              # Defines no. of regions
fem define coordinate;r;cube1;example reg 1        # Defines coordinates for solid
fem define node;r;cube1;example reg 1              # Reads in nodal information  
fem define base;r;cube1;example reg 1              # Defines basis functions (05) 
fem define element;r;cube1;example reg 1           # Reads in element information
fem define fibre;r;cube1;example reg 1             # Reads in fibre information for solid 
fem define element;r;cube1;example fib reg 1       # Defines fibre interpolation      
fem define equation;r;cube1;example reg 1 class 1  # Defines equations for solid
fem define material;r;cube1;example reg 1 class 1  # Reads in mech properties of solid
fem define initial;r;cube1;example reg 1 class 1   # Defines initial & boundary conds.
fem define solve;r;cube1;example reg 1 class 1     # Defines solution parameters
 
fem define coordinates;r;coro;example reg 2        # Defines coordinates for fluid flow
fem define node;r;coro;example reg 2               # Reads in nodal information 
fem define;add base;r;coro;example                 # Defines basis functions (01) 
fem define field;r;coro;example reg 2              # Reads in field information 
fem define element;r;coro;example reg 2            # Reads in element information

fem define xi;r;coro;example node reg 2            # Reads in xi positions of reg 2 wrt reg 1
fem update mesh position of 2 in 1 undeform        # Updates relative positions

fem define element;r;coro;example field reg 2      # Defines field interpolation     
fem define grid;r;coro;example coronary reg 2      # Reads in spatial discretisation    
fem update grid geometry reg 2                     # Updates grid point coordinates 
fem define equation;r;coro;example reg 2 class 2   # Defines equations for fluid flow 
fem define material;r;coro;example reg 2 class 2   # Reads in phy. & mech. properties   
fem update material field material 4 no_field 1 reg 2 class 2
#                                                  # Updates unstressed radii
fem define time;r;coro;example                     # Defines boundary conditions       
fem define initial;r;coro;example reg 2 class 2    # Defines initial conditions           
fem define solve;r;coro;example reg 2 class 2      # Defines solution parameters
fem define coup;r;coro;example                     # Defines solid & fluid mech. coupling

$time_initial=0.000;
$time_final=0.005;
$time_step=0.0000008;   
$write_out_period=0.0001;

fem export nodes;solid${time_initial} as undeformed offset 500 reg 1 class 1
fem export elements;solid_undef as undeformed offset_elem 500 offset_node 500 reg 1 class 1 
fem export elements;solid_def as deformed offset_elem 600 offset_node 600 reg 1 class 1 
fem export grid;fluid as flow offset_elem 0 offset_node 0 elastic_tube reg 2 class 2
fem export points;fluid${time_initial} grid as flow offset 0 reg 2 class 2
 
for ($time=$time_initial;$time<$time_final;$time=$time+$write_out_period)      
{
$w_time=$time+$write_out_period;
fem solve restart to $write_out_period delta_t $time_step reg 1,2 class 1,2;
fem update geometry from solution reg 1;
fem export nodes;solid${w_time} as deformed offset 600 reg 1 class 1;
fem update nodes derivative 1 linear reg 2 
fem update mesh position of 2 in 1 deform
fem update grid geometry reg 2
fem export points;fluid${w_time} grid as flow offset 0 reg 2 class 2
}

Additional testing commands:

fem list nodes reg 2
fem list elem reg 2

Files used by this example are:

Name                Modified     Size

example_847.com     04-May-2005  3.5k
coro.ipbase         26-Sep-2001  2.0k
coro.ipcoor         26-Sep-2001  601 
coro.ipcoup         21-Aug-2002  769 
coro.ipcoup.old     26-Sep-2001  743 
coro.ipelem         26-Sep-2001  1.5k
coro.ipelfd         26-Sep-2001  1.1k
coro.ipequa         22-Nov-2005  1.3k
coro.ipfiel         26-Sep-2001  1.6k
coro.ipgrid         04-May-2005  250 
coro.ipinit         01-Dec-2005  723 
coro.ipmate         26-Sep-2001  6.8k
coro.ipnode         26-Sep-2001  2.3k
coro.ipsolv         27-Aug-2002  799 
coro.iptime         26-Sep-2001  1.4k
coro.ipxi           26-Sep-2001  325 
cube1.ipbase        26-Sep-2001  3.4k
cube1.ipcoor        26-Sep-2001  602 
cube1.ipelem        26-Sep-2001  433 
cube1.ipelfb        26-Sep-2001  435 
cube1.ipequa        02-May-2004  2.1k
cube1.ipfibr        26-Sep-2001  2.0k
cube1.ipinit        12-Dec-2002  1.5k
cube1.ipmate        12-Dec-2002  2.3k
cube1.ipnode        26-Sep-2001  1.9k
cube1.ipsolv        16-Aug-2010  2.4k
cube1.ipsolv.old    13-Apr-2007  2.2k
cube_coro.ippara    12-Nov-2002  5.9k
cube_coro.ipregi    26-Sep-2001  136 
fluid.exelem        04-May-2005  3.4k
fluid0.0005.exnode  04-May-2005  1.9k
fluid0.001.exnode   04-May-2005  1.9k
fluid0.0015.exnode  04-May-2005  1.9k
fluid0.002.exnode   04-May-2005  1.9k
fluid0.0025.exnode  04-May-2005  1.9k
fluid0.003.exnode   04-May-2005  1.9k
fluid0.0035.exnode  04-May-2005  1.9k
fluid0.004.exnode   04-May-2005  1.9k
fluid0.0045.exnode  04-May-2005  1.9k
fluid0.005.exnode   04-May-2005  1.9k
fluid0.exnode       04-May-2005  1.9k
solid0.0005.exnode  04-May-2005  1.9k
solid0.001.exnode   04-May-2005  1.9k
solid0.0015.exnode  04-May-2005  1.9k
solid0.002.exnode   04-May-2005  1.9k
solid0.0025.exnode  04-May-2005  1.9k
solid0.003.exnode   04-May-2005  1.9k
solid0.0035.exnode  04-May-2005  1.9k
solid0.004.exnode   04-May-2005  1.9k
solid0.0045.exnode  04-May-2005  1.9k
solid0.005.exnode   04-May-2005  1.9k
solid0.exnode       26-Sep-2001  1.9k
solid_def.exelem    26-Sep-2001  5.8k
solid_undef.exelem  26-Sep-2001  5.8k
test_output.com     27-Sep-2001   41

Download the entire example:

Name                      Modified     Size

examples_8_84_847.tar.gz  17-Aug-2010  259k

Testing status by version:

	Status	Tested	Real time (s)
i686-linux
cm	Success	Sun Mar 6 00:02:07 2016	2
cm-debug	Success	Sat Mar 5 00:08:01 2016	3
mips-irix
cm	Success	Sun Aug 19 01:48:48 2007	18
cm-debug	Success	Wed Aug 15 01:44:20 2007	45
cm-debug-clear-malloc	Success	Sat Aug 18 01:58:39 2007	69
cm-debug-clear-malloc7	Success	Mon Aug 20 01:52:45 2007	56
cm64	Success	Sun Aug 19 01:48:57 2007	18
cm64-debug	Success	Tue Aug 21 01:43:15 2007	44
cm64-debug-clear-malloc	Success	Thu Apr 1 10:51:42 2004	12
rs6000-aix
cm	Success	Wed Mar 4 01:09:40 2009	3
cm-debug	Success	Mon Mar 2 01:11:01 2009	9
cm64	Success	Wed Mar 4 01:09:33 2009	4
cm64-debug	Success	Tue Mar 3 01:16:19 2009	10
x86_64-linux
cm	Success	Sun Mar 6 00:01:04 2016	1
cm-debug	Success	Sat Mar 5 00:01:47 2016	1

Testing status by file:

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.

cmiss_test.out

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.
Success	cm-debug-clear-malloc:	ndiff test: no significant differences with 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.

fluid0.005.exnode

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.
Success	cm-debug-clear-malloc:	ndiff test: no significant differences with 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:34 2016

Input last modified: Mon Aug 16 11:22:34 2010

CMISS Help / Examples / 8 / 84 / 847