the code is still in draft form, but it works. The displayed strains/stresses are elemental values.

srcs/FEM/my_geo.geo

 Lx = 5;
 Ly = 2;
-nx = 5;
-ny = 1;
+nx = 50;
+ny = 20;
 
 Point(1) = {0, 0, 0, 1.0};
 Point(2) = {Lx, 0, 0, 1.0};
@@ -21,12 +21,18 @@ Recombine Surface {1}; // quads instead of triangles
 
 Physical Curve("left_edge", 5) = {4};
 Physical Surface("domain", 6) = {1};
-Physical Curve("top_edge",7) = {3};
+Physical Curve("top_edge", 7) = {3};
+Physical Curve("right_edge", 8) = {2};
 
 // additional parameters given to the solver
-SetNumber("Boundary Conditions/left_edge/ux", 0.);
-SetNumber("Boundary Conditions/left_edge/uy", 0.);
+SetNumber("Boundary Conditions/left_edge/ux", 0.); // HERE YOU DO NOT HAVE TO IMPOSE BOTH ux and uy simultaneously ! (permet aussi de simuler appuis à roulettes)
+//SetNumber("Boundary Conditions/left_edge/uy", 0.);
 SetNumber("Materials/domain/Young", 210e3);
 SetNumber("Materials/domain/Poisson", 0.3);
-SetNumber("Boundary Conditions/top_edge/tx", 0.);
-SetNumber("Boundary Conditions/top_edge/ty", -100.);
\ No newline at end of file
+SetNumber("Materials/domain/rho",7800); //volumic mass of acier
+SetNumber("Boundary Conditions/top_edge/tx", 0.); // ALWAYS NEED TO IMPOSE BOTH tx AND ty ON A GIVEN EDGE (realiste, OK) !
+SetNumber("Boundary Conditions/top_edge/ty", 0.); //set to some non-zero value to induce vertical deflection
+SetNumber("Boundary Conditions/right_edge/tx", 21e3); // for simple tension conditions
+SetNumber("Boundary Conditions/right_edge/ty", 0.);
+SetNumber("Volumic Forces/domain/bx",0.);
+SetNumber("Volumic Forces/domain/by",0.); //set to -9.81 for gravity