Example 527: 10% uniaxial extension of a unit cube with fibres (90 degrees)

A 3D unit cube of incompressible transversely isotropic material (exponential material law) where the fibre direction is at a constant 90 degree angle wrt the X axis, is displacement loaded on the X=1 face 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_527  10% uniaxial extension of a unit cube with fibers (90 degrees)
 
fem
fem define coordinates;r;unifib90;example #Defines 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 define fibre;r;;example         #Define fibre orientations in the 1-2
#                               !  plane for all 8 nodes in DEGREES wrt
#                               !  Xi 1 direction as 90 degrees.
fem define element;r;;example fibre
fem define window
fem draw lines
fem draw fibres
fem define equation;r;;example      #Defines equation information.
fem define material;r;;example      #Defines material information.
fem define initial;r;;example       #Defines initial conditions.
fem define solve;r;;example         #Defines solution information.
fem solve step 1 error 1.0D-10  #Solve the problem.
fem draw line deformed dotted
fem list strain at 1 ref        #List strain information wrt reference
#                               !  coordinates at 1st gauss point.
fem list stress at 1 ref        #List stress information wrt reference
#                               !  coordinates at 1st gauss point.
 
              #For this 90 degree orientation, the stress in the
              #  X-direction is the lowest out of all 4 directions. The
              #  amount of shortening in the Y-direction (fibre
              #  direction) is not the same as in the Z-direction
              #  because the stiffnesses are different.

Additional testing commands:

# There are repeated eigenvalues so Euler angles and eigenvectors are not unique.
fem list stress at 1 ref
fem list strain at 1 ref

Files used by this example are:

Name                 Modified     Size

example_527.com 10-Apr-2000 1.5k test_output.com 10-Apr-2000 132 unifib90.ipbase 10-Apr-2000 3.4k unifib90.ipcoor 10-Apr-2000 608 unifib90.ipelem 10-Apr-2000 446 unifib90.ipelfb 10-Apr-2000 275 unifib90.ipequa 02-May-2004 2.1k unifib90.ipfibr 30-Jan-2001 1.1k unifib90.ipinit 05-Dec-2002 2.9k unifib90.ipmate 13-Oct-2003 3.0k unifib90.ipnode 10-Apr-2000 1.9k unifib90.ipsolv 16-Aug-2010 2.3k unifib90.ipsolv.old 13-Apr-2007 2.2k

Download the entire example:

Name                      Modified     Size

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

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:01:58 20160
cm-debugSuccessSat Mar 5 00:02:15 20160
mips-irix
cmSuccessSun Aug 19 01:24:58 20073
cm-debugSuccessWed Aug 15 01:21:12 20076
cm-debug-clear-mallocSuccessSat Aug 18 01:24:28 20079
cm-debug-clear-malloc7SuccessMon Aug 20 01:22:41 20079
cm64SuccessSun Aug 19 01:25:01 20073
cm64-debugSuccessTue Aug 21 01:27:02 20076
cm64-debug-clear-mallocSuccessTue Feb 1 09:44:23 20053
rs6000-aix
cmSuccessWed Mar 4 01:07:21 20091
cm-debugSuccessMon Mar 2 01:06:10 20091
cm64SuccessWed Mar 4 01:07:21 20091
cm64-debugSuccessTue Mar 3 01:11:44 20091
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:22 2016

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


CMISS Help / Examples / 5 / 52 / 527