( 2D contact. Geometrical target. Friction.
  Node 1,3 penetrates geometrical point 1 with 0.01.
  After the contact algorithm has been applied,
  node 1,3 should get coordinate 1.
  This gives sigxx=-0.01.
  Because the friction coefficient is 0.1
  we expect about sigxy=-0.1*0.01=-0.001.
  Full integration is needed, otherwise the hourglass
  mode of the element would be strongly activated.
  Bacause of stick the sigxy is more high then for full slip. )

echo -no
number_of_space_dimensions 2
materi_velocity
materi_strain_total
materi_stress
end_initia

node 1  0.00 0.
node 2  1.01 0.
node 3  0.00 1.
node 4  1.01 1.
element 1  -quad4 1 2 3 4

group_type 0  -materi
group_materi_elasti_young 0  1.0
group_materi_memory 0  -updated

bounda_unknown 0  -ra 1 3 -ra -velx
bounda_time    0  0. 0. 100. 0.
bounda_unknown 1  -ra 1 3 -ra -vely
bounda_time    1  0. 1. 100. 1.

geometry_line 1  1. -1. 1. 2. 1.e-1
contact_friction  0.1
contact_penalty_velocity 100.
contact_geometry  1 -geometry_line 1

print_filter 1  -node_dof -all -velx -vely -sigxx -sigxy

control_timestep             10  1.e-3 2.e-3
control_timestep_iterations  10  10
(
control_print                10  -time_current -node_dof
control_print_history        10  -node_dof 4 -sigxy
)

target_item 1  -node_dof 4 -sigxy
target_value 1  -0.001 5.e-4
end_data
