(
Two axisymm elements with interference  2D
inner radius shaft=1.0
outer radius shaft=2.0
inner radius hub=1.998
outer radius hub=3.0

E=1.e7

The contact algorithm displaces the nodes
to get rid of the penetration.
Node 2 should go about 1.e-3 to the left.

Test supplied by: Osman F Buyukisik, osman@fuse.net
)

echo -no
number_of_space_dimensions 2
materi_velocity
materi_strain_total
materi_stress
end_initia

node 1  1.      0.0
node 2  2.000   0.0
node 3  1.998   0.0
node 4  3.0     0.0
node 5  1.      1.0
node 6  2.      1.0
node 7  1.998   1.0
node 8  3.0     1.0

element 1  -quad4 1 2 5 6
element 2  -quad4 3 4 7 8

node_boundary   2  -yes   
node_boundary   6  -yes   

bounda_unknown 0 -ra 1 2 3 4 -ra  -vely
bounda_time    0 0.0 0.0 100.0 0.0

contact_penalty_velocity  1.0e12

group_type 0  -materi
group_materi_elasti_young 0  1.0e7
group_materi_elasti_poisson 0 0.3
group_materi_memory 0  -updated
group_integration_points 0  -maximal
group_axisymmetric 0 -yes

control_timestep    10  1. 1.
(control_print      10 -node_dof )

target_item    0  -node_dof 2 -velx
target_value   0  -1.e-3 1.e-4
end_data
