Example 5513: Inflation of a trilinear prolate element of incomp + fluid shift pole-zero material

Finite elasticity analysis: 2kPa inflation of a prolate element with rotating fibre and sheet fields and an incompressible + fluid shift pole-zero material law. Geometry is interpolated isochorically (Focus^3.Cosh[lambda].Sinh^2[lambda]).Hydrostatic pressure varies quartically in the transmural direction. Uses Newton's method with GMRES to solve the sparsely stored linear equations.


The comfile run by this example is as follows:

#Example_5513  Inflation of a prolate element incomp + fluid shift pole-zero
#             material

fem
fem define coordinates;r;prolfluid;example #3D prolate spheroidal
fem define nodes;r;;example                #4 nodes of previous example
fem define bases;r;;example                #trilinear + quartic hydr. pressure
fem define element;r;;example              #single axisymmetric element
fem define fibre;r;;example                #fibre angles of prev example
fem define element;r;;example fibre        #fibre elements
fem define window
fem draw lines
fem define equation;r;;example             #3D, finite elasticity, FEM
fem define material;r;;example             #incomp+fluid; pole-zero law
fem define initial;r;;example              #inflate to 2kPa
fem define solve;r;;example                #Newton with GMRES
fem solve increment 0.1 iterate 3
fem solve increment 0.4 iterate 3
fem solve increment 0.5 iterate 12         #should take 10 iterations here 
fem draw lines def dot
fem list strain xi 3 p 3                   #transmural strains
fem list stress xi 3 p 3                   #transmural stresses

Files used by this example are:

Name                  Modified     Size

example_5513.com 10-Apr-2000 1.1k prolfluid.ipbase 10-Apr-2000 4.3k prolfluid.ipcoor 10-Apr-2000 707 prolfluid.ipelem 10-Apr-2000 409 prolfluid.ipequa 02-May-2004 2.2k prolfluid.ipfibr 30-Jan-2001 1.3k prolfluid.ipinit 12-Dec-2002 2.0k prolfluid.ipmate 12-Dec-2002 5.9k prolfluid.ipnode 10-Apr-2000 1.2k prolfluid.ippara 17-Dec-2002 5.0k prolfluid.ipshee 10-Apr-2000 781 prolfluid.ipsolv 16-Aug-2010 2.6k prolfluid.ipsolv.old 13-Apr-2007 2.4k test_output_old.com 10-Apr-2000 17

Download the entire example:

Name                           Modified     Size

examples_5_55_551_5513.tar.gz 17-Aug-2010 8.6k

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

Input last modified: Mon Aug 16 11:19:21 2010


CMISS Help / Examples / 5 / 55 / 551 / 5513