femagtools 1.7.3__py3-none-any.whl → 1.7.5__py3-none-any.whl

femagtools/__init__.py

@@ -2,7 +2,7 @@
 
 """
 __title__ = 'femagtools'
-__version__ = '1.7.3'
+__version__ = '1.7.5'
 __author__ = 'Ronald Tanner'
 __license__ = 'BSD'
 __copyright__ = 'Copyright 2023-2024 Gamma Technology'

femagtools/airgap.py

@@ -52,3 +52,14 @@ def read(filename, pmod=0):
         return(dict())
 
     return fft(bag[0], bag[1], pmod)
+
+if __name__ == '__main__':
+    import sys
+    import matplotlib.pyplot as plt
+    import femagtools.plot.fluxdens
+    ag = read(sys.argv[1])
+    fig, axs = plt.subplots(nrows=2)
+    femagtools.plot.fluxdens.airgap(ag, ax=axs[0])
+    femagtools.plot.fluxdens.airgap_fft(ag, ax=axs[1])
+    fig.tight_layout()
+    plt.show()

femagtools/bch.py

@@ -70,7 +70,7 @@ def sttylosses(losses):
         except KeyError:
             pass
         return d
-    
+
     if sregs == {'StYoke', 'StTeeth'}:
         return hysteddy('StYoke', 'StTeeth', losses)
     if sregs == {'StJo', 'StZa'}:
@@ -95,9 +95,9 @@ def sttylosses(losses):
         return l
     return {}
 
-def losses_mapping_external_rotor(losses): 
+def losses_mapping_external_rotor(losses):
     styoke = 'rotor'
-    rotor = 'Iron'    
+    rotor = 'Iron'
     d = {}
     try:
         d['styoke'] = losses[styoke]
@@ -613,6 +613,7 @@ class Reader:
 
     def __read_general_machine_data(self, content):
         mcfiles = []
+        slotsides = 1
         for l in content:
             try:
                 if l.find('Armature Length') > -1:
@@ -631,8 +632,10 @@ class Reader:
                     self.machine['p_sim'] = int(l.split()[-1])
                 elif l.find('Total Number of Slots') > -1:
                     self.machine['Q'] = int(l.split()[-1])
+                elif l.find('Total Number of Slot-Sides') > -1:
+                    slotsides = int(l.split()[-1])
                 elif l.find('Number of Slot-Sides sim.') > -1:
-                    self.machine['qs_sim'] = int(l.split()[-1])
+                    self.machine['qs_sim'] = int(l.split()[-1])*self.machine['Q']//slotsides
                 elif l.find('POC-File used in calculation') > -1:
                     self.machine['pocfile'] = l.split(
                         ':')[-1].strip().replace('\\', '\\\\')
@@ -1046,7 +1049,7 @@ class Reader:
         '''
         diff = np.floor(np.abs(np.diff(idList)))
         if idList[-1] == 0 and len(idList) > 2 and \
-           not np.all(diff == diff[0]):
+           diff[-1] == 0:
             idList = idList[:-1]
         return idList
 
@@ -1062,7 +1065,8 @@ class Reader:
                 m.append([floatnan(x) for x in rec])
 
         m = np.array(m).T
-        ncols = np.argmax(np.abs(m[1][1:]-m[1][:-1]))+1
+        d = np.diff(m[1])
+        ncols = (len(d)+1)//(len(d[d < 0])+1)
         if ncols == 1 and len(m[1]) > 1 and m[1][0] != m[1][1]:  # simple correction
             ncols = 2
         iq = np.reshape(m[1], (-1, ncols))[0]
@@ -1071,7 +1075,8 @@ class Reader:
             ncols = ncols-1
 
         id = np.reshape(m[0], (-1, ncols)).T[0]
-        id = self.__removeTrailingZero(id)
+        if id[0] >= 0:
+            id = self.__removeTrailingZero(id)
         nrows = len(id)
         if nrows > 1 and id[nrows-1] < id[nrows-2]:
             nrows = nrows-1
@@ -1096,7 +1101,8 @@ class Reader:
                 m.append([floatnan(x) for x in rec])
 
         m = np.array(m).T
-        ncols = np.argmax(np.abs(m[1][1:]-m[1][:-1]))+1
+        d = np.diff(m[1])
+        ncols = (len(d)+1)//(len(d[d < 0])+1)
         if ncols == 1 and len(m[1]) > 1 and m[1][0] != m[1][1]:  # simple correction
             ncols = 2
         iq = np.linspace(np.min(m[1]), np.max(m[1]), ncols)
@@ -1105,7 +1111,8 @@ class Reader:
             ncols = ncols-1
 
         id = np.reshape(m[0], (-1, ncols)).T[0]
-        id = self.__removeTrailingZero(id)
+        if id[0] >= 0:
+            id = self.__removeTrailingZero(id)
         nrows = len(id)
         if nrows > 1 and id[nrows-1] < id[nrows-2]:
             nrows = nrows-1
@@ -1168,7 +1175,7 @@ class Reader:
                    for name in content[2].split('\t') if name][2:]
         # outer rotor motor
         if 'Iron' in subregs and \
-            'Rotor' in subregs: 
+            'Rotor' in subregs:
             self.external_rotor = True
         logger.info("Stator Subregions: %s", subregs)
         speed = float(content[0].split()[-1])/60.
@@ -1184,7 +1191,11 @@ class Reader:
         if not m:
             return
         m = np.array(m).T
-        ncols = np.argmax(np.abs(m[1][1:]-m[1][:-1]))+1
+        d = np.diff(m[1])
+        if self.ldq:
+            ncols = (len(d)+1)//(len(d[d > 0])+1)
+        else:
+            ncols = (len(d)+1)//(len(d[d < 0])+1)
         if ncols == 1 and len(m[1]) > 1 and m[1][0] != m[1][1]:
             ncols = 2
         id = np.reshape(m[0], (-1, ncols)).T[0]
@@ -1229,9 +1240,9 @@ class Reader:
                                          (nrows, ncols)).T.tolist()
                            for k, v in zip(cols, m[2:])})
         ls['speed'] = speed
-        if self.external_rotor: 
+        if self.external_rotor:
             ls.update(losses_mapping_external_rotor(ls))
-        else: 
+        else:
             ls.update(sttylosses(ls))
         if self.ldq:
             self.ldq['losses'] = ls
@@ -1487,9 +1498,9 @@ class Reader:
     def __read_losses(self, content):
         losses = {}
         i = 0
-        # check if external rotor 
+        # check if external rotor
         if self.weights[0][1] == 0.0 and \
-            self.weights[0][-1] > 0: 
+            self.weights[0][-1] > 0:
             self.external_rotor = True
         # find results for angle:
         while True:
@@ -1548,8 +1559,8 @@ class Reader:
                     losses['total'] += losses['staza']+losses['stajo']
                 elif len(rec) == 1:
                     t = l.split(':')[-1].strip()
-                    if t == 'Iron': 
-                        if self.external_rotor: 
+                    if t == 'Iron':
+                        if self.external_rotor:
                             losses['rotfe'] = floatnan(rec[0])
                             losses['total'] += losses['rotfe']
                         else:
@@ -1561,14 +1572,14 @@ class Reader:
                 continue
 
             if _rotloss.search(l):
-                if l.find('StZa') > -1: 
+                if l.find('StZa') > -1:
                     self.external_rotor = True
                 rec = self.__findNums(content[i+2])
                 if len(rec) == 1:
-                    if self.external_rotor: 
+                    if self.external_rotor:
                         losses['staza'] = floatnan(rec[0])
                         losses['total'] += losses['staza']
-                    else: 
+                    else:
                         rotfe = floatnan(rec[0])
                         losses['rotfe'] += rotfe
                         losses['total'] += rotfe
@@ -1602,7 +1613,7 @@ class Reader:
                     if content[i+1].split() == ['rotf', '----']:
                         losses['rotfe'] = sum([floatnan(x) for x in rec])
                         losses['total'] += losses['rotfe']
-                    
+
                     if content[i+1].split() == ['Iron', '----']: # external rotor
                         losses['rotfe'] = sum([floatnan(x) for x in rec])
                         losses['total'] += losses['rotfe']

femagtools/dxfsl/area.py

@@ -84,6 +84,7 @@ class Area(object):
         self.is_child = False
         self.areas_inside = {}
         self.areas_of_group = []
+        self.group_is_inside = False
 
     def identifier(self):
         return "{}-{}".format(self.id, self.type)

femagtools/dxfsl/fslrenderer.py

@@ -189,7 +189,7 @@ class FslRenderer(object):
                          '-- min_corner = {}, {}'.format(
                              geom.start_min_corner(0),
                              geom.start_min_corner(1)),
-                         '-- max_corner = {}, {}'.format(
+                         '-- max_corner start = {}, {}'.format(
                              geom.start_max_corner(0),
                              geom.start_max_corner(1)),
                          '\n']
@@ -211,11 +211,12 @@ class FslRenderer(object):
                 'r, phi = c2pr(x0, y0)',
                 'x1, y1 = pr2c(r1, phi)',
                 'x2, y2 = pr2c(r1, {}*math.pi/parts)'.format(slice),
+                f'r = {geom.dist_end_max_corner()}',
                 'x3, y3 = pr2c(r, {}*math.pi/parts)'.format(slice),
                 'nc_line(x0, y0, x1, y1, 0)',
                 'nc_circle_m(x1, y1, x2, y2, 0.0, 0.0, 0)',
                 'nc_line(x2, y2, x3, y3, 0)',
-                'x0, y0 = pr2c(r1 - hair/2, math.pi/parts)',
+                'x0, y0 = pr2c(r1 - hair/2, math.pi/parts/2)',
                 'create_mesh_se(x0, y0)',
                 '\n',
                 'outer_da_start = {}'.format(

femagtools/dxfsl/geom.py

@@ -2478,7 +2478,8 @@ class Geometry(object):
             groups_inside = {g.id: g for g in grouplist}
             area.areas_inside = groups_inside
             areas_outside.append(area)
-            for g in  grouplist:
+            for g in grouplist:
+                g.group_is_inside = True
                 alist = groups.get(g.id, [])
                 alist.append(area)
                 groups[g.id] = alist
@@ -2532,6 +2533,11 @@ class Geometry(object):
             for area in area_list:
                 if self.create_aux_lines(area, rightangle, leftangle):
                     done = True
+            main_groups = [g for g in self.areagroup_list
+                           if not g.group_is_inside]
+            if len(main_groups) > 1:
+                if self.create_outside_aux_lines(main_groups):
+                    done = True
         else:
             logger.debug("-> start create_auxiliary_lines")
             self.set_areas_inside_for_all_areas()
@@ -2856,6 +2862,38 @@ class Geometry(object):
                      area.get_id())
         return done
 
+    def create_outside_aux_lines(self, grouplist):
+        done = False
+        aux_color = 'blue'
+        aux_linestyle = 'dotted'
+
+        i = 0
+        gaps = []
+        for group in grouplist:
+            i += 1
+            for g in grouplist[i:]:
+                gaps += group.get_lowest_gap_list(g,
+                                                  self.center, self.max_radius,
+                                                  None, None)
+        gaps.sort()
+
+        l = len(grouplist) -1
+        for d, points, token, id in gaps[:l]:
+            logger.info("Token %s", token)
+            line = Line(Element(start=points[0],
+                                end=points[1]),
+                                color=aux_color,
+                                linestyle=aux_linestyle,
+                                attr='outside_auxline')
+            n1 = self.find_the_node(line.node1(ndec))
+            n2 = self.find_the_node(line.node2(ndec))
+            if n1 and n2:
+                self.add_element(line,
+                                 rtol=self.rtol,
+                                 atol=self.atol)
+                done = True
+        return done
+
     def set_rotor(self):
         self.sym_counterpart = 1
         self.sym_part = 0

femagtools/fsl.py

@@ -106,6 +106,7 @@ class Builder:
         params['show_plots'] = model.stator[templ].get('plot', False)
         params['write_fsl'] = True
         params['airgap'] = -1.0
+        params['nodedist'] = model.stator.get('nodedist', 1)
         pos = 'in' if model.external_rotor else 'out'
         params['part'] = ('stator', pos)
         conv = convert(model.stator['dxffile']['name'], **params)
@@ -491,6 +492,11 @@ class Builder:
         if not hasattr(model, 'stator'):
             setattr(model, 'stator', {})
         model.stator['num_slots'] = conv.get('tot_num_slot')
+        if model.get('num_agnodes', 0) == 0:
+            model.set_value('agndst', conv['agndst']*1e-3)
+            logger.info("num poles %d num slots %d outer diameter %.4f m agndst %.4f mm",
+                        model.poles, model.stator['num_slots'],
+                        model.outer_diam, model.agndst*1e3)
         if params['full_model']:
             model.stator['num_slots_gen'] = model.stator['num_slots']
         else:
@@ -531,21 +537,6 @@ class Builder:
             logger.info("create new model '%s'", model.name)
             if model.is_dxffile():
                 self.prepare_model_with_dxf(model)
-                if model.get('num_agnodes', 0) == 0:
-                    from femagtools.dxfsl.fslrenderer import agndst
-                    ag = model.get('airgap')
-                    model.set_value(
-                        'agndst',
-                        agndst(model.get('bore_diam'),
-                               model.get('bore_diam') - 2*ag,
-                               model.stator.get('num_slots'),
-                               model.get('poles'),
-                               model.dxffile.get('nodedist')))
-
-                logger.info(" num poles %d num slots %d outer diameter %.4f m",
-                            model.poles, model.stator['num_slots'],
-                            model.outer_diam)
-
             else:
                 self.prepare_stator(model)
                 if hasattr(model, 'magnet'):
@@ -553,9 +544,9 @@ class Builder:
                 self.prepare_diameter(model)
                 if self.fsl_stator and model.get('num_agnodes', 0) == 0:
                     from femagtools.dxfsl.fslrenderer import agndst
-                    if model.get('agndst'): 
+                    if model.get('agndst',0):
                         pass
-                    else: 
+                    else:
                         ag = model.get('airgap')
                         model.set_value(
                             'agndst',