32-python…

flat_2d.geo
View Options

			a = 1.0;
			r = a;
			gap = 0;
			nbpoints = 1.0;
			cl1 = a/(2.0*nbpoints);
			cl2 = r/(nbpoints);

			Point(1) = {-a/2, 0 + gap, 0, cl1};
			Point(2) = {a/2, 0 + gap, 0, cl1};
			Point(3) = {a/2, a/2., 0, cl1};
			Point(4) = {-a/2, a/2., 0, cl1};
			Point(5) = {r/2, 0, 0, cl2};
			Point(6) = {-r/2, 0, 0, cl2};
			Point(7) = {-r/2, -r/2, 0, cl2};
			Point(8) = { r/2, -r/2, 0, cl2};
			Line(2) = {1, 2};
			Line(3) = {2, 3};
			Line(4) = {3, 4};
			Line(5) = {4, 1};
			Line(6) = {7, 8};
			Line(7) = {8, 5};
			Line(8) = {5, 6};
			Line(9) = {6, 7};
			Line Loop(1) = {2, 3, 4, 5};
			Line Loop(2) = {6, 7, 8, 9};
			Plane Surface(1) = {1};
			Plane Surface(2) = {2};
			Physical Line("contact_top") = {2};
			Physical Line("contact_bottom") = {8};
			Physical Line("top") = {4};
			Physical Line("bottom") = {6};
			Physical Surface("bot_body") = {2};
			Physical Surface("top_body") = {1};