Example 523: Uniaxial force for 10% elongation on a unit cube

A 3D unit cube of incompressible Mooney-Rivlin material is force loaded on the X=1 face, so that it moves by 10% in the +X direction, while keeping the X=0 face fixed (uniaxial extension).


The comfile run by this example is as follows:

#Example_523  Uniaxial force for 10% elongation on a unit cube

              #This problem is set up the same as in the uniaxial 
              #extension of the unit cube except that a force is
              #applied to certain nodes in the X-direction instead of 
              #a 10% extension.

if (!defined $example) {
    $example = "./";
}
# Drop off the trailing / in the example path
$chopped = chop $example;
if ($chopped ne "/") {
    $example .= $chopped;
}
fem define coordinates;r;uniforce;example #Define 3D rectangular cartesian. 
fem define node;r;;example         #Defines 8 nodes.
fem define base;r;;example	       #Defines 3 basis functions. 
fem define element;r;;example      #Defines 1 element.            
fem def fibr;r;;example
fem def elem;r;;example fibre
fem define window
fem draw lines
fem define equation;r;;example     #Read in equation information.
fem define material;r;;example     #Read in material information.
fem define initial;r;;example      #For the initial conditions, set up the
#                              !  same as in the first uniaxial problem
#                              !  except instead of fixing nodes 2,4,6,8
#                              !  for equation 1 to 0.1, apply a force
#                              !  boundary condition to these nodes of
#                              !  magnitude that will give the same 10% 
#                              !  elongation, ie Cauchy stress*deformed
#                              !  area/# nodes or 4.486*0.9535*0.9535/4
#                              !  =1.0196.
fem define solve;r;;example        #Read in solution information.
fem solve step 1 error 10D-10  #Solve the problem.  Results should be 
#                              !  the same as in the 10% uniaxial 
#                              !  extension problem.
fem draw line deformed dotted
fem list strain at 1 ref       #List strain information.
fem list stress at 1 ref       #List stress information.

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_523.com 24-May-2004 1.9k test_output.com 24-May-2004 50 uniforce.ipbase 10-Apr-2000 3.4k uniforce.ipcoor 10-Apr-2000 619 uniforce.ipelem 10-Apr-2000 457 uniforce.ipelfb 24-May-2004 398 uniforce.ipequa 02-May-2004 2.1k uniforce.ipfibr 24-May-2004 2.1k uniforce.ipinit 26-May-2004 1.5k uniforce.ipmate 05-Dec-2002 2.3k uniforce.ipnode 10-Apr-2000 1.9k uniforce.ipsolv 16-Aug-2010 2.4k uniforce.ipsolv.old 13-Apr-2007 2.2k

Download the entire example:

Name                      Modified     Size

examples_5_52_523.tar.gz 14-Aug-2014 7.1k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:01:57 20160
cm-debugSuccessSat Mar 5 00:02:14 20161
mips-irix
cmSuccessSun Aug 19 01:24:53 20073
cm-debugSuccessWed Aug 15 01:21:27 20075
cm-debug-clear-mallocSuccessSat Aug 18 01:25:09 20079
cm-debug-clear-malloc7SuccessMon Aug 20 01:21:19 20079
cm64SuccessSun Aug 19 01:22:56 20073
cm64-debugSuccessTue Aug 21 01:20:03 20076
cm64-debug-clear-mallocSuccessTue Feb 1 09:44:09 20053
rs6000-aix
cmSuccessWed Mar 4 01:07:19 20091
cm-debugSuccessMon Mar 2 01:06:20 20091
cm64SuccessWed Mar 4 01:07:19 20091
cm64-debugSuccessTue Mar 3 01:11:43 20091
x86_64-linux
cmSuccessSun Mar 6 00:01:03 20160
cm-debugSuccessSat Mar 5 00:01:13 20160

Testing status by file:


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

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


CMISS Help / Examples / 5 / 52 / 523