Example 848: Inflation followed by contraction of a prolate spheroid with coronaries

A single prolate spheroidal element is defined with fibres oriented 65 degrees wrt the horizontal on the inside and -55 degrees wrt the horizontal on the outside. Sheet angles are defined as -45 degrees and +45 degrees at inner and outer nodes respectively. Geometry is interpolated isochorically (Focus^3.Cosh[lambda].Sinh^2[lambda]). Hydrostatic pressure is interpolated quadratically in Xi3 (element based) and the pole-zero law defines the incompressible material behaviour. The prolate is inflated with a pressure of 2kPa on the inside face while the the pressure is held at zero on the outside face. It is then contracted with active fibre stresses (steady state Calcium-length-tension relation). the solution starts from pressure inflated solution of 2kPa. Newton's method is used to solve the sparsely stored linear equations. Note this is based on example 545 with the maxium Ca concentration reduced to maintain stability. Coronary vessels are added in region 2.


The comfile run by this example is as follows:

#Example_848  Inflation followed by contraction of an orthotropic prolate spheroid
# with fibres and sheets in region 1 and coronaries region 2
 
fem
fem de para;r;coro;example
fem reallocate
fem de reg;r;coro;example
fem define coordinate;r;profib;example reg 1
fem define coordinate;r;profib;example reg 2
fem define node;r;profib;example  reg 1
fem define base;r;profib;example
fem de;add base;r;prolisot2;example
fem define element;r;profib;example reg 1
fem define fibre;r;profib;example reg 1
fem define element;r;profib;example fibre  reg 1
fem de no;r;coronary;example reg 2
fem de ele;r;;example reg 2
fem de xi;r;;example  node reg 2
fem de fiel;r;;example reg 2
fem de ele;r;;example field reg 2
FEM update mesh position of 2 in 1 undeform
fem export node;rest as time reg 2
fem export ele;rest as rest reg 2
#
fem define equation;r;profib;example reg 1 class 1
fem define material;r;profib;example reg 1 class 1
# Inflation
fem define initial;r;profib;example  reg 1 class 1
fem define solve;r;profib;example   reg 1 class 1
fem solve increment 0.2 iter 3  class 1
fem solve increment 0.3 iter 3  class 1
fem solve increment 0.5         class 1
FEM update mesh position of 2 in 1 deform
fem export node;inflated as time reg 2
fem export ele;inflated as inf reg 2
#  Contraction
fem define material;r;prolcontract;example class 1
fem define active;r;prolcontract;example class 1
fem define solve;r;prolcontract;example class 1
fem solve increment 0.0 iterate 10 class 1
FEM update mesh position of 2 in 1 deform
fem export node;contract as time reg 2
fem export ele;contract as con reg 2
 
 
 
 

Files used by this example are:

Name                       Modified     Size

example_848.com 10-Apr-2000 1.6k contract.exelem 10-Apr-2000 475k contract.exnode 10-Apr-2000 47k coro.ippara 12-Nov-2002 5.9k coro.ipregi 10-Apr-2000 128 coronary.ipelem 10-Apr-2000 119k coronary.ipelfd 10-Apr-2000 53k coronary.ipfiel 10-Apr-2000 34k coronary.ipgrid 27-Nov-2001 216 coronary.ipnode 10-Apr-2000 83k coronary.ipxi 10-Apr-2000 24k cylinder_flow_mech.ippara 12-Nov-2002 5.9k inflated.exelem 10-Apr-2000 45k inflated.exnode 10-Apr-2000 36k profib.ipbase 10-Apr-2000 4.2k profib.ipcoor 10-Apr-2000 958 profib.ipelem 10-Apr-2000 409 profib.ipelfb 10-Apr-2000 353 profib.ipequa 02-May-2004 2.1k profib.ipfibr 30-Jan-2001 1.3k profib.ipinit 10-Apr-2000 1.7k profib.ipmate 12-Dec-2002 5.7k profib.ipnode 10-Apr-2000 1.2k profib.ipshee 10-Apr-2000 781 profib.ipsolv 16-Aug-2010 2.3k profib.ipsolv.old 13-Apr-2007 2.2k profib.ipwind 10-Apr-2000 262 prolcontract.ipacti 03-Mar-2004 727 prolcontract.ipmate 12-Dec-2002 6.2k prolcontract.ipsolv 16-Aug-2010 2.3k prolcontract.ipsolv.old 13-Apr-2007 2.1k prolisot2.ipbase 10-Apr-2000 893 rest.exelem 10-Apr-2000 45k rest.exnode 10-Apr-2000 36k test.movie 10-Apr-2000 367k test2.ipelem 10-Apr-2000 119k test2.ipelfd 10-Apr-2000 53k test2.ipfiel 10-Apr-2000 34k test2.ipgrid 27-Nov-2001 216 test2.ipnode 10-Apr-2000 83k test2.ippara 12-Nov-2002 5.9k test2.ipxi 10-Apr-2000 24k test_output.comold 10-Apr-2000 50 view.com 10-Apr-2000 758

Download the entire example:

Name                      Modified     Size

examples_8_84_848.tar.gz 17-Aug-2010 297k

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

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


CMISS Help / Examples / 8 / 84 / 848