femagtools 1.7.8__py3-none-any.whl → 1.8.0__py3-none-any.whl

femagtools/templates/magnetSector.mako

@@ -68,55 +68,40 @@ for i = 0, m.npols_gen-1 do
     end
 end
 %endif
-%if model.get('thcond', 0) and model.get('thcap', 0):
-beta = math.pi/m.num_poles + m.zeroangl/180*math.pi
-xrb,yrb = pr2c(m.yoke_rad+0.1, beta)  -- rotor lamination
-thcond = 24 -- ${model['thcond']}
-thcap = 480 -- ${model['thcap']}
-def_mat_therm(xrb,yrb,'blue',7700,thcond,thcap,1)
---def_mat_therm(m.yoke_rad/2,0.1,'blue',7700,thcond,thcap,1) -- Shaft
+%if model.get('thcond', 0) and model.get('thcap', 0) and model.get('density', 0):
+if m.shaft_rad == nil then 
+    m.shaft_rad = m.condshaft_r
+end 
+if m.shaft_rad > dy2/2 then 
+    m.shaft_rad = dy2/2
+end 
+beta = math.pi/m.num_poles
+rotor_thcond = ${model['thcond']}
+rotor_thcap = ${model['thcap']}
+rotor_density = ${model.get('density')}
 
-rm = m.rotor_rad - m.magn_height/2
-thcond = 8
-thcap = 440
-for i = 1,m.npols_gen do -- Magnets
-  alfa = (2*i-1) * beta
-  xmx,ymx = pr2c(rm,alfa-beta/2)
-  def_mat_therm(xmx,ymx,darkgreen-i%2,7500,thcond,thcap,1)
+%if model.get('thcond_shaft', 0) and model.get('thcap_shaft', 0):
+if m.shaft_rad < m.yoke_rad then
+   shaft_thcond = ${model['thcond_shaft']}
+   shaft_thcap = ${model['thcap_shaft']}
+   shaft_density = ${model['spmaweight_shaft']*1e3}
+   r_shaft = (m.shaft_rad + m.yoke_rad)/2
+   x0_shaft, y0_shaft = pd2c(r_shaft, beta/2)
 end
+%endif
+if x0_shaft == nil then
+   -- add air layer (inside) for heat transfer
+   h = 3.8
+   beta = 360*m.npols_gen/m.num_poles
 
-thcond = 0.026   -- air
-thcap = 1007
-get_spel_keys("sekeys")      -- Get all subregions of the model
-for i=1, #sekeys do
-  srkey =  get_spel_data("srkey", sekeys[i])
-  if srkey > 0 then
-    srname = get_sreg_data("name",srkey)
-    if srname == '    ' then
-      srkey = 0
-    end
-  end
-  if srkey == 0 then
-    elkeys = get_spel_data("elkeys", sekeys[i])
-    Ex, Ey = get_elem_data("xycp", elkeys[1])
-    r, phi = c2pr(Ex, Ey)
-    if r < m.rotor_rad then -- rotor only
-      def_mat_therm(Ex,Ey,'skyblue',1.12,thcond,thcap,1)
-    end
-  end
+   x0, y0 = pd2c(m.shaft_rad, m.zeroangl)
+   x1, y1 = pd2c(m.shaft_rad-h, m.zeroangl)
+   x2, y2 = pd2c(m.shaft_rad-h, beta+m.zeroangl)
+   x3, y3 = pd2c(m.shaft_rad, beta+m.zeroangl)
+   nc_line(x0, y0, x1, y1, 0)
+   nc_circle(x1, y1, x2, y2, 0)
+   nc_line(x2, y2, x3, y3, 0)
+   x0, y0 = pd2c(m.shaft_rad-h/2, beta/2+m.zeroangl)
+   create_mesh_se(x0, y0)
 end
---[[
--- add air layer (inside) for heat transfer
-h = 3.8
-beta = 360*m.npols_gen/m.num_poles
-x0, y0 = pd2c(dy2/2, m.zeroangl)
-x1, y1 = pd2c(dy2/2-h, m.zeroangl)
-x2, y2 = pd2c(dy2/2-h, beta+m.zeroangl)
-x3, y3 = pd2c(dy2/2, beta+m.zeroangl)
-nc_line(x0, y0, x1, y1, 0)
-nc_circle(x1, y1, x2, y2, 0)
-nc_line(x2, y2, x3, y3, 0)
-x0, y0 = pd2c(dy2/2-h/2, beta/2+m.zeroangl)
-create_mesh_se(x0, y0)
---]]
 %endif

