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