Example 536: Pressure inflation of a cylinder with free ends using r2 interpolation

A trilinear cylinder (r^2 interpolation) of incompressible Mooney-Rivlin material is pressure inflated by 1kPa while leaving the ends free.


The comfile run by this example is as follows:

#Example_536  Pressure inflation of a 2 element cylinder with free ends
#             using r^2 interpolation

#
#  ERROR: There is a problem using r^2 interpolation and pressure boundary
#         conditions - produces singular arrays and needs fixing
#

fem 
fem define coordinate;r;cylinp;example #Defines 3 cylindrical polar coords,
#                              !  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 the window dimensions.
fem draw lines                 #Draws undeformed mesh.
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.
fem define solve;r;;example        #Residual and intermediate output.
fem solve step 1
fem draw lines deformed dotted #Draws 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.

Files used by this example are:

Name                 Modified     Size

example_536.com 10-Apr-2000 1.3k cylinp.ipbase 10-Apr-2000 4.2k cylinp.ipcoor 10-Apr-2000 646 cylinp.ipelem 10-Apr-2000 2.0k cylinp.ipequa 02-May-2004 2.9k cylinp.ipfibr 30-Jan-2001 1.9k cylinp.ipinit 05-Dec-2002 1.4k cylinp.ipmate 05-Dec-2002 2.3k cylinp.ipnode 10-Apr-2000 4.3k cylinp.ipsolv 16-Aug-2010 2.2k cylinp.ipsolv.old 13-Apr-2007 2.0k cylinp.ipwind 10-Apr-2000 262 test.ipequa 02-May-2004 2.6k test_output_old.com 10-Apr-2000 200

Download the entire example:

Name                      Modified     Size

examples_5_53_536.tar.gz 17-Aug-2010 6.8k

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 / 536