*heading
normal impact of a sphere with a rigid plane
**
*preprint, model=yes, history=yes
**
**NODES
**
*node
1, -0.25, 0.0, 0.0
2, 0.25, 0.0, 0.0
3, 0.5, 0.5, -0.1000001
4, -0.5, 0.5, -0.1000001
5, -0.5, -0.5, -0.1000001
6, 0.5, -0.5, -0.1000001
7, 0.0, 0.0, 0.0
**
**NODE SETS
**
*Nset, nset=Al_Alloy
1,
*Nset, nset=Mg_Alloy
2,
*Nset, nset=dem
1,2
*nset,nset=refnode
7,
**
**ELEMENTS and ELEMENT SETS
**
*element, type=pd3d, elset=Al_alloy
1, 1
*element, type=pd3d, elset=Mg_alloy
2, 2
*element, type=s4r, elset=rigid_plane
3, 3, 4, 5, 6
*Elset, elset=dem
1, 2
**
**SECTIONS
**
*shell section, elset=rigid_plane, material=mat-1, offset=SPOS
0.1, 5
*rigidbody, refnode=refnode, elset=rigid_plane
*discrete section, elset=dem,shape=sphere,density=distDens
distRad
*distribution table, name=distDensTable
density
*distribution, name=distDens, table=distDensTable, location=element
 ,1.0
1, 2699
2, 1800
*distribution table, name=distRadTable
length
*distribution, name=distRad, table=distRadTable, location=element
 ,1.0
1, 0.1
2, 0.1
**
**SURFACES
**
*surface, type=element, name=Al_alloy_surf
Al_alloy,
*surface, type=element, name=Mg_alloy_surf
Mg_Alloy,
*surface, type=element, name=rigid_plane_surf
rigid_plane, SPOS
**
**SURFACE INTERACTION 
**
*surface interaction, name=Al_alloy_int
*friction
0.0
*surface behavior, pressure-overclosure=tabular
 0, 0
