39-the-python…

cohesive-contact.geo
View Options

			h = 0.5;

			y = 2;
			x = 1;

			Point(1) = { x/2., y/2, 0, h};
			Point(2) = {-x/2., y/2, 0, h};
			Point(3) = {-x/2.,-y/2, 0, h};
			Point(4) = { x/2.,-y/2, 0, h};
			Point(5) = {-x/2., 0, 0, h};
			Point(6) = { x/2., 0, 0, h};

			Line(1) = {1, 2};
			Line(2) = {2, 5};
			Line(3) = {5, 6};
			Line(4) = {6, 1};

			Line(5) = {5, 3};
			Line(6) = {3, 4};
			Line(7) = {4, 6};

			Line Loop(1) = {1, 2, 3, 4};
			Line Loop(2) = {-3, 5, 6, 7};
			Plane Surface(1) = {1};
			Plane Surface(2) = {2};

			Physical Line("fixed") = {6};
			Physical Line("loading") = {1};
			Physical Line("interface") = {3};
			Physical Line("sides") = {2, 5, 7, 4};

			Physical Surface("body") = {1, 2};

			Recombine Surface "*";
			Transfinite Surface "*";
			Mesh.SecondOrderIncomplete = 1;