( 
 2D  beam 200mmx2.5mmx2.5mm
 Both ends are built-in
 2.0 N/mm distributed force along the top
                  4              4
              w  L          2 200
 max defl = ----------- = ---------------- = 0.213  mm
             384 E I       384 30E6 3.255
 40x6 mesh
 force   N
 length  mm  --- actually units could be any consistent set.
 Modeled 1/2 of the length. Symm BC at the middle
 Good results with 2nd order elements
 Osman Buyukisik      osman@fuse.net
)

echo -no
number_of_space_dimensions 2
materi_velocity
materi_stress
materi_strain_total
end_initia

(
options_convection  -no
options_inertia     -no
)
options_mesh  -fixed_in_space -fixed_in_space

(point at lower left edge)
geometry_point   0  0.0 0.0               0.01
(point at middle of right edge)
geometry_point   1  100.0 1.25              0.01
(left edge)
geometry_line    1  0.0  0.0   0.0  2.5   0.01
(right edge)
geometry_line    2  100.0  0.0   100.0  2.5   0.01
(top edge)
geometry_line    3  0.0 2.5  100.0 2.5   0.01

(left edge fixed)
bounda_unknown   0 -geometry_line 1 -vely
bounda_time      0 0.0 0.0 100.00  0.0
bounda_unknown   1 -geometry_line 1 -velx
bounda_time      1 0.0 0.0 100.00  0.0
(right edge symm BC)
bounda_unknown   2 -geometry_line 2 -velx
bounda_time      2 0.0 0.0 100.00  0.0
(pressure force at the top of beam)
(
force_element_edge_geometry 4 -geometry_line 3
force_element_edge          4  0.0   -5.00  
)

force_element_edge_normal_geometry 4 -geometry_line 3
force_element_edge_normal          4  -5.00   

( material props  )
group_type 			            0  -materi
group_materi_elasti_young   0  30.0e6
group_materi_elasti_poisson 0  0.3
group_materi_membrane       0 -yes
group_materi_memory         0 -updated_without_rotation
(thickness)
volume_factor                   2.50
(this is for checking the results at x=0)
post_point   0    0.010 0.010
post_point   1    100.0 2.499

(create the the beam)
control_mesh_macro              0 -rectangle 0 20 6
control_mesh_macro_parameters   0  50. 1.25  100.0 2.50
(better for bending )

control_mesh_refine_globally    5   -p_refinement

control_timestep                      10  1  1.0
control_timestep_iterations           10  3

control_print  10 -post_line_dof

target_item 1  -post_point_dof 1 -vely
target_value 1 -0.213 0.01

end_data
