Example 528: Shear of unit cube in xz-plane

A 3D unit cube of incompressible myocardium (pole-zero behaviour) is defined to vary linearly in the X and Y directions anc cubic Hermite in Z and is sheared in the XZ plane (Z=1 face is sheared in -X direction). The first principal stress is fitted to a bilinear-cubic field and exported to CMGUI on the deformed mesh (which has been transferred into the geometry).


The comfile run by this example is as follows:

#Example_528  Shear of unit cube in xz-plane; first principal stress fitted 
#             and exported to cmgui

#            Note that principal stress is fitted and exported 
#              to frontend.

fem
fem define coordinate;r;sheared_cube;example #Defines 3D rectangular
#                              !  cartesian coordinates.
fem define node;r;;example         #Defines 8 nodes in a unit cube.
fem define base;r;;example         #Defines 6 basis functions: three 3D 
#                              !  and three 2D.
fem define element;r;;example      #Defines 1 element: x bilinear-cubic; 
#                              !  y,z trilinear.
fem define window;r;;example on 3  #Applied shear will be in xz plane
fem draw lines
fem define fibre;d             #  default fibre angle
fem define element;r;;example fibre
fem define equation;r;;example lock  #Defines 3D finite elasticity.
fem define material;r;;example     #Pole-zero.
fem define initial;r;;example      #Defines simple shear in xz plane.
fem define solve;r;;example                 
fem solve
fem draw line deformed dotted
fem update gauss stress principal #Principal stresses in YG gauss point
#                              !  array.
fem list gauss yg
fem update geometry from solution #Transfer deformed coordinates to 
#                              !  geometry. 
fem define field;r;;example        #Sets up one field component.
fem define element;r;;example field
fem define fit;r;;example gauss    #Gauss point fit.
fem fit                        #Fit first principal stress.
fem export element;sheared_cube #For use with cmgui.
fem export node;sheared_cube   #Note: nodes are in deformed position.

Additional testing commands:

fem list node

Files used by this example are:

Name                     Modified     Size

example_528.com 10-Apr-2000 1.6k sheared_cube.exelem 16-Dec-2002 5.7k sheared_cube.exnode 16-Dec-2002 2.1k sheared_cube.iod 10-Apr-2000 70k sheared_cube.ipbase 10-Apr-2000 7.8k sheared_cube.ipcoor 10-Apr-2000 570 sheared_cube.ipelem 10-Apr-2000 466 sheared_cube.ipelfb 09-Jul-2001 237 sheared_cube.ipelfd 10-Apr-2000 238 sheared_cube.ipequa 02-May-2004 2.0k sheared_cube.ipfibr 30-Jan-2001 1.0k sheared_cube.ipfiel 10-Apr-2000 1.4k sheared_cube.ipfit 13-Apr-2007 1.6k sheared_cube.ipinit 05-Oct-2005 1.5k sheared_cube.ipmate 05-Dec-2002 5.6k sheared_cube.ipnode 10-Apr-2000 2.3k sheared_cube.ipsolv 16-Aug-2010 2.3k sheared_cube.ipsolv.old 13-Apr-2007 2.1k sheared_cube.ipwind 10-Apr-2000 289 temp.iod 10-Apr-2000 56k test_output.com 10-Apr-2000 14

Download the entire example:

Name                      Modified     Size

examples_5_52_528.tar.gz 17-Aug-2010 19k

Testing status by version:

StatusTestedReal time (s)
i686-linux
cmSuccessSun Mar 6 00:01:58 20160
cm-debugSuccessSat Mar 5 00:02:16 20161
mips-irix
cmSuccessSun Aug 19 01:28:28 20073
cm-debugSuccessWed Aug 15 01:25:45 20078
cm-debug-clear-mallocSuccessSat Aug 18 01:27:55 200712
cm-debug-clear-malloc7SuccessMon Aug 20 01:26:17 200711
cm64SuccessSun Aug 19 01:29:33 20074
cm64-debugSuccessTue Aug 21 01:23:16 20079
cm64-debug-clear-mallocSuccessTue Feb 1 09:44:27 20052
rs6000-aix
cmSuccessWed Mar 4 01:07:22 20091
cm-debugSuccessMon Mar 2 01:06:58 20092
cm64SuccessWed Mar 4 01:07:22 20091
cm64-debugSuccessTue Mar 3 01:11:56 20092
x86_64-linux
cmSuccessSun Mar 6 00:01:03 20161
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: Mon Aug 16 11:19:13 2010


CMISS Help / Examples / 5 / 52 / 528