diff --git a/Coil_geometry/13_Properties_new_HH_coil_incr_N.py b/Coil_geometry/13_Properties_new_HH_coil_incr_N.py index 54efb6d..a487ab7 100644 --- a/Coil_geometry/13_Properties_new_HH_coil_incr_N.py +++ b/Coil_geometry/13_Properties_new_HH_coil_incr_N.py @@ -15,12 +15,43 @@ from IPython import get_ipython I = 1.25 -HH_Coil = BC.BCoil(HH = 1, distance = 54 ,radius = 48 , layers = 8, windings = 8, wire_height = 0.45, wire_width = 0.45,insulation_thickness= 0.09,is_round = True, winding_offset= True) -HH_Coil.set_R_outer(49.3) -HH_Coil.set_d_min(49.8) -HH_Coil.print_info() +Wire_1 = BC.BCoil(HH = 1, distance = 54 ,radius = 48 , layers = 8, windings = 8, wire_height = 0.425, wire_width = 0.425,insulation_thickness= 0.09,is_round = True, winding_offset= False) +Wire_1.set_R_outer(49.3) +Wire_1.set_d_min(49.8) +Wire_1.print_info() +print(Wire_1.get_coil_width() * 1e3 * Wire_1.get_coil_height() * 1e3) -HH_Coil.cooling(I,22.5) +Wire_1.cooling(I,22.5) -HH_Coil.plot_raster(30) +Wire_1.plot_raster(30) + +Wire_2 = BC.BCoil(HH = 1, distance = 54 ,radius = 48 , layers = 8, windings = 8, wire_height = 0.425, wire_width = 0.425,insulation_thickness= 0.09,is_round = False, winding_offset= False) +Wire_2.set_R_outer(49.3) +Wire_2.set_d_min(49.8) +Wire_2.print_info() +print(Wire_2.get_coil_width() * 1e3 * Wire_2.get_coil_height() * 1e3) + +Wire_2.cooling(I,22.5) + +Wire_2.plot_raster(30) + +Wire_2 = BC.BCoil(HH = 1, distance = 54 ,radius = 48 , layers = 8, windings = 8, wire_height = 0.55, wire_width = 0.55,insulation_thickness= 0.09,is_round = True, winding_offset= True) +Wire_2.set_R_outer(49.3) +Wire_2.set_d_min(49.8) +Wire_2.print_info() +print(Wire_2.get_coil_width() * 1e3 * Wire_2.get_coil_height() * 1e3) +Wire_2.cooling(I,22.5) + +Wire_2.plot_raster(30) + +I = 64/42 * 1.25 +Wire_2 = BC.BCoil(HH = 1, distance = 54 ,radius = 48 , layers = 6, windings = 7, wire_height = 0.55, wire_width = 0.55,insulation_thickness= 0.09,is_round = True, winding_offset= True) +Wire_2.set_R_outer(49.3) +Wire_2.set_d_min(49.8) +Wire_2.print_info() +print(Wire_2.get_coil_width() * 1e3 * Wire_2.get_coil_height() * 1e3) +Wire_2.cooling(I,22.5) +print(I) + +Wire_2.plot_raster(30) diff --git a/src/coil_class.py b/src/coil_class.py index b5dc789..0d1a2d5 100644 --- a/src/coil_class.py +++ b/src/coil_class.py @@ -199,12 +199,19 @@ class BCoil: return full_ras def plot_raster(self, raster_value): - full_structure = self.full_raster(raster_value) * 1e3 - - plt.figure(77) - plt.scatter(full_structure[:,:,1],full_structure[:,:,0]) + full_structure = self.full_raster(100) * 1e3 + if self.get_coil_width() > self.get_coil_height(): + extension = self.get_coil_width() + else: + extension = self.get_coil_height() + extension *= 1e3 + plt.figure(77,figsize=(5,5)) + plt.scatter(full_structure[:,:,1],full_structure[:,:,0], linewidths= 0.01) plt.xlabel("radius [mm]") plt.ylabel("z position [mm]") + plt.xlim(1e3*self.get_R_inner()-0.5, 1e3*self.get_R_inner() + extension + 0.5) + plt.ylim(1e3*self.get_zmin()- 0.5,1e3*self.get_zmin() + extension + 0.5) + plt.show() plt.close(77) @@ -583,12 +590,11 @@ class BCoil: return BCoil.resistivity_copper(T) * self.get_wire_length() / self.get_wire_area() -HH_Coil = BCoil(HH=1, distance=50, radius=40, layers=4, windings=4, wire_height=2, wire_width=1.5, - insulation_thickness=0.1, is_round=True, winding_offset=False) -HH_Coil.set_R_outer(70) -HH_Coil.print_info() +#HH_Coil = BCoil(HH=1, distance=50, radius=40, layers=4, windings=4, wire_height=0.4, wire_width=0.4,insulation_thickness=0.1, is_round=False, winding_offset=False) +#HH_Coil.set_R_outer(70) +#HH_Coil.print_info() -raster = HH_Coil.full_raster(10) +#raster = HH_Coil.plot_raster(10) #print(raster[:,:,0])