***********************************************************************73
** COPYRIGHT 2011,2012 ANTON SHTERENLIKHT, UNIVERSITY OF BRISTOL, UK.
** ALL RIGHTS RESERVED.
**
** REDISTRIBUTION AND USE IN SOURCE AND BINARY FORMS, WITH OR
** WITHOUT MODIFICATION, ARE PERMITTED PROVIDED THAT THE FOLLOWING
** CONDITIONS ARE MET:
**  1. REDISTRIBUTIONS OF SOURCE CODE MUST RETAIN THE ABOVE COPYRIGHT
**     NOTICE, THIS LIST OF CONDITIONS AND THE FOLLOWING DISCLAIMER.
**  2. REDISTRIBUTIONS IN BINARY FORM MUST REPRODUCE THE ABOVE
**     COPYRIGHT NOTICE, THIS LIST OF CONDITIONS AND THE FOLLOWING
**     DISCLAIMER IN THE DOCUMENTATION AND/OR OTHER MATERIALS
**     PROVIDED WITH THE DISTRIBUTION.
**
** THIS SOFTWARE IS PROVIDED BY THE AUTHOR ``AS IS'' AND ANY
** EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED
** TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS
** FOR A PARTICULAR PURPOSE ARE DISCLAIMED. IN NO EVENT SHALL
** THE AUTHOR OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT,
** INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
** (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE
** GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS
** INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF LIABILITY,
** WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT (INCLUDING
** NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF
** THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
***********************************************************************73
** 1 FE IN TENSION.
*HEADING
1 FE, ROTATION+TENSION, TO CHECK MATERIAL STRESSES, NO UMAT
*PREPRINT,MODEL=YES,HISTORY=YES
*NODE
*NSET,NSET=NALL
*NSET,NSET=DOWN
*NSET,NSET=UP
*NSET,NSET=NOUT
*NSET,NSET=FIXED
*NSET,NSET=CONSTR
*NSET,NSET=MOVE
*ELEMENT,TYPE=C3D8,ELSET=ALL
************************************************************************
*SOLID SECTION,ELSET=ALL,MATERIAL=ZAPPA
*MATERIAL,NAME=ZAPPA
*DENSITY
7.87D-9
*ELASTIC,TYPE=ISOTROPIC
2.0D5, 0.3
*PLASTIC,HARDENING=ISOTROPIC
** THE DATA IS FROM HARD
 384.59 , 0
 414.665 , 0.0056083
 432.623 , 0.0084694
 448.568 , 0.0127718
 464.232 , 0.0211876
 476.584 , 0.0255166
 489.561 , 0.0340914
 501.66 , 0.0427065
*POROUS METAL PLASTICITY, RELATIVE DENSITY=0.99
*VOID NUCLEATION
** 1 - EPSILON_N
** 2 - S_N
** 3 - F_N
************************************************************************
*BOUNDARY
FIXED,1,3
CONSTR,2,3
************************************************************************
*STEP,NLGEOM=YES
*STATIC
*AMPLITUDE,NAME=TENSILE,DEFINITION=TABULAR,TIME=TOTAL TIME,VALUE=RELATIVE,SMOOTH=.05
0.0,0.0, 1.0,1.0
************************************************************************
*BOUNDARY,TYPE=DISPLACEMENT,AMPLITUDE=TENSILE
MOVE,2,2,2.0
************************************************************************
*OUTPUT,HISTORY,FREQUENCY=1
*ELEMENT OUTPUT,ELSET=ALL
LE,PE,ER,ERV,PEEQ
S,MISES,PRESS
ELEN
VVF,VVFG,VVFN
*NODE OUTPUT,NSET=NOUT
U2
*OUTPUT,FIELD,NUMBER INTERVAL=10
*NODE OUTPUT
U
**
** OUTPUT TO DAT FILE
**
*EL PRINT,ELSET=ALL,SUMMARY=NO,TOTALS=NO,POSITION=CENTROIDAL
LE
*EL PRINT,ELSET=ALL,SUMMARY=NO,TOTALS=NO,POSITION=CENTROIDAL
PE
*EL PRINT,ELSET=ALL,SUMMARY=NO,TOTALS=NO,POSITION=CENTROIDAL
ER
*EL PRINT,ELSET=ALL,SUMMARY=NO,TOTALS=NO,POSITION=CENTROIDAL
S
*EL PRINT,ELSET=ALL,SUMMARY=NO,TOTALS=NO,POSITION=CENTROIDAL
MISES,PRESS
*EL PRINT,ELSET=ALL,SUMMARY=NO,TOTALS=NO,POSITION=CENTROIDAL
ELEN
*EL PRINT,ELSET=ALL,SUMMARY=NO,TOTALS=NO,POSITION=CENTROIDAL
ENER
*EL PRINT,ELSET=ALL,SUMMARY=NO,TOTALS=NO,POSITION=CENTROIDAL
VVF,VVFG,VVFN
*NODE PRINT,NSET=NOUT
U2
**
** OUTPUT TO FIL FILE
**
*EL FILE,ELSET=ALL,POSITION=CENTROIDAL
S,SINV,LE,EE,PE,VVF,VVFG,VVFN
*EL FILE,ELSET=ALL,POSITION=CENTROIDAL
ELEN
*EL FILE,ELSET=ALL,POSITION=CENTROIDAL
ENER
************************************************************************
*END STEP