*PARAMETER
tend = 0.05  , "Termination time"
#
L    = 0.4   , "Length of specimen"
t    = 0.05  , "Thickness and width of specimen"
dx   = 1.0e-2, "Element size"
#
R    = 0.03  , "Radius of punch and supports"
disp = 1e-3  , "Displacement of punch at termination"
#
*UNIT_SYSTEM
SI
*TIME
%tend
*COMPONENT_BOX
"Specimen"
1, 1, %L/%dx, %t/%dx, %t/%dx
-%L/2, -%t/2, 0, %L/2, %t/2, %t
*COMPONENT_CYLINDER
"Support X-"
2, 2, 5, 3
-%L/2+%R, -%t, -%R, -%L/2+%R, %t, -%R, %R
*COMPONENT_CYLINDER
"Support X+"
2, 2, 5, 3
%L/2-%R, -%t, -%R, %L/2-%R, %t, -%R, %R
*COMPONENT_CYLINDER
"Punch X-"
3, 3, 5, 3
-%L/4+%R, -%t, %t+%R, -%L/4+%R, %t, %t+%R, %R
*COMPONENT_CYLINDER
"Punch X+"
3, 3, 5, 3
%L/4-%R, -%t, %t+%R, %L/4-%R, %t, %t+%R, %R
*CHANGE_P-ORDER
P, 2, 2
P, 3, 2
*PARAMETER
E        = 30e9
nu       = 0.2
#
rho      = 2300
G        = %E/(2*(1+%nu))
A        = 0.10
B        = 1.0
n        = 1.0
f_c      = 50e6
T        = 5e6
epsf_min = 2e-3
p_c      = 100e6
p_l      = 500e6
xmu_l    = 0.05
D1       = 0.02
D2       = 0.5
K        = %E/(3*(1-2*%nu))
*MAT_USER_2
"IMPETUS example user material"
1, %rho, %G
%A, %B, %n, %f_c, %T, %epsf_min, %p_c, %p_l
%xmu_l, %D1, %D2, %K
*MAT_RIGID
23, 7800
*PART
1, 1
2, 23, 0, 0, 30 
3, 23, 0, 0, 30
*CONTACT
1
ALL, 0, ALL, 0, 0.2
*BC_MOTION
2
P, 2, XYZ
*BC_MOTION
3
P, 3, XY
D, Z, 3, -1
*FUNCTION
"Displacement"
3, 0, 0, 0, time, disp
smooth_d(%disp, 0, %tend)
*END
