reset 
create vertex 0 190.786 0 
create vertex -75.4151 131.1413 0 
create vertex 0 113.2860 0 
create curve arc center vertex 3 2 1  
delete vertex 3 
create vertex -85.7842 87.3459 0 
create curve vertex 2  vertex 5 
delete vertex 1
#{a = 22.1}
#{b = 35.1}
#{centx = -106.4700}
#{centy = 64.6 + b}
create vertex {centx} {centy-b} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (0.25*0.25) )}
create vertex {centx+0.25} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (0.5*0.5) )}
create vertex {centx+0.5} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (1*1) )}
create vertex {centx+1} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (2*2) )}
create vertex {centx+2} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (3*3) )}
create vertex {centx+3} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (4*4) )}
create vertex {centx+4} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (5*5) )}
create vertex {centx+5} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (7.5*7.5) )}
create vertex {centx+7.5} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (10*10) )}
create vertex {centx+10} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (12.5*12.5) )}
create vertex {centx+12.5} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (15*15) )}
create vertex {centx+15} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (16*16) )}
create vertex {centx+16} {-sqrt(temp)+centy} 0

create vertex on curve 2 fraction 0.2
create vertex on curve 2 fraction 0.4
create vertex on curve 2 fraction 0.6
create vertex on curve 2 fraction 0.8
create curve spline vertex 7 8 9 10 11 12 13 14 15 16 17 18 19 5 23 22 21 20 6 delete
vertex vis on
delete curve 2

curve all copy reflect X 
curve all copy move x -212.94
curve all copy move x 425.88
delete curve 12 13
curve 8 9 copy move -212.94

#{A = 74.45}
#{B = 83.27}
#{CENTX = -425.88}
#{CENTY = 107.515994}
#{temp = B*B - ( ((B*B)/(A*A)) * ((-427-CENTX)*(-427-CENTX)) )}
create vertex -427 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((-428-CENTX)*(-428-CENTX)) )}
create vertex -428 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((-430-CENTX)*(-430-CENTX)) )}
create vertex -430 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((-435-CENTX)*(-435-CENTX)) )}
create vertex -435 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((-440-CENTX)*(-440-CENTX)) )}
create vertex -440 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((-445-CENTX)*(-445-CENTX)) )}
create vertex -445 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((-450-CENTX)*(-450-CENTX)) )}
create vertex -450 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((-455-CENTX)*(-455-CENTX)) )}
create vertex -455 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((-460-CENTX)*(-460-CENTX)) )}
create vertex -460 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((-465-CENTX)*(-465-CENTX)) )}
create vertex -465 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((-470-CENTX)*(-470-CENTX)) )}
create vertex -470 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((-475-CENTX)*(-475-CENTX)) )}
create vertex -475 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((-480-CENTX)*(-480-CENTX)) )}
create vertex -480 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((-485-CENTX)*(-485-CENTX)) )}
create vertex -485 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((-490-CENTX)*(-490-CENTX)) )}
create vertex -490 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((-495-CENTX)*(-495-CENTX)) )}
create vertex -495 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((-497-CENTX)*(-497-CENTX)) )}
create vertex -497 {sqrt(temp)+CENTY} 0

#{a = 18.5}
#{b = 24.9}
#{centx = CENTX-103.07}
#{centy = 65+b}
create vertex {centx} {centy-b} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (0.25*0.25) )}
create vertex {centx+0.25} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (0.5*0.5) )}
create vertex {centx+0.5} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (1*1) )}
create vertex {centx+1} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (2*2) )}
create vertex {centx+2} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (3*3) )}
create vertex {centx+3} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (4*4) )}
create vertex {centx+4} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (5*5) )}
create vertex {centx+5} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (7.5*7.5) )}
create vertex {centx+7.5} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (10*10) )}
create vertex {centx+10} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (12.5*12.5) )}
create vertex {centx+12.5} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (15*15) )}
create vertex {centx+15} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (16*16) )}
create vertex {centx+16} {-sqrt(temp)+centy} 0

create curve spline vertex 74 75 76 77 78 79 80 81 82 83 84 85 86 73 72 71 70 69 68 67 66 65 64 63 62 61 60 59 58 57 54 delete
vertex vis on

