Example 5515: Inflation of an isotropic (Mooney-Rivlin) bilinearCubic prolate spheroid with linear interpolated hydrostatic pressure

A single prolate spheroidal element: lambda - bilinar in Xi1 and Xi2, and cubic in Xi3; mu and theta - trilinear. Hydrostatic pressure is interpolated quadratically in Xi3 (element based). The material is isotropic (Mooney-Rivlin) and incompressible. The prolate is inflated with a pressure of 5kPa on the inside face while the the pressure is held at zero on the outside face.

The comfile run by this example is as follows:

#Example_5515:  inflation of a prolate (lambda bilinearCubic and mu,theta trilinear.)                       

fem define para;r;ellipsoid

fem define coor;r;
fem define nodes;r;
#fem define bases;r;3CubicHermite
fem define bases;r;2Linear-CubicHermite
fem define ;add bases;r;3Linear
fem define ;add bases;r;pressure
fem define ;add bases;r;2Linear
fem define ;add bases;r;Linear-CubicHermite
fem define ;add bases;r;CubicHermite-Linear
fem define elements;r;ellipsoid
#
#
# Define element groups
fem group elem 1 as apexEndocardElem
fem group elem 1 as apexEpicardElem
fem group elem 1 as apexElem
fem group elem 1 as baseEndocardElem
fem group elem 1 as baseEpicardElem
fem group elem 1 as endocardElem
fem group elem 1 as epicardElem
#
#
# Define node groups
#
fem group node xi2=1 xi3=1 external elem baseEpicardElem as baseEpicardNodes
fem group node xi2=0 external elem apexElem as apexNodes
fem group node xi3=0 external elem endocardElem as endocardNodes
fem group node xi3=1 external elem epicardElem as epicardNodes
fem group node xi1=0 xi2=1 xi3=1 external elem baseEpicardElem as constrainedThetaNode
fem list nodes group
#
#
fem define fibre;d
fem define element;r;ellipsoid fibre
fem define equation;r;ellipsoid
fem define material;r;
fem define initial;r;
fem define solve;r;ellipsoid
read com;solveInflation

# Run "cmgui view" to view the result

Files used by this example are:

Name                         Modified     Size

example_5515.com             11-Sep-2000  1.3k
2CubicHermite.ipbase         08-Sep-2000  1.5k
2Linear-CubicHermite.ipbase  08-Sep-2000  1.7k
2Linear.ipbase               08-Sep-2000  1.1k
3CubicHermite.ipbase         08-Sep-2000  1.9k
3Linear.ipbase               08-Sep-2000  1.4k
CubicHermite-Linear.ipbase   08-Sep-2000  1.5k
Linear-CubicHermite.ipbase   08-Sep-2000  1.5k
ellipsoid.com                08-Sep-2000  1.2k
ellipsoid.ipcoor             08-Sep-2000  725 
ellipsoid.ipelem             08-Sep-2000  598 
ellipsoid.ipelfb             08-Sep-2000  273 
ellipsoid.ipequa             02-May-2004  2.1k
ellipsoid.ipfibr             30-Jan-2001  5.1k
ellipsoid.ipinit             12-Dec-2002  1.8k
ellipsoid.ipmate             12-Dec-2002  2.2k
ellipsoid.ipnode             08-Sep-2000  1.5k
ellipsoid.ippara             12-Nov-2002  5.9k
ellipsoid.ipsolv             16-Aug-2010  2.7k
ellipsoid.ipsolv.old         13-Apr-2007  2.6k
pressure.ipbase              08-Sep-2000  1.9k
solveInflation.com           11-Sep-2000  1.0k
view.com                     17-Sep-2000  2.0k

Download the entire example:

Name                           Modified     Size

examples_5_55_551_5515.tar.gz  17-Aug-2010  9.8k

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

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