Inflation of an isotropic (Mooney-Rivlin) trilinear prolate spheroid

Example 5511: Inflation of an isotropic (Mooney-Rivlin) trilinear prolate spheroid

A single prolate spheroidal element is defined to be incompressible and isotropic (Mooney-Rivlin material law). Geometry is interpolated isochorically (Focus^3.Cosh[lambda].Sinh^2[lambda]).Hydrostatic pressure is interpolated quadratically in Xi3 (element based). 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.

The comfile run by this example is as follows:

#Example_5511  Inflation of an isotropic Mooney-Rivlin prolate spheroid 

fem
fem define coordinate;r;prolisot;example #Defines 3 GLOBAL PROLATE SPHEROIDAL
#                              !  coordinates. 
#                              !  Transmural interpolation is in
#                              !  Focus^3.Cosh(Lamda).Sinh^2(Lamda).
fem define node;r;;example         #Focus = 37.5. Reads 4 NODES, 3 COORDS,
#                              !  1 VERSION per node per nj, 0 
#                              !  DERIVATIVES for all coords. Nodal 
#                              !  coordinates, Xj, are (lamda,mu,theta):
#                              !  1=(0.43,0,0);2=(0.43,120,0);
#                              !  3=(0.72,0,0); 4=(0.72,120,0).
fem define base;r;;example         #Defines 3 basis functions: 
#                              !  1) Lagrange/Hermite with 3 Xi coords,
#                              !     Linear Lagrange interpolant and 1
#                              !     Gauss point in Xi1 dirctn (theta-
#                              !     circumferential) but 3 in Xi2 (mu-
#                              !     longitudinal) and Xi3 (lamda-
#                              !     transmural);
#                              !  2) An auxiliary basis with 5 auxiliary
#                              !     element parameters, pressure basis
#                              !     same number of Gauss points as for
#                              !     basis function 1, and polynomial
#                              !     degrees of: 0 for Xi1 and Xi2 for
#                              !     all 5 parameters; 0 in Xi3 for 
#                              !     parameter 1; 1 for parameter 2; 2 
#                              !     for parameter 3; -1 for parameter
#                              !     4 (Xi3=0 face pressure bc); & -2
#                              !     for parameter 5 (Xi3=1 face 
#                              !     pressure bc).
#                              !  3) 2 Xi coords with linear Lagrange
#                              !     interpolant, and 3 Gauss points
#                              !     in Xi1 and Xi2 directions.
fem define element;r;;example      #Defines 1 prolate spheroidal element.
fem define fibr;d
fem define elem;r;;example fibre
fem define window;r;;example       #Defines window dimensions in X,Y,Z
#                              !  directions as (min,max):(-5,5),(-5,5),
#                              !  (-5,5).
fem draw lines                 #Makes line segments visible on window.
fem define equation;r;;example     #Static 3D finite elasticity problem with
#                              !  geometric coords as the dependent 
#                              !  variables. Solution is by Galerkin
#                              !  finite elements in a nonlinear
#                              !  analysis of a nonlinear equation.
#                              !  Incompressible material with basis 1
#                              !  in all elements for dep vars 1,2,3
#                              !  and basis 2 in all elements for dep
#                              !  var 4 (hydrostatic pressure).   
fem define material;r;;example     #Stresses in constitutive law are
#                              !  referred to body coords.  Defines a
#                              !  hyperelastic constitutive law for
#                              !  which the strain energy function, W,
#                              !  is a polynomial function of principal 
#                              !  strain invariants. The Mooney-Rivlin 
#                              !  material law is epressed in terms of
#                              !  I1 and I2 to the maximum power of 1.
#                              !  The 4 constant parameters are 0,2,6,0
#                              !  defining the strain energy function as 
#                              !  W = 2*(I1-3) + 6*(I2-3)
fem define initial;r;;example      #Boundary pressure increments are entered
#                              !  and hydrostatic pressure is matched
#                              !  across boundaries with adjacent Xi3
#                              !  faces. Initial displacements are all
#                              !  zero. For dependent variable 2
#                              !  (mu-dirn) fix all nodes (1..4);
#                              !  Dependent variable 3 (theta-dirn), 
#                              !  fix node 2. No force bcs.
#                              !  Dependent variable 4 (hyd press) fix
#                              !  element 1, parameter 4 (Xi3=0 face
#                              !  inside prolate) with increment 
#                              !  magnitude 2, and parameter 5 (Xi3=1
#                              !  face outside prolate) with increment  
#                              !  magnitude 0. No force bcs.
fem define solve;r;;example        #  Read in solution information.
fem solve increment 0.2 iter 3 #Solve the problem.
fem solve increment 0.3 iter 3 #Incrementing of the ventricular 
#                              !  pressure is neccessary to avoid
fem solve increment 0.5        #  unstable solutions. This depends on
#                              !  the constitutive law and pressure
#                              !  magnitude, as the wall is gradually
#                              !  stiffened.
fem draw lines def dotted      #Draw deformed mesh on window.
fem list strain at 1 ref       #List strain information wrt reference
#                              !  coordinates.
fem list stress at 1 ref       #List stress information wrt reference
#                              !  coordinates.
#                                    
#              !Note that if node 1, rather than node 2, is fixed in the
#              !  theta-direction, unpredictable results will follow 
#              !  because the moment needed to fix a node in this
#              !  direction has no moment arm (on axis).
#
#              !Note that volume is not conserved even with a quadratic
#              !  hydrostatic pressure (I3 not equal to 1).

