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


CMISS Help / Examples / 5 / 55 / 551 / 5515