Example 532: Pressure inflation and axial extension of a cylinder using r2 interpolation

A trilinear cylinder (r^2 interpolation) of incompressible Mooney-Rivlin material is inflated to 1kPa.


The comfile run by this example is as follows:

#Example_532 Inflation of a 2 element cylinder using r2 interpolation
#
#
#           N   N  OO  TTT EEE
#           NN  N O  O  T  E
#           N N N O  O  T  EEE
#           N  NN O  O  T  E
#           N   N  OO   T  EEE
#
#  ERROR: There is a problem using r^2 interpolation and pressure boundary
#         conditions - produces singular arrays and needs fixing
#
#
#
#
 
fem
fem define coordinates;r;cylind2;example #3D cylindrical polar coordinates,
#                              !  with r^2 interpolation.
fem define node;r;;example         #Defines 6 nodes.
fem define base;r;;example         #Defines 3 basis functions.
fem define element;r;;example      #Defines 2 elements.
fem define window;r;;example       #Defines window dimensions.
fem draw lines                 #Draws undeformed mesh on window.
fem define fibre;d                 #Default fibre angles
fem define element;r;;example fibre
fem define equation;r;;example     #Defines finite elasticity equations.
fem define material;r;;example     #Incompressible Mooney-Rivlin material.
fem define initial;r;;example      #Internal pressure of 1, external 0
#                              !  constrain axial movement at centre
#                              !  nodes.
fem define solve;r;;example        #Residual and intermediate output.
fem solve step 1
fem draw lines def dot         #Draw deformed mesh.
fem list strain at 1 ref       #List strain wrt reference coordinates
#                              !  at 1st Gauss point.
fem list stress at 1 ref       #List stress wrt reference coordinates
#                              !  at 1st Gauss point.
fem list elem tot
fem list elem def tot
             #Note the differences in strain and stress results in
             #  comparison to the previous problem.  The third strain
             #  invariant indicates that volume is conserved (I3=1);
             #  the inner radial stress is close to -1 while the outer
             #  radial stress is close to zero (when all Gauss points
             #  results are displayed).

Additional testing commands:

#ERROR: This example doesn't work properly (check r^2 interpolation stuff)
# list solve is temporary test output until problem fixed
fem list solve
#fem list stress at 1 ref
#fem list strain at 1 ref

Files used by this example are:

Name                Modified     Size

example_532.com 10-Apr-2000 2.0k cylind2.ipbase 10-Apr-2000 4.2k cylind2.ipcoor 10-Apr-2000 647 cylind2.ipelem 10-Apr-2000 2.0k cylind2.ipelfb 10-Apr-2000 382 cylind2.ipequa 02-May-2004 2.0k cylind2.ipfibr 30-Jan-2001 534 cylind2.ipinit 05-Dec-2002 1.7k cylind2.ipmate 05-Dec-2002 2.2k cylind2.ipnode 10-Apr-2000 1.5k cylind2.ipsolv 16-Aug-2010 2.3k cylind2.ipsolv.old 13-Apr-2007 2.1k cylind2.ipwind 10-Apr-2000 263 test_output.com 10-Apr-2000 200

Download the entire example:

Name                      Modified     Size

examples_5_53_532.tar.gz 17-Aug-2010 7.0k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:01:59 20160
cm-debugSuccessSat Mar 5 00:02:17 20161
mips-irix
cmSuccessSun Aug 19 01:31:13 20075
cm-debugSuccessWed Aug 15 01:29:51 200713
cm-debug-clear-mallocSuccessSat Aug 18 01:36:01 200718
cm-debug-clear-malloc7SuccessMon Aug 20 01:33:31 200719
cm64SuccessSun Aug 19 01:31:40 20075
cm64-debugSuccessTue Aug 21 01:28:02 200713
cm64-debug-clear-mallocSuccessTue Feb 1 09:45:26 20058
rs6000-aix
cmSuccessWed Mar 4 01:07:23 20091
cm-debugSuccessMon Mar 2 01:07:26 20094
cm64SuccessWed Mar 4 01:07:23 20091
cm64-debugSuccessTue Mar 3 01:12:58 20094
x86_64-linux
cmSuccessSun Mar 6 00:01:03 20160
cm-debugSuccessSat Mar 5 00:01:14 20161

Testing status by file:


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

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


CMISS Help / Examples / 5 / 53 / 532