femagtools/templates/mesh-airgap.mako

@@ -34,20 +34,26 @@ if not airgap_created then
     if inner_da_start == nil then
       inner_da_start = da2/2
     end
-    nc_line(r1, 0.0, r2, 0.0, 2)
+    x1, y1 = pr2c(inner_da_start, 0.0)
+    nc_line(x1, y1, r1, 0.0, 0.0)
 
     if outer_da_start == nil then
       outer_da_start = da1/2
     end
-    nc_line(r2, 0.0, outer_da_start, 0.0, 0)
+    x2, y2 = pr2c(outer_da_start, 0.0)
+    nc_line(r2, 0.0, x2, y2, 0.0)
 
     if m.tot_num_slot > m.num_sl_gen then
-      x4, y4 = pr2c(outer_da_start, alfa)
-      nc_line(x1, y1, x2, y2, 2)
-      nc_line(x4, y4, x2, y2, 0)
+      x3, y3 = pr2c(inner_da_end, alfa)
+      x4, y4 = pr2c(r1, alfa)
+      nc_line(x3, y3, x4, y4, 0, 0)
+
+      x3, y3 = pr2c(outer_da_end, alfa)
+      x4, y4 = pr2c(r2, alfa)
+      nc_line(x3, y3, x4, y4, 0, 0)
     end
 
-    x0, y0 = pr2c(r2-ag/6, alfa/2)
+    x0, y0 = pr2c(r1-ag/6, alfa/2)
     create_mesh_se(x0, y0)
     x0, y0 = pr2c(r2+ag/6, alfa/2)
     create_mesh_se(x0, y0)

femagtools/templates/prepare_thermal.mako

