femagtools 1.7.6__py3-none-any.whl → 1.7.7__py3-none-any.whl

femagtools/fsl.py

@@ -42,7 +42,7 @@ class Builder:
             default_filters=['decode.utf8'])
 
         self.fsl_stator = False
-        self.fsl_magnet = False
+        self.fsl_rotor = False
 
     def create_wdg_def(self, model):
         name = 'winding'
@@ -216,7 +216,7 @@ class Builder:
                .format(model.magnet.get('mcvkey_shaft', 'dummy'))]
 
         if 'magnetFsl' in model.magnet:
-            self.fsl_magnet = True
+            self.fsl_rotor = True
             #  obsolete
             if 'parameter' in model.magnet['magnetFsl']:
                 return mcv + self.render_template(
@@ -242,10 +242,10 @@ class Builder:
         magmodel['mcvkey_magnet'] = model.get_mcvkey_magnet()
         if templ == 'dxf':
             return mcv + [
-                u'xmag = {}',
-                u'ymag = {}',
-                u'mag_orient = {}',
-                u'ndt(agndst)'] + model.magnet['dxf']['fsl']
+                'xmag = {}',
+                'ymag = {}',
+                'mag_orient = {}',
+                'ndt(agndst)'] + model.magnet['dxf']['fsl']
 
         return mcv + self.render_rotor(magmodel, templ)
 
@@ -267,12 +267,16 @@ class Builder:
                 rotmodel['winding_inside'] = not model.external_rotor
                 culosses = self.create_cu_losses(rotmodel, condMat)
 
-        rotmodel.update(model.rotor[templ])
         rotmodel['is_rotor'] = True  # just in case for the template
+        if templ == 'EESM':
+            if 'dxf' in rotmodel:
+                return mcv + ['ndt(agndst)'] + rotmodel['dxf']['fsl']
+            templ = 'rot_hsm'
+        rotmodel.update(model.rotor[templ])
         return mcv + culosses + self.render_rotor(rotmodel, templ)
 
     def create_rotor_winding(self, model):
-        if hasattr(model, 'rotor') and model.rotortype() == 'rot_hsm':
+        if hasattr(model, 'rotor') and model.rotortype() == 'EESM':
             return self.render_rotor(model.rotor, 'rotor_winding')
         return []
 
@@ -306,21 +310,22 @@ class Builder:
 
         if templ == 'dxf':
             # reuse dxfsl model
-            self.fsl_magnet = True
+            self.fsl_rotor = True
         if templ != 'dxffile':
             return
 
         from femagtools.dxfsl.converter import convert
         params = {}
-        params['split'] = model.magnet[templ].get('split', False)
-        params['show_plots'] = model.magnet[templ].get('plot', False)
+        rotor = model.magnet
+        params['split'] = rotor[templ].get('split', False)
+        params['show_plots'] = rotor[templ].get('plot', False)
         params['write_fsl'] = True
         params['airgap'] = -1.0
         pos = 'out' if model.external_rotor else 'in'
         params['part'] = ('rotor', pos)
         logger.info("Conv rotor from %s",
-                    model.magnet[templ]['name'])
-        conv = convert(model.magnet[templ]['name'], **params)
+                    rotor[templ]['name'])
+        conv = convert(rotor[templ]['name'], **params)
         model.set_value('poles', int(conv.get('num_poles')))
         self.set_diameter_parameter(model, conv)
         if model.get('da2'):
@@ -328,9 +333,9 @@ class Builder:
             ag = (model.get('bore_diam') - model.get('da2')/1e3)/2
             model.set_value('airgap', ag)
 
-        model.magnet['dxf'] = dict(fsl=conv['fsl'])
-        self.fsl_magnet = True
-        del model.magnet[templ]
+        rotor['dxf'] = dict(fsl=conv['fsl'])
+        self.fsl_rotor = True
+        del rotor[templ]
 
     def render_rotor(self, magmodel, templ):
         fslcode = self.__render(magmodel, templ, magnet=True)
