( 
 A cyldrical structure loaded with hydrostatic pressure in +z dir
 The cylinder is clamped at the bottom. 
 Does not have end cover plates.
 Osman Buyukisik   osman@fuse.net
)

echo -no
number_of_space_dimensions 3
materi_velocity
materi_strain_total
materi_stress
end_initia

( define some macros  )
start_define
  lower_edge geometry_quadrilateral 0
end_define
start_define
  upper_edge geometry_quadrilateral 1
end_define
start_define
  inner_edge geometry_cylinder 0
end_define
start_define
  outer_edge geometry_cylinder 1
end_define

( define geometry )
lower_edge   -10.0 -10.0 0.0   10. -10.0 0.0  -10.0 10. 0.0  10.0 10. 0.0 1.e-3
upper_edge   -10.0 -10.0 1.0   10. -10.0 1.0  -10.0 10. 1.0  10.0 10. 1.0 1.e-3
inner_edge   0. 0. 0.   0. 0. 1.   0.90   1.e-3
outer_edge   0. 0. 0.   0. 0. 1.   1.00   1.e-3

(Steel AISI-304 material properties, SI units)
group_type 0  -materi 
group_integration_points 0  -minimal 
group_materi_density 0  8030.0
group_materi_elasti_young   0  1.93e11
group_materi_elasti_poisson 0  0.29
group_materi_memory 0  -updated_without_rotation

(fix the lower edge )
bounda_unknown 0  -lower_edge -velx -vely -velz
bounda_time    0  0. 0. 1.e6 0.

( apply the pressure to inner edge )
force_element_edge_normal  5 -1e6
force_element_edge_normal_factor  5 3 -1
force_element_edge_normal_geometry 5 -inner_edge

(here we create the mesh)
control_mesh_macro             10  -cylinder_hollow 0 20 5 20
control_mesh_macro_parameters  10  0. 0. 0.  0. 0. 1.  0.95  0.1  0. 360.

(just to check stresses )
post_point     1  1.0 0.0 0.0

control_timestep               30  1 1.
control_timestep_iterations    30  3

end_data

