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:

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:	diff test: no differences with generic answer.
  Success	cm-debug:	diff test: no differences with generic answer.
  mips-irix
  Success	cm:	diff test: no differences with generic answer.
  Success	cm-debug:	diff test: no differences with generic answer.
  Success	cm-debug-clear-malloc:	diff test: no differences with generic answer.
  Success	cm-debug-clear-malloc7:	diff test: no differences with generic answer.
  Success	cm64:	diff test: no differences with generic answer.
  Success	cm64-debug:	diff test: no differences with generic answer.
  Success	cm64-debug-clear-malloc:	diff test: no differences with generic answer.
  rs6000-aix
  Success	cm:	diff test: no differences with generic answer.
  Success	cm-debug:	diff test: no differences with generic answer.
  Success	cm64:	diff test: no differences with generic answer.
  Success	cm64-debug:	diff test: no differences with generic answer.
  x86_64-linux
  Success	cm:	diff test: no differences with generic answer.
  Success	cm-debug:	diff test: no differences with generic answer.

  
  

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

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