femagtools 1.6.8__py3-none-any.whl → 1.7.0__py3-none-any.whl

femagtools/plot/char.py

@@ -3,15 +3,17 @@
 """
 import numpy as np
 import matplotlib.pyplot as plt
+from ..machine.utils import iqd
 import logging
 
 
-def mtpa(pmrel, i1max, title='', projection='', ax=0):
+def mtpa(pmrel, i1max, u1max=0, title='', projection='', ax=0):
     """create a line or surface plot with torque and mtpa curve
 
     Args:
         pmrel: pm machine object
         i1max: maximum phase current / A
+        u1max: maximum phase voltage / V (optional)
         title: plot title
         projection: projection to be used (surface if '3d')
         ax: axis to be used
@@ -46,19 +48,29 @@ def mtpa(pmrel, i1max, title='', projection='', ax=0):
 
         ax.set_aspect('equal')
         x, y = np.meshgrid(id, iq)
-        CS = ax.contour(x, y, torque_iqd, 6, colors='k')
-        ax.clabel(CS, fmt='%d', inline=1)
+        nlevels = 6
+        CST = ax.contour(x, y, torque_iqd, nlevels, colors='k')
+        ax.clabel(CST, fmt='%d Nm', inline=1)
+        if u1max > 0:
+            n_iqd = np.array([
+                [pmrel.w1_umax(u1max, iqx, idx)/np.pi/2/pmrel.p*60
+                 for idx in id] for iqx in iq])
+            CSN = ax.contour(x, y, n_iqd, nlevels, colors='k',
+                             linestyles='dotted')  # linestyles='dashed')
+            ax.clabel(CSN, fmt='%d rpm', inline=1)
 
         ax.set_xlabel('Id/A')
         ax.set_ylabel('Iq/A')
         ax.plot(iopt[1], iopt[0],
-                color='red', linewidth=2, label='MTPA: {0:5.0f} Nm'.format(
+                color='red', linewidth=2, label='MTPA: {0:.0f} Nm'.format(
                     np.max(iopt[2][-1])))
+        iqx, idx = np.array([iqd(b, i1max) for b in np.linspace(-np.pi/2, 0)]).T
+        ax.plot(idx, iqx, color='blue', label='I1max: {0:.0f} A'.format(i1max))
         ax.grid()
 
     if title:
         ax.set_title(title)
-    ax.legend()
+    ax.legend(loc='upper left', framealpha=1)
 
 
 def mtpv(pmrel, u1max, i1max, title='', projection='', ax=0):
@@ -176,7 +188,12 @@ def characteristics(char, title=''):
     axs[1, 0].set_xlabel("Speed / rpm")
     axs[1, 0].set_ylabel("Current / A")
     axs[1, 0].legend(loc='center left')
-    if 'beta' in char:
+    if 'iex' in char:
+        ax3 = axs[1, 0].twinx()
+        ax3.plot(n, char['iex'], 'C3-', label='Iexc')
+        ax3.set_ylabel("Iexc / A")
+        ax3.legend(loc='center right')
+    elif 'beta' in char:
         ax3 = axs[1, 0].twinx()
         ax3.plot(n, char['beta'], 'C3-', label='Beta')
         ax3.set_ylabel("Beta / °")
@@ -337,7 +354,7 @@ def _plot_contour(speed, torque, z, ax, title='', levels=[],
     return contf, xscale, yscale
 
 
-def efficiency_map(rmap, ax=0, title='Efficiency Map', clabel=True,
+def efficiency_map(rmap, ax=0, title='', clabel=True,
                    cmap='YlOrRd', levels=None, cbar=False):
     if ax == 0:
         fig, ax = plt.subplots(figsize=(12, 12))

femagtools/plot/forcedens.py

@@ -8,7 +8,7 @@
 """
 import numpy as np
 import matplotlib.pyplot as plt
-from matplotlib import cm
+from matplotlib import colormaps
 
 def _create_3d_axis():
     """creates a subplot with 3d projection if one does not already exist"""
@@ -24,7 +24,7 @@ def _create_3d_axis():
         plt.subplot(111, projection='3d')
 
 
-def _plot_surface(ax, x, y, z, labels, azim=None):
+def _plot_surface(ax, x, y, z, labels, azim=None, cmap=colormaps['viridis']):
     """helper function for surface plots"""
     # ax.tick_params(axis='both', which='major', pad=-3)
     assert np.size(x) > 1 and np.size(y) > 1 and np.size(z) > 1
@@ -34,7 +34,7 @@ def _plot_surface(ax, x, y, z, labels, azim=None):
     Z = np.ma.masked_invalid(z)
     ax.plot_surface(X, Y, Z,
                     rstride=1, cstride=1,
-                    cmap=cm.viridis, alpha=0.85,
+                    cmap=cmap,
                     vmin=np.nanmin(z), vmax=np.nanmax(z),
                     linewidth=0, antialiased=True)
 #                    edgecolor=(0, 0, 0, 0))
@@ -60,43 +60,70 @@ def forcedens(title, pos, fdens, ax=0):
     ax.set_ylabel('Force Density / kN/m²')
 
 
-def forcedens_surface(fdens, ax=0):
+def forcedens_surface(fdens, attr='FN', ax=0, cmap=colormaps['jet']):
     if ax == 0:
         _create_3d_axis()
         ax = plt.gca()
-    xpos = [p for p in fdens.positions[0]['X']]
-    ypos = [p['position'] for p in fdens.positions]
-    z = 1e-3*np.array([p['FN']
-                       for p in fdens.positions])
-    _plot_surface(ax, xpos, ypos, z,
-                  ('Rotor pos/°', 'Pos/°', 'F N / kN/m²'))
-
-
-def forcedens_fft(title, fdens, nharm=40, ax=0):
+    z = 1e-3*np.array(getattr(fdens, attr))
+    n, m = z.shape
+    x = np.linspace(0, 360, m)
+    y = np.linspace(0, 360, n)
+    _plot_surface(ax, x, y, z,
+                  ('Rotor pos/°', 'Pos/°', f'{attr} / kN/m²'),
+                  cmap=cmap)
+
+def forcedens_contour(fdens, attr='FN', ax=0, cmap=colormaps['jet']):
+    if ax == 0:
+        ax = plt.gca()
+    z = 1e-3*np.array(getattr(fdens, attr))
+    n, m = z.shape
+    x = np.linspace(0, 360, m)
+    y = np.linspace(0, 360, n)
+    cs = ax.contourf(x, y, z, cmap=cmap)
+    ax.set_xlabel('Rotor pos/°')
+    ax.set_ylabel('Pos/°')
+    if attr == 'FN':
+        ax.set_title('Radial Force Density')
+    elif attr == 'FT':
+        ax.set_title('Tangential Force Density')
+    cfig = ax.get_figure()
+    cbar = cfig.colorbar(cs, ax=ax)
+    cbar.ax.set_ylabel('kN/m²')
+
+def forcedens_fft(title, fdens, harmmax=(), #(200, 40),
+                  cmap=colormaps['YlOrBr'],
+                  ax=0):
     """plot force densities FFT
     Args:
       title: plot title
       fdens: force density object
-      nharm: (int) num harmonics
+      cmap: colormap
     """
     if ax == 0:
         ax = plt.axes(projection="3d")
 
-    F = 1e-3*fdens.fft(nharm)['fn_harm']['amplitude']
-    fmin = 0.2
-    num_bars = F.shape[0] + 1
-    _xx, _yy = np.meshgrid(np.arange(1, num_bars),
-                           np.arange(1, num_bars))
-    z_size = F[F > fmin]
-    x_pos, y_pos = _xx[F > fmin], _yy[F > fmin]
-    z_pos = np.zeros_like(z_size)
-    x_size = 2
-    y_size = 2
-
-    ax.bar3d(x_pos, y_pos, z_pos, x_size, y_size, z_size)
-    #ax.view_init(azim=120)
-    ax.set_xlim(0, num_bars+1)
-    ax.set_ylim(0, num_bars+1)
+    FN = 1e-3*fdens.fft()['fn_harm']['amplitude']
+    if harmmax:
+        shape = np.min((harmmax, np.shape(FN)), axis=0)
+    else:
+        shape = np.shape(FN)
+
+    _x = np.arange(shape[0])+1
+    _y = np.arange(shape[1])+1
+    _xx, _yy = np.meshgrid(_x, _y)
+    x, y = _xx.ravel(), _yy.ravel()
+    top = np.ravel(FN[:shape[0],:shape[1]])
+    bottom = np.zeros_like(top)
+    maxy = max(x[-1], y[-1])
+    width, depth = 5*x[-1]/maxy, 5*y[-1]/maxy
+
+    min_height = np.min(top)
+    max_height = np.max(top)
+    rgba = [cmap((k-min_height)/max_height) for k in top]
+    ax.bar3d(x, y, bottom, width, depth, top, color=rgba)
+
+    ax.set_xlim(0, shape[0]+1)
+    ax.set_ylim(0, shape[1]+1)
     ax.set_title(title)
     ax.set_xlabel('M')
     ax.set_ylabel('N')

femagtools/plot/mcv.py

@@ -70,7 +70,8 @@ def felosses(losses, coeffs, title='', log=True, ax=0):
     Args:
       losses: dict with f, B, pfe values
       coeffs: list with steinmetz (cw, alpha, beta) or
-              jordan (cw, alpha, ch, beta, gamma) coeffs
+              jordan (cw, alpha, ch, beta, gamma)  or
+              bertotti (cw, alpha, cw, ce) coeffs
       title: title string
       log: log scale for x and y axes if True
 
@@ -89,6 +90,8 @@ def felosses(losses, coeffs, title='', log=True, ax=0):
         if f > 0:
             if len(coeffs) == 5:
                 ax.plot(B, lc.pfe_jordan(f, B, *coeffs, fo=fo, Bo=Bo))
+            if len(coeffs) == 4:
+                ax.plot(B, lc.pfe_bertotti(f, B, *coeffs))
             elif len(coeffs) == 3:
                 ax.plot(B, lc.pfe_steinmetz(f, B, *coeffs, fo=fo, Bo=Bo))
             ax.plot(losses['B'][:len(pfe)], pfe,

femagtools/plot/phasor.py

@@ -141,6 +141,11 @@ def phasor(bch, ax=0):
     xd = w1*bch.dqPar['ld'][-1]
     xq = w1*bch.dqPar['lq'][-1]
     r1 = bch.machine['r1']
-    i1beta_phasor(bch.dqPar['up'][-1],
+    if np.isscalar(bch.dqPar['up']):
+        up = bch.dqPar['up']
+    else:
+        up = bch.dqPar['up'][-1]
+
+    i1beta_phasor(up,
                   bch.dqPar['i1'][-1], bch.dqPar['beta'][-1],
                   r1, xd, xq, ax)

femagtools/poc.py

@@ -11,11 +11,12 @@
 
 
 class Poc:
-    def __init__(self, arg, parameters=dict()):
+    def __init__(self, arg=999, parameters={}, current_angles=[]):
         """initialize this object from a set of parameters or a file
 
         Args:
-            arg filename or pole_pitch
+            arg: filename or pole_pitch
+            current_angles: list of current angles
         """
         self.skew_angle = 0.0
         self.num_skew_steps = 0
@@ -33,15 +34,21 @@ class Poc:
                                               [str(i+1) for i in range(3)])
             b = parameters.get('offset', 0)
             num_winding = len(self.key_winding)
-            self.phi_voltage_winding = parameters.get('phi_voltage_winding',
-                                                      [b+i*360/num_winding
-                                                       for i in range(num_winding)])
+            if current_angles:
+                self.phi_voltage_winding = current_angles
+            else:
+                self.phi_voltage_winding = parameters.get('phi_voltage_winding',
+                                                          [b+i*360/num_winding
+                                                           for i in range(num_winding)])
 
     def __setattr__(self, name, val):
         self.__dict__[name] = val  # this will create the attribute name
 
-    def filename(self):
-        prefix = self.pocType
+    def filename(self, modelName=''):
+        if modelName:
+            prefix = modelName
+        else:
+            prefix = self.pocType
         if prefix == 'Function':
             prefix = self.shape_current
         return '{}_{}p.poc'.format(
@@ -125,7 +132,7 @@ class Poc:
         else:
             self.shape_current=self.pocType
             self.pocType='Function'
-        
+
         try:
             self.skew_angle=float(pocfile.readline())
             self.num_skew_steps=int(pocfile.readline())
@@ -155,7 +162,7 @@ class HspPoc(Poc):
             pocType='hsp', harmonic_id=harm,
             func_current=amp, func_phi=phi))
 
-        
+
 def curr_har(x, n, A, phi ):
     "return fourier sum"
     import numpy as np
@@ -169,7 +176,7 @@ def curr_fun(phi, A, phi_wdg ):
     "return curr"
     import numpy as np
     if np.asarray(A).shape == 1:
-        return [[np.asarray(phi)+x, A] for x in phi_wdg] 
+        return [[np.asarray(phi)+x, A] for x in phi_wdg]
     return [phi, A]
 
 

femagtools/templates/cogg_calc.mako

@@ -18,7 +18,6 @@ m.line            =    0
 m.two_pole_wi     =    2*m.pole_width
 m.range_x         =    m.two_pole_wi
 m.range_y         =    0.0
-
 m.fc_force_points =  5
 m.fcpx_mm1        =   m.npols_gen*m.pole_width +1.0
 m.fcpy_mm1        =   -ag/2
@@ -39,7 +38,14 @@ m.fcpy_mm5        =   m.fcpy_mm1
 % endif
 
 m.nu_skew_steps   =    ${model.get('num_skew_steps', 0)}
+% if model.get('noload_ex_cur', 0): 
+m.nloa_ex_cur = ${model.get('noload_ex_cur',0)}
+% endif 
+% if model.get("noload_ex_cur", 0): 
+m.magn_temp = 20   -- dummy parameter
+% else: 
 m.magn_temp       =    ${model.get('magn_temp')}
+% endif
 m.fc_mult_move_type =  1.0 --  Type of move path in air gap
 m.nu_move_steps   =    ${model.get('num_move_steps', 49)}
 

femagtools/templates/displ_stator_rotor.mako

@@ -0,0 +1,33 @@
+-- calculate static or dynamic eccentricity
+-- NOTE: requires full model
+--
+static = 0
+dynamic = 1
+m.disp_info   =  ${model.get('type', 'static')}  -- Displace: Stator: 0, Rotor: 1
+
+da1 = ${model.get('bore_diam',0)*1e3}
+ag = ${model.get('airgap',0)*1e3}
+
+if m.b_min == 0 then -- move inside: outer stator
+  if m.disp_info == static then
+    m.fc_radius = da1/2 - ag/6
+  else
+    m.fc_radius = da1/2 - 5*ag/6
+  end
+else
+  if m.disp_info == static then
+    m.fc_radius = da1/2 + ag/6
+  else
+    m.fc_radius = da1/2 + 5*ag/6
+  end
+end
+m.disp_radius =  ${model.get('ecc', 0)*1e3}  -- Displacement: r [mm]
+m.disp_phi    =  0.0         -- Displacement: phi [Degr]
+
+pre_models("Displ_Stat/Rot")
+-- must correct fc_radius
+if m.b_min == 0 then -- move inside: outer stator
+  m.fc_radius = da1/2 - ag/2
+else
+  m.fc_radius = da1/2 + ag/2
+end

femagtools/templates/fieldcalc.mako

@@ -6,29 +6,23 @@ du_u0=1e-3
 
 get_wdg_keys("wkeys")
 
--- check the machine type PSM/ESM
-magn = get_dev_data("magn_remanence")
-
-if magn > 0 then -- IPM
-  -- set cur to zero 
-  for i = 1, #wkeys do   
-    def_curr_wdg(wkeys[i],0, 0) 
-  end 
-else -- ESM
-  for i = 1, #wkeys - 1 do   
-    def_curr_wdg(wkeys[i],0, 0) 
-  end 
+-- set cur to zero 
+for i = 1, #wkeys do   
+  def_curr_wdg(wkeys[i],0, 0) 
 end 
-   
+% if model.get('noload_ex_cur', 0):
+def_curr_wdg(wkeys[#wkeys], ${model['noload_ex_cur']}, 0) 
+% endif
+
 calc_field_single(maxit, reset, du_u0)
 
-post_models("induct(x)","b")    -- Calculate field distribution
+post_models("induct(x)","b_airgap")    -- Calculate field distribution
 
   data=io.open("bag.dat","w")              -- Output in data file
-  N = table.getn(b)                             -- Number of elements in array
+  N = table.getn(b_airgap)                             -- Number of elements in array
   i = 1
   repeat
-    data:write(string.format("%g %g %g\n",b[i],b[i+1],b[i+2]))
+    data:write(string.format("%g %g %g\n",b_airgap[i],b_airgap[i+1],b_airgap[i+2]))
     i = i+3
   until i>=N
   io.close(data)

femagtools/templates/pm_sym_f_cur.mako

@@ -1,7 +1,7 @@
 --
 -- SM Simulation with excitation current
 --
-m.wdcon           = ${model.get('wdgcon', 1)}   --    Connection: 0=open, 1=star, 2=delta
+m.wdcon           = ${model.get('wdgcon', 0)}   --    Connection: 0=open, 1=star, 2=delta
 m.move_action     = ${model.get('move_action', 0)}
 m.speed           = ${model.get('speed')*60}
 m.skew_angle      = ${model.get('skew_angle',0)}

femagtools/tks.py

@@ -56,7 +56,6 @@ class Reader(object):
         self.rho = 7.6
         self.losses = dict(f=[], B=[])
         pfe = []
-        b_pfe = []
 
         # filecontent is used
         if filecontent:
@@ -92,7 +91,6 @@ class Reader(object):
                     b, p = readlist(content[i+1:])
                     self.losses['f'].append(fxref)
                     self.losses['B'].append(b)
-                    b_pfe.append(b)
                     pfe.append(p)
         logger.info("%s Bmax %3.2f", filename, max(self.curve[0]['bi']))
 
@@ -126,19 +124,15 @@ class Reader(object):
 
             self.losses['Bo'] = self.Bo
             self.losses['fo'] = self.fo
-
-            z = lc.fit_bertotti0(self.losses['f'], b_pfe, pfe)
-            if z.shape[0] == 4:
-                self.bertotti = {'ch': z[0], 'alpha': z[1], 'cw': z[2], 'ce': z[3]}
-            else:
-                self.bertotti = {'ch': z[0], 'alpha': 2.0, 'cw': z[1], 'ce': z[2]}
+            z = lc.fit_bertotti(self.losses['f'],
+                                self.losses['B'], pfe)
+            self.bertotti = {'ch': z[0], 'alpha': 2.0, 'cw': z[1], 'ce': z[2]}
             logger.info("Bertotti loss coeffs %s", z)
 
             # must normalize pfe matrix:
             (self.losses['B'],
              self.losses['pfe']) = femagtools.mcv.norm_pfe(
                  self.losses['B'], pfe)
-        self.b_pfe = b_pfe
         self.pfe = pfe
 
 

{femagtools-1.6.8.dist-info → femagtools-1.7.0.dist-info}/LICENSE

@@ -1,4 +1,5 @@
 Copyright (c) 2016-2023, Semafor Informatik & Energie AG, Basel
+Copyright (c) 2023-2024, Gamma Technology LLC
 All rights reserved.
 
 Redistribution and use in source and binary forms, with or without

{femagtools-1.6.8.dist-info → femagtools-1.7.0.dist-info}/METADATA

@@ -1,9 +1,10 @@
 Metadata-Version: 2.1
 Name: femagtools
-Version: 1.6.8
+Version: 1.7.0
 Summary: Python API for FEMAG
 Author-email: Ronald Tanner <tar@semafor.ch>, Dapu Zhang <dzhang@gtisoft.com>, Beat Holm <hob@semafor.ch>, Günther Amsler <amg@semafor.ch>, Nicolas Mauchle <mau@semafor.ch>
 License: Copyright (c) 2016-2023, Semafor Informatik & Energie AG, Basel
+        Copyright (c) 2023-2024, Gamma Technology LLC
         All rights reserved.
         
         Redistribution and use in source and binary forms, with or without
@@ -42,6 +43,8 @@ Requires-Dist: mako
 Requires-Dist: six
 Requires-Dist: lmfit
 Requires-Dist: netCDF4 >=1.6.5
+Provides-Extra: all
+Requires-Dist: femagtools[dxfsl,meshio,mplot,svgfsl,test,vtk,zmq] ; extra == 'all'
 Provides-Extra: dxfsl
 Requires-Dist: dxfgrabber ; extra == 'dxfsl'
 Requires-Dist: networkx ; extra == 'dxfsl'
@@ -81,14 +84,14 @@ Femagtools is an Open-Source Python-API for FEMAG offering following features:
 * calculate machine characteristics by using analytic machine models
 * execute parameter studies and multi-objective optimization
 
-The package can be used with Python 3.x on Linux, MacOS or Windows and is hosted on github: <https://github.com/SEMAFORInformatik/femagtools/>` where also many examples can be found in the examples directory. Contributions and feedback to this project are highly welcome.
+The package can be used with Python 3.x on Linux, MacOS or Windows and is hosted on github: <https://github.com/SEMAFORInformatik/femagtools/> where also many examples can be found in the examples directory. Contributions and feedback to this project are highly welcome.
 
 The installation can be done in the usual ways with pip:
 
 ```
-pip install femagtools[all]
+pip install 'femagtools[all]'
 ```
-([all] pulls in all optional dependencies.
+`[all]` pulls in all optional dependencies. Up-to-date information about optional dependencies can be found in the [pyproject.toml](pyproject.toml) file under `[project.optional-dependencies]`.
 
 For details see the documentation <http://docs.semafor.ch/femagtools>