femagtools 1.7.9__py3-none-any.whl → 1.8.1__py3-none-any.whl

femagtools/templates/prepare_thermal.mako

@@ -1,7 +1,7 @@
 -- prepare thermal model
 
 save_model('cont')
-
+state_of_problem("therm_static")
 function get_xy(elkey)
   -- get element xy
   local x, y
@@ -40,12 +40,13 @@ shaft_spel = {} ---- lightgrey, 16
 magnet_spel = {} -- black, 8
 outer_air = {} -- darkred, 10
 inner_air = {} -- lightgrey, 16
-rotor_air = {} 
+rotor_air = {}
 stator_air = {} -- black, 8
 airgap = {} -- white, 8
 magnet_pocket_spel = {}  -- red, 1
 non_uniform_airgap = {} -- white, 8
-
+slot_opening = {}
+slot_sreg_exist = 0
 -- search all components
 for i =1,  #spel_keys do
   -- loop over all spel keys
@@ -73,7 +74,7 @@ for i =1,  #spel_keys do
 %if model.get('htc_outer', 0):
               def_heat_transfer(xc,yc,yellow,${model['htc_outer']}, 1.0)
 %endif
-              color_spel(spel_keys[i], 10) 
+              color_spel(spel_keys[i], 10)
 
            elseif rc > dy2/2 and rc < da2/2 then
               -- rotor air
@@ -86,11 +87,18 @@ for i =1,  #spel_keys do
               -- stator air
               table.insert(stator_air, spel_keys[i])
               def_mat_therm(xc, yc, 'yellow', 1.19,0.15,1007, 1)
-              color_spel(spel_keys[i], 8) 
-
+              color_spel(spel_keys[i], 8)
+              key_exist_sl = get_sreg_key(xc, yc)
+              if key_exist_sl <= 0 and slot_sreg_exist == 0 then
+                  def_new_sreg(xc, yc, 'Slot', 'yellow')
+                  slot_sreg_exist = 1
+              else
+                  add_to_sreg(xc, yc, 'Slot')
+              end
            elseif rc > da2/2 and rc < da1/2 then
               table.insert(airgap, spel_keys[i])
-              def_mat_therm(xc, yc, 'yellow', 1.19,0.063,1007, 1)
+              def_mat_therm(xc, yc, 'yellow', 1.19,0.15,1007, 1)
+              add_to_sreg(xc, yc, 'Slot')
               color_spel(spel_keys[i], 7) -- white air
 
 --           elseif rc > dy2/2 then
@@ -104,14 +112,6 @@ for i =1,  #spel_keys do
               if x0_shaft ~= nil and x0_shaft ~= 0.0 then
                   table.insert(shaft_spel, spel_keys[i])
                   def_mat_therm(xc, yc, 'lightgrey', shaft_density, shaft_thcond, shaft_thcap, 1)
-              else
-                  -- is inner air
-                  table.insert(inner_air, spel_keys[i])
-%if model.get('htc_inner', 0):
-                  def_heat_transfer(xc,yc,yellow,${model['htc_inner']}, 1.0)
-%endif
-                  def_new_sreg(xc,yc, 'INAR', 'yellow')
-                  color_spel(spel_keys[i], 16) -- light grey inner air
               end
            end
         end
@@ -139,81 +139,219 @@ end
 
 function is_exist(arr, val)
    local i
-   if #arr == 0 then 
-      return false 
-   end 
-   for i = 1, #arr do 
-      if arr[i] == val then 
+   if #arr == 0 then
+      return false
+   end
+   for i = 1, #arr do
+      if arr[i] == val then
          return true
-      end 
-   end 
+      end
+   end
    return false
-end 
--- identify magnet air pocket 
+end
+-- identify magnet air pocket
 
-for i = 1, #magnet_spel do 
+for i = 1, #magnet_spel do
     bkeys = get_spel_data ( "bndkeys", magnet_spel[i])
-    for j = 1, #bkeys do 
+    for j = 1, #bkeys do
          elks = get_node_data ( "elkeys", bkeys[j] )
-         for k = 1, #elks do 
+         for k = 1, #elks do
             p1, p2 = get_elem_data("perm", elks[k])
             sek = get_elem_data("sekey", elks[k])
-            if p1 == 1.0 then 
-                  if (not is_exist(magnet_pocket_spel, sek)) and (not is_exist(airgap, sek)) then  
+            if p1 == 1.0 then
+                  if (not is_exist(magnet_pocket_spel, sek)) and (not is_exist(airgap, sek)) then
                      table.insert(magnet_pocket_spel, sek)
                      els = get_spel_data('elkeys', sek)
                      xc, yc = get_xy(els[1])
                      def_mat_therm(xc, yc, 'red', 1.12,0.026,1007, 1)
                      --table.remove(rotor_air, sek)
                      color_spel(sek, 1)