Additional testing commands:

fem list stress at 1 ref
fem list strain at 1 ref

Files used by this example are:

Name                 Modified     Size

example_5511.com     17-Dec-2002  6.0k
prolisot.ipbase      10-Apr-2000  4.2k
prolisot.ipcoor      10-Apr-2000  707 
prolisot.ipelem      10-Apr-2000  409 
prolisot.ipelfb      10-Apr-2000  232 
prolisot.ipequa      02-May-2004  2.0k
prolisot.ipfibr      30-Jan-2001  738 
prolisot.ipinit      12-Dec-2002  1.7k
prolisot.ipmate      12-Dec-2002  2.3k
prolisot.ipnode      10-Apr-2000  1.2k
prolisot.ipsolv      16-Aug-2010  2.3k
prolisot.ipsolv.old  13-Apr-2007  2.1k
prolisot.ipwind      10-Apr-2000  262 
test_output.com      17-Dec-2002   50

Download the entire example:

Name                           Modified     Size

examples_5_55_551_5511.tar.gz  14-Aug-2014  8.2k

Testing status by version:

	Status	Tested	Real time (s)
i686-linux
cm	Success	Sun Mar 6 00:02:01 2016	1
cm-debug	Success	Sat Mar 5 00:02:19 2016	1
mips-irix
cm	Success	Sun Aug 19 01:30:08 2007	5
cm-debug	Success	Wed Aug 15 01:27:02 2007	10
cm-debug-clear-malloc	Success	Sat Aug 18 01:31:59 2007	14
cm-debug-clear-malloc7	Success	Mon Aug 20 01:29:54 2007	14
cm64	Success	Sun Aug 19 01:30:50 2007	5
cm64-debug	Success	Tue Aug 21 01:25:29 2007	11
cm64-debug-clear-malloc	Success	Tue Feb 1 09:45:17 2005	6
rs6000-aix
cm	Success	Wed Mar 4 01:07:25 2009	1
cm-debug	Success	Mon Mar 2 01:07:27 2009	2
cm64	Success	Wed Mar 4 01:07:25 2009	1
cm64-debug	Success	Tue Mar 3 01:12:58 2009	3
x86_64-linux
cm	Success	Sun Mar 6 00:01:03 2016	0
cm-debug	Success	Sat Mar 5 00:01:14 2016	0

Testing status by file:

cmiss_test.log

i686-linux
Success	cm:	cmiss_test.log.retain.
Success	cm-debug:	cmiss_test.log.retain.
mips-irix
Success	cm:	cmiss_test.log.retain.
Success	cm-debug:	cmiss_test.log.retain.
Success	cm-debug-clear-malloc:	cmiss_test.log.retain.
Success	cm-debug-clear-malloc7:	cmiss_test.log.retain.
Success	cm64:	cmiss_test.log.retain.
Success	cm64-debug:	cmiss_test.log.retain.
Success	cm64-debug-clear-malloc:	cmiss_test.log.retain.
rs6000-aix
Success	cm:	cmiss_test.log.retain.
Success	cm-debug:	cmiss_test.log.retain.
Success	cm64:	cmiss_test.log.retain.
Success	cm64-debug:	cmiss_test.log.retain.
x86_64-linux
Success	cm:	cmiss_test.log.retain.
Success	cm-debug:	cmiss_test.log.retain.

cmiss_test.out

i686-linux
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
mips-irix
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
Success	cm-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
Success	cm-debug-clear-malloc7:	ndiff test: no significant differences with generic answer.
Success	cm64:	ndiff test: no significant differences with generic answer.
Success	cm64-debug:	ndiff test: no significant differences with generic answer.
Success	cm64-debug-clear-malloc:	ndiff test: no significant differences with generic answer.
rs6000-aix
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.
Success	cm64:	ndiff test: no significant differences with generic answer.
Success	cm64-debug:	ndiff test: no significant differences with generic answer.
x86_64-linux
Success	cm:	ndiff test: no significant differences with generic answer.
Success	cm-debug:	ndiff test: no significant differences with generic answer.

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

Input last modified: Thu Aug 14 10:03:14 2014

CMISS Help / Examples / 5 / 55 / 551 / 5511