@@ -351,23 +356,25 @@ class Builder:
             if 'thcond' in model.stator:
                 fslcmds += [
                     '-- thermal properties in airgap',
-                    'ag_cond = 0.063',
-                    'thcap = 1007',
-                    'beta = math.pi*m.npols_gen/m.num_poles + m.zeroangl/180*math.pi',
-                    'xai, yai = pr2c((da1+da2)/4, beta)',
-                    'def_mat_therm(xai,yai,"cyan",1.19,ag_cond,thcap,1)',
-                    'xai, yai = pr2c((da1+da2)/4-ag/4, beta)',
-                    'def_mat_therm(xai,yai,"cyan",1.19,ag_cond,thcap,1)',
-                    'xai, yai = pr2c((da1+da2)/4+ag/4, beta)',
-                    'def_mat_therm(xai,yai,"cyan",1.19,ag_cond,thcap,1)',
+                    'if m.zeroangl ~= nil then',
+                    '   ag_cond = 0.063',
+                    '   thcap = 1007',
+                    '   beta = math.pi*m.npols_gen/m.num_poles + m.zeroangl/180*math.pi',
+                    '   xai, yai = pr2c((da1+da2)/4, beta)',
+                    '   def_mat_therm(xai,yai,"cyan",1.19,ag_cond,thcap,1)',
+                    '   xai, yai = pr2c((da1+da2)/4-ag/4, beta)',
+                    '   def_mat_therm(xai,yai,"cyan",1.19,ag_cond,thcap,1)',
+                    '   xai, yai = pr2c((da1+da2)/4+ag/4, beta)',
+                    '   def_mat_therm(xai,yai,"cyan",1.19,ag_cond,thcap,1)',
                     '',
-                    'state_of_problem("therm_static") -- thermic boundary conditions',
-                    'x1,y1 = pd2c(dy2/2,m.zeroangl)',
-                    'x2,y2 = pd2c(dy1/2,m.zeroangl)',
-                    'beta = 360*m.npols_gen/m.num_poles',
-                    'x3,y3 = pd2c(dy1/2,beta+m.zeroangl)',
-                    'x4,y4 = pd2c(dy2/2,beta+m.zeroangl)',
-                    'def_bcond_tp(x1,y1,x2,y2,x3,y3,x4,y4, 4)',
+                    '   state_of_problem("therm_static") -- thermic boundary conditions',
+                    '   x1,y1 = pd2c(dy2/2,m.zeroangl)',
+                    '   x2,y2 = pd2c(dy1/2,m.zeroangl)',
+                    '   beta = 360*m.npols_gen/m.num_poles',
+                    '   x3,y3 = pd2c(dy1/2,beta+m.zeroangl)',
+                    '   x4,y4 = pd2c(dy2/2,beta+m.zeroangl)',
+                    '   def_bcond_tp(x1,y1,x2,y2,x3,y3,x4,y4, 4)',
+                    'end',
                     'state_of_problem("mag_static")']
             return fslcmds
         return []
@@ -409,10 +416,10 @@ class Builder:
                 cond = condMat.find(windings['material'])
             if not cond:
                 raise FslBuilderError(
-                    'conductor material {} not found'.format(
-                        windings['material']))
+                    'conductor material {} not found in {}'.format(
+                        windings['material'], condMat))
             windings['cuconduct'] = cond['elconduct']
-            for k in ('thcond', 'thcap'):
+            for k in ('thcond', 'thcap', 'spmaweight'):
                 if k in cond:
                     windings[k] = cond[k]
 
@@ -433,7 +440,14 @@ class Builder:
         return []
 
     def create_gen_winding(self, model):
-        genwdg = self.__render(model, 'gen_winding')
+        try:
+            if model.winding['wire']['name'] == 'hairpin_winding':
+                model.winding['wire'].update(
+                    {"num_layers": model.winding["num_layers"]})
+                genwdg = self.__render(model.winding,
+                                       'gen_' + model.winding['wire'].get('name'))
+        except KeyError:  # not hairpin_winding
+            genwdg = self.__render(model, 'gen_winding')
         k = list({'leak_dist_wind',
                   'leak_evol_wind',
                   'leak_tooth_wind'}.intersection(model.winding))
@@ -480,6 +494,9 @@ class Builder:
         params['airgap'] = model.dxffile.get('airgap', 0.0)
         params['nodedist'] = model.dxffile.get('nodedist', 1)
         params['full_model'] = model.dxffile.get('full_model', False)
+        params['EESM'] = model.dxffile.get('type', 'PMSM') == 'EESM'
+        if params['EESM']:
+            model.rotor['EESM'] = {}
         conv = convert(dxfname, **params)
 
         model.set_value('poles', conv.get('num_poles'))
@@ -492,6 +509,7 @@ class Builder:
         if not hasattr(model, 'stator'):
             setattr(model, 'stator', {})
         model.stator['num_slots'] = conv.get('tot_num_slot')
+        model.stator['slot_area'] = conv.get('slot_area')
         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",
@@ -509,11 +527,17 @@ class Builder:
         if 'fsl_stator' in conv:
             self.fsl_stator = True
             model.stator['dxf'] = dict(fsl=conv['fsl_stator'])
-        if not hasattr(model, 'magnet'):
-            setattr(model, 'magnet', {})
-        if 'fsl_magnet' in conv:
-            self.fsl_magnet = True
-            model.magnet['dxf'] = dict(fsl=conv['fsl_magnet'])
+        if not (hasattr(model, 'magnet') or hasattr(model, 'rotor')):
+            if params['EESM']:
+                setattr(model, 'rotor', {})
+            else:
+                setattr(model, 'magnet', {})
+        if 'fsl_rotor' in conv:
+            self.fsl_rotor = True
+            if hasattr(model, 'magnet'):
+                model.magnet['dxf'] = dict(fsl=conv['fsl_rotor'])
+            if hasattr(model, 'rotor'):
+                model.rotor['dxf'] = dict(fsl=conv['fsl_rotor'])
 
     def create_model(self, model, magnets=[], condMat=[], ignore_material=False):
         magnetMat = {}
@@ -563,10 +587,9 @@ class Builder:
                         'remanenc', 1.2)
                     model['magnet']['relperm'] = magnetMat.get('relperm', 1.05)
                     model['magnet']['rlen'] = magnetMat.get('rlen', 1.0)