-                  end 
-            end 
-         end 
-    end 
-end 
+                  end
+            end
+         end
+    end
+end
 
 
 
-for i = 1, #rotor_air do 
-  if rotor_air[i] ~= nil then 
+for i = 1, #rotor_air do
+  if rotor_air[i] ~= nil then
     bkeys = get_spel_data ( "bndkeys", rotor_air[i])
-    for j = 1, #bkeys do 
+    for j = 1, #bkeys do
       x, y = get_node_data ( "xy", bkeys[j] )
-      r, phi = c2pd(x, y) 
-      if math.abs(r - da2/2) < 1e-5 then 
-        if not is_exist(non_uniform_airgap,rotor_air[i]) then 
+      r, phi = c2pd(x, y)
+      if math.abs(r - da2/2) < 1e-5 then
+        if not is_exist(non_uniform_airgap,rotor_air[i]) then
           table.insert(non_uniform_airgap, rotor_air[i])
           els = get_spel_data('elkeys', rotor_air[i])
           xc, yc = get_xy(els[1])
-          def_mat_therm(xc, yc, 'yellow', 1.19,0.063,1007, 1)
+          def_mat_therm(xc, yc, 'yellow', 1.19,1.15,1007, 1)
+          add_to_sreg(xc, yc, 'Slot')
           color_spel(rotor_air[i], 7)
-        end 
-      end 
-    end 
-  end 
-end 
+        end
+      end
+    end
+  end
+end
 
-for i = 1, #rotor_air do 
-  if rotor_air[i] ~= nil then 
-    if is_exist(non_uniform_airgap, rotor_air[i]) then 
+for i = 1, #rotor_air do
+  if rotor_air[i] ~= nil then
+    if is_exist(non_uniform_airgap, rotor_air[i]) then
       rotor_air[i] = nil
-    end    
-  end 
-end 
+    end
+  end
+end
 
-for i = 1, #rotor_air do 
-   if is_exist(magnet_pocket_spel, rotor_air[i]) then 
+for i = 1, #rotor_air do
+   if is_exist(magnet_pocket_spel, rotor_air[i]) then
      rotor_air[i] = nil
-   end 
- end 
- 
+   end
+ end
+
 
-for i = 1, #rotor_air do 
-   if rotor_air[i] ~= nil then 
+for i = 1, #rotor_air do
+   if rotor_air[i] ~= nil then
       els = get_spel_data('elkeys', rotor_air[i])
       xc, yc = get_xy(els[1])
       def_mat_therm(xc, yc, 'red', 1.12,0.026,1007, 1)
-   end 
-end 
+   end
+end
+
 
 save_metafile(model..'.ps')
