( Homepage example: Pressed rubber cylinder.
  Ref. Sussman T.,Bathe K.J. Computers&Structures Vol 26, No:1&2 1987
  rubber cylinder (0.4m diam, plane strain) compressed by 0.2m
  force is 1.40 +- 0.05 MN  
  This is a 1/4 model.
  Mooney-Rivlin parameters are given.
  C1=0.293 C2=0.177
  Small strain ELastic Modulus = 2.82, Poisson's ratio 0.49967

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

echo -no
number_of_space_dimensions 2
materi_velocity
materi_displacement
materi_stress
materi_strain_elasti
materi_strain_total
end_initia

node  0  0. -0.2
node  1  0. -0.1
node  2  0.  0.0
node  3  .1  0.0
node  4  .2  0.0
node  5  .1414 -0.14142
node  6  .0805 -0.06667

element 0 -quad4 0 5 1 6
element 1 -quad4 5 4 6 3
element 2 -quad4 1 6 2 3

group_type 0  -materi
group_materi_memory 0  -total
group_materi_elasti_young 0  2.82
group_materi_elasti_poisson 0  0.49967
group_materi_hyper_mooney_rivlin 0 .293 .177

(left edge)
geometry_line 0  0.0 -0.2 0.0 0.0   1.e-3

(upper edge)
geometry_line 1  0.0  0.0 0.2 0.0  1.e-3

(contact line)
geometry_line 2  0.0 -0.2 1.0 -0.2 1.e-2

(circle edge)
geometry_circle 0  0.0 0.0 0.2 1.e-3

(all edges)
geometry_set 0  -geometry_line 0 -geometry_line 1 -geometry_circle 0

bounda_unknown 0  -geometry_line 0 -velx
bounda_time    0  0. 0. 100. 0.
bounda_unknown 1  -geometry_line 1 -vely
bounda_time    1  0.0 -0.1 100.0 -0.1

contact_penalty_velocity  1.e5
contact_geometry  1 -geometry_line 2

post_node 1 -node_rhside -sum -geometry_line 1

control_mesh_refine_globally           7  -h_refinement
control_mesh_refine_globally_geometry  7  -geometry_set 0
control_mesh_refine_globally           8  -h_refinement
control_mesh_refine_globally_geometry  8  -geometry_set 0

control_timestep                      10  0.01 1.
control_timestep_iterations           10  10
(
control_print                         10  -time_current 
control_print_history                 10  -post_node_result 1 -vely
)

target_item 1  -post_node_result 1 -vely 
target_value 1  1.41727 0.02
end_data

