A trilinear cylinder (r^2 interpolation) of incompressible Mooney-Rivlin material is pressure inflated by 1kPa while leaving the ends free.
#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.
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
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