*HEADING
ADAPTIVE MESHING EXAMPLE
SECTION ROLLING
Units - N, m, seconds
*RESTART, WRITE,NUMBER=10
*SYSTEM
  -0.0949794,          0.,          0.,  -0.0949794,          1.,          0.
  -0.0949794,          0.,          1.
*INCLUDE, INPUT=ale_rolling_sectionnode.inp
*INCLUDE, INPUT=ale_rolling_sectionelem.inp
*INCLUDE, INPUT=ale_rolling_sectionnelset.inp
** SECTION: METAL
*SOLID SECTION, ELSET=METAL, MATERIAL=STEEL
1., 
*MATERIAL, NAME=STEEL
*DENSITY
       7833.,
*ELASTIC
    2.12E+11,       0.281
*PLASTIC
       8e+07,          0.,   
    2.35e+08,         0.6
*ELEMENT, TYPE=ROTARYI, ELSET=ROTI
50000,10117
*ROTARY INERTIA, ELSET=ROTI
1.E-4,11.05,1.E-4
*ELEMENT, TYPE=MASS, ELSET=MASS
50001,10117
50002,10116
*MASS, ELSET=MASS
10.,
** INITIAL CONDITION: VELOCITY
*INITIAL CONDITIONS, TYPE=VELOCITY
ROLNODES, 1, -4.187
ROLNODES, 2, 0.
ROLNODES, 3, 0.
** STEP: STEP-1
** 
*SURFACE, TYPE=REVOLUTION, NAME=VROLSURF
-0.1800, 0., 0.6221, -0.1800, 100., 0.6221
START,          0.2,          0.2
 LINE,       0.4354,          0.2
 LINE,      0.45182,     0.035821
CIRCL,      0.45182,    -0.035821,     0.093604,   6.9389E-18
 LINE,       0.4354,         -0.2
 LINE,          0.2,         -0.2
*SURFACE, TYPE=REVOLUTION, NAME=HROLSURF
-0.1800, 0.760, 0., -0.180, 0.760, 100.
START,        0.547,        0.137
 LINE,      0.71051,      0.12261
CIRCL,        0.747,     0.082765,        0.707,     0.082765
 LINE,        0.747,    -0.082765
CIRCL,      0.71051,     -0.12261,        0.707,    -0.082765
 LINE,        0.547,       -0.137
*INCLUDE, INPUT=ale_rolling_sectionsurf.inp
*RIGID BODY, REF NODE=10116, ANALYTICAL SURFACE =HROLSURF
*RIGID BODY, REF NODE=10117, ANALYTICAL SURFACE =VROLSURF
*SURFACE INTERACTION, NAME=PRO-1
*FRICTION
0.3,
*STEP
*DYNAMIC, EXPLICIT
, 0.227
*BULK VISCOSITY
0.06, 1.2
** BC: H_ROL_HLD TYPE: DISPLACEMENT
*BOUNDARY, OP=NEW
HROLREF, 1, 1
HROLREF, 2, 2
HROLREF, 3, 3
HROLREF, 4, 4
HROLREF, 5, 5
** BC: H_ROLL_V TYPE: VELOCITY
*BOUNDARY, OP=NEW, TYPE=VELOCITY
HROLREF, 6, 6,-5.
** BC: V_ROL_HLD TYPE: DISPLACEMENT
*BOUNDARY, OP=NEW
VROLREF, 1, 1
VROLREF, 2, 2
VROLREF, 3, 3
VROLREF, 4, 4
VROLREF, 5, 5
VROLREF, 6, 6
** BC: Y_SYMM TYPE: TYPE_CONSTRAINT
*BOUNDARY, OP=NEW
YSYM, YSYMM
** BC: Z_SYMM TYPE: TYPE_CONSTRAINT
*BOUNDARY, OP=NEW
ZSYM, ZSYMM
*CONTACT 
*CONTACT INCLUSIONS
,
*CONTACT PROPERTY ASSIGNMENT
,, PRO-1
*FILE OUTPUT, NUMBER INTERVAL=1, TIME MARKS=YES
*NODE FILE, NSET=HROLREF
RF, 
*EL FILE, ELSET=QA_TEST
PEEQ,
*****
*NSET, NSET=QA_TEST_REFN
HROLREF,VROLREF
*ELSET, ELSET=QA_TEST
7659, 7660, 7664, 7665, 7666, 7683, 7726, 7718, 7719, 7715
*NSET, NSET=QA_TEST
ZSYM,
*OUTPUT, FIELD, TIME MARKS=YES, NUMBER INTERVAL=1
*ELEMENT OUTPUT, ELSET=QA_TEST

PEEQ,
*NODE OUTPUT, NSET=QA_TEST
U,
*NODE OUTPUT, NSET=QA_TEST_REFN
RF,
**** Output below for generating the figures.
*NODE OUTPUT
U,
*ELEMENT OUTPUT

PEEQ,
*OUTPUT, HIST, FREQ=9999
*ENERGY OUTPUT, VAR=PRESELECT
******
*** Following history output for generating figures.
*OUTPUT,HISTORY,TIME INTERVAL=0.0026
*NODE OUTPUT,NSET=HROLREF
RF2,RF3,RM2,RM3
*NODE OUTPUT,NSET=VROLREF
RF2,RF3,RM3,RM2
***
**MASSSCALING:WHOLEMODEL1
**ROLLING
*VARIABLE MASS SCALING,TYPE=BELOW MIN,FREQUENCY =50,ELSET=METAL,DT=2.E-5
*END STEP