( Home page example: Elasto-plastic buckling of heated hollow cylinder.
 
 A hollow cylinder is heated from inside, 
 and convects heat from the outer edge to the environment.
 The cylinder expands because of the heat.
 The cylinder is clamped at the bottom and top, 
 so that it buckles because of the thermal expansion.

 An updated Lagrange large deformation theory is 
 used to model the large strains and rotations in the cylinder.

 Plasticity shows up because of the high deformations.
)

echo -no
number_of_space_dimensions 3
condif_temperature
materi_velocity
materi_strain_total
materi_stress
materi_plasti_kappa
end_initia

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

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 -condif
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
group_materi_plasti_vonmises 0  24.3e6
group_materi_expansion_linear 0  1.78e-5
group_condif_conductivity     0  16.3   
(
group_condif_capacity         0  502.0   
)

condif_convection  0   10.0 0.0
condif_convection_geometry 0  -outer_edge

bounda_unknown 0  -lower_edge -velx -vely -velz
bounda_time    0  0. 0. 1.e6 0.
bounda_unknown 1  -upper_edge -velx -vely -velz
bounda_time    1  0. 0. 1.e6 0.
bounda_unknown 2  -inner_edge -temp
bounda_time    2  0. 0. 1. 2.e3

post_point     1  1.0 0.0 0.0

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.

control_timestep               30  0.1 1.
control_timestep_iterations    30  5
(
control_print                  30  -time_current -post_node_rhside_ratio
)

target_item 0  -post_point_dof 1 -kap
target_value 0   0.0747218 0.01
end_data