3.243362e+001, 0.000001
9.173612e+001, 0.000002
1.685300e+002, 0.000003
2.594689e+002, 0.000004
3.626189e+002, 0.000005
4.766749e+002, 0.000006
6.006790e+002, 0.000007
7.338890e+002, 0.000008
8.757077e+002, 0.000009
1.025641e+003, 0.000010
1.183272e+003, 0.000011
1.348240e+003, 0.000012
1.520234e+003, 0.000013
1.698977e+003, 0.000014
1.884223e+003, 0.000015
2.075751e+003, 0.000016
2.273363e+003, 0.000017
2.476875e+003, 0.000018
2.686122e+003, 0.000019
2.900951e+003, 0.000020
3.121220e+003, 0.000021
3.346797e+003, 0.000022
3.577562e+003, 0.000023
3.813399e+003, 0.000024
4.054202e+003, 0.000025
4.299871e+003, 0.000026
4.550310e+003, 0.000027
4.805432e+003, 0.000028
5.065151e+003, 0.000029
5.329387e+003, 0.000030
5.598065e+003, 0.000031
5.871112e+003, 0.000032
6.148459e+003, 0.000033
6.430041e+003, 0.000034
6.715795e+003, 0.000035
7.005661e+003, 0.000036
7.299582e+003, 0.000037
7.597501e+003, 0.000038
7.899367e+003, 0.000039
8.205128e+003, 0.000040
8.514736e+003, 0.000041
8.828142e+003, 0.000042
9.145302e+003, 0.000043
9.466172e+003, 0.000044
9.790709e+003, 0.000045
1.011887e+004, 0.000046
1.045062e+004, 0.000047
1.078592e+004, 0.000048
1.112473e+004, 0.000049
1.146702e+004, 0.000050
1.181274e+004, 0.000051
1.216187e+004, 0.000052
1.251438e+004, 0.000053
1.287022e+004, 0.000054
1.322938e+004, 0.000055
1.359181e+004, 0.000056
1.395750e+004, 0.000057
1.432641e+004, 0.000058
1.469851e+004, 0.000059
1.507378e+004, 0.000060
1.545219e+004, 0.000061
1.583372e+004, 0.000062
1.621833e+004, 0.000063
1.660601e+004, 0.000064
1.699673e+004, 0.000065
*contact damping, definition=damping coefficient
0.0
*surface interaction, name=Mg_alloy_int
*friction
0.0
*surface behavior, pressure-overclosure=tabular
0, 0
1.921992e+001, 0.000001
5.436215e+001, 0.000002
9.986964e+001, 0.000003
1.537594e+002, 0.000004
2.148853e+002, 0.000005
2.824740e+002, 0.000006
3.559579e+002, 0.000007
4.348972e+002, 0.000008
5.189379e+002, 0.000009
6.077873e+002, 0.000010
7.011979e+002, 0.000011
7.989571e+002, 0.000012
9.008793e+002, 0.000013
1.006801e+003, 0.000014
1.116577e+003, 0.000015
1.230075e+003, 0.000016
1.347178e+003, 0.000017
1.467778e+003, 0.000018
1.591776e+003, 0.000019
1.719082e+003, 0.000020
1.849612e+003, 0.000021
1.983287e+003, 0.000022
2.120037e+003, 0.000023
2.259792e+003, 0.000024
2.402490e+003, 0.000025
2.548072e+003, 0.000026
2.696480e+003, 0.000027
2.847663e+003, 0.000028
3.001571e+003, 0.000029
3.158155e+003, 0.000030
3.317372e+003, 0.000031
3.479177e+003, 0.000032
3.643531e+003, 0.000033
3.810395e+003, 0.000034
3.979731e+003, 0.000035
4.151503e+003, 0.000036
4.325678e+003, 0.000037
4.502223e+003, 0.000038
4.681106e+003, 0.000039
4.862298e+003, 0.000040
5.045769e+003, 0.000041
5.231492e+003, 0.000042
5.419438e+003, 0.000043
5.609584e+003, 0.000044
5.801902e+003, 0.000045
5.996369e+003, 0.000046
6.192962e+003, 0.000047
6.391657e+003, 0.000048
6.592433e+003, 0.000049
6.795268e+003, 0.000050
7.000142e+003, 0.000051
7.207035e+003, 0.000052
7.415926e+003, 0.000053
7.626798e+003, 0.000054
7.839631e+003, 0.000055
8.054408e+003, 0.000056
8.271112e+003, 0.000057
8.489724e+003, 0.000058
8.710230e+003, 0.000059
8.932612e+003, 0.000060
9.156855e+003, 0.000061
9.382944e+003, 0.000062
9.610864e+003, 0.000063
9.840600e+003, 0.000064
1.007214e+004, 0.000065
*contact damping, definition=damping coefficient
0.0
**
**INITIAL CONDITIONS
**
*initial conditions, type=velocity
1, 3, -0.2
2, 3, -0.2
**
** MATERIAL
**
*material,name=mat-1
*elastic
2e11,0.3
*density
8000
**
**BOUNDARY CONDITIONS
**
*boundary
7, 1, 6
**
**STEP
**
*step
*dynamic, explicit, direct
1e-6, 8e-4
**
**CONTACT DEFINITIONS
**
*contact
*contact inclusions
Al_alloy_surf, rigid_plane_surf
Mg_alloy_surf, rigid_plane_surf
*contact property assignment
Al_alloy_surf, rigid_plane_surf, Al_alloy_int
Mg_alloy_surf, rigid_plane_surf, Mg_alloy_int
*contact controls assignment, rotational terms=structural
**
**HISTORY OUTPUT
**
*output, history, freq=1
*node output, nset=Al_alloy
U, V, VR
*node output, nset=Mg_alloy
U, V, VR
*energy output, variable=all
**
**FIELD OUTPUT
**
*output, field, timeinterval=1e-6
*node output, nset=dem
U, V
*contact output
CFORCE
*end step