( 
 Bimetallic beam under thermal load.
 10x0.05x1 one metal and 10x0.05x1 another.
 Modeled only 1/2 of the beam, used symmetry at x=5.0
 Using Roark & Young artcl 7.2 The linear solution is:
 max deflection: 0.9 in  max stress : 12000 psi
 If the nonlinear solution is used, there is no hand calc
 to check but it turns out to lower both the deflection and stress. 
 And it seems to match other FE program results.
 Osman Buyukisik      osman@fuse.net
)

echo -no
number_of_space_dimensions 2
materi_velocity
materi_strain_total
materi_stress
condif_temperature
end_initia

( comment out the next 3 lines for nonlinear solution )
options_convection  -no
options_inertia     -no
options_mesh  -fixed_in_space -fixed_in_space

geometry_point   1  0.0 0.1               0.01
geometry_line    1  0.0  0.0   0.0  0.1   0.01
geometry_line    2  5.0  0.0   5.0  0.1   0.01

bounda_unknown   0 -geometry_point 1 -vely
bounda_time      0 0.0 0.0 100.00  0.0
bounda_unknown   1 -geometry_line 2 -velx
bounda_time      1 0.0 0.0 100.00  0.0
bounda_unknown   2 -all -temp
bounda_time      2 0.0 400.0 1.0 400.0  100.00  400.0

( material props  )
group_type 			0  -materi -condif
group_materi_elasti_young 	0  1e7
group_materi_elasti_poisson 	0  0.3
group_condif_conductivity       0 1.
group_materi_expansion_linear   0 14.5e-6
group_materi_membrane           0 -yes
( comment out the next line if nonlinear solution is required )
group_materi_memory             0  -updated_without_rotation

group_type 			1  -materi -condif
group_materi_elasti_young 	1  1e7
group_materi_elasti_poisson 	1  0.3
group_condif_conductivity       1 1.
group_materi_expansion_linear   1 2.5e-6
( comment out the next line if nonlinear solution is required )
group_materi_memory             1  -updated_without_rotation
group_materi_membrane           1 -yes


(this is for checking the results at x=5)
post_point   0    5.0 0.000001
post_point   1    5.0 0.0500
post_point   2    5.0 0.099999

(create the the first beam)
control_mesh_macro              0 -rectangle 0 5 3
control_mesh_macro_parameters   0  2.5 0.025  5.0 0.05

(create the the second beam)
control_mesh_macro              1 -rectangle 1 5 3
control_mesh_macro_parameters   1  2.5 0.075  5.0 0.05

( merge these two beams )
control_mesh_merge              2  -yes

( WE REALLY NEED THIS!!  linear quads lock quite badly.  shear lock)
control_mesh_refine_globally    3   -p_refinement

control_timestep                       5  1  1.0
control_timestep_iterations            5  3
(
control_print                          5 -time_current
)

( the following checks are for linear solution only !! )
target_item 1  -post_point_dof 0 -vely
target_value 1 -0.9  0.004
target_item 2  -post_point_dof 0 -sigxx
target_value 2 12000.0 1000.0
target_item 3  -post_point_dof 2 -sigxx
target_value 3 -12000.0 1000.0

end_data
