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:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:02:07 20162
cm-debugSuccessSat Mar 5 00:08:01 20163
mips-irix
cmSuccessSun Aug 19 01:48:48 200718
cm-debugSuccessWed Aug 15 01:44:20 200745
cm-debug-clear-mallocSuccessSat Aug 18 01:58:39 200769
cm-debug-clear-malloc7SuccessMon Aug 20 01:52:45 200756
cm64SuccessSun Aug 19 01:48:57 200718
cm64-debugSuccessTue Aug 21 01:43:15 200744
cm64-debug-clear-mallocSuccessThu Apr 1 10:51:42 200412
rs6000-aix
cmSuccessWed Mar 4 01:09:40 20093
cm-debugSuccessMon Mar 2 01:11:01 20099
cm64SuccessWed Mar 4 01:09:33 20094
cm64-debugSuccessTue Mar 3 01:16:19 200910
x86_64-linux
cmSuccessSun Mar 6 00:01:04 20161
cm-debugSuccessSat Mar 5 00:01:47 20161

Testing status by file:


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