SetFactory("OpenCASCADE");
h0 = 0.15e-3;
h1 = 2e-3;
D = 15e-3;
d = 13e-3;
di = 0.01e-3;
r = 0.5 * (D - d);
L = 2 * D;
t = 0.5 * (D - d);
w = 2 * Sqrt(r^2 - (r - t)^2);

Mesh.SecondOrderIncomplete=1;
Point(1) = {-0.5 * L, 0.0, 0.5 * di, h1};
Point(2) = {-0.5 * L, 0.0, 0.5 * D, h1};
Point(3) = {-0.5 * w, 0.0, 0.5 * D, h0};
Point(4) = {     0.0, 0.0, 0.5 * d, h0};
Point(5) = {     0.0, 0.0, 0.5 * d + r, h0};
Point(6) = { 0.5 * w, 0.0, 0.5 * D, h0};
Point(7) = { 0.5 * L, 0.0, 0.5 * D, h1};
Point(8) = { 0.5 * L, 0.0, 0.5 * di, h1};
Line(1) = {1,2};
Line(2) = {2,3};
Circle(3) = {3,5,4};
Circle(4) = {4,5,6};
Line(5) = {6,7};
Line(6) = {7,8};
Line(7) = {8,1};
Line Loop(9) = {1,2,3,4,5,6,7};
Plane Surface(11) = {9};
v[] = Extrude {{1.0,0.0,0.0},{0.0,0.0,0.0},-Pi/2} {
  Surface{11}; Layers{Ceil(d * Pi / 4 / h0)}; Recombine;
};
Recombine Surface{11,v[0]};
Physical Volume("volume") = {v[1]};
Physical Surface("bottom") = {v[0]};
Physical Surface("front") = {11};
Physical Surface("clamp") = {v[2]};
Physical Surface("load") = {v[7]};