@@ -0,0 +1,354 @@
+-- prepare thermal model
+
+save_model('cont')
+state_of_problem("therm_static")
+function get_xy(elkey)
+  -- get element xy
+  local x, y
+  x, y = get_elem_data ( "xycp", elkey)
+  return x, y
+end
+
+function get_perm(elkey)
+  -- get element permeability
+  local x, y
+  x, y = get_elem_data ( "perm", elkey)
+  return x, y
+end
+
+function get_mag(elkey)
+  -- get element rem
+  local x, y
+  x, y = get_elem_data ("mag", elkey)
+  return x, y
+end
+
+function color_spel(spel_key, color)
+  local el
+  el = get_spel_data('elkeys', spel_key)
+  for i = 1, #el do
+     draw_elem(el[i], color, 0)
+  end
+end
+
+get_spel_keys("spel_keys")      -- Get all superelements of the model
+
+copper_spel = {}  -- red, 1
+stator_lam_spel = {} -- darkblue, 9
+rotor_lam_spel = {} -- darkblue, 9
+shaft_spel = {} ---- lightgrey, 16
+magnet_spel = {} -- black, 8
+outer_air = {} -- darkred, 10
+inner_air = {} -- lightgrey, 16
+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
+  is_winding = get_spel_data("condtyp", spel_keys[i])
+  els = get_spel_data('elkeys', spel_keys[i])
+  mag1, mag2 = get_mag(els[1])
+  perm1, perm2 = get_perm(els[1])
+  xc, yc = get_xy(els[1])
+  rc, pc = c2pd(xc, yc)
+
+  if (is_winding == 0) then
+     mctype = get_spel_data("mcvtyp", spel_keys[i])
+     if (mctype == 0) then -- no curve
+        -- check if it is air or magnet
+        if math.abs(mag1) > 0 or math.abs(mag2) > 0 then
+           -- is magnet
+           table.insert(magnet_spel, spel_keys[i])
+           def_mat_therm(xc, yc,'red', magn_density, magn_thcond, magn_thcap, 1)
+           color_spel(spel_keys[i], 8) -- black magnet
+        else
+           -- is air
+           if rc > dy1/2 then
+              -- outer air
+              table.insert(outer_air, spel_keys[i])
+%if model.get('htc_outer', 0):
+              def_heat_transfer(xc,yc,yellow,${model['htc_outer']}, 1.0)
+%endif
+              color_spel(spel_keys[i], 10) 
+
+           elseif rc > dy2/2 and rc < da2/2 then
+              -- rotor air
+              -- check is magnet pocket/air pocket/non_uniform airgap
+              table.insert(rotor_air, spel_keys[i])
+             --def_mat_therm(xc, yc, 'yellow', 1.12,0.026,1007, 1)
+              --color_spel(spel_keys[i], 3) -- yellow air
+
+           elseif rc > da1/2 and rc < dy1/2 then
+              -- 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) 
+              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.15,1007, 1)
+              add_to_sreg(xc, yc, 'Slot')
+              color_spel(spel_keys[i], 7) -- white air
+
+--           elseif rc > dy2/2 then
+              -- airgap
+--              table.insert(airgap, spel_keys[i])
+--              color_spel(spel_keys[i], 3) -- yellow air
+--              def_mat_therm(xc, yc, 'yellow', 1.19,0.15,1007, 1)
+
+            elseif rc < dy2/2 then
+              -- check if shaft or inner air
+              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)
+              end
+           end
+        end
+     else
+        -- check if it is stator / rotor
+        if rc > m.fc_radius then
+           color_spel(spel_keys[i], 9) -- magenta stator
+           table.insert(stator_lam_spel, spel_keys[i])
+           def_mat_therm(xc, yc, 'darkblue', stator_density, stator_thcond, stator_thcap, 1)
+        else
+           color_spel(spel_keys[i], 9) -- violet rotor
+           table.insert(rotor_lam_spel, spel_keys[i])
+           def_mat_therm(xc, yc, 'darkblue', rotor_density, rotor_thcond, rotor_thcap, 1)
+        end
+     end
+
+  else
+     table.insert(copper_spel, spel_keys[i])
+     def_mat_therm(xc, yc, 'green', conductor_density, conductor_thcond, conductor_thcap, 1)
+     color_spel(spel_keys[i], 1) --
+
+  end
+
+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 
+         return true
+      end 
+   end 
+   return false
+end 
+-- identify magnet air pocket 
+
+for i = 1, #magnet_spel do 
+    bkeys = get_spel_data ( "bndkeys", magnet_spel[i])
+    for j = 1, #bkeys do 
+         elks = get_node_data ( "elkeys", bkeys[j] )
+         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  
+                     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 
+
+
+
+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 
+      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 
+          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,1.15,1007, 1)
+          add_to_sreg(xc, yc, 'Slot')
+          color_spel(rotor_air[i], 7)
+        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 
+      rotor_air[i] = nil
+    end    
+  end 
+end 
+
+for i = 1, #rotor_air do 
+   if is_exist(magnet_pocket_spel, rotor_air[i]) then 
+     rotor_air[i] = nil
+   end 
+ end 
+ 
+
+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 
+
+
+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.get('slot_insul', 0.15)}
+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/templates/statorRotor3.mako

@@ -148,26 +148,7 @@ end
 %endif
 
 %if model.get('thcond', 0) and model.get('thcap', 0):
-thcond = ${model.get('thcond', 24)}
-thcap = ${model.get('thcap', 480)}
-density = ${model.get('density', 7700)}
--- Thermal properties: iron in yoke and tooth
-xst,yst = pd2c(m.inside_diam/2+1,m.zeroangl+0.1)    -- stator tooth
-xsj,ysj = pd2c(m.yoke_diam/2-1,m.zeroangl+0.1)  -- stator yoke
-def_mat_therm(xst,yst,'darkblue',density,thcond,thcap,1)
-def_mat_therm(xsj,ysj,'darkblue',density,thcond,thcap,1)
--- insulation
-if m.slot_indul > 0 then
-  thickness = 0.25 -- mm
-  thcond = 0.31
-  thcap = 1100
-  xip, yip = pr2c(m.slot_height + m.inside_diam/2-m.slot_indul/2,
-                  math.pi/m.tot_num_slot+m.zeroangl/180*math.pi)
-  def_insulation(xip,yip,'darkred',thickness,1340,thcond,thcap)
-end
--- air in slot opening
-xnl, ynl = pr2c(m.inside_diam/2+m.slot_h1/2, math.pi/m.tot_num_slot+m.zeroangl/180*math.pi)
-sl_cond = 0.15
-sl_cap = 1007
-def_mat_therm(xnl,ynl,'cyan',1.19,sl_cond,sl_cap,1)
+stator_thcond = ${model.get('thcond', 24)}
+stator_thcap = ${model.get('thcap', 480)}
+stator_density = ${model.get('density', 7700)}
 %endif