This example tests the pole-zero stress/strain relationship. A unit cube is first compressed in the x-direction (fiber direction). The calculated strain and stress at a Gauss point is found. Note that it is the 2nd Piola-Kirchhoff stress before hysdrostatic pressure that is found. Using the extracted strain, a Perl routine calculates the stress for this strain to compare with the cmiss calculated stress. The three values are written in as three coloumns of a .dat file. There might be a difference in the last decimals between the two stresses as the extracted strain value has fewer decimals than what is actually used to calculate the cmiss stress. The cube is incrementally deformed from a very compressed state along the x-direction until it reaches a very elongated state (stretched beyond the fibre pole). At each increment the procedure for calculating the strain and stresses is carried out and written into the file.

Afterwards a simple shear deformation is applied to the unit cube. The top surface (z=1) is incrementally moved from a negative shearing angle (beyond the shear pole) to a positive shear angle. The shear strains and stresses are obtained as described before and written to a separate .dat file. A Gnuplot is used for visualization of the results.

Created by Espen Remme November 2003.