From 36ca5ef693f39ba75fb396a6ac89f5c7a939c436 Mon Sep 17 00:00:00 2001
From: Louis Denis <louis.denis@student.uliege.be>
Date: Tue, 10 May 2022 14:20:54 +0200
Subject: [PATCH] large rotation

---
 srcs/FEM/large_rotation_validation.geo | 64 ++++++++++++++++++++++++++
 1 file changed, 64 insertions(+)
 create mode 100644 srcs/FEM/large_rotation_validation.geo

diff --git a/srcs/FEM/large_rotation_validation.geo b/srcs/FEM/large_rotation_validation.geo
new file mode 100644
index 0000000..d10805c
--- /dev/null
+++ b/srcs/FEM/large_rotation_validation.geo
@@ -0,0 +1,64 @@
+
+h = 1;
+H = 10;
+
+n = 20;
+
+Point(1) = {0, 0, 0, 0.5};
+Point(2) = {0.9*H, 0, 0, 0.5};
+Point(3) = {0.9*H, h, 0, 0.5};
+Point(4) = {0, h, 0, 0.5};
+
+Point(5) = {H, 0, 0, 0.5};
+Point(6) = {H, h, 0, 0.5};
+
+Point(7) = {H, H, 0, 0.5};
+Point(8) = {0.9*H, H, 0, 0.5};
+
+Line(1) = {1, 2};
+Line(2) = {2, 3};
+Line(3) = {3, 4};
+Line(4) = {4, 1};
+Curve Loop(1) = {1, 2, 3, 4};
+Plane Surface(1) = {1};
+
+Line(5) = {2, 5};
+Line(6) = {5, 6};
+Line(7) = {6, 3};
+Curve Loop(2) = {5, 6, 7, -2};
+Plane Surface(2) = {2};
+
+Line(8) = {6, 7};
+Line(9) = {7, 8};
+Line(10) = {8, 3};
+Curve Loop(3) = {-7, 8, 9, 10};
+Plane Surface(3) = {3};
+
+Transfinite Curve {1, 3, 8, 10} = 9*n+1 Using Progression 1;
+Transfinite Curve {2, 4, 5, 6, 7, 9} = n+1 Using Progression 1;
+Transfinite Surface {1};
+Transfinite Surface {2};
+Transfinite Surface {3};
+
+Recombine Surface {1};
+Recombine Surface {2};
+Recombine Surface {3};
+
+Physical Curve("left_edge", 1) = {4};
+Physical Surface("FEM_domain", 2) = {1, 2, 3}; // the trick is to include both plane surfaces in one single domain
+Physical Curve("top_edge", 3) = {9};
+
+F = 500;
+
+// additional parameters given to the solver
+SetNumber("Boundary Conditions/left_edge/ux", 0.); // ALWAYS NEED TO IMPOSE BOTH ux AND uy ON A GIVEN EDGE !! (pas très réaliste, faut y réfléchir)
+SetNumber("Boundary Conditions/left_edge/uy", 0.);
+SetNumber("Materials/FEM_domain/Young", 3e7);
+SetNumber("Materials/FEM_domain/Poisson", 0.3);
+SetNumber("Materials/FEM_domain/rho",7800); //volumic mass of acier
+SetNumber("Boundary Conditions/top_edge/tx", F); // ALWAYS NEED TO IMPOSE BOTH tx AND ty ON A GIVEN EDGE (realiste, OK) !
+SetNumber("Boundary Conditions/top_edge/ty", 0); //set to some other value for vertical deflection
+SetNumber("Volumic Forces/FEM_domain/bx",0.);
+SetNumber("Volumic Forces/FEM_domain/by",0.); //set to -9.81 for gravity
+
+Physical Curve("BEM_FEM_boundary", 4) = {4};
\ No newline at end of file
-- 
GitLab