( 
 Disk spinning at 1000 rad/s =  9549 RPM  1/2 symmetry model

 Maximum tangential (hoop) stress is at the ID center and about 51ksi.
 Maximum radial stress is at the fillet and about 34ksi.
 Maximum radial deflection at od is about 5.6x10-3 in.

 Units: lbf, lbf s^2/in, in , s, psi
 Osman Buyukisik      osman@fuse.net
)

echo -no
number_of_space_dimensions 2
materi_velocity
materi_stress
materi_strain_total
end_initia

start_arithmetic
   ri     2.0
end_arithmetic
start_arithmetic
   ro     10.0
end_arithmetic
start_arithmetic
   r1     5.5 
end_arithmetic
start_arithmetic
   r2     7.5
end_arithmetic
start_arithmetic
   w1     4.0
end_arithmetic
start_arithmetic
   w2     1.0
end_arithmetic
start_arithmetic
   omega  1000.0 
end_arithmetic
start_arithmetic
   rho    7.4e-4
end_arithmetic
start_arithmetic
   fbody  rho MULTIPLY omega MULTIPLY omega 
end_arithmetic
start_define
   disk_rim  geometry_line 1
end_define
( due to symmetry only 1/2 of the weight )
start_arithmetic
   blade_wgt   0.300
end_arithmetic
start_arithmetic
   number_of_blades 80.0
end_arithmetic
start_arithmetic
   blade_cg    12.0
end_arithmetic
start_arithmetic
   gc  386.4
end_arithmetic
start_arithmetic
   pi2  6.2831853
end_arithmetic
start_arithmetic 
 blade_cf blade_wgt MULTIPLY number_of_blades MULTIPLY omega MULTIPLY omega MULTIPLY blade_cg DIVIDE gc DIVIDE pi2 DIVIDE ro
end_arithmetic

node         0     ri   w1
node         1     ri   0.0
node         2     r1   w1
node         3     r1   0.
node         4     7.203  1.446
node         5     7.203  0.0
node         6     7.563  1.118
node         7     7.563  0.0
node         8     8.035  1.0
node         9     8.035  0.0
node        10     ro   w2
node        11     ro   0.

element 0 -quad4  0 1 2 3
element 1 -quad4  2 3 4 5
element 2 -quad4  4 5 6 7
element 3 -quad4  6 7 8 9
element 4 -quad4  8 9 10 11

geometry_line    0  0.0  0.0   ro   0.0   0.01
disk_rim            ro   w2    ro   0.0   0.01

( circle for the fillet )
geometry_circle  2  8.035 2.0 1.0         0.01

( symmetry plane at Y=0 )
bounda_unknown   0 -geometry_line 0 -vely
bounda_time      0 0.0

( apply the blade rim load )
force_element_edge 0          blade_cf 0. 
force_element_edge_geometry 0 -disk_rim

( apply body load due to the rotational speed  )
force_element_volume            0  1. 0.
force_element_volume_factor     0  0. fbody

( material props  )
group_type 			0  -materi
group_materi_elasti_young 	0  30e6
group_materi_elasti_poisson 0  0.3
group_materi_density        0  rho
group_axisymmetric          0 -yes
group_materi_memory 		0  -updated_without_rotation

post_point 0    ri 0.0
post_point 1    ri 4.0
post_point 2    ro 0.0
post_point 3    ro 1.0

control_mesh_adjust_geometry    0 -geometry_circle  2  -geometry_circle  2
control_mesh_refine_globally    1 -h_refinement
control_mesh_adjust_geometry    1 -geometry_circle  2  -geometry_circle  2
control_mesh_refine_globally    2 -h_refinement
control_mesh_adjust_geometry    2 -geometry_circle  2  -geometry_circle  2

control_timestep                     		5  1. 1.
control_timestep_iterations_automatic           5  5.0e-2  1
(
control_print                                   5 -time_current
)

target_item 1  -post_point_dof 0 -sigzz
target_value 1  51000.0 20000.0
target_item 2  -post_point_dof 2 -sigxx
target_value 2  13000.0 1000.0
end_data