#{A = 74.45}
#{B = 76.89}
#{CENTX = 425.88}
#{CENTY = 190.785994-B}
#{temp = B*B - ( ((B*B)/(A*A)) * ((427-CENTX)*(427-CENTX)) )}
create vertex 427 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((428-CENTX)*(428-CENTX)) )}
create vertex 428 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((430-CENTX)*(430-CENTX)) )}
create vertex 430 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((435-CENTX)*(435-CENTX)) )}
create vertex 435 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((440-CENTX)*(440-CENTX)) )}
create vertex 440 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((445-CENTX)*(445-CENTX)) )}
create vertex 445 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((450-CENTX)*(450-CENTX)) )}
create vertex 450 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((455-CENTX)*(455-CENTX)) )}
create vertex 455 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((460-CENTX)*(460-CENTX)) )}
create vertex 460 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((465-CENTX)*(465-CENTX)) )}
create vertex 465 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((470-CENTX)*(470-CENTX)) )}
create vertex 470 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((475-CENTX)*(475-CENTX)) )}
create vertex 475 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((480-CENTX)*(480-CENTX)) )}
create vertex 480 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((485-CENTX)*(485-CENTX)) )}
create vertex 485 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((490-CENTX)*(490-CENTX)) )}
create vertex 490 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((495-CENTX)*(495-CENTX)) )}
create vertex 495 {sqrt(temp)+CENTY} 0
#{temp = B*B - ( ((B*B)/(A*A)) * ((497-CENTX)*(497-CENTX)) )}
create vertex 497 {sqrt(temp)+CENTY} 0

#{a = 18.5}
#{b = 24.9}
#{centx = CENTX+103.07}
#{centy = 70+b}
create vertex {centx} {centy-b} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (0.25*0.25) )}
create vertex {centx-0.25} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (0.5*0.5) )}
create vertex {centx-0.5} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (1*1) )}
create vertex {centx-1} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (2*2) )}
create vertex {centx-2} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (3*3) )}
create vertex {centx-3} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (4*4) )}
create vertex {centx-4} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (5*5) )}
create vertex {centx-5} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (7.5*7.5) )}
create vertex {centx-7.5} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (10*10) )}
create vertex {centx-10} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (12.5*12.5) )}
create vertex {centx-12.5} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (15*15) )}
create vertex {centx-15} {-sqrt(temp)+centy} 0
#{temp = b*b - ( ((b*b)/(a*a)) * (16*16) )}
create vertex {centx-16} {-sqrt(temp)+centy} 0

create curve spline vertex 105 106 107 108 109 110 111 112 113 114 115 116 117 104 103 102 101 100 99 98 97 96 95 94 93 92 91 90 89 88 38 delete
vertex vis on

create vertex -613.95 65 0
create vertex -641.95 65 0
create vertex -671.95 65 0
create vertex -671.95 0 0

create vertex 672.95 70 0
create vertex 706.95 65 0
create vertex 736.95 65 0
create vertex 736.95 0 0 

create curve vertex 74 vertex 119
create curve vertex 119 vertex 120
create curve vertex 120 vertex 121
create curve vertex 121 vertex 122
create curve vertex 122 vertex 126
create curve vertex 126 vertex 125
create curve vertex 125 vertex 124
create curve vertex 124 vertex 123
create curve vertex 123 vertex 105

create surface curve all
sweep surface 1 xaxis angle 360

create cylinder height 300 radius 50
rotate volume 2 about X angle 90
volume 2 move x 614.58 y 130
unite volume 1 2
modify curve 84 blend radius 7

create cylinder height 300 radius 21.75
rotate volume 3 about x angle 90
volume 3 move x 614.58 y 210
modify curve 94 blend radius 5
subtract volume 3 from volume 1

webcut volume 1 with plane yplane offset 200 noimprint nomerge
delete volume 4

volume 1 move x 671.95
volume 1 copy move x -1408.9
volume 1 copy move x -2816.8
volume 1 copy move x 1408.9
unite volume all

#{factor = 1 - 0/212.94}
volume 1 scale 0.001
volume all scale X {factor} Y 1.0 Z 1.0

rotate volume 1 about y angle 270

volume 1 sizing function type skeleton scale 5 time_accuracy_level 2 min_size auto max_size 0.016 max_gradient 1.5
volume 1 scheme tetmesh
mesh volume 1

block 1 volume all
block 1 element type TETRA10

sideset 6 surface all

quality volume 1 shape global

export mesh "beta1.gen" overwrite