( 
 Rotating circular plate - constant thickness
 R = 24 mm
 thick = 4 mm
 density = 2.7e-8  ton/mm^3
 Rotational speed = 150000 RPM = 15708 rad/s
 Distributed force is = rho*omega^2*r
 This force is linearly dependent on r (which is the X coord).
 So rho*omega^2 = 6.662  <-- force_element_volume X coeff

 sigxx at x=0 = sigzz = 1/8 rho*omega^2*(3+nu)*R^2 
                      = 1583 MPa
 Units :
 length: mm
 Force : N
 Mass  : ton
 time  : sec
 Stress: MPa

 Osman 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

geometry_line    1  0.0  -10.0   0.0  50.0   0.1

bounda_unknown   0 -geometry_line 1 -velx
bounda_time      0 0.0 0.0 100.00  0.0

(Aluminum material props  )
group_type 			0  -materi
group_materi_elasti_young 	0  70000.00
group_materi_elasti_poisson 	0  0.3
group_materi_density         	0 2.7e-8
group_materi_memory 		0  -updated_without_rotation
group_axisymmetric 0 -yes

(due to the model starting from x=0, we need -minimal )
group_integration_points        0 -minimal

( 
  Body force 150,000 RPM --> rho*omega^2 = 6.662
  force/volume = 0.0 + 6.662*X 
)
force_element_volume            1  1. 0.
force_element_volume_factor     1  0. 6.662

(this is for checking the results at x=0)
post_line 1  0. 0. 0. 4.

control_mesh_macro              0 -rectangle 0 13 5
control_mesh_macro_parameters   0  12.000 2.000  24.00 4.00 

control_timestep                5  1. 1.

target_item 0 -post_line_dof 1 -sigxx

target_value 0 1583.0 10.0
end_data