-                    for k in ('thcond', 'thcap'):
+                    for k in ('spmaweight', 'thcond', 'thcap'):
                         if k in magnetMat:
-                            model['magnet'][k] = magnetMat[k]
-
+                            model['magnet'][k+'_magnet'] = magnetMat[k]
                 rotor = (self.create_magnet(model) +
                          self.create_magnet_model(model))
                 if magnetMat:
@@ -714,11 +737,10 @@ class Builder:
         return self.__render(model, 'colorgrad')
 
     def mesh_airgap(self, model):
-        if ((self.fsl_stator and self.fsl_magnet) or
+        if ((self.fsl_stator and self.fsl_rotor) or
                 model.get('num_agnodes', 0)):
             return self.__render(model, 'mesh-airgap')
-        else:
-            return []
+        return []
 
     def create(self, model, sim, magnets=None, condMat=[]):
         "create model and analysis function"
@@ -768,10 +790,10 @@ class Builder:
 
         return (fslmodel + self.create_analysis(sim) +
                 ['save_model("close")'])
-    
-    def create_detailed_wire(self, params, templ): 
+
+    def create_detailed_wire(self, params, templ):
         return self.__render(params, templ)
-    
+
     def __render(self, model, templ, stator=False, magnet=False):
         if templ.split('.')[-1] in ('fsl', 'mako'):
             try:
@@ -791,7 +813,7 @@ class Builder:
             if stator:
                 self.fsl_stator = True
             if magnet:
-                self.fsl_magnet = True
+                self.fsl_rotor = True
 
         return template.render_unicode(model=model).split('\n')
 

femagtools/isa7.py

@@ -528,7 +528,7 @@ class Isa7(object):
 
     color = {1: [1.0, 0.0, 0.0],  # RED
              2: [0.0, 1.0, 0.0],  # GREEN
-             3: [1.0, 1.0, 0.0],  # DARKYELLOW
+             3: [1.0, 0.8, 0.0],  # DARKYELLOW
              4: [0.0, 0.5019607843137255, 1.0], # BLUE
              5: [0.9803921568627451, 0.0, 1.0], # MAGENTA
              6: [0.0, 1.0, 0.8235294117647058], # CYAN
@@ -749,7 +749,7 @@ class Isa7(object):
                                      for e in self.elements])
 
         for a in ('FC_RADIUS', 'pole_pairs', 'poles_sim',
-                  'delta_node_angle', 'speed',
+                  'layers', 'coil_span', 'delta_node_angle', 'speed',
                   'MAGN_TEMPERATURE', 'BR_TEMP_COEF',
                   'MA_SPEZ_WEIGHT', 'CU_SPEZ_WEIGHT'):
             v = getattr(reader, a, '')

femagtools/machine/effloss.py

@@ -323,6 +323,8 @@ def efficiency_losses_map(eecpars, u1, T, temp, n, npoints=(60, 40),
         plcu1 = m.iqd_plcu1(iqd[0], iqd[1], 2*np.pi*f1)
         plcu2 = m.iqd_plcu2(*iqd)
         tfric = m.tfric
+        logger.info("Iex %f %f",
+                    np.min(iqd[2]), np.max(iqd[2]))
     else:
         plfe1 = np.array(r['plfe1'])
         plfe2 = np.zeros(ntmesh.shape[1])

femagtools/machine/pm.py

@@ -141,6 +141,14 @@ class PmRelMachine(object):
     def rstat(self, w):
         """stator resistance
         """
+        if isinstance(self.zeta1, list):
+            logger.info("setup ac loss parameters...")
+            # polyfit from ac loss calculation
+            freq = w/2/np.pi
+            kr = self.zeta1[0]*freq**3 + self.zeta1[1]*freq**2 + \
+                self.zeta1[2]*freq + self.zeta1[3]
+            kr[kr<1] = 1.
+            return self.r1*(1 - self.kth1*(self.tcu1 - 20))*kr  # ref 20°C
         if self.skin_resistance is not None:
             return self.skin_resistance(self.r1, w, self.tcu1, kth=self.kth1)
 
@@ -947,7 +955,7 @@ class PmRelMachine(object):
         i1max: max. phase current (RMS)
         """
         r = dict(id=[], iq=[], uq=[], ud=[], u1=[], i1=[], T=[],
-                 beta=[], gamma=[], phi=[], cosphi=[], pmech=[], n=[])
+                 beta=[], gamma=[], phi=[], cosphi=[], pmech=[], n=[], type_op=[])
 
         if kwargs.get('i1max', 0):
             w1type, T = self.w1_imax_umax(kwargs['i1max'], u1max)
@@ -1032,6 +1040,7 @@ class PmRelMachine(object):
                 r['n'].append(nx)
                 r['T'].append(Tf)
 
+            r['type_op'] = list(betai1(iq, id))
             Pmax = 2*np.pi*n1*Tf
             for ns, nu, iv in zip(nstab[1:], speedrange[2:], interv):
                 # find id, iq, torque in fieldweakening range