+
+
+
+---------------------------------------------
+-- material slot insulation
+---------------------------------------------
+
+function get_boundary_node(num_slots)
+
+   local ctr
+   -- get the number of winding
+   get_wdg_keys("wkeys")
+   -- search subregions
+   local srkey_sl = {}
+   ctr = 1
+   for i = 1, #wkeys do
+      tmpksl = get_wdg_data("srkeys", wkeys[i])
+      for j = 1, #tmpksl do
+         srkey_sl[ctr] = tmpksl[j]
+         ctr = ctr + 1
+      end
+   end
+   -- search superelements
+   local sekey_sl = {}
+   ctr = 1
+   for k = 1, #srkey_sl do
+      sl_se = get_sreg_data("sekeys", srkey_sl[k])
+      for kk =1, #sl_se do
+         sl_nd = get_spel_data("bndkeys", sl_se[kk])
+         x, y = get_node_data("xy", sl_nd[1])
+         r, phi = c2pd(x, y)
+         if phi < 360/num_slots then
+            sekey_sl[ctr] = sl_se[kk]
+            ctr = ctr + 1
+         end
+      end
+   end
+
+   local xn  = {}
+   local yn = {}
+   local node = {}
+   local bndnodes = {}
+   local bnd_unique = {}
+
+   ctr = 1
+
+   for i = 1, #sekey_sl do
+      bnd = get_spel_data("bndkeys", sekey_sl[i])
+
+      for j = 1, #bnd do
+         bndnodes[ctr] = bnd[j]
+         ctr = ctr + 1
+         if bnd_unique[bnd[j]] == true then
+            bnd_unique[bnd[j]] = false
+         else
+            bnd_unique[bnd[j]] = true
+         end
+      end
+   end
+
+   ctr =  1
+   for j = 1, #bndnodes do
+     x, y = get_node_data("xy", bndnodes[j])
+
+     r, phi = c2pd(x, y)
+     if (phi < 360/num_slots/2 - 0.05) and (bnd_unique[bndnodes[j]] == true)  then
+
+        node[ctr] = bndnodes[j]
+        xn[ctr] = x
+        yn[ctr] = y
+        ctr = ctr + 1
+       end
+   end
+
+
+   local indx = {1, math.floor(#node/4), math.floor(#node/2), math.floor(#node/4*3), #node}
+
+   local x_new = {}
+   local y_new = {}
+
+   for i = 1, 10 do
+      x_new[i] = 0
+      y_new[i] = 0
+   end
+
+   for i = 1, #indx do
+      x_new[i] = xn[indx[i]]
+      y_new[i] = yn[indx[i]]
+      r, phi = c2pd(x_new[i], y_new[i])
+      x_new[10 - (i-1)], y_new[10 - (i-1)] = pd2c(r, 360/num_slots - phi)
+   end
+
+   local rn = {}
+   local phin = {}
+
+   for i = 1, #x_new do
+      rn[i], phin[i] = c2pd(x_new[i], y_new[i])
+   end
+   return rn, phin
+end
+
+
+rn, phin = get_boundary_node(m.tot_num_slot)
+thickness = ${model.windings.get('slot_insulation_thickness', 0.15e-3)*1e3}
+if thickness == 0.0 then 
+   thickness = 0.15
+end 
+conductivity = ${model.get('slot_insul_cond', 0.31)}
+
+for i = 1,m.num_sl_gen do
+x = {}
+y = {}
+
+for j = 1, #rn do
+   x[j], y[j] = pd2c(rn[j], phin[j] + (i-1)*360/m.tot_num_slot)
+   point(x[j], y[j],"black","x")
+end
+
+def_insulation_by_nodechain(thickness,conductivity,
+      x[1],y[1],
+      x[2],y[2],
+      x[3],y[3],
+      x[4],y[4],
+      x[5],y[5],
+      x[6],y[6],
+      x[7],y[7],
+      x[8],y[8],
+      x[9],y[9],
+      x[10],y[10],
+      x[1], y[1]
+      )
+
+end
+save_metafile('insulation.ps')
+
+state_of_problem("mag_static")

femagtools/windings.py

@@ -32,18 +32,18 @@ def q1q2yk(Q, p, m, l=1):
     """returns q1, q2, Yk, Qb"""
     t = num_basic_windings(Q, p, l)
     Qb = Q//t
-    qqb = Qb if l == 2 else Qb//2
+    qbb = Qb if l==2 else Qb//2
     pb = p//t
-    if qqb//m % 2:  # odd
-        q2 = (qqb + m)//(2*m) - 1
+    if qbb//m % 2:  # odd
+        q2 = (qbb + m)//(2*m) - 1
         q1 = q2 + 1
     else:
-        q2 = (qqb)//(2*m)
+        q2 = (qbb)//(2*m)
         q1 = q2
     n = 1
-    while (n*qqb + 1) % pb:
+    while (n*qbb + 1) % pb:
         n += 1
-    Yk = (n*qqb + 1)//pb
+    Yk = (n*qbb + 1)//pb
     return q1, q2, Yk, Qb
 
 
@@ -118,22 +118,22 @@ class Winding(object):
         self.yd = coilwidth
 
         q1, q2, Yk, Qb = q1q2yk(self.Q, self.p, self.m, self.l)
+        j = 2 if layers == 1 else 1
         k1 = [(q1 + q2)*i for i in range(self.m)]
         k2 = (q1*(self.m+1) + q2*(self.m-1))//2
-        j = 2 if layers == 1 else 1
         pos = [[(j*Yk*(k + n)) % Qb
                 for n in range(q1)] for k in k1]
         neg = [[j*Yk*(k + n + k2) % Qb
                 for n in range(q2)] for k in k1]
-        if layers > 1:
+        if self.l > 1:
             slots = [sorted([(k, 1, 1)
                              for k in p] + [(k, -1, 1)
                                             for k in n])
                      for n, p in zip(neg, pos)]
             for i, p in enumerate(slots):
                 slots[i] = sorted(slots[i] +
-                                  [((k[0]+coilwidth) % Qb, -k[1], 0)
-                                   for k in slots[i]], key=lambda s: s[0])
+                                  [((k[0]+self.yd) % Qb, -k[1], 0)
+                                   for k in slots[i]])
         else:
             if (coilwidth + 1) % 2:
                 coilwidth += 1
@@ -147,10 +147,6 @@ class Winding(object):
             slots = [sorted([(k, 1, 1) for k in p] + [(k, -1, 1) for k in n])
                      for n, p in zip(xneg, xpos)]
 
-        if self.m > 3:  # TODO: check this hack
-            slots = slots[:self.m//2] + [[((k[0]+1) % self.Q, k[1], k[2])
-                                          for k in s]
-                                         for s in slots[:self.m//2]]
         taus = 360/self.Q
         self.windings = {
             i+1:  dict(dir=[k[1] for k in s],
@@ -183,6 +179,10 @@ class Winding(object):
             nue = n
         else:
             nue = self.kw_order(n)
+        #if q1 == q2:  # integral slot winding
+        #    q = self.Q/2/self.m/self.p
+        #    nuep = nue/self.p
+        #    return np.sin(nuep*np.pi/2/self.m)/q/np.sin(nuep*np.pi/2/self.m/q)
         k = 2 if self.l == 1 else 1
         a = nue*k*np.pi/self.Q*Yk
         t = self.Q//Qb
@@ -244,14 +244,19 @@ class Winding(object):
 
     def slots(self, key):
         """returns slot indexes of winding key"""
-        ngen = self.m*self.Q//np.gcd(self.Q, self.m*2*self.p)
         taus = 360/self.Q
-
-        dim = int(self.l*ngen/self.m)
-        slots = [(round((x-taus/2)/taus) + ngen*n) % self.Q + 1
-                 for n in range(self.Q//ngen)
+        t = num_basic_windings(self.Q, self.p, self.l)
+        Qb = self.Q//t
+        dim = self.l*Qb//self.m
+        ngen = t
+        qgen = Qb
+        if dim > len(self.windings[key]['PHI']):
+            qgen = qgen//2
+            ngen = 2*t
+        slots = [(round((x-taus/2)/taus) + qgen*n) % self.Q + 1
+                 for n in range(ngen)
                  for x in self.windings[key]['PHI'][:dim]]
-        return np.array(slots).reshape((np.gcd(self.Q, self.p), -1))
+        return np.array(slots).reshape(t, -1)
 
     def axis(self, k=1):
         """returns axis angle of winding 1 in mechanical system"""

{femagtools-1.7.9.dist-info → femagtools-1.8.1.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: femagtools
-Version: 1.7.9
+Version: 1.8.1
 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
@@ -37,33 +37,33 @@ Classifier: Topic :: Scientific/Engineering
 Requires-Python: >=3.7
 Description-Content-Type: text/markdown
 License-File: LICENSE
-Requires-Dist: numpy <=1.26.4
+Requires-Dist: numpy<=1.26.4
 Requires-Dist: scipy
 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'
+Requires-Dist: netCDF4>=1.6.5
+Provides-Extra: rdp
+Requires-Dist: rdp; extra == "rdp"
 Provides-Extra: dxfsl
-Requires-Dist: dxfgrabber ; extra == 'dxfsl'
-Requires-Dist: networkx ; extra == 'dxfsl'
-Requires-Dist: rdp ; extra == 'dxfsl'
-Provides-Extra: meshio
-Requires-Dist: meshio ; extra == 'meshio'
-Provides-Extra: mplot
-Requires-Dist: matplotlib ; extra == 'mplot'
+Requires-Dist: dxfgrabber; extra == "dxfsl"
+Requires-Dist: networkx; extra == "dxfsl"
 Provides-Extra: svgfsl
-Requires-Dist: dxfgrabber ; extra == 'svgfsl'
-Requires-Dist: networkx ; extra == 'svgfsl'
-Requires-Dist: lxml ; extra == 'svgfsl'
-Requires-Dist: rdp ; extra == 'svgfsl'
-Provides-Extra: test
-Requires-Dist: pytest ; extra == 'test'
+Requires-Dist: dxfgrabber; extra == "svgfsl"
+Requires-Dist: networkx; extra == "svgfsl"
+Requires-Dist: lxml; extra == "svgfsl"
+Provides-Extra: mplot
+Requires-Dist: matplotlib; extra == "mplot"
+Provides-Extra: meshio
+Requires-Dist: meshio; extra == "meshio"
 Provides-Extra: vtk
-Requires-Dist: vtk ; extra == 'vtk'
+Requires-Dist: vtk; extra == "vtk"
 Provides-Extra: zmq
-Requires-Dist: pyzmq ; extra == 'zmq'
+Requires-Dist: pyzmq; extra == "zmq"
+Provides-Extra: test
+Requires-Dist: pytest; extra == "test"
+Provides-Extra: all
+Requires-Dist: femagtools[dxfsl,meshio,mplot,svgfsl,test,vtk,zmq]; extra == "all"
 
 
 # Introduction to Femagtools