L = 0.1; H = 0.6; y = 0; lc = 40*1e-3; lc_min = 20*1e-3; // Defining the rotation alpha = 30*Pi/180; a1 = -L; b1 = -2*H; a2 = L; b2 = -2*H; a3 = L; b3 = 0; a4 = -L; b4 = 0; a5 = -L; b5 = -H; a6 = L; b6 = -H; s = Sin(alpha); c = Cos(alpha); mx = 0; my = -0.52; Point(1) = {c*a4-s*b4, s*a4+c*b4, 0,lc_min}; Point(2) = {c*a1-s*b1, s*a1+c*b1, 0,lc_min}; // Point(3) = {c*a2-s*b2, s*a2+c*b2, 0,lc_min}; // Point(4) = {c*a3-s*b3, s*a3+c*b3, 0,lc_min}; Point(5) = {c*a5-s*b5, s*a5+c*b5, 0,lc_min}; Point(6) = {c*a6-s*b6, s*a6+c*b6, 0,lc_min}; Line(1) = {1, 5}; // Lateral left Line(2) = {5, 2}; // Lateral left 2 Line(3) = {2, 3}; // Bottom Line(4) = {3, 6}; // Lateral right Line(5) = {6, 4}; // Lateral right 2 Line(6) = {4, 1}; // Top Point(7) = {c*mx-s*my,s*mx+c*my,0}; Line Loop(1) = {1:6}; Plane Surface(1) = {1}; Physical Line("Bottom") = {3}; Physical Line("Lateral") = {1,2, 4, 5}; Physical Line("Top") = {6}; Physical Surface("Domain") = {1}; Physical Point("PtFix") = {1}; Mesh.Algorithm = 5; Mesh.MshFileVersion = 2;