*PARAMETER
# 
# Test settings
tend  = 0.1  , "Termination time"
disp  = 0.008, "Displacement at termination"
Lg    = 0.05 , "Gauge length"
#
# Material data
rho   = 7800 , "Density"
E     = 210e9, "Young's modulus" 
nu    = 0.3  , "Poisson's ratio"
A     = 500e6, "Yield strength"
B     = 300e6, "Strain hardening parameter"
n     = 0.25 , "Strain hardening parameter"
c     = 0.003, "Strain rate coefficient"
rate0 = 1.0  , "Reference strain rate"
#
# Simulation settings
*UNIT_SYSTEM
SI
*TIME
%tend
*FREQUENCY_CUTOFF
1
ALL, 0, 2.0e-6
#
# Mesh
*INCLUDE
specimen.k
*CHANGE_P-ORDER
ALL, 0, 3
#
# Materials
*MAT_USER_1
"IMPETUS example user material"
1, %rho, %E, %nu
1, 0, 0, %c, %rate0
*FUNCTION
"Flow stress vs effective plastic strain"
1
%A + %B*epsp^%n
*MAT_RIGID
23, %rho
#
# Parts
*PART
"Deformable section"
1, 1
"Rigid part"
2, 23
"Rigid part"
3, 23
#
# Boundary conditions
*BC_MOTION
"Prescribed displacement"
2
P, 2, YZ
D, X, 23, 0.5
*BC_MOTION
"Prescribed displacement"
3
P, 3, YZ
D, X, 23, -0.5
*FUNCTION
"Smoothed displacement vs time"
23
smooth_d(%disp, 0, %tend)
#
# Output
*OUTPUT_SENSOR
1, 1, -%Lg/2, 0, 0
2, 1, %Lg/2, 0, 0
*FUNCTION
"Nominal strain"
11
(dxs(2) - dxs(1))/%Lg
*FUNCTION
"Nominal stress"
12
fxr(3)/(12.5e-3*2.0e-3)
*OUTPUT_USER
"Stress vs strain"
10, "Nominal strain", 11, strain
"Nominal stress", 12, stress
*END
