modularize…

hertz_2D.geo
View Options

			cl1 = 0.01;
			cl2 = 0.1;
			cl3 = 1;
			Dy = 0.1;
			Point(1) = {0, 0.1-Dy, 0, cl1};
			Point(2) = {0.5, 0.6-Dy, 0, cl2};
			Point(3) = {-0.5, 0.6-Dy, 0, cl2};
			Point(4) = {0.5, 0, 0, cl3};
			Point(5) = {-0.5, 0, 0, cl3};
			Point(6) = {-0.5, -0.25, 0, cl3};
			Point(7) = {0.5, -0.25, 0, cl3};
			Point(8) = {0, 0.6-Dy, 0, cl2};
			Circle(1) = {3, 8, 1};
			Circle(2) = {1, 8, 2};
			Line(3) = {2, 8};
			Line(13) = {8, 3};
			Line(4) = {6, 7};
			Line(5) = {7, 4};
			Line(6) = {4, 5};
			Line(7) = {5, 6};
			Line Loop(9) = {2, 3, 13, 1};
			Plane Surface(9) = {9};
			Line Loop(11) = {6, 7, 4, 5};
			Plane Surface(11) = {11};
			Physical Line("rigid") = {6};
			Physical Line("contact_surface") = {1, 2};
			Physical Line("top") = {3, 13};
			Physical Surface("top_body") = {9};
			Physical Surface("bottom_body") = {11};