(
Single axisymmetric element test
inner radius=10.0
outer radius=11.0
length = 1
innner surface has 1000 psi pressure
E=1.e7, nu=0.3
from J. Akin's 1995 text p.131-132 solution is:
ur = 9.9642E-3   9.5309E-3  deflections at ID and OD resp.
e_r=-4.33E-4  ID  strains radial and tangential
e_t=9.964E-4  ID
S_t=10910  at ID  tangential stress

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

options_mesh -fixed_in_space -fixed_in_space
options_convection -no

node 1  10. 0.
node 2  11.0 0.
node 3  10.0 1.
node 4  11.0 1.
element 1  -quad4 1 2 3 4

group_type 0  -materi
group_materi_elasti_young 0  1.0e7
group_materi_memory 0  -updated_without_rotation
group_axisymmetric 0 -yes

(inner radius)
geometry_line  1  10. 0. 10. 1. 1.e-4
(lower edge)
geometry_line  2  10. 0. 11. 0. 1.e-4

bounda_unknown 2  -geometry_line 2 -vely 
bounda_time    2  0. 0. 100. 0.

force_element_edge 0  1000. 0.
force_element_edge_geometry 0  -geometry_line 1
force_element_edge_time 0  0.0 1.0 1.0 1.0

post_point 1  10. 0.
print_filter 1  -post_point_dof -all  -velx -sigxx -sigzz

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

target_item 1  -post_point_dof 1 -sigzz
target_value 1  10910.0 2000.0
end_data
