##PARAMETERS##
reset

#{air_X = 1}
#{magnet_Y = 0.2}
#{magnet_X = 0.14}

#{core_Y = 0.16}
#{core_X = 0.2}
#{coil_Y = 0.05}
#{coil_X = 0.05}

#{gap_size_Y = 0.1}
#{gap_size_X = 0.3}

#{mesh_size = 0.02}

##GEOMETRY##

create surface rectangle width {air_X} height {air_X} zplane 
create surface rectangle width {magnet_X} height {magnet_Y} zplane 
move surface 2 y {magnet_Y/2+gap_size_Y/2} include_merged
create surface rectangle width {core_X} height {core_Y} zplane
move surface 3 y {-core_Y/2-gap_size_Y/2} include_merged
create surface rectangle width {coil_X} height {core_Y} zplane
move surface 4 x {core_X/2+coil_X/2} y {-core_Y/2-gap_size_Y/2} include_merged
create surface rectangle width {coil_X} height {core_Y} zplane
move surface 5 x {-core_X/2-coil_X/2} y {-core_Y/2-gap_size_Y/2} include_merged
subtract surface 2 3 4 5  from surface 1  keep_tool 
webcut body 1 with plane normal to curve 8  fraction 0 from start 
webcut body 6  with plane normal to curve 16  fraction 0 from end 
section volume all with xplane offset 0 normal
delete body 5
sweep surface all perpendicular distance {mesh_size}  keep 
delete body 1 to 4, 6, 7
unite volume 12 9
imprint volume all 
merge volume all 

##REGIONS##

set duplicate block elements off
block 1 add volume 10 
block 1 name "V_Magnet"
block 2 add volume 11
block 2 name "V_Core"
block 3 add volume 8
block 3 name "V_Coil"
block 4 add volume 12
block 4 name "V_Air_up"
block 5 add volume 13 
block 5 name "V_Air_down"

block all element type HEX20

nodeset 1 add surface 43 51 45 33 27 53 39
nodeset 1 name "S_Flux_y"
nodeset 2 add surface 22 52
nodeset 2 name "S_Flux_x"
nodeset 3 add surface 22 32 43 52 51 25 26
nodeset 3 name "S_Smooth_fix"
nodeset 4 add surface 36 19 16 17 38
nodeset 4 name "S_Smooth_c"
nodeset 6 add surface 39  
nodeset 6 name "S_fix_xz"
nodeset 7 add surface 45 53 27
nodeset 7 name "S_smooth_fix_xz"
nodeset 8 add surface 40 54 20 46 34  
nodeset 8 name "S_Flux_front"
nodeset 9 add surface 35 47 21 41 55  
nodeset 9 name "S_Flux_back"


## MESH ##

volume all size {mesh_size}
volume all scheme auto
mesh volume all

nodeset 10 add node 51  
nodeset 10 name "P_Spring"

export ansys cdb "magnet_coil_mesh.cdb" geometry overwrite