def _make_pf_coils(self): list_of_components = [] if (self.pf_coil_vertical_thicknesses, self.pf_coil_radial_thicknesses, self.pf_coil_to_rear_blanket_radial_gap) != (None, None, None): self._pf_coil = paramak.PoloidalFieldCoilSet( heights=self.pf_coil_vertical_thicknesses, widths=self.pf_coil_radial_thicknesses, center_points=self._pf_coils_xy_values, rotation_angle=self.rotation_angle, stp_filename='pf_coils.stp', stl_filename='pf_coils.stl', name="pf_coil", material_tag="pf_coil_mat", ) self._pf_coils_casing = paramak.PoloidalFieldCoilCaseSetFC( pf_coils=self._pf_coil, casing_thicknesses=self.pf_coil_case_thickness, rotation_angle=self.rotation_angle, stp_filename='pf_coil_cases.stp', stl_filename='pf_coil_cases.stl', name="pf_coil_case", material_tag="pf_coil_case_mat", ) list_of_components = [self._pf_coils_casing, self._pf_coil] return list_of_components
def test_PoloidalFieldCoilCaseSetFC_with_number_thickness(self): """Creates a set of PF coil cases from a list of PF coils with a single numerical thicknesses.""" pf_coils_1 = paramak.PoloidalFieldCoil(height=10, width=10, center_point=(100, 100)) pf_coils_2 = paramak.PoloidalFieldCoil(height=10, width=10, center_point=(100, 150)) pf_coils_3 = paramak.PoloidalFieldCoil(height=20, width=20, center_point=(50, 200)) pf_coils_4 = paramak.PoloidalFieldCoil(height=20, width=40, center_point=(50, 50)) test_shape = paramak.PoloidalFieldCoilCaseSetFC( pf_coils=[pf_coils_1, pf_coils_2, pf_coils_3, pf_coils_4], casing_thicknesses=10, ) assert test_shape.casing_thicknesses == 10 assert test_shape.solid is not None assert len(test_shape.solid.Solids()) == 4
def _make_coils(self): list_of_components = [] self._inboard_tf_coils = paramak.CenterColumnShieldCylinder( height=self._blanket_rear_wall_end_height * 2, inner_radius=self._inboard_tf_coils_start_radius, outer_radius=self._inboard_tf_coils_end_radius, rotation_angle=self.rotation_angle, stp_filename="inboard_tf_coils.stp", stl_filename="inboard_tf_coils.stl", name="inboard_tf_coils", material_tag="inboard_tf_coils_mat", ) list_of_components.append(self._inboard_tf_coils) if self._tf_info_provided: self._tf_coil = paramak.ToroidalFieldCoilCoatHanger( with_inner_leg=False, horizontal_start_point=( self._inboard_tf_coils_start_radius, self._blanket_rear_wall_end_height, ), vertical_mid_point=(self._outboard_tf_coil_start_radius, 0), thickness=self.outboard_tf_coil_radial_thickness, number_of_coils=self.number_of_tf_coils, distance=self.outboard_tf_coil_poloidal_thickness, stp_filename="outboard_tf_coil.stp", stl_filename="outboard_tf_coil.stl", rotation_angle=self.rotation_angle, horizontal_length=self._outboard_tf_coils_horizontal_length, vertical_length=self._outboard_tf_coils_vertical_height) list_of_components.append(self._tf_coil) if self._pf_info_provided: self._pf_coil = paramak.PoloidalFieldCoilSet( heights=self.pf_coil_vertical_thicknesses, widths=self.pf_coil_radial_thicknesses, center_points=self._pf_coils_xy_values, rotation_angle=self.rotation_angle, stp_filename='pf_coils.stp', stl_filename='pf_coils.stl', name="pf_coil", material_tag="pf_coil_mat", ) list_of_components.append(self._pf_coil) self._pf_coils_casing = paramak.PoloidalFieldCoilCaseSetFC( pf_coils=self._pf_coil, casing_thicknesses=self.pf_coil_case_thickness, rotation_angle=self.rotation_angle, stp_filename='pf_coil_cases.stp', stl_filename='pf_coil_cases.stl', name="pf_coil_case", material_tag="pf_coil_case_mat", ) list_of_components.append(self._pf_coils_casing) return list_of_components
def test_PoloidalFieldCoilCaseSetFC_from_list(self): """Creates a set of PF coil cases from a list of PF coils""" pf_coils_1 = paramak.PoloidalFieldCoil(height=10, width=10, center_point=(100, 100), rotation_angle=180) pf_coils_2 = paramak.PoloidalFieldCoil(height=10, width=10, center_point=(100, 150), rotation_angle=180) pf_coils_3 = paramak.PoloidalFieldCoil(height=20, width=20, center_point=(50, 200), rotation_angle=180) pf_coils_4 = paramak.PoloidalFieldCoil(height=20, width=40, center_point=(50, 50), rotation_angle=180) test_shape = paramak.PoloidalFieldCoilCaseSetFC( pf_coils=[pf_coils_1, pf_coils_2, pf_coils_3, pf_coils_4], casing_thicknesses=[5, 5, 10, 10], rotation_angle=180) assert test_shape.solid is not None assert len(test_shape.solid.Solids()) == 4
def test_PoloidalFieldCoilSet_incorrect_lengths(): """Checks PoloidalFieldCoilSet with the wrong number of casing thicknesses using a list""" pf_coils_1 = paramak.PoloidalFieldCoil(height=10, width=10, center_point=(100, 100), rotation_angle=180) paramak.PoloidalFieldCoilCaseSetFC( pf_coils=[pf_coils_1], casing_thicknesses=[5, 5, 10, 10], rotation_angle=180).solid
def test_PoloidalFieldCoilSet_incorrect_lengths_FC(): """Checks PoloidalFieldCoilSet with the wrong number of casing thicknesses using a coil set object""" pf_coils_set = paramak.PoloidalFieldCoilSet( heights=[10, 10, 20, 20], widths=[10, 10, 20, 40], center_points=[(100, 100), (100, 150), (50, 200), (50, 50)], rotation_angle=180) paramak.PoloidalFieldCoilCaseSetFC(pf_coils=pf_coils_set, casing_thicknesses=[5, 5, 10], rotation_angle=180).solid
def setUp(self): self.pf_coils_set = paramak.PoloidalFieldCoilSet( heights=[10, 10, 20, 20], widths=[10, 10, 20, 40], center_points=[(100, 100), (100, 150), (50, 200), (50, 50)], ) self.test_shape = paramak.PoloidalFieldCoilCaseSetFC( pf_coils=self.pf_coils_set, casing_thicknesses=[5, 10, 5, 10], )
def main(): number_of_toroidal_field_coils = 8 angle_offset = (360 / number_of_toroidal_field_coils) / 2. tf_coil_thickness = 50 tf_coil_distance = 130 tf_coil_casing_distance = tf_coil_distance + 40 tf_coil_casing_thickness = 20 inner_tf_case = paramak.InnerTfCoilsFlat( height=1800, inner_radius=330, outer_radius=430, number_of_coils=number_of_toroidal_field_coils, gap_size=5, rotation_angle=180, azimuth_start_angle=angle_offset) tf_coils = paramak.ToroidalFieldCoilPrincetonD( R1=400, R2=1500, # height thickness=tf_coil_thickness, distance=tf_coil_distance, number_of_coils=number_of_toroidal_field_coils, rotation_angle=180, ) tf_coil_casing = paramak.TFCoilCasing( magnet=tf_coils, distance=tf_coil_casing_distance, inner_offset=tf_coil_casing_thickness, outer_offset=tf_coil_casing_thickness, vertical_section_offset=tf_coil_casing_thickness, # rotation_angle=180, # producing occational errors with this arg ) pf_coils = paramak.PoloidalFieldCoilSet( heights=[100, 120, 80, 80, 120, 180], widths=[100, 120, 80, 80, 120, 180], center_points=[(530, 1030), (1370, 790), (1740, 250), (1750, -250), (1360, -780), (680, -1050)], rotation_angle=180) pf_coils_casing = paramak.PoloidalFieldCoilCaseSetFC( pf_coils=pf_coils, casing_thicknesses=[10] * 6, rotation_angle=180) pf_coils.export_stp('pf_coils.stp') pf_coils_casing.export_stp('pf_coils_case.stp') tf_coils.export_stp('tf_coil.stp') tf_coil_casing.export_stp('tf_coil_casing.stp') inner_tf_case.export_stp('inner_tf_case.stp')
def test_PoloidalFieldCoilCaseSetFC_from_pf_coil_set(self): """Creates a set of PF coil cases from a PF coils object""" pf_coils_set = paramak.PoloidalFieldCoilSet(heights=[10, 10, 20, 20], widths=[10, 10, 20, 40], center_points=[(100, 100), (100, 150), (50, 200), (50, 50)], rotation_angle=180) test_shape = paramak.PoloidalFieldCoilCaseSetFC( pf_coils=pf_coils_set, casing_thicknesses=[5, 5, 10, 10], rotation_angle=180) assert test_shape.solid is not None assert len(test_shape.solid.Solids()) == 4 assert len(pf_coils_set.solid.Solids()) == 4
def main(): rot_angle = 180 all_components = [] plasma = paramak.Plasma( # default parameters rotation_angle=rot_angle, stp_filename="plasma_shape.stp", ) all_components.append(plasma) component = paramak.BlanketFP( plasma=plasma, thickness=100, stop_angle=90, start_angle=-90, offset_from_plasma=30, rotation_angle=rot_angle, stp_filename="blanket_constant_thickness_outboard_plasma.stp", ) all_components.append(component) component = paramak.BlanketCutterStar(height=2000, width=2000, distance=100) all_components.append(component) component = paramak.BlanketFP( plasma=plasma, thickness=100, stop_angle=90, start_angle=250, offset_from_plasma=30, rotation_angle=rot_angle, stp_filename="blanket_constant_thickness_inboard_plasma.stp", ) all_components.append(component) component = paramak.BlanketFP( plasma=plasma, thickness=100, stop_angle=250, start_angle=-90, offset_from_plasma=30, rotation_angle=rot_angle, stp_filename="blanket_constant_thickness_plasma.stp", ) all_components.append(component) CenterColumnShieldCylinder = paramak.CenterColumnShieldCylinder( inner_radius=80, outer_radius=100, height=300, rotation_angle=rot_angle, stp_filename="center_column_shield_cylinder.stp", ) all_components.append(CenterColumnShieldCylinder) component = paramak.InboardFirstwallFCCS( central_column_shield=CenterColumnShieldCylinder, thickness=50, rotation_angle=rot_angle, stp_filename="firstwall_from_center_column_shield_cylinder.stp", ) all_components.append(component) CenterColumnShieldHyperbola = paramak.CenterColumnShieldHyperbola( inner_radius=50, mid_radius=75, outer_radius=100, height=300, rotation_angle=rot_angle, stp_filename="center_column_shield_hyperbola.stp", ) all_components.append(CenterColumnShieldHyperbola) component = paramak.InboardFirstwallFCCS( central_column_shield=CenterColumnShieldHyperbola, thickness=50, rotation_angle=rot_angle, stp_filename="firstwall_from_center_column_shield_hyperbola.stp", ) all_components.append(component) CenterColumnShieldCircular = paramak.CenterColumnShieldCircular( inner_radius=50, mid_radius=75, outer_radius=100, height=300, rotation_angle=rot_angle, stp_filename="center_column_shield_circular.stp", ) all_components.append(CenterColumnShieldCircular) component = paramak.InboardFirstwallFCCS( central_column_shield=CenterColumnShieldCircular, thickness=50, rotation_angle=rot_angle, stp_filename="firstwall_from_center_column_shield_circular.stp", ) all_components.append(component) CenterColumnShieldFlatTopHyperbola = paramak.CenterColumnShieldFlatTopHyperbola( inner_radius=50, mid_radius=75, outer_radius=100, arc_height=220, height=300, rotation_angle=rot_angle, stp_filename="center_column_shield_flat_top_hyperbola.stp", ) all_components.append(CenterColumnShieldFlatTopHyperbola) component = paramak.InboardFirstwallFCCS( central_column_shield=CenterColumnShieldFlatTopHyperbola, thickness=50, rotation_angle=rot_angle, stp_filename= "firstwall_from_center_column_shield_flat_top_hyperbola.stp", ) all_components.append(component) CenterColumnShieldFlatTopCircular = paramak.CenterColumnShieldFlatTopCircular( inner_radius=50, mid_radius=75, outer_radius=100, arc_height=220, height=300, rotation_angle=rot_angle, stp_filename="center_column_shield_flat_top_Circular.stp", ) all_components.append(CenterColumnShieldFlatTopCircular) component = paramak.InboardFirstwallFCCS( central_column_shield=CenterColumnShieldFlatTopCircular, thickness=50, rotation_angle=rot_angle, stp_filename= "firstwall_from_center_column_shield_flat_top_Circular.stp", ) all_components.append(component) CenterColumnShieldPlasmaHyperbola = paramak.CenterColumnShieldPlasmaHyperbola( inner_radius=150, mid_offset=50, edge_offset=40, height=800, rotation_angle=rot_angle, stp_filename="center_column_shield_plasma_hyperbola.stp", ) all_components.append(CenterColumnShieldPlasmaHyperbola) component = paramak.InboardFirstwallFCCS( central_column_shield=CenterColumnShieldPlasmaHyperbola, thickness=50, rotation_angle=rot_angle, stp_filename="firstwall_from_center_column_shield_plasma_hyperbola.stp", ) all_components.append(component) component = paramak.InnerTfCoilsCircular( inner_radius=25, outer_radius=100, number_of_coils=10, gap_size=5, height=300, stp_filename="inner_tf_coils_circular.stp", ) all_components.append(component) component = paramak.InnerTfCoilsFlat( inner_radius=25, outer_radius=100, number_of_coils=10, gap_size=5, height=300, stp_filename="inner_tf_coils_flat.stp", ) all_components.append(component) # this makes 4 pf coil cases pf_coil_set = paramak.PoloidalFieldCoilCaseSet( heights=[10, 10, 20, 20], widths=[10, 10, 20, 40], casing_thicknesses=[5, 5, 10, 10], center_points=[(100, 100), (100, 150), (50, 200), (50, 50)], rotation_angle=rot_angle, stp_filename="pf_coil_case_set.stp") all_components.append(pf_coil_set) # this makes 4 pf coils pf_coil_set = paramak.PoloidalFieldCoilSet(heights=[10, 10, 20, 20], widths=[10, 10, 20, 40], center_points=[(100, 100), (100, 150), (50, 200), (50, 50)], rotation_angle=rot_angle, stp_filename="pf_coil_set.stp") all_components.append(pf_coil_set) # this makes 4 pf coil cases for the 4 pf coils made above component = paramak.PoloidalFieldCoilCaseSetFC( pf_coils=pf_coil_set, casing_thicknesses=[5, 5, 10, 10], rotation_angle=rot_angle, stp_filename="pf_coil_cases_set.stp") all_components.append(component) # this makes 1 pf coils pf_coil = paramak.PoloidalFieldCoil(center_point=(100, 100), height=20, width=20, rotation_angle=rot_angle, stp_filename="poloidal_field_coil.stp") all_components.append(pf_coil) # this makes one PF coil case for the provided pf coil component = paramak.PoloidalFieldCoilCaseSetFC( pf_coils=[pf_coil], casing_thicknesses=[10], rotation_angle=rot_angle, stp_filename="pf_coil_cases_set_fc.stp") all_components.append(component) component = paramak.PoloidalFieldCoilCaseFC( pf_coil=pf_coil, casing_thickness=10, rotation_angle=rot_angle, stp_filename="poloidal_field_coil_case_fc.stp", ) all_components.append(component) component = paramak.PoloidalFieldCoilCase( center_point=(100, 100), coil_height=20, coil_width=20, casing_thickness=10, rotation_angle=rot_angle, stp_filename="poloidal_field_coil_case.stp", ) all_components.append(component) component = paramak.BlanketConstantThicknessArcV( inner_lower_point=(300, -200), inner_mid_point=(500, 0), inner_upper_point=(300, 200), thickness=100, rotation_angle=rot_angle, stp_filename="blanket_arc_v.stp", ) all_components.append(component) component = paramak.BlanketConstantThicknessArcH( inner_lower_point=(300, -200), inner_mid_point=(400, 0), inner_upper_point=(300, 200), thickness=100, rotation_angle=rot_angle, stp_filename="blanket_arc_h.stp", ) all_components.append(component) component = paramak.ToroidalFieldCoilRectangle( horizontal_start_point=(100, 700), vertical_mid_point=(800, 0), thickness=150, distance=60, stp_filename="tf_coil_rectangle.stp", number_of_coils=1, ) all_components.append(component) component = paramak.ToroidalFieldCoilCoatHanger( horizontal_start_point=(200, 500), horizontal_length=400, vertical_mid_point=(700, 50), vertical_length=500, thickness=50, distance=50, stp_filename="toroidal_field_coil_coat_hanger.stp", number_of_coils=1, ) all_components.append(component) component = paramak.ToroidalFieldCoilTripleArc( R1=80, h=200, radii=(70, 100), coverages=(60, 60), thickness=30, distance=30, number_of_coils=1, stp_filename="toroidal_field_coil_triple_arc.stp") all_components.append(component) magnet = paramak.ToroidalFieldCoilPrincetonD( R1=80, R2=300, thickness=30, distance=30, number_of_coils=1, stp_filename="toroidal_field_coil_princeton_d.stp") all_components.append(magnet) component = paramak.ITERtypeDivertor( # default parameters rotation_angle=rot_angle, stp_filename="ITER_type_divertor.stp", ) all_components.append(component) component = paramak.PortCutterRotated(center_point=(450, 0), polar_coverage_angle=20, rotation_angle=10, polar_placement_angle=45, azimuth_placement_angle=0) all_components.append(component) component = paramak.PortCutterRectangular( distance=3, z_pos=0, height=0.2, width=0.4, fillet_radius=0.02, azimuth_placement_angle=[0, 45, 90, 180]) all_components.append(component) component = paramak.PortCutterCircular( distance=3, z_pos=0.25, radius=0.1, # azimuth_placement_angle=[0, 45, 90, 180], # TODO: fix issue #548 azimuth_placement_angle=[0, 45, 90], ) all_components.append(component) component = paramak.VacuumVessel(height=2, inner_radius=1, thickness=0.2, rotation_angle=270) all_components.append(component) component = paramak.CoolantChannelRingStraight( height=200, channel_radius=10, ring_radius=70, number_of_coolant_channels=8, workplane="XY", rotation_axis="Z", stp_filename="coolant_channel_ring_straight.stp", ) all_components.append(component) component = paramak.CoolantChannelRingCurved( height=200, channel_radius=10, ring_radius=70, mid_offset=-20, number_of_coolant_channels=8, workplane="XY", path_workplane="XZ", stp_filename="coolant_channel_ring_curved.stp", force_cross_section=True) all_components.append(component) component = paramak.RotatedIsoscelesTriangle( height=20, base_length=15, pivot_angle=0, pivot_point=(100, 50), rotation_angle=70, workplane='XY', stp_filename='rotated_isosceles_triangle.stp') all_components.append(component) component = paramak.RotatedTrapezoid(length_1=10, length_2=20, length_3=30, pivot_angle=0, pivot_point=(100, 50), rotation_angle=45, stp_filename='rotated_trapezoid.stp') all_components.append(component) component = paramak.PoloidalSegments(number_of_segments=5, center_point=(400, 50)) all_components.append(component) component = paramak.TFCoilCasing(magnet=magnet, inner_offset=10, outer_offset=15, vertical_section_offset=20, distance=40) all_components.append(component) component = paramak.HexagonPin(length_of_side=5, distance=10, center_point=(10, 20)) all_components.append(component) return all_components
def test_PoloidalFieldCoilSet_incorrect_pf_coil(): """Checks PoloidalFieldCoilSet with the pf_coils as an incorrect entry""" paramak.PoloidalFieldCoilCaseSetFC( pf_coils=20, casing_thicknesses=[5, 5, 10, 10], rotation_angle